G4StatMFMacroMultiNucleon Class Reference

#include <G4StatMFMacroMultiNucleon.hh>

Inheritance diagram for G4StatMFMacroMultiNucleon:

Public Member Functions
	G4StatMFMacroMultiNucleon (const G4int Size)
	~G4StatMFMacroMultiNucleon ()
G4double	CalcMeanMultiplicity (const G4double FreeVol, const G4double mu, const G4double nu, const G4double T)
G4double	CalcZARatio (const G4double nu)
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 G4StatMFMacroMultiNucleon.hh.

Constructor & Destructor Documentation

G4StatMFMacroMultiNucleon::G4StatMFMacroMultiNucleon ( const G4int  Size )

inline

Definition at line 43 of file G4StatMFMacroMultiNucleon.hh.

: G4VStatMFMacroCluster(Size) {};

G4StatMFMacroMultiNucleon::~G4StatMFMacroMultiNucleon ( )

inline

Definition at line 46 of file G4StatMFMacroMultiNucleon.hh.

{};

Member Function Documentation

G4double G4StatMFMacroMultiNucleon::CalcEnergy ( const G4double  T )

virtual

Implements G4VStatMFMacroCluster.

Definition at line 126 of file G4StatMFMacroMultiNucleon.cc.

References G4VStatMFMacroCluster::_Energy, G4VStatMFMacroCluster::_InvLevelDensity, A23, G4StatMFParameters::Beta(), G4StatMFParameters::DBetaDT(), python.hepunit::elm_coupling, G4StatMFParameters::GetE0(), G4StatMFParameters::GetGamma0(), G4StatMFParameters::GetKappaCoulomb(), G4StatMFParameters::Getr0(), G4VStatMFMacroCluster::theA, and G4VStatMFMacroCluster::theZARatio.

{
    const G4double Coulomb = (3./5.)*(elm_coupling/G4StatMFParameters::Getr0())*
    (1.0 - 1.0/std::pow(1.0+G4StatMFParameters::GetKappaCoulomb(),1./3.));
    
    const G4double A23 = std::pow(static_cast<G4double>(theA),2./3.);

    // Volume term 
    G4double EVol = static_cast<G4double>(theA) * (T*T/_InvLevelDensity - G4StatMFParameters::GetE0());
    
    // Symmetry term
    G4double ESym = static_cast<G4double>(theA) * G4StatMFParameters::GetGamma0() *(1. - 2.* theZARatio) * (1. - 2.* theZARatio);
    
    // Surface term
    G4double ESurf = A23*(G4StatMFParameters::Beta(T) - T*G4StatMFParameters::DBetaDT(T));
 
    // Coulomb term
    G4double ECoul = Coulomb*A23*static_cast<G4double>(theA)*theZARatio*theZARatio;
    
    // Translational term
    G4double ETrans = (3./2.)*T;
    
       
    return _Energy = EVol + ESurf + ECoul + ETrans + ESym;
}

G4double G4StatMFMacroMultiNucleon::CalcEntropy ( const G4double  T, const G4double  FreeVol  )

virtual

Implements G4VStatMFMacroCluster.

Definition at line 153 of file G4StatMFMacroMultiNucleon.cc.

References G4VStatMFMacroCluster::_InvLevelDensity, G4VStatMFMacroCluster::_MeanMultiplicity, G4StatMFParameters::DBetaDT(), python.hepunit::fermi, and G4VStatMFMacroCluster::theA.

{
    const G4double ThermalWaveLenght = 16.15*fermi/std::sqrt(T);
    const G4double lambda3 = ThermalWaveLenght*ThermalWaveLenght*ThermalWaveLenght;

    G4double Entropy = 0.0;
    if (_MeanMultiplicity > 0.0) {
    // Volume term
    G4double SV = 2.0*static_cast<G4double>(theA)*T/_InvLevelDensity;
        
    // Surface term
    G4double SS = -G4StatMFParameters::DBetaDT(T)*std::pow(static_cast<G4double>(theA),2./3.);
        
    // Translational term
    G4double ST = (5./2.)+std::log(FreeVol * std::sqrt(static_cast<G4double>(theA)) * 
                  static_cast<G4double>(theA)/(lambda3*_MeanMultiplicity));
        
        
    Entropy = _MeanMultiplicity*(SV + SS + ST);
    }
                                
                                
    return Entropy;
}

G4double G4StatMFMacroMultiNucleon::CalcMeanMultiplicity ( const G4double  FreeVol, const G4double  mu, const G4double  nu, const G4double  T  )

virtual

Implements G4VStatMFMacroCluster.

Definition at line 83 of file G4StatMFMacroMultiNucleon.cc.

References G4VStatMFMacroCluster::_InvLevelDensity, G4VStatMFMacroCluster::_MeanMultiplicity, A23, G4StatMFParameters::Beta(), python.hepunit::elm_coupling, python.hepunit::fermi, G4StatMFParameters::GetE0(), G4StatMFParameters::GetGamma0(), G4StatMFParameters::GetKappaCoulomb(), G4StatMFParameters::Getr0(), G4INCL::Math::max(), G4VStatMFMacroCluster::theA, and G4VStatMFMacroCluster::theZARatio.

{
    const G4double ThermalWaveLenght = 16.15*fermi/std::sqrt(T);
    
    const G4double lambda3 = ThermalWaveLenght*ThermalWaveLenght*ThermalWaveLenght;
    
    const G4double A23 = std::pow(static_cast<G4double>(theA),2./3.);
    
    const G4double Coulomb = (3./5.)*(elm_coupling/G4StatMFParameters::Getr0())*
    (1.0 - 1.0/std::pow(1.0+G4StatMFParameters::GetKappaCoulomb(),1./3.));
    
    G4double exponent = (mu + nu*theZARatio+ G4StatMFParameters::GetE0() + T*T/_InvLevelDensity 
             - G4StatMFParameters::GetGamma0()*(1.0 - 2.0*theZARatio)*
             (1.0 - 2.0*theZARatio))*theA
    - G4StatMFParameters::Beta(T)*A23 - Coulomb*theZARatio*theZARatio*A23*theA;
    
    exponent /= T;
    
    if (exponent > 30.0) exponent = 30.0;
    
    _MeanMultiplicity = std::max((FreeVol * static_cast<G4double>(theA) * 
                    std::sqrt(static_cast<G4double>(theA))/lambda3) *
                   std::exp(exponent),1.0e-30);
    return _MeanMultiplicity;   
}

G4double G4StatMFMacroMultiNucleon::CalcZARatio ( const G4double  nu )

virtual

Implements G4VStatMFMacroCluster.

Definition at line 111 of file G4StatMFMacroMultiNucleon.cc.

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

{
    const G4double Coulomb = (3./5.)*(elm_coupling/G4StatMFParameters::Getr0())*
    (1.0 - 1.0/std::pow(1.0+G4StatMFParameters::GetKappaCoulomb(),1./3.));

    G4double den = 8.0*G4StatMFParameters::GetGamma0()+2.0*Coulomb*std::pow(static_cast<G4double>(theA),2./3.);
    G4double num = 4.0*G4StatMFParameters::GetGamma0()+nu;
    
    return theZARatio = num/den;
    

}

---

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

• G4StatMFMacroMultiNucleon.hh
• G4StatMFMacroMultiNucleon.cc