#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 43 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 47 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 61 of file G4StatMFMacroChemicalPotential.hh.

61{};

Member Function Documentation

◆ CalcChemicalPotentialNu()

G4double G4StatMFMacroChemicalPotential::CalcChemicalPotentialNu ( void )

Definition at line 56 of file G4StatMFMacroChemicalPotential.cc.

{
  G4Pow* g4calc = G4Pow::GetInstance();
  G4double CP = G4StatMFParameters::GetCoulomb();
 
  // Initial value for _ChemPotentialNu 
  _ChemPotentialNu = (theZ/theA)*(8.0*G4StatMFParameters::GetGamma0()
                  +2.0*CP*g4calc->Z23(theA))
    - 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);   
    // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
      } while (fChemPb < 0.0);
    } else {
      do {
    ChemPb += 1.5*std::abs(ChemPb-ChemPa);
    fChemPb = this->operator()(ChemPb);
    // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
      } 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)){
    G4cout <<"G4StatMFMacroChemicalPotential:"<<" ChemPa="<<ChemPa
       <<" ChemPb="<<ChemPb<< G4endl;
    G4cout <<"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;
}