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

#include <G4ProtonField.hh>

Inheritance diagram for G4ProtonField:
G4VNuclearField

Public Member Functions

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

Protected Attributes

const G4double radius
 
G4V3DNucleus * theNucleus
 

Private Member Functions

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

Private Attributes

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

Detailed Description

Definition at line 48 of file G4ProtonField.hh.

Constructor & Destructor Documentation

◆ G4ProtonField() [1/2]

G4ProtonField::G4ProtonField ( G4V3DNucleus *  nucleus)

Definition at line 47 of file G4ProtonField.cc.

47 :
48 G4VNuclearField(aNucleus), theDensity(theNucleus->GetNuclearDensity())
49{
50 theA = theNucleus->GetMassNumber();
51 theZ = theNucleus->GetCharge();
52 theBarrier = GetBarrier();
53 theRadius = 2.*theNucleus->GetOuterRadius();
54 theFermi.Init(theA, theZ);
55 for (G4double aR=0.;aR<theRadius; aR+=0.3*fermi)
56 {
57 G4ThreeVector aPosition(0,0,aR);
58 G4double density = GetDensity(aPosition);
59 G4double fermiMom = GetFermiMomentum(density);
60 theFermiMomBuffer.push_back(fermiMom);
61 }
62 {
63 G4ThreeVector aPosition(0,0,theRadius);
64 G4double density = GetDensity(aPosition);
65 G4double fermiMom = GetFermiMomentum(density);
66 theFermiMomBuffer.push_back(fermiMom);
67 }
68 {
69 G4ThreeVector aPosition(0,0,theRadius+0.001*fermi);
70 theFermiMomBuffer.push_back(0);
71 }
72 {
73 G4ThreeVector aPosition(0,0,1.*m);
74 theFermiMomBuffer.push_back(0);
75 }
76}
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
G4ProtonField::GetFermiMomentum
G4double GetFermiMomentum(const G4double aDensity)
Definition: G4ProtonField.hh:72
G4ProtonField::GetDensity
G4double GetDensity(const G4ThreeVector &aPosition)
Definition: G4ProtonField.hh:68
G4ProtonField::theDensity
const G4VNuclearDensity * theDensity
Definition: G4ProtonField.hh:82
G4ProtonField::theBarrier
G4double theBarrier
Definition: G4ProtonField.hh:79
G4ProtonField::GetBarrier
virtual G4double GetBarrier()
Definition: G4ProtonField.cc:98
G4ProtonField::theFermiMomBuffer
std::vector< G4double > theFermiMomBuffer
Definition: G4ProtonField.hh:84
G4ProtonField::theFermi
G4FermiMomentum theFermi
Definition: G4ProtonField.hh:81
G4ProtonField::theRadius
G4double theRadius
Definition: G4ProtonField.hh:80
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, GetBarrier(), G4V3DNucleus::GetCharge(), GetDensity(), GetFermiMomentum(), G4V3DNucleus::GetMassNumber(), G4V3DNucleus::GetOuterRadius(), G4FermiMomentum::Init(), m, theA, theBarrier, theFermi, theFermiMomBuffer, G4VNuclearField::theNucleus, theRadius, and theZ.

◆ ~G4ProtonField()

G4ProtonField::~G4ProtonField ( )
virtual

Definition at line 79 of file G4ProtonField.cc.

80{ }

◆ G4ProtonField() [2/2]

G4ProtonField::G4ProtonField ( const G4ProtonField &  right)
private

Member Function Documentation

◆ GetBarrier()

G4double G4ProtonField::GetBarrier ( )
virtual

Implements G4VNuclearField.

Definition at line 98 of file G4ProtonField.cc.

