G4NuclearFermiDensity Class Reference

#include <G4NuclearFermiDensity.hh>

Inheritance diagram for G4NuclearFermiDensity:

Public Member Functions
	G4NuclearFermiDensity (G4int anA, G4int aZ)
	~G4NuclearFermiDensity ()
G4double	GetRelativeDensity (const G4ThreeVector &aPosition) const
G4double	GetRadius (const G4double maxRelativeDenisty) const
G4double	GetDeriv (const G4ThreeVector &aPosition) const
Public Member Functions inherited from G4VNuclearDensity
	G4VNuclearDensity ()
virtual	~G4VNuclearDensity ()
G4double	GetDensity (const G4ThreeVector &aPosition) const

Additional Inherited Members
Protected Member Functions inherited from G4VNuclearDensity
void	Setrho0 (G4double arho0)
G4double	Getrho0 () const

Detailed Description

Definition at line 40 of file G4NuclearFermiDensity.hh.

Constructor & Destructor Documentation

G4NuclearFermiDensity()

G4NuclearFermiDensity::G4NuclearFermiDensity	(	G4int	anA,
		G4int	aZ )

Definition at line 33 of file G4NuclearFermiDensity.cc.

 : theA(anA), a(0.545 * fermi) 
{
  G4double a13 = G4Pow::GetInstance()->Z13(anA);
  const G4double r0 = 1.16 * (1. - 1.16/(a13*a13)) * fermi;
  theR = r0 * a13;
  Setrho0(3./ (4.*pi *r0*r0*r0 * theA * (1. + sqr(a/theR)*pi2 )));
}

~G4NuclearFermiDensity()

G4NuclearFermiDensity::~G4NuclearFermiDensity

(

)

Definition at line 42 of file G4NuclearFermiDensity.cc.

{}

Member Function Documentation

GetDeriv()

G4double G4NuclearFermiDensity::GetDeriv ( const G4ThreeVector & aPosition ) const

inlinevirtual

Implements G4VNuclearDensity.

Definition at line 57 of file G4NuclearFermiDensity.hh.

    {
      G4double currentR=aPosition.mag();
      if (currentR > 40*theR  ) {return 0;}
      else return -G4Exp((currentR-theR)/a) * sqr(GetDensity(aPosition)) / (a*Getrho0());
    }   

GetRadius()

G4double G4NuclearFermiDensity::GetRadius ( const G4double maxRelativeDenisty ) const

inlinevirtual

Implements G4VNuclearDensity.

Definition at line 51 of file G4NuclearFermiDensity.hh.

    {
      return (maxRelativeDenisty>0 && maxRelativeDenisty <= 1 ) ?
             (theR + a*G4Log((1-maxRelativeDenisty+G4Exp(-1*theR/a))/maxRelativeDenisty))  : DBL_MAX;
    }

GetRelativeDensity()

G4double G4NuclearFermiDensity::GetRelativeDensity ( const G4ThreeVector & aPosition ) const

inlinevirtual

Implements G4VNuclearDensity.

Definition at line 46 of file G4NuclearFermiDensity.hh.

    {
      return 1./(1.+G4Exp((aPosition.mag()-theR)/a));
    }

---

The documentation for this class was generated from the following files:

• G4NuclearFermiDensity.hh
• G4NuclearFermiDensity.cc