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 // =========================================== 29 // GEANT4 class header file 30 // 31 // Class: G4IonDEDXScalingICRU7 32 // 33 // Base class: G4VIonDEDXScalingAlgo 34 // 35 // Author: Anton Lechner (Anton. 36 // 37 // First implementation: 10. 05. 2009 38 // 39 // Modifications: 12. 11. 2009 - Moved all dec 40 // scaling for h 41 // Adapting Scal 42 // to changes in 43 // 44 // Class description: 45 // dE/dx scaling algorithm applied on top o 46 // covered by the ICRU 73 report) 47 // 48 // Comments: 49 // 50 // =========================================== 51 52 #ifndef G4IONDEDXSCALINGICRU73_HH 53 #define G4IONDEDXSCALINGICRU73_HH 54 55 #include "globals.hh" 56 #include "G4VIonDEDXScalingAlgorithm.hh" 57 #include "G4Material.hh" 58 #include "G4ParticleDefinition.hh" 59 #include <vector> 60 #include "G4Exp.hh" 61 62 class G4IonDEDXScalingICRU73 : public G4VIonDE 63 64 public: 65 explicit G4IonDEDXScalingICRU73(G4int minAto 66 G4int maxAtomicNumberIon = 102); 67 ~G4IonDEDXScalingICRU73(); 68 69 // Function for scaling the kinetic energy ( 70 // Returns scaling factor for a given ion. 71 G4double ScalingFactorEnergy( 72 const G4ParticleDefinition* parti 73 const G4Material* material) overr 74 75 76 // Function for scaling the dE/dx value (no 77 // Returns scaling factor for a given ion-m 78 // a given kinetic energy. 79 G4double ScalingFactorDEDX( 80 const G4ParticleDefinition* parti 81 const G4Material*, 82 G4double kineticEnergy) override; 83 84 85 // Function for defining a base particle fo 86 // (no base particle by default). Returns a 87 // particle. 88 G4int AtomicNumberBaseIon( 89 G4int atomicNumberIon, 90 const G4Material*) override; 91 92 private: 93 void UpdateCacheParticle( 94 const G4ParticleDefinition* parti 95 96 void UpdateCacheMaterial( 97 const G4Material* material); 98 99 void CreateReferenceParticles(); 100 101 G4double EquilibriumCharge( 102 G4double mass, 103 G4double charge, 104 G4double atomicNumberPow, 105 G4double kineticEnergy); 106 107 // Scaling is only applied for ions with at 108 // defined by the following parameters: 109 G4int minAtomicNumber; 110 G4int maxAtomicNumber; 111 112 G4bool referencePrepared; 113 114 // Some properties of reference particle (F 115 ///////////////////////////G4ParticleDefini 116 G4int atomicNumberRefFe; 117 G4int massNumberRefFe; 118 G4double atomicNumberRefPow23Fe; 119 G4double chargeRefFe; 120 G4double massRefFe; 121 122 // Some properties of reference particle (A 123 ///////////////////////////G4ParticleDefini 124 G4int atomicNumberRefAr; 125 G4int massNumberRefAr; 126 G4double atomicNumberRefPow23Ar; 127 G4double chargeRefAr; 128 G4double massRefAr; 129 130 // Flag indicating the use of Fe ions as re 131 G4bool useFe; 132 133 // Some properties of projectiles are store 134 const G4ParticleDefinition* cacheParticle; 135 G4int cacheMassNumber; 136 G4int cacheAtomicNumber; 137 G4double cacheAtomicNumberPow23; 138 G4double cacheCharge; 139 G4double cacheMass; 140 141 // Material pointer 142 const G4Material* cacheMaterial; 143 }; 144 145 // ########################################### 146 147 inline void G4IonDEDXScalingICRU73::UpdateCach 148 const G4ParticleDefinition* partic 149 150 if(particle != cacheParticle) { 151 152 cacheParticle = particle; 153 cacheAtomicNumber = particle -> GetAtomic 154 cacheMassNumber = particle -> GetAtomicMa 155 cacheCharge = particle -> GetPDGCharge(); 156 cacheMass = particle -> GetPDGMass(); 157 cacheAtomicNumberPow23 = std::pow(G4doubl 158 } 159 } 160 161 // ########################################### 162 163 inline void G4IonDEDXScalingICRU73::UpdateCach 164 const G4Material* material) { 165 166 if(cacheMaterial != material) { 167 168 cacheMaterial = material; 169 170 useFe = true; 171 172 size_t nmbElements = material -> GetNumbe 173 if( nmbElements > 1 ) useFe = false; 174 175 if( material -> GetName() == "G4_WATER" ) 176 } 177 } 178 179 // ########################################### 180 181 inline G4double G4IonDEDXScalingICRU73::Equili 182 G4double m 183 G4double c 184 G4double a 185 G4double k 186 187 G4double totalEnergy = kineticEnergy + mass 188 G4double betaSquared = kineticEnergy * 189 (totalEnergy + mass) / (tota 190 191 G4double beta = std::sqrt( betaSquared ); 192 193 G4double velOverBohrVel = beta / CLHEP::fine 194 195 G4double q1 = 1.0 - G4Exp(-velOverBohrVel / 196 197 return q1 * charge; 198 } 199 200 // ########################################### 201 202 #endif 203