G4ProtonField Class Reference

#include <G4ProtonField.hh>

Inheritance diagram for G4ProtonField:

Collaboration diagram for G4ProtonField:

Public Member Functions
	G4ProtonField (G4V3DNucleus *nucleus)
virtual	~G4ProtonField ()
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
	G4ProtonField (const G4ProtonField &right)
const G4ProtonField &	operator= (const G4ProtonField &right)
int	operator== (const G4ProtonField &right) const
int	operator!= (const G4ProtonField &right) const
G4double	GetDensity (const G4ThreeVector &aPosition)
G4double	GetFermiMomentum (const G4double aDensity)

Private Attributes
G4int	theA
G4int	theZ
G4double	theBarrier
G4double	theRadius
G4FermiMomentum	theFermi
const G4VNuclearDensity *	theDensity
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 G4ProtonField.hh.

Constructor & Destructor Documentation

G4ProtonField() [1/2]

G4ProtonField::G4ProtonField

(

G4V3DNucleus *

nucleus

)

Definition at line 47 of file G4ProtonField.cc.

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

~G4ProtonField()

G4ProtonField::~G4ProtonField ( )

virtual

Definition at line 79 of file G4ProtonField.cc.

{ }

G4ProtonField() [2/2]

G4ProtonField::G4ProtonField ( const G4ProtonField &  right )

private

Member Function Documentation

GetBarrier()

G4double G4ProtonField::GetBarrier ( )

virtual

Implements G4VNuclearField.

Definition at line 98 of file G4ProtonField.cc.

{
  G4double coulombBarrier = (1.44/1.14) * MeV * theZ / (1.0 + G4Pow::GetInstance()->Z13(theA));
//GF   G4double bindingEnergy = G4NucleiPropertiesTable::GetBindingEnergy(Z, A);
  G4double bindingEnergy =0;
/*
 *   G4cout << " coulombBarrier/bindingEnergy : " 
 *       << coulombBarrier << " /" << bindingEnergy << G4endl;
 */
  return bindingEnergy/theA+coulombBarrier;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetDensity()

G4double G4ProtonField::GetDensity ( const G4ThreeVector &  aPosition )

inlineprivate

Definition at line 68 of file G4ProtonField.hh.

  {
    return theDensity->GetDensity(aPosition);
  }

Here is the call graph for this function:

Here is the caller graph for this function:

GetFermiMomentum()

G4double G4ProtonField::GetFermiMomentum ( const G4double  aDensity )

inlineprivate

Definition at line 72 of file G4ProtonField.hh.

  {
    return theFermi.GetFermiMomentum(aDensity);
  }

Here is the call graph for this function:

Here is the caller graph for this function:

GetField()

G4double G4ProtonField::GetField ( const G4ThreeVector &  aPosition )

virtual

Implements G4VNuclearField.

Definition at line 82 of file G4ProtonField.cc.

{
//G4cout << " Fermi Potential " << (fermiMom*fermiMom)/(2*proton_mass_c2) <<G4endl;
  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);
  G4double y = -1*(fermiMom*fermiMom)/(2*proton_mass_c2)+theBarrier;
//  G4cout <<" Protonfield test "<<index<<" "<< x1<<" "<<y1<<" "<<x2<<" "<<y2<<" "<<x<<" "<<y<<" "<<theBarrier<<G4endl;
  return y;
}

Here is the call graph for this function:

operator!=()

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

private

operator=()

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

private

operator==()

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

private

Member Data Documentation

theA

G4int G4ProtonField::theA

private

Definition at line 77 of file G4ProtonField.hh.

theBarrier

G4double G4ProtonField::theBarrier

private

Definition at line 79 of file G4ProtonField.hh.

theDensity

const G4VNuclearDensity* G4ProtonField::theDensity

private

Definition at line 82 of file G4ProtonField.hh.

theFermi

G4FermiMomentum G4ProtonField::theFermi

private

Definition at line 81 of file G4ProtonField.hh.

theFermiMomBuffer

std::vector<G4double> G4ProtonField::theFermiMomBuffer

private

Definition at line 84 of file G4ProtonField.hh.

theRadius

G4double G4ProtonField::theRadius

private

Definition at line 80 of file G4ProtonField.hh.

theZ

G4int G4ProtonField::theZ

private

Definition at line 78 of file G4ProtonField.hh.

---

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

• G4ProtonField.hh
• G4ProtonField.cc