| 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: G4eBremsstrahlungRelModel 32 // File name: G4eBremsstrahlungRelModel
33 // extention of standard G4eBre 33 // extention of standard G4eBremsstrahlungModel
34 // 34 //
35 // Author: Andreas Schaelicke 35 // Author: Andreas Schaelicke
36 // 36 //
37 // Creation date: 28.03.2008 37 // Creation date: 28.03.2008
38 // 38 //
39 // Modifications: 39 // Modifications:
40 // 40 //
41 // 15.07.18 introduced data structures to sto 41 // 15.07.18 introduced data structures to store LPM functions and element depen
42 // dent data for faster run-time com 42 // dent data for faster run-time computation (see more in .cc M.Novak)
43 // 43 //
44 // Class Description: 44 // Class Description:
45 // 45 //
46 // Implementation of energy loss for gamma emi 46 // Implementation of energy loss for gamma emission by electrons and
47 // positrons including an improved version of 47 // positrons including an improved version of the LPM effect
48 48
49 // ------------------------------------------- 49 // -------------------------------------------------------------------
50 // 50 //
51 51
52 #ifndef G4eBremsstrahlungRelModel_h 52 #ifndef G4eBremsstrahlungRelModel_h
53 #define G4eBremsstrahlungRelModel_h 1 53 #define G4eBremsstrahlungRelModel_h 1
54 54
55 #include "G4VEmModel.hh" 55 #include "G4VEmModel.hh"
56 #include <memory> <<
57 56
58 class G4ParticleChangeForLoss; 57 class G4ParticleChangeForLoss;
59 58
60 class G4eBremsstrahlungRelModel : public G4VEm 59 class G4eBremsstrahlungRelModel : public G4VEmModel {
61 60
62 public: 61 public:
63 62
64 explicit G4eBremsstrahlungRelModel(const G4P << 63 explicit G4eBremsstrahlungRelModel(const G4ParticleDefinition* p=0,
65 const G4S << 64 const G4String& nam="eBremLPM");
66 65
67 ~G4eBremsstrahlungRelModel() override; << 66 virtual ~G4eBremsstrahlungRelModel();
68 67
69 void Initialise(const G4ParticleDefinition*, << 68 virtual void Initialise(const G4ParticleDefinition*,
>> 69 const G4DataVector&) override;
70 70
71 void InitialiseLocal(const G4ParticleDefinit << 71 virtual void InitialiseLocal(const G4ParticleDefinition*,
72 G4VEmModel* masterModel << 72 G4VEmModel* masterModel) override;
73 << 73
74 G4double ComputeDEDXPerVolume(const G4Materi << 74 virtual G4double ComputeDEDXPerVolume(const G4Material*,
75 const G4Partic << 75 const G4ParticleDefinition*,
76 G4double ekin, << 76 G4double ekin,
77 G4double cutEn << 77 G4double cutEnergy) override;
78 << 78
79 G4double ComputeCrossSectionPerAtom(const G4 << 79 virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition*,
80 G4double << 80 G4double ekin,
81 G4double << 81 G4double zet,
82 G4double << 82 G4double,
83 G4double << 83 G4double cutEnergy,
84 G4double << 84 G4double maxEnergy = DBL_MAX) override;
85 << 85
86 void SampleSecondaries(std::vector<G4Dynamic << 86 virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
87 const G4MaterialCutsC << 87 const G4MaterialCutsCouple*,
88 const G4DynamicPartic << 88 const G4DynamicParticle*,
89 G4double cutEnergy, << 89 G4double cutEnergy,
90 G4double maxEnergy) o << 90 G4double maxEnergy) override;
91 << 91
92 void SetupForMaterial(const G4ParticleDefini << 92 virtual void SetupForMaterial(const G4ParticleDefinition*,
93 const G4Material*, G4d << 93 const G4Material*,
94 << 94 G4double) override;
95 G4double MinPrimaryEnergy(const G4Material*, << 95
96 const G4ParticleDe << 96 virtual G4double MinPrimaryEnergy(const G4Material*,
97 G4double cutEnergy << 97 const G4ParticleDefinition*,
>> 98 G4double cutEnergy) override;
98 99
99 protected: 100 protected:
100 101
101 virtual G4double ComputeDXSectionPerAtom(G4d 102 virtual G4double ComputeDXSectionPerAtom(G4double gammaEnergy);
102 103
103 void SetParticle(const G4ParticleDefinition* << 104 void SetParticle(const G4ParticleDefinition* p);
104 105
105 private: 106 private:
106 107
107 G4double ComputeBremLoss(G4double cutEnergy) 108 G4double ComputeBremLoss(G4double cutEnergy);
108 109
109 G4double ComputeXSectionPerAtom(G4double cut 110 G4double ComputeXSectionPerAtom(G4double cutEnergy);
110 111
111 G4double ComputeRelDXSectionPerAtom(G4double 112 G4double ComputeRelDXSectionPerAtom(G4double gammaEnergy);
112 113
113 // init special data per element i.e. per Z 114 // init special data per element i.e. per Z
114 void InitialiseElementData(); << 115 void InitialiseElementData();
115 116
116 // methods for initialisation and run-time e 117 // methods for initialisation and run-time evaluation of LPM functions:
117 void InitLPMFunctions(); << 118 void InitLPMFunctions();
118 <<
119 void ComputeLPMfunctions(G4double& funcXiS, <<
120 G4double& funcGS, <<
121 G4double& funcPhiS, <<
122 const G4double egam <<
123 <<
124 void GetLPMFunctions(G4double& lpmGs, <<
125 G4double& lpmPhis, <<
126 const G4double ss); <<
127 <<
128 void ComputeLPMGsPhis(G4double& funcGS, <<
129 G4double& funcPhiS, <<
130 const G4double varShat <<
131 119
>> 120 void ComputeLPMfunctions(G4double& funcXiS,
>> 121 G4double& funcGS,
>> 122 G4double& funcPhiS,
>> 123 const G4double egamma);
>> 124
>> 125 void GetLPMFunctions(G4double& lpmGs,
>> 126 G4double& lpmPhis,
>> 127 const G4double s);
>> 128
>> 129 void ComputeLPMGsPhis(G4double& funcGS,
>> 130 G4double& funcPhiS,
>> 131 const G4double varShat);
>> 132 //
132 // for evaluating screening related function 133 // for evaluating screening related functions
133 void ComputeScreeningFunctions(G4double& phi << 134 void ComputeScreeningFunctions(G4double& phi1,
134 G4double& phi << 135 G4double& phi1m2,
135 G4double& psi << 136 G4double& psi1,
136 G4double& psi << 137 G4double& psi1m2,
137 const G4doubl << 138 const G4double gam,
138 const G4doubl << 139 const G4double eps);
139 // hide assignment operator and cctr 140 // hide assignment operator and cctr
140 G4eBremsstrahlungRelModel& operator= << 141 G4eBremsstrahlungRelModel& operator=(const G4eBremsstrahlungRelModel& right) = delete;
141 (const G4eBremsstrahlungRelModel& right) = d <<
142 G4eBremsstrahlungRelModel(const G4eBremsstr 142 G4eBremsstrahlungRelModel(const G4eBremsstrahlungRelModel&) = delete;
143 143
144 private: <<
145 <<
146 G4bool fIsUseCompleteSc <<
147 G4bool fIsInitializer = <<
148 G4bool fUseLPM = true; <<
149 <<
150 protected: 144 protected:
151 145
152 G4bool fIsElectron = tr << 146 G4bool fIsElectron;
153 G4bool fIsScatOffElectr << 147 G4bool fIsScatOffElectron;
154 G4bool fIsLPMActive = f << 148 G4bool fIsLPMActive;
155 // << 149 //
156 G4int fCurrentIZ = 0; << 150 G4int fCurrentIZ;
157 const G4ParticleDefinition* fPrimaryParticle <<
158 G4ParticleDefinition* fGammaParticle = <<
159 G4ParticleChangeForLoss* fParticleChange <<
160 // cash 151 // cash
161 G4double fPrimaryParticle << 152 G4double fPrimaryParticleMass;
162 G4double fPrimaryKinEnerg << 153 G4double fPrimaryKinEnergy;
163 G4double fPrimaryTotalEne << 154 G4double fPrimaryTotalEnergy;
164 G4double fDensityFactor = << 155 G4double fDensityFactor;
165 G4double fDensityCorr = 0 << 156 G4double fDensityCorr;
166 G4double fLowestKinEnergy 157 G4double fLowestKinEnergy;
167 // scattering off electrons 158 // scattering off electrons
168 G4double fNucTerm = 0.; << 159 G4double fNucTerm;
169 G4double fSumTerm = 0.; << 160 G4double fSumTerm;
170 << 161 //
171 private: <<
172 <<
173 // LPM related members <<
174 G4double fLPMEnergyThresh <<
175 G4double fLPMEnergy; <<
176 <<
177 protected: <<
178 <<
179 static const G4double gBremFactor; 162 static const G4double gBremFactor;
180 static const G4double gMigdalConstant; 163 static const G4double gMigdalConstant;
>> 164 //
>> 165 const G4ParticleDefinition* fPrimaryParticle;
>> 166 G4ParticleDefinition* fGammaParticle;
>> 167 G4ParticleChangeForLoss* fParticleChange;
181 168
182 private: 169 private:
183 <<
184 static const G4int gMaxZet; 170 static const G4int gMaxZet;
185 // 171 //
186 static const G4double gLPMconstant; 172 static const G4double gLPMconstant;
187 // 173 //
188 static const G4double gXGL[8]; 174 static const G4double gXGL[8];
189 static const G4double gWGL[8]; 175 static const G4double gWGL[8];
190 static const G4double gFelLowZet[8]; 176 static const G4double gFelLowZet[8];
191 static const G4double gFinelLowZet[8]; 177 static const G4double gFinelLowZet[8];
192 // 178 //
193 struct ElementData { 179 struct ElementData {
194 /** @brief \f$ \ln(Z) \f$ */ 180 /** @brief \f$ \ln(Z) \f$ */
195 G4double fLogZ; 181 G4double fLogZ;
196 /** @brief \f$ \ln(Z)/3 + f_c \f$ */ 182 /** @brief \f$ \ln(Z)/3 + f_c \f$ */
197 G4double fFz; 183 G4double fFz;
198 /** @brief \f$ ((Fel-fc)+Finel*invZ)\f$ * 184 /** @brief \f$ ((Fel-fc)+Finel*invZ)\f$ */
199 G4double fZFactor1; 185 G4double fZFactor1;
200 /** @brief \f$ (Fel-fc)\f$ */ 186 /** @brief \f$ (Fel-fc)\f$ */
201 G4double fZFactor11; 187 G4double fZFactor11;
202 /** @brief \f$ (1.0+invZ)/12 \f$ */ 188 /** @brief \f$ (1.0+invZ)/12 \f$ */
203 G4double fZFactor2; 189 G4double fZFactor2;
204 // LPM variables 190 // LPM variables
205 G4double fVarS1; 191 G4double fVarS1;
206 G4double fILVarS1; 192 G4double fILVarS1;
207 G4double fILVarS1Cond; 193 G4double fILVarS1Cond;
208 // constant factors to the screening funct 194 // constant factors to the screening function evaluations
209 G4double fGammaFactor; 195 G4double fGammaFactor;
210 G4double fEpsilonFactor; 196 G4double fEpsilonFactor;
211 }; 197 };
212 // 198 //
213 struct LPMFuncs { 199 struct LPMFuncs {
214 LPMFuncs() : fIsInitialized(false), fISDel 200 LPMFuncs() : fIsInitialized(false), fISDelta(100.), fSLimit(2.) {}
215 G4bool fIsInitialized; 201 G4bool fIsInitialized;
216 G4double fISDelta; 202 G4double fISDelta;
217 G4double fSLimit; 203 G4double fSLimit;
218 std::vector<G4double> fLPMFuncG; 204 std::vector<G4double> fLPMFuncG;
219 std::vector<G4double> fLPMFuncPhi; 205 std::vector<G4double> fLPMFuncPhi;
220 }; 206 };
221 // 207 //
222 static std::shared_ptr<LPMFuncs> gLPMFuncs() << 208 static LPMFuncs gLPMFuncs;
223 static std::shared_ptr<std::vector<ElementDa << 209 static std::vector<ElementData*> gElementData;
224 std::shared_ptr<LPMFuncs> fLPMFuncs; << 210 //
225 std::shared_ptr<std::vector<ElementData*>> f << 211 G4bool fIsUseCompleteScreening;
>> 212 // LPM related members
>> 213 G4double fLPMEnergyThreshold;
>> 214 G4double fLPMEnergy;
>> 215
226 }; 216 };
227 217
228 #endif 218 #endif
229 219