Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/materials/include/G4IonisParamMat.hh

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /materials/include/G4IonisParamMat.hh (Version 11.3.0) and /materials/include/G4IonisParamMat.hh (Version 10.5.p1)


  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 //
                                                   >>  26 //
 25                                                    27 
 26 // class description                               28 // class description
 27 //                                                 29 //
 28 // The class contains few (physical) quantitie     30 // The class contains few (physical) quantities related to the Ionisation
 29 // process, for a material defined by its poin     31 // process, for a material defined by its pointer G4Material*
 30 //                                                 32 //
                                                   >>  33 
                                                   >>  34 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
                                                   >>  35 
 31 // 09-07-98: data moved from G4Material (mma)      36 // 09-07-98: data moved from G4Material (mma)
 32 // 09-03-01: copy constructor and assignement      37 // 09-03-01: copy constructor and assignement operator in public (mma)
 33 // 28-10-02: add setMeanExcitationEnergy (V.Iv     38 // 28-10-02: add setMeanExcitationEnergy (V.Ivanchenko)
 34 // 27-09-07: add computation of parameters for     39 // 27-09-07: add computation of parameters for ions (V.Ivanchenko)
 35 // 04-03-08: add fBirks constant (mma)             40 // 04-03-08: add fBirks constant (mma)
 36 // 16-01-19, add exact computation of the dens <<  41 
                                                   >>  42 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
 37                                                    43 
 38 #ifndef G4IonisParamMat_HH                         44 #ifndef G4IonisParamMat_HH
 39 #define G4IonisParamMat_HH                         45 #define G4IonisParamMat_HH
 40                                                    46 
 41 #include "G4DensityEffectCalculator.hh"        << 
 42 #include "G4ios.hh"                                47 #include "G4ios.hh"
 43 #include "globals.hh"                              48 #include "globals.hh"
                                                   >>  49 #include "G4Log.hh"
                                                   >>  50 #include "G4Exp.hh"
                                                   >>  51 #include "G4Threading.hh"
 44                                                    52 
 45 class G4Material;                              <<  53 class G4Material;                        // forward declaration
 46 class G4DensityEffectData;                         54 class G4DensityEffectData;
 47                                                    55 
 48 class G4IonisParamMat                          <<  56 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
                                                   >>  57 
                                                   >>  58 class G4IonisParamMat  // with description
 49 {                                                  59 {
 50  public:                                       <<  60 public:
 51   G4IonisParamMat(const G4Material*);          <<  61 
                                                   >>  62   G4IonisParamMat(const G4Material*); 
 52   ~G4IonisParamMat();                              63   ~G4IonisParamMat();
 53   G4IonisParamMat& operator=(const G4IonisPara << 
 54   G4IonisParamMat(const G4IonisParamMat&) = de << 
 55                                                    64 
 56   // parameters for mean energy loss calculati     65   // parameters for mean energy loss calculation:
 57   inline G4double GetMeanExcitationEnergy() co <<  66   inline
                                                   >>  67   G4double  GetMeanExcitationEnergy()   const {return fMeanExcitationEnergy;};
 58                                                    68 
 59   void SetMeanExcitationEnergy(G4double value) <<  69   void      SetMeanExcitationEnergy(G4double value);
 60   G4double FindMeanExcitationEnergy(const G4Ma <<  70   G4double  FindMeanExcitationEnergy(const G4Material*) const;
 61                                                << 
 62   inline G4double GetLogMeanExcEnergy() const  << 
 63   inline G4double* GetShellCorrectionVector()  << 
 64   inline G4double GetTaul() const { return fTa << 
 65                                                    71 
                                                   >>  72   inline
                                                   >>  73   G4double  GetLogMeanExcEnergy()       const {return fLogMeanExcEnergy;};
                                                   >>  74   inline
                                                   >>  75   G4double* GetShellCorrectionVector()  const {return fShellCorrectionVector;};
                                                   >>  76   inline
                                                   >>  77   G4double  GetTaul()                   const {return fTaul;};
                                                   >>  78     
 66   // parameters of the density correction:         79   // parameters of the density correction:
 67   inline G4double GetPlasmaEnergy() const { re <<  80   inline
 68   inline G4double GetAdjustmentFactor() const  <<  81   G4double  GetPlasmaEnergy()           const {return fPlasmaEnergy;};
 69   inline G4double GetCdensity() const { return <<  82   inline
 70   inline G4double GetMdensity() const { return <<  83   G4double  GetAdjustmentFactor()       const {return fAdjustmentFactor;};
 71   inline G4double GetAdensity() const { return <<  84   inline
 72   inline G4double GetX0density() const { retur <<  85   G4double  GetCdensity()               const {return fCdensity;};
 73   inline G4double GetX1density() const { retur <<  86   inline
 74   inline G4double GetD0density() const { retur <<  87   G4double  GetMdensity()               const {return fMdensity;};
                                                   >>  88   inline
                                                   >>  89   G4double  GetAdensity()               const {return fAdensity;};
                                                   >>  90   inline
                                                   >>  91   G4double  GetX0density()              const {return fX0density;};
                                                   >>  92   inline
                                                   >>  93   G4double  GetX1density()              const {return fX1density;};
                                                   >>  94   inline
                                                   >>  95   G4double  GetD0density()              const {return fD0density;};
 75                                                    96 
 76   // user defined density correction parameter     97   // user defined density correction parameterisation
 77   void SetDensityEffectParameters(             <<  98   void SetDensityEffectParameters(G4double cd, G4double md, G4double ad,
 78     G4double cd, G4double md, G4double ad, G4d <<  99                                   G4double x0, G4double x1, G4double d0);
 79                                                   100 
 80   // defined density correction parameterisati    101   // defined density correction parameterisation via base material
 81   void SetDensityEffectParameters(const G4Mate    102   void SetDensityEffectParameters(const G4Material* bmat);
 82                                                << 103     
 83   void ComputeDensityEffectOnFly(G4bool);      << 
 84                                                << 
 85   inline G4DensityEffectCalculator* GetDensity << 
 86   {                                            << 
 87     return fDensityEffectCalc;                 << 
 88   }                                            << 
 89                                                << 
 90   // compute density correction as a function     104   // compute density correction as a function of the kinematic variable
 91   // x = log10(beta*gamma) using parameterisat << 105   // x = log10(beta*gamma)  
 92   inline G4double DensityCorrection(G4double x << 106   inline G4double DensityCorrection(G4double x);
 93   {                                            << 
 94     return (nullptr == fDensityEffectCalc) ? G << 
 95                                            : f << 
 96   }                                            << 
 97                                                << 
 98   // use parameterisation                      << 
 99   G4double GetDensityCorrection(G4double x) co << 
100                                                   107 
101   static G4DensityEffectData* GetDensityEffect    108   static G4DensityEffectData* GetDensityEffectData();
102                                                   109 
103   // parameters of the energy loss fluctuation    110   // parameters of the energy loss fluctuation model:
104   inline G4double GetF1fluct() const { return  << 111   inline
105   inline G4double GetF2fluct() const { return  << 112   G4double  GetF1fluct()                const {return fF1fluct;};
106   inline G4double GetEnergy1fluct() const { re << 113   inline
107   inline G4double GetLogEnergy1fluct() const { << 114   G4double  GetF2fluct()                const {return fF2fluct;};
108   inline G4double GetEnergy2fluct() const { re << 115   inline
109   inline G4double GetLogEnergy2fluct() const { << 116   G4double  GetEnergy1fluct()           const {return fEnergy1fluct;};
110   inline G4double GetEnergy0fluct() const { re << 117   inline
111   inline G4double GetRateionexcfluct() const { << 118   G4double  GetLogEnergy1fluct()        const {return fLogEnergy1fluct;};
                                                   >> 119   inline
                                                   >> 120   G4double  GetEnergy2fluct()           const {return fEnergy2fluct;};
                                                   >> 121   inline
                                                   >> 122   G4double  GetLogEnergy2fluct()        const {return fLogEnergy2fluct;};
                                                   >> 123   inline
                                                   >> 124   G4double  GetEnergy0fluct()           const {return fEnergy0fluct;};
                                                   >> 125   inline
                                                   >> 126   G4double  GetRateionexcfluct()        const {return fRateionexcfluct;};
112                                                   127 
113   // parameters for ion corrections computatio    128   // parameters for ion corrections computations
114   inline G4double GetZeffective() const { retu << 129   inline
115   inline G4double GetFermiEnergy() const { ret << 130   G4double  GetZeffective()             const {return fZeff;};
116   inline G4double GetLFactor() const { return  << 131   inline
117   inline G4double GetInvA23() const { return f << 132   G4double  GetFermiEnergy()            const {return fFermiEnergy;};
118                                                << 133   inline
                                                   >> 134   G4double  GetLFactor()                const {return fLfactor;};
                                                   >> 135   inline
                                                   >> 136   G4double  GetInvA23()                 const {return fInvA23;};
                                                   >> 137     
119   // parameters for Birks attenuation:            138   // parameters for Birks attenuation:
120   inline void SetBirksConstant(G4double value) << 139   inline
121   inline G4double GetBirksConstant() const { r << 140   void      SetBirksConstant(G4double value) {fBirks = value;}; 
                                                   >> 141   inline
                                                   >> 142   G4double  GetBirksConstant()         const {return fBirks;};
                                                   >> 143 
                                                   >> 144   // parameters for average energy per ion 
                                                   >> 145   inline
                                                   >> 146   void      SetMeanEnergyPerIonPair(G4double value) {fMeanEnergyPerIon = value;}; 
                                                   >> 147   inline
                                                   >> 148   G4double  GetMeanEnergyPerIonPair()         const {return fMeanEnergyPerIon;};
                                                   >> 149       
                                                   >> 150   G4IonisParamMat(__void__&);
                                                   >> 151   // Fake default constructor for usage restricted to direct object
                                                   >> 152   // persistency for clients requiring preallocation of memory for
                                                   >> 153   // persistifiable objects.
122                                                   154 
123   // parameters for average energy per ion     << 155 private:
124   inline void SetMeanEnergyPerIonPair(G4double << 156     
125   inline G4double GetMeanEnergyPerIonPair() co << 
126                                                << 
127   // parameter for sampling of positron annihi << 
128   inline void SetOrtoPositroniumFraction(G4dou << 
129   inline G4double GetOrtoPositroniumFraction() << 
130                                                << 
131                                                << 
132   // operators                                 << 
133   G4bool operator==(const G4IonisParamMat&) co << 
134   G4bool operator!=(const G4IonisParamMat&) co << 
135                                                << 
136  private:                                      << 
137   // Compute mean parameters : ExcitationEnerg    157   // Compute mean parameters : ExcitationEnergy,Shell corretion vector ...
138   void ComputeMeanParameters();                   158   void ComputeMeanParameters();
139                                                   159 
140   // Compute parameters for the density effect    160   // Compute parameters for the density effect
141   void ComputeDensityEffectParameters();       << 161   void ComputeDensityEffect();
142                                                   162 
143   // Compute parameters for the energy fluctua    163   // Compute parameters for the energy fluctuation model
144   void ComputeFluctModel();                       164   void ComputeFluctModel();
145                                                   165 
146   // Compute parameters for ion parameterizati    166   // Compute parameters for ion parameterizations
147   void ComputeIonParameters();                    167   void ComputeIonParameters();
148                                                   168 
                                                   >> 169   // operators
                                                   >> 170   G4IonisParamMat& operator=(const G4IonisParamMat&) = delete;
                                                   >> 171   G4bool operator==(const G4IonisParamMat&) const = delete;
                                                   >> 172   G4bool operator!=(const G4IonisParamMat&) const = delete;
                                                   >> 173   G4IonisParamMat(const G4IonisParamMat&) = delete;
                                                   >> 174 
149   //                                              175   //
150   // data members                                 176   // data members
151   //                                              177   //
152   const G4Material* fMaterial;  // this materi << 178   const G4Material* fMaterial;       // this material
153                                                << 
154   G4DensityEffectCalculator* fDensityEffectCal << 
155   G4double* fShellCorrectionVector;  // shell  << 
156                                                   179 
157   // parameters for mean energy loss calculati    180   // parameters for mean energy loss calculation
158   G4double fMeanExcitationEnergy;  //          << 181   G4double  fMeanExcitationEnergy;          // 
159   G4double fLogMeanExcEnergy;  //              << 182   G4double  fLogMeanExcEnergy;              // 
160   G4double fTaul;  // lower limit of Bethe-Blo << 183   G4double* fShellCorrectionVector;         // shell correction coefficients
                                                   >> 184   G4double  fTaul;                          // lower limit of Bethe-Bloch formula
161                                                   185 
162   // parameters of the density correction         186   // parameters of the density correction
163   G4double fCdensity;  // mat.constant         << 187   G4double fCdensity;                      // mat.constant
164   G4double fMdensity;  // exponent             << 188   G4double fMdensity;                      // exponent
165   G4double fAdensity;  //                      << 189   G4double fAdensity;                      //
166   G4double fX0density;  //                     << 190   G4double fX0density;                     //
167   G4double fX1density;  //                     << 191   G4double fX1density;                     //
168   G4double fD0density;                            192   G4double fD0density;
169                                                   193 
170   G4double fPlasmaEnergy;                         194   G4double fPlasmaEnergy;
171   G4double fAdjustmentFactor;                     195   G4double fAdjustmentFactor;
172                                                   196 
173   // parameters of the energy loss fluctuation    197   // parameters of the energy loss fluctuation model
174   G4double fF1fluct;                           << 198   G4double fF1fluct;                       
175   G4double fF2fluct;                           << 199   G4double fF2fluct;                       
176   G4double fEnergy1fluct;                         200   G4double fEnergy1fluct;
177   G4double fLogEnergy1fluct;                      201   G4double fLogEnergy1fluct;
178   G4double fEnergy2fluct;                         202   G4double fEnergy2fluct;
179   G4double fLogEnergy2fluct;                      203   G4double fLogEnergy2fluct;
180   G4double fEnergy0fluct;                         204   G4double fEnergy0fluct;
181   G4double fRateionexcfluct;                      205   G4double fRateionexcfluct;
182                                                   206 
183   // parameters for ion corrections computatio    207   // parameters for ion corrections computations
184   G4double fZeff;                                 208   G4double fZeff;
185   G4double fFermiEnergy;                          209   G4double fFermiEnergy;
186   G4double fLfactor;                              210   G4double fLfactor;
187   G4double fInvA23;                               211   G4double fInvA23;
188                                                << 212     
189   // parameter for Birks attenuation              213   // parameter for Birks attenuation
190   G4double fBirks;                                214   G4double fBirks;
191   // average energy per ion pair                  215   // average energy per ion pair
192   G4double fMeanEnergyPerIon;                     216   G4double fMeanEnergyPerIon;
193   G4double twoln10;                            << 
194   // parameter for sampling of positron annihi << 
195   G4double fOrtoPositroniumFraction{0.035};    << 
196                                                   217 
197   // static data created only once                218   // static data created only once
198   static G4DensityEffectData* fDensityData;       219   static G4DensityEffectData* fDensityData;
                                                   >> 220   G4double twoln10;
                                                   >> 221 #ifdef G4MULTITHREADED
                                                   >> 222   static G4Mutex ionisMutex;
                                                   >> 223 #endif
199 };                                                224 };
                                                   >> 225 
                                                   >> 226 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... ....oooOO0OOooo....
                                                   >> 227 
                                                   >> 228 inline G4double G4IonisParamMat::DensityCorrection(G4double x)
                                                   >> 229 {
                                                   >> 230   // x = log10(beta*gamma)  
                                                   >> 231   G4double y = 0.0;
                                                   >> 232   if(x < fX0density) {
                                                   >> 233     if(fD0density > 0.0) { y = fD0density*G4Exp(twoln10*(x - fX0density)); }
                                                   >> 234   } else if(x >= fX1density) { y = twoln10*x - fCdensity; }
                                                   >> 235   else {y = twoln10*x - fCdensity + fAdensity*G4Exp(G4Log(fX1density - x)*fMdensity);}
                                                   >> 236   return y;
                                                   >> 237 }
200                                                   238 
201 #endif                                            239 #endif
202                                                   240