45 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMacroMultiplicity::operator= meant to not be accessible");
51 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMacroMultiplicity::operator== meant to not be accessible");
58 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMacroMultiplicity::operator!= meant to not be accessible");
78 (5.0/3.0)*CP*ZA5*ZA5*g4calc->
Z23(5) -
90 const G4double intervalWidth = 1.e-4;
95 while (fChemPa*fChemPb > 0.0 && iterations < 100)
100 ChemPa += 0.6*(ChemPa-ChemPb);
105 ChemPb += 0.6*(ChemPb-ChemPa);
110 if (fChemPa*fChemPb > 0.0)
112 G4cout <<
"G4StatMFMacroMultiplicity:"<<
" ChemPa="<<ChemPa
113 <<
" ChemPb="<<ChemPb<<
G4endl;
114 G4cout <<
"G4StatMFMacroMultiplicity:"<<
" fChemPa="<<fChemPa
115 <<
" fChemPb="<<fChemPb<<
G4endl;
116 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMacroMultiplicity::CalcChemicalPotentialMu: I couldn't bracket the root.");
118 else if (fChemPa*fChemPb < 0.0 &&
std::abs(ChemPa-ChemPb) > intervalWidth)
124 if (!theSolver->
Brent(*
this))
126 G4cout <<
"G4StatMFMacroMultiplicity:"<<
" ChemPa="<<ChemPa
127 <<
" ChemPb="<<ChemPb<<
G4endl;
128 throw G4HadronicException(__FILE__, __LINE__,
"G4StatMFMacroMultiplicity::CalcChemicalPotentialMu: I couldn't find the root.");
151 for (std::vector<G4VStatMFMacroCluster*>::iterator i =
_theClusters->begin();
156 MeanA += multip*(n++);