Loading...
Searching...
No Matches
G4TransparentRegXTRadiator Class Reference
#include <G4TransparentRegXTRadiator.hh>
Inheritance diagram for G4TransparentRegXTRadiator:
Detailed Description
Definition at line 51 of file G4TransparentRegXTRadiator.hh.
Constructor & Destructor Documentation
◆ G4TransparentRegXTRadiator()
| G4TransparentRegXTRadiator::G4TransparentRegXTRadiator | ( | G4LogicalVolume * | anEnvelope, |
| G4Material * | foilMat, | ||
| G4Material * | gasMat, | ||
| G4double | a, | ||
| G4double | b, | ||
| G4int | n, | ||
| const G4String & | processName = "TransparentRegXTRadiator" |
||
| ) |
Definition at line 42 of file G4TransparentRegXTRadiator.cc.
45 :
46 G4VXTRenergyLoss(anEnvelope,foilMat,gasMat,a,b,n,processName)
47{
50
51 // Build energy and angular integral spectra of X-ray TR photons from
52 // a radiator
53
54 fAlphaPlate = 10000;
55 fAlphaGas = 1000;
56
57 // BuildTable();
58}
G4DLLIMPORT std::ostream G4cout
Definition: G4VXTRenergyLoss.hh:75
◆ ~G4TransparentRegXTRadiator()
| G4TransparentRegXTRadiator::~G4TransparentRegXTRadiator | ( | ) |
Definition at line 62 of file G4TransparentRegXTRadiator.cc.
63{
64 ;
65}
Member Function Documentation
◆ GetStackFactor()
|
virtual |
Reimplemented from G4VXTRenergyLoss.
Definition at line 169 of file G4TransparentRegXTRadiator.cc.
171{
172 /*
173 G4double result, Za, Zb, Ma, Mb, sigma;
174
175 Za = GetPlateFormationZone(energy,gamma,varAngle);
176 Zb = GetGasFormationZone(energy,gamma,varAngle);
177 Ma = GetPlateLinearPhotoAbs(energy);
178 Mb = GetGasLinearPhotoAbs(energy);
179 sigma = Ma*fPlateThick + Mb*fGasThick;
180
181 G4complex Ca(1.0+0.5*fPlateThick*Ma/fAlphaPlate,fPlateThick/Za/fAlphaPlate);
182 G4complex Cb(1.0+0.5*fGasThick*Mb/fAlphaGas,fGasThick/Zb/fAlphaGas);
183
184 G4complex Ha = std::pow(Ca,-fAlphaPlate);
185 G4complex Hb = std::pow(Cb,-fAlphaGas);
186 G4complex H = Ha*Hb;
187 G4complex F1 = (1.0 - Ha)*(1.0 - Hb )/(1.0 - H)
188 * G4double(fPlateNumber) ;
189 G4complex F2 = (1.0-Ha)*(1.0-Ha)*Hb/(1.0-H)/(1.0-H)
190 * (1.0 - std::exp(-0.5*fPlateNumber*sigma)) ;
191 // *(1.0 - std::pow(H,fPlateNumber)) ;
192 G4complex R = (F1 + F2)*OneInterfaceXTRdEdx(energy,gamma,varAngle);
193 // G4complex R = F2*OneInterfaceXTRdEdx(energy,gamma,varAngle);
194 result = 2.0*std::real(R);
195 return result;
196 */
197 // numerically unstable result
198
199 G4double result, Qa, Qb, Q, aZa, bZb, aMa, bMb, D, sigma;
200
206 Qa = std::exp(-0.5*aMa);
207 Qb = std::exp(-0.5*bMb);
208 Q = Qa*Qb;
209
210 G4complex Ha( Qa*std::cos(aZa), -Qa*std::sin(aZa) );
211 G4complex Hb( Qb*std::cos(bZb), -Qb*std::sin(bZb) );
212 G4complex H = Ha*Hb;
213 G4complex Hs = conj(H);
214 D = 1.0 /( (1 - Q)*(1 - Q) +
215 4*Q*std::sin(0.5*(aZa + bZb))*std::sin(0.5*(aZa + bZb)) );
216 G4complex F1 = (1.0 - Ha)*(1.0 - Hb)*(1.0 - Hs)
218 G4complex F2 = (1.0 - Ha)*(1.0 - Ha)*Hb*(1.0 - Hs)*(1.0 - Hs)
219 // * (1.0 - std::pow(H,fPlateNumber)) * D*D;
220 * (1.0 - std::exp(-0.5*fPlateNumber*sigma)) * D*D;
222 result = 2.0*std::real(R);
223 return result;
224
225}
G4double GetPlateLinearPhotoAbs(G4double)
Definition: G4VXTRenergyLoss.cc:1085
G4double GetGasFormationZone(G4double, G4double, G4double)
Definition: G4VXTRenergyLoss.cc:1105
G4complex OneInterfaceXTRdEdx(G4double energy, G4double gamma, G4double varAngle)
Definition: G4VXTRenergyLoss.cc:858
G4double GetPlateFormationZone(G4double, G4double, G4double)
Definition: G4VXTRenergyLoss.cc:1033
G4double GetGasLinearPhotoAbs(G4double)
Definition: G4VXTRenergyLoss.cc:1157
◆ SpectralXTRdEdx()
Reimplemented from G4VXTRenergyLoss.
Definition at line 71 of file G4TransparentRegXTRadiator.cc.
72{
73 G4double result, sum = 0., tmp, cof1, cof2, cofMin, cofPHC, theta2, theta2k /*, aMa, bMb ,sigma*/;
74 G4int k, kMax, kMin;
75
76 //aMa = fPlateThick*GetPlateLinearPhotoAbs(energy);
77 //bMb = fGasThick*GetGasLinearPhotoAbs(energy);
78 //sigma = aMa + bMb;
79
80 cofPHC = 4*pi*hbarc;
82 cof1 = fPlateThick*tmp;
83 cof2 = fGasThick*tmp;
84
87 cofMin /= cofPHC;
88
90
91 // if (fGamma < 1200) kMin = G4int(cofMin); // 1200 ?
92 // else kMin = 1;
93
94
95 kMin = G4int(cofMin);
96 if (cofMin > kMin) kMin++;
97
98 // tmp = (fPlateThick + fGasThick)*energy*fMaxThetaTR;
99 // tmp /= cofPHC;
100 // kMax = G4int(tmp);
101 // if(kMax < 0) kMax = 0;
102 // kMax += kMin;
103
104
105 kMax = kMin + 49; // 19; // kMin + G4int(tmp);
106
107 // tmp /= fGamma;
108 // if( G4int(tmp) < kMin ) kMin = G4int(tmp);
109
111 {
114 }
115 for( k = kMin; k <= kMax; k++ )
116 {
118 result = (k - cof1)*(k - cof1)*(k + cof2)*(k + cof2);
119 // tmp = std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/result;
121 {
122 sum += 0.5*std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/result;
123 }
124 else
125 {
126 sum += std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/result;
127 }
128 theta2k = std::sqrt(theta2*std::abs(k-cofMin));
129
131 {
132 // G4cout<<"k = "<<k<<"; sqrt(theta2k) = "<<theta2k<<"; tmp = "<<std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/result
133 // <<"; sum = "<<sum<<G4endl;
136 }
137 }
138 result = 4*( cof1 + cof2 )*( cof1 + cof2 )*sum/energy;
139 // result *= ( 1 - std::exp(-0.5*fPlateNumber*sigma) )/( 1 - std::exp(-0.5*sigma) );
140 // fPlateNumber;
142 // +1-std::exp(-0.5*fPlateNumber*sigma);
143 /*
144 fEnergy = energy;
145 // G4Integrator<G4VXTRenergyLoss,G4double(G4VXTRenergyLoss::*)(G4double)> integral;
146 G4Integrator<G4TransparentRegXTRadiator,G4double(G4VXTRenergyLoss::*)(G4double)> integral;
147
148 tmp = integral.Legendre96(this,&G4VXTRenergyLoss::SpectralAngleXTRdEdx,
149 0.0,0.3*fMaxThetaTR) +
150 integral.Legendre96(this,&G4VXTRenergyLoss::SpectralAngleXTRdEdx,
151 0.3*fMaxThetaTR,0.6*fMaxThetaTR) +
152 integral.Legendre96(this,&G4VXTRenergyLoss::SpectralAngleXTRdEdx,
153 0.6*fMaxThetaTR,fMaxThetaTR) ;
154 result += tmp;
155 */
156 return result;
157}
The documentation for this class was generated from the following files:
