#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::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::~G4StatMFMacroChemicalPotential ( )

inline

Definition at line 62 of file G4StatMFMacroChemicalPotential.hh.

62{};

Member Function Documentation

◆ CalcChemicalPotentialNu()

G4double G4StatMFMacroChemicalPotential::CalcChemicalPotentialNu ( void )

Definition at line 59 of file G4StatMFMacroChemicalPotential.cc.

{
    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;
}