45 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMacroMultiplicity::operator= meant to not be accessable");
51 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMacroMultiplicity::operator== meant to not be accessable");
58 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMacroMultiplicity::operator!= meant to not be accessable");
81 (5.0/3.0)*CP*ZA5*ZA5*std::pow(5.,2./3.) -
88 G4double ChemPb = ChemPa - 0.5*std::abs(ChemPa);
100 G4int iterations = 0;
101 while (fChemPa*fChemPb > 0.0 && iterations < 100)
103 if (std::abs(fChemPa) <= std::abs(fChemPb))
105 ChemPa += 0.6*(ChemPa-ChemPb);
111 ChemPb += 0.6*(ChemPb-ChemPa);
117 if (fChemPa*fChemPb > 0.0)
119 G4cerr <<
"G4StatMFMacroMultiplicity:"<<
" ChemPa="<<ChemPa<<
" ChemPb="<<ChemPb<<
G4endl;
120 G4cerr <<
"G4StatMFMacroMultiplicity:"<<
" fChemPa="<<fChemPa<<
" fChemPb="<<fChemPb<<
G4endl;
121 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMacroMultiplicity::CalcChemicalPotentialMu: I couldn't bracket the root.");
123 else if (fChemPa*fChemPb < 0.0 && std::abs(ChemPa-ChemPb) > intervalWidth)
128 if (!theSolver->
Brent(*
this))
130 G4cerr <<
"G4StatMFMacroMultiplicity:"<<
" ChemPa="<<ChemPa<<
" ChemPb="<<ChemPb<<
G4endl;
131 G4cerr <<
"G4StatMFMacroMultiplicity:"<<
" fChemPa="<<fChemPa<<
" fChemPb="<<fChemPb<<
G4endl;
132 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMacroMultiplicity::CalcChemicalPotentialMu: I couldn't find the root.");
158 for (std::vector<G4VStatMFMacroCluster*>::iterator i =
_theClusters->begin();
162 MeanA += multip*
static_cast<G4double>(n++);
static G4double GetGamma0()
G4double _ChemPotentialNu
static G4double GetKappaCoulomb()
G4StatMFMacroMultiplicity & operator=(const G4StatMFMacroMultiplicity &right)
G4double _MeanTemperature
void SetIntervalLimits(const G4double Limit1, const G4double Limit2)
G4double _MeanMultiplicity
G4double CalcChemicalPotentialMu(void)
G4bool Brent(Function &theFunction)
G4bool operator!=(const G4StatMFMacroMultiplicity &right) const
G4bool operator==(const G4StatMFMacroMultiplicity &right) const
std::vector< G4VStatMFMacroCluster * > * _theClusters
G4double GetRoot(void) const
G4double CalcMeanA(const G4double mu)
static G4double Beta(G4double T)
G4double _ChemPotentialMu
G4GLOB_DLL std::ostream G4cerr
G4double operator()(const G4double mu)