| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // 26 //
27 // ------------------------------------------- 27 // -------------------------------------------------------------------
28 // 28 //
29 // GEANT4 Class header file 29 // GEANT4 Class header file
30 // 30 //
31 // 31 //
32 // File name: G4EmCorrections 32 // File name: G4EmCorrections
33 // 33 //
34 // Author: Vladimir Ivanchenko 34 // Author: Vladimir Ivanchenko
35 // 35 //
36 // Creation date: 13.01.2005 36 // Creation date: 13.01.2005
37 // 37 //
38 // Modifications: 38 // Modifications:
39 // 28.04.2006 General cleanup, add finite size 39 // 28.04.2006 General cleanup, add finite size corrections (V.Ivanchenko)
40 // 13.05.2006 Add corrections for ion stopping 40 // 13.05.2006 Add corrections for ion stopping (V.Ivanhcenko)
41 // 20.05.2008 Removed Finite Size correction ( 41 // 20.05.2008 Removed Finite Size correction (V.Ivanchenko)
42 // 12.09.2008 Added inlined interfaces to effe 42 // 12.09.2008 Added inlined interfaces to effective charge (V.Ivanchenko)
43 // 19.04.2012 Fix reproducibility problem (A.R 43 // 19.04.2012 Fix reproducibility problem (A.Ribon)
44 // 44 //
45 // Class Description: 45 // Class Description:
46 // 46 //
47 // This class provides calculation of EM corre 47 // This class provides calculation of EM corrections to ionisation
48 // 48 //
49 49
50 // ------------------------------------------- 50 // -------------------------------------------------------------------
51 // 51 //
52 52
53 #ifndef G4EmCorrections_h 53 #ifndef G4EmCorrections_h
54 #define G4EmCorrections_h 1 54 #define G4EmCorrections_h 1
55 55
56 #include "globals.hh" 56 #include "globals.hh"
57 #include "G4ionEffectiveCharge.hh" 57 #include "G4ionEffectiveCharge.hh"
58 #include "G4Material.hh" 58 #include "G4Material.hh"
59 #include "G4ParticleDefinition.hh" 59 #include "G4ParticleDefinition.hh"
60 60
61 class G4VEmModel; 61 class G4VEmModel;
62 class G4PhysicsVector; 62 class G4PhysicsVector;
63 class G4IonTable; 63 class G4IonTable;
64 class G4MaterialCutsCouple; 64 class G4MaterialCutsCouple;
65 class G4PhysicsFreeVector; 65 class G4PhysicsFreeVector;
66 class G4Pow; 66 class G4Pow;
67 67
68 class G4EmCorrections 68 class G4EmCorrections
69 { 69 {
70 70
71 public: 71 public:
72 72
73 explicit G4EmCorrections(G4int verb); 73 explicit G4EmCorrections(G4int verb);
74 74
75 ~G4EmCorrections(); 75 ~G4EmCorrections();
76 76
77 G4double HighOrderCorrections(const G4Partic 77 G4double HighOrderCorrections(const G4ParticleDefinition*,
78 const G4Materi 78 const G4Material*,
79 const G4double 79 const G4double kineticEnergy,
80 const G4double 80 const G4double cutEnergy);
81 81
82 G4double IonHighOrderCorrections(const G4Par 82 G4double IonHighOrderCorrections(const G4ParticleDefinition*,
83 const G4Mat 83 const G4MaterialCutsCouple*,
84 const G4dou 84 const G4double kineticEnergy);
85 85
86 G4double ComputeIonCorrections(const G4Parti 86 G4double ComputeIonCorrections(const G4ParticleDefinition*,
87 const G4Mater 87 const G4Material*,
88 const G4doubl 88 const G4double kineticEnergy);
89 89
90 G4double IonBarkasCorrection(const G4Particl 90 G4double IonBarkasCorrection(const G4ParticleDefinition*,
91 const G4Materia 91 const G4Material*,
92 const G4double 92 const G4double kineticEnergy);
93 93
94 G4double Bethe(const G4ParticleDefinition*, 94 G4double Bethe(const G4ParticleDefinition*,
95 const G4Material*, 95 const G4Material*,
96 const G4double kineticEnergy) 96 const G4double kineticEnergy);
97 97
98 G4double SpinCorrection(const G4ParticleDefi 98 G4double SpinCorrection(const G4ParticleDefinition*,
99 const G4Material*, 99 const G4Material*,
100 const G4double kinet 100 const G4double kineticEnergy);
101 101
102 G4double KShellCorrection(const G4ParticleDe 102 G4double KShellCorrection(const G4ParticleDefinition*,
103 const G4Material*, 103 const G4Material*,
104 const G4double kin 104 const G4double kineticEnergy);
105 105
106 G4double LShellCorrection(const G4ParticleDe 106 G4double LShellCorrection(const G4ParticleDefinition*,
107 const G4Material*, 107 const G4Material*,
108 const G4double kin 108 const G4double kineticEnergy);
109 109
110 G4double ShellCorrection(const G4ParticleDef 110 G4double ShellCorrection(const G4ParticleDefinition*,
111 const G4Material*, 111 const G4Material*,
112 const G4double kine 112 const G4double kineticEnergy);
113 113
114 G4double ShellCorrectionSTD(const G4Particle 114 G4double ShellCorrectionSTD(const G4ParticleDefinition*,
115 const G4Material 115 const G4Material*,
116 const G4double k 116 const G4double kineticEnergy);
117 117
118 G4double DensityCorrection(const G4ParticleD 118 G4double DensityCorrection(const G4ParticleDefinition*,
119 const G4Material* 119 const G4Material*,
120 const G4double ki 120 const G4double kineticEnergy);
121 121
122 G4double BarkasCorrection(const G4ParticleDe 122 G4double BarkasCorrection(const G4ParticleDefinition*,
123 const G4Material*, 123 const G4Material*,
124 const G4double kin 124 const G4double kineticEnergy,
125 const G4bool isIni 125 const G4bool isInitialized = false);
126 126
127 G4double BlochCorrection(const G4ParticleDef 127 G4double BlochCorrection(const G4ParticleDefinition*,
128 const G4Material*, 128 const G4Material*,
129 const G4double kine 129 const G4double kineticEnergy,
130 const G4bool isInit 130 const G4bool isInitialized = false);
131 131
132 G4double MottCorrection(const G4ParticleDefi 132 G4double MottCorrection(const G4ParticleDefinition*,
133 const G4Material*, 133 const G4Material*,
134 const G4double kinet 134 const G4double kineticEnergy,
135 const G4bool isIniti 135 const G4bool isInitialized = false);
136 136
137 void AddStoppingData(const G4int Z, const G4 137 void AddStoppingData(const G4int Z, const G4int A,
138 const G4String& materia 138 const G4String& materialName,
139 G4PhysicsVector* dVecto 139 G4PhysicsVector* dVector);
140 140
141 void InitialiseForNewRun(); 141 void InitialiseForNewRun();
142 142
143 // effective charge correction using stoppin 143 // effective charge correction using stopping power data
144 G4double EffectiveChargeCorrection(const G4P 144 G4double EffectiveChargeCorrection(const G4ParticleDefinition*,
145 const G4M 145 const G4Material*,
146 const G4d 146 const G4double kineticEnergy);
147 147
148 // effective charge of an ion 148 // effective charge of an ion
149 inline G4double GetParticleCharge(const G4Pa 149 inline G4double GetParticleCharge(const G4ParticleDefinition*,
150 const G4Ma 150 const G4Material*,
151 const G4do 151 const G4double kineticEnergy);
152 152
153 inline 153 inline
154 G4double EffectiveChargeSquareRatio(const G4 154 G4double EffectiveChargeSquareRatio(const G4ParticleDefinition*,
155 const G4 155 const G4Material*,
156 const G4 156 const G4double kineticEnergy);
157 157
158 // ionisation models for ions 158 // ionisation models for ions
159 inline void SetIonisationModels(G4VEmModel* 159 inline void SetIonisationModels(G4VEmModel* mod1 = nullptr,
160 G4VEmModel* 160 G4VEmModel* mod2 = nullptr);
161 161
162 inline G4int GetNumberOfStoppingVectors() co 162 inline G4int GetNumberOfStoppingVectors() const;
163 163
164 inline void SetVerbose(G4int verb); 164 inline void SetVerbose(G4int verb);
165 165
166 // hide assignment operator 166 // hide assignment operator
167 G4EmCorrections & operator=(const G4EmCorrec 167 G4EmCorrections & operator=(const G4EmCorrections &right) = delete;
168 G4EmCorrections(const G4EmCorrections&) = de 168 G4EmCorrections(const G4EmCorrections&) = delete;
169 169
170 private: 170 private:
171 171
172 void Initialise(); 172 void Initialise();
173 173
174 void BuildCorrectionVector(); 174 void BuildCorrectionVector();
175 175
176 void SetupKinematics(const G4ParticleDefinit 176 void SetupKinematics(const G4ParticleDefinition*,
177 const G4Material*, 177 const G4Material*,
178 const G4double kineticE 178 const G4double kineticEnergy);
179 179
180 G4double KShell(const G4double theta, const 180 G4double KShell(const G4double theta, const G4double eta);
181 181
182 G4double LShell(const G4double theta, const 182 G4double LShell(const G4double theta, const G4double eta);
183 183
184 G4int Index(const G4double x, const G4double 184 G4int Index(const G4double x, const G4double* y, const G4int n) const;
185 185
186 G4double Value(const G4double xv, const G4do 186 G4double Value(const G4double xv, const G4double x1, const G4double x2,
187 const G4double y1, const G4do 187 const G4double y1, const G4double y2) const;
188 188
189 G4double Value2(const G4double xv, const G4d 189 G4double Value2(const G4double xv, const G4double yv,
190 const G4double x1, const G4d 190 const G4double x1, const G4double x2,
191 const G4double y1, const G4d 191 const G4double y1, const G4double y2,
192 const G4double z11, const G4 192 const G4double z11, const G4double z21,
193 const G4double z12, const G4double z22) 193 const G4double z12, const G4double z22) const;
194 194
195 G4Pow* g4calc; 195 G4Pow* g4calc;
196 G4IonTable* ionTable; 196 G4IonTable* ionTable;
197 197
198 const G4ParticleDefinition* particle = nullp 198 const G4ParticleDefinition* particle = nullptr;
199 const G4ParticleDefinition* curParticle = nu 199 const G4ParticleDefinition* curParticle = nullptr;
200 const G4Material* material = nullptr; 200 const G4Material* material = nullptr;
201 const G4Material* curMaterial = nullptr; 201 const G4Material* curMaterial = nullptr;
202 const G4ElementVector* theElementVector = nu 202 const G4ElementVector* theElementVector = nullptr;
203 const G4double* atomDensity = nullptr; 203 const G4double* atomDensity = nullptr;
204 204
205 G4PhysicsVector* curVector = nullptr; 205 G4PhysicsVector* curVector = nullptr;
206 206
207 G4VEmModel* ionLEModel = nullptr; 207 G4VEmModel* ionLEModel = nullptr;
208 G4VEmModel* ionHEModel = nullptr; 208 G4VEmModel* ionHEModel = nullptr;
209 209
210 G4double kinEnergy = 0.0; 210 G4double kinEnergy = 0.0;
211 G4double mass = 0.0; 211 G4double mass = 0.0;
212 G4double massFactor = 1.0; 212 G4double massFactor = 1.0;
213 G4double tau = 0.0; 213 G4double tau = 0.0;
214 G4double gamma = 1.0; 214 G4double gamma = 1.0;
215 G4double bg2 = 0.0; 215 G4double bg2 = 0.0;
216 G4double beta2 = 0.0; 216 G4double beta2 = 0.0;
217 G4double beta = 0.0; 217 G4double beta = 0.0;
218 G4double ba2 = 0.0; 218 G4double ba2 = 0.0;
219 G4double tmax = 0.0; 219 G4double tmax = 0.0;
220 G4double charge = 0.0; 220 G4double charge = 0.0;
221 G4double q2 = 0.0; 221 G4double q2 = 0.0;
222 G4double eth; 222 G4double eth;
223 G4double eCorrMin; 223 G4double eCorrMin;
224 G4double eCorrMax; 224 G4double eCorrMax;
225 225
226 std::size_t ncouples = 0; 226 std::size_t ncouples = 0;
227 std::size_t idxBarkas = 0; 227 std::size_t idxBarkas = 0;
228 G4int nK = 20; 228 G4int nK = 20;
229 G4int nL = 26; 229 G4int nL = 26;
230 G4int nEtaK = 29; 230 G4int nEtaK = 29;
231 G4int nEtaL = 28; 231 G4int nEtaL = 28;
232 G4int nbinCorr = 52; 232 G4int nbinCorr = 52;
233 G4int numberOfElements = 0; 233 G4int numberOfElements = 0;
234 234
235 // Ion stopping data 235 // Ion stopping data
236 G4int nIons = 0; 236 G4int nIons = 0;
237 G4int idx = 0; 237 G4int idx = 0;
238 G4int currentZ = 0; 238 G4int currentZ = 0;
239 239
240 G4int verbose; 240 G4int verbose;
241 G4bool isInitializer = false; 241 G4bool isInitializer = false;
242 242
243 std::vector<G4int> Zion; 243 std::vector<G4int> Zion;
244 std::vector<G4int> Aion; 244 std::vector<G4int> Aion;
245 std::vector<G4String> materialName; 245 std::vector<G4String> materialName;
246 std::vector<const G4ParticleDefinition*> ion 246 std::vector<const G4ParticleDefinition*> ionList;
247 247
248 std::map< G4int, std::vector<G4double> > thc 248 std::map< G4int, std::vector<G4double> > thcorr;
249 249
250 std::vector<const G4Material*> currmat; 250 std::vector<const G4Material*> currmat;
251 std::vector<const G4Material*> materialList; 251 std::vector<const G4Material*> materialList;
252 std::vector<G4PhysicsVector*> stopData; 252 std::vector<G4PhysicsVector*> stopData;
253 253
254 G4ionEffectiveCharge effCharge; 254 G4ionEffectiveCharge effCharge;
255 255
256 static const G4double ZD[11]; 256 static const G4double ZD[11];
257 static const G4double UK[20]; 257 static const G4double UK[20];
258 static const G4double VK[20]; 258 static const G4double VK[20];
259 static G4double ZK[20]; 259 static G4double ZK[20];
260 static const G4double Eta[29]; 260 static const G4double Eta[29];
261 static G4double CK[20][29]; 261 static G4double CK[20][29];
262 static G4double CL[26][28]; 262 static G4double CL[26][28];
263 static const G4double UL[26]; 263 static const G4double UL[26];
264 static G4double VL[26]; 264 static G4double VL[26];
265 265
266 static G4double sWmaxBarkas; 266 static G4double sWmaxBarkas;
267 static G4PhysicsFreeVector* sBarkasCorr; 267 static G4PhysicsFreeVector* sBarkasCorr;
268 static G4PhysicsFreeVector* sThetaK; 268 static G4PhysicsFreeVector* sThetaK;
269 static G4PhysicsFreeVector* sThetaL; 269 static G4PhysicsFreeVector* sThetaL;
270 }; 270 };
271 271
272 inline G4int 272 inline G4int
273 G4EmCorrections::Index(const G4double x, const 273 G4EmCorrections::Index(const G4double x, const G4double* y, const G4int n) const
274 { 274 {
275 G4int iddd = n-1; 275 G4int iddd = n-1;
276 // Loop checking, 03-Aug-2015, Vladimir Ivan 276 // Loop checking, 03-Aug-2015, Vladimir Ivanchenko
277 do {--iddd;} while (iddd>0 && x<y[iddd]); 277 do {--iddd;} while (iddd>0 && x<y[iddd]);
278 return iddd; 278 return iddd;
279 } 279 }
280 280
281 inline G4double G4EmCorrections::Value(const G 281 inline G4double G4EmCorrections::Value(const G4double xv, const G4double x1,
282 const G 282 const G4double x2,
283 const G 283 const G4double y1, const G4double y2) const
284 { 284 {
285 return y1 + (y2 - y1)*(xv - x1)/(x2 - x1); 285 return y1 + (y2 - y1)*(xv - x1)/(x2 - x1);
286 } 286 }
287 287
288 inline G4double G4EmCorrections::Value2(const 288 inline G4double G4EmCorrections::Value2(const G4double xv, const G4double yv,
289 const 289 const G4double x1, const G4double x2,
290 const 290 const G4double y1, const G4double y2,
291 const 291 const G4double z11, const G4double z21,
292 const 292 const G4double z12, const G4double z22) const
293 { 293 {
294 return ( z11*(x2-xv)*(y2-yv) + z22*(xv-x1)*( 294 return ( z11*(x2-xv)*(y2-yv) + z22*(xv-x1)*(yv-y1) +
295 z12*(x2-xv)*(yv-y1) + z21*(xv-x1)*(y2-yv) 295 z12*(x2-xv)*(yv-y1) + z21*(xv-x1)*(y2-yv) )
296 / ((x2-x1)*(y2-y1)); 296 / ((x2-x1)*(y2-y1));
297 } 297 }
298 298
299 inline void 299 inline void
300 G4EmCorrections::SetIonisationModels(G4VEmMode 300 G4EmCorrections::SetIonisationModels(G4VEmModel* mod1, G4VEmModel* mod2)
301 { 301 {
302 if(nullptr != mod1) { ionLEModel = mod1; } 302 if(nullptr != mod1) { ionLEModel = mod1; }
303 if(nullptr != mod2) { ionHEModel = mod2; } 303 if(nullptr != mod2) { ionHEModel = mod2; }
304 } 304 }
305 305
306 inline G4int G4EmCorrections::GetNumberOfStopp 306 inline G4int G4EmCorrections::GetNumberOfStoppingVectors() const
307 { 307 {
308 return nIons; 308 return nIons;
309 } 309 }
310 310
311 inline G4double 311 inline G4double
312 G4EmCorrections::GetParticleCharge(const G4Par 312 G4EmCorrections::GetParticleCharge(const G4ParticleDefinition* p,
313 const G4Mat 313 const G4Material* mat,
314 const G4dou 314 const G4double kineticEnergy)
315 { 315 {
316 return effCharge.EffectiveCharge(p,mat,kinet 316 return effCharge.EffectiveCharge(p,mat,kineticEnergy);
317 } 317 }
318 318
319 inline G4double 319 inline G4double
320 G4EmCorrections::EffectiveChargeSquareRatio(co 320 G4EmCorrections::EffectiveChargeSquareRatio(const G4ParticleDefinition* p,
321 co 321 const G4Material* mat,
322 co 322 const G4double kineticEnergy)
323 { 323 {
324 return effCharge.EffectiveChargeSquareRatio( 324 return effCharge.EffectiveChargeSquareRatio(p,mat,kineticEnergy);
325 } 325 }
326 326
327 inline void G4EmCorrections::SetVerbose(G4int 327 inline void G4EmCorrections::SetVerbose(G4int verb)
328 { 328 {
329 verbose = verb; 329 verbose = verb;
330 } 330 }
331 331
332 //....oooOO0OOooo........oooOO0OOooo........oo 332 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
333 333
334 #endif 334 #endif
335 335