86 unsigned int index =
static_cast<unsigned int>(x/(0.3*
fermi));
107 return bindingEnergy/
theA+coulombBarrier;
void Init(G4int anA, G4int aZ)
G4double Z13(G4int Z) const
static constexpr double MeV
Float_t y1[n_points_granero]
Float_t x1[n_points_granero]
std::vector< G4double > theFermiMomBuffer
Float_t y2[n_points_geant4]
static constexpr double proton_mass_c2
static constexpr double m
static constexpr double fermi
static G4Pow * GetInstance()
G4ProtonField(G4V3DNucleus *nucleus)
virtual G4double GetBarrier()
virtual G4int GetCharge()=0
virtual G4int GetMassNumber()=0
G4double GetFermiMomentum(const G4double aDensity)
G4V3DNucleus * theNucleus
virtual G4double GetField(const G4ThreeVector &aPosition)
Float_t x2[n_points_geant4]
G4double GetDensity(const G4ThreeVector &aPosition)
G4double bindingEnergy(G4int A, G4int Z)
virtual G4double GetOuterRadius()=0