| 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 // GEANT4 Class header file 28 // GEANT4 Class header file
29 // 29 //
30 // 30 //
31 // File name: G4VEmModel 31 // File name: G4VEmModel
32 // 32 //
33 // Author: Vladimir Ivanchenko 33 // Author: Vladimir Ivanchenko
34 // 34 //
35 // Creation date: 03.01.2002 35 // Creation date: 03.01.2002
36 // 36 //
37 // Modifications: 37 // Modifications:
38 // 38 //
39 // 23-12-02 V.Ivanchenko change interface befo 39 // 23-12-02 V.Ivanchenko change interface before move to cut per region
40 // 24-01-03 Cut per region (V.Ivanchenko) 40 // 24-01-03 Cut per region (V.Ivanchenko)
41 // 13-02-03 Add name (V.Ivanchenko) 41 // 13-02-03 Add name (V.Ivanchenko)
42 // 25-02-03 Add sample theta and displacement 42 // 25-02-03 Add sample theta and displacement (V.Ivanchenko)
43 // 23-07-03 Replace G4Material by G4MaterialCu 43 // 23-07-03 Replace G4Material by G4MaterialCutCouple in dE/dx and CrossSection
44 // calculation (V.Ivanchenko) 44 // calculation (V.Ivanchenko)
45 // 01-03-04 L.Urban signature changed in Sampl 45 // 01-03-04 L.Urban signature changed in SampleCosineTheta
46 // 23-04-04 L.urban signature of SampleCosineT 46 // 23-04-04 L.urban signature of SampleCosineTheta changed back
47 // 17-11-04 Add method CrossSectionPerAtom (V. 47 // 17-11-04 Add method CrossSectionPerAtom (V.Ivanchenko)
48 // 14-03-05 Reduce number of pure virtual meth 48 // 14-03-05 Reduce number of pure virtual methods and make inline part
49 // separate (V.Ivanchenko) 49 // separate (V.Ivanchenko)
50 // 24-03-05 Remove IsInCharge and add G4VParti 50 // 24-03-05 Remove IsInCharge and add G4VParticleChange in the constructor (VI)
51 // 08-04-05 Major optimisation of internal int 51 // 08-04-05 Major optimisation of internal interfaces (V.Ivantchenko)
52 // 15-04-05 optimize internal interface for ms 52 // 15-04-05 optimize internal interface for msc (V.Ivanchenko)
53 // 08-05-05 A -> N (V.Ivanchenko) 53 // 08-05-05 A -> N (V.Ivanchenko)
54 // 25-07-05 Move constructor and destructor to 54 // 25-07-05 Move constructor and destructor to the body (V.Ivanchenko)
55 // 02-02-06 ComputeCrossSectionPerAtom: defaul 55 // 02-02-06 ComputeCrossSectionPerAtom: default value A=0. (mma)
56 // 06-02-06 add method ComputeMeanFreePath() ( 56 // 06-02-06 add method ComputeMeanFreePath() (mma)
57 // 07-03-06 Optimize msc methods (V.Ivanchenko 57 // 07-03-06 Optimize msc methods (V.Ivanchenko)
58 // 29-06-06 Add member currentElement and Get/ 58 // 29-06-06 Add member currentElement and Get/Set methods (V.Ivanchenko)
59 // 29-10-07 Added SampleScattering (V.Ivanchen 59 // 29-10-07 Added SampleScattering (V.Ivanchenko)
60 // 15-07-08 Reorder class members and improve 60 // 15-07-08 Reorder class members and improve comments (VI)
61 // 21-07-08 Added vector of G4ElementSelector 61 // 21-07-08 Added vector of G4ElementSelector and methods to use it (VI)
62 // 12-09-08 Added methods GetParticleCharge, G 62 // 12-09-08 Added methods GetParticleCharge, GetChargeSquareRatio,
63 // CorrectionsAlongStep, ActivateNucl 63 // CorrectionsAlongStep, ActivateNuclearStopping (VI)
64 // 16-02-09 Moved implementations of virtual m 64 // 16-02-09 Moved implementations of virtual methods to source (VI)
65 // 07-04-09 Moved msc methods from G4VEmModel 65 // 07-04-09 Moved msc methods from G4VEmModel to G4VMscModel (VI)
66 // 13-10-10 Added G4VEmAngularDistribution (VI 66 // 13-10-10 Added G4VEmAngularDistribution (VI)
67 // 67 //
68 // Class Description: 68 // Class Description:
69 // 69 //
70 // Abstract interface to energy loss models 70 // Abstract interface to energy loss models
71 71
72 // ------------------------------------------- 72 // -------------------------------------------------------------------
73 // 73 //
74 74
75 #ifndef G4VEmModel_h 75 #ifndef G4VEmModel_h
76 #define G4VEmModel_h 1 76 #define G4VEmModel_h 1
77 77
78 #include "globals.hh" 78 #include "globals.hh"
79 #include "G4DynamicParticle.hh" 79 #include "G4DynamicParticle.hh"
80 #include "G4ParticleDefinition.hh" 80 #include "G4ParticleDefinition.hh"
81 #include "G4MaterialCutsCouple.hh" 81 #include "G4MaterialCutsCouple.hh"
82 #include "G4Material.hh" 82 #include "G4Material.hh"
83 #include "G4Element.hh" 83 #include "G4Element.hh"
84 #include "G4ElementVector.hh" 84 #include "G4ElementVector.hh"
85 #include "G4Isotope.hh" 85 #include "G4Isotope.hh"
86 #include "G4DataVector.hh" 86 #include "G4DataVector.hh"
87 #include "G4VEmFluctuationModel.hh" 87 #include "G4VEmFluctuationModel.hh"
88 #include "G4VEmAngularDistribution.hh" 88 #include "G4VEmAngularDistribution.hh"
89 #include "G4EmElementSelector.hh" 89 #include "G4EmElementSelector.hh"
90 #include <CLHEP/Random/RandomEngine.h> 90 #include <CLHEP/Random/RandomEngine.h>
91 #include <vector> 91 #include <vector>
92 92
93 class G4ElementData; 93 class G4ElementData;
94 class G4PhysicsTable; 94 class G4PhysicsTable;
95 class G4Region; 95 class G4Region;
96 class G4VParticleChange; 96 class G4VParticleChange;
97 class G4ParticleChangeForLoss; 97 class G4ParticleChangeForLoss;
98 class G4ParticleChangeForGamma; 98 class G4ParticleChangeForGamma;
99 class G4Track; 99 class G4Track;
100 class G4LossTableManager; 100 class G4LossTableManager;
101 101
102 class G4VEmModel 102 class G4VEmModel
103 { 103 {
104 104
105 public: 105 public:
106 106
107 explicit G4VEmModel(const G4String& nam); 107 explicit G4VEmModel(const G4String& nam);
108 108
109 virtual ~G4VEmModel(); 109 virtual ~G4VEmModel();
110 110
111 //------------------------------------------ 111 //------------------------------------------------------------------------
112 // Virtual methods to be implemented for any 112 // Virtual methods to be implemented for any concrete model
113 //------------------------------------------ 113 //------------------------------------------------------------------------
114 114
115 virtual void Initialise(const G4ParticleDefi 115 virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&) = 0;
116 116
117 virtual void SampleSecondaries(std::vector<G 117 virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
118 const G4Mater 118 const G4MaterialCutsCouple*,
119 const G4Dynam 119 const G4DynamicParticle*,
120 G4double tmin 120 G4double tmin = 0.0,
121 G4double tmax 121 G4double tmax = DBL_MAX) = 0;
122 122
123 //------------------------------------------ 123 //------------------------------------------------------------------------
124 // Methods for initialisation of MT; may be 124 // Methods for initialisation of MT; may be overwritten if needed
125 //------------------------------------------ 125 //------------------------------------------------------------------------
126 126
127 // initialisation in local thread << 127 // initilisation in local thread
128 virtual void InitialiseLocal(const G4Particl 128 virtual void InitialiseLocal(const G4ParticleDefinition*,
129 G4VEmModel* mas 129 G4VEmModel* masterModel);
130 130
131 // initialisation of a new material at run t << 131 // initilisation of a new material at run time
132 virtual void InitialiseForMaterial(const G4P 132 virtual void InitialiseForMaterial(const G4ParticleDefinition*,
133 const G4M 133 const G4Material*);
134 134
135 // initialisation of a new element at run ti << 135 // initilisation of a new element at run time
136 virtual void InitialiseForElement(const G4Pa 136 virtual void InitialiseForElement(const G4ParticleDefinition*,
137 G4int Z); 137 G4int Z);
138 138
139 //------------------------------------------ 139 //------------------------------------------------------------------------
140 // Methods with standard implementation; may 140 // Methods with standard implementation; may be overwritten if needed
141 //------------------------------------------ 141 //------------------------------------------------------------------------
142 142
143 // main method to compute dEdx 143 // main method to compute dEdx
144 virtual G4double ComputeDEDXPerVolume(const 144 virtual G4double ComputeDEDXPerVolume(const G4Material*,
145 const 145 const G4ParticleDefinition*,
146 G4doub 146 G4double kineticEnergy,
147 G4doub 147 G4double cutEnergy = DBL_MAX);
148 148
149 // main method to compute cross section per 149 // main method to compute cross section per Volume
150 virtual G4double CrossSectionPerVolume(const 150 virtual G4double CrossSectionPerVolume(const G4Material*,
151 const 151 const G4ParticleDefinition*,
152 G4dou 152 G4double kineticEnergy,
153 G4dou 153 G4double cutEnergy = 0.0,
154 G4dou 154 G4double maxEnergy = DBL_MAX);
155 155
156 // method to get partial cross section 156 // method to get partial cross section
157 virtual G4double GetPartialCrossSection(cons 157 virtual G4double GetPartialCrossSection(const G4Material*,
158 G4in 158 G4int level,
159 cons 159 const G4ParticleDefinition*,
160 G4do 160 G4double kineticEnergy);
161 161
162 // main method to compute cross section per 162 // main method to compute cross section per atom
163 virtual G4double ComputeCrossSectionPerAtom( 163 virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition*,
164 164 G4double kinEnergy,
165 165 G4double Z,
166 166 G4double A = 0., /* amu */
167 167 G4double cutEnergy = 0.0,
168 168 G4double maxEnergy = DBL_MAX);
169 169
170 // main method to compute cross section per 170 // main method to compute cross section per atomic shell
171 virtual G4double ComputeCrossSectionPerShell 171 virtual G4double ComputeCrossSectionPerShell(const G4ParticleDefinition*,
172 172 G4int Z, G4int shellIdx,
173 173 G4double kinEnergy,
174 174 G4double cutEnergy = 0.0,
175 175 G4double maxEnergy = DBL_MAX);
176 176
177 // Compute effective ion charge square 177 // Compute effective ion charge square
178 virtual G4double ChargeSquareRatio(const G4T 178 virtual G4double ChargeSquareRatio(const G4Track&);
179 179
180 // Compute effective ion charge square 180 // Compute effective ion charge square
181 virtual G4double GetChargeSquareRatio(const 181 virtual G4double GetChargeSquareRatio(const G4ParticleDefinition*,
182 const 182 const G4Material*,
183 G4doub 183 G4double kineticEnergy);
184 184
185 // Compute ion charge 185 // Compute ion charge
186 virtual G4double GetParticleCharge(const G4P 186 virtual G4double GetParticleCharge(const G4ParticleDefinition*,
187 const G4M 187 const G4Material*,
188 G4double 188 G4double kineticEnergy);
189 189
190 // Initialisation for a new track 190 // Initialisation for a new track
191 virtual void StartTracking(G4Track*); 191 virtual void StartTracking(G4Track*);
192 192
193 // add correction to energy loss and compute 193 // add correction to energy loss and compute non-ionizing energy loss
194 virtual void CorrectionsAlongStep(const G4Ma 194 virtual void CorrectionsAlongStep(const G4MaterialCutsCouple*,
195 const G4Dy 195 const G4DynamicParticle*,
196 const G4do << 196 G4double& eloss,
197 G4double& << 197 G4double& niel,
>> 198 G4double length);
198 199
199 // value which may be tabulated (by default 200 // value which may be tabulated (by default cross section)
200 virtual G4double Value(const G4MaterialCutsC 201 virtual G4double Value(const G4MaterialCutsCouple*,
201 const G4ParticleDefin 202 const G4ParticleDefinition*,
202 G4double kineticEnerg 203 G4double kineticEnergy);
203 204
204 // threshold for zero value 205 // threshold for zero value
205 virtual G4double MinPrimaryEnergy(const G4Ma 206 virtual G4double MinPrimaryEnergy(const G4Material*,
206 const G4Pa 207 const G4ParticleDefinition*,
207 G4double c 208 G4double cut = 0.0);
208 209
209 // model can define low-energy limit for the 210 // model can define low-energy limit for the cut
210 virtual G4double MinEnergyCut(const G4Partic 211 virtual G4double MinEnergyCut(const G4ParticleDefinition*,
211 const G4Materi 212 const G4MaterialCutsCouple*);
212 213
213 // initialisation at run time for a given ma << 214 // initilisation at run time for a given material
214 virtual void SetupForMaterial(const G4Partic 215 virtual void SetupForMaterial(const G4ParticleDefinition*,
215 const G4Materi 216 const G4Material*,
216 G4double kinet 217 G4double kineticEnergy);
217 218
218 // add a region for the model 219 // add a region for the model
219 virtual void DefineForRegion(const G4Region* 220 virtual void DefineForRegion(const G4Region*);
220 221
221 // fill number of different type of secondar <<
222 virtual void FillNumberOfSecondaries(G4int& <<
223 G4int& <<
224 <<
225 // for automatic documentation 222 // for automatic documentation
226 virtual void ModelDescription(std::ostream& 223 virtual void ModelDescription(std::ostream& outFile) const;
227 224
228 protected: 225 protected:
229 226
230 // initialisation of the ParticleChange for 227 // initialisation of the ParticleChange for the model
231 G4ParticleChangeForLoss* GetParticleChangeFo 228 G4ParticleChangeForLoss* GetParticleChangeForLoss();
232 229
233 // initialisation of the ParticleChange for 230 // initialisation of the ParticleChange for the model
234 G4ParticleChangeForGamma* GetParticleChangeF 231 G4ParticleChangeForGamma* GetParticleChangeForGamma();
235 232
236 // kinematically allowed max kinetic energy 233 // kinematically allowed max kinetic energy of a secondary
237 virtual G4double MaxSecondaryEnergy(const G4 234 virtual G4double MaxSecondaryEnergy(const G4ParticleDefinition*,
238 G4double 235 G4double kineticEnergy);
239 236
240 public: 237 public:
241 238
242 //------------------------------------------ 239 //------------------------------------------------------------------------
243 // Generic methods common to all models 240 // Generic methods common to all models
244 //------------------------------------------ 241 //------------------------------------------------------------------------
245 242
246 // should be called at initialisation to bui 243 // should be called at initialisation to build element selectors
247 void InitialiseElementSelectors(const G4Part 244 void InitialiseElementSelectors(const G4ParticleDefinition*,
248 const G4Data 245 const G4DataVector&);
249 246
250 // should be called at initialisation to acc 247 // should be called at initialisation to access element selectors
251 inline std::vector<G4EmElementSelector*>* Ge 248 inline std::vector<G4EmElementSelector*>* GetElementSelectors();
252 249
253 // should be called at initialisation to set 250 // should be called at initialisation to set element selectors
254 inline void SetElementSelectors(std::vector< 251 inline void SetElementSelectors(std::vector<G4EmElementSelector*>*);
255 252
256 // dEdx per unit length, base material appro 253 // dEdx per unit length, base material approach may be used
257 inline G4double ComputeDEDX( const G4Materia << 254 virtual inline G4double ComputeDEDX(const G4MaterialCutsCouple*,
258 const G4Particl 255 const G4ParticleDefinition*,
259 G4double kineti 256 G4double kineticEnergy,
260 G4double cutEne 257 G4double cutEnergy = DBL_MAX);
261 258
262 // cross section per volume, base material a 259 // cross section per volume, base material approach may be used
263 inline G4double CrossSection(const G4Materia 260 inline G4double CrossSection(const G4MaterialCutsCouple*,
264 const G4Particl 261 const G4ParticleDefinition*,
265 G4double kineti 262 G4double kineticEnergy,
266 G4double cutEne 263 G4double cutEnergy = 0.0,
267 G4double maxEne 264 G4double maxEnergy = DBL_MAX);
268 265
269 // compute mean free path via cross section 266 // compute mean free path via cross section per volume
270 inline G4double ComputeMeanFreePath(const G4 267 inline G4double ComputeMeanFreePath(const G4ParticleDefinition*,
271 G4double 268 G4double kineticEnergy,
272 const G4 269 const G4Material*,
273 G4double 270 G4double cutEnergy = 0.0,
274 G4double 271 G4double maxEnergy = DBL_MAX);
275 272
276 // generic cross section per element 273 // generic cross section per element
277 inline G4double ComputeCrossSectionPerAtom(c 274 inline G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition*,
278 c 275 const G4Element*,
279 G 276 G4double kinEnergy,
280 G 277 G4double cutEnergy = 0.0,
281 G 278 G4double maxEnergy = DBL_MAX);
282 279
283 // atom can be selected effitiantly if eleme 280 // atom can be selected effitiantly if element selectors are initialised
284 inline const G4Element* SelectRandomAtom(con 281 inline const G4Element* SelectRandomAtom(const G4MaterialCutsCouple*,
285 con 282 const G4ParticleDefinition*,
286 G4d 283 G4double kineticEnergy,
287 G4d 284 G4double cutEnergy = 0.0,
288 G4d 285 G4double maxEnergy = DBL_MAX);
289 // same as SelectRandomAtom above but more e 286 // same as SelectRandomAtom above but more efficient since log-ekin is known
290 inline const G4Element* SelectTargetAtom(con 287 inline const G4Element* SelectTargetAtom(const G4MaterialCutsCouple*,
291 con 288 const G4ParticleDefinition*,
292 G4d 289 G4double kineticEnergy,
293 G4d 290 G4double logKineticEnergy,
294 G4d 291 G4double cutEnergy = 0.0,
295 G4d 292 G4double maxEnergy = DBL_MAX);
296 293
>> 294
297 // to select atom cross section per volume i 295 // to select atom cross section per volume is recomputed for each element
298 const G4Element* SelectRandomAtom(const G4Ma 296 const G4Element* SelectRandomAtom(const G4Material*,
299 const G4Pa 297 const G4ParticleDefinition*,
300 G4double k 298 G4double kineticEnergy,
301 G4double c 299 G4double cutEnergy = 0.0,
302 G4double m 300 G4double maxEnergy = DBL_MAX);
303 301
304 // to select atom if cross section is propor 302 // to select atom if cross section is proportional number of electrons
305 const G4Element* GetCurrentElement(const G4M << 303 G4int SelectRandomAtomNumber(const G4Material*);
306 G4int SelectRandomAtomNumber(const G4Materia <<
307 304
308 // select isotope in order to have precise m 305 // select isotope in order to have precise mass of the nucleus
309 const G4Isotope* GetCurrentIsotope(const G4E << 306 G4int SelectIsotopeNumber(const G4Element*);
310 G4int SelectIsotopeNumber(const G4Element*) <<
311 307
312 //------------------------------------------ 308 //------------------------------------------------------------------------
313 // Get/Set methods 309 // Get/Set methods
314 //------------------------------------------ 310 //------------------------------------------------------------------------
315 311
316 void SetParticleChange(G4VParticleChange*, G 312 void SetParticleChange(G4VParticleChange*, G4VEmFluctuationModel* f=nullptr);
317 313
318 void SetCrossSectionTable(G4PhysicsTable*, G 314 void SetCrossSectionTable(G4PhysicsTable*, G4bool isLocal);
319 315
320 inline G4ElementData* GetElementData(); 316 inline G4ElementData* GetElementData();
321 317
322 inline G4PhysicsTable* GetCrossSectionTable( 318 inline G4PhysicsTable* GetCrossSectionTable();
323 319
324 inline G4VEmFluctuationModel* GetModelOfFluc 320 inline G4VEmFluctuationModel* GetModelOfFluctuations();
325 321
326 inline G4VEmAngularDistribution* GetAngularD 322 inline G4VEmAngularDistribution* GetAngularDistribution();
327 323
328 inline G4VEmModel* GetTripletModel(); 324 inline G4VEmModel* GetTripletModel();
329 325
330 inline void SetTripletModel(G4VEmModel*); 326 inline void SetTripletModel(G4VEmModel*);
331 327
332 inline void SetAngularDistribution(G4VEmAngu 328 inline void SetAngularDistribution(G4VEmAngularDistribution*);
333 329
334 inline G4double HighEnergyLimit() const; 330 inline G4double HighEnergyLimit() const;
335 331
336 inline G4double LowEnergyLimit() const; 332 inline G4double LowEnergyLimit() const;
337 333
338 inline G4double HighEnergyActivationLimit() 334 inline G4double HighEnergyActivationLimit() const;
339 335
340 inline G4double LowEnergyActivationLimit() c 336 inline G4double LowEnergyActivationLimit() const;
341 337
342 inline G4double PolarAngleLimit() const; 338 inline G4double PolarAngleLimit() const;
343 339
344 inline G4double SecondaryThreshold() const; 340 inline G4double SecondaryThreshold() const;
345 341
>> 342 inline G4bool LPMFlag() const;
>> 343
346 inline G4bool DeexcitationFlag() const; 344 inline G4bool DeexcitationFlag() const;
347 345
348 inline G4bool ForceBuildTableFlag() const; 346 inline G4bool ForceBuildTableFlag() const;
349 347
350 inline G4bool UseAngularGeneratorFlag() cons 348 inline G4bool UseAngularGeneratorFlag() const;
351 349
352 inline void SetAngularGeneratorFlag(G4bool); 350 inline void SetAngularGeneratorFlag(G4bool);
353 351
354 inline void SetHighEnergyLimit(G4double); 352 inline void SetHighEnergyLimit(G4double);
355 353
356 inline void SetLowEnergyLimit(G4double); 354 inline void SetLowEnergyLimit(G4double);
357 355
358 inline void SetActivationHighEnergyLimit(G4d 356 inline void SetActivationHighEnergyLimit(G4double);
359 357
360 inline void SetActivationLowEnergyLimit(G4do 358 inline void SetActivationLowEnergyLimit(G4double);
361 359
362 inline G4bool IsActive(G4double kinEnergy) c 360 inline G4bool IsActive(G4double kinEnergy) const;
363 361
364 inline void SetPolarAngleLimit(G4double); 362 inline void SetPolarAngleLimit(G4double);
365 363
366 inline void SetSecondaryThreshold(G4double); 364 inline void SetSecondaryThreshold(G4double);
367 365
>> 366 inline void SetLPMFlag(G4bool val);
>> 367
368 inline void SetDeexcitationFlag(G4bool val); 368 inline void SetDeexcitationFlag(G4bool val);
369 369
370 inline void SetForceBuildTable(G4bool val); 370 inline void SetForceBuildTable(G4bool val);
371 371
372 inline void SetFluctuationFlag(G4bool val); 372 inline void SetFluctuationFlag(G4bool val);
373 373
>> 374 inline void SetMasterThread(G4bool val);
>> 375
374 inline G4bool IsMaster() const; 376 inline G4bool IsMaster() const;
375 377
376 inline void SetUseBaseMaterials(G4bool val); 378 inline void SetUseBaseMaterials(G4bool val);
377 379
378 inline G4bool UseBaseMaterials() const; 380 inline G4bool UseBaseMaterials() const;
379 381
380 inline G4double MaxSecondaryKinEnergy(const 382 inline G4double MaxSecondaryKinEnergy(const G4DynamicParticle* dynParticle);
381 383
382 inline const G4String& GetName() const; 384 inline const G4String& GetName() const;
383 385
384 inline void SetCurrentCouple(const G4Materia 386 inline void SetCurrentCouple(const G4MaterialCutsCouple*);
385 387
386 inline G4bool IsLocked() const; << 388 inline const G4Element* GetCurrentElement() const;
387 389
388 inline void SetLocked(G4bool); << 390 inline const G4Isotope* GetCurrentIsotope() const;
389 391
390 // obsolete methods << 392 inline G4bool IsLocked() const;
391 [[deprecated("Use G4EmParameters::Instance() <<
392 void SetLPMFlag(G4bool); <<
393 <<
394 void SetMasterThread(G4bool); <<
395 393
396 // hide assignment operator << 394 inline void SetLocked(G4bool);
397 G4VEmModel & operator=(const G4VEmModel &ri <<
398 G4VEmModel(const G4VEmModel&) = delete; <<
399 395
400 protected: 396 protected:
401 397
402 inline const G4MaterialCutsCouple* CurrentCo 398 inline const G4MaterialCutsCouple* CurrentCouple() const;
403 399
404 inline void SetCurrentElement(const G4Elemen 400 inline void SetCurrentElement(const G4Element*);
405 401
406 private: 402 private:
407 403
>> 404 // hide assignment operator
>> 405 G4VEmModel & operator=(const G4VEmModel &right);
>> 406 G4VEmModel(const G4VEmModel&);
>> 407
408 // ======== Parameters of the class fixed at 408 // ======== Parameters of the class fixed at construction =========
409 409
410 G4VEmFluctuationModel* flucModel = null << 410 G4VEmFluctuationModel* flucModel;
411 G4VEmAngularDistribution* anglModel = null << 411 G4VEmAngularDistribution* anglModel;
412 G4VEmModel* fTripletModel = << 412 const G4String name;
413 const G4MaterialCutsCouple* fCurrentCouple = <<
414 const G4Element* fCurrentElement <<
415 std::vector<G4EmElementSelector*>* elmSelect <<
416 G4LossTableManager* fEmManager; <<
417 <<
418 protected: <<
419 <<
420 G4ElementData* fElementData = <<
421 G4VParticleChange* pParticleChange <<
422 G4PhysicsTable* xSectionTable = <<
423 const G4Material* pBaseMaterial = <<
424 const std::vector<G4double>* theDensityFacto <<
425 const std::vector<G4int>* theDensityIdx = <<
426 <<
427 G4double inveplus; <<
428 G4double pFactor = 1.0; <<
429 <<
430 private: <<
431 413
432 G4double lowLimit; << 414 // ======== Parameters of the class fixed at initialisation =======
433 G4double highLimit; <<
434 G4double eMinActive = 0.0; <<
435 G4double eMaxActive = DBL_MAX; <<
436 G4double secondaryThreshold = DBL_MAX; <<
437 G4double polarAngleLimit; <<
438 415
439 G4int nSelectors = 0; << 416 G4double lowLimit;
440 G4int nsec = 5; << 417 G4double highLimit;
>> 418 G4double eMinActive;
>> 419 G4double eMaxActive;
>> 420 G4double polarAngleLimit;
>> 421 G4double secondaryThreshold;
>> 422 G4bool theLPMflag;
>> 423 G4bool flagDeexcitation;
>> 424 G4bool flagForceBuildTable;
>> 425 G4bool isMaster;
>> 426
>> 427 G4bool localTable;
>> 428 G4bool localElmSelectors;
>> 429 G4bool useAngularGenerator;
>> 430 G4bool useBaseMaterials;
>> 431 G4bool isLocked;
>> 432 G4int nSelectors;
>> 433 std::vector<G4EmElementSelector*>* elmSelectors;
>> 434 G4LossTableManager* fEmManager;
441 435
442 protected: 436 protected:
443 437
444 std::size_t currentCoupleIndex = 0; << 438 G4ElementData* fElementData;
445 std::size_t basedCoupleIndex = 0; << 439 G4VParticleChange* pParticleChange;
446 G4bool lossFlucFlag = true; << 440 G4PhysicsTable* xSectionTable;
>> 441 const G4Material* pBaseMaterial;
>> 442 const std::vector<G4double>* theDensityFactor;
>> 443 const std::vector<G4int>* theDensityIdx;
>> 444 size_t idxTable;
>> 445 G4bool lossFlucFlag;
>> 446 G4double inveplus;
>> 447 G4double pFactor;
>> 448
>> 449 // ======== Cached values - may be state dependent ================
447 450
448 private: 451 private:
449 452
450 G4bool flagDeexcitation = false; << 453 const G4MaterialCutsCouple* fCurrentCouple;
451 G4bool flagForceBuildTable = false; << 454 const G4Element* fCurrentElement;
452 G4bool isMaster = true; << 455 const G4Isotope* fCurrentIsotope;
453 << 456 G4VEmModel* fTripletModel;
454 G4bool localTable = true; <<
455 G4bool localElmSelectors = true; <<
456 G4bool useAngularGenerator = false; <<
457 G4bool useBaseMaterials = false; <<
458 G4bool isLocked = false; <<
459 457
460 const G4String name; << 458 G4int nsec;
461 std::vector<G4double> xsec; 459 std::vector<G4double> xsec;
462 460
463 }; 461 };
464 462
465 // ======== Run time inline methods ========== 463 // ======== Run time inline methods ================
466 464
467 inline void G4VEmModel::SetCurrentCouple(const 465 inline void G4VEmModel::SetCurrentCouple(const G4MaterialCutsCouple* ptr)
468 { 466 {
469 if(fCurrentCouple != ptr) { 467 if(fCurrentCouple != ptr) {
470 fCurrentCouple = ptr; 468 fCurrentCouple = ptr;
471 basedCoupleIndex = currentCoupleIndex = pt <<
472 pBaseMaterial = ptr->GetMaterial(); 469 pBaseMaterial = ptr->GetMaterial();
473 pFactor = 1.0; 470 pFactor = 1.0;
474 if(useBaseMaterials) { << 471 if(useBaseMaterials && pBaseMaterial->GetBaseMaterial()) {
475 basedCoupleIndex = (*theDensityIdx)[curr << 472 pBaseMaterial = pBaseMaterial->GetBaseMaterial();
476 if(nullptr != pBaseMaterial->GetBaseMate << 473 pFactor = (*theDensityFactor)[(*theDensityIdx)[ptr->GetIndex()]];
477 pBaseMaterial = pBaseMaterial->GetBaseMateri <<
478 pFactor = (*theDensityFactor)[currentCou <<
479 } 474 }
480 } 475 }
481 } 476 }
482 477
483 //....oooOO0OOooo........oooOO0OOooo........oo 478 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
484 479
485 inline const G4MaterialCutsCouple* G4VEmModel: 480 inline const G4MaterialCutsCouple* G4VEmModel::CurrentCouple() const
486 { 481 {
487 return fCurrentCouple; 482 return fCurrentCouple;
488 } 483 }
489 484
490 //....oooOO0OOooo........oooOO0OOooo........oo 485 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
491 486
492 inline void G4VEmModel::SetCurrentElement(cons 487 inline void G4VEmModel::SetCurrentElement(const G4Element* elm)
493 { 488 {
494 fCurrentElement = elm; 489 fCurrentElement = elm;
>> 490 fCurrentIsotope = nullptr;
>> 491 }
>> 492
>> 493 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
>> 494
>> 495 inline const G4Element* G4VEmModel::GetCurrentElement() const
>> 496 {
>> 497 return fCurrentElement;
>> 498 }
>> 499
>> 500 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
>> 501
>> 502 inline const G4Isotope* G4VEmModel::GetCurrentIsotope() const
>> 503 {
>> 504 return fCurrentIsotope;
495 } 505 }
496 506
497 //....oooOO0OOooo........oooOO0OOooo........oo 507 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
498 508
499 inline 509 inline
500 G4double G4VEmModel::MaxSecondaryKinEnergy(con 510 G4double G4VEmModel::MaxSecondaryKinEnergy(const G4DynamicParticle* dynPart)
501 { 511 {
502 return MaxSecondaryEnergy(dynPart->GetPartic 512 return MaxSecondaryEnergy(dynPart->GetParticleDefinition(),
503 dynPart->GetKineti 513 dynPart->GetKineticEnergy());
504 } 514 }
505 515
506 //....oooOO0OOooo........oooOO0OOooo........oo 516 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
507 517
508 inline G4double G4VEmModel::ComputeDEDX(const 518 inline G4double G4VEmModel::ComputeDEDX(const G4MaterialCutsCouple* couple,
509 const 519 const G4ParticleDefinition* part,
510 G4doub 520 G4double kinEnergy,
511 G4doub 521 G4double cutEnergy)
512 { 522 {
513 SetCurrentCouple(couple); 523 SetCurrentCouple(couple);
514 return pFactor*ComputeDEDXPerVolume(pBaseMat 524 return pFactor*ComputeDEDXPerVolume(pBaseMaterial,part,kinEnergy,cutEnergy);
515 } 525 }
516 526
517 //....oooOO0OOooo........oooOO0OOooo........oo 527 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
518 528
519 inline G4double G4VEmModel::CrossSection(const 529 inline G4double G4VEmModel::CrossSection(const G4MaterialCutsCouple* couple,
520 const 530 const G4ParticleDefinition* part,
521 G4dou 531 G4double kinEnergy,
522 G4dou 532 G4double cutEnergy,
523 G4dou 533 G4double maxEnergy)
524 { 534 {
525 SetCurrentCouple(couple); 535 SetCurrentCouple(couple);
526 return pFactor*CrossSectionPerVolume(pBaseMa 536 return pFactor*CrossSectionPerVolume(pBaseMaterial,part,kinEnergy,
527 cutEner 537 cutEnergy,maxEnergy);
528 } 538 }
529 539
530 //....oooOO0OOooo........oooOO0OOooo........oo 540 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
531 541
532 inline 542 inline
533 G4double G4VEmModel::ComputeMeanFreePath(const 543 G4double G4VEmModel::ComputeMeanFreePath(const G4ParticleDefinition* part,
534 G4dou 544 G4double ekin,
535 const 545 const G4Material* material,
536 G4dou 546 G4double emin,
537 G4dou 547 G4double emax)
538 { 548 {
539 G4double cross = CrossSectionPerVolume(mater 549 G4double cross = CrossSectionPerVolume(material,part,ekin,emin,emax);
540 return (cross > 0.0) ? 1./cross : DBL_MAX; 550 return (cross > 0.0) ? 1./cross : DBL_MAX;
541 } 551 }
542 552
543 //....oooOO0OOooo........oooOO0OOooo........oo 553 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
544 554
545 inline G4double 555 inline G4double
546 G4VEmModel::ComputeCrossSectionPerAtom(const G 556 G4VEmModel::ComputeCrossSectionPerAtom(const G4ParticleDefinition* part,
547 const G 557 const G4Element* elm,
548 G4doubl 558 G4double kinEnergy,
549 G4doubl 559 G4double cutEnergy,
550 G4doubl 560 G4double maxEnergy)
551 { 561 {
552 fCurrentElement = elm; << 562 SetCurrentElement(elm);
553 return ComputeCrossSectionPerAtom(part,kinEn 563 return ComputeCrossSectionPerAtom(part,kinEnergy,elm->GetZ(),elm->GetN(),
554 cutEnergy, 564 cutEnergy,maxEnergy);
555 } 565 }
556 566
557 //....oooOO0OOooo........oooOO0OOooo........oo 567 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
558 568
559 inline const G4Element* 569 inline const G4Element*
560 G4VEmModel::SelectRandomAtom(const G4MaterialC 570 G4VEmModel::SelectRandomAtom(const G4MaterialCutsCouple* couple,
561 const G4ParticleD 571 const G4ParticleDefinition* part,
562 G4double kinEnerg 572 G4double kinEnergy,
563 G4double cutEnerg 573 G4double cutEnergy,
564 G4double maxEnerg 574 G4double maxEnergy)
565 { 575 {
566 SetCurrentCouple(couple); 576 SetCurrentCouple(couple);
567 fCurrentElement = (nSelectors > 0) ? 577 fCurrentElement = (nSelectors > 0) ?
568 ((*elmSelectors)[couple->GetIndex()])->Sel 578 ((*elmSelectors)[couple->GetIndex()])->SelectRandomAtom(kinEnergy) :
569 SelectRandomAtom(pBaseMaterial,part,kinEne 579 SelectRandomAtom(pBaseMaterial,part,kinEnergy,cutEnergy,maxEnergy);
>> 580 fCurrentIsotope = nullptr;
570 return fCurrentElement; 581 return fCurrentElement;
571 } 582 }
572 583
573 //....oooOO0OOooo........oooOO0OOooo........oo 584 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
574 585
575 inline const G4Element* 586 inline const G4Element*
576 G4VEmModel::SelectTargetAtom(const G4MaterialC 587 G4VEmModel::SelectTargetAtom(const G4MaterialCutsCouple* couple,
577 const G4ParticleD 588 const G4ParticleDefinition* part,
578 G4double kinEnerg 589 G4double kinEnergy,
579 G4double logKinE, 590 G4double logKinE,
580 G4double cutEnerg 591 G4double cutEnergy,
581 G4double maxEnerg 592 G4double maxEnergy)
582 { 593 {
583 SetCurrentCouple(couple); 594 SetCurrentCouple(couple);
584 fCurrentElement = (nSelectors > 0) 595 fCurrentElement = (nSelectors > 0)
585 ? ((*elmSelectors)[couple->GetIndex()])->Se 596 ? ((*elmSelectors)[couple->GetIndex()])->SelectRandomAtom(kinEnergy,logKinE)
586 : SelectRandomAtom(pBaseMaterial,part,kinEn 597 : SelectRandomAtom(pBaseMaterial,part,kinEnergy,cutEnergy,maxEnergy);
>> 598 fCurrentIsotope = nullptr;
587 return fCurrentElement; 599 return fCurrentElement;
588 } 600 }
589 601
590 // ======== Get/Set inline methods used at ini 602 // ======== Get/Set inline methods used at initialisation ================
591 603
592 inline G4VEmFluctuationModel* G4VEmModel::GetM 604 inline G4VEmFluctuationModel* G4VEmModel::GetModelOfFluctuations()
593 { 605 {
594 return flucModel; 606 return flucModel;
595 } 607 }
596 608
597 //....oooOO0OOooo........oooOO0OOooo........oo 609 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
598 610
599 inline G4VEmAngularDistribution* G4VEmModel::G 611 inline G4VEmAngularDistribution* G4VEmModel::GetAngularDistribution()
600 { 612 {
601 return anglModel; 613 return anglModel;
602 } 614 }
603 615
604 //....oooOO0OOooo........oooOO0OOooo........oo 616 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
605 617
606 inline void G4VEmModel::SetAngularDistribution 618 inline void G4VEmModel::SetAngularDistribution(G4VEmAngularDistribution* p)
607 { 619 {
608 if(p != anglModel) { 620 if(p != anglModel) {
609 delete anglModel; 621 delete anglModel;
610 anglModel = p; 622 anglModel = p;
611 } 623 }
612 } 624 }
613 625
614 //....oooOO0OOooo........oooOO0OOooo........oo 626 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
615 627
616 inline G4VEmModel* G4VEmModel::GetTripletModel 628 inline G4VEmModel* G4VEmModel::GetTripletModel()
617 { 629 {
618 return fTripletModel; 630 return fTripletModel;
619 } 631 }
620 632
621 //....oooOO0OOooo........oooOO0OOooo........oo 633 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
622 634
623 inline void G4VEmModel::SetTripletModel(G4VEmM 635 inline void G4VEmModel::SetTripletModel(G4VEmModel* p)
624 { 636 {
625 if(p != fTripletModel) { 637 if(p != fTripletModel) {
626 delete fTripletModel; 638 delete fTripletModel;
627 fTripletModel = p; 639 fTripletModel = p;
628 } 640 }
629 } 641 }
630 642
631 //....oooOO0OOooo........oooOO0OOooo........oo 643 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
632 644
633 inline G4double G4VEmModel::HighEnergyLimit() 645 inline G4double G4VEmModel::HighEnergyLimit() const
634 { 646 {
635 return highLimit; 647 return highLimit;
636 } 648 }
637 649
638 //....oooOO0OOooo........oooOO0OOooo........oo 650 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
639 651
640 inline G4double G4VEmModel::LowEnergyLimit() c 652 inline G4double G4VEmModel::LowEnergyLimit() const
641 { 653 {
642 return lowLimit; 654 return lowLimit;
643 } 655 }
644 656
645 //....oooOO0OOooo........oooOO0OOooo........oo 657 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
646 658
647 inline G4double G4VEmModel::HighEnergyActivati 659 inline G4double G4VEmModel::HighEnergyActivationLimit() const
648 { 660 {
649 return eMaxActive; 661 return eMaxActive;
650 } 662 }
651 663
652 //....oooOO0OOooo........oooOO0OOooo........oo 664 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
653 665
654 inline G4double G4VEmModel::LowEnergyActivatio 666 inline G4double G4VEmModel::LowEnergyActivationLimit() const
655 { 667 {
656 return eMinActive; 668 return eMinActive;
657 } 669 }
658 670
659 //....oooOO0OOooo........oooOO0OOooo........oo 671 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
660 672
661 inline G4double G4VEmModel::PolarAngleLimit() 673 inline G4double G4VEmModel::PolarAngleLimit() const
662 { 674 {
663 return polarAngleLimit; 675 return polarAngleLimit;
664 } 676 }
665 677
666 //....oooOO0OOooo........oooOO0OOooo........oo 678 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
667 679
668 inline G4double G4VEmModel::SecondaryThreshold 680 inline G4double G4VEmModel::SecondaryThreshold() const
669 { 681 {
670 return secondaryThreshold; 682 return secondaryThreshold;
671 } 683 }
672 684
673 //....oooOO0OOooo........oooOO0OOooo........oo 685 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
674 686
>> 687 inline G4bool G4VEmModel::LPMFlag() const
>> 688 {
>> 689 return theLPMflag;
>> 690 }
>> 691
>> 692 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
>> 693
675 inline G4bool G4VEmModel::DeexcitationFlag() c 694 inline G4bool G4VEmModel::DeexcitationFlag() const
676 { 695 {
677 return flagDeexcitation; 696 return flagDeexcitation;
678 } 697 }
679 698
680 //....oooOO0OOooo........oooOO0OOooo........oo 699 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
681 700
682 inline G4bool G4VEmModel::ForceBuildTableFlag( 701 inline G4bool G4VEmModel::ForceBuildTableFlag() const
683 { 702 {
684 return flagForceBuildTable; 703 return flagForceBuildTable;
685 } 704 }
686 705
687 //....oooOO0OOooo........oooOO0OOooo........oo 706 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
688 707
689 inline G4bool G4VEmModel::UseAngularGeneratorF 708 inline G4bool G4VEmModel::UseAngularGeneratorFlag() const
690 { 709 {
691 return useAngularGenerator; 710 return useAngularGenerator;
692 } 711 }
693 712
694 //....oooOO0OOooo........oooOO0OOooo........oo 713 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
695 714
696 inline void G4VEmModel::SetAngularGeneratorFla 715 inline void G4VEmModel::SetAngularGeneratorFlag(G4bool val)
697 { 716 {
698 useAngularGenerator = val; 717 useAngularGenerator = val;
699 } 718 }
700 719
701 //....oooOO0OOooo........oooOO0OOooo........oo 720 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
702 721
703 inline void G4VEmModel::SetFluctuationFlag(G4b 722 inline void G4VEmModel::SetFluctuationFlag(G4bool val)
704 { 723 {
705 lossFlucFlag = val; 724 lossFlucFlag = val;
706 } 725 }
707 726
708 //....oooOO0OOooo........oooOO0OOooo........oo 727 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
709 728
>> 729 inline void G4VEmModel::SetMasterThread(G4bool val)
>> 730 {
>> 731 isMaster = val;
>> 732 }
>> 733
>> 734 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
>> 735
710 inline G4bool G4VEmModel::IsMaster() const 736 inline G4bool G4VEmModel::IsMaster() const
711 { 737 {
712 return isMaster; 738 return isMaster;
713 } 739 }
714 740
715 //....oooOO0OOooo........oooOO0OOooo........oo 741 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
716 742
717 inline void G4VEmModel::SetUseBaseMaterials(G4 743 inline void G4VEmModel::SetUseBaseMaterials(G4bool val)
718 { 744 {
719 useBaseMaterials = val; 745 useBaseMaterials = val;
720 } 746 }
721 747
722 //....oooOO0OOooo........oooOO0OOooo........oo 748 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
723 749
724 inline G4bool G4VEmModel::UseBaseMaterials() c 750 inline G4bool G4VEmModel::UseBaseMaterials() const
725 { 751 {
726 return useBaseMaterials; 752 return useBaseMaterials;
727 } 753 }
728 754
729 //....oooOO0OOooo........oooOO0OOooo........oo 755 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
730 756
731 inline void G4VEmModel::SetHighEnergyLimit(G4d 757 inline void G4VEmModel::SetHighEnergyLimit(G4double val)
732 { 758 {
733 highLimit = val; 759 highLimit = val;
734 } 760 }
735 761
736 //....oooOO0OOooo........oooOO0OOooo........oo 762 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
737 763
738 inline void G4VEmModel::SetLowEnergyLimit(G4do 764 inline void G4VEmModel::SetLowEnergyLimit(G4double val)
739 { 765 {
740 lowLimit = val; 766 lowLimit = val;
741 } 767 }
742 768
743 //....oooOO0OOooo........oooOO0OOooo........oo 769 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
744 770
745 inline void G4VEmModel::SetActivationHighEnerg 771 inline void G4VEmModel::SetActivationHighEnergyLimit(G4double val)
746 { 772 {
747 eMaxActive = val; 773 eMaxActive = val;
748 } 774 }
749 775
750 //....oooOO0OOooo........oooOO0OOooo........oo 776 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
751 777
752 inline void G4VEmModel::SetActivationLowEnergy 778 inline void G4VEmModel::SetActivationLowEnergyLimit(G4double val)
753 { 779 {
754 eMinActive = val; 780 eMinActive = val;
755 } 781 }
756 782
757 //....oooOO0OOooo........oooOO0OOooo........oo 783 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
758 784
759 inline G4bool G4VEmModel::IsActive(G4double ki 785 inline G4bool G4VEmModel::IsActive(G4double kinEnergy) const
760 { 786 {
761 return (kinEnergy >= eMinActive && kinEnergy 787 return (kinEnergy >= eMinActive && kinEnergy <= eMaxActive);
762 } 788 }
763 789
764 //....oooOO0OOooo........oooOO0OOooo........oo 790 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
765 791
766 inline void G4VEmModel::SetPolarAngleLimit(G4d 792 inline void G4VEmModel::SetPolarAngleLimit(G4double val)
767 { 793 {
768 if(!isLocked) { polarAngleLimit = val; } 794 if(!isLocked) { polarAngleLimit = val; }
769 } 795 }
770 796
771 //....oooOO0OOooo........oooOO0OOooo........oo 797 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
772 798
773 inline void G4VEmModel::SetSecondaryThreshold( 799 inline void G4VEmModel::SetSecondaryThreshold(G4double val)
774 { 800 {
775 secondaryThreshold = val; 801 secondaryThreshold = val;
776 } 802 }
777 803
778 //....oooOO0OOooo........oooOO0OOooo........oo 804 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
779 805
>> 806 inline void G4VEmModel::SetLPMFlag(G4bool val)
>> 807 {
>> 808 theLPMflag = val;
>> 809 }
>> 810
>> 811 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
>> 812
780 inline void G4VEmModel::SetDeexcitationFlag(G4 813 inline void G4VEmModel::SetDeexcitationFlag(G4bool val)
781 { 814 {
782 flagDeexcitation = val; 815 flagDeexcitation = val;
783 } 816 }
784 817
785 //....oooOO0OOooo........oooOO0OOooo........oo 818 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
786 819
787 inline void G4VEmModel::SetForceBuildTable(G4b 820 inline void G4VEmModel::SetForceBuildTable(G4bool val)
788 { 821 {
789 flagForceBuildTable = val; 822 flagForceBuildTable = val;
790 } 823 }
791 824
792 //....oooOO0OOooo........oooOO0OOooo........oo 825 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
793 826
794 inline const G4String& G4VEmModel::GetName() c 827 inline const G4String& G4VEmModel::GetName() const
795 { 828 {
796 return name; 829 return name;
797 } 830 }
798 831
799 //....oooOO0OOooo........oooOO0OOooo........oo 832 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
800 833
801 inline std::vector<G4EmElementSelector*>* G4VE 834 inline std::vector<G4EmElementSelector*>* G4VEmModel::GetElementSelectors()
802 { 835 {
803 return elmSelectors; 836 return elmSelectors;
804 } 837 }
805 838
806 //....oooOO0OOooo........oooOO0OOooo........oo 839 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
807 840
808 inline void 841 inline void
809 G4VEmModel::SetElementSelectors(std::vector<G4 842 G4VEmModel::SetElementSelectors(std::vector<G4EmElementSelector*>* p)
810 { 843 {
811 if(p != elmSelectors) { 844 if(p != elmSelectors) {
812 elmSelectors = p; 845 elmSelectors = p;
813 nSelectors = (nullptr != elmSelectors) ? G << 846 nSelectors = (elmSelectors) ? G4int(elmSelectors->size()) : 0;
814 localElmSelectors = false; 847 localElmSelectors = false;
815 } 848 }
816 } 849 }
817 850
818 //....oooOO0OOooo........oooOO0OOooo........oo 851 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
819 852
820 inline G4ElementData* G4VEmModel::GetElementDa 853 inline G4ElementData* G4VEmModel::GetElementData()
821 { 854 {
822 return fElementData; 855 return fElementData;
823 } 856 }
824 857
825 //....oooOO0OOooo........oooOO0OOooo........oo 858 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
826 859
827 inline G4PhysicsTable* G4VEmModel::GetCrossSec 860 inline G4PhysicsTable* G4VEmModel::GetCrossSectionTable()
828 { 861 {
829 return xSectionTable; 862 return xSectionTable;
830 } 863 }
831 864
832 //....oooOO0OOooo........oooOO0OOooo........oo 865 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
833 866
834 inline G4bool G4VEmModel::IsLocked() const 867 inline G4bool G4VEmModel::IsLocked() const
835 { 868 {
836 return isLocked; 869 return isLocked;
837 } 870 }
838 871
839 //....oooOO0OOooo........oooOO0OOooo........oo 872 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
840 873
841 inline void G4VEmModel::SetLocked(G4bool val) 874 inline void G4VEmModel::SetLocked(G4bool val)
842 { 875 {
843 isLocked = val; 876 isLocked = val;
844 } 877 }
845 878
846 //....oooOO0OOooo........oooOO0OOooo........oo 879 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
847 880
848 #endif 881 #endif
849 882