Geant4.10
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions
G4StatMFMacroBiNucleon Class Reference

#include <G4StatMFMacroBiNucleon.hh>

Inheritance diagram for G4StatMFMacroBiNucleon:
G4VStatMFMacroCluster

Public Member Functions

 G4StatMFMacroBiNucleon ()
 
 ~G4StatMFMacroBiNucleon ()
 
G4double CalcMeanMultiplicity (const G4double FreeVol, const G4double mu, const G4double nu, const G4double T)
 
G4double CalcZARatio (const G4double)
 
G4double CalcEnergy (const G4double T)
 
G4double CalcEntropy (const G4double T, const G4double FreeVol)
 
- Public Member Functions inherited from G4VStatMFMacroCluster
 G4VStatMFMacroCluster (const G4int Size)
 
virtual ~G4VStatMFMacroCluster ()
 
G4bool operator== (const G4VStatMFMacroCluster &right) const
 
G4bool operator!= (const G4VStatMFMacroCluster &right) const
 
G4double GetMeanMultiplicity (void) const
 
G4double GetInvLevelDensity (void) const
 
void SetZARatio (const G4double value)
 
G4double GetZARatio (void) const
 
void SetSize (const G4double value)
 
G4double GetSize (void) const
 

Additional Inherited Members

- Protected Attributes inherited from G4VStatMFMacroCluster
G4int theA
 
G4double _InvLevelDensity
 
G4double _Entropy
 
G4double theZARatio
 
G4double _MeanMultiplicity
 
G4double _Energy
 

Detailed Description

Definition at line 38 of file G4StatMFMacroBiNucleon.hh.

Constructor & Destructor Documentation

G4StatMFMacroBiNucleon::G4StatMFMacroBiNucleon ( )
inline
G4StatMFMacroBiNucleon::~G4StatMFMacroBiNucleon ( )
inline

Definition at line 46 of file G4StatMFMacroBiNucleon.hh.

46 {};

Member Function Documentation

G4double G4StatMFMacroBiNucleon::CalcEnergy ( const G4double  T)
virtual

Implements G4VStatMFMacroCluster.

Definition at line 87 of file G4StatMFMacroBiNucleon.cc.

References G4VStatMFMacroCluster::_Energy, python.hepunit::elm_coupling, G4NucleiProperties::GetBindingEnergy(), G4StatMFParameters::GetKappaCoulomb(), G4StatMFParameters::Getr0(), G4VStatMFMacroCluster::theA, and G4VStatMFMacroCluster::theZARatio.

88 {
89  const G4double Coulomb = (3./5.)*(elm_coupling/G4StatMFParameters::Getr0())*
90  (1.0 - 1.0/std::pow(1.0+G4StatMFParameters::GetKappaCoulomb(),1./3.));
91 
93  Coulomb * theZARatio * theZARatio * std::pow(G4double(theA),5./3.) +
94  (3./2.) * T;
95 
96  return _Energy;
97 }
static G4double GetKappaCoulomb()
tuple elm_coupling
Definition: hepunit.py:286
static G4double Getr0()
static G4double GetBindingEnergy(const G4int A, const G4int Z)
double G4double
Definition: G4Types.hh:76
G4double G4StatMFMacroBiNucleon::CalcEntropy ( const G4double  T,
const G4double  FreeVol 
)
virtual

Implements G4VStatMFMacroCluster.

Definition at line 101 of file G4StatMFMacroBiNucleon.cc.

References G4VStatMFMacroCluster::_MeanMultiplicity, python.hepunit::fermi, and G4VStatMFMacroCluster::theA.

102 {
103  const G4double ThermalWaveLenght = 16.15*fermi/std::sqrt(T);
104  const G4double lambda3 = ThermalWaveLenght*ThermalWaveLenght*ThermalWaveLenght;
105 
106  G4double Entropy = 0.0;
107  if (_MeanMultiplicity > 0.0)
108  // Is this formula correct?
109  Entropy = _MeanMultiplicity*(5./2.+
110  std::log(3.0*static_cast<G4double>(theA)*
111  std::sqrt(static_cast<G4double>(theA))*FreeVol/
112  (lambda3*_MeanMultiplicity)));
113 
114 
115  return Entropy;
116 }
double G4double
Definition: G4Types.hh:76
G4double G4StatMFMacroBiNucleon::CalcMeanMultiplicity ( const G4double  FreeVol,
const G4double  mu,
const G4double  nu,
const G4double  T 
)
virtual

Implements G4VStatMFMacroCluster.

Definition at line 60 of file G4StatMFMacroBiNucleon.cc.

References G4VStatMFMacroCluster::_MeanMultiplicity, python.hepunit::elm_coupling, python.hepunit::fermi, G4NucleiProperties::GetBindingEnergy(), G4StatMFParameters::GetKappaCoulomb(), G4StatMFParameters::Getr0(), G4VStatMFMacroCluster::theA, and G4VStatMFMacroCluster::theZARatio.

62 {
63  const G4double ThermalWaveLenght = 16.15*fermi/std::sqrt(T);
64 
65  const G4double lambda3 = ThermalWaveLenght*ThermalWaveLenght*ThermalWaveLenght;
66 
67  const G4double degeneracy = 3.0;
68 
69  const G4double Coulomb = (3./5.)*(elm_coupling/G4StatMFParameters::Getr0())*
70  (1.0 - 1.0/std::pow(1.0+G4StatMFParameters::GetKappaCoulomb(),1./3.));
71 
72  const G4double BindingE = G4NucleiProperties::GetBindingEnergy(theA,1); //old value was 2.796*MeV
73  G4double exponent = (BindingE + theA*(mu+nu*theZARatio) -
74  Coulomb*theZARatio*theZARatio*std::pow(G4double(theA),5./3.))/T;
75 
76  // To avoid numerical problems
77  if (exponent < -700.0) exponent = -700.0;
78  else if (exponent > 700.0) exponent = 700.0;
79 
80  _MeanMultiplicity = (degeneracy*FreeVol*static_cast<G4double>(theA)*std::sqrt(static_cast<G4double>(theA))/lambda3)*
81  std::exp(exponent);
82 
83  return _MeanMultiplicity;
84 }
static G4double GetKappaCoulomb()
tuple elm_coupling
Definition: hepunit.py:286
static G4double Getr0()
static G4double GetBindingEnergy(const G4int A, const G4int Z)
double G4double
Definition: G4Types.hh:76
G4double G4StatMFMacroBiNucleon::CalcZARatio ( const G4double  )
inlinevirtual

Implements G4VStatMFMacroCluster.

Definition at line 66 of file G4StatMFMacroBiNucleon.hh.

References G4VStatMFMacroCluster::theZARatio.

66 {return theZARatio = 0.5;}

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