G4SigmaZeroField Class Reference

#include <G4SigmaZeroField.hh>

Inheritance diagram for G4SigmaZeroField:

Public Member Functions
	G4SigmaZeroField (G4V3DNucleus *nucleus, G4double coeff=0.36 *CLHEP::fermi)
virtual	~G4SigmaZeroField ()
virtual G4double	GetField (const G4ThreeVector &aPosition)
virtual G4double	GetBarrier ()
virtual G4double	GetCoeff ()
Public Member Functions inherited from G4VNuclearField
	G4VNuclearField (G4V3DNucleus *aNucleus=0)
virtual	~G4VNuclearField ()
void	SetNucleus (G4V3DNucleus *aNucleus)

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

Detailed Description

Definition at line 46 of file G4SigmaZeroField.hh.

Constructor & Destructor Documentation

G4SigmaZeroField::G4SigmaZeroField	(	G4V3DNucleus *	nucleus,
		G4double	coeff = 0.36*CLHEP::fermi
	)

Definition at line 47 of file G4SigmaZeroField.cc.

  : G4VNuclearField(nucleus)
{
  theCoeff = coeff;
}

G4SigmaZeroField::~G4SigmaZeroField ( )

virtual

Definition at line 53 of file G4SigmaZeroField.cc.

{ }

Member Function Documentation

G4double G4SigmaZeroField::GetBarrier ( )

virtual

Implements G4VNuclearField.

Definition at line 76 of file G4SigmaZeroField.cc.

{
  return 0.;
}

virtual G4double G4SigmaZeroField::GetCoeff ( )

inlinevirtual

Reimplemented from G4VNuclearField.

Definition at line 61 of file G4SigmaZeroField.hh.

{ return theCoeff; }

G4double G4SigmaZeroField::GetField ( const G4ThreeVector &  aPosition )

virtual

Implements G4VNuclearField.

Definition at line 56 of file G4SigmaZeroField.cc.

References G4InuclSpecialFunctions::bindingEnergy(), density, G4NucleiProperties::GetBindingEnergy(), G4V3DNucleus::GetCharge(), G4VNuclearDensity::GetDensity(), G4V3DNucleus::GetMassNumber(), G4V3DNucleus::GetNuclearDensity(), G4ParticleDefinition::GetPDGMass(), python.hepunit::hbarc, CLHEP::Hep3Vector::mag(), python.hepunit::neutron_mass_c2, python.hepunit::pi, python.hepunit::proton_mass_c2, G4VNuclearField::radius, G4SigmaZero::SigmaZero(), and G4VNuclearField::theNucleus.

{
// Field is 0 out of the nucleus!
  if(aPosition.mag() >= radius) return 0.0;

  G4double sigmaZeroMass = G4SigmaZero::SigmaZero()->GetPDGMass();

  G4int A = theNucleus->GetMassNumber();
  G4int Z = theNucleus->GetCharge();
  G4double bindingEnergy = G4NucleiProperties::GetBindingEnergy(A, Z);
  G4double nucleusMass = Z*proton_mass_c2+(A-Z)*neutron_mass_c2+bindingEnergy;
  G4double reducedMass = sigmaZeroMass*nucleusMass/(sigmaZeroMass+nucleusMass);

  const G4VNuclearDensity * nuclearDensity=theNucleus->GetNuclearDensity();
  G4double density = nuclearDensity->GetDensity(aPosition);

  return -2.*pi*hbarc*hbarc/reducedMass*(2.0)*theCoeff*density;
}

---

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

• G4SigmaZeroField.hh
• G4SigmaZeroField.cc