28 #ifndef G4NuclearFermiDensity_h
29 #define G4NuclearFermiDensity_h 1
35 #include <CLHEP/Units/PhysicalConstants.h>
48 return 1./(1.+
G4Exp((aPosition.mag()-
theR)/
a));
53 return (maxRelativeDenisty>0 && maxRelativeDenisty <= 1 ) ?
60 if (currentR > 40*
theR ) {
return 0;}
CLHEP::Hep3Vector G4ThreeVector
G4double GetDeriv(const G4ThreeVector &aPosition) const
G4double GetRelativeDensity(const G4ThreeVector &aPosition) const
G4NuclearFermiDensity(G4int anA, G4int aZ)
G4double G4Log(G4double x)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
G4double GetDensity(const G4ThreeVector &aPosition) const
G4double GetRadius(const G4double maxRelativeDenisty) const