Geant4  10.03.p01
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
G4VEmModel.hh
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // $Id: G4VEmModel.hh 97742 2016-06-08 09:24:54Z gcosmo $
27 //
28 // -------------------------------------------------------------------
29 //
30 // GEANT4 Class header file
31 //
32 //
33 // File name: G4VEmModel
34 //
35 // Author: Vladimir Ivanchenko
36 //
37 // Creation date: 03.01.2002
38 //
39 // Modifications:
40 //
41 // 23-12-02 V.Ivanchenko change interface before move to cut per region
42 // 24-01-03 Cut per region (V.Ivanchenko)
43 // 13-02-03 Add name (V.Ivanchenko)
44 // 25-02-03 Add sample theta and displacement (V.Ivanchenko)
45 // 23-07-03 Replace G4Material by G4MaterialCutCouple in dE/dx and CrossSection
46 // calculation (V.Ivanchenko)
47 // 01-03-04 L.Urban signature changed in SampleCosineTheta
48 // 23-04-04 L.urban signature of SampleCosineTheta changed back
49 // 17-11-04 Add method CrossSectionPerAtom (V.Ivanchenko)
50 // 14-03-05 Reduce number of pure virtual methods and make inline part
51 // separate (V.Ivanchenko)
52 // 24-03-05 Remove IsInCharge and add G4VParticleChange in the constructor (VI)
53 // 08-04-05 Major optimisation of internal interfaces (V.Ivantchenko)
54 // 15-04-05 optimize internal interface for msc (V.Ivanchenko)
55 // 08-05-05 A -> N (V.Ivanchenko)
56 // 25-07-05 Move constructor and destructor to the body (V.Ivanchenko)
57 // 02-02-06 ComputeCrossSectionPerAtom: default value A=0. (mma)
58 // 06-02-06 add method ComputeMeanFreePath() (mma)
59 // 07-03-06 Optimize msc methods (V.Ivanchenko)
60 // 29-06-06 Add member currentElement and Get/Set methods (V.Ivanchenko)
61 // 29-10-07 Added SampleScattering (V.Ivanchenko)
62 // 15-07-08 Reorder class members and improve comments (VI)
63 // 21-07-08 Added vector of G4ElementSelector and methods to use it (VI)
64 // 12-09-08 Added methods GetParticleCharge, GetChargeSquareRatio,
65 // CorrectionsAlongStep, ActivateNuclearStopping (VI)
66 // 16-02-09 Moved implementations of virtual methods to source (VI)
67 // 07-04-09 Moved msc methods from G4VEmModel to G4VMscModel (VI)
68 // 13-10-10 Added G4VEmAngularDistribution (VI)
69 //
70 // Class Description:
71 //
72 // Abstract interface to energy loss models
73 
74 // -------------------------------------------------------------------
75 //
76 
77 #ifndef G4VEmModel_h
78 #define G4VEmModel_h 1
79 
80 #include "globals.hh"
81 #include "G4DynamicParticle.hh"
82 #include "G4ParticleDefinition.hh"
83 #include "G4MaterialCutsCouple.hh"
84 #include "G4Material.hh"
85 #include "G4Element.hh"
86 #include "G4ElementVector.hh"
87 #include "G4Isotope.hh"
88 #include "G4DataVector.hh"
89 #include "G4VEmFluctuationModel.hh"
90 #include "G4VEmAngularDistribution.hh"
91 #include "G4EmElementSelector.hh"
92 #include <CLHEP/Random/RandomEngine.h>
93 #include <vector>
94 
95 class G4ElementData;
96 class G4PhysicsTable;
97 class G4Region;
98 class G4VParticleChange;
99 class G4ParticleChangeForLoss;
100 class G4ParticleChangeForGamma;
101 class G4Track;
102 class G4LossTableManager;
103 
104 class G4VEmModel
105 {
106 
107 public:
108 
109  explicit G4VEmModel(const G4String& nam);
110 
111  virtual ~G4VEmModel();
112 
113  //------------------------------------------------------------------------
114  // Virtual methods to be implemented for any concrete model
115  //------------------------------------------------------------------------
116 
117  virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&) = 0;
118 
119  virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
120  const G4MaterialCutsCouple*,
121  const G4DynamicParticle*,
122  G4double tmin = 0.0,
123  G4double tmax = DBL_MAX) = 0;
124 
125  //------------------------------------------------------------------------
126  // Methods for initialisation of MT; may be overwritten if needed
127  //------------------------------------------------------------------------
128 
129  // initilisation in local thread
130  virtual void InitialiseLocal(const G4ParticleDefinition*,
131  G4VEmModel* masterModel);
132 
133  // initilisation of a new material at run time
134  virtual void InitialiseForMaterial(const G4ParticleDefinition*,
135  const G4Material*);
136 
137  // initilisation of a new element at run time
138  virtual void InitialiseForElement(const G4ParticleDefinition*,
139  G4int Z);
140 
141  //------------------------------------------------------------------------
142  // Methods with standard implementation; may be overwritten if needed
143  //------------------------------------------------------------------------
144 
145  // main method to compute dEdx
146  virtual G4double ComputeDEDXPerVolume(const G4Material*,
147  const G4ParticleDefinition*,
148  G4double kineticEnergy,
149  G4double cutEnergy = DBL_MAX);
150 
151  // main method to compute cross section per Volume
152  virtual G4double CrossSectionPerVolume(const G4Material*,
153  const G4ParticleDefinition*,
154  G4double kineticEnergy,
155  G4double cutEnergy = 0.0,
156  G4double maxEnergy = DBL_MAX);
157 
158  // method to get partial cross section
159  virtual G4double GetPartialCrossSection(const G4Material*,
160  G4int level,
161  const G4ParticleDefinition*,
162  G4double kineticEnergy);
163 
164  // main method to compute cross section per atom
165  virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition*,
166  G4double kinEnergy,
167  G4double Z,
168  G4double A = 0., /* amu */
169  G4double cutEnergy = 0.0,
170  G4double maxEnergy = DBL_MAX);
171 
172  // main method to compute cross section per atomic shell
173  virtual G4double ComputeCrossSectionPerShell(const G4ParticleDefinition*,
174  G4int Z, G4int shellIdx,
175  G4double kinEnergy,
176  G4double cutEnergy = 0.0,
177  G4double maxEnergy = DBL_MAX);
178 
179  // Compute effective ion charge square
180  virtual G4double ChargeSquareRatio(const G4Track&);
181 
182  // Compute effective ion charge square
183  virtual G4double GetChargeSquareRatio(const G4ParticleDefinition*,
184  const G4Material*,
185  G4double kineticEnergy);
186 
187  // Compute ion charge
188  virtual G4double GetParticleCharge(const G4ParticleDefinition*,
189  const G4Material*,
190  G4double kineticEnergy);
191 
192  // Initialisation for a new track
193  virtual void StartTracking(G4Track*);
194 
195  // add correction to energy loss and compute non-ionizing energy loss
196  virtual void CorrectionsAlongStep(const G4MaterialCutsCouple*,
197  const G4DynamicParticle*,
198  G4double& eloss,
199  G4double& niel,
200  G4double length);
201 
202  // value which may be tabulated (by default cross section)
203  virtual G4double Value(const G4MaterialCutsCouple*,
204  const G4ParticleDefinition*,
205  G4double kineticEnergy);
206 
207  // threshold for zero value
208  virtual G4double MinPrimaryEnergy(const G4Material*,
209  const G4ParticleDefinition*,
210  G4double cut = 0.0);
211 
212  // model can define low-energy limit for the cut
213  virtual G4double MinEnergyCut(const G4ParticleDefinition*,
214  const G4MaterialCutsCouple*);
215 
216  // initilisation at run time for a given material
217  virtual void SetupForMaterial(const G4ParticleDefinition*,
218  const G4Material*,
219  G4double kineticEnergy);
220 
221  // add a region for the model
222  virtual void DefineForRegion(const G4Region*);
223 
224  // for automatic documentation
225  virtual void ModelDescription(std::ostream& outFile) const;
226 
227 protected:
228 
229  // initialisation of the ParticleChange for the model
230  G4ParticleChangeForLoss* GetParticleChangeForLoss();
231 
232  // initialisation of the ParticleChange for the model
233  G4ParticleChangeForGamma* GetParticleChangeForGamma();
234 
235  // kinematically allowed max kinetic energy of a secondary
236  virtual G4double MaxSecondaryEnergy(const G4ParticleDefinition*,
237  G4double kineticEnergy);
238 
239 public:
240 
241  //------------------------------------------------------------------------
242  // Generic methods common to all models
243  //------------------------------------------------------------------------
244 
245  // should be called at initialisation to build element selectors
246  void InitialiseElementSelectors(const G4ParticleDefinition*,
247  const G4DataVector&);
248 
249  // should be called at initialisation to access element selectors
250  inline std::vector<G4EmElementSelector*>* GetElementSelectors();
251 
252  // should be called at initialisation to set element selectors
253  inline void SetElementSelectors(std::vector<G4EmElementSelector*>*);
254 
255  // dEdx per unit length
256  virtual inline G4double ComputeDEDX(const G4MaterialCutsCouple*,
257  const G4ParticleDefinition*,
258  G4double kineticEnergy,
259  G4double cutEnergy = DBL_MAX);
260 
261  // cross section per volume
262  inline G4double CrossSection(const G4MaterialCutsCouple*,
263  const G4ParticleDefinition*,
264  G4double kineticEnergy,
265  G4double cutEnergy = 0.0,
266  G4double maxEnergy = DBL_MAX);
267 
268  // compute mean free path via cross section per volume
269  inline G4double ComputeMeanFreePath(const G4ParticleDefinition*,
270  G4double kineticEnergy,
271  const G4Material*,
272  G4double cutEnergy = 0.0,
273  G4double maxEnergy = DBL_MAX);
274 
275  // generic cross section per element
276  inline G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition*,
277  const G4Element*,
278  G4double kinEnergy,
279  G4double cutEnergy = 0.0,
280  G4double maxEnergy = DBL_MAX);
281 
282  // select isotope in order to have precise mass of the nucleus
283  inline G4int SelectIsotopeNumber(const G4Element*);
284 
285  // atom can be selected effitiantly if element selectors are initialised
286  inline const G4Element* SelectRandomAtom(const G4MaterialCutsCouple*,
287  const G4ParticleDefinition*,
288  G4double kineticEnergy,
289  G4double cutEnergy = 0.0,
290  G4double maxEnergy = DBL_MAX);
291 
292  // to select atom cross section per volume is recomputed for each element
293  const G4Element* SelectRandomAtom(const G4Material*,
294  const G4ParticleDefinition*,
295  G4double kineticEnergy,
296  G4double cutEnergy = 0.0,
297  G4double maxEnergy = DBL_MAX);
298 
299  // to select atom if cross section is proportional number of electrons
300  inline G4int SelectRandomAtomNumber(const G4Material*);
301 
302  //------------------------------------------------------------------------
303  // Get/Set methods
304  //------------------------------------------------------------------------
305 
306  void SetParticleChange(G4VParticleChange*, G4VEmFluctuationModel* f=nullptr);
307 
308  void SetCrossSectionTable(G4PhysicsTable*, G4bool isLocal);
309 
310  inline G4ElementData* GetElementData();
311 
312  inline G4PhysicsTable* GetCrossSectionTable();
313 
314  inline G4VEmFluctuationModel* GetModelOfFluctuations();
315 
316  inline G4VEmAngularDistribution* GetAngularDistribution();
317 
318  inline void SetAngularDistribution(G4VEmAngularDistribution*);
319 
320  inline G4double HighEnergyLimit() const;
321 
322  inline G4double LowEnergyLimit() const;
323 
324  inline G4double HighEnergyActivationLimit() const;
325 
326  inline G4double LowEnergyActivationLimit() const;
327 
328  inline G4double PolarAngleLimit() const;
329 
330  inline G4double SecondaryThreshold() const;
331 
332  inline G4bool LPMFlag() const;
333 
334  inline G4bool DeexcitationFlag() const;
335 
336  inline G4bool ForceBuildTableFlag() const;
337 
338  inline G4bool UseAngularGeneratorFlag() const;
339 
340  inline void SetAngularGeneratorFlag(G4bool);
341 
342  inline void SetHighEnergyLimit(G4double);
343 
344  inline void SetLowEnergyLimit(G4double);
345 
346  inline void SetActivationHighEnergyLimit(G4double);
347 
348  inline void SetActivationLowEnergyLimit(G4double);
349 
350  inline G4bool IsActive(G4double kinEnergy);
351 
352  inline void SetPolarAngleLimit(G4double);
353 
354  inline void SetSecondaryThreshold(G4double);
355 
356  inline void SetLPMFlag(G4bool val);
357 
358  inline void SetDeexcitationFlag(G4bool val);
359 
360  inline void SetForceBuildTable(G4bool val);
361 
362  inline void SetMasterThread(G4bool val);
363 
364  inline G4bool IsMaster() const;
365 
366  inline G4double MaxSecondaryKinEnergy(const G4DynamicParticle* dynParticle);
367 
368  inline const G4String& GetName() const;
369 
370  inline void SetCurrentCouple(const G4MaterialCutsCouple*);
371 
372  inline const G4Element* GetCurrentElement() const;
373 
374  inline const G4Isotope* GetCurrentIsotope() const;
375 
376  inline G4bool IsLocked() const;
377 
378  inline void SetLocked(G4bool);
379 
380 protected:
381 
382  inline const G4MaterialCutsCouple* CurrentCouple() const;
383 
384  inline void SetCurrentElement(const G4Element*);
385 
386 private:
387 
388  // hide assignment operator
389  G4VEmModel & operator=(const G4VEmModel &right) = delete;
390  G4VEmModel(const G4VEmModel&) = delete;
391 
392  // ======== Parameters of the class fixed at construction =========
393 
394  G4VEmFluctuationModel* flucModel;
395  G4VEmAngularDistribution* anglModel;
396  const G4String name;
397 
398  // ======== Parameters of the class fixed at initialisation =======
399 
400  G4double lowLimit;
401  G4double highLimit;
402  G4double eMinActive;
403  G4double eMaxActive;
404  G4double polarAngleLimit;
405  G4double secondaryThreshold;
406  G4bool theLPMflag;
407  G4bool flagDeexcitation;
408  G4bool flagForceBuildTable;
409  G4bool isMaster;
410 
411  G4bool localTable;
412  G4bool localElmSelectors;
413  G4bool useAngularGenerator;
414  G4bool isLocked;
415  G4int nSelectors;
416  std::vector<G4EmElementSelector*>* elmSelectors;
417  G4LossTableManager* fEmManager;
418 
419 protected:
420 
421  G4ElementData* fElementData;
422  G4VParticleChange* pParticleChange;
423  G4PhysicsTable* xSectionTable;
424  const std::vector<G4double>* theDensityFactor;
425  const std::vector<G4int>* theDensityIdx;
426  size_t idxTable;
427  const static G4double inveplus;
428 
429  // ======== Cached values - may be state dependent ================
430 
431 private:
432 
433  const G4MaterialCutsCouple* fCurrentCouple;
434  const G4Element* fCurrentElement;
435  const G4Isotope* fCurrentIsotope;
436 
437  G4int nsec;
438  std::vector<G4double> xsec;
439 
440 };
441 
442 // ======== Run time inline methods ================
443 
444 inline void G4VEmModel::SetCurrentCouple(const G4MaterialCutsCouple* p)
445 {
446  fCurrentCouple = p;
447 }
448 
449 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
450 
451 inline const G4MaterialCutsCouple* G4VEmModel::CurrentCouple() const
452 {
453  return fCurrentCouple;
454 }
455 
456 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
457 
458 inline void G4VEmModel::SetCurrentElement(const G4Element* elm)
459 {
460  fCurrentElement = elm;
461  fCurrentIsotope = nullptr;
462 }
463 
464 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
465 
466 inline const G4Element* G4VEmModel::GetCurrentElement() const
467 {
468  return fCurrentElement;
469 }
470 
471 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
472 
473 inline const G4Isotope* G4VEmModel::GetCurrentIsotope() const
474 {
475  return fCurrentIsotope;
476 }
477 
478 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
479 
480 inline
481 G4double G4VEmModel::MaxSecondaryKinEnergy(const G4DynamicParticle* dynPart)
482 {
483  return MaxSecondaryEnergy(dynPart->GetParticleDefinition(),
484  dynPart->GetKineticEnergy());
485 }
486 
487 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
488 
489 inline G4double G4VEmModel::ComputeDEDX(const G4MaterialCutsCouple* couple,
490  const G4ParticleDefinition* part,
491  G4double kinEnergy,
492  G4double cutEnergy)
493 {
494  SetCurrentCouple(couple);
495  return ComputeDEDXPerVolume(couple->GetMaterial(),part,kinEnergy,cutEnergy);
496 }
497 
498 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
499 
500 inline G4double G4VEmModel::CrossSection(const G4MaterialCutsCouple* couple,
501  const G4ParticleDefinition* part,
502  G4double kinEnergy,
503  G4double cutEnergy,
504  G4double maxEnergy)
505 {
506  SetCurrentCouple(couple);
507  return CrossSectionPerVolume(couple->GetMaterial(),part,kinEnergy,
508  cutEnergy,maxEnergy);
509 }
510 
511 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
512 
513 inline
514 G4double G4VEmModel::ComputeMeanFreePath(const G4ParticleDefinition* part,
515  G4double ekin,
516  const G4Material* material,
517  G4double emin,
518  G4double emax)
519 {
520  G4double cross = CrossSectionPerVolume(material,part,ekin,emin,emax);
521  return cross > 0.0 ? 1./cross : DBL_MAX;
522 }
523 
524 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
525 
526 inline G4double
527 G4VEmModel::ComputeCrossSectionPerAtom(const G4ParticleDefinition* part,
528  const G4Element* elm,
529  G4double kinEnergy,
530  G4double cutEnergy,
531  G4double maxEnergy)
532 {
533  SetCurrentElement(elm);
534  return ComputeCrossSectionPerAtom(part,kinEnergy,elm->GetZ(),elm->GetN(),
535  cutEnergy,maxEnergy);
536 }
537 
538 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
539 
540 inline const G4Element*
541 G4VEmModel::SelectRandomAtom(const G4MaterialCutsCouple* couple,
542  const G4ParticleDefinition* part,
543  G4double kinEnergy,
544  G4double cutEnergy,
545  G4double maxEnergy)
546 {
547  fCurrentCouple = couple;
548  if(nSelectors > 0) {
549  fCurrentElement =
550  ((*elmSelectors)[couple->GetIndex()])->SelectRandomAtom(kinEnergy);
551  } else {
552  fCurrentElement = SelectRandomAtom(couple->GetMaterial(),part,kinEnergy,
553  cutEnergy,maxEnergy);
554  }
555  fCurrentIsotope = nullptr;
556  return fCurrentElement;
557 }
558 
559 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
560 
561 inline G4int G4VEmModel::SelectRandomAtomNumber(const G4Material* mat)
562 {
563  // this algorith assumes that cross section is proportional to
564  // number electrons multiplied by number of atoms
565  size_t nn = mat->GetNumberOfElements();
566  const G4ElementVector* elmv = mat->GetElementVector();
567  G4int Z = (*elmv)[0]->GetZasInt();
568  if(1 < nn) {
569  const G4double* at = mat->GetVecNbOfAtomsPerVolume();
570  G4double tot = mat->GetTotNbOfAtomsPerVolume()*G4UniformRand();
571  for( size_t i=0; i<nn; ++i) {
572  tot -= at[i];
573  if(tot <= 0.0) {
574  Z = (*elmv)[i]->GetZasInt();
575  break;
576  }
577  }
578  }
579  return Z;
580 }
581 
582 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
583 
584 inline G4int G4VEmModel::SelectIsotopeNumber(const G4Element* elm)
585 {
586  SetCurrentElement(elm);
587  G4int N = G4lrint(elm->GetN());
588  G4int ni = elm->GetNumberOfIsotopes();
589  fCurrentIsotope = nullptr;
590  if(ni > 0) {
591  G4int idx = 0;
592  if(ni > 1) {
593  G4double* ab = elm->GetRelativeAbundanceVector();
594  G4double x = G4UniformRand();
595  for(; idx<ni; ++idx) {
596  x -= ab[idx];
597  if (x <= 0.0) { break; }
598  }
599  if(idx >= ni) { idx = ni - 1; }
600  }
601  fCurrentIsotope = elm->GetIsotope(idx);
602  N = fCurrentIsotope->GetN();
603  }
604  return N;
605 }
606 
607 // ======== Get/Set inline methods used at initialisation ================
608 
609 inline G4VEmFluctuationModel* G4VEmModel::GetModelOfFluctuations()
610 {
611  return flucModel;
612 }
613 
614 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
615 
616 inline G4VEmAngularDistribution* G4VEmModel::GetAngularDistribution()
617 {
618  return anglModel;
619 }
620 
621 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
622 
623 inline void G4VEmModel::SetAngularDistribution(G4VEmAngularDistribution* p)
624 {
625  if(p != anglModel) {
626  delete anglModel;
627  anglModel = p;
628  }
629 }
630 
631 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
632 
633 inline G4double G4VEmModel::HighEnergyLimit() const
634 {
635  return highLimit;
636 }
637 
638 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
639 
640 inline G4double G4VEmModel::LowEnergyLimit() const
641 {
642  return lowLimit;
643 }
644 
645 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
646 
647 inline G4double G4VEmModel::HighEnergyActivationLimit() const
648 {
649  return eMaxActive;
650 }
651 
652 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
653 
654 inline G4double G4VEmModel::LowEnergyActivationLimit() const
655 {
656  return eMinActive;
657 }
658 
659 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
660 
661 inline G4double G4VEmModel::PolarAngleLimit() const
662 {
663  return polarAngleLimit;
664 }
665 
666 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
667 
668 inline G4double G4VEmModel::SecondaryThreshold() const
669 {
670  return secondaryThreshold;
671 }
672 
673 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
674 
675 inline G4bool G4VEmModel::LPMFlag() const
676 {
677  return theLPMflag;
678 }
679 
680 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
681 
682 inline G4bool G4VEmModel::DeexcitationFlag() const
683 {
684  return flagDeexcitation;
685 }
686 
687 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
688 
689 inline G4bool G4VEmModel::ForceBuildTableFlag() const
690 {
691  return flagForceBuildTable;
692 }
693 
694 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
695 
696 inline G4bool G4VEmModel::UseAngularGeneratorFlag() const
697 {
698  return useAngularGenerator;
699 }
700 
701 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
702 
703 inline void G4VEmModel::SetAngularGeneratorFlag(G4bool val)
704 {
705  useAngularGenerator = val;
706 }
707 
708 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
709 
710 inline void G4VEmModel::SetMasterThread(G4bool val)
711 {
712  isMaster = val;
713 }
714 
715 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
716 
717 inline G4bool G4VEmModel::IsMaster() const
718 {
719  return isMaster;
720 }
721 
722 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
723 
724 inline void G4VEmModel::SetHighEnergyLimit(G4double val)
725 {
726  highLimit = val;
727 }
728 
729 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
730 
731 inline void G4VEmModel::SetLowEnergyLimit(G4double val)
732 {
733  lowLimit = val;
734 }
735 
736 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
737 
738 inline void G4VEmModel::SetActivationHighEnergyLimit(G4double val)
739 {
740  eMaxActive = val;
741 }
742 
743 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
744 
745 inline void G4VEmModel::SetActivationLowEnergyLimit(G4double val)
746 {
747  eMinActive = val;
748 }
749 
750 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
751 
752 inline G4bool G4VEmModel::IsActive(G4double kinEnergy)
753 {
754  return (kinEnergy >= eMinActive && kinEnergy <= eMaxActive);
755 }
756 
757 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
758 
759 inline void G4VEmModel::SetPolarAngleLimit(G4double val)
760 {
761  if(!isLocked) { polarAngleLimit = val; }
762 }
763 
764 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
765 
766 inline void G4VEmModel::SetSecondaryThreshold(G4double val)
767 {
768  secondaryThreshold = val;
769 }
770 
771 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
772 
773 inline void G4VEmModel::SetLPMFlag(G4bool val)
774 {
775  theLPMflag = val;
776 }
777 
778 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
779 
780 inline void G4VEmModel::SetDeexcitationFlag(G4bool val)
781 {
782  flagDeexcitation = val;
783 }
784 
785 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
786 
787 inline void G4VEmModel::SetForceBuildTable(G4bool val)
788 {
789  flagForceBuildTable = val;
790 }
791 
792 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
793 
794 inline const G4String& G4VEmModel::GetName() const
795 {
796  return name;
797 }
798 
799 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
800 
801 inline std::vector<G4EmElementSelector*>* G4VEmModel::GetElementSelectors()
802 {
803  return elmSelectors;
804 }
805 
806 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
807 
808 inline void
809 G4VEmModel::SetElementSelectors(std::vector<G4EmElementSelector*>* p)
810 {
811  elmSelectors = p;
812  if(elmSelectors) { nSelectors = elmSelectors->size(); }
813  localElmSelectors = false;
814 }
815 
816 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
817 
818 inline G4ElementData* G4VEmModel::GetElementData()
819 {
820  return fElementData;
821 }
822 
823 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
824 
825 inline G4PhysicsTable* G4VEmModel::GetCrossSectionTable()
826 {
827  return xSectionTable;
828 }
829 
830 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
831 
832 inline G4bool G4VEmModel::IsLocked() const
833 {
834  return isLocked;
835 }
836 
837 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
838 
839 inline void G4VEmModel::SetLocked(G4bool val)
840 {
841  isLocked = val;
842 }
843 
844 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
845 
846 #endif
847 
G4Element::GetNumberOfIsotopes
size_t GetNumberOfIsotopes() const
Definition: G4Element.hh:159
G4VEmModel::LowEnergyActivationLimit
G4double LowEnergyActivationLimit() const
Definition: G4VEmModel.hh:654
G4VEmModel::HighEnergyActivationLimit
G4double HighEnergyActivationLimit() const
Definition: G4VEmModel.hh:647
G4Isotope.hh
CLHEP::N
const int N
Definition: mixmax.h:43
G4VEmModel::MaxSecondaryKinEnergy
G4double MaxSecondaryKinEnergy(const G4DynamicParticle *dynParticle)
Definition: G4VEmModel.hh:481
G4VEmModel::idxTable
size_t idxTable
Definition: G4VEmModel.hh:426
name
const XML_Char * name
Definition: expat.h:151
G4VEmModel::SetActivationHighEnergyLimit
void SetActivationHighEnergyLimit(G4double)
Definition: G4VEmModel.hh:738
G4VEmModel::LowEnergyLimit
G4double LowEnergyLimit() const
Definition: G4VEmModel.hh:640
G4VEmModel::CrossSectionPerVolume
virtual G4double CrossSectionPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
Definition: G4VEmModel.cc:257
G4VEmModel::InitialiseForElement
virtual void InitialiseForElement(const G4ParticleDefinition *, G4int Z)
Definition: G4VEmModel.cc:243
G4VEmModel::StartTracking
virtual void StartTracking(G4Track *)
Definition: G4VEmModel.cc:292
G4VEmModel::GetParticleChangeForLoss
G4ParticleChangeForLoss * GetParticleChangeForLoss()
Definition: G4VEmModel.cc:118
G4VEmModel::SecondaryThreshold
G4double SecondaryThreshold() const
Definition: G4VEmModel.hh:668
G4ElementVector
std::vector< G4Element * > G4ElementVector
Definition: G4ElementVector.hh:44
G4VEmModel::ForceBuildTableFlag
G4bool ForceBuildTableFlag() const
Definition: G4VEmModel.hh:689
G4VEmModel::CorrectionsAlongStep
virtual void CorrectionsAlongStep(const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double &eloss, G4double &niel, G4double length)
Definition: G4VEmModel.cc:370
G4VEmModel::SampleSecondaries
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin=0.0, G4double tmax=DBL_MAX)=0
G4VEmAngularDistribution.hh
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
G4VEmModel::InitialiseElementSelectors
void InitialiseElementSelectors(const G4ParticleDefinition *, const G4DataVector &)
Definition: G4VEmModel.cc:146
G4VEmModel::HighEnergyLimit
G4double HighEnergyLimit() const
Definition: G4VEmModel.hh:633
G4Element::GetN
G4double GetN() const
Definition: G4Element.hh:135
G4VEmModel::SetParticleChange
void SetParticleChange(G4VParticleChange *, G4VEmFluctuationModel *f=nullptr)
Definition: G4VEmModel.cc:418
G4Element.hh
G4MaterialCutsCouple.hh
G4VEmModel::DefineForRegion
virtual void DefineForRegion(const G4Region *)
Definition: G4VEmModel.cc:340
p
const char * p
Definition: xmltok.h:285
G4VEmModel::MinPrimaryEnergy
virtual G4double MinPrimaryEnergy(const G4Material *, const G4ParticleDefinition *, G4double cut=0.0)
Definition: G4VEmModel.cc:386
G4VEmModel::InitialiseLocal
virtual void InitialiseLocal(const G4ParticleDefinition *, G4VEmModel *masterModel)
Definition: G4VEmModel.cc:218
G4VEmModel::SetupForMaterial
virtual void SetupForMaterial(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
Definition: G4VEmModel.cc:411
G4Element::GetZ
G4double GetZ() const
Definition: G4Element.hh:131
G4VEmModel::GetAngularDistribution
G4VEmAngularDistribution * GetAngularDistribution()
Definition: G4VEmModel.hh:616
G4DynamicParticle.hh
G4VEmModel::SetSecondaryThreshold
void SetSecondaryThreshold(G4double)
Definition: G4VEmModel.hh:766
G4VEmModel::LPMFlag
G4bool LPMFlag() const
Definition: G4VEmModel.hh:675
G4VEmModel::IsMaster
G4bool IsMaster() const
Definition: G4VEmModel.hh:717
G4VEmModel::GetCrossSectionTable
G4PhysicsTable * GetCrossSectionTable()
Definition: G4VEmModel.hh:825
G4EmElementSelector.hh
G4VEmModel::GetModelOfFluctuations
G4VEmFluctuationModel * GetModelOfFluctuations()
Definition: G4VEmModel.hh:609
G4VEmModel::MinEnergyCut
virtual G4double MinEnergyCut(const G4ParticleDefinition *, const G4MaterialCutsCouple *)
Definition: G4VEmModel.cc:395
G4VEmModel::GetElementData
G4ElementData * GetElementData()
Definition: G4VEmModel.hh:818
G4Material::GetElementVector
const G4ElementVector * GetElementVector() const
Definition: G4Material.hh:190
G4int
int G4int
Definition: G4Types.hh:78
G4VEmModel::G4VEmModel
G4VEmModel(const G4String &nam)
Definition: G4VEmModel.cc:68
G4VEmModel::SetHighEnergyLimit
void SetHighEnergyLimit(G4double)
Definition: G4VEmModel.hh:724
G4VEmModel::Initialise
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)=0
G4VEmModel::SetAngularGeneratorFlag
void SetAngularGeneratorFlag(G4bool)
Definition: G4VEmModel.hh:703
G4VEmModel::theDensityIdx
const std::vector< G4int > * theDensityIdx
Definition: G4VEmModel.hh:425
G4Material::GetVecNbOfAtomsPerVolume
const G4double * GetVecNbOfAtomsPerVolume() const
Definition: G4Material.hh:206
G4VEmModel::CurrentCouple
const G4MaterialCutsCouple * CurrentCouple() const
Definition: G4VEmModel.hh:451
G4VEmModel::InitialiseForMaterial
virtual void InitialiseForMaterial(const G4ParticleDefinition *, const G4Material *)
Definition: G4VEmModel.cc:224
G4Isotope::GetN
G4int GetN() const
Definition: G4Isotope.hh:94
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:97
A
double A(double temperature)
Definition: G4DNAElectronHoleRecombination.cc:60
G4VEmModel::ModelDescription
virtual void ModelDescription(std::ostream &outFile) const
Definition: G4VEmModel.cc:440
G4VEmModel::UseAngularGeneratorFlag
G4bool UseAngularGeneratorFlag() const
Definition: G4VEmModel.hh:696
G4VEmModel::SetLocked
void SetLocked(G4bool)
Definition: G4VEmModel.hh:839
G4bool
bool G4bool
Definition: G4Types.hh:79
G4Material.hh
G4VEmModel::CrossSection
G4double CrossSection(const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
Definition: G4VEmModel.hh:500
G4VEmModel::inveplus
static const G4double inveplus
Definition: G4VEmModel.hh:427
G4Element::GetRelativeAbundanceVector
G4double * GetRelativeAbundanceVector() const
Definition: G4Element.hh:167
G4VEmModel::~G4VEmModel
virtual ~G4VEmModel()
Definition: G4VEmModel.cc:92
G4ElementVector.hh
G4ParticleDefinition.hh
G4VEmModel::SetCrossSectionTable
void SetCrossSectionTable(G4PhysicsTable *, G4bool isLocal)
Definition: G4VEmModel.cc:426
G4VEmModel::SelectIsotopeNumber
G4int SelectIsotopeNumber(const G4Element *)
Definition: G4VEmModel.hh:584
G4DynamicParticle::GetParticleDefinition
const G4ParticleDefinition * GetParticleDefinition() const
G4VEmModel::GetPartialCrossSection
virtual G4double GetPartialCrossSection(const G4Material *, G4int level, const G4ParticleDefinition *, G4double kineticEnergy)
Definition: G4VEmModel.cc:283
G4VEmModel::GetElementSelectors
std::vector< G4EmElementSelector * > * GetElementSelectors()
Definition: G4VEmModel.hh:801
G4VEmModel::ComputeCrossSectionPerAtom
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A=0., G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
Definition: G4VEmModel.cc:321
G4VEmModel::ChargeSquareRatio
virtual G4double ChargeSquareRatio(const G4Track &)
Definition: G4VEmModel.cc:345
G4VEmModel::theDensityFactor
const std::vector< G4double > * theDensityFactor
Definition: G4VEmModel.hh:424
G4VEmModel::IsActive
G4bool IsActive(G4double kinEnergy)
Definition: G4VEmModel.hh:752
G4VEmModel::SetActivationLowEnergyLimit
void SetActivationLowEnergyLimit(G4double)
Definition: G4VEmModel.hh:745
emax
static const G4double emax
Definition: G4ChatterjeeCrossSection.hh:39
G4VEmModel::SetMasterThread
void SetMasterThread(G4bool val)
Definition: G4VEmModel.hh:710
G4Material::GetTotNbOfAtomsPerVolume
G4double GetTotNbOfAtomsPerVolume() const
Definition: G4Material.hh:209
G4VEmModel::GetChargeSquareRatio
virtual G4double GetChargeSquareRatio(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
Definition: G4VEmModel.cc:353
G4VEmModel::GetParticleCharge
virtual G4double GetParticleCharge(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
Definition: G4VEmModel.cc:362
G4VEmModel::DeexcitationFlag
G4bool DeexcitationFlag() const
Definition: G4VEmModel.hh:682
G4VEmModel::ComputeDEDXPerVolume
virtual G4double ComputeDEDXPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=DBL_MAX)
Definition: G4VEmModel.cc:248
G4VEmModel::SetElementSelectors
void SetElementSelectors(std::vector< G4EmElementSelector * > *)
Definition: G4VEmModel.hh:809
ab
static const G4double ab
Definition: G4ElectroNuclearCrossSection.cc:87
G4VEmModel::ComputeDEDX
virtual G4double ComputeDEDX(const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=DBL_MAX)
Definition: G4VEmModel.hh:489
G4lrint
int G4lrint(double ad)
Definition: templates.hh:163
G4VEmModel::SetLPMFlag
void SetLPMFlag(G4bool val)
Definition: G4VEmModel.hh:773
G4VEmFluctuationModel.hh
G4VEmModel::SetCurrentCouple
void SetCurrentCouple(const G4MaterialCutsCouple *)
Definition: G4VEmModel.hh:444
G4VEmModel::SetAngularDistribution
void SetAngularDistribution(G4VEmAngularDistribution *)
Definition: G4VEmModel.hh:623
G4VEmModel::PolarAngleLimit
G4double PolarAngleLimit() const
Definition: G4VEmModel.hh:661
G4VEmModel::ComputeCrossSectionPerShell
virtual G4double ComputeCrossSectionPerShell(const G4ParticleDefinition *, G4int Z, G4int shellIdx, G4double kinEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
Definition: G4VEmModel.cc:331
G4VEmModel::SetForceBuildTable
void SetForceBuildTable(G4bool val)
Definition: G4VEmModel.hh:787
G4Element::GetIsotope
const G4Isotope * GetIsotope(G4int iso) const
Definition: G4Element.hh:170
G4VEmModel::pParticleChange
G4VParticleChange * pParticleChange
Definition: G4VEmModel.hh:422
G4VEmModel::GetName
const G4String & GetName() const
Definition: G4VEmModel.hh:794
G4VEmModel::MaxSecondaryEnergy
virtual G4double MaxSecondaryEnergy(const G4ParticleDefinition *, G4double kineticEnergy)
Definition: G4VEmModel.cc:403
G4Material::GetNumberOfElements
size_t GetNumberOfElements() const
Definition: G4Material.hh:186
G4VEmModel::IsLocked
G4bool IsLocked() const
Definition: G4VEmModel.hh:832
RandomEngine.h
G4VEmModel::xSectionTable
G4PhysicsTable * xSectionTable
Definition: G4VEmModel.hh:423
G4double
double G4double
Definition: G4Types.hh:76
G4VEmModel::SetLowEnergyLimit
void SetLowEnergyLimit(G4double)
Definition: G4VEmModel.hh:731
G4VEmModel::SetDeexcitationFlag
void SetDeexcitationFlag(G4bool val)
Definition: G4VEmModel.hh:780
G4DataVector.hh
G4VEmModel::ComputeMeanFreePath
G4double ComputeMeanFreePath(const G4ParticleDefinition *, G4double kineticEnergy, const G4Material *, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
Definition: G4VEmModel.hh:514
G4VEmModel::SetCurrentElement
void SetCurrentElement(const G4Element *)
Definition: G4VEmModel.hh:458
DBL_MAX
#define DBL_MAX
Definition: templates.hh:83
G4VEmModel::fElementData
G4ElementData * fElementData
Definition: G4VEmModel.hh:421
G4VEmModel::SelectRandomAtom
const G4Element * SelectRandomAtom(const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
Definition: G4VEmModel.hh:541
G4VEmModel::SetPolarAngleLimit
void SetPolarAngleLimit(G4double)
Definition: G4VEmModel.hh:759
G4VEmModel::GetParticleChangeForGamma
G4ParticleChangeForGamma * GetParticleChangeForGamma()
Definition: G4VEmModel.cc:132
Z
G4int Z
Definition: G4ParticleHPJENDLHEData.cc:262
G4VEmModel::Value
virtual G4double Value(const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy)
Definition: G4VEmModel.cc:377
G4MaterialCutsCouple::GetMaterial
const G4Material * GetMaterial() const
Definition: G4MaterialCutsCouple.hh:150
G4VEmModel::SelectRandomAtomNumber
G4int SelectRandomAtomNumber(const G4Material *)
Definition: G4VEmModel.hh:561
G4VEmModel::GetCurrentIsotope
const G4Isotope * GetCurrentIsotope() const
Definition: G4VEmModel.hh:473
G4VEmModel::GetCurrentElement
const G4Element * GetCurrentElement() const
Definition: G4VEmModel.hh:466