Geant4-11
Public Member Functions | Protected Attributes | Private Member Functions | Private Attributes
G4NeutronField Class Reference

#include <G4NeutronField.hh>

Inheritance diagram for G4NeutronField:
G4VNuclearField

Public Member Functions

 G4NeutronField (G4V3DNucleus *nucleus)
 
virtual G4double GetBarrier ()
 
virtual G4double GetCoeff ()
 
virtual G4double GetField (const G4ThreeVector &aPosition)
 
void SetNucleus (G4V3DNucleus *aNucleus)
 
virtual ~G4NeutronField ()
 

Protected Attributes

const G4double radius
 
G4V3DNucleus * theNucleus
 

Private Member Functions

 G4NeutronField (const G4NeutronField &right)
 
G4double GetDensity (const G4ThreeVector &aPosition)
 
G4double GetFermiMomentum (const G4double aDensity)
 
G4bool operator!= (const G4NeutronField &right) const
 
const G4NeutronField & operator= (const G4NeutronField &right)
 
G4bool operator== (const G4NeutronField &right) const
 

Private Attributes

G4int theA
 
const G4VNuclearDensity * theDensity
 
G4FermiMomentum theFermi
 
std::vector< G4double > theFermiMomBuffer
 
G4double theR
 
G4int theZ
 

Detailed Description

Definition at line 48 of file G4NeutronField.hh.

Constructor & Destructor Documentation

◆ G4NeutronField() [1/2]

G4NeutronField::G4NeutronField ( G4V3DNucleus *  nucleus)

Definition at line 45 of file G4NeutronField.cc.

45 :
46 G4VNuclearField(aNucleus), theDensity(theNucleus->GetNuclearDensity())
47{
48 theA = theNucleus->GetMassNumber();
49 theZ = theNucleus->GetCharge();
50 theFermi.Init(theA, theZ);
51 theR = 2.*theNucleus->GetOuterRadius();
52 for (G4double aR=0.;aR<theR; aR+=0.3*fermi)
53 {
54 G4ThreeVector aPosition(0,0,aR);
55 G4double density = GetDensity(aPosition);
56 G4double fermiMom = GetFermiMomentum(density);
57 theFermiMomBuffer.push_back(fermiMom);
58 }
59 {
60 G4ThreeVector aPosition(0,0,theR);
61 G4double density = GetDensity(aPosition);
62 G4double fermiMom = GetFermiMomentum(density);
63 theFermiMomBuffer.push_back(fermiMom);
64 }
65 {
66 G4ThreeVector aPosition(0,0,theR+0.001*fermi);
67 theFermiMomBuffer.push_back(0);
68 }
69 {
70 G4ThreeVector aPosition(0,0,1.*m);
71 theFermiMomBuffer.push_back(0);
72 }
73}
m
static constexpr double m
Definition: G4SIunits.hh:109
fermi
static constexpr double fermi
Definition: G4SIunits.hh:83
G4double
double G4double
Definition: G4Types.hh:83
G4FermiMomentum::Init
void Init(G4int anA, G4int aZ)
Definition: G4FermiMomentum.hh:43
G4NeutronField::theDensity
const G4VNuclearDensity * theDensity
Definition: G4NeutronField.hh:81
G4NeutronField::GetFermiMomentum
G4double GetFermiMomentum(const G4double aDensity)
Definition: G4NeutronField.hh:73
G4NeutronField::theFermiMomBuffer
std::vector< G4double > theFermiMomBuffer
Definition: G4NeutronField.hh:84
G4NeutronField::theFermi
G4FermiMomentum theFermi
Definition: G4NeutronField.hh:78
G4NeutronField::GetDensity
G4double GetDensity(const G4ThreeVector &aPosition)
Definition: G4NeutronField.hh:68
G4V3DNucleus::GetOuterRadius
virtual G4double GetOuterRadius()=0
G4V3DNucleus::GetNuclearDensity
virtual const G4VNuclearDensity * GetNuclearDensity() const =0
G4V3DNucleus::GetCharge
virtual G4int GetCharge()=0
G4V3DNucleus::GetMassNumber
virtual G4int GetMassNumber()=0
G4VNuclearField::theNucleus
G4V3DNucleus * theNucleus
Definition: G4VNuclearField.hh:47
G4VNuclearField::G4VNuclearField
G4VNuclearField(G4V3DNucleus *aNucleus=0)
Definition: G4VNuclearField.cc:42

References fermi, G4V3DNucleus::GetCharge(), GetDensity(), GetFermiMomentum(), G4V3DNucleus::GetMassNumber(), G4V3DNucleus::GetOuterRadius(), G4FermiMomentum::Init(), m, theA, theFermi, theFermiMomBuffer, G4VNuclearField::theNucleus, theR, and theZ.

◆ ~G4NeutronField()

G4NeutronField::~G4NeutronField ( )
virtual

Definition at line 75 of file G4NeutronField.cc.

76{ }

◆ G4NeutronField() [2/2]

G4NeutronField::G4NeutronField ( const G4NeutronField &  right)
private

Member Function Documentation

◆ GetBarrier()

G4double G4NeutronField::GetBarrier ( )
virtual

Implements G4VNuclearField.

Definition at line 91 of file G4NeutronField.cc.

92{
93/*
94 * G4double A = theNucleus->GetMassNumber();
95 * G4double Z = theNucleus->GetCharge();
96 *
97 * return G4NucleiPropertiesTable::GetBindingEnergy(Z, A)/A;
98 */
99 return 0.;
100}

◆ GetCoeff()

virtual G4double G4VNuclearField::GetCoeff ( )
inlinevirtualinherited

