61 res = (aZ >= 70) ? 0.80 :
62 (((0.2357e-5*aZ) - 0.42679e-3)*aZ + 0.27035e-1)*aZ + 0.19025;
63 res += 0.06*(
GetA() - 1);
65 }
else if(
GetZ() == 2 &&
GetA() <= 4) {
66 res = (aZ >= 70) ? 0.98 :
67 (((0.23684e-5*aZ) - 0.42143e-3)*aZ + 0.25222e-1)*aZ + 0.46699;
68 res += 0.12*(4 -
GetA());
virtual G4double BarrierPenetrationFactor(G4int Eexc) const
virtual ~G4CoulombBarrier()
G4CoulombBarrier(G4int anA, G4int aZ)
G4double GetCoulombBarrier(G4int ARes, G4int ZRes, G4double U) const
static G4double RadiusCB(G4int Z, G4int A)
G4double GetRho(void) const
void SetParameters(G4double rho, G4double r0)
static constexpr double elm_coupling
static constexpr double fermi
std::map< G4String, G4AttDef > * GetInstance(const G4String &storeKey, G4bool &isNew)