Loading...
Searching...
No Matches
G4XTRTransparentRegRadModel Class Reference
#include <G4XTRTransparentRegRadModel.hh>
Inheritance diagram for G4XTRTransparentRegRadModel:
Detailed Description
Definition at line 48 of file G4XTRTransparentRegRadModel.hh.
Constructor & Destructor Documentation
◆ G4XTRTransparentRegRadModel()
| G4XTRTransparentRegRadModel::G4XTRTransparentRegRadModel | ( | G4LogicalVolume * | anEnvelope, |
| G4Material * | foilMat, | ||
| G4Material * | gasMat, | ||
| G4double | a, | ||
| G4double | b, | ||
| G4int | n, | ||
| const G4String & | processName = "XTRTransparentRegRadModel" |
||
| ) |
Definition at line 40 of file G4XTRTransparentRegRadModel.cc.
43 :
44 G4VXTRenergyLoss(anEnvelope,foilMat,gasMat,a,b,n,processName)
45{
47
48 // Build energy and angular integral spectra of X-ray TR photons from
49 // a radiator
51 fAlphaPlate = 10000;
52 fAlphaGas = 1000;
53
54 // BuildTable();
55}
G4DLLIMPORT std::ostream G4cout
Definition: G4VXTRenergyLoss.hh:75
◆ ~G4XTRTransparentRegRadModel()
| G4XTRTransparentRegRadModel::~G4XTRTransparentRegRadModel | ( | ) |
Definition at line 59 of file G4XTRTransparentRegRadModel.cc.
60{
61 ;
62}
Member Function Documentation
◆ GetStackFactor()
|
virtual |
Reimplemented from G4VXTRenergyLoss.
Definition at line 145 of file G4XTRTransparentRegRadModel.cc.
147{
148 /*
149 G4double result, Za, Zb, Ma, Mb, sigma;
150
151 Za = GetPlateFormationZone(energy,gamma,varAngle);
152 Zb = GetGasFormationZone(energy,gamma,varAngle);
153 Ma = GetPlateLinearPhotoAbs(energy);
154 Mb = GetGasLinearPhotoAbs(energy);
155 sigma = Ma*fPlateThick + Mb*fGasThick;
156
157 G4complex Ca(1.0+0.5*fPlateThick*Ma/fAlphaPlate,fPlateThick/Za/fAlphaPlate);
158 G4complex Cb(1.0+0.5*fGasThick*Mb/fAlphaGas,fGasThick/Zb/fAlphaGas);
159
160 G4complex Ha = std::pow(Ca,-fAlphaPlate);
161 G4complex Hb = std::pow(Cb,-fAlphaGas);
162 G4complex H = Ha*Hb;
163 G4complex F1 = (1.0 - Ha)*(1.0 - Hb )/(1.0 - H)
164 * G4double(fPlateNumber) ;
165 G4complex F2 = (1.0-Ha)*(1.0-Ha)*Hb/(1.0-H)/(1.0-H)
166 * (1.0 - std::exp(-0.5*fPlateNumber*sigma)) ;
167 // *(1.0 - std::pow(H,fPlateNumber)) ;
168 G4complex R = (F1 + F2)*OneInterfaceXTRdEdx(energy,gamma,varAngle);
169 // G4complex R = F2*OneInterfaceXTRdEdx(energy,gamma,varAngle);
170 result = 2.0*std::real(R);
171 return result;
172 */
173 // numerically unstable result
174
175 G4double result, Qa, Qb, Q, aZa, bZb, aMa, bMb, D, sigma;
176
182 Qa = std::exp(-0.5*aMa);
183 Qb = std::exp(-0.5*bMb);
184 Q = Qa*Qb;
185
186 G4complex Ha( Qa*std::cos(aZa), -Qa*std::sin(aZa) );
187 G4complex Hb( Qb*std::cos(bZb), -Qb*std::sin(bZb) );
188 G4complex H = Ha*Hb;
189 G4complex Hs = conj(H);
190 D = 1.0 /( (1 - Q)*(1 - Q) +
191 4*Q*std::sin(0.5*(aZa + bZb))*std::sin(0.5*(aZa + bZb)) );
192 G4complex F1 = (1.0 - Ha)*(1.0 - Hb)*(1.0 - Hs)
194 G4complex F2 = (1.0 - Ha)*(1.0 - Ha)*Hb*(1.0 - Hs)*(1.0 - Hs)
195 // * (1.0 - std::pow(H,fPlateNumber)) * D*D;
196 * (1.0 - std::exp(-0.5*fPlateNumber*sigma)) * D*D;
198 result = 2.0*std::real(R);
199 return result;
200
201}
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 68 of file G4XTRTransparentRegRadModel.cc.
69{
70 G4double result, sum = 0., tmp, cof1, cof2, cofMin, cofPHC,aMa, bMb, sigma;
71 G4int k, kMax, kMin;
72
73 aMa = GetPlateLinearPhotoAbs(energy);
74 bMb = GetGasLinearPhotoAbs(energy);
75
77 {
78 aMa += GetPlateCompton(energy);
79 bMb += GetGasCompton(energy);
80 }
81 aMa *= fPlateThick;
82 bMb *= fGasThick;
83
84 sigma = aMa + bMb;
85
86 cofPHC = 4*pi*hbarc;
88 cof1 = fPlateThick*tmp;
89 cof2 = fGasThick*tmp;
90
93 cofMin /= cofPHC;
94
95 // if (fGamma < 1200) kMin = G4int(cofMin); // 1200 ?
96 // else kMin = 1;
97
98
99 kMin = G4int(cofMin);
100 if (cofMin > kMin) kMin++;
101
102 // tmp = (fPlateThick + fGasThick)*energy*fMaxThetaTR;
103 // tmp /= cofPHC;
104 // kMax = G4int(tmp);
105 // if(kMax < 0) kMax = 0;
106 // kMax += kMin;
107
108
109 kMax = kMin + 19; // 5; // 9; // kMin + G4int(tmp);
110
111 // tmp /= fGamma;
112 // if( G4int(tmp) < kMin ) kMin = G4int(tmp);
113 // G4cout<<"kMin = "<<kMin<<"; kMax = "<<kMax<<G4endl;
114
115 for( k = kMin; k <= kMax; k++ )
116 {
118 result = (k - cof1)*(k - cof1)*(k + cof2)*(k + cof2);
119
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 // G4cout<<"k = "<<k<<"; sum = "<<sum<<G4endl;
129 }
130 result = 4.*( cof1 + cof2 )*( cof1 + cof2 )*sum/energy;
131 result *= ( 1. - std::exp(-fPlateNumber*sigma) )/( 1. - std::exp(-sigma) );
132 return result;
133}
G4double GetPlateCompton(G4double)
Definition: G4VXTRenergyLoss.cc:1252
The documentation for this class was generated from the following files:
