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.

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

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(), python.hepunit::elm_coupling, 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