Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // 27 // ------------------------------------------- 28 // GEANT 4 class header file 29 // 30 // History: New Implementation 31 // 32 // ---------- G4QAOLowEnergyLoss physics 33 // by Stephane Chauvie, 21 Ma 34 // 35 // Modified: 36 // 16/09/2000 S. Chauvie Oscillator for all m 37 // 23/05/2000 MGP Made compliant to design 38 // 01/06/2001 V.Ivanchenko replace names by Z 39 // 40 // Class description: 41 // Quantal Harmonic Oscillator Model for energ 42 // Further documentation available from http:/ 43 44 // ------------------------------------------- 45 46 47 #ifndef G4QAOLowEnergyLoss_hh 48 #define G4QAOLowEnergyLoss_hh 1 49 50 #include "G4VLowEnergyModel.hh" 51 #include "globals.hh" 52 53 class G4QAOLowEnergyLoss : public G4VLowEnergy 54 { 55 public: 56 explicit G4QAOLowEnergyLoss(const G4String& 57 ~G4QAOLowEnergyLoss(); 58 59 G4double HighEnergyLimit(const G4ParticleDef 60 const G4Material* m 61 // returns the higher limit for model validi 62 63 G4double LowEnergyLimit(const G4ParticleDefi 64 const G4Material* ma 65 // returns the lower limit for model validit 66 67 G4double HighEnergyLimit(const G4ParticleDef 68 // returns the higher limit for model validi 69 70 G4double LowEnergyLimit(const G4ParticleDefi 71 // returns the lower limit for model validit 72 73 G4bool IsInCharge(const G4DynamicParticle* p 74 const G4Material* material) const over 75 // returns true if the model is applicable a 76 // that particle for that material 77 78 G4bool IsInCharge(const G4ParticleDefinition 79 const G4Material* material) const over 80 // returns true if the model is applicable a 81 // that particle for that material 82 83 G4double TheValue(const G4DynamicParticle* p 84 const G4Material* material) ove 85 // returns the energy loss via the quantal h 86 87 G4double TheValue(const G4ParticleDefinition 88 const G4Material* material 89 G4double kinet 90 // returns the energy loss via the quantal h 91 92 private: 93 G4double EnergyLoss(const G4Material* materi 94 G4double kineticEn 95 G4double zParticle 96 // returns the energy loss via the quantal h 97 98 // get number of shell, energy and oscillato 99 G4int GetNumberOfShell(const G4Material* mat 100 101 G4double GetShellEnergy(const G4Material* ma 102 G4double GetOscillatorEnergy(const G4Materia 103 G4double GetShellStrength(const G4Material* 104 G4double GetOccupationNumber(G4int Z, G4int 105 106 // calculate stopping number for L's term 107 G4double GetL0(G4double normEnergy) const; 108 // terms in Z^2 109 G4double GetL1(G4double normEnergy) const; 110 // terms in Z^3 111 G4double GetL2(G4double normEnergy) const; 112 // terms in Z^4 113 114 // number, energy and oscillator strengths 115 // for an harmonic oscillator model of mater 116 static const G4int nbofShellForMaterial[6]; 117 static const G4double alShellEnergy[3]; 118 static const G4double alShellStrength[3]; 119 static const G4double siShellEnergy[3]; 120 static const G4double siShellStrength[3]; 121 static const G4double cuShellEnergy[4]; 122 static const G4double cuShellStrength[4]; 123 static const G4double taShellEnergy[6]; 124 static const G4double taShellStrength[6]; 125 static const G4double auShellEnergy[6]; 126 static const G4double auShellStrength[6]; 127 static const G4double ptShellEnergy[6]; 128 static const G4double ptShellStrength[6]; 129 // variable for calculation of stopping num 130 static const G4double L0[67][2]; 131 static const G4double L1[22][2]; 132 static const G4double L2[14][2]; 133 static const G4int nbOfElectronPerSubShell[1 134 static const G4int fNumberOfShells[101]; 135 136 // Z of element at now avaliable for the mod 137 static const G4int materialAvailable[6]; 138 139 G4int numberOfMaterials; 140 G4int sizeL0; 141 G4int sizeL1; 142 G4int sizeL2; 143 144 }; 145 146 #endif 147