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