G4ProtonField Class Reference

#include <G4ProtonField.hh>

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

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

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);  
  }
}

~G4ProtonField()

G4ProtonField::~G4ProtonField ( )

virtual

Definition at line 79 of file G4ProtonField.cc.

{ }

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;
}

Referenced by G4ProtonField().

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;
}

---

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

• G4ProtonField.hh
• G4ProtonField.cc