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

(

G4V3DNucleus *

nucleus

)

Definition at line 46 of file G4ProtonField.cc.

References density, python.hepunit::fermi, GetBarrier(), 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();
  theBarrier = GetBarrier();
  theRadius = 2.*theNucleus->GetOuterRadius();
  theFermi.Init(theA, theZ);
  G4double aR=0;
  while(aR<theRadius)
  {
    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,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 ( )

virtual

Definition at line 80 of file G4ProtonField.cc.

{ }

Member Function Documentation

G4double G4ProtonField::GetBarrier ( )

virtual

Implements G4VNuclearField.

Definition at line 99 of file G4ProtonField.cc.

References G4InuclSpecialFunctions::bindingEnergy(), and python.hepunit::MeV.

Referenced by G4ProtonField().

{
  G4double coulombBarrier = (1.44/1.14) * MeV * theZ / (1.0 + std::pow(theA,1./3.));
//GF   G4double bindingEnergy = G4NucleiPropertiesTable::GetBindingEnergy(Z, A);
  G4double bindingEnergy =0;
/*
 *   G4cout << " coulombBarrier/bindingEnergy : " 
 *       << coulombBarrier << " /" << bindingEnergy << G4endl;
 */
  return bindingEnergy/theA+coulombBarrier;
}

G4double G4ProtonField::GetField ( const G4ThreeVector &  aPosition )

virtual

Implements G4VNuclearField.

Definition at line 83 of file G4ProtonField.cc.

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

{
//G4cout << " Fermi Potential " << (fermiMom*fermiMom)/(2*proton_mass_c2) <<G4endl;
  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);
  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