Loading...
Searching...
No Matches
G4RegularXTRadiator Class Reference
#include <G4RegularXTRadiator.hh>
Inheritance diagram for G4RegularXTRadiator:
Detailed Description
Definition at line 50 of file G4RegularXTRadiator.hh.
Constructor & Destructor Documentation
◆ G4RegularXTRadiator()
| G4RegularXTRadiator::G4RegularXTRadiator | ( | G4LogicalVolume * | anEnvelope, |
| G4Material * | foilMat, | ||
| G4Material * | gasMat, | ||
| G4double | a, | ||
| G4double | b, | ||
| G4int | n, | ||
| const G4String & | processName = "XTRegularRadiator" |
||
| ) |
Definition at line 42 of file G4RegularXTRadiator.cc.
45 :
46 G4VXTRenergyLoss(anEnvelope,foilMat,gasMat,a,b,n,processName)
47{
49
50 // Build energy and angular integral spectra of X-ray TR photons from
51 // a radiator
52
53 fAlphaPlate = 10000;
54 fAlphaGas = 1000;
56
57 // BuildTable() ;
58}
G4DLLIMPORT std::ostream G4cout
Definition: G4VXTRenergyLoss.hh:75
◆ ~G4RegularXTRadiator()
| G4RegularXTRadiator::~G4RegularXTRadiator | ( | ) |
Definition at line 62 of file G4RegularXTRadiator.cc.
63{
64 ;
65}
Member Function Documentation
◆ GetStackFactor()
|
virtual |
Reimplemented from G4VXTRenergyLoss.
Definition at line 171 of file G4RegularXTRadiator.cc.
173{
174
175 // some gamma (10000/1000) like algorithm
176
177 G4double result, Za, Zb, Ma, Mb;
178
179 Za = GetPlateFormationZone(energy,gamma,varAngle);
180 Zb = GetGasFormationZone(energy,gamma,varAngle);
181
182 Ma = GetPlateLinearPhotoAbs(energy);
183 Mb = GetGasLinearPhotoAbs(energy);
184
185
188
191 G4complex H = Ha*Hb;
192
193 G4complex F1 = (1.0 - Ha)*(1.0 - Hb )/(1.0 - H)
195
196 G4complex F2 = (1.0-Ha)*(1.0-Ha)*Hb/(1.0-H)/(1.0-H)
197 * (1.0 - std::pow(H,fPlateNumber));
198
200
201 result = 2.0*std::real(R);
202
203 return result;
204
205 /*
206 // numerically stable but slow algorithm
207
208 G4double result, Qa, Qb, Q, aZa, bZb, aMa, bMb; // , D;
209
210 aZa = fPlateThick/GetPlateFormationZone(energy,gamma,varAngle);
211 bZb = fGasThick/GetGasFormationZone(energy,gamma,varAngle);
212 aMa = fPlateThick*GetPlateLinearPhotoAbs(energy);
213 bMb = fGasThick*GetGasLinearPhotoAbs(energy);
214 Qa = std::exp(-aMa);
215 Qb = std::exp(-bMb);
216 Q = Qa*Qb;
217 G4complex Ha( std::exp(-0.5*aMa)*std::cos(aZa),
218 -std::exp(-0.5*aMa)*std::sin(aZa) );
219 G4complex Hb( std::exp(-0.5*bMb)*std::cos(bZb),
220 -std::exp(-0.5*bMb)*std::sin(bZb) );
221 G4complex H = Ha*Hb;
222
223 G4complex Hs = conj(H);
224 D = 1.0 /( (1 - std::sqrt(Q))*(1 - std::sqrt(Q)) +
225 4*std::sqrt(Q)*std::sin(0.5*(aZa+bZb))*std::sin(0.5*(aZa+bZb)) );
226 G4complex F1 = (1.0 - Ha)*(1.0 - Hb)*(1.0 - Hs)
227 * G4double(fPlateNumber)*D;
228 G4complex F2 = (1.0-Ha)*(1.0-Ha)*Hb*(1.0-Hs)*(1.0-Hs)
229 * (1.0 - std::pow(H,fPlateNumber)) * D*D;
230 G4complex R = (F1 + F2)*OneInterfaceXTRdEdx(energy,gamma,varAngle);
231
232
233 G4complex S(0.,0.), c(1.,0.);
234 G4int k;
235 for(k = 1; k < fPlateNumber; k++)
236 {
237 c *= H;
238 S += ( G4double(fPlateNumber) - G4double(k) )*c;
239 }
240 G4complex R = (2.- Ha - 1./Ha)*S + (1. - Ha)*G4double(fPlateNumber);
241 R *= OneInterfaceXTRdEdx(energy,gamma,varAngle);
242 result = 2.0*std::real(R);
243 return result;
244 */
245}
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 G4RegularXTRadiator.cc.
72{
73 G4double result, sum = 0., tmp, cof1, cof2, cofMin, cofPHC, theta2, theta2k;
74 G4double aMa, bMb ,sigma, dump;
75 G4int k, kMax, kMin;
76
79 sigma = 0.5*(aMa + bMb);
80 dump = std::exp(-fPlateNumber*sigma);
82 cofPHC = 4*pi*hbarc;
84 cof1 = fPlateThick*tmp;
85 cof2 = fGasThick*tmp;
86
89 cofMin /= cofPHC;
90
92
93 // if (fGamma < 1200) kMin = G4int(cofMin); // 1200 ?
94 // else kMin = 1;
95
96
97 kMin = G4int(cofMin);
98 if (cofMin > kMin) kMin++;
99
100 // tmp = (fPlateThick + fGasThick)*energy*fMaxThetaTR;
101 // tmp /= cofPHC;
102 // kMax = G4int(tmp);
103 // if(kMax < 0) kMax = 0;
104 // kMax += kMin;
105
106
107 kMax = kMin + 49; // 19; // kMin + G4int(tmp);
108
109 // tmp /= fGamma;
110 // if( G4int(tmp) < kMin ) kMin = G4int(tmp);
111
113 {
116 }
117 for( k = kMin; k <= kMax; k++ )
118 {
120 result = (k - cof1)*(k - cof1)*(k + cof2)*(k + cof2);
121 // tmp = std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/result;
123 {
124 sum += 0.5*std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/result;
125 }
126 else
127 {
128 sum += std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/result;
129 }
130 theta2k = std::sqrt(theta2*std::abs(k-cofMin));
131
133 {
134 // G4cout<<"k = "<<k<<"; sqrt(theta2k) = "<<theta2k<<"; tmp = "<<std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/result
135 // <<"; sum = "<<sum<<G4endl;
138 }
139 }
140 result = 2*( cof1 + cof2 )*( cof1 + cof2 )*sum/energy;
141 // result *= ( 1 - std::exp(-0.5*fPlateNumber*sigma) )/( 1 - std::exp(-0.5*sigma) );
142 // fPlateNumber;
143 result *= ( 1 - dump + 2*dump*fPlateNumber );
144 /*
145 fEnergy = energy;
146 // G4Integrator<G4VXTRenergyLoss,G4double(G4VXTRenergyLoss::*)(G4double)> integral;
147 G4Integrator<G4TransparentRegXTRadiator,G4double(G4VXTRenergyLoss::*)(G4double)> integral;
148
149 tmp = integral.Legendre96(this,&G4VXTRenergyLoss::SpectralAngleXTRdEdx,
150 0.0,0.3*fMaxThetaTR) +
151 integral.Legendre96(this,&G4VXTRenergyLoss::SpectralAngleXTRdEdx,
152 0.3*fMaxThetaTR,0.6*fMaxThetaTR) +
153 integral.Legendre96(this,&G4VXTRenergyLoss::SpectralAngleXTRdEdx,
154 0.6*fMaxThetaTR,fMaxThetaTR) ;
155 result += tmp;
156 */
157 return result;
158}
The documentation for this class was generated from the following files:
