126 {0.2865,0.1266,-0.001429,0.02402,-0.01135,0.001475};
131 for (
G4int i=1; i<6; ++i) {
136 if(x < 0.2) { ex = x * (1 - 0.5*
x); }
137 else { ex = 1. -
G4Exp(-x); }
142 if(tq2 < 0.2) { tt *= (1.0 - tq2 + 0.5*tq2*tq2); }
143 else { tt *=
G4Exp(-tq2); }
145 effCharge = charge*(1.0 + tt) * std::sqrt(ex);
162 y = vF * std::sqrt(v1sq) * ( 1.0 + 0.2/v1sq ) / zi23 ;
166 y = 0.692308 * vF * (1.0 + 0.666666*v1sq + v1sq*v1sq/15.0) / zi23 ;
172 q = 1.0 -
G4Exp( 0.803*y3 - 1.3167*y3*y3 - 0.38157*y - 0.008983*y*y);
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4double A13(G4double A) const
static constexpr double MeV
static constexpr double keV
G4double A23(G4double A) const
G4double chargeCorrection
const G4Material * lastMat
G4double GetPDGCharge() const
G4double GetPDGMass() const
G4double G4Log(G4double x)
virtual ~G4ionEffectiveCharge()
static constexpr double proton_mass_c2
static constexpr double amu_c2
static G4Pow * GetInstance()
G4double GetFermiEnergy() const
static constexpr double eplus
G4IonisParamMat * GetIonisation() const
G4double GetZeffective() const
G4double EffectiveCharge(const G4ParticleDefinition *p, const G4Material *material, G4double kineticEnergy)
const G4ParticleDefinition * lastPart
static constexpr double eplus