G4SigmaZeroField Class Reference

#include <G4SigmaZeroField.hh>

Inheritance diagram for G4SigmaZeroField:

G4VNuclearField

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

Detailed Description

Definition at line 46 of file G4SigmaZeroField.hh.


Constructor & Destructor Documentation

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]

Definition at line 53 of file G4SigmaZeroField.cc.

00054 { }


Member Function Documentation

G4double G4SigmaZeroField::GetBarrier (  )  [virtual]

Implements G4VNuclearField.

Definition at line 76 of file G4SigmaZeroField.cc.

00077 {
00078   return 0.;
00079 }

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

Reimplemented from G4VNuclearField.

Definition at line 61 of file G4SigmaZeroField.hh.

00061 { return theCoeff; }

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 }


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