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

#include <G4StatMFMacroNucleon.hh>

Inheritance diagram for G4StatMFMacroNucleon:
G4VStatMFMacroCluster

Public Member Functions

 G4StatMFMacroNucleon ()
 
 ~G4StatMFMacroNucleon ()
 
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 37 of file G4StatMFMacroNucleon.hh.

Constructor & Destructor Documentation

G4StatMFMacroNucleon::G4StatMFMacroNucleon ( )

Definition at line 36 of file G4StatMFMacroNucleon.cc.

37  : G4VStatMFMacroCluster(1), _NeutronMeanMultiplicity(0.0),
38  _ProtonMeanMultiplicity(0.0)
39 {}
G4StatMFMacroNucleon::~G4StatMFMacroNucleon ( )

Definition at line 41 of file G4StatMFMacroNucleon.cc.

42 {}

Member Function Documentation

G4double G4StatMFMacroNucleon::CalcEnergy ( const G4double  T)
virtual

Implements G4VStatMFMacroCluster.

Definition at line 79 of file G4StatMFMacroNucleon.cc.

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

80 {
81  G4double Coulomb = (3./5.)*(elm_coupling/G4StatMFParameters::Getr0())*
82  (1.0 - 1.0/std::pow(1.0+G4StatMFParameters::GetKappaCoulomb(),1./3.));
83 
84  return _Energy = Coulomb * theZARatio * theZARatio + (3./2.) * T;
85 }
static G4double GetKappaCoulomb()
tuple elm_coupling
Definition: hepunit.py:286
static G4double Getr0()
double G4double
Definition: G4Types.hh:76
G4double G4StatMFMacroNucleon::CalcEntropy ( const G4double  T,
const G4double  FreeVol 
)
virtual

Implements G4VStatMFMacroCluster.

Definition at line 88 of file G4StatMFMacroNucleon.cc.

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

89 {
90  G4double ThermalWaveLenght = 16.15*fermi/std::sqrt(T);
91  G4double lambda3 = ThermalWaveLenght*ThermalWaveLenght*ThermalWaveLenght;
92 
93  G4double NeutronEntropy = 0.0;
94  if (_NeutronMeanMultiplicity > 0.0)
95  NeutronEntropy = _NeutronMeanMultiplicity*(5./2.+
96  std::log(2.0*static_cast<G4double>(theA)*FreeVol/
97  (lambda3*_NeutronMeanMultiplicity)));
98 
99  G4double ProtonEntropy = 0.0;
100  if (_ProtonMeanMultiplicity > 0.0)
101  ProtonEntropy = _ProtonMeanMultiplicity*(5./2.+
102  std::log(2.0*static_cast<G4double>(theA)*FreeVol/
103  (lambda3*_ProtonMeanMultiplicity)));
104 
105  return NeutronEntropy+ProtonEntropy;
106 }
double G4double
Definition: G4Types.hh:76
G4double G4StatMFMacroNucleon::CalcMeanMultiplicity ( const G4double  FreeVol,
const G4double  mu,
const G4double  nu,
const G4double  T 
)
virtual

Implements G4VStatMFMacroCluster.

Definition at line 45 of file G4StatMFMacroNucleon.cc.

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

48 {
49  if (T <= 0.0) {
50  throw G4HadronicException(__FILE__, __LINE__,
51  "G4StatMFMacroNucleon::CalcMeanMultiplicity: Temperature less or equal 0");
52  }
53 
54  G4double ThermalWaveLenght = 16.15*fermi/std::sqrt(T);
55 
56  G4double lambda3 = ThermalWaveLenght*ThermalWaveLenght*ThermalWaveLenght;
57 
58  static const G4double degeneracy = 2.0;
59 
60  G4double Coulomb = (3./5.)*(elm_coupling/G4StatMFParameters::Getr0())*
61  (1.0 - 1.0/std::pow(1.0+G4StatMFParameters::GetKappaCoulomb(),1./3.));
62 
63  G4double exponent_proton = (mu+nu-Coulomb)/T;
64  G4double exponent_neutron = mu/T;
65 
66  if (exponent_neutron > 700.0) exponent_neutron = 700.0;
67  if (exponent_proton > 700.0) exponent_proton = 700.0;
68 
69  _NeutronMeanMultiplicity =
70  (degeneracy*FreeVol/lambda3)*std::exp(exponent_neutron);
71 
72  _ProtonMeanMultiplicity =
73  (degeneracy*FreeVol/lambda3)*std::exp(exponent_proton);
74 
75  return _MeanMultiplicity = _NeutronMeanMultiplicity + _ProtonMeanMultiplicity;
76 }
static G4double GetKappaCoulomb()
tuple elm_coupling
Definition: hepunit.py:286
static G4double Getr0()
double G4double
Definition: G4Types.hh:76
G4double G4StatMFMacroNucleon::CalcZARatio ( const G4double  )
inlinevirtual

Implements G4VStatMFMacroCluster.

Definition at line 48 of file G4StatMFMacroNucleon.hh.

References G4VStatMFMacroCluster::theZARatio.

49  {
50  theZARatio = 0.0;
51  if (_ProtonMeanMultiplicity+_NeutronMeanMultiplicity > 0.0) {
52  theZARatio = _ProtonMeanMultiplicity/
53  (_ProtonMeanMultiplicity+_NeutronMeanMultiplicity);
54  }
55  return theZARatio;
56  }

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