include/Geant4/G4ICRU73QOModel.hh

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 //
0027 // -------------------------------------------------------------------
0028 //
0029 // GEANT4 Class header file
0030 //
0031 //
0032 // File name:     G4ICRU73QOModel
0033 //
0034 // Author:        Alexander Bagulya
0035 //
0036 // Creation date: 21.05.2010
0037 //
0038 // Modifications:
0039 //
0040 //
0041 // Class Description:
0042 //
0043 // Quantum Harmonic Oscillator Model for energy loss using atomic shell
0044 // structure of atoms taking into account Q^2 (main for projectile charge Q),
0045 // Q^3 and Q^4 terms for computation of energy loss due to binary collisions.
0046 // Can be applied on heavy negatively charged particles for the energy interval
0047 // 10 keV - 10 MeV scaled to the proton mass.
0048 //
0049 // Used data and formula of
0050 // 1. G4QAOLowEnergyLoss class, S.Chauvie, P.Nieminen, M.G.Pia. IEEE Trans.
0051 //    Nucl. Sci. 54 (2007) 578.
0052 // 2. ShellStrength and ShellEnergy from ICRU'73 Report 2005,
0053 // 3. Data for Ta (Z=73) from P.Sigmund, A.Shinner. Eur. Phys. J. D15 (2001)
0054 //    165-172
0055 //
0056 // -------------------------------------------------------------------
0057 //
0058
0059 #ifndef G4ICRU73QOModel_h
0060 #define G4ICRU73QOModel_h 1
0061
0062 #include <CLHEP/Units/PhysicalConstants.h>
0063
0064 #include "G4VEmModel.hh"
0065 #include "G4AtomicShells.hh"
0066 #include "G4DensityEffectData.hh"
0067
0068 class G4ParticleChangeForLoss;
0069
0070 class G4ICRU73QOModel : public G4VEmModel
0071 {
0072
0073 public:
0074
0075   explicit G4ICRU73QOModel(const G4ParticleDefinition* p = nullptr,
0076                            const G4String& nam = "ICRU73QO");
0077
0078   ~G4ICRU73QOModel() = default;
0079
0080   void Initialise(const G4ParticleDefinition*, const G4DataVector&) override;
0081
0082   G4double ComputeCrossSectionPerElectron(
0083                                  const G4ParticleDefinition*,
0084                                  G4double kineticEnergy,
0085                                  G4double cutEnergy,
0086                                  G4double maxEnergy);
0087
0088   G4double ComputeCrossSectionPerAtom(
0089                                  const G4ParticleDefinition*,
0090                                  G4double kineticEnergy,
0091                                  G4double Z, G4double A,
0092                                  G4double cutEnergy,
0093                                  G4double maxEnergy) override;
0094
0095   G4double CrossSectionPerVolume(const G4Material*,
0096                                  const G4ParticleDefinition*,
0097                                  G4double kineticEnergy,
0098                                  G4double cutEnergy,
0099                                  G4double maxEnergy) override;
0100
0101   G4double ComputeDEDXPerVolume(const G4Material*,
0102                 const G4ParticleDefinition*,
0103                 G4double kineticEnergy,
0104                 G4double) override;
0105
0106   void SampleSecondaries(std::vector<G4DynamicParticle*>*,
0107              const G4MaterialCutsCouple*,
0108              const G4DynamicParticle*,
0109              G4double tmin,
0110              G4double maxEnergy) override;
0111
0112   // hide assignment operator
0113   G4ICRU73QOModel & operator=(const  G4ICRU73QOModel &right) = delete;
0114   G4ICRU73QOModel(const  G4ICRU73QOModel&) = delete;
0115
0116 protected:
0117
0118   G4double MaxSecondaryEnergy(const G4ParticleDefinition*,
0119                   G4double kinEnergy) final;
0120
0121 private:
0122
0123   inline void SetParticle(const G4ParticleDefinition* p);
0124   inline void SetLowestKinEnergy(G4double val);
0125
0126   G4double DEDX(const G4Material* material, G4double kineticEnergy);
0127
0128   G4double DEDXPerElement(G4int Z, G4double kineticEnergy);
0129
0130   // get number of shell, energy and oscillator strengths for material
0131   G4int GetNumberOfShells(G4int Z) const;
0132
0133   G4double GetShellEnergy(G4int Z, G4int nbOfTheShell) const;
0134   G4double GetOscillatorEnergy(G4int Z, G4int nbOfTheShell) const;
0135   G4double GetShellStrength(G4int Z, G4int nbOfTheShell) const;
0136
0137   // calculate stopping number for L's term
0138   G4double GetL0(G4double normEnergy) const;
0139   // terms in Z^2
0140   G4double GetL1(G4double normEnergy) const;
0141   // terms in Z^3
0142   G4double GetL2(G4double normEnergy) const;
0143   // terms in Z^4
0144
0145   const G4ParticleDefinition* particle;
0146   G4ParticleDefinition*       theElectron;
0147   G4ParticleChangeForLoss*    fParticleChange;
0148   G4DensityEffectData*        denEffData;
0149
0150   G4double mass;
0151   G4double charge;
0152   G4double chargeSquare;
0153   G4double massRate;
0154   G4double ratio;
0155   G4double lowestKinEnergy;
0156
0157   G4bool   isInitialised;
0158
0159   // Z of element at now avaliable for the model
0160   static const G4int NQOELEM  = 26;
0161   static const G4int NQODATA  = 130;
0162   static const G4int ZElementAvailable[NQOELEM];
0163
0164   // number, energy and oscillator strengths
0165   // for an harmonic oscillator model of material
0166   static const G4int startElemIndex[NQOELEM];
0167   static const G4int nbofShellsForElement[NQOELEM];
0168   static const G4double ShellEnergy[NQODATA];
0169   static const G4double SubShellOccupation[NQODATA];  // Z * ShellStrength
0170
0171   G4int indexZ[100];
0172
0173   //  variable for calculation of stopping number of L's term
0174   static const G4double L0[67][2];
0175   static const G4double L1[22][2];
0176   static const G4double L2[14][2];
0177
0178   G4int sizeL0;
0179   G4int sizeL1;
0180   G4int sizeL2;
0181
0182   static const G4double factorBethe[99];
0183
0184 };
0185
0186 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0187
0188 inline void G4ICRU73QOModel::SetParticle(const G4ParticleDefinition* p)
0189 {
0190   particle = p;
0191   mass = particle->GetPDGMass();
0192   charge = particle->GetPDGCharge()/CLHEP::eplus;
0193   chargeSquare = charge*charge;
0194   massRate     = mass/CLHEP::proton_mass_c2;
0195   ratio = CLHEP::electron_mass_c2/mass;
0196 }
0197
0198 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0199
0200 inline void G4ICRU73QOModel::SetLowestKinEnergy(G4double val)
0201 {
0202   lowestKinEnergy = val;
0203 }
0204
0205 #endif