Reimplemented in G4AntiProtonField, G4KaonMinusField, G4KaonPlusField, G4KaonZeroField, G4PionMinusField, G4PionPlusField, G4PionZeroField, G4SigmaMinusField, G4SigmaPlusField, and G4SigmaZeroField.

Definition at line 44 of file G4VNuclearField.hh.

44{ return 0; }

◆ GetDensity()

G4double G4NeutronField::GetDensity ( const G4ThreeVector &  aPosition)
inlineprivate

Definition at line 68 of file G4NeutronField.hh.

69 {
70 return theDensity->GetDensity(aPosition);
71 }
G4VNuclearDensity::GetDensity
G4double GetDensity(const G4ThreeVector &aPosition) const
Definition: G4VNuclearDensity.hh:42

References G4VNuclearDensity::GetDensity(), and theDensity.

Referenced by G4NeutronField().

◆ GetFermiMomentum()

G4double G4NeutronField::GetFermiMomentum ( const G4double  aDensity)
inlineprivate

Definition at line 73 of file G4NeutronField.hh.

74 {
75 return theFermi.GetFermiMomentum(aDensity);
76 }
G4FermiMomentum::GetFermiMomentum
G4double GetFermiMomentum(G4double density)
Definition: G4FermiMomentum.hh:45

References G4FermiMomentum::GetFermiMomentum(), and theFermi.

Referenced by G4NeutronField().

◆ GetField()

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

Implements G4VNuclearField.

Definition at line 78 of file G4NeutronField.cc.

79{
80 G4double x = aPosition.mag();
81 unsigned int index = static_cast<unsigned int>(x/(0.3*fermi));
82 if( (index+2) > theFermiMomBuffer.size()) return theFermiMomBuffer.back();
83 G4double y1 = theFermiMomBuffer[index];
84 G4double y2 = theFermiMomBuffer[index+1];
85 G4double x1 = (0.3*fermi)*index;
86 G4double x2 = (0.3*fermi)*(index+1);
87 G4double fermiMom = y1 + (x-x1)*(y2-y1)/(x2-x1);
88 return -1*(fermiMom*fermiMom)/(2*neutron_mass_c2);
89}
CLHEP::Hep3Vector::mag
double mag() const
source.hepunit.neutron_mass_c2
float neutron_mass_c2
Definition: hepunit.py:275

References fermi, CLHEP::Hep3Vector::mag(), source.hepunit::neutron_mass_c2, and theFermiMomBuffer.

◆ operator!=()

G4bool G4NeutronField::operator!= ( const G4NeutronField &  right) const
private

◆ operator=()

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

◆ operator==()

G4bool G4NeutronField::operator== ( const G4NeutronField &  right) const
private

◆ SetNucleus()

void G4VNuclearField::SetNucleus ( G4V3DNucleus *  aNucleus)
inlineinherited

Definition at line 61 of file G4VNuclearField.hh.

62{
63 theNucleus = aNucleus;
64}

References G4VNuclearField::theNucleus.

Field Documentation

◆ radius

const G4double G4VNuclearField::radius
protectedinherited

Definition at line 48 of file G4VNuclearField.hh.

Referenced by G4AntiProtonField::GetField(), G4KaonMinusField::GetField(), G4KaonPlusField::GetField(), G4KaonZeroField::GetField(), G4PionMinusField::GetField(), G4PionPlusField::GetField(), G4PionZeroField::GetField(), G4SigmaMinusField::GetField(), G4SigmaPlusField::GetField(), and G4SigmaZeroField::GetField().

◆ theA

G4int G4NeutronField::theA
private

Definition at line 79 of file G4NeutronField.hh.

Referenced by G4NeutronField().

◆ theDensity

const G4VNuclearDensity* G4NeutronField::theDensity
private

Definition at line 81 of file G4NeutronField.hh.

Referenced by GetDensity().

◆ theFermi

G4FermiMomentum G4NeutronField::theFermi
private

Definition at line 78 of file G4NeutronField.hh.

Referenced by G4NeutronField(), and GetFermiMomentum().

◆ theFermiMomBuffer

std::vector<G4double> G4NeutronField::theFermiMomBuffer
private

Definition at line 84 of file G4NeutronField.hh.

Referenced by G4NeutronField(), and GetField().

◆ theNucleus

G4V3DNucleus* G4VNuclearField::theNucleus
protectedinherited

Definition at line 47 of file G4VNuclearField.hh.

Referenced by G4NeutronField(), G4ProtonField::G4ProtonField(), G4AntiProtonField::GetBarrier(), G4KaonMinusField::GetBarrier(), G4KaonPlusField::GetBarrier(), G4PionMinusField::GetBarrier(), G4PionPlusField::GetBarrier(), G4SigmaMinusField::GetBarrier(), G4SigmaPlusField::GetBarrier(), G4AntiProtonField::GetField(), G4KaonMinusField::GetField(), G4KaonPlusField::GetField(), G4KaonZeroField::GetField(), G4PionMinusField::GetField(), G4PionPlusField::GetField(), G4PionZeroField::GetField(), G4SigmaMinusField::GetField(), G4SigmaPlusField::GetField(), G4SigmaZeroField::GetField(), and G4VNuclearField::SetNucleus().

◆ theR

G4double G4NeutronField::theR
private

Definition at line 82 of file G4NeutronField.hh.

Referenced by G4NeutronField().

◆ theZ

G4int G4NeutronField::theZ
private

Definition at line 80 of file G4NeutronField.hh.

Referenced by G4NeutronField().

The documentation for this class was generated from the following files:
Generated by doxygen 1.9.3