G4PionMinusField Class Reference

#include <G4PionMinusField.hh>

Inheritance diagram for G4PionMinusField:

G4VNuclearField

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

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 48 of file G4PionMinusField.cc.

00049   : G4VNuclearField(nucleus)
00050 {
00051   theCoeff = coeff;
00052 }

G4PionMinusField::~G4PionMinusField (  )  [virtual]

Definition at line 55 of file G4PionMinusField.cc.

00056 { }


Member Function Documentation

G4double G4PionMinusField::GetBarrier (  )  [virtual]

Implements G4VNuclearField.

Definition at line 77 of file G4PionMinusField.cc.

References G4V3DNucleus::GetCharge(), G4V3DNucleus::GetMassNumber(), and G4VNuclearField::theNucleus.

Referenced by GetField().

00078 {
00079   G4int A = theNucleus->GetMassNumber();
00080   G4int Z = theNucleus->GetCharge();
00081   G4double coulombBarrier = (1.44/1.14) * MeV * Z / (1.0 + std::pow(A,1./3.));
00082   return -coulombBarrier;
00083 }

virtual G4double G4PionMinusField::GetCoeff (  )  [inline, virtual]

Reimplemented from G4VNuclearField.

Definition at line 64 of file G4PionMinusField.hh.

00064 { return theCoeff; }

G4double G4PionMinusField::GetField ( const G4ThreeVector aPosition  )  [virtual]

Implements G4VNuclearField.

Definition at line 58 of file G4PionMinusField.cc.

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

00059 {
00060 // Field is 0 out of the nucleus!
00061   if(aPosition.mag() >= radius) return 0.0;
00062 
00063   G4int A = theNucleus->GetMassNumber();
00064   G4int Z = theNucleus->GetCharge();
00065   G4double pionMinusMass = G4PionMinus::PionMinus()->GetPDGMass();
00066   G4double bindingEnergy = G4NucleiProperties::GetBindingEnergy(G4lrint(A), G4lrint(Z));
00067   G4double nucleusMass = Z*proton_mass_c2+(A-Z)*neutron_mass_c2+bindingEnergy;
00068   G4double reducedMass = pionMinusMass*nucleusMass/(pionMinusMass+nucleusMass);
00069 
00070   G4double density = A*theNucleus->GetNuclearDensity()->GetDensity(aPosition);
00071   G4double nucleonMass = (proton_mass_c2+neutron_mass_c2)/2;
00072 
00073   return 2.*pi*hbarc*hbarc/reducedMass*(1+pionMinusMass/nucleonMass)*theCoeff*density + GetBarrier();
00074 }


The documentation for this class was generated from the following files:
Generated on Mon May 27 17:52:55 2013 for Geant4 by  doxygen 1.4.7