65 lowEnergyLimit(1.*
MeV),
66 highEnergyLimit(100.*
GeV),
67 twoln10(2.*std::log(10.)),
167 G4double tau = kineticEnergy/particleMass ;
171 if ( tau < taul ) tau = taul ;
188 bg2 = tau*(tau+2.0) ;
189 beta2 = bg2/(gamma*gamma) ;
201 if ( x < x1den ) delta += aden*std::pow((x1den-x),mden) ;
211 sh += shellCorrectionVector[
k]/
x;
217 sh += shellCorrectionVector[
k]/
x;
219 sh *= std::log(tau/taul)/std::log(
taulim/taul) ;
224 ionloss -= delta + sh ;
227 if ( ionloss < 0.0) ionloss = 0.0 ;