G4StatMFMacroChemicalPotential Class Reference

#include <G4StatMFMacroChemicalPotential.hh>

Public Member Functions
G4double	CalcChemicalPotentialNu (void)
	G4StatMFMacroChemicalPotential (const G4double anA, const G4double aZ, const G4double kappa, const G4double temp, std::vector< G4VStatMFMacroCluster * > *ClusterVector)
G4double	GetChemicalPotentialMu (void) const
G4double	GetChemicalPotentialNu (void) const
G4double	GetMeanMultiplicity (void) const
G4double	operator() (const G4double nu)
	~G4StatMFMacroChemicalPotential ()

Private Member Functions
void	CalcChemicalPotentialMu (const G4double nu)
G4double	CalcMeanZ (const G4double nu)
	G4StatMFMacroChemicalPotential ()
	G4StatMFMacroChemicalPotential (const G4StatMFMacroChemicalPotential &)
G4bool	operator!= (const G4StatMFMacroChemicalPotential &right) const
G4StatMFMacroChemicalPotential &	operator= (const G4StatMFMacroChemicalPotential &right)
G4bool	operator== (const G4StatMFMacroChemicalPotential &right) const

Private Attributes
G4double	_ChemPotentialMu
G4double	_ChemPotentialNu
G4double	_Kappa
G4double	_MeanMultiplicity
G4double	_MeanTemperature
std::vector< G4VStatMFMacroCluster * > *	_theClusters
G4double	theA
G4double	theZ

Detailed Description

Definition at line 43 of file G4StatMFMacroChemicalPotential.hh.

Constructor & Destructor Documentation

G4StatMFMacroChemicalPotential() [1/3]

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.

{};

G4StatMFMacroChemicalPotential() [2/3]

G4StatMFMacroChemicalPotential::G4StatMFMacroChemicalPotential ( )

inlineprivate

Definition at line 68 of file G4StatMFMacroChemicalPotential.hh.

{};

G4StatMFMacroChemicalPotential() [3/3]

G4StatMFMacroChemicalPotential::G4StatMFMacroChemicalPotential ( const G4StatMFMacroChemicalPotential &  )

inlineprivate

Definition at line 71 of file G4StatMFMacroChemicalPotential.hh.

{};

Member Function Documentation

CalcChemicalPotentialMu()

void G4StatMFMacroChemicalPotential::CalcChemicalPotentialMu ( const G4double  nu )

private

Definition at line 130 of file G4StatMFMacroChemicalPotential.cc.

{
  G4StatMFMacroMultiplicity * theMultip = new
    G4StatMFMacroMultiplicity(theA,_Kappa,_MeanTemperature,nu,_theClusters);
  
  _ChemPotentialMu = theMultip->CalcChemicalPotentialMu();
  _MeanMultiplicity = theMultip->GetMeanMultiplicity();
  
  delete theMultip;
  
  return; 
}

References _ChemPotentialMu, _Kappa, _MeanMultiplicity, _MeanTemperature, _theClusters, G4StatMFMacroMultiplicity::CalcChemicalPotentialMu(), G4StatMFMacroMultiplicity::GetMeanMultiplicity(), and theA.

Referenced by CalcMeanZ().

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

References _ChemPotentialNu, G4Solver< Function >::Brent(), anonymous_namespace{paraMaker.cc}::CP, G4cout, G4endl, G4StatMFParameters::GetCoulomb(), G4StatMFParameters::GetGamma0(), G4Pow::GetInstance(), G4Solver< Function >::GetRoot(), operator()(), G4Solver< Function >::SetIntervalLimits(), theA, theZ, and G4Pow::Z23().

Referenced by G4StatMFMacroTemperature::CalcChemicalPotentialNu().

CalcMeanZ()

G4double G4StatMFMacroChemicalPotential::CalcMeanZ ( const G4double  nu )

private

Definition at line 108 of file G4StatMFMacroChemicalPotential.cc.

{
  std::vector<G4VStatMFMacroCluster*>::iterator i;
  for (i= _theClusters->begin()+1; i != _theClusters->end(); ++i) 
    { 
      (*i)->CalcZARatio(nu);
    }
  CalcChemicalPotentialMu(nu);
  // This is important, the Z over A ratio for proton and neutron depends on the 
  // chemical potential Mu, while for the first guess for Chemical potential mu 
  // some values of Z over A ratio. This is the reason for that.
  (*_theClusters->begin())->CalcZARatio(nu);
  
  G4double MeanZ = 0.0;
  G4int n = 1;
  for (i = _theClusters->begin(); i != _theClusters->end(); ++i)
    {
      MeanZ += (n++) * (*i)->GetZARatio() * (*i)->GetMeanMultiplicity(); 
    }
  return MeanZ;
}

References _theClusters, CalcChemicalPotentialMu(), and CLHEP::detail::n.

Referenced by operator()().

GetChemicalPotentialMu()

G4double G4StatMFMacroChemicalPotential::GetChemicalPotentialMu ( void  ) const

inline

Definition at line 83 of file G4StatMFMacroChemicalPotential.hh.

{return _ChemPotentialMu;}

References _ChemPotentialMu.

Referenced by G4StatMFMacroTemperature::CalcChemicalPotentialNu().

GetChemicalPotentialNu()

G4double G4StatMFMacroChemicalPotential::GetChemicalPotentialNu ( void  ) const

inline

Definition at line 85 of file G4StatMFMacroChemicalPotential.hh.

{return _ChemPotentialNu;}

References _ChemPotentialNu.

GetMeanMultiplicity()

G4double G4StatMFMacroChemicalPotential::GetMeanMultiplicity ( void  ) const

inline

Definition at line 81 of file G4StatMFMacroChemicalPotential.hh.

{return _MeanMultiplicity;}

References _MeanMultiplicity.

Referenced by G4StatMFMacroTemperature::CalcChemicalPotentialNu().

operator!=()

G4bool G4StatMFMacroChemicalPotential::operator!= ( const G4StatMFMacroChemicalPotential &  right ) const

private

Definition at line 50 of file G4StatMFMacroChemicalPotential.cc.

{
    throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMacroChemicalPotential::operator!= meant to not be accessible");
    return true;
}

operator()()

G4double G4StatMFMacroChemicalPotential::operator() ( const G4double  nu )

inline

Definition at line 63 of file G4StatMFMacroChemicalPotential.hh.

    { return (theZ - this->CalcMeanZ(nu))/theZ; }   

References CalcMeanZ(), and theZ.

Referenced by CalcChemicalPotentialNu().

operator=()

G4StatMFMacroChemicalPotential & G4StatMFMacroChemicalPotential::operator= ( const G4StatMFMacroChemicalPotential &  right )

private

Definition at line 37 of file G4StatMFMacroChemicalPotential.cc.

{
    throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMacroChemicalPotential::operator= meant to not be accessible");
    return *this;
}

operator==()

G4bool G4StatMFMacroChemicalPotential::operator== ( const G4StatMFMacroChemicalPotential &  right ) const

private

Definition at line 43 of file G4StatMFMacroChemicalPotential.cc.

{
    throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMacroChemicalPotential::operator== meant to not be accessible");
    return false;
}

Field Documentation

_ChemPotentialMu

G4double G4StatMFMacroChemicalPotential::_ChemPotentialMu

private

Definition at line 107 of file G4StatMFMacroChemicalPotential.hh.

Referenced by CalcChemicalPotentialMu(), and GetChemicalPotentialMu().

_ChemPotentialNu

G4double G4StatMFMacroChemicalPotential::_ChemPotentialNu

private

Definition at line 109 of file G4StatMFMacroChemicalPotential.hh.

Referenced by CalcChemicalPotentialNu(), and GetChemicalPotentialNu().

_Kappa

G4double G4StatMFMacroChemicalPotential::_Kappa

private

Definition at line 101 of file G4StatMFMacroChemicalPotential.hh.

Referenced by CalcChemicalPotentialMu().

_MeanMultiplicity

G4double G4StatMFMacroChemicalPotential::_MeanMultiplicity

private

Definition at line 103 of file G4StatMFMacroChemicalPotential.hh.

Referenced by CalcChemicalPotentialMu(), and GetMeanMultiplicity().

_MeanTemperature

G4double G4StatMFMacroChemicalPotential::_MeanTemperature

private

Definition at line 105 of file G4StatMFMacroChemicalPotential.hh.

Referenced by CalcChemicalPotentialMu().

_theClusters

std::vector<G4VStatMFMacroCluster*>* G4StatMFMacroChemicalPotential::_theClusters

private

Definition at line 111 of file G4StatMFMacroChemicalPotential.hh.

Referenced by CalcChemicalPotentialMu(), and CalcMeanZ().

theA

G4double G4StatMFMacroChemicalPotential::theA

private

Definition at line 97 of file G4StatMFMacroChemicalPotential.hh.

Referenced by CalcChemicalPotentialMu(), and CalcChemicalPotentialNu().

theZ

G4double G4StatMFMacroChemicalPotential::theZ

private

Definition at line 99 of file G4StatMFMacroChemicalPotential.hh.

Referenced by CalcChemicalPotentialNu(), and operator()().

---

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

• source/processes/hadronic/models/de_excitation/multifragmentation/include/G4StatMFMacroChemicalPotential.hh
• source/processes/hadronic/models/de_excitation/multifragmentation/src/G4StatMFMacroChemicalPotential.cc