#include <G4PionMinusField.hh>
Inheritance 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 () |
Definition at line 47 of file G4PionMinusField.hh.
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] |
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] |
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 }