61 if (ZRes > ARes || ARes < 1) {
62 std::ostringstream errOs;
63 errOs <<
"G4CoulombBarrier::GetCoulombBarrier: ";
64 errOs <<
"Wrong values for ";
65 errOs <<
"residual nucleus A = " << ARes <<
" ";
66 errOs <<
"and residual nucleus Z = " << ZRes <<
G4endl;
83 if(
GetA()==1 &&
GetZ()==1){ rho=0.0;}
static G4Pow * GetInstance()
virtual ~G4CoulombBarrier()
G4double GetCoulombBarrier(G4int ARes, G4int ZRes, G4double U) const
G4double Z13(G4int Z) const
virtual G4double BarrierPenetrationFactor(G4double) const
static const double fermi