G4PionZeroField Class Reference

#include <G4PionZeroField.hh>

Inheritance diagram for G4PionZeroField:

Public Member Functions
	G4PionZeroField (G4V3DNucleus *nucleus, G4double coeff=0.042 *CLHEP::fermi)
virtual	~G4PionZeroField ()
virtual G4double	GetField (const G4ThreeVector &aPosition)
virtual G4double	GetBarrier ()
virtual G4double	GetCoeff ()

---

Detailed Description

Definition at line 47 of file G4PionZeroField.hh.

---

Constructor & Destructor Documentation

G4PionZeroField::G4PionZeroField	(	G4V3DNucleus *	nucleus,
		G4double	coeff = 0.042 *CLHEP::fermi
	)

Definition at line 47 of file G4PionZeroField.cc.

  : G4VNuclearField(nucleus)
{
  theCoeff = coeff;
}

G4PionZeroField::~G4PionZeroField

(

)

[virtual]

Definition at line 54 of file G4PionZeroField.cc.

{ }

---

Member Function Documentation

G4double G4PionZeroField::GetBarrier

(

)

[virtual]

Implements G4VNuclearField.

Definition at line 77 of file G4PionZeroField.cc.

{
  return 0;
}

virtual G4double G4PionZeroField::GetCoeff

(

)

[inline, virtual]

Reimplemented from G4VNuclearField.

Definition at line 62 of file G4PionZeroField.hh.

{ return theCoeff; }

G4double G4PionZeroField::GetField

(

const G4ThreeVector &

aPosition

)

[virtual]

Implements G4VNuclearField.

Definition at line 57 of file G4PionZeroField.cc.

References G4InuclSpecialFunctions::bindingEnergy(), G4NucleiProperties::GetBindingEnergy(), G4V3DNucleus::GetCharge(), G4VNuclearDensity::GetDensity(), G4V3DNucleus::GetMassNumber(), G4V3DNucleus::GetNuclearDensity(), G4ParticleDefinition::GetPDGMass(), G4INCL::Math::pi, G4PionZero::PionZero(), G4VNuclearField::radius, and G4VNuclearField::theNucleus.

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

  G4double pionZeroMass = G4PionZero::PionZero()->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 = pionZeroMass*nucleusMass/(pionZeroMass+nucleusMass);


  G4double density = A*theNucleus->GetNuclearDensity()->GetDensity(aPosition);
  G4double nucleonMass = (proton_mass_c2+neutron_mass_c2)/2;

  return 2.*pi*hbarc*hbarc/reducedMass*(1+pionZeroMass/nucleonMass)*theCoeff*density;
}

---

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

• G4PionZeroField.hh
• G4PionZeroField.cc

---