G4StatMFMacroMultiplicity Class Reference

#include <G4StatMFMacroMultiplicity.hh>

Public Member Functions
	G4StatMFMacroMultiplicity (const G4double anA, const G4double kappa, const G4double temp, const G4double nu, std::vector< G4VStatMFMacroCluster * > *ClusterVector)
	~G4StatMFMacroMultiplicity ()
G4double	operator() (const G4double mu)
G4double	GetMeanMultiplicity (void) const
G4double	GetChemicalPotentialMu (void) const
G4double	CalcChemicalPotentialMu (void)

---

Detailed Description

Definition at line 43 of file G4StatMFMacroMultiplicity.hh.

---

Constructor & Destructor Documentation

G4StatMFMacroMultiplicity::G4StatMFMacroMultiplicity	(	const G4double	anA,
		const G4double	kappa,
		const G4double	temp,
		const G4double	nu,
		std::vector< G4VStatMFMacroCluster * > *	ClusterVector
	)			[inline]

Definition at line 47 of file G4StatMFMacroMultiplicity.hh.

00051                                                                                  :
00052         theA(anA),
00053         _Kappa(kappa),
00054         _MeanMultiplicity(0.0),
00055         _MeanTemperature(temp),
00056         _ChemPotentialMu(0.0),
00057         _ChemPotentialNu(nu),
00058         _theClusters(ClusterVector) 
00059         {};
        

G4StatMFMacroMultiplicity::~G4StatMFMacroMultiplicity

(

)

[inline]

Definition at line 61 of file G4StatMFMacroMultiplicity.hh.

{};

---

Member Function Documentation

G4double G4StatMFMacroMultiplicity::CalcChemicalPotentialMu

(

void

)

Definition at line 65 of file G4StatMFMacroMultiplicity.cc.

References G4StatMFParameters::Beta(), G4Solver< Function >::Brent(), G4cerr, G4endl, G4StatMFParameters::GetE0(), 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));

    // starting value for chemical potential \mu
    // it is the derivative of F(T,V)-\nu*Z w.r.t. Af in Af=5
    G4double ZA5 = _theClusters->operator[](4)->GetZARatio();
    G4double ILD5 = _theClusters->operator[](4)->GetInvLevelDensity();
    _ChemPotentialMu = -G4StatMFParameters::GetE0()-
        _MeanTemperature*_MeanTemperature/ILD5 -
        _ChemPotentialNu*ZA5 + 
        G4StatMFParameters::GetGamma0()*(1.0-2.0*ZA5)*(1.0-2.0*ZA5) +
        (2.0/3.0)*G4StatMFParameters::Beta(_MeanTemperature)/std::pow(5.,1./3.) +
        (5.0/3.0)*CP*ZA5*ZA5*std::pow(5.,2./3.) -
        1.5*_MeanTemperature/5.0;
                


    G4double ChemPa = _ChemPotentialMu;
    if (ChemPa/_MeanTemperature > 10.0) ChemPa = 10.0*_MeanTemperature;
    G4double ChemPb = ChemPa - 0.5*std::abs(ChemPa);
    
    
    G4double fChemPa = this->operator()(ChemPa); 
    G4double fChemPb = this->operator()(ChemPb); 
    

    // Set the precision level for locating the root. 
    // If the root is inside this interval, then it's done! 
    G4double intervalWidth = 1.e-4;

    // bracketing the solution
    G4int iterations = 0;
    while (fChemPa*fChemPb > 0.0 && iterations < 100) 
    {
        if (std::abs(fChemPa) <= std::abs(fChemPb)) 
        {
            ChemPa += 0.6*(ChemPa-ChemPb);
            fChemPa = this->operator()(ChemPa);
            iterations++;
        } 
        else 
        {
            ChemPb += 0.6*(ChemPb-ChemPa);
            fChemPb = this->operator()(ChemPb);
            iterations++;
        }
    }

    if (fChemPa*fChemPb > 0.0) // the bracketing failed, complain 
    {
      G4cerr <<"G4StatMFMacroMultiplicity:"<<" ChemPa="<<ChemPa<<" ChemPb="<<ChemPb<< G4endl;
      G4cerr <<"G4StatMFMacroMultiplicity:"<<" fChemPa="<<fChemPa<<" fChemPb="<<fChemPb<< G4endl;
        throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMacroMultiplicity::CalcChemicalPotentialMu: I couldn't bracket the root.");
    }
    else if (fChemPa*fChemPb < 0.0 && std::abs(ChemPa-ChemPb) > intervalWidth) // the bracketing was OK, try to locate the root
    {   
    G4Solver<G4StatMFMacroMultiplicity> * theSolver = new G4Solver<G4StatMFMacroMultiplicity>(100,intervalWidth);
    theSolver->SetIntervalLimits(ChemPa,ChemPb);
    //    if (!theSolver->Crenshaw(*this)) 
    if (!theSolver->Brent(*this)) 
    {
      G4cerr <<"G4StatMFMacroMultiplicity:"<<" ChemPa="<<ChemPa<<" ChemPb="<<ChemPb<< G4endl;
      G4cerr <<"G4StatMFMacroMultiplicity:"<<" fChemPa="<<fChemPa<<" fChemPb="<<fChemPb<< G4endl;
        throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMacroMultiplicity::CalcChemicalPotentialMu: I couldn't find the root.");
    }
    _ChemPotentialMu = theSolver->GetRoot();
    delete theSolver;
    }
    else // the root is within the interval, which is shorter then the precision level - all done 
    {
     _ChemPotentialMu = ChemPa;
    }

    return _ChemPotentialMu;
}

G4double G4StatMFMacroMultiplicity::GetChemicalPotentialMu

(

void

)

const [inline]

Definition at line 83 of file G4StatMFMacroMultiplicity.hh.

{return _ChemPotentialMu;}

G4double G4StatMFMacroMultiplicity::GetMeanMultiplicity

(

void

)

const [inline]

Definition at line 81 of file G4StatMFMacroMultiplicity.hh.

{return _MeanMultiplicity;}

G4double G4StatMFMacroMultiplicity::operator()

(

const G4double

mu

)

[inline]

Definition at line 63 of file G4StatMFMacroMultiplicity.hh.

Referenced by CalcChemicalPotentialMu().

        { return (theA - this->CalcMeanA(mu))/theA; }   

---

The documentation for this class was generated from the following files:

• G4StatMFMacroMultiplicity.hh
• G4StatMFMacroMultiplicity.cc

---