G4NeutronField Class Reference

#include <G4NeutronField.hh>

Inheritance 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 ()

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::G4NeutronField

(

G4V3DNucleus *

nucleus

)

Definition at line 45 of file G4NeutronField.cc.

References density, python.hepunit::fermi, G4V3DNucleus::GetCharge(), G4V3DNucleus::GetMassNumber(), G4V3DNucleus::GetOuterRadius(), G4FermiMomentum::Init(), python.hepunit::m, and G4VNuclearField::theNucleus.

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

G4NeutronField::~G4NeutronField ( )

virtual

Definition at line 77 of file G4NeutronField.cc.

{ }

Member Function Documentation

G4double G4NeutronField::GetBarrier ( )

virtual

Implements G4VNuclearField.

Definition at line 93 of file G4NeutronField.cc.

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

G4double G4NeutronField::GetField ( const G4ThreeVector &  aPosition )

virtual

Implements G4VNuclearField.

Definition at line 80 of file G4NeutronField.cc.

References python.hepunit::fermi, CLHEP::Hep3Vector::mag(), python.hepunit::neutron_mass_c2, and test::x.

{
  G4double x = aPosition.mag();
  G4int index = static_cast<G4int>(x/(0.3*fermi) );
  if(index+2> static_cast<G4int>(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);
}

---

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

• G4NeutronField.hh
• G4NeutronField.cc