52 #ifndef G4IONDEDXSCALINGICRU73_HH
53 #define G4IONDEDXSCALINGICRU73_HH
66 G4int maxAtomicNumberIon = 102);
89 G4int atomicNumberIon,
172 size_t nmbElements = material -> GetNumberOfElements();
173 if( nmbElements > 1 )
useFe =
false;
175 if( material -> GetName() ==
"G4_WATER" )
useFe =
true;
187 G4double totalEnergy = kineticEnergy + mass;
188 G4double betaSquared = kineticEnergy *
189 (totalEnergy + mass) / (totalEnergy * totalEnergy);
191 G4double beta = std::sqrt( betaSquared );
193 G4double velOverBohrVel = beta / CLHEP::fine_structure_const;
195 G4double q1 = 1.0 - std::exp(-velOverBohrVel / atomicNumberPow);
void UpdateCacheParticle(const G4ParticleDefinition *particle)
G4double ScalingFactorDEDX(const G4ParticleDefinition *particle, const G4Material *, G4double kineticEnergy)
const G4Material * cacheMaterial
G4double ScalingFactorEnergy(const G4ParticleDefinition *particle, const G4Material *material)
G4int AtomicNumberBaseIon(G4int atomicNumberIon, const G4Material *)
G4double EquilibriumCharge(G4double mass, G4double charge, G4double atomicNumberPow, G4double kineticEnergy)
~G4IonDEDXScalingICRU73()
void CreateReferenceParticles()
const G4ParticleDefinition * cacheParticle
G4double atomicNumberRefPow23Ar
void UpdateCacheMaterial(const G4Material *material)
G4double cacheAtomicNumberPow23
G4IonDEDXScalingICRU73(G4int minAtomicNumberIon=19, G4int maxAtomicNumberIon=102)
G4double atomicNumberRefPow23Fe