80 G4double rm = (m1 + m2Local) * ( std::pow(z1, .23) + std::pow(z2, .23)) ;
82 G4double er = 32.536 * m2Local * energy / ( z1 * z2 * rm ) ;
85 ionloss = 0.5*std::log(1+1.1383*er)/
86 (er+0.01312*std::pow(er,0.21226)+0.19593*std::sqrt(er)) ;
89 ionloss = 0.5*std::log(er)/er ;
94 G4double sig = 4.0 * m1 * m2Local / ((m1 + m2Local)*(m1 + m2Local)*
95 (4.0 + 0.197*std::pow(er,-1.6991)+6.584*std::pow(er,-1.0494))) ;
101 ionloss *= 8.462 * z1 * z2 * m1 / rm ;
103 if ( ionloss < 0.0) ionloss = 0.0 ;