Geant4  10.00.p01
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 76333 2013-11-08 14:31:50Z 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 "G4DataVector.hh"
88 #include "G4VEmFluctuationModel.hh"
89 #include "G4VEmAngularDistribution.hh"
90 #include "G4EmElementSelector.hh"
91 #include "Randomize.hh"
92 #include <vector>
93 
94 class G4ElementData;
95 class G4PhysicsTable;
96 class G4Region;
97 class G4VParticleChange;
98 class G4ParticleChangeForLoss;
99 class G4ParticleChangeForGamma;
100 class G4Track;
101 class G4LossTableManager;
102 
103 class G4VEmModel
104 {
105 
106 public:
107 
108  G4VEmModel(const G4String& nam);
109 
110  virtual ~G4VEmModel();
111 
112  //------------------------------------------------------------------------
113  // Virtual methods to be implemented for any concrete model
114  //------------------------------------------------------------------------
115 
116  virtual void Initialise(const G4ParticleDefinition*,
117  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  // main method to compute cross section per atom
159  virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition*,
160  G4double kinEnergy,
161  G4double Z,
162  G4double A = 0., /* amu */
163  G4double cutEnergy = 0.0,
164  G4double maxEnergy = DBL_MAX);
165 
166  // Compute effective ion charge square
167  virtual G4double ChargeSquareRatio(const G4Track&);
168 
169  // Compute effective ion charge square
170  virtual G4double GetChargeSquareRatio(const G4ParticleDefinition*,
171  const G4Material*,
172  G4double kineticEnergy);
173 
174  // Compute ion charge
175  virtual G4double GetParticleCharge(const G4ParticleDefinition*,
176  const G4Material*,
177  G4double kineticEnergy);
178 
179  // Initialisation for a new track
180  virtual void StartTracking(G4Track*);
181 
182  // add correction to energy loss and compute non-ionizing energy loss
183  virtual void CorrectionsAlongStep(const G4MaterialCutsCouple*,
184  const G4DynamicParticle*,
185  G4double& eloss,
186  G4double& niel,
187  G4double length);
188 
189  // value which may be tabulated (by default cross section)
190  virtual G4double Value(const G4MaterialCutsCouple*,
191  const G4ParticleDefinition*,
192  G4double kineticEnergy);
193 
194  // threshold for zero value
195  virtual G4double MinPrimaryEnergy(const G4Material*,
196  const G4ParticleDefinition*,
197  G4double cut = 0.0);
198 
199  // model can define low-energy limit for the cut
200  virtual G4double MinEnergyCut(const G4ParticleDefinition*,
201  const G4MaterialCutsCouple*);
202 
203  // initilisation at run time for a given material
204  virtual void SetupForMaterial(const G4ParticleDefinition*,
205  const G4Material*,
206  G4double kineticEnergy);
207 
208  // add a region for the model
209  virtual void DefineForRegion(const G4Region*);
210 
211 protected:
212 
213  // initialisation of the ParticleChange for the model
214  G4ParticleChangeForLoss* GetParticleChangeForLoss();
215 
216  // initialisation of the ParticleChange for the model
217  G4ParticleChangeForGamma* GetParticleChangeForGamma();
218 
219  // kinematically allowed max kinetic energy of a secondary
220  virtual G4double MaxSecondaryEnergy(const G4ParticleDefinition*,
221  G4double kineticEnergy);
222 
223 public:
224 
225  //------------------------------------------------------------------------
226  // Generic methods common to all models
227  //------------------------------------------------------------------------
228 
229  // should be called at initialisation to build element selectors
230  void InitialiseElementSelectors(const G4ParticleDefinition*,
231  const G4DataVector&);
232 
233  // should be called at initialisation to access element selectors
234  inline std::vector<G4EmElementSelector*>* GetElementSelectors();
235 
236  // should be called at initialisation to set element selectors
237  inline void SetElementSelectors(std::vector<G4EmElementSelector*>*);
238 
239  // dEdx per unit length
240  inline G4double ComputeDEDX(const G4MaterialCutsCouple*,
241  const G4ParticleDefinition*,
242  G4double kineticEnergy,
243  G4double cutEnergy = DBL_MAX);
244 
245  // cross section per volume
246  inline G4double CrossSection(const G4MaterialCutsCouple*,
247  const G4ParticleDefinition*,
248  G4double kineticEnergy,
249  G4double cutEnergy = 0.0,
250  G4double maxEnergy = DBL_MAX);
251 
252  // compute mean free path via cross section per volume
253  inline G4double ComputeMeanFreePath(const G4ParticleDefinition*,
254  G4double kineticEnergy,
255  const G4Material*,
256  G4double cutEnergy = 0.0,
257  G4double maxEnergy = DBL_MAX);
258 
259  // generic cross section per element
260  inline G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition*,
261  const G4Element*,
262  G4double kinEnergy,
263  G4double cutEnergy = 0.0,
264  G4double maxEnergy = DBL_MAX);
265 
266  // select isotope in order to have precise mass of the nucleus
267  inline G4int SelectIsotopeNumber(const G4Element*);
268 
269  // atom can be selected effitiantly if element selectors are initialised
270  inline const G4Element* SelectRandomAtom(const G4MaterialCutsCouple*,
271  const G4ParticleDefinition*,
272  G4double kineticEnergy,
273  G4double cutEnergy = 0.0,
274  G4double maxEnergy = DBL_MAX);
275 
276  // to select atom cross section per volume is recomputed for each element
277  const G4Element* SelectRandomAtom(const G4Material*,
278  const G4ParticleDefinition*,
279  G4double kineticEnergy,
280  G4double cutEnergy = 0.0,
281  G4double maxEnergy = DBL_MAX);
282 
283  // to select atom if cross section is proportional number of electrons
284  inline G4int SelectRandomAtomNumber(const G4Material*);
285 
286  //------------------------------------------------------------------------
287  // Get/Set methods
288  //------------------------------------------------------------------------
289 
290  void SetParticleChange(G4VParticleChange*, G4VEmFluctuationModel* f=0);
291 
292  void SetCrossSectionTable(G4PhysicsTable*, G4bool isLocal);
293 
294  inline G4ElementData* GetElementData();
295 
296  inline G4PhysicsTable* GetCrossSectionTable();
297 
298  inline G4VEmFluctuationModel* GetModelOfFluctuations();
299 
300  inline G4VEmAngularDistribution* GetAngularDistribution();
301 
302  inline void SetAngularDistribution(G4VEmAngularDistribution*);
303 
304  inline G4double HighEnergyLimit() const;
305 
306  inline G4double LowEnergyLimit() const;
307 
308  inline G4double HighEnergyActivationLimit() const;
309 
310  inline G4double LowEnergyActivationLimit() const;
311 
312  inline G4double PolarAngleLimit() const;
313 
314  inline G4double SecondaryThreshold() const;
315 
316  inline G4bool LPMFlag() const;
317 
318  inline G4bool DeexcitationFlag() const;
319 
320  inline G4bool ForceBuildTableFlag() const;
321 
322  inline G4bool UseAngularGeneratorFlag() const;
323 
324  inline void SetAngularGeneratorFlag(G4bool);
325 
326  inline void SetHighEnergyLimit(G4double);
327 
328  inline void SetLowEnergyLimit(G4double);
329 
330  inline void SetActivationHighEnergyLimit(G4double);
331 
332  inline void SetActivationLowEnergyLimit(G4double);
333 
334  inline G4bool IsActive(G4double kinEnergy);
335 
336  inline void SetPolarAngleLimit(G4double);
337 
338  inline void SetSecondaryThreshold(G4double);
339 
340  inline void SetLPMFlag(G4bool val);
341 
342  inline void SetDeexcitationFlag(G4bool val);
343 
344  inline void SetForceBuildTable(G4bool val);
345 
346  inline void SetMasterThread(G4bool val);
347 
348  inline G4bool IsMaster() const;
349 
350  inline G4double MaxSecondaryKinEnergy(const G4DynamicParticle* dynParticle);
351 
352  inline const G4String& GetName() const;
353 
354  inline void SetCurrentCouple(const G4MaterialCutsCouple*);
355 
356  inline const G4Element* GetCurrentElement() const;
357 
358 protected:
359 
360  inline const G4MaterialCutsCouple* CurrentCouple() const;
361 
362  inline void SetCurrentElement(const G4Element*);
363 
364 private:
365 
366  // hide assignment operator
367  G4VEmModel & operator=(const G4VEmModel &right);
368  G4VEmModel(const G4VEmModel&);
369 
370  // ======== Parameters of the class fixed at construction =========
371 
372  G4VEmFluctuationModel* flucModel;
373  G4VEmAngularDistribution* anglModel;
374  const G4String name;
375 
376  // ======== Parameters of the class fixed at initialisation =======
377 
378  G4double lowLimit;
379  G4double highLimit;
380  G4double eMinActive;
381  G4double eMaxActive;
382  G4double polarAngleLimit;
383  G4double secondaryThreshold;
384  G4bool theLPMflag;
385  G4bool flagDeexcitation;
386  G4bool flagForceBuildTable;
387  G4bool isMaster;
388 
389  G4bool localTable;
390  G4bool localElmSelectors;
391  G4bool useAngularGenerator;
392  G4int nSelectors;
393  std::vector<G4EmElementSelector*>* elmSelectors;
394 
395 protected:
396 
397  G4ElementData* fElementData;
398  G4VParticleChange* pParticleChange;
399  G4PhysicsTable* xSectionTable;
400  const std::vector<G4double>* theDensityFactor;
401  const std::vector<G4int>* theDensityIdx;
402  size_t idxTable;
403 
404  // ======== Cashed values - may be state dependent ================
405 
406 private:
407 
408  G4LossTableManager* fManager;
409  const G4MaterialCutsCouple* fCurrentCouple;
410  const G4Element* fCurrentElement;
411 
412  G4int nsec;
413  std::vector<G4double> xsec;
414 
415 };
416 
417 // ======== Run time inline methods ================
418 
419 inline void G4VEmModel::SetCurrentCouple(const G4MaterialCutsCouple* p)
420 {
421  fCurrentCouple = p;
422 }
423 
424 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
425 
426 inline const G4MaterialCutsCouple* G4VEmModel::CurrentCouple() const
427 {
428  return fCurrentCouple;
429 }
430 
431 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
432 
433 inline void G4VEmModel::SetCurrentElement(const G4Element* elm)
434 {
435  fCurrentElement = elm;
436 }
437 
438 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
439 
440 inline const G4Element* G4VEmModel::GetCurrentElement() const
441 {
442  return fCurrentElement;
443 }
444 
445 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
446 
447 inline
448 G4double G4VEmModel::MaxSecondaryKinEnergy(const G4DynamicParticle* dynPart)
449 {
450  return MaxSecondaryEnergy(dynPart->GetParticleDefinition(),
451  dynPart->GetKineticEnergy());
452 }
453 
454 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
455 
456 inline G4double G4VEmModel::ComputeDEDX(const G4MaterialCutsCouple* couple,
457  const G4ParticleDefinition* part,
458  G4double kinEnergy,
459  G4double cutEnergy)
460 {
461  SetCurrentCouple(couple);
462  return ComputeDEDXPerVolume(couple->GetMaterial(),part,kinEnergy,cutEnergy);
463 }
464 
465 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
466 
467 inline G4double G4VEmModel::CrossSection(const G4MaterialCutsCouple* couple,
468  const G4ParticleDefinition* part,
469  G4double kinEnergy,
470  G4double cutEnergy,
471  G4double maxEnergy)
472 {
473  SetCurrentCouple(couple);
474  return CrossSectionPerVolume(couple->GetMaterial(),part,kinEnergy,
475  cutEnergy,maxEnergy);
476 }
477 
478 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
479 
480 inline
481 G4double G4VEmModel::ComputeMeanFreePath(const G4ParticleDefinition* part,
482  G4double ekin,
483  const G4Material* material,
484  G4double emin,
485  G4double emax)
486 {
487  G4double mfp = DBL_MAX;
488  G4double cross = CrossSectionPerVolume(material,part,ekin,emin,emax);
489  if (cross > DBL_MIN) { mfp = 1./cross; }
490  return mfp;
491 }
492 
493 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
494 
495 inline G4double G4VEmModel::ComputeCrossSectionPerAtom(
496  const G4ParticleDefinition* part,
497  const G4Element* elm,
498  G4double kinEnergy,
499  G4double cutEnergy,
500  G4double maxEnergy)
501 {
502  SetCurrentElement(elm);
503  return ComputeCrossSectionPerAtom(part,kinEnergy,elm->GetZ(),elm->GetN(),
504  cutEnergy,maxEnergy);
505 }
506 
507 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
508 
509 inline const G4Element*
510 G4VEmModel::SelectRandomAtom(const G4MaterialCutsCouple* couple,
511  const G4ParticleDefinition* part,
512  G4double kinEnergy,
513  G4double cutEnergy,
514  G4double maxEnergy)
515 {
516  fCurrentCouple = couple;
517  if(nSelectors > 0) {
518  fCurrentElement =
519  ((*elmSelectors)[couple->GetIndex()])->SelectRandomAtom(kinEnergy);
520  } else {
521  fCurrentElement = SelectRandomAtom(couple->GetMaterial(),part,kinEnergy,
522  cutEnergy,maxEnergy);
523  }
524  return fCurrentElement;
525 }
526 
527 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
528 
529 inline G4int G4VEmModel::SelectRandomAtomNumber(const G4Material* mat)
530 {
531  size_t nn = mat->GetNumberOfElements();
532  const G4ElementVector* elmv = mat->GetElementVector();
533  G4int Z = G4int((*elmv)[0]->GetZ());
534  if(1 < nn) {
535  const G4double* at = mat->GetVecNbOfAtomsPerVolume();
536  G4double tot = mat->GetTotNbOfAtomsPerVolume()*G4UniformRand();
537  for( size_t i=0; i<nn; ++i) {
538  Z = G4int((*elmv)[0]->GetZ());
539  tot -= Z*at[i];
540  if(tot <= 0.0) { break; }
541  }
542  }
543  return Z;
544 }
545 
546 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
547 
548 inline G4int G4VEmModel::SelectIsotopeNumber(const G4Element* elm)
549 {
550  SetCurrentElement(elm);
551  G4int N = G4int(elm->GetN() + 0.5);
552  G4int ni = elm->GetNumberOfIsotopes();
553  if(ni > 0) {
554  G4int idx = 0;
555  if(ni > 1) {
556  G4double* ab = elm->GetRelativeAbundanceVector();
557  G4double x = G4UniformRand();
558  for(; idx<ni; ++idx) {
559  x -= ab[idx];
560  if (x <= 0.0) { break; }
561  }
562  if(idx >= ni) { idx = ni - 1; }
563  }
564  N = elm->GetIsotope(idx)->GetN();
565  }
566  return N;
567 }
568 
569 // ======== Get/Set inline methods used at initialisation ================
570 
571 inline G4VEmFluctuationModel* G4VEmModel::GetModelOfFluctuations()
572 {
573  return flucModel;
574 }
575 
576 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
577 
578 inline G4VEmAngularDistribution* G4VEmModel::GetAngularDistribution()
579 {
580  return anglModel;
581 }
582 
583 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
584 
585 inline void G4VEmModel::SetAngularDistribution(G4VEmAngularDistribution* p)
586 {
587  anglModel = p;
588 }
589 
590 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
591 
592 inline G4double G4VEmModel::HighEnergyLimit() const
593 {
594  return highLimit;
595 }
596 
597 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
598 
599 inline G4double G4VEmModel::LowEnergyLimit() const
600 {
601  return lowLimit;
602 }
603 
604 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
605 
606 inline G4double G4VEmModel::HighEnergyActivationLimit() const
607 {
608  return eMaxActive;
609 }
610 
611 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
612 
613 inline G4double G4VEmModel::LowEnergyActivationLimit() const
614 {
615  return eMinActive;
616 }
617 
618 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
619 
620 inline G4double G4VEmModel::PolarAngleLimit() const
621 {
622  return polarAngleLimit;
623 }
624 
625 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
626 
627 inline G4double G4VEmModel::SecondaryThreshold() const
628 {
629  return secondaryThreshold;
630 }
631 
632 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
633 
634 inline G4bool G4VEmModel::LPMFlag() const
635 {
636  return theLPMflag;
637 }
638 
639 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
640 
641 inline G4bool G4VEmModel::DeexcitationFlag() const
642 {
643  return flagDeexcitation;
644 }
645 
646 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
647 
648 inline G4bool G4VEmModel::ForceBuildTableFlag() const
649 {
650  return flagForceBuildTable;
651 }
652 
653 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
654 
655 inline G4bool G4VEmModel::UseAngularGeneratorFlag() const
656 {
657  return useAngularGenerator;
658 }
659 
660 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
661 
662 inline void G4VEmModel::SetAngularGeneratorFlag(G4bool val)
663 {
664  useAngularGenerator = val;
665 }
666 
667 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
668 
669 inline void G4VEmModel::SetMasterThread(G4bool val)
670 {
671  isMaster = val;
672 }
673 
674 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
675 
676 inline G4bool G4VEmModel::IsMaster() const
677 {
678  return isMaster;
679 }
680 
681 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
682 
683 inline void G4VEmModel::SetHighEnergyLimit(G4double val)
684 {
685  highLimit = val;
686 }
687 
688 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
689 
690 inline void G4VEmModel::SetLowEnergyLimit(G4double val)
691 {
692  lowLimit = val;
693 }
694 
695 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
696 
697 inline void G4VEmModel::SetActivationHighEnergyLimit(G4double val)
698 {
699  eMaxActive = val;
700 }
701 
702 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
703 
704 inline void G4VEmModel::SetActivationLowEnergyLimit(G4double val)
705 {
706  eMinActive = val;
707 }
708 
709 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
710 
711 inline G4bool G4VEmModel::IsActive(G4double kinEnergy)
712 {
713  return (kinEnergy >= eMinActive && kinEnergy <= eMaxActive);
714 }
715 
716 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
717 
718 inline void G4VEmModel::SetPolarAngleLimit(G4double val)
719 {
720  polarAngleLimit = val;
721 }
722 
723 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
724 
725 inline void G4VEmModel::SetSecondaryThreshold(G4double val)
726 {
727  secondaryThreshold = val;
728 }
729 
730 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
731 
732 inline void G4VEmModel::SetLPMFlag(G4bool val)
733 {
734  theLPMflag = val;
735 }
736 
737 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
738 
739 inline void G4VEmModel::SetDeexcitationFlag(G4bool val)
740 {
741  flagDeexcitation = val;
742 }
743 
744 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
745 
746 inline void G4VEmModel::SetForceBuildTable(G4bool val)
747 {
748  flagForceBuildTable = val;
749 }
750 
751 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
752 
753 inline const G4String& G4VEmModel::GetName() const
754 {
755  return name;
756 }
757 
758 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
759 
760 inline std::vector<G4EmElementSelector*>* G4VEmModel::GetElementSelectors()
761 {
762  return elmSelectors;
763 }
764 
765 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
766 
767 inline void
768 G4VEmModel::SetElementSelectors(std::vector<G4EmElementSelector*>* p)
769 {
770  elmSelectors = p;
771  if(elmSelectors) { nSelectors = elmSelectors->size(); }
772  localElmSelectors = false;
773 }
774 
775 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
776 
777 inline G4ElementData* G4VEmModel::GetElementData()
778 {
779  return fElementData;
780 }
781 
782 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
783 
784 inline G4PhysicsTable* G4VEmModel::GetCrossSectionTable()
785 {
786  return xSectionTable;
787 }
788 
789 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
790 
791 #endif
792 
G4Element::GetNumberOfIsotopes
size_t GetNumberOfIsotopes() const
Definition: G4Element.hh:158
G4VEmModel::LowEnergyActivationLimit
G4double LowEnergyActivationLimit() const
Definition: G4VEmModel.hh:613
G4VEmModel::HighEnergyActivationLimit
G4double HighEnergyActivationLimit() const
Definition: G4VEmModel.hh:606
G4VEmModel::localTable
G4bool localTable
Definition: G4VEmModel.hh:389
G4VEmModel::MaxSecondaryKinEnergy
G4double MaxSecondaryKinEnergy(const G4DynamicParticle *dynParticle)
Definition: G4VEmModel.hh:448
G4VEmModel::idxTable
size_t idxTable
Definition: G4VEmModel.hh:402
G4VEmModel::SetActivationHighEnergyLimit
void SetActivationHighEnergyLimit(G4double)
Definition: G4VEmModel.hh:697
G4VEmModel::flagForceBuildTable
G4bool flagForceBuildTable
Definition: G4VEmModel.hh:386
G4VEmModel::LowEnergyLimit
G4double LowEnergyLimit() const
Definition: G4VEmModel.hh:599
G4VEmModel::CrossSectionPerVolume
virtual G4double CrossSectionPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
Definition: G4VEmModel.cc:245
G4VEmModel::InitialiseForElement
virtual void InitialiseForElement(const G4ParticleDefinition *, G4int Z)
Definition: G4VEmModel.cc:231
G4VEmModel::StartTracking
virtual void StartTracking(G4Track *)
Definition: G4VEmModel.cc:271
G4VEmModel::GetParticleChangeForLoss
G4ParticleChangeForLoss * GetParticleChangeForLoss()
Definition: G4VEmModel.cc:107
G4VEmModel::SecondaryThreshold
G4double SecondaryThreshold() const
Definition: G4VEmModel.hh:627
right
Definition: F04UserTrackInformation.hh:37
G4ElementVector
std::vector< G4Element * > G4ElementVector
Definition: G4ElementVector.hh:44
G4VEmModel::ForceBuildTableFlag
G4bool ForceBuildTableFlag() const
Definition: G4VEmModel.hh:648
G4VEmModel::CorrectionsAlongStep
virtual void CorrectionsAlongStep(const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double &eloss, G4double &niel, G4double length)
Definition: G4VEmModel.cc:339
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:135
G4VEmModel::HighEnergyLimit
G4double HighEnergyLimit() const
Definition: G4VEmModel.hh:592
G4VEmModel::flagDeexcitation
G4bool flagDeexcitation
Definition: G4VEmModel.hh:385
G4Element::GetN
G4double GetN() const
Definition: G4Element.hh:134
G4Element.hh
G4VEmModel::theLPMflag
G4bool theLPMflag
Definition: G4VEmModel.hh:384
G4MaterialCutsCouple.hh
G4VEmModel::DefineForRegion
virtual void DefineForRegion(const G4Region *)
Definition: G4VEmModel.cc:309
G4VEmModel::isMaster
G4bool isMaster
Definition: G4VEmModel.hh:387
G4VEmModel::MinPrimaryEnergy
virtual G4double MinPrimaryEnergy(const G4Material *, const G4ParticleDefinition *, G4double cut=0.0)
Definition: G4VEmModel.cc:355
G4VEmModel::InitialiseLocal
virtual void InitialiseLocal(const G4ParticleDefinition *, G4VEmModel *masterModel)
Definition: G4VEmModel.cc:206
G4VEmModel::SetupForMaterial
virtual void SetupForMaterial(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
Definition: G4VEmModel.cc:380
G4Element::GetZ
G4double GetZ() const
Definition: G4Element.hh:131
G4VEmModel::GetAngularDistribution
G4VEmAngularDistribution * GetAngularDistribution()
Definition: G4VEmModel.hh:578
G4DynamicParticle.hh
G4VEmModel::SetSecondaryThreshold
void SetSecondaryThreshold(G4double)
Definition: G4VEmModel.hh:725
G4VEmModel::LPMFlag
G4bool LPMFlag() const
Definition: G4VEmModel.hh:634
G4VEmModel::IsMaster
G4bool IsMaster() const
Definition: G4VEmModel.hh:676
G4VEmModel::GetCrossSectionTable
G4PhysicsTable * GetCrossSectionTable()
Definition: G4VEmModel.hh:784
G4EmElementSelector.hh
G4VEmModel::name
const G4String name
Definition: G4VEmModel.hh:374
G4VEmModel::GetModelOfFluctuations
G4VEmFluctuationModel * GetModelOfFluctuations()
Definition: G4VEmModel.hh:571
G4VEmModel::MinEnergyCut
virtual G4double MinEnergyCut(const G4ParticleDefinition *, const G4MaterialCutsCouple *)
Definition: G4VEmModel.cc:364
G4VEmModel::GetElementData
G4ElementData * GetElementData()
Definition: G4VEmModel.hh:777
G4Material::GetElementVector
const G4ElementVector * GetElementVector() const
Definition: G4Material.hh:188
G4int
int G4int
Definition: G4Types.hh:78
G4VEmModel::fCurrentElement
const G4Element * fCurrentElement
Definition: G4VEmModel.hh:410
G4VEmModel::flucModel
G4VEmFluctuationModel * flucModel
Definition: G4VEmModel.hh:372
G4VEmModel::G4VEmModel
G4VEmModel(const G4String &nam)
Definition: G4VEmModel.cc:65
G4VEmModel::nSelectors
G4int nSelectors
Definition: G4VEmModel.hh:392
G4VEmModel::fManager
G4LossTableManager * fManager
Definition: G4VEmModel.hh:408
G4VEmModel::SetHighEnergyLimit
void SetHighEnergyLimit(G4double)
Definition: G4VEmModel.hh:683
G4VEmModel::Initialise
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)=0
G4VEmModel::SetAngularGeneratorFlag
void SetAngularGeneratorFlag(G4bool)
Definition: G4VEmModel.hh:662
G4VEmModel::theDensityIdx
const std::vector< G4int > * theDensityIdx
Definition: G4VEmModel.hh:401
G4Material::GetVecNbOfAtomsPerVolume
const G4double * GetVecNbOfAtomsPerVolume() const
Definition: G4Material.hh:204
G4VEmModel::eMaxActive
G4double eMaxActive
Definition: G4VEmModel.hh:381
G4VEmModel::CurrentCouple
const G4MaterialCutsCouple * CurrentCouple() const
Definition: G4VEmModel.hh:426
G4VEmModel::InitialiseForMaterial
virtual void InitialiseForMaterial(const G4ParticleDefinition *, const G4Material *)
Definition: G4VEmModel.cc:212
G4Isotope::GetN
G4int GetN() const
Definition: G4Isotope.hh:94
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:87
G4VEmModel::secondaryThreshold
G4double secondaryThreshold
Definition: G4VEmModel.hh:383
G4VEmModel::SetParticleChange
void SetParticleChange(G4VParticleChange *, G4VEmFluctuationModel *f=0)
Definition: G4VEmModel.cc:387
G4VEmModel::fCurrentCouple
const G4MaterialCutsCouple * fCurrentCouple
Definition: G4VEmModel.hh:409
G4VEmModel::polarAngleLimit
G4double polarAngleLimit
Definition: G4VEmModel.hh:382
G4VEmModel::UseAngularGeneratorFlag
G4bool UseAngularGeneratorFlag() const
Definition: G4VEmModel.hh:655
G4VEmModel::elmSelectors
std::vector< G4EmElementSelector * > * elmSelectors
Definition: G4VEmModel.hh:393
G4bool
bool G4bool
Definition: G4Types.hh:79
G4Material.hh
G4VEmModel::operator=
G4VEmModel & operator=(const G4VEmModel &right)
G4VEmModel::CrossSection
G4double CrossSection(const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
Definition: G4VEmModel.hh:467
G4Element::GetRelativeAbundanceVector
G4double * GetRelativeAbundanceVector() const
Definition: G4Element.hh:166
G4VEmModel::~G4VEmModel
virtual ~G4VEmModel()
Definition: G4VEmModel.cc:88
G4ElementVector.hh
G4ParticleDefinition.hh
G4VEmModel::anglModel
G4VEmAngularDistribution * anglModel
Definition: G4VEmModel.hh:373
G4VEmModel::SetCrossSectionTable
void SetCrossSectionTable(G4PhysicsTable *, G4bool isLocal)
Definition: G4VEmModel.cc:395
Randomize.hh
G4VEmModel::SelectIsotopeNumber
G4int SelectIsotopeNumber(const G4Element *)
Definition: G4VEmModel.hh:548
G4DynamicParticle::GetParticleDefinition
const G4ParticleDefinition * GetParticleDefinition() const
G4VEmModel::GetElementSelectors
std::vector< G4EmElementSelector * > * GetElementSelectors()
Definition: G4VEmModel.hh:760
A
static const G4double A[nN]
Definition: G4ElectroNuclearCrossSection.cc:115
G4VEmModel::nsec
G4int nsec
Definition: G4VEmModel.hh:412
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:300
G4VEmModel::ChargeSquareRatio
virtual G4double ChargeSquareRatio(const G4Track &)
Definition: G4VEmModel.cc:314
G4VEmModel::theDensityFactor
const std::vector< G4double > * theDensityFactor
Definition: G4VEmModel.hh:400
G4VEmModel::IsActive
G4bool IsActive(G4double kinEnergy)
Definition: G4VEmModel.hh:711
G4VEmModel::SetActivationLowEnergyLimit
void SetActivationLowEnergyLimit(G4double)
Definition: G4VEmModel.hh:704
G4VEmModel::SetMasterThread
void SetMasterThread(G4bool val)
Definition: G4VEmModel.hh:669
G4VEmModel::eMinActive
G4double eMinActive
Definition: G4VEmModel.hh:380
G4VEmModel::xsec
std::vector< G4double > xsec
Definition: G4VEmModel.hh:413
G4Material::GetTotNbOfAtomsPerVolume
G4double GetTotNbOfAtomsPerVolume() const
Definition: G4Material.hh:207
G4VEmModel::GetChargeSquareRatio
virtual G4double GetChargeSquareRatio(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
Definition: G4VEmModel.cc:322
G4VEmModel::GetParticleCharge
virtual G4double GetParticleCharge(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
Definition: G4VEmModel.cc:331
G4VEmModel::DeexcitationFlag
G4bool DeexcitationFlag() const
Definition: G4VEmModel.hh:641
G4VEmModel::ComputeDEDXPerVolume
virtual G4double ComputeDEDXPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=DBL_MAX)
Definition: G4VEmModel.cc:236
G4VEmModel::SetElementSelectors
void SetElementSelectors(std::vector< G4EmElementSelector * > *)
Definition: G4VEmModel.hh:768
ab
static const G4double ab
Definition: G4ElectroNuclearCrossSection.cc:87
G4VEmModel::ComputeDEDX
G4double ComputeDEDX(const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=DBL_MAX)
Definition: G4VEmModel.hh:456
G4VEmModel::localElmSelectors
G4bool localElmSelectors
Definition: G4VEmModel.hh:390
G4VEmModel::SetLPMFlag
void SetLPMFlag(G4bool val)
Definition: G4VEmModel.hh:732
G4VEmFluctuationModel.hh
G4VEmModel::SetCurrentCouple
void SetCurrentCouple(const G4MaterialCutsCouple *)
Definition: G4VEmModel.hh:419
G4VEmModel::SetAngularDistribution
void SetAngularDistribution(G4VEmAngularDistribution *)
Definition: G4VEmModel.hh:585
G4VEmModel::highLimit
G4double highLimit
Definition: G4VEmModel.hh:379
G4VEmModel::useAngularGenerator
G4bool useAngularGenerator
Definition: G4VEmModel.hh:391
G4VEmModel::PolarAngleLimit
G4double PolarAngleLimit() const
Definition: G4VEmModel.hh:620
DBL_MIN
#define DBL_MIN
Definition: templates.hh:75
G4VEmModel::lowLimit
G4double lowLimit
Definition: G4VEmModel.hh:378
G4VEmModel::SetForceBuildTable
void SetForceBuildTable(G4bool val)
Definition: G4VEmModel.hh:746
G4Element::GetIsotope
const G4Isotope * GetIsotope(G4int iso) const
Definition: G4Element.hh:169
G4VEmModel::pParticleChange
G4VParticleChange * pParticleChange
Definition: G4VEmModel.hh:398
G4VEmModel::GetName
const G4String & GetName() const
Definition: G4VEmModel.hh:753
G4VEmModel::MaxSecondaryEnergy
virtual G4double MaxSecondaryEnergy(const G4ParticleDefinition *, G4double kineticEnergy)
Definition: G4VEmModel.cc:372
G4Material::GetNumberOfElements
size_t GetNumberOfElements() const
Definition: G4Material.hh:184
G4VEmModel::xSectionTable
G4PhysicsTable * xSectionTable
Definition: G4VEmModel.hh:399
G4double
double G4double
Definition: G4Types.hh:76
G4VEmModel::SetLowEnergyLimit
void SetLowEnergyLimit(G4double)
Definition: G4VEmModel.hh:690
G4VEmModel::SetDeexcitationFlag
void SetDeexcitationFlag(G4bool val)
Definition: G4VEmModel.hh:739
G4DataVector.hh
G4VEmModel::ComputeMeanFreePath
G4double ComputeMeanFreePath(const G4ParticleDefinition *, G4double kineticEnergy, const G4Material *, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
Definition: G4VEmModel.hh:481
G4VEmModel::SetCurrentElement
void SetCurrentElement(const G4Element *)
Definition: G4VEmModel.hh:433
DBL_MAX
#define DBL_MAX
Definition: templates.hh:83
G4VEmModel::fElementData
G4ElementData * fElementData
Definition: G4VEmModel.hh:397
G4VEmModel::SelectRandomAtom
const G4Element * SelectRandomAtom(const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
Definition: G4VEmModel.hh:510
G4VEmModel::SetPolarAngleLimit
void SetPolarAngleLimit(G4double)
Definition: G4VEmModel.hh:718
G4VEmModel::GetParticleChangeForGamma
G4ParticleChangeForGamma * GetParticleChangeForGamma()
Definition: G4VEmModel.cc:121
G4VEmModel::Value
virtual G4double Value(const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy)
Definition: G4VEmModel.cc:346
G4MaterialCutsCouple::GetMaterial
const G4Material * GetMaterial() const
Definition: G4MaterialCutsCouple.hh:150
G4VEmModel::SelectRandomAtomNumber
G4int SelectRandomAtomNumber(const G4Material *)
Definition: G4VEmModel.hh:529
G4VEmModel::GetCurrentElement
const G4Element * GetCurrentElement() const
Definition: G4VEmModel.hh:440