#include <G4StatMFMacroTemperature.hh>

Public Member Functions
G4double	CalcTemperature (void)

	G4StatMFMacroTemperature (const G4double anA, const G4double aZ, const G4double ExEnergy, const G4double FreeE0, const G4double kappa, std::vector< G4VStatMFMacroCluster * > *ClusterVector)

G4double	GetChemicalPotentialMu (void) const

G4double	GetChemicalPotentialNu (void) const

G4double	GetEntropy (void) const

G4double	GetMeanMultiplicity (void) const

G4double	GetTemperature (void) const

G4double	operator() (const G4double T)

	~G4StatMFMacroTemperature ()

Private Member Functions
void	CalcChemicalPotentialNu (const G4double T)

G4double	FragsExcitEnergy (const G4double T)

	G4StatMFMacroTemperature ()

	G4StatMFMacroTemperature (const G4StatMFMacroTemperature &)

G4bool	operator!= (const G4StatMFMacroTemperature &right) const

G4StatMFMacroTemperature &	operator= (const G4StatMFMacroTemperature &right)

G4bool	operator== (const G4StatMFMacroTemperature &right) const

Private Attributes
G4double	_ChemPotentialMu

G4double	_ChemPotentialNu

G4double	_ExEnergy

G4double	_FreeInternalE0

G4double	_Kappa

G4double	_MeanEntropy

G4double	_MeanMultiplicity

G4double	_MeanTemperature

std::vector< G4VStatMFMacroCluster * > *	_theClusters

G4double	theA

G4double	theZ

Detailed Description

Definition at line 41 of file G4StatMFMacroTemperature.hh.

Constructor & Destructor Documentation

◆ G4StatMFMacroTemperature() [1/3]

G4StatMFMacroTemperature::G4StatMFMacroTemperature	(	const G4double	anA,
		const G4double	aZ,
		const G4double	ExEnergy,
		const G4double	FreeE0,
		const G4double	kappa,
		std::vector< G4VStatMFMacroCluster * > *	ClusterVector
	)

Definition at line 46 of file G4StatMFMacroTemperature.cc.

                                                     :
  theA(anA),
  theZ(aZ),
  _ExEnergy(ExEnergy),
  _FreeInternalE0(FreeE0),
  _Kappa(kappa),
  _MeanMultiplicity(0.0),
  _MeanTemperature(0.0),
  _ChemPotentialMu(0.0),
  _ChemPotentialNu(0.0),
  _MeanEntropy(0.0),
  _theClusters(ClusterVector) 
{}

◆ ~G4StatMFMacroTemperature()

G4StatMFMacroTemperature::~G4StatMFMacroTemperature ( )

Definition at line 62 of file G4StatMFMacroTemperature.cc.

63{}

◆ G4StatMFMacroTemperature() [2/3]

G4StatMFMacroTemperature::G4StatMFMacroTemperature ( )

private

◆ G4StatMFMacroTemperature() [3/3]

G4StatMFMacroTemperature::G4StatMFMacroTemperature ( const G4StatMFMacroTemperature & )

inlineprivate

Definition at line 61 of file G4StatMFMacroTemperature.hh.

61{};

Member Function Documentation

◆ CalcChemicalPotentialNu()

void G4StatMFMacroTemperature::CalcChemicalPotentialNu ( const G4double T )

private

Definition at line 185 of file G4StatMFMacroTemperature.cc.

{
  G4StatMFMacroChemicalPotential * theChemPot = new
    G4StatMFMacroChemicalPotential(theA,theZ,_Kappa,T,_theClusters);
 
  _ChemPotentialNu = theChemPot->CalcChemicalPotentialNu();
  _ChemPotentialMu = theChemPot->GetChemicalPotentialMu();
  _MeanMultiplicity = theChemPot->GetMeanMultiplicity();        
  delete theChemPot;
            
  return;
}

References _ChemPotentialMu, _ChemPotentialNu, _Kappa, _MeanMultiplicity, _theClusters, G4StatMFMacroChemicalPotential::CalcChemicalPotentialNu(), G4StatMFMacroChemicalPotential::GetChemicalPotentialMu(), G4StatMFMacroChemicalPotential::GetMeanMultiplicity(), theA, and theZ.

