Geant4 11.2.2
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4PionZeroField Class Reference

#include <G4PionZeroField.hh>

+ Inheritance diagram for G4PionZeroField:

Public Member Functions

 G4PionZeroField (G4V3DNucleus *nucleus, G4double coeff=0.042 *CLHEP::fermi)
 
virtual ~G4PionZeroField ()
 
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
G4V3DNucleustheNucleus
 
const G4double radius
 

Detailed Description

Definition at line 47 of file G4PionZeroField.hh.

Constructor & Destructor Documentation

◆ G4PionZeroField()

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

Definition at line 47 of file G4PionZeroField.cc.

48 : G4VNuclearField(nucleus)
49{
50 theCoeff = coeff;
51}
G4VNuclearField(G4V3DNucleus *aNucleus=0)

◆ ~G4PionZeroField()

G4PionZeroField::~G4PionZeroField ( )
virtual

Definition at line 54 of file G4PionZeroField.cc.

55{ }

Member Function Documentation

◆ GetBarrier()

G4double G4PionZeroField::GetBarrier ( )
virtual

Implements G4VNuclearField.

Definition at line 77 of file G4PionZeroField.cc.

78{
79 return 0;
80}

◆ GetCoeff()

virtual G4double G4PionZeroField::GetCoeff ( )
inlinevirtual

Reimplemented from G4VNuclearField.

Definition at line 62 of file G4PionZeroField.hh.

62{ return theCoeff; }

◆ GetField()

G4double G4PionZeroField::GetField ( const G4ThreeVector & aPosition)
virtual

Implements G4VNuclearField.

Definition at line 57 of file G4PionZeroField.cc.

58{
59// Field is 0 out of the nucleus!
60 if(aPosition.mag() >= radius) return 0.0;
61
62 G4double pionZeroMass = G4PionZero::PionZero()->GetPDGMass();
65
67 G4double nucleusMass = Z*proton_mass_c2+(A-Z)*neutron_mass_c2+bindingEnergy;
68 G4double reducedMass = pionZeroMass*nucleusMass/(pionZeroMass+nucleusMass);
69
70
71 G4double density = A*theNucleus->GetNuclearDensity()->GetDensity(aPosition);
72 G4double nucleonMass = (proton_mass_c2+neutron_mass_c2)/2;
73
74 return 2.*pi*hbarc*hbarc/reducedMass*(1+pionZeroMass/nucleonMass)*theCoeff*density;
75}
double G4double
Definition G4Types.hh:83
int G4int
Definition G4Types.hh:85
const G4double A[17]
double mag() const
static G4double GetBindingEnergy(const G4int A, const G4int Z)
static G4PionZero * PionZero()
virtual const G4VNuclearDensity * GetNuclearDensity() const =0
virtual G4int GetCharge()=0
virtual G4int GetMassNumber()=0
G4double GetDensity(const G4ThreeVector &aPosition) const
G4V3DNucleus * theNucleus
const G4double radius
G4double bindingEnergy(G4int A, G4int Z)

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