G4StatMFMacroChemicalPotential Class Reference

#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 44 of file G4StatMFMacroChemicalPotential.hh.

Constructor & Destructor Documentation

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 ( )

inline

Definition at line 62 of file G4StatMFMacroChemicalPotential.hh.

{};

Member Function Documentation

G4double G4StatMFMacroChemicalPotential::CalcChemicalPotentialNu

(

void

)

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

Here is the call graph for this function:

G4double G4StatMFMacroChemicalPotential::GetChemicalPotentialMu ( void  ) const

inline

Definition at line 84 of file G4StatMFMacroChemicalPotential.hh.

{return _ChemPotentialMu;}

G4double G4StatMFMacroChemicalPotential::GetChemicalPotentialNu ( void  ) const

inline

Definition at line 86 of file G4StatMFMacroChemicalPotential.hh.

{return _ChemPotentialNu;}

G4double G4StatMFMacroChemicalPotential::GetMeanMultiplicity ( void  ) const

inline

Definition at line 82 of file G4StatMFMacroChemicalPotential.hh.

{return _MeanMultiplicity;}

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

inline

Definition at line 64 of file G4StatMFMacroChemicalPotential.hh.

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

Here is the caller graph for this function:

---

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

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