Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/standard/include/G4ICRU73QOModel.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/G4ICRU73QOModel.hh (Version 11.3.0) and /processes/electromagnetic/standard/include/G4ICRU73QOModel.hh (Version 11.0)


  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 //                                                 26 //
 27 // -------------------------------------------     27 // -------------------------------------------------------------------
 28 //                                                 28 //
 29 // GEANT4 Class header file                        29 // GEANT4 Class header file
 30 //                                                 30 //
 31 //                                                 31 //
 32 // File name:     G4ICRU73QOModel                  32 // File name:     G4ICRU73QOModel
 33 //                                                 33 //
 34 // Author:        Alexander Bagulya                34 // Author:        Alexander Bagulya
 35 //                                                 35 //
 36 // Creation date: 21.05.2010                       36 // Creation date: 21.05.2010
 37 //                                                 37 //
 38 // Modifications:                                  38 // Modifications:
 39 //                                                 39 //
 40 //                                                 40 //
 41 // Class Description:                              41 // Class Description:
 42 //                                                 42 //
 43 // Quantum Harmonic Oscillator Model for energ     43 // Quantum Harmonic Oscillator Model for energy loss using atomic shell 
 44 // structure of atoms taking into account Q^2      44 // structure of atoms taking into account Q^2 (main for projectile charge Q), 
 45 // Q^3 and Q^4 terms for computation of energy     45 // Q^3 and Q^4 terms for computation of energy loss due to binary collisions. 
 46 // Can be applied on heavy negatively charged      46 // Can be applied on heavy negatively charged particles for the energy interval 
 47 // 10 keV - 10 MeV scaled to the proton mass.      47 // 10 keV - 10 MeV scaled to the proton mass.
 48 //                                                 48 //
 49 // Used data and formula of                        49 // Used data and formula of 
 50 // 1. G4QAOLowEnergyLoss class, S.Chauvie, P.N     50 // 1. G4QAOLowEnergyLoss class, S.Chauvie, P.Nieminen, M.G.Pia. IEEE Trans. 
 51 //    Nucl. Sci. 54 (2007) 578.                    51 //    Nucl. Sci. 54 (2007) 578.
 52 // 2. ShellStrength and ShellEnergy from ICRU'     52 // 2. ShellStrength and ShellEnergy from ICRU'73 Report 2005,
 53 // 3. Data for Ta (Z=73) from P.Sigmund, A.Shi     53 // 3. Data for Ta (Z=73) from P.Sigmund, A.Shinner. Eur. Phys. J. D15 (2001) 
 54 //    165-172                                      54 //    165-172
 55 //                                                 55 //
 56 // -------------------------------------------     56 // -------------------------------------------------------------------
 57 //                                                 57 //
 58                                                    58 
 59 #ifndef G4ICRU73QOModel_h                          59 #ifndef G4ICRU73QOModel_h
 60 #define G4ICRU73QOModel_h 1                        60 #define G4ICRU73QOModel_h 1
 61                                                    61 
 62 #include <CLHEP/Units/PhysicalConstants.h>         62 #include <CLHEP/Units/PhysicalConstants.h>
 63                                                    63 
 64 #include "G4VEmModel.hh"                           64 #include "G4VEmModel.hh"
 65 #include "G4AtomicShells.hh"                       65 #include "G4AtomicShells.hh"
 66 #include "G4DensityEffectData.hh"                  66 #include "G4DensityEffectData.hh"
 67                                                    67 
 68 class G4ParticleChangeForLoss;                     68 class G4ParticleChangeForLoss;
 69                                                    69 
 70 class G4ICRU73QOModel : public G4VEmModel          70 class G4ICRU73QOModel : public G4VEmModel
 71 {                                                  71 {
 72                                                    72 
 73 public:                                            73 public:
 74                                                    74 
 75   explicit G4ICRU73QOModel(const G4ParticleDef <<  75   explicit G4ICRU73QOModel(const G4ParticleDefinition* p = 0,
 76                            const G4String& nam     76                            const G4String& nam = "ICRU73QO");
 77                                                    77 
 78   ~G4ICRU73QOModel() = default;                    78   ~G4ICRU73QOModel() = default;
 79                                                    79 
 80   void Initialise(const G4ParticleDefinition*,     80   void Initialise(const G4ParticleDefinition*, const G4DataVector&) override;
 81                                                    81 
 82   G4double ComputeCrossSectionPerElectron(         82   G4double ComputeCrossSectionPerElectron(
 83                                  const G4Parti     83                                  const G4ParticleDefinition*,
 84                                  G4double kine     84                                  G4double kineticEnergy,
 85                                  G4double cutE     85                                  G4double cutEnergy,
 86                                  G4double maxE     86                                  G4double maxEnergy);
 87                                                    87                                  
 88   G4double ComputeCrossSectionPerAtom(             88   G4double ComputeCrossSectionPerAtom(
 89                                  const G4Parti     89                                  const G4ParticleDefinition*,
 90                                  G4double kine     90                                  G4double kineticEnergy,
 91                                  G4double Z, G     91                                  G4double Z, G4double A,
 92                                  G4double cutE     92                                  G4double cutEnergy,
 93                                  G4double maxE     93                                  G4double maxEnergy) override;
 94                                                    94                                                                   
 95   G4double CrossSectionPerVolume(const G4Mater     95   G4double CrossSectionPerVolume(const G4Material*,
 96                                  const G4Parti     96                                  const G4ParticleDefinition*,
 97                                  G4double kine     97                                  G4double kineticEnergy,
 98                                  G4double cutE     98                                  G4double cutEnergy,
 99                                  G4double maxE     99                                  G4double maxEnergy) override;
100                                                   100                                  
101   G4double ComputeDEDXPerVolume(const G4Materi    101   G4double ComputeDEDXPerVolume(const G4Material*,
102         const G4ParticleDefinition*,              102         const G4ParticleDefinition*,
103         G4double kineticEnergy,                   103         G4double kineticEnergy,
104         G4double) override;                       104         G4double) override;
105                                                   105 
106   void SampleSecondaries(std::vector<G4Dynamic    106   void SampleSecondaries(std::vector<G4DynamicParticle*>*,
107        const G4MaterialCutsCouple*,               107        const G4MaterialCutsCouple*,
108        const G4DynamicParticle*,                  108        const G4DynamicParticle*,
109        G4double tmin,                             109        G4double tmin,
110        G4double maxEnergy) override;              110        G4double maxEnergy) override;
111                                                   111 
112   // add correction to energy loss and compute    112   // add correction to energy loss and compute non-ionizing energy loss
113   void CorrectionsAlongStep(const G4MaterialCu    113   void CorrectionsAlongStep(const G4MaterialCutsCouple*,
114           const G4DynamicParticle*,               114           const G4DynamicParticle*,
115           const G4double& length,                 115           const G4double& length,
116           G4double& eloss) override;              116           G4double& eloss) override;
117                                                   117 
118   // hide assignment operator                     118   // hide assignment operator
119   G4ICRU73QOModel & operator=(const  G4ICRU73Q    119   G4ICRU73QOModel & operator=(const  G4ICRU73QOModel &right) = delete;
120   G4ICRU73QOModel(const  G4ICRU73QOModel&) = d    120   G4ICRU73QOModel(const  G4ICRU73QOModel&) = delete;
121                                                   121  
122 protected:                                        122 protected:
123                                                   123 
124   G4double MaxSecondaryEnergy(const G4Particle    124   G4double MaxSecondaryEnergy(const G4ParticleDefinition*,
125             G4double kinEnergy) final;            125             G4double kinEnergy) final;
126                                                   126 
127 private:                                          127 private:
128                                                   128 
129   inline void SetParticle(const G4ParticleDefi    129   inline void SetParticle(const G4ParticleDefinition* p);
130   inline void SetLowestKinEnergy(G4double val)    130   inline void SetLowestKinEnergy(G4double val);
131                                                   131 
132   G4double DEDX(const G4Material* material, G4    132   G4double DEDX(const G4Material* material, G4double kineticEnergy);
133                                                   133 
134   G4double DEDXPerElement(G4int Z, G4double ki    134   G4double DEDXPerElement(G4int Z, G4double kineticEnergy);
135                                                   135 
136   // get number of shell, energy and oscillato    136   // get number of shell, energy and oscillator strengths for material
137   G4int GetNumberOfShells(G4int Z) const;         137   G4int GetNumberOfShells(G4int Z) const;
138                                                   138 
139   G4double GetShellEnergy(G4int Z, G4int nbOfT    139   G4double GetShellEnergy(G4int Z, G4int nbOfTheShell) const; 
140   G4double GetOscillatorEnergy(G4int Z, G4int     140   G4double GetOscillatorEnergy(G4int Z, G4int nbOfTheShell) const; 
141   G4double GetShellStrength(G4int Z, G4int nbO    141   G4double GetShellStrength(G4int Z, G4int nbOfTheShell) const;
142                                                   142 
143   // calculate stopping number for L's term       143   // calculate stopping number for L's term
144   G4double GetL0(G4double normEnergy) const;      144   G4double GetL0(G4double normEnergy) const;
145   // terms in Z^2                                 145   // terms in Z^2
146   G4double GetL1(G4double normEnergy) const;      146   G4double GetL1(G4double normEnergy) const;
147   // terms in Z^3                                 147   // terms in Z^3
148   G4double GetL2(G4double normEnergy) const;      148   G4double GetL2(G4double normEnergy) const;
149   // terms in Z^4                                 149   // terms in Z^4
150                                                   150   
151   const G4ParticleDefinition* particle;           151   const G4ParticleDefinition* particle;
152   G4ParticleDefinition*       theElectron;        152   G4ParticleDefinition*       theElectron;   
153   G4ParticleChangeForLoss*    fParticleChange;    153   G4ParticleChangeForLoss*    fParticleChange;
154   G4DensityEffectData*        denEffData;         154   G4DensityEffectData*        denEffData;
155                                                   155 
156   G4double mass;                                  156   G4double mass;
157   G4double charge;                                157   G4double charge;
158   G4double chargeSquare;                          158   G4double chargeSquare;
159   G4double massRate;                              159   G4double massRate;
160   G4double ratio;                                 160   G4double ratio;
161   G4double lowestKinEnergy;                       161   G4double lowestKinEnergy;
162                                                   162 
163   G4bool   isInitialised;                         163   G4bool   isInitialised;
164                                                   164 
165   // Z of element at now avaliable for the mod    165   // Z of element at now avaliable for the model
166   static const G4int NQOELEM  = 26;               166   static const G4int NQOELEM  = 26;
167   static const G4int NQODATA  = 130;              167   static const G4int NQODATA  = 130;
168   static const G4int ZElementAvailable[NQOELEM    168   static const G4int ZElementAvailable[NQOELEM];
169                                                   169   
170   // number, energy and oscillator strengths      170   // number, energy and oscillator strengths
171   // for an harmonic oscillator model of mater    171   // for an harmonic oscillator model of material
172   static const G4int startElemIndex[NQOELEM];     172   static const G4int startElemIndex[NQOELEM];
173   static const G4int nbofShellsForElement[NQOE    173   static const G4int nbofShellsForElement[NQOELEM];
174   static const G4double ShellEnergy[NQODATA];     174   static const G4double ShellEnergy[NQODATA];
175   static const G4double SubShellOccupation[NQO    175   static const G4double SubShellOccupation[NQODATA];  // Z * ShellStrength
176                                                   176 
177   G4int indexZ[100];                              177   G4int indexZ[100];
178                                                   178 
179   //  variable for calculation of stopping num    179   //  variable for calculation of stopping number of L's term
180   static const G4double L0[67][2];                180   static const G4double L0[67][2];
181   static const G4double L1[22][2];                181   static const G4double L1[22][2];
182   static const G4double L2[14][2];                182   static const G4double L2[14][2];
183                                                   183   
184   G4int sizeL0;                                   184   G4int sizeL0;
185   G4int sizeL1;                                   185   G4int sizeL1;
186   G4int sizeL2;                                   186   G4int sizeL2;
187                                                   187 
188   static const G4double factorBethe[99];          188   static const G4double factorBethe[99];
189                                                   189   
190 };                                                190 };
191                                                   191 
192 //....oooOO0OOooo........oooOO0OOooo........oo    192 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
193                                                   193 
194 inline void G4ICRU73QOModel::SetParticle(const    194 inline void G4ICRU73QOModel::SetParticle(const G4ParticleDefinition* p)
195 {                                                 195 {
196   particle = p;                                   196   particle = p;
197   mass = particle->GetPDGMass();                  197   mass = particle->GetPDGMass();
198   charge = particle->GetPDGCharge()/CLHEP::epl    198   charge = particle->GetPDGCharge()/CLHEP::eplus;
199   chargeSquare = charge*charge;                   199   chargeSquare = charge*charge;
200   massRate     = mass/CLHEP::proton_mass_c2;      200   massRate     = mass/CLHEP::proton_mass_c2;
201   ratio = CLHEP::electron_mass_c2/mass;           201   ratio = CLHEP::electron_mass_c2/mass;
202 }                                                 202 }
203                                                   203 
204 //....oooOO0OOooo........oooOO0OOooo........oo    204 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
205                                                   205 
206 inline void G4ICRU73QOModel::SetLowestKinEnerg    206 inline void G4ICRU73QOModel::SetLowestKinEnergy(G4double val)
207 {                                                 207 {
208   lowestKinEnergy = val;                          208   lowestKinEnergy = val;
209 }                                                 209 }
210                                                   210 
211 #endif                                            211 #endif
212                                                   212