#include <G4StatMFMacroChemicalPotential.hh>

Public Member Functions
	G4StatMFMacroChemicalPotential (const G4double anA, const G4double aZ, const G4double kappa, const G4double temp, std::vector< G4VStatMFMacroCluster * > *ClusterVector)

	~G4StatMFMacroChemicalPotential ()

G4double	operator() (const G4double nu)

G4double	GetMeanMultiplicity (void) const

G4double	GetChemicalPotentialMu (void) const

G4double	GetChemicalPotentialNu (void) const

G4double	CalcChemicalPotentialNu (void)

Detailed Description

Definition at line 44 of file G4StatMFMacroChemicalPotential.hh.

Constructor & Destructor Documentation

G4StatMFMacroChemicalPotential::G4StatMFMacroChemicalPotential	(	const G4double	anA,
		const G4double	aZ,
		const G4double	kappa,
		const G4double	temp,
		std::vector< G4VStatMFMacroCluster * > *	ClusterVector
	)

inline

Definition at line 48 of file G4StatMFMacroChemicalPotential.hh.

                                                                       :
     theA(anA),
     theZ(aZ),
     _Kappa(kappa),
     _MeanMultiplicity(0.0),
     _MeanTemperature(temp),
     _ChemPotentialMu(0.0),
     _ChemPotentialNu(0.0),
     _theClusters(ClusterVector) 
     {};

G4StatMFMacroChemicalPotential::~G4StatMFMacroChemicalPotential ( )

inline

Definition at line 62 of file G4StatMFMacroChemicalPotential.hh.

62 {};

Member Function Documentation

G4double G4StatMFMacroChemicalPotential::CalcChemicalPotentialNu ( void )

Definition at line 59 of file G4StatMFMacroChemicalPotential.cc.

References G4Solver< Function >::Brent(), python.hepunit::elm_coupling, G4cerr, G4endl, G4StatMFParameters::GetGamma0(), G4StatMFParameters::GetKappaCoulomb(), G4StatMFParameters::Getr0(), G4Solver< Function >::GetRoot(), operator()(), and G4Solver< Function >::SetIntervalLimits().

 {
     G4double CP = ((3./5.)*elm_coupling/G4StatMFParameters::Getr0())*
     (1.0-1.0/std::pow(1.0+G4StatMFParameters::GetKappaCoulomb(),1.0/3.0));
 
     // Initial value for _ChemPotentialNu   
     _ChemPotentialNu = (theZ/theA)*(8.0*G4StatMFParameters::GetGamma0()+2.0*CP*std::pow(theA,2./3.)) -
     4.0*G4StatMFParameters::GetGamma0();
         
 
     G4double ChemPa = _ChemPotentialNu;
     G4double ChemPb = 0.5*_ChemPotentialNu;
     
     G4double fChemPa = this->operator()(ChemPa); 
     G4double fChemPb = this->operator()(ChemPb); 
 
     if (fChemPa*fChemPb > 0.0) {    
     // bracketing the solution
     if (fChemPa < 0.0) {
         do {
         ChemPb -= 1.5*std::abs(ChemPb-ChemPa);
         fChemPb = this->operator()(ChemPb);   
         } while (fChemPb < 0.0);
     } else {
         do {
         ChemPb += 1.5*std::abs(ChemPb-ChemPa);
         fChemPb = this->operator()(ChemPb);
         } while (fChemPb > 0.0);
     }
     }
 
     G4Solver<G4StatMFMacroChemicalPotential> * theSolver =
       new G4Solver<G4StatMFMacroChemicalPotential>(100,1.e-4);
     theSolver->SetIntervalLimits(ChemPa,ChemPb);
     //    if (!theSolver->Crenshaw(*this)) 
     if (!theSolver->Brent(*this)){
       G4cerr <<"G4StatMFMacroChemicalPotential:"<<" ChemPa="<<ChemPa<<" ChemPb="<<ChemPb<< G4endl;
       G4cerr <<"G4StatMFMacroChemicalPotential:"<<" fChemPa="<<fChemPa<<" fChemPb="<<fChemPb<< G4endl;
       throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMacroChemicalPotential::CalcChemicalPotentialNu: I couldn't find the root.");
     }
     _ChemPotentialNu = theSolver->GetRoot();
     delete theSolver;
     return _ChemPotentialNu;
 }