| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26
27 #include "G4XTRRegularRadModel.hh"
28
29 #include "G4PhysicalConstants.hh"
30
31 //////////////////////////////////////////////
32 // Constructor, destructor
33 G4XTRRegularRadModel::G4XTRRegularRadModel(G4L
34 G4M
35 G4M
36 G4d
37 con
38 : G4VXTRenergyLoss(anEnvelope, foilMat, gasM
39 {
40 G4cout << " XTR Regular discrete radiator mo
41
42 fExitFlux = true;
43 }
44
45 //////////////////////////////////////////////
46 G4XTRRegularRadModel::~G4XTRRegularRadModel()
47
48 //////////////////////////////////////////////
49 void G4XTRRegularRadModel::ProcessDescription(
50 {
51 out << "Describes X-ray transition radiation
52 "plates\n"
53 "fixed.\n";
54 }
55
56 //////////////////////////////////////////////
57 G4double G4XTRRegularRadModel::SpectralXTRdEdx
58 {
59 static constexpr G4double cofPHC = 4. * pi *
60 G4double result, sum = 0., tmp, cof1, cof2,
61 G4double aMa, bMb, sigma, dump;
62 G4int k, kMax, kMin;
63
64 aMa = fPlateThick * GetPlateLinearPhotoAbs
65 bMb = fGasThick * GetGasLinearPhotoAbs(ene
66 sigma = 0.5 * (aMa + bMb);
67 dump = std::exp(-fPlateNumber * sigma);
68 if(verboseLevel > 2)
69 G4cout << " dump = " << dump << G4endl;
70 tmp = (fSigma1 - fSigma2) / cofPHC / energy
71 cof1 = fPlateThick * tmp;
72 cof2 = fGasThick * tmp;
73
74 cofMin = energy * (fPlateThick + fGasThick)
75 cofMin += (fPlateThick * fSigma1 + fGasThick
76 cofMin /= cofPHC;
77
78 theta2 = cofPHC / (energy * (fPlateThick + f
79
80 kMin = G4int(cofMin);
81 if(cofMin > kMin)
82 kMin++;
83
84 kMax = kMin + 49;
85
86 if(verboseLevel > 2)
87 {
88 G4cout << cof1 << " " << cof2 << "
89 G4cout << "kMin = " << kMin << "; kMax
90 }
91 for(k = kMin; k <= kMax; ++k)
92 {
93 tmp = pi * fPlateThick * (k + cof2) / (
94 result = (k - cof1) * (k - cof1) * (k + co
95 if(k == kMin && kMin == G4int(cofMin))
96 {
97 sum +=
98 0.5 * std::sin(tmp) * std::sin(tmp) *
99 }
100 else
101 {
102 sum += std::sin(tmp) * std::sin(tmp) * s
103 }
104 theta2k = std::sqrt(theta2 * std::abs(k -
105
106 if(verboseLevel > 2)
107 {
108 G4cout << k << " " << theta2k << "
109 << std::sin(tmp) * std::sin(tmp)
110 << " " << sum << G4endl;
111 }
112 }
113 result = 2 * (cof1 + cof2) * (cof1 + cof2) *
114 result *= dump * (-1 + dump + 2 * fPlateNumb
115
116 return result;
117 }
118
119 //////////////////////////////////////////////
120 // Approximation for radiator interference fac
121 // fully Regular radiator. The plate and gas g
122 // The mean values of the plate and gas gap th
123 // are supposed to be about XTR formation zone
124 // mean absorption length of XTR photons in co
125 G4double G4XTRRegularRadModel::GetStackFactor(
126
127 {
128 G4double aZa = fPlateThick / GetPlateFormati
129 G4double bZb = fGasThick / GetGasFormationZo
130
131 G4double aMa = fPlateThick * GetPlateLinearP
132 G4double bMb = fGasThick * GetGasLinearPhoto
133
134 G4double Qa = std::exp(-aMa);
135 G4double Qb = std::exp(-bMb);
136 G4double Q = Qa * Qb;
137
138 G4complex Ha(std::exp(-0.5 * aMa) * std::cos
139 -std::exp(-0.5 * aMa) * std::si
140
141 G4complex Hb(std::exp(-0.5 * bMb) * std::cos
142 -std::exp(-0.5 * bMb) * std::si
143
144 G4complex H = Ha * Hb;
145 G4complex Hs = std::conj(H);
146
147 G4complex F2 = (1.0 - Ha) * (Qa - Ha) * Hb *
148 F2 *= std::pow(Q, G4double(fPlateNumber)) -
149
150 G4double result = (1. - std::pow(Q, G4double
151 result *= (1. - Qa) * (1. + Qa - 2. * std::s
152 result /= (1. - std::sqrt(Q)) * (1. - std::s
153 4. * std::sqrt(Q) * std::sin(0.5 *
154 std::sin(0.5 * (aZa + bZb));
155
156 G4double I2 = 1.;
157 I2 /= (1. - std::sqrt(Q)) * (1. - std::sqrt(
158 4. * std::sqrt(Q) * std::sin(0.5 * (aZ
159 std::sin(0.5 * (aZa + bZb));
160
161 I2 /= Q * ((std::sqrt(Q) - std::cos(aZa + bZ
162 (std::sqrt(Q) - std::cos(aZa +
163 std::sin(aZa + bZb) * std::sin(aZ
164
165 G4complex stack = 2. * I2 * F2;
166 stack += result;
167 stack *= OneInterfaceXTRdEdx(energy, gamma,
168
169 return std::real(stack);
170 }
171