Geant4  10.02.p03
G4SigmaPlusField Class Reference

#include <G4SigmaPlusField.hh>

Inheritance diagram for G4SigmaPlusField:
Collaboration diagram for G4SigmaPlusField:

Public Member Functions

 G4SigmaPlusField (G4V3DNucleus *nucleus, G4double coeff=0.36 *CLHEP::fermi)
 
virtual ~G4SigmaPlusField ()
 
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)
 

Private Member Functions

 G4SigmaPlusField (const G4SigmaPlusField &right)
 
const G4SigmaPlusFieldoperator= (const G4SigmaPlusField &right)
 
int operator== (const G4SigmaPlusField &right) const
 
int operator!= (const G4SigmaPlusField &right) const
 

Private Attributes

G4double theCoeff
 

Additional Inherited Members

- Protected Attributes inherited from G4VNuclearField
G4V3DNucleustheNucleus
 
const G4double radius
 

Detailed Description

Definition at line 35 of file G4SigmaPlusField.hh.

Constructor & Destructor Documentation

◆ G4SigmaPlusField() [1/2]

G4SigmaPlusField::G4SigmaPlusField ( G4V3DNucleus nucleus,
G4double  coeff = 0.36*CLHEP::fermi 
)

Definition at line 49 of file G4SigmaPlusField.cc.

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

◆ ~G4SigmaPlusField()

G4SigmaPlusField::~G4SigmaPlusField ( )
virtual

Definition at line 56 of file G4SigmaPlusField.cc.

57 { }

◆ G4SigmaPlusField() [2/2]

G4SigmaPlusField::G4SigmaPlusField ( const G4SigmaPlusField right)
private

Member Function Documentation

◆ GetBarrier()

G4double G4SigmaPlusField::GetBarrier ( )
virtual

Implements G4VNuclearField.

Definition at line 78 of file G4SigmaPlusField.cc.

79 {
82  G4double coulombBarrier = (1.44/1.14) * MeV * Z / (1.0 + G4Pow::GetInstance()->A13(A));
83  return coulombBarrier;
84 }
static G4Pow * GetInstance()
Definition: G4Pow.cc:55
static const double MeV
Definition: G4SIunits.hh:211
virtual G4int GetCharge()=0
virtual G4int GetMassNumber()=0
double A(double temperature)
Float_t Z
G4double A13(G4double A) const
Definition: G4Pow.hh:132
G4V3DNucleus * theNucleus
double G4double
Definition: G4Types.hh:76
Here is the call graph for this function:
Here is the caller graph for this function:

◆ GetCoeff()

virtual G4double G4SigmaPlusField::GetCoeff ( )
inlinevirtual

Reimplemented from G4VNuclearField.

Definition at line 50 of file G4SigmaPlusField.hh.

50 { return theCoeff; }

◆ GetField()

G4double G4SigmaPlusField::GetField ( const G4ThreeVector aPosition)
virtual

Implements G4VNuclearField.

Definition at line 59 of file G4SigmaPlusField.cc.

60 {
61 // Field is 0 out of the nucleus!
62  if(aPosition.mag() >= radius) return 0.0;
63 
64  G4double sigmaPlusMass = G4SigmaPlus::SigmaPlus()->GetPDGMass();
65 
68  G4double bindingEnergy = G4NucleiProperties::GetBindingEnergy(A, Z);
69  G4double nucleusMass = Z*proton_mass_c2+(A-Z)*neutron_mass_c2+bindingEnergy;
70  G4double reducedMass = sigmaPlusMass*nucleusMass/(sigmaPlusMass+nucleusMass);
71 
72  const G4VNuclearDensity * nuclearDensity=theNucleus->GetNuclearDensity();
73  G4double density = nuclearDensity->GetDensity(aPosition);
74 
75  return -2.*pi*hbarc*hbarc/reducedMass*(2.0)*theCoeff*density+GetBarrier();
76 }
virtual G4int GetCharge()=0
virtual const G4VNuclearDensity * GetNuclearDensity() const =0
virtual G4int GetMassNumber()=0
G4double GetDensity(const G4ThreeVector &aPosition) const
int G4int
Definition: G4Types.hh:78
G4double density
Definition: TRTMaterials.hh:39
double A(double temperature)
float hbarc
Definition: hepunit.py:265
Float_t Z
double mag() const
float proton_mass_c2
Definition: hepunit.py:275
float neutron_mass_c2
Definition: hepunit.py:276
const G4double radius
virtual G4double GetBarrier()
static const double pi
Definition: G4SIunits.hh:74
static G4double GetBindingEnergy(const G4int A, const G4int Z)
G4V3DNucleus * theNucleus
static G4SigmaPlus * SigmaPlus()
Definition: G4SigmaPlus.cc:108
double G4double
Definition: G4Types.hh:76
Here is the call graph for this function:

◆ operator!=()

int G4SigmaPlusField::operator!= ( const G4SigmaPlusField right) const
private

◆ operator=()

const G4SigmaPlusField& G4SigmaPlusField::operator= ( const G4SigmaPlusField right)
private

◆ operator==()

int G4SigmaPlusField::operator== ( const G4SigmaPlusField right) const
private

Member Data Documentation

◆ theCoeff

G4double G4SigmaPlusField::theCoeff
private

Definition at line 53 of file G4SigmaPlusField.hh.


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