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