G4NeutronField Class Reference

#include <G4NeutronField.hh>

Inheritance diagram for G4NeutronField:

Collaboration diagram for G4NeutronField:

Public Member Functions
	G4NeutronField (G4V3DNucleus *nucleus)
virtual	~G4NeutronField ()
virtual G4double	GetField (const G4ThreeVector &aPosition)
virtual G4double	GetBarrier ()
Public Member Functions inherited from G4VNuclearField
	G4VNuclearField (G4V3DNucleus *aNucleus=0)
virtual	~G4VNuclearField ()
void	SetNucleus (G4V3DNucleus *aNucleus)
virtual G4double	GetCoeff ()

Private Member Functions
	G4NeutronField (const G4NeutronField &right)
const G4NeutronField &	operator= (const G4NeutronField &right)
int	operator== (const G4NeutronField &right) const
int	operator!= (const G4NeutronField &right) const
G4double	GetDensity (const G4ThreeVector &aPosition)
G4double	GetFermiMomentum (const G4double aDensity)

Private Attributes
G4FermiMomentum	theFermi
G4int	theA
G4int	theZ
const G4VNuclearDensity *	theDensity
G4double	theR
std::vector< G4double >	theFermiMomBuffer

Additional Inherited Members
Protected Attributes inherited from G4VNuclearField
G4V3DNucleus *	theNucleus
const G4double	radius

Detailed Description

Definition at line 48 of file G4NeutronField.hh.

Constructor & Destructor Documentation

G4NeutronField() [1/2]

G4NeutronField::G4NeutronField

(

G4V3DNucleus *

nucleus

)

Definition at line 45 of file G4NeutronField.cc.

                                                      : 
    G4VNuclearField(aNucleus), theDensity(theNucleus->GetNuclearDensity())
{ 
  theA = theNucleus->GetMassNumber();
  theZ = theNucleus->GetCharge();
  theFermi.Init(theA, theZ);
  theR = 2.*theNucleus->GetOuterRadius();
  for (G4double aR=0.;aR<theR; aR+=0.3*fermi)
  {
    G4ThreeVector aPosition(0,0,aR);
    G4double density = GetDensity(aPosition);
    G4double fermiMom = GetFermiMomentum(density);
    theFermiMomBuffer.push_back(fermiMom);
  }
  {
  G4ThreeVector aPosition(0,0,theR);
  G4double density = GetDensity(aPosition);
  G4double fermiMom = GetFermiMomentum(density);
  theFermiMomBuffer.push_back(fermiMom);  
  }
  {
  G4ThreeVector aPosition(0,0,theR+0.001*fermi);
  theFermiMomBuffer.push_back(0);  
  }
  {
  G4ThreeVector aPosition(0,0,1.*m);
  theFermiMomBuffer.push_back(0);  
  }
}

Here is the call graph for this function:

~G4NeutronField()

G4NeutronField::~G4NeutronField ( )

virtual

Definition at line 75 of file G4NeutronField.cc.

{ }

G4NeutronField() [2/2]

G4NeutronField::G4NeutronField ( const G4NeutronField &  right )

private

Member Function Documentation

GetBarrier()

G4double G4NeutronField::GetBarrier ( )

virtual

Implements G4VNuclearField.

Definition at line 91 of file G4NeutronField.cc.

{
/*
 *   G4double A = theNucleus->GetMassNumber();
 *   G4double Z = theNucleus->GetCharge();
 * 
 *   return G4NucleiPropertiesTable::GetBindingEnergy(Z, A)/A;
 */
    return 0.;
}

GetDensity()

G4double G4NeutronField::GetDensity ( const G4ThreeVector &  aPosition )

inlineprivate

Definition at line 68 of file G4NeutronField.hh.

  {
    return theDensity->GetDensity(aPosition);
  }

Here is the call graph for this function:

Here is the caller graph for this function:

GetFermiMomentum()

G4double G4NeutronField::GetFermiMomentum ( const G4double  aDensity )

inlineprivate

Definition at line 73 of file G4NeutronField.hh.

  { 
    return theFermi.GetFermiMomentum(aDensity);
  }

Here is the call graph for this function:

Here is the caller graph for this function:

GetField()

G4double G4NeutronField::GetField ( const G4ThreeVector &  aPosition )

virtual

Implements G4VNuclearField.

Definition at line 78 of file G4NeutronField.cc.

{
  G4double x = aPosition.mag();
  unsigned int index = static_cast<unsigned int>(x/(0.3*fermi));
  if( (index+2) > theFermiMomBuffer.size()) return theFermiMomBuffer.back();
  G4double y1 = theFermiMomBuffer[index];
  G4double y2 = theFermiMomBuffer[index+1];
  G4double x1 = (0.3*fermi)*index;
  G4double x2 = (0.3*fermi)*(index+1);
  G4double fermiMom = y1 + (x-x1)*(y2-y1)/(x2-x1);
  return -1*(fermiMom*fermiMom)/(2*neutron_mass_c2);
}

Here is the call graph for this function:

operator!=()

int G4NeutronField::operator!= ( const G4NeutronField &  right ) const

private

operator=()

const G4NeutronField& G4NeutronField::operator= ( const G4NeutronField &  right )

private

operator==()

int G4NeutronField::operator== ( const G4NeutronField &  right ) const

private

Member Data Documentation

theA

G4int G4NeutronField::theA

private

Definition at line 79 of file G4NeutronField.hh.

theDensity

const G4VNuclearDensity* G4NeutronField::theDensity

private

Definition at line 81 of file G4NeutronField.hh.

theFermi

G4FermiMomentum G4NeutronField::theFermi

private

Definition at line 78 of file G4NeutronField.hh.

theFermiMomBuffer

std::vector<G4double> G4NeutronField::theFermiMomBuffer

private

Definition at line 84 of file G4NeutronField.hh.

theR

G4double G4NeutronField::theR

private

Definition at line 82 of file G4NeutronField.hh.

theZ

G4int G4NeutronField::theZ

private

Definition at line 80 of file G4NeutronField.hh.

---

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

• G4NeutronField.hh
• G4NeutronField.cc