#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 () |
Definition at line 46 of file G4SigmaZeroField.hh.
G4SigmaZeroField::G4SigmaZeroField | ( | G4V3DNucleus * | nucleus, | |
G4double | coeff = 0.36 *CLHEP::fermi | |||
) |
Definition at line 47 of file G4SigmaZeroField.cc.
00048 : G4VNuclearField(nucleus) 00049 { 00050 theCoeff = coeff; 00051 }
G4SigmaZeroField::~G4SigmaZeroField | ( | ) | [virtual] |
G4double G4SigmaZeroField::GetBarrier | ( | ) | [virtual] |
virtual G4double G4SigmaZeroField::GetCoeff | ( | ) | [inline, virtual] |
G4double G4SigmaZeroField::GetField | ( | const G4ThreeVector & | aPosition | ) | [virtual] |
Implements G4VNuclearField.
Definition at line 56 of file G4SigmaZeroField.cc.
References G4InuclSpecialFunctions::bindingEnergy(), G4NucleiProperties::GetBindingEnergy(), G4V3DNucleus::GetCharge(), G4VNuclearDensity::GetDensity(), G4V3DNucleus::GetMassNumber(), G4V3DNucleus::GetNuclearDensity(), G4ParticleDefinition::GetPDGMass(), G4INCL::Math::pi, G4VNuclearField::radius, G4SigmaZero::SigmaZero(), and G4VNuclearField::theNucleus.
00057 { 00058 // Field is 0 out of the nucleus! 00059 if(aPosition.mag() >= radius) return 0.0; 00060 00061 G4double sigmaZeroMass = G4SigmaZero::SigmaZero()->GetPDGMass(); 00062 00063 G4int A = theNucleus->GetMassNumber(); 00064 G4int Z = theNucleus->GetCharge(); 00065 G4double bindingEnergy = G4NucleiProperties::GetBindingEnergy(A, Z); 00066 G4double nucleusMass = Z*proton_mass_c2+(A-Z)*neutron_mass_c2+bindingEnergy; 00067 G4double reducedMass = sigmaZeroMass*nucleusMass/(sigmaZeroMass+nucleusMass); 00068 00069 const G4VNuclearDensity * nuclearDensity=theNucleus->GetNuclearDensity(); 00070 G4double density = nuclearDensity->GetDensity(aPosition); 00071 00072 return -2.*pi*hbarc*hbarc/reducedMass*(2.0)*theCoeff*density; 00073 }