99{
100 G4double coulombBarrier = (1.44/1.14) * MeV * theZ / (1.0 + G4Pow::GetInstance()->Z13(theA));
101//GF G4double bindingEnergy = G4NucleiPropertiesTable::GetBindingEnergy(Z, A);
102 G4double bindingEnergy =0;
103/*
104 * G4cout << " coulombBarrier/bindingEnergy : "
105 * << coulombBarrier << " /" << bindingEnergy << G4endl;
106 */
107 return bindingEnergy/theA+coulombBarrier;
108}
MeV
static constexpr double MeV
Definition: G4SIunits.hh:200
G4Pow::GetInstance
static G4Pow * GetInstance()
Definition: G4Pow.cc:41
G4Pow::Z13
G4double Z13(G4int Z) const
Definition: G4Pow.hh:123
G4InuclSpecialFunctions::bindingEnergy
G4double bindingEnergy(G4int A, G4int Z)
Definition: bindingEnergy.cc:36

References G4InuclSpecialFunctions::bindingEnergy(), G4Pow::GetInstance(), MeV, theA, theZ, and G4Pow::Z13().

Referenced by G4ProtonField().

◆ 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 G4ProtonField::GetDensity ( const G4ThreeVector &  aPosition)
inlineprivate

Definition at line 68 of file G4ProtonField.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 G4ProtonField().

◆ GetFermiMomentum()

G4double G4ProtonField::GetFermiMomentum ( const G4double  aDensity)
inlineprivate

Definition at line 72 of file G4ProtonField.hh.

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

References G4FermiMomentum::GetFermiMomentum(), and theFermi.

Referenced by G4ProtonField().

◆ GetField()

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

Implements G4VNuclearField.

Definition at line 82 of file G4ProtonField.cc.

83{
84//G4cout << " Fermi Potential " << (fermiMom*fermiMom)/(2*proton_mass_c2) <<G4endl;
85 G4double x = aPosition.mag();
86 unsigned int index = static_cast<unsigned int>(x/(0.3*fermi));
87 if((index+2) > theFermiMomBuffer.size()) return theFermiMomBuffer.back();
88 G4double y1 = theFermiMomBuffer[index];
89 G4double y2 = theFermiMomBuffer[index+1];
90 G4double x1 = (0.3*fermi)*index;
91 G4double x2 = (0.3*fermi)*(index+1);
92 G4double fermiMom = y1 + (x-x1)*(y2-y1)/(x2-x1);
93 G4double y = -1*(fermiMom*fermiMom)/(2*proton_mass_c2)+theBarrier;
94// G4cout <<" Protonfield test "<<index<<" "<< x1<<" "<<y1<<" "<<x2<<" "<<y2<<" "<<x<<" "<<y<<" "<<theBarrier<<G4endl;
95 return y;
96}
CLHEP::Hep3Vector::mag
double mag() const
source.hepunit.proton_mass_c2
float proton_mass_c2
Definition: hepunit.py:274

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

◆ operator!=()

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

◆ operator=()

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

◆ operator==()

G4bool G4ProtonField::operator== ( const G4ProtonField &  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 G4ProtonField::theA
private

Definition at line 77 of file G4ProtonField.hh.

Referenced by G4ProtonField(), and GetBarrier().

◆ theBarrier

G4double G4ProtonField::theBarrier
private

Definition at line 79 of file G4ProtonField.hh.

Referenced by G4ProtonField(), and GetField().

◆ theDensity

const G4VNuclearDensity* G4ProtonField::theDensity
private

Definition at line 82 of file G4ProtonField.hh.

Referenced by GetDensity().

◆ theFermi

G4FermiMomentum G4ProtonField::theFermi
private

Definition at line 81 of file G4ProtonField.hh.

Referenced by G4ProtonField(), and GetFermiMomentum().

◆ theFermiMomBuffer

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

Definition at line 84 of file G4ProtonField.hh.

Referenced by G4ProtonField(), and GetField().

◆ theNucleus

G4V3DNucleus* G4VNuclearField::theNucleus
protectedinherited

Definition at line 47 of file G4VNuclearField.hh.

Referenced by G4NeutronField::G4NeutronField(), 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().

◆ theRadius

G4double G4ProtonField::theRadius
private

Definition at line 80 of file G4ProtonField.hh.

Referenced by G4ProtonField().

◆ theZ

G4int G4ProtonField::theZ
private

Definition at line 78 of file G4ProtonField.hh.

Referenced by G4ProtonField(), and GetBarrier().

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