48 if (ZRes > ARes || ARes < 1) {
49 G4cout <<
"G4GEMCoulombBarrierHE::GetCoulombBarrier: "
50 <<
"Wrong values for "
51 <<
"residual nucleus A = " << ARes <<
" "
52 <<
"and residual nucleus Z = " << ZRes <<
G4endl;
58 G4double CompoundRadius = CalcCompoundRadius(ARes);
65 Barrier /= (1.0 + std::sqrt(U/static_cast<G4double>(2*ARes)));
77 G4double Result = 1.12*(AresOneThird + AejectOneThird) -
78 0.86*(AresOneThird+AejectOneThird)/(AresOneThird*AejectOneThird);