55 return (
A < 65) ? blimit :
56 BarashenkovFissionBarrier(
A,
Z)/(1.0 + std::sqrt(U/(
G4double)(2*
A)));
72 G4double x = (aCoul/(2.0*aSurf))*(Z*
Z)/static_cast<G4double>(
A);
73 x /= (1.0 - k*(N-
Z)*(N-Z)/
static_cast<G4double>(
A*
A));
77 if (x <= 2./3.) { BF0 *= 0.38*(3./4.-
x); }
78 else { BF0 *= 0.83*(1. -
x)*(1. - x)*(1. -
x); }
80 G4int d = N - 2*(N/2) + Z - 2*(Z/2);
83 SPtr->GetPairingCorrection(N,Z,res);
86 return BF0 + D*d - res;
G4double BarashenkovFissionBarrier(G4int A, G4int Z) const
virtual G4double FissionBarrier(G4int A, G4int Z, G4double U) const final
static G4NuclearLevelData * GetInstance()
static G4Pow * GetInstance()
virtual ~G4FissionBarrier()
static constexpr double MeV
double A(double temperature)
G4ShellCorrection * GetShellCorrection()
const G4CameronShellPlusPairingCorrections * GetCameronShellPlusPairingCorrections() const
static constexpr double GeV
const G4CameronShellPlusPairingCorrections * SPtr
G4double Z23(G4int Z) const