39 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMacroChemicalPotential::operator= meant to not be accessable");
45 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMacroChemicalPotential::operator== meant to not be accessable");
52 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMacroChemicalPotential::operator!= meant to not be accessable");
76 if (fChemPa*fChemPb > 0.0) {
80 ChemPb -= 1.5*std::abs(ChemPb-ChemPa);
82 }
while (fChemPb < 0.0);
85 ChemPb += 1.5*std::abs(ChemPb-ChemPa);
87 }
while (fChemPb > 0.0);
95 if (!theSolver->
Brent(*
this)){
96 G4cerr <<
"G4StatMFMacroChemicalPotential:"<<
" ChemPa="<<ChemPa<<
" ChemPb="<<ChemPb<<
G4endl;
97 G4cerr <<
"G4StatMFMacroChemicalPotential:"<<
" fChemPa="<<fChemPa<<
" fChemPb="<<fChemPb<<
G4endl;
98 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMacroChemicalPotential::CalcChemicalPotentialNu: I couldn't find the root.");
109 std::vector<G4VStatMFMacroCluster*>::iterator i;
112 (*i)->CalcZARatio(nu);
124 MeanZ +=
static_cast<G4double>(n++) *
126 (*i)->GetMeanMultiplicity();
static G4double GetGamma0()
G4double operator()(const G4double nu)
static G4double GetKappaCoulomb()
G4double _MeanTemperature
void SetIntervalLimits(const G4double Limit1, const G4double Limit2)
G4double CalcChemicalPotentialMu(void)
G4double _ChemPotentialMu
std::vector< G4VStatMFMacroCluster * > * _theClusters
G4bool Brent(Function &theFunction)
G4bool operator!=(const G4StatMFMacroChemicalPotential &right) const
G4double CalcMeanZ(const G4double nu)
G4bool operator==(const G4StatMFMacroChemicalPotential &right) const
G4double GetMeanMultiplicity(void) const
G4double GetRoot(void) const
void CalcChemicalPotentialMu(const G4double nu)
G4double _MeanMultiplicity
G4double _ChemPotentialNu
G4double CalcChemicalPotentialNu(void)
G4StatMFMacroChemicalPotential & operator=(const G4StatMFMacroChemicalPotential &right)
G4GLOB_DLL std::ostream G4cerr