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.1.p2)


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