G4PionMinusField Class Reference

#include <G4PionMinusField.hh>

Inheritance diagram for G4PionMinusField:

Collaboration diagram for G4PionMinusField:

Public Member Functions
	G4PionMinusField (G4V3DNucleus *nucleus, G4double coeff=0.042 *CLHEP::fermi)
virtual	~G4PionMinusField ()
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 47 of file G4PionMinusField.hh.

Constructor & Destructor Documentation

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

Definition at line 49 of file G4PionMinusField.cc.

  : G4VNuclearField(nucleus)
{
  theCoeff = coeff;
}

G4PionMinusField::~G4PionMinusField ( )

virtual

Definition at line 56 of file G4PionMinusField.cc.

{ }

Member Function Documentation

G4double G4PionMinusField::GetBarrier ( )

virtual

Implements G4VNuclearField.

Definition at line 78 of file G4PionMinusField.cc.

{
  G4int A = theNucleus->GetMassNumber();
  G4int Z = theNucleus->GetCharge();
  G4double coulombBarrier = (1.44/1.14) * MeV * Z / (1.0 + G4Pow::GetInstance()->Z13(A));
  return -coulombBarrier;
}

Here is the call graph for this function:

Here is the caller graph for this function:

virtual G4double G4PionMinusField::GetCoeff ( )

inlinevirtual

Reimplemented from G4VNuclearField.

Definition at line 64 of file G4PionMinusField.hh.

{ return theCoeff; }

G4double G4PionMinusField::GetField ( const G4ThreeVector &  aPosition )

virtual

Implements G4VNuclearField.

Definition at line 59 of file G4PionMinusField.cc.

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

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

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

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

Here is the call graph for this function:

---

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

• source/geant4.10.03.p03/source/processes/hadronic/models/binary_cascade/include/G4PionMinusField.hh
• source/geant4.10.03.p03/source/processes/hadronic/models/binary_cascade/src/G4PionMinusField.cc