Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/standard/include/G4eBremsstrahlungRelModel.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 /processes/electromagnetic/standard/include/G4eBremsstrahlungRelModel.hh (Version 11.3.0) and /processes/electromagnetic/standard/include/G4eBremsstrahlungRelModel.hh (Version 10.0.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 //                                                 25 //
                                                   >>  26 // $Id: G4eBremsstrahlungRelModel.hh 75582 2013-11-04 12:13:01Z gcosmo $
 26 //                                                 27 //
 27 // -------------------------------------------     28 // -------------------------------------------------------------------
 28 //                                                 29 //
 29 // GEANT4 Class header file                        30 // GEANT4 Class header file
 30 //                                                 31 //
 31 //                                                 32 //
 32 // File name:     G4eBremsstrahlungRelModel        33 // File name:     G4eBremsstrahlungRelModel
 33 //                extention of standard G4eBre     34 //                extention of standard G4eBremsstrahlungModel
 34 //                                                 35 //
 35 // Author:        Andreas Schaelicke           <<  36 // Author:        Andreas Schaelicke 
 36 //                                                 37 //
 37 // Creation date: 28.03.2008                       38 // Creation date: 28.03.2008
 38 //                                                 39 //
 39 // Modifications:                                  40 // Modifications:
 40 //                                                 41 //
 41 // 15.07.18  introduced data structures to sto << 
 42 //           dent data for faster run-time com << 
 43 //                                                 42 //
 44 // Class Description:                              43 // Class Description:
 45 //                                                 44 //
 46 // Implementation of energy loss for gamma emi     45 // Implementation of energy loss for gamma emission by electrons and
 47 // positrons including an improved version of      46 // positrons including an improved version of the LPM effect
 48                                                    47 
 49 // -------------------------------------------     48 // -------------------------------------------------------------------
 50 //                                                 49 //
 51                                                    50 
 52 #ifndef G4eBremsstrahlungRelModel_h                51 #ifndef G4eBremsstrahlungRelModel_h
 53 #define G4eBremsstrahlungRelModel_h 1              52 #define G4eBremsstrahlungRelModel_h 1
 54                                                    53 
 55 #include "G4VEmModel.hh"                           54 #include "G4VEmModel.hh"
 56 #include <memory>                              <<  55 #include "G4NistManager.hh"
                                                   >>  56 #include "G4Exp.hh"
                                                   >>  57 #include "G4Log.hh"
 57                                                    58 
 58 class G4ParticleChangeForLoss;                     59 class G4ParticleChangeForLoss;
                                                   >>  60 class G4PhysicsVector;
 59                                                    61 
 60 class G4eBremsstrahlungRelModel : public G4VEm <<  62 class G4eBremsstrahlungRelModel : public G4VEmModel
                                                   >>  63 {
 61                                                    64 
 62 public:                                            65 public:
 63                                                    66 
 64   explicit G4eBremsstrahlungRelModel(const G4P <<  67   G4eBremsstrahlungRelModel(const G4ParticleDefinition* p = 0, 
 65                                      const G4S <<  68           const G4String& nam = "eBremLPM");
 66                                                    69 
 67   ~G4eBremsstrahlungRelModel() override;       <<  70   virtual ~G4eBremsstrahlungRelModel();
 68                                                    71 
 69   void Initialise(const G4ParticleDefinition*, <<  72   virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
                                                   >>  73 
                                                   >>  74   virtual void InitialiseLocal(const G4ParticleDefinition*,
                                                   >>  75              G4VEmModel* masterModel);
                                                   >>  76 
                                                   >>  77   virtual G4double ComputeDEDXPerVolume(const G4Material*,
                                                   >>  78           const G4ParticleDefinition*,
                                                   >>  79           G4double kineticEnergy,
                                                   >>  80           G4double cutEnergy);
                                                   >>  81           
                                                   >>  82   virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition*,
                                                   >>  83                 G4double tkin, 
                                                   >>  84                 G4double Z,   G4double,
                                                   >>  85                 G4double cutEnergy,
                                                   >>  86                 G4double maxEnergy = DBL_MAX);
                                                   >>  87   
                                                   >>  88   virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
                                                   >>  89          const G4MaterialCutsCouple*,
                                                   >>  90          const G4DynamicParticle*,
                                                   >>  91          G4double cutEnergy,
                                                   >>  92          G4double maxEnergy);
                                                   >>  93 
                                                   >>  94   virtual void SetupForMaterial(const G4ParticleDefinition*,
                                                   >>  95                                 const G4Material*,G4double);
                                                   >>  96 
                                                   >>  97   virtual G4double MinPrimaryEnergy(const G4Material*,
                                                   >>  98             const G4ParticleDefinition*,
                                                   >>  99             G4double cut);
                                                   >> 100 
                                                   >> 101   inline void SetLPMconstant(G4double val);
                                                   >> 102   inline G4double LPMconstant() const;
                                                   >> 103 
                                                   >> 104   inline void SetLowestKinEnergy(G4double);
                                                   >> 105   inline G4double LowestKinEnergy() const;
 70                                                   106 
 71   void InitialiseLocal(const G4ParticleDefinit << 
 72                        G4VEmModel* masterModel << 
 73                                                << 
 74   G4double ComputeDEDXPerVolume(const G4Materi << 
 75                                 const G4Partic << 
 76                                 G4double ekin, << 
 77                                 G4double cutEn << 
 78                                                << 
 79   G4double ComputeCrossSectionPerAtom(const G4 << 
 80                                       G4double << 
 81                                       G4double << 
 82                                       G4double << 
 83                                       G4double << 
 84                                       G4double << 
 85                                                << 
 86   void SampleSecondaries(std::vector<G4Dynamic << 
 87                          const G4MaterialCutsC << 
 88                          const G4DynamicPartic << 
 89                          G4double cutEnergy,   << 
 90                          G4double maxEnergy) o << 
 91                                                << 
 92   void SetupForMaterial(const G4ParticleDefini << 
 93                         const G4Material*, G4d << 
 94                                                << 
 95   G4double MinPrimaryEnergy(const G4Material*, << 
 96                             const G4ParticleDe << 
 97                             G4double cutEnergy << 
 98                                                   107 
 99 protected:                                        108 protected:
100                                                   109 
101   virtual G4double ComputeDXSectionPerAtom(G4d    110   virtual G4double ComputeDXSectionPerAtom(G4double gammaEnergy);
102                                                   111 
103   void SetParticle(const G4ParticleDefinition* << 112   // * fast inline functions *
                                                   >> 113   inline void SetCurrentElement(const G4double);
104                                                   114 
105 private:                                          115 private:
106                                                   116 
                                                   >> 117   void InitialiseConstants();
                                                   >> 118 
                                                   >> 119   void CalcLPMFunctions(G4double gammaEnergy);
                                                   >> 120 
107   G4double ComputeBremLoss(G4double cutEnergy)    121   G4double ComputeBremLoss(G4double cutEnergy);
108                                                   122 
109   G4double ComputeXSectionPerAtom(G4double cut    123   G4double ComputeXSectionPerAtom(G4double cutEnergy);
110                                                   124 
111   G4double ComputeRelDXSectionPerAtom(G4double    125   G4double ComputeRelDXSectionPerAtom(G4double gammaEnergy);
112                                                   126 
113   // init special data per element i.e. per Z  << 127   void SetParticle(const G4ParticleDefinition* p);
114   void InitialiseElementData();                << 
115                                                   128 
116   // methods for initialisation and run-time e << 129   inline G4double Phi1(G4double,G4double);
117   void InitLPMFunctions();                     << 130   inline G4double Phi1M2(G4double,G4double);
                                                   >> 131   inline G4double Psi1(G4double,G4double);
                                                   >> 132   inline G4double Psi1M2(G4double,G4double);
                                                   >> 133 
                                                   >> 134   // hide assignment operator
                                                   >> 135   G4eBremsstrahlungRelModel & operator=(const  G4eBremsstrahlungRelModel &right);
                                                   >> 136   G4eBremsstrahlungRelModel(const  G4eBremsstrahlungRelModel&);
118                                                   137 
119   void ComputeLPMfunctions(G4double& funcXiS,  << 138 protected:
120                            G4double& funcGS,   << 
121                            G4double& funcPhiS, << 
122                            const G4double egam << 
123                                                << 
124   void GetLPMFunctions(G4double& lpmGs,        << 
125                        G4double& lpmPhis,      << 
126                        const G4double ss);     << 
127                                                << 
128   void ComputeLPMGsPhis(G4double& funcGS,      << 
129                         G4double& funcPhiS,    << 
130                         const G4double varShat << 
131                                                << 
132   // for evaluating screening related function << 
133   void ComputeScreeningFunctions(G4double& phi << 
134                                  G4double& phi << 
135                                  G4double& psi << 
136                                  G4double& psi << 
137                                  const G4doubl << 
138                                  const G4doubl << 
139   // hide assignment operator and cctr         << 
140   G4eBremsstrahlungRelModel& operator=         << 
141   (const G4eBremsstrahlungRelModel& right) = d << 
142   G4eBremsstrahlungRelModel(const  G4eBremsstr << 
143                                                   139 
144 private:                                       << 140   G4NistManager*              nist;
                                                   >> 141   const G4ParticleDefinition* particle;
                                                   >> 142   G4ParticleDefinition*       theGamma;
                                                   >> 143   G4ParticleChangeForLoss*    fParticleChange;
145                                                   144 
146   G4bool                      fIsUseCompleteSc << 145   G4double bremFactor;
147   G4bool                      fIsInitializer = << 
148   G4bool                      fUseLPM = true;  << 
149                                                   146 
150 protected:                                     << 
151                                                << 
152   G4bool                      fIsElectron = tr << 
153   G4bool                      fIsScatOffElectr << 
154   G4bool                      fIsLPMActive = f << 
155   //                                           << 
156   G4int                       fCurrentIZ = 0;  << 
157   const G4ParticleDefinition* fPrimaryParticle << 
158   G4ParticleDefinition*       fGammaParticle = << 
159   G4ParticleChangeForLoss*    fParticleChange  << 
160   // cash                                         147   // cash
161   G4double                    fPrimaryParticle << 148   G4double particleMass;
162   G4double                    fPrimaryKinEnerg << 149   G4double kinEnergy;
163   G4double                    fPrimaryTotalEne << 150   G4double totalEnergy;
164   G4double                    fDensityFactor = << 151   G4double currentZ;
165   G4double                    fDensityCorr = 0 << 152   G4double densityFactor;
166   G4double                    fLowestKinEnergy << 153   G4double densityCorr;
167   // scattering off electrons                  << 154 
168   G4double                    fNucTerm = 0.;   << 155   G4bool   isElectron;
169   G4double                    fSumTerm = 0.;   << 
170                                                   156 
171 private:                                          157 private:
172                                                   158 
173   // LPM related members                       << 159   static const G4double xgi[8], wgi[8];
174   G4double                    fLPMEnergyThresh << 160   static const G4double Fel_light[5];
175   G4double                    fLPMEnergy;      << 161   static const G4double Finel_light[5];
                                                   >> 162 
                                                   >> 163   // consts
                                                   >> 164   G4double lowestKinEnergy;
                                                   >> 165   G4double fMigdalConstant;
                                                   >> 166   G4double fLPMconstant;
                                                   >> 167   G4double energyThresholdLPM;
                                                   >> 168   G4double facFel, facFinel;
                                                   >> 169   G4double preS1,logTwo;
176                                                   170 
177 protected:                                     << 171   // cash
                                                   >> 172   G4double z13, z23, lnZ;
                                                   >> 173   G4double Fel, Finel, fCoulomb, fMax; 
178                                                   174 
179   static const G4double       gBremFactor;     << 175   // LPM effect
180   static const G4double       gMigdalConstant; << 176   G4double lpmEnergy;
                                                   >> 177   G4PhysicsVector  *fXiLPM, *fPhiLPM, *fGLPM;
                                                   >> 178   G4double xiLPM, phiLPM, gLPM;
                                                   >> 179 
                                                   >> 180   // critical gamma energies
                                                   >> 181   G4double klpm, kp;
                                                   >> 182 
                                                   >> 183   // flags
                                                   >> 184   G4bool   use_completescreening;
                                                   >> 185   G4bool   isInitialised;
                                                   >> 186 };
181                                                   187 
182 private:                                       << 188 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >> 189 
                                                   >> 190 inline void G4eBremsstrahlungRelModel::SetCurrentElement(const G4double Z)
                                                   >> 191 {
                                                   >> 192   if(Z != currentZ) {
                                                   >> 193     currentZ = Z;
                                                   >> 194 
                                                   >> 195     G4int iz = G4int(Z);
                                                   >> 196     z13 = nist->GetZ13(iz);
                                                   >> 197     z23 = z13*z13;
                                                   >> 198     lnZ = nist->GetLOGZ(iz);
                                                   >> 199 
                                                   >> 200     if (iz <= 4) {
                                                   >> 201       Fel = Fel_light[iz];  
                                                   >> 202       Finel = Finel_light[iz] ; 
                                                   >> 203     }
                                                   >> 204     else {
                                                   >> 205       Fel = facFel - lnZ/3. ;
                                                   >> 206       Finel = facFinel - 2.*lnZ/3. ;
                                                   >> 207     }
                                                   >> 208 
                                                   >> 209     fCoulomb = GetCurrentElement()->GetfCoulomb();
                                                   >> 210     fMax = Fel-fCoulomb + Finel/currentZ  +  (1.+1./currentZ)/12.;
                                                   >> 211   }
                                                   >> 212 }
                                                   >> 213 
                                                   >> 214 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
                                                   >> 215 
                                                   >> 216 
                                                   >> 217 inline G4double G4eBremsstrahlungRelModel::Phi1(G4double gg, G4double)
                                                   >> 218 {
                                                   >> 219   //       Thomas-Fermi FF from Tsai, eq.(3.38) for Z>=5
                                                   >> 220   return 20.863 - 2.*G4Log(1. + sqr(0.55846*gg) )
                                                   >> 221     - 4.*( 1. - 0.6*G4Exp(-0.9*gg) - 0.4*G4Exp(-1.5*gg) );
                                                   >> 222 }
                                                   >> 223 
                                                   >> 224 inline G4double G4eBremsstrahlungRelModel::Phi1M2(G4double gg, G4double)
                                                   >> 225 {
                                                   >> 226   //       Thomas-Fermi FF from Tsai, eq. (3.39) for Z>=5
                                                   >> 227   // return Phi1(gg,Z) - 
                                                   >> 228   return 2./(3.*(1. + 6.5*gg +6.*gg*gg) );
                                                   >> 229 }
                                                   >> 230 
                                                   >> 231 inline G4double G4eBremsstrahlungRelModel::Psi1(G4double eps, G4double)
                                                   >> 232 {
                                                   >> 233   //       Thomas-Fermi FF from Tsai, eq.(3.40) for Z>=5 
                                                   >> 234   return 28.340 - 2.*G4Log(1. + sqr(3.621*eps) )
                                                   >> 235     - 4.*( 1. - 0.7*G4Exp(-8*eps) - 0.3*G4Exp(-29.*eps) );
                                                   >> 236 }
                                                   >> 237 
                                                   >> 238 inline G4double G4eBremsstrahlungRelModel::Psi1M2(G4double eps, G4double)
                                                   >> 239 {
                                                   >> 240   //       Thomas-Fermi FF from Tsai, eq. (3.41) for Z>=5
                                                   >> 241   return  2./(3.*(1. + 40.*eps +400.*eps*eps) );
                                                   >> 242 }
                                                   >> 243 
                                                   >> 244 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
                                                   >> 245 
                                                   >> 246 inline 
                                                   >> 247 void G4eBremsstrahlungRelModel::SetLPMconstant(G4double val) 
                                                   >> 248 {
                                                   >> 249   fLPMconstant = val;
                                                   >> 250 }
                                                   >> 251 
                                                   >> 252 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
                                                   >> 253 
                                                   >> 254 inline 
                                                   >> 255 G4double G4eBremsstrahlungRelModel::LPMconstant() const 
                                                   >> 256 {
                                                   >> 257   return fLPMconstant;
                                                   >> 258 }
                                                   >> 259 
                                                   >> 260 inline void G4eBremsstrahlungRelModel::SetLowestKinEnergy(G4double val)
                                                   >> 261 {
                                                   >> 262   lowestKinEnergy = val;
                                                   >> 263 }
                                                   >> 264 
                                                   >> 265 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
                                                   >> 266 
                                                   >> 267 inline G4double G4eBremsstrahlungRelModel::LowestKinEnergy() const
                                                   >> 268 {
                                                   >> 269   return lowestKinEnergy;
                                                   >> 270 }
                                                   >> 271 
                                                   >> 272 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
183                                                   273 
184   static const G4int          gMaxZet;         << 
185   //                                           << 
186   static const G4double       gLPMconstant;    << 
187   //                                           << 
188   static const G4double       gXGL[8];         << 
189   static const G4double       gWGL[8];         << 
190   static const G4double       gFelLowZet[8];   << 
191   static const G4double       gFinelLowZet[8]; << 
192   //                                           << 
193   struct ElementData {                         << 
194     /** @brief \f$ \ln(Z) \f$  */              << 
195     G4double  fLogZ;                           << 
196     /** @brief \f$ \ln(Z)/3 + f_c \f$  */      << 
197     G4double  fFz;                             << 
198     /** @brief \f$ ((Fel-fc)+Finel*invZ)\f$  * << 
199     G4double  fZFactor1;                       << 
200     /** @brief \f$ (Fel-fc)\f$  */             << 
201     G4double  fZFactor11;                      << 
202     /** @brief \f$ (1.0+invZ)/12  \f$  */      << 
203     G4double  fZFactor2;                       << 
204     // LPM variables                           << 
205     G4double  fVarS1;                          << 
206     G4double  fILVarS1;                        << 
207     G4double  fILVarS1Cond;                    << 
208     // constant factors to the screening funct << 
209     G4double  fGammaFactor;                    << 
210     G4double  fEpsilonFactor;                  << 
211   };                                           << 
212   //                                           << 
213   struct LPMFuncs {                            << 
214     LPMFuncs() : fIsInitialized(false), fISDel << 
215     G4bool                 fIsInitialized;     << 
216     G4double               fISDelta;           << 
217     G4double               fSLimit;            << 
218     std::vector<G4double>  fLPMFuncG;          << 
219     std::vector<G4double>  fLPMFuncPhi;        << 
220   };                                           << 
221   //                                           << 
222   static std::shared_ptr<LPMFuncs> gLPMFuncs() << 
223   static std::shared_ptr<std::vector<ElementDa << 
224   std::shared_ptr<LPMFuncs> fLPMFuncs;         << 
225   std::shared_ptr<std::vector<ElementData*>> f << 
226 };                                             << 
227                                                   274 
228 #endif                                            275 #endif
229                                                   276