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 ]

  1 //
  2 // ********************************************************************
  3 // * License and Disclaimer                                           *
  4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.                             *
 10 // *                                                                  *
 11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                                                  *
 18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // ********************************************************************
 25 //
 26 //
 27 // -------------------------------------------------------------------
 28 //
 29 // GEANT4 Class header file
 30 //
 31 //
 32 // File name:     G4ICRU73QOModel
 33 //
 34 // Author:        Alexander Bagulya
 35 //
 36 // Creation date: 21.05.2010
 37 //
 38 // Modifications:
 39 //
 40 //
 41 // Class Description:
 42 //
 43 // Quantum Harmonic Oscillator Model for energy loss using atomic shell 
 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 loss due to binary collisions. 
 46 // Can be applied on heavy negatively charged particles for the energy interval 
 47 // 10 keV - 10 MeV scaled to the proton mass.
 48 //
 49 // Used data and formula of 
 50 // 1. G4QAOLowEnergyLoss class, S.Chauvie, P.Nieminen, M.G.Pia. IEEE Trans. 
 51 //    Nucl. Sci. 54 (2007) 578.
 52 // 2. ShellStrength and ShellEnergy from ICRU'73 Report 2005,
 53 // 3. Data for Ta (Z=73) from P.Sigmund, A.Shinner. Eur. Phys. J. D15 (2001) 
 54 //    165-172
 55 //
 56 // -------------------------------------------------------------------
 57 //
 58 
 59 #ifndef G4ICRU73QOModel_h
 60 #define G4ICRU73QOModel_h 1
 61 
 62 #include <CLHEP/Units/PhysicalConstants.h>
 63 
 64 #include "G4VEmModel.hh"
 65 #include "G4AtomicShells.hh"
 66 #include "G4DensityEffectData.hh"
 67 
 68 class G4ParticleChangeForLoss;
 69 
 70 class G4ICRU73QOModel : public G4VEmModel
 71 {
 72 
 73 public:
 74 
 75   explicit G4ICRU73QOModel(const G4ParticleDefinition* p = nullptr,
 76                            const G4String& nam = "ICRU73QO");
 77 
 78   ~G4ICRU73QOModel() = default;
 79 
 80   void Initialise(const G4ParticleDefinition*, const G4DataVector&) override;
 81 
 82   G4double ComputeCrossSectionPerElectron(
 83                                  const G4ParticleDefinition*,
 84                                  G4double kineticEnergy,
 85                                  G4double cutEnergy,
 86                                  G4double maxEnergy);
 87                                  
 88   G4double ComputeCrossSectionPerAtom(
 89                                  const G4ParticleDefinition*,
 90                                  G4double kineticEnergy,
 91                                  G4double Z, G4double A,
 92                                  G4double cutEnergy,
 93                                  G4double maxEnergy) override;
 94                                                                   
 95   G4double CrossSectionPerVolume(const G4Material*,
 96                                  const G4ParticleDefinition*,
 97                                  G4double kineticEnergy,
 98                                  G4double cutEnergy,
 99                                  G4double maxEnergy) override;
100                                  
101   G4double ComputeDEDXPerVolume(const G4Material*,
102         const G4ParticleDefinition*,
103         G4double kineticEnergy,
104         G4double) override;
105 
106   void SampleSecondaries(std::vector<G4DynamicParticle*>*,
107        const G4MaterialCutsCouple*,
108        const G4DynamicParticle*,
109        G4double tmin,
110        G4double maxEnergy) override;
111 
112   // add correction to energy loss and compute non-ionizing energy loss
113   void CorrectionsAlongStep(const G4MaterialCutsCouple*,
114           const G4DynamicParticle*,
115           const G4double& length,
116           G4double& eloss) override;
117 
118   // hide assignment operator
119   G4ICRU73QOModel & operator=(const  G4ICRU73QOModel &right) = delete;
120   G4ICRU73QOModel(const  G4ICRU73QOModel&) = delete;
121  
122 protected:
123 
124   G4double MaxSecondaryEnergy(const G4ParticleDefinition*,
125             G4double kinEnergy) final;
126 
127 private:
128 
129   inline void SetParticle(const G4ParticleDefinition* p);
130   inline void SetLowestKinEnergy(G4double val);
131 
132   G4double DEDX(const G4Material* material, G4double kineticEnergy);
133 
134   G4double DEDXPerElement(G4int Z, G4double kineticEnergy);
135 
136   // get number of shell, energy and oscillator strengths for material
137   G4int GetNumberOfShells(G4int Z) const;
138 
139   G4double GetShellEnergy(G4int Z, G4int nbOfTheShell) const; 
140   G4double GetOscillatorEnergy(G4int Z, G4int nbOfTheShell) const; 
141   G4double GetShellStrength(G4int Z, G4int nbOfTheShell) const;
142 
143   // calculate stopping number for L's term
144   G4double GetL0(G4double normEnergy) const;
145   // terms in Z^2
146   G4double GetL1(G4double normEnergy) const;
147   // terms in Z^3
148   G4double GetL2(G4double normEnergy) const;
149   // terms in Z^4
150   
151   const G4ParticleDefinition* particle;
152   G4ParticleDefinition*       theElectron;   
153   G4ParticleChangeForLoss*    fParticleChange;
154   G4DensityEffectData*        denEffData;
155 
156   G4double mass;
157   G4double charge;
158   G4double chargeSquare;
159   G4double massRate;
160   G4double ratio;
161   G4double lowestKinEnergy;
162 
163   G4bool   isInitialised;
164 
165   // Z of element at now avaliable for the model
166   static const G4int NQOELEM  = 26;
167   static const G4int NQODATA  = 130;
168   static const G4int ZElementAvailable[NQOELEM];
169   
170   // number, energy and oscillator strengths
171   // for an harmonic oscillator model of material
172   static const G4int startElemIndex[NQOELEM];
173   static const G4int nbofShellsForElement[NQOELEM];
174   static const G4double ShellEnergy[NQODATA];
175   static const G4double SubShellOccupation[NQODATA];  // Z * ShellStrength
176 
177   G4int indexZ[100];
178 
179   //  variable for calculation of stopping number of L's term
180   static const G4double L0[67][2];
181   static const G4double L1[22][2];
182   static const G4double L2[14][2];
183   
184   G4int sizeL0;
185   G4int sizeL1;
186   G4int sizeL2;
187 
188   static const G4double factorBethe[99];
189   
190 };
191 
192 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
193 
194 inline void G4ICRU73QOModel::SetParticle(const G4ParticleDefinition* p)
195 {
196   particle = p;
197   mass = particle->GetPDGMass();
198   charge = particle->GetPDGCharge()/CLHEP::eplus;
199   chargeSquare = charge*charge;
200   massRate     = mass/CLHEP::proton_mass_c2;
201   ratio = CLHEP::electron_mass_c2/mass;
202 }
203 
204 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
205 
206 inline void G4ICRU73QOModel::SetLowestKinEnergy(G4double val)
207 {
208   lowestKinEnergy = val;
209 }
210 
211 #endif
212