82 G4double rm = (m1 + m2Local) * ( std::pow(z1, .23) + std::pow(z2, .23) ) ;
84 G4double er = 32.536 * m2Local * energy / ( z1 * z2 * rm ) ;
205 for (
G4int i=1; i<104; i++)
209 (a[i][1]-a[i-1][1])*(er-a[i-1][0])/(a[i][0]-a[i-1][0])+a[i-1][1];
216 G4double sig = 4.0 * m1 * m2Local / ((m1 + m2Local)*(m1 + m2Local)*
217 (4.0 + 0.197*std::pow(er,-1.6991)+6.584*std::pow(er,-1.0494))) ;
222 ionloss *= 8.462 * z1 * z2 * m1 / rm ;
224 if ( ionloss < 0.0) ionloss = 0.0 ;
std::vector< ExP01TrackerHit * > a
G4double NuclearStoppingPower(G4double kineticEnergy, G4double z1, G4double z2, G4double m1, G4double m2) const
static constexpr double keV
ThreeVector shoot(const G4int Ap, const G4int Af)