38 _ProtonMeanMultiplicity(0.0)
51 "G4StatMFMacroNucleon::CalcMeanMultiplicity: Temperature less or equal 0");
56 G4double lambda3 = ThermalWaveLenght*ThermalWaveLenght*ThermalWaveLenght;
58 static const G4double degeneracy = 2.0;
63 G4double exponent_proton = (mu+nu-Coulomb)/T;
66 if (exponent_neutron > 700.0) exponent_neutron = 700.0;
67 if (exponent_proton > 700.0) exponent_proton = 700.0;
69 _NeutronMeanMultiplicity =
70 (degeneracy*FreeVol/lambda3)*std::exp(exponent_neutron);
72 _ProtonMeanMultiplicity =
73 (degeneracy*FreeVol/lambda3)*std::exp(exponent_proton);
91 G4double lambda3 = ThermalWaveLenght*ThermalWaveLenght*ThermalWaveLenght;
94 if (_NeutronMeanMultiplicity > 0.0)
95 NeutronEntropy = _NeutronMeanMultiplicity*(5./2.+
96 std::log(2.0*static_cast<G4double>(
theA)*FreeVol/
97 (lambda3*_NeutronMeanMultiplicity)));
100 if (_ProtonMeanMultiplicity > 0.0)
101 ProtonEntropy = _ProtonMeanMultiplicity*(5./2.+
102 std::log(2.0*static_cast<G4double>(
theA)*FreeVol/
103 (lambda3*_ProtonMeanMultiplicity)));
105 return NeutronEntropy+ProtonEntropy;
G4double CalcMeanMultiplicity(const G4double FreeVol, const G4double mu, const G4double nu, const G4double T)
static G4double GetKappaCoulomb()
G4double CalcEnergy(const G4double T)
G4double CalcEntropy(const G4double T, const G4double FreeVol)
G4double _MeanMultiplicity