Loading...
Searching...
No Matches
G4RegularXTRadiator Class Reference
#include <G4RegularXTRadiator.hh>
Inheritance diagram for G4RegularXTRadiator:
Detailed Description
Definition at line 49 of file G4RegularXTRadiator.hh.
Constructor & Destructor Documentation
◆ G4RegularXTRadiator()
|
explicit |
Definition at line 41 of file G4RegularXTRadiator.cc.
44 :
45 G4VXTRenergyLoss(anEnvelope,foilMat,gasMat,a,b,n,processName)
46{
48
49 // Build energy and angular integral spectra of X-ray TR photons from
50 // a radiator
51
52 fAlphaPlate = 10000;
53 fAlphaGas = 1000;
55
56 // BuildTable() ;
57}
G4GLOB_DLL std::ostream G4cout
Definition: G4VXTRenergyLoss.hh:74
◆ ~G4RegularXTRadiator()
| G4RegularXTRadiator::~G4RegularXTRadiator | ( | ) |
Definition at line 61 of file G4RegularXTRadiator.cc.
62{
63 ;
64}
Member Function Documentation
◆ GetStackFactor()
|
overridevirtual |
Reimplemented from G4VXTRenergyLoss.
Definition at line 170 of file G4RegularXTRadiator.cc.
172{
173
174 // some gamma (10000/1000) like algorithm
175
176 G4double result, Za, Zb, Ma, Mb;
177
178 Za = GetPlateFormationZone(energy,gamma,varAngle);
179 Zb = GetGasFormationZone(energy,gamma,varAngle);
180
181 Ma = GetPlateLinearPhotoAbs(energy);
182 Mb = GetGasLinearPhotoAbs(energy);
183
184
187
190 G4complex H = Ha*Hb;
191
192 G4complex F1 = (1.0 - Ha)*(1.0 - Hb )/(1.0 - H)
194
195 G4complex F2 = (1.0-Ha)*(1.0-Ha)*Hb/(1.0-H)/(1.0-H)
196 * (1.0 - std::pow(H,fPlateNumber));
197
199
200 result = 2.0*std::real(R);
201
202 return result;
203
204 /*
205 // numerically stable but slow algorithm
206
207 G4double result, Qa, Qb, Q, aZa, bZb, aMa, bMb; // , D;
208
209 aZa = fPlateThick/GetPlateFormationZone(energy,gamma,varAngle);
210 bZb = fGasThick/GetGasFormationZone(energy,gamma,varAngle);
211 aMa = fPlateThick*GetPlateLinearPhotoAbs(energy);
212 bMb = fGasThick*GetGasLinearPhotoAbs(energy);
213 Qa = std::exp(-aMa);
214 Qb = std::exp(-bMb);
215 Q = Qa*Qb;
216 G4complex Ha( std::exp(-0.5*aMa)*std::cos(aZa),
217 -std::exp(-0.5*aMa)*std::sin(aZa) );
218 G4complex Hb( std::exp(-0.5*bMb)*std::cos(bZb),
219 -std::exp(-0.5*bMb)*std::sin(bZb) );
220 G4complex H = Ha*Hb;
221
222 G4complex Hs = conj(H);
223 D = 1.0 /( (1 - std::sqrt(Q))*(1 - std::sqrt(Q)) +
224 4*std::sqrt(Q)*std::sin(0.5*(aZa+bZb))*std::sin(0.5*(aZa+bZb)) );
225 G4complex F1 = (1.0 - Ha)*(1.0 - Hb)*(1.0 - Hs)
226 * G4double(fPlateNumber)*D;
227 G4complex F2 = (1.0-Ha)*(1.0-Ha)*Hb*(1.0-Hs)*(1.0-Hs)
228 * (1.0 - std::pow(H,fPlateNumber)) * D*D;
229 G4complex R = (F1 + F2)*OneInterfaceXTRdEdx(energy,gamma,varAngle);
230
231
232 G4complex S(0.,0.), c(1.,0.);
233 G4int k;
234 for(k = 1; k < fPlateNumber; k++)
235 {
236 c *= H;
237 S += ( G4double(fPlateNumber) - G4double(k) )*c;
238 }
239 G4complex R = (2.- Ha - 1./Ha)*S + (1. - Ha)*G4double(fPlateNumber);
240 R *= OneInterfaceXTRdEdx(energy,gamma,varAngle);
241 result = 2.0*std::real(R);
242 return result;
243 */
244}
G4double GetPlateLinearPhotoAbs(G4double)
Definition: G4VXTRenergyLoss.cc:1102
G4double GetGasFormationZone(G4double, G4double, G4double)
Definition: G4VXTRenergyLoss.cc:1122
G4complex OneInterfaceXTRdEdx(G4double energy, G4double gamma, G4double varAngle)
Definition: G4VXTRenergyLoss.cc:871
G4double GetPlateFormationZone(G4double, G4double, G4double)
Definition: G4VXTRenergyLoss.cc:1050
G4double GetGasLinearPhotoAbs(G4double)
Definition: G4VXTRenergyLoss.cc:1174
◆ SpectralXTRdEdx()
Reimplemented from G4VXTRenergyLoss.
Definition at line 70 of file G4RegularXTRadiator.cc.
71{
72 G4double result, sum = 0., tmp, cof1, cof2, cofMin, cofPHC, theta2, theta2k;
73 G4double aMa, bMb ,sigma, dump;
75
78 sigma = 0.5*(aMa + bMb);
79 dump = std::exp(-fPlateNumber*sigma);
81 cofPHC = 4*pi*hbarc;
83 cof1 = fPlateThick*tmp;
84 cof2 = fGasThick*tmp;
85
88 cofMin /= cofPHC;
89
91
92 // if (fGamma < 1200) kMin = G4int(cofMin); // 1200 ?
93 // else kMin = 1;
94
95
98
99 // tmp = (fPlateThick + fGasThick)*energy*fMaxThetaTR;
100 // tmp /= cofPHC;
101 // kMax = G4int(tmp);
102 // if(kMax < 0) kMax = 0;
103 // kMax += kMin;
104
105
107
108 // tmp /= fGamma;
109 // if( G4int(tmp) < kMin ) kMin = G4int(tmp);
110
112 {
115 }
117 {
119 result = (k - cof1)*(k - cof1)*(k + cof2)*(k + cof2);
120 // tmp = std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/result;
122 {
123 sum += 0.5*std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/result;
124 }
125 else
126 {
127 sum += std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/result;
128 }
129 theta2k = std::sqrt(theta2*std::abs(k-cofMin));
130
132 {
133 // G4cout<<"k = "<<k<<"; sqrt(theta2k) = "<<theta2k<<"; tmp = "<<std::sin(tmp)*std::sin(tmp)*std::abs(k-cofMin)/result
134 // <<"; sum = "<<sum<<G4endl;
137 }
138 }
139 result = 2*( cof1 + cof2 )*( cof1 + cof2 )*sum/energy;
140 // result *= ( 1 - std::exp(-0.5*fPlateNumber*sigma) )/( 1 - std::exp(-0.5*sigma) );
141 // fPlateNumber;
142 result *= ( 1 - dump + 2*dump*fPlateNumber );
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}
G4double energy(const ThreeVector &p, const G4double m)
Definition: G4INCLKinematicsUtils.cc:158
The documentation for this class was generated from the following files:
