G4StatMFMacroChemicalPotential Class Reference

#include <G4StatMFMacroChemicalPotential.hh>

Collaboration diagram for G4StatMFMacroChemicalPotential:

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)

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

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

Detailed Description

Definition at line 44 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 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.

{};

G4StatMFMacroChemicalPotential() [2/3]

G4StatMFMacroChemicalPotential::G4StatMFMacroChemicalPotential ( )

inlineprivate

Definition at line 69 of file G4StatMFMacroChemicalPotential.hh.

{};

G4StatMFMacroChemicalPotential() [3/3]

G4StatMFMacroChemicalPotential::G4StatMFMacroChemicalPotential ( const G4StatMFMacroChemicalPotential &  )

inlineprivate

Definition at line 72 of file G4StatMFMacroChemicalPotential.hh.

{};

Here is the call graph for this function:

Member Function Documentation

CalcChemicalPotentialMu()

void G4StatMFMacroChemicalPotential::CalcChemicalPotentialMu ( const G4double  nu )

private

Definition at line 131 of file G4StatMFMacroChemicalPotential.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

CalcChemicalPotentialNu()

G4double G4StatMFMacroChemicalPotential::CalcChemicalPotentialNu

(

void

)

Definition at line 57 of file G4StatMFMacroChemicalPotential.cc.

{
  G4Pow* g4pow = G4Pow::GetInstance();
  G4double CP = G4StatMFParameters::GetCoulomb();

  // Initial value for _ChemPotentialNu 
  _ChemPotentialNu = (theZ/theA)*(8.0*G4StatMFParameters::GetGamma0()
                  +2.0*CP*g4pow->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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

CalcMeanZ()

G4double G4StatMFMacroChemicalPotential::CalcMeanZ ( const G4double  nu )

private

Definition at line 109 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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetChemicalPotentialMu()

G4double G4StatMFMacroChemicalPotential::GetChemicalPotentialMu ( void  ) const

inline

Definition at line 84 of file G4StatMFMacroChemicalPotential.hh.

{return _ChemPotentialMu;}

Here is the caller graph for this function:

GetChemicalPotentialNu()

G4double G4StatMFMacroChemicalPotential::GetChemicalPotentialNu ( void  ) const

inline

Definition at line 86 of file G4StatMFMacroChemicalPotential.hh.

{return _ChemPotentialNu;}

Here is the call graph for this function:

GetMeanMultiplicity()

G4double G4StatMFMacroChemicalPotential::GetMeanMultiplicity ( void  ) const

inline

Definition at line 82 of file G4StatMFMacroChemicalPotential.hh.

{return _MeanMultiplicity;}

Here is the caller graph for this function:

operator!=()

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

private

Definition at line 51 of file G4StatMFMacroChemicalPotential.cc.

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

Here is the caller graph for this function:

operator()()

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

inline

Definition at line 64 of file G4StatMFMacroChemicalPotential.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

operator=()

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

private

Definition at line 38 of file G4StatMFMacroChemicalPotential.cc.

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

Here is the caller graph for this function:

operator==()

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

private

Definition at line 44 of file G4StatMFMacroChemicalPotential.cc.

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

Here is the caller graph for this function:

Member Data Documentation

_ChemPotentialMu

G4double G4StatMFMacroChemicalPotential::_ChemPotentialMu

private

Definition at line 108 of file G4StatMFMacroChemicalPotential.hh.

_ChemPotentialNu

G4double G4StatMFMacroChemicalPotential::_ChemPotentialNu

private

Definition at line 110 of file G4StatMFMacroChemicalPotential.hh.

_Kappa

G4double G4StatMFMacroChemicalPotential::_Kappa

private

Definition at line 102 of file G4StatMFMacroChemicalPotential.hh.

_MeanMultiplicity

G4double G4StatMFMacroChemicalPotential::_MeanMultiplicity

private

Definition at line 104 of file G4StatMFMacroChemicalPotential.hh.

_MeanTemperature

G4double G4StatMFMacroChemicalPotential::_MeanTemperature

private

Definition at line 106 of file G4StatMFMacroChemicalPotential.hh.

_theClusters

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

private

Definition at line 112 of file G4StatMFMacroChemicalPotential.hh.

theA

G4double G4StatMFMacroChemicalPotential::theA

private

Definition at line 98 of file G4StatMFMacroChemicalPotential.hh.

theZ

G4double G4StatMFMacroChemicalPotential::theZ

private

Definition at line 100 of file G4StatMFMacroChemicalPotential.hh.

---

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

• G4StatMFMacroChemicalPotential.hh
• G4StatMFMacroChemicalPotential.cc