Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4PionMinusField Class Reference

#include <G4PionMinusField.hh>

Inheritance diagram for G4PionMinusField:
Collaboration 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 ()
 
- 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 G4PionMinusField.hh.

Constructor & Destructor Documentation

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

Definition at line 49 of file G4PionMinusField.cc.

50  : G4VNuclearField(nucleus)
51 {
52  theCoeff = coeff;
53 }
G4VNuclearField(G4V3DNucleus *aNucleus=0)
G4PionMinusField::~G4PionMinusField ( )
virtual

Definition at line 56 of file G4PionMinusField.cc.

57 { }

Member Function Documentation

G4double G4PionMinusField::GetBarrier ( )
virtual

Implements G4VNuclearField.

Definition at line 78 of file G4PionMinusField.cc.

79 {
82  G4double coulombBarrier = (1.44/1.14) * MeV * Z / (1.0 + G4Pow::GetInstance()->Z13(A));
83  return -coulombBarrier;
84 }
static G4Pow * GetInstance()
Definition: G4Pow.cc:55
virtual G4int GetCharge()=0
virtual G4int GetMassNumber()=0
int G4int
Definition: G4Types.hh:78
G4double Z13(G4int Z) const
Definition: G4Pow.hh:127
double A(double temperature)
G4V3DNucleus * theNucleus
static constexpr double MeV
Definition: G4SIunits.hh:214
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

Here is the caller graph for this function:

virtual G4double G4PionMinusField::GetCoeff ( )
inlinevirtual

Reimplemented from G4VNuclearField.

Definition at line 64 of file G4PionMinusField.hh.

64 { return theCoeff; }
G4double G4PionMinusField::GetField ( const G4ThreeVector aPosition)
virtual

Implements G4VNuclearField.

Definition at line 59 of file G4PionMinusField.cc.

60 {
61 // Field is 0 out of the nucleus!
62  if(aPosition.mag() >= radius) return 0.0;
63 
66  G4double pionMinusMass = G4PionMinus::PionMinus()->GetPDGMass();
68  G4double nucleusMass = Z*proton_mass_c2+(A-Z)*neutron_mass_c2+bindingEnergy;
69  G4double reducedMass = pionMinusMass*nucleusMass/(pionMinusMass+nucleusMass);
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+pionMinusMass/nucleonMass)*theCoeff*density + GetBarrier();
75 }
virtual G4int GetCharge()=0
virtual const G4VNuclearDensity * GetNuclearDensity() const =0
virtual G4int GetMassNumber()=0
int G4int
Definition: G4Types.hh:78
double A(double temperature)
float proton_mass_c2
Definition: hepunit.py:275
float neutron_mass_c2
Definition: hepunit.py:276
const G4double radius
G4double GetPDGMass() const
G4double GetDensity(const G4ThreeVector &aPosition) const
static G4double GetBindingEnergy(const G4int A, const G4int Z)
int G4lrint(double ad)
Definition: templates.hh:163
static G4PionMinus * PionMinus()
Definition: G4PionMinus.cc:98
G4V3DNucleus * theNucleus
virtual G4double GetBarrier()
static constexpr double pi
Definition: G4SIunits.hh:75
double G4double
Definition: G4Types.hh:76
G4double bindingEnergy(G4int A, G4int Z)
double mag() const

Here is the call graph for this function:


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