Referenced by FragsExcitEnergy().

◆ CalcTemperature()

G4double G4StatMFMacroTemperature::CalcTemperature ( void )

Definition at line 65 of file G4StatMFMacroTemperature.cc.

{
  // Inital guess for the interval of the ensemble temperature values
  G4double Ta = 0.5; 
  G4double Tb = std::max(std::sqrt(_ExEnergy/(theA*0.12)),0.01*MeV);
    
  G4double fTa = this->operator()(Ta); 
  G4double fTb = this->operator()(Tb); 
 
  // Bracketing the solution
  // T should be greater than 0.
  // The interval is [Ta,Tb]
  // We start with a value for Ta = 0.5 MeV
  // it should be enough to have fTa > 0 If it isn't 
  // the case, we decrease Ta. But carefully, because 
  // fTa growes very fast when Ta is near 0 and we could have
  // an overflow.
 
  G4int iterations = 0;  
  // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
  while (fTa < 0.0 && ++iterations < 10) {
    Ta -= 0.5*Ta;
    fTa = this->operator()(Ta);
  }
  // Usually, fTb will be less than 0, but if it is not the case: 
  iterations = 0;  
  // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
  while (fTa*fTb > 0.0 && iterations++ < 10) {
    Tb += 2.*std::fabs(Tb-Ta);
    fTb = this->operator()(Tb);
  }
    
  if (fTa*fTb > 0.0) {
    G4cerr <<"G4StatMFMacroTemperature:"<<" Ta="<<Ta<<" Tb="<<Tb<< G4endl;
    G4cerr <<"G4StatMFMacroTemperature:"<<" fTa="<<fTa<<" fTb="<<fTb<< G4endl;
    throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMacroTemperature::CalcTemperature: I couldn't bracket the solution.");
  }
 
  G4Solver<G4StatMFMacroTemperature> * theSolver = 
    new G4Solver<G4StatMFMacroTemperature>(100,1.e-4);
  theSolver->SetIntervalLimits(Ta,Tb);
  if (!theSolver->Crenshaw(*this)){ 
    G4cout <<"G4StatMFMacroTemperature, Crenshaw method failed:"<<" Ta="
       <<Ta<<" Tb="<<Tb<< G4endl;
    G4cout <<"G4StatMFMacroTemperature, Crenshaw method failed:"<<" fTa="
       <<fTa<<" fTb="<<fTb<< G4endl;
  }
  _MeanTemperature = theSolver->GetRoot();
  G4double FunctionValureAtRoot =  this->operator()(_MeanTemperature);
  delete  theSolver;
 
  // Verify if the root is found and it is indeed within the physical domain, 
  // say, between 1 and 50 MeV, otherwise try Brent method:
  if (std::fabs(FunctionValureAtRoot) > 5.e-2) {
    if (_MeanTemperature < 1. || _MeanTemperature > 50.) {
      G4cout << "Crenshaw method failed; function = " << FunctionValureAtRoot 
         << " solution? = " << _MeanTemperature << " MeV " << G4endl;
      G4Solver<G4StatMFMacroTemperature> * theSolverBrent = 
    new G4Solver<G4StatMFMacroTemperature>(200,1.e-3);
      theSolverBrent->SetIntervalLimits(Ta,Tb);
      if (!theSolverBrent->Brent(*this)){
    G4cout <<"G4StatMFMacroTemperature, Brent method failed:"
           <<" Ta="<<Ta<<" Tb="<<Tb<< G4endl;
    G4cout <<"G4StatMFMacroTemperature, Brent method failed:"
           <<" fTa="<<fTa<<" fTb="<<fTb<< G4endl; 
    throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMacroTemperature::CalcTemperature: I couldn't find the root with any method.");
      }
 
      _MeanTemperature = theSolverBrent->GetRoot();
      FunctionValureAtRoot =  this->operator()(_MeanTemperature);
      delete theSolverBrent;
    }
    if (std::abs(FunctionValureAtRoot) > 5.e-2) {
      G4cout << "Brent method failed; function = " << FunctionValureAtRoot 
         << " solution? = " << _MeanTemperature << " MeV " << G4endl;
      throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMacroTemperature::CalcTemperature: I couldn't find the root with any method.");
    }
  }
  //G4cout << "G4StatMFMacroTemperature::CalcTemperature: function = " 
  //<< FunctionValureAtRoot 
  //       << " T(MeV)= " << _MeanTemperature << G4endl;
  return _MeanTemperature;
}

References _ExEnergy, _MeanTemperature, G4Solver< Function >::Brent(), G4Solver< Function >::Crenshaw(), G4cerr, G4cout, G4endl, G4Solver< Function >::GetRoot(), G4INCL::Math::max(), MeV, operator()(), G4Solver< Function >::SetIntervalLimits(), and theA.

Referenced by G4StatMFMacroCanonical::CalculateTemperature().

◆ FragsExcitEnergy()

G4double G4StatMFMacroTemperature::FragsExcitEnergy ( const G4double T )

private

Definition at line 149 of file G4StatMFMacroTemperature.cc.

{
  // Model Parameters
  G4Pow* g4calc = G4Pow::GetInstance();
  G4double R0 = G4StatMFParameters::Getr0()*g4calc->Z13(theA);
  G4double R = R0*g4calc->A13(1.0+G4StatMFParameters::GetKappaCoulomb());
  G4double FreeVol = _Kappa*(4.*pi/3.)*R0*R0*R0; 
 
  // Calculate Chemical potentials
  CalcChemicalPotentialNu(T);
 
 
  // Average total fragment energy
  G4double AverageEnergy = 0.0;
  std::vector<G4VStatMFMacroCluster*>::iterator i;
  for (i =  _theClusters->begin(); i != _theClusters->end(); ++i) 
    {
      AverageEnergy += (*i)->GetMeanMultiplicity() * (*i)->CalcEnergy(T);
    }
    
  // Add Coulomb energy         
  AverageEnergy += 0.6*elm_coupling*theZ*theZ/R;        
    
  // Calculate mean entropy
  _MeanEntropy = 0.0;
  for (i = _theClusters->begin(); i != _theClusters->end(); ++i) 
    {
      _MeanEntropy += (*i)->CalcEntropy(T,FreeVol); 
    }
 
  // Excitation energy per nucleon
  return AverageEnergy - _FreeInternalE0;
}

References _FreeInternalE0, _Kappa, _MeanEntropy, _theClusters, G4Pow::A13(), CalcChemicalPotentialNu(), source.hepunit::elm_coupling, G4Pow::GetInstance(), G4StatMFParameters::GetKappaCoulomb(), G4StatMFParameters::Getr0(), pi, theA, theZ, and G4Pow::Z13().

Referenced by operator()().

◆ GetChemicalPotentialMu()

G4double G4StatMFMacroTemperature::GetChemicalPotentialMu ( void ) const

inline

Definition at line 73 of file G4StatMFMacroTemperature.hh.

73{return _ChemPotentialMu;}

References _ChemPotentialMu.

Referenced by G4StatMFMacroCanonical::CalculateTemperature().

◆ GetChemicalPotentialNu()

G4double G4StatMFMacroTemperature::GetChemicalPotentialNu ( void ) const

inline

Definition at line 75 of file G4StatMFMacroTemperature.hh.

75{return _ChemPotentialNu;}

References _ChemPotentialNu.

Referenced by G4StatMFMacroCanonical::CalculateTemperature().

◆ GetEntropy()

G4double G4StatMFMacroTemperature::GetEntropy ( void ) const

inline

Definition at line 79 of file G4StatMFMacroTemperature.hh.

79{return _MeanEntropy;}

References _MeanEntropy.

Referenced by G4StatMFMacroCanonical::CalculateTemperature().

◆ GetMeanMultiplicity()

G4double G4StatMFMacroTemperature::GetMeanMultiplicity ( void ) const

inline

Definition at line 71 of file G4StatMFMacroTemperature.hh.

71{return _MeanMultiplicity;}

References _MeanMultiplicity.

Referenced by G4StatMFMacroCanonical::CalculateTemperature().

◆ GetTemperature()

G4double G4StatMFMacroTemperature::GetTemperature ( void ) const

inline

Definition at line 77 of file G4StatMFMacroTemperature.hh.

77{return _MeanTemperature;}

References _MeanTemperature.

◆ operator!=()

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

private

◆ operator()()

G4double G4StatMFMacroTemperature::operator() ( const G4double T )

inline

Definition at line 52 of file G4StatMFMacroTemperature.hh.

53 { return (_ExEnergy - this->FragsExcitEnergy(T))/_ExEnergy; }

G4StatMFMacroTemperature::FragsExcitEnergy

G4double FragsExcitEnergy(const G4double T)

Definition: G4StatMFMacroTemperature.cc:149

References _ExEnergy, and FragsExcitEnergy().

Referenced by CalcTemperature().

◆ operator=()

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

private

◆ operator==()

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

private

Field Documentation

◆ _ChemPotentialMu

G4double G4StatMFMacroTemperature::_ChemPotentialMu

private

Definition at line 105 of file G4StatMFMacroTemperature.hh.

Referenced by CalcChemicalPotentialNu(), and GetChemicalPotentialMu().

◆ _ChemPotentialNu

G4double G4StatMFMacroTemperature::_ChemPotentialNu

private

Definition at line 107 of file G4StatMFMacroTemperature.hh.

Referenced by CalcChemicalPotentialNu(), and GetChemicalPotentialNu().

◆ _ExEnergy

G4double G4StatMFMacroTemperature::_ExEnergy

private

Definition at line 95 of file G4StatMFMacroTemperature.hh.

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

◆ _FreeInternalE0

G4double G4StatMFMacroTemperature::_FreeInternalE0

private

Definition at line 97 of file G4StatMFMacroTemperature.hh.

Referenced by FragsExcitEnergy().

◆ _Kappa

G4double G4StatMFMacroTemperature::_Kappa

private

Definition at line 99 of file G4StatMFMacroTemperature.hh.

Referenced by CalcChemicalPotentialNu(), and FragsExcitEnergy().

◆ _MeanEntropy

G4double G4StatMFMacroTemperature::_MeanEntropy

private

Definition at line 109 of file G4StatMFMacroTemperature.hh.

Referenced by FragsExcitEnergy(), and GetEntropy().

◆ _MeanMultiplicity

G4double G4StatMFMacroTemperature::_MeanMultiplicity

private

Definition at line 101 of file G4StatMFMacroTemperature.hh.

Referenced by CalcChemicalPotentialNu(), and GetMeanMultiplicity().

◆ _MeanTemperature

G4double G4StatMFMacroTemperature::_MeanTemperature

private

Definition at line 103 of file G4StatMFMacroTemperature.hh.

Referenced by CalcTemperature(), and GetTemperature().

◆ _theClusters

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

private

Definition at line 111 of file G4StatMFMacroTemperature.hh.

Referenced by CalcChemicalPotentialNu(), and FragsExcitEnergy().

◆ theA

G4double G4StatMFMacroTemperature::theA

private

Definition at line 91 of file G4StatMFMacroTemperature.hh.

Referenced by CalcChemicalPotentialNu(), CalcTemperature(), and FragsExcitEnergy().

◆ theZ

G4double G4StatMFMacroTemperature::theZ

private

Definition at line 93 of file G4StatMFMacroTemperature.hh.

Referenced by CalcChemicalPotentialNu(), and FragsExcitEnergy().

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

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

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4StatMFMacroTemperature() [1/3]

◆ ~G4StatMFMacroTemperature()

◆ G4StatMFMacroTemperature() [2/3]

◆ G4StatMFMacroTemperature() [3/3]

Member Function Documentation

◆ CalcChemicalPotentialNu()

◆ CalcTemperature()

◆ FragsExcitEnergy()

◆ GetChemicalPotentialMu()

◆ GetChemicalPotentialNu()

◆ GetEntropy()

◆ GetMeanMultiplicity()

◆ GetTemperature()

◆ operator!=()

◆ operator()()

◆ operator=()

◆ operator==()

Field Documentation

◆ _ChemPotentialMu

◆ _ChemPotentialNu

◆ _ExEnergy

◆ _FreeInternalE0

◆ _Kappa

◆ _MeanEntropy

◆ _MeanMultiplicity

◆ _MeanTemperature

◆ _theClusters

◆ theA

◆ theZ