G4StatMFMicroCanonical Class Reference

#include <G4StatMFMicroCanonical.hh>

Inheritance diagram for G4StatMFMicroCanonical:

Data Structures
struct	DeleteFragment

Public Types
enum	{ MaxAllowedMultiplicity = 4 }

Public Member Functions
G4StatMFChannel *	ChooseAandZ (const G4Fragment &theFragment)
	G4StatMFMicroCanonical (const G4Fragment &theFragment)
G4double	GetMeanMultiplicity (void) const
G4double	GetMeanTemperature (void) const
	~G4StatMFMicroCanonical ()

Protected Attributes
G4double	__FreeInternalE0
G4double	__MeanEntropy
G4double	__MeanMultiplicity
G4double	__MeanTemperature

Private Member Functions
G4double	CalcEntropyOfCompoundNucleus (const G4Fragment &theFragment, G4double &TConf)
G4double	CalcFreeInternalEnergy (const G4Fragment &theFragment, G4double T)
G4double	CalcInvLevelDensity (G4int anA)
	G4StatMFMicroCanonical ()
	G4StatMFMicroCanonical (const G4StatMFMicroCanonical &right)
void	Initialize (const G4Fragment &theFragment)
G4bool	operator!= (const G4StatMFMicroCanonical &right) const
G4StatMFMicroCanonical &	operator= (const G4StatMFMicroCanonical &right)
G4bool	operator== (const G4StatMFMicroCanonical &right) const

Private Attributes
std::vector< G4StatMFMicroManager * >	_ThePartitionManagerVector
G4double	_WCompoundNucleus

Detailed Description

Definition at line 46 of file G4StatMFMicroCanonical.hh.

Member Enumeration Documentation

anonymous enum

anonymous enum

Enumerator
MaxAllowedMultiplicity

Definition at line 76 of file G4StatMFMicroCanonical.hh.

{MaxAllowedMultiplicity = 4};

Constructor & Destructor Documentation

G4StatMFMicroCanonical() [1/3]

G4StatMFMicroCanonical::G4StatMFMicroCanonical

(

const G4Fragment &

theFragment

)

Definition at line 40 of file G4StatMFMicroCanonical.cc.

{
  // Perform class initialization
  Initialize(theFragment);
}

References Initialize().

~G4StatMFMicroCanonical()

G4StatMFMicroCanonical::~G4StatMFMicroCanonical

(

)

Definition at line 47 of file G4StatMFMicroCanonical.cc.

{
  // garbage collection
  if (!_ThePartitionManagerVector.empty()) {
    std::for_each(_ThePartitionManagerVector.begin(),
            _ThePartitionManagerVector.end(),
            DeleteFragment());
  }
}

References _ThePartitionManagerVector.

G4StatMFMicroCanonical() [2/3]

G4StatMFMicroCanonical::G4StatMFMicroCanonical ( )

inlineprivate

Definition at line 58 of file G4StatMFMicroCanonical.hh.

{};

G4StatMFMicroCanonical() [3/3]

G4StatMFMicroCanonical::G4StatMFMicroCanonical ( const G4StatMFMicroCanonical &  right )

private

Member Function Documentation

CalcEntropyOfCompoundNucleus()

G4double G4StatMFMicroCanonical::CalcEntropyOfCompoundNucleus ( const G4Fragment &  theFragment, G4double &  TConf  )

private

Definition at line 170 of file G4StatMFMicroCanonical.cc.

{
  G4int A = theFragment.GetA_asInt();
  G4double U = theFragment.GetExcitationEnergy();
  G4double A13 = G4Pow::GetInstance()->Z13(A);
  
  G4double Ta = std::max(std::sqrt(U/(0.125*A)),0.0012*MeV); 
  G4double Tb = Ta;
  
  G4double ECompoundNucleus = CalcFreeInternalEnergy(theFragment,Ta);
  G4double Da = (U+__FreeInternalE0-ECompoundNucleus)/U;
  G4double Db = 0.0;
    
  G4double InvLevelDensity = CalcInvLevelDensity(A);
  
  // bracketing the solution
  if (Da == 0.0) {
    TConf = Ta;
    return 2*Ta*A/InvLevelDensity - G4StatMFParameters::DBetaDT(Ta)*A13*A13;
  } else if (Da < 0.0) {
    do {
      Tb -= 0.5*Tb;
      ECompoundNucleus = CalcFreeInternalEnergy(theFragment,Tb);
      Db = (U+__FreeInternalE0-ECompoundNucleus)/U;
    } while (Db < 0.0);
  } else {
    do {
      Tb += 0.5*Tb;
      ECompoundNucleus = CalcFreeInternalEnergy(theFragment,Tb);
      Db = (U+__FreeInternalE0-ECompoundNucleus)/U;
    } while (Db > 0.0);
  }
  
  G4double eps = 1.0e-14 * std::abs(Tb-Ta);
  
  for (G4int i = 0; i < 1000; i++) {
    G4double Tc = (Ta+Tb)*0.5;
    if (std::abs(Ta-Tb) <= eps) {
      TConf = Tc;
      return 2*Tc*A/InvLevelDensity - G4StatMFParameters::DBetaDT(Tc)*A13*A13;
    }
    ECompoundNucleus = CalcFreeInternalEnergy(theFragment,Tc);
    G4double Dc = (U+__FreeInternalE0-ECompoundNucleus)/U;
    
    if (Dc == 0.0) {
      TConf = Tc;
      return 2*Tc*A/InvLevelDensity - G4StatMFParameters::DBetaDT(Tc)*A13*A13;
    }
    
    if (Da*Dc < 0.0) {
      Tb = Tc;
      Db = Dc;
    } else {
      Ta = Tc;
      Da = Dc;
    } 
  }
  
  G4cout << 
    "G4StatMFMicrocanoncal::CalcEntropyOfCompoundNucleus: I can't calculate the temperature" 
     << G4endl;
  
  return 0.0;
}

References G4VStatMFEnsemble::__FreeInternalE0, A, A13, CalcFreeInternalEnergy(), CalcInvLevelDensity(), G4StatMFParameters::DBetaDT(), eps, G4cout, G4endl, G4Fragment::GetA_asInt(), G4Fragment::GetExcitationEnergy(), G4Pow::GetInstance(), G4INCL::Math::max(), MeV, and G4Pow::Z13().

Referenced by Initialize().

CalcFreeInternalEnergy()

G4double G4StatMFMicroCanonical::CalcFreeInternalEnergy ( const G4Fragment &  theFragment, G4double  T  )

private

Definition at line 146 of file G4StatMFMicroCanonical.cc.

{
  G4int A = theFragment.GetA_asInt();
  G4int Z = theFragment.GetZ_asInt();
  G4double A13 = G4Pow::GetInstance()->Z13(A);
  
  G4double InvLevelDensityPar = G4StatMFParameters::GetEpsilon0()
    *(1.0 + 3.0/G4double(A-1));
  
  G4double VolumeTerm = (-G4StatMFParameters::GetE0()+T*T/InvLevelDensityPar)*A;
  
  G4double SymmetryTerm = G4StatMFParameters::GetGamma0()
    *(A - 2*Z)*(A - 2*Z)/G4double(A);
  
  G4double SurfaceTerm = (G4StatMFParameters::Beta(T)
              - T*G4StatMFParameters::DBetaDT(T))*A13*A13;
  
  G4double CoulombTerm = elm_coupling*0.6*Z*Z/(G4StatMFParameters::Getr0()*A13);
  
  return VolumeTerm + SymmetryTerm + SurfaceTerm + CoulombTerm;
}

References A, A13, G4StatMFParameters::Beta(), G4StatMFParameters::DBetaDT(), source.hepunit::elm_coupling, G4Fragment::GetA_asInt(), G4StatMFParameters::GetE0(), G4StatMFParameters::GetEpsilon0(), G4StatMFParameters::GetGamma0(), G4Pow::GetInstance(), G4StatMFParameters::Getr0(), G4Fragment::GetZ_asInt(), Z, and G4Pow::Z13().

Referenced by CalcEntropyOfCompoundNucleus().

CalcInvLevelDensity()

G4double G4StatMFMicroCanonical::CalcInvLevelDensity ( G4int  anA )

private

Definition at line 265 of file G4StatMFMicroCanonical.cc.

{
  G4double res = 0.0;
  if (anA > 1) {
    res = G4StatMFParameters::GetEpsilon0()*(1.0+3.0/(anA - 1.0));
  }
  return res;
}

References G4StatMFParameters::GetEpsilon0().

Referenced by CalcEntropyOfCompoundNucleus().

ChooseAandZ()

G4StatMFChannel * G4StatMFMicroCanonical::ChooseAandZ ( const G4Fragment &  theFragment )

virtual

Implements G4VStatMFEnsemble.

Definition at line 237 of file G4StatMFMicroCanonical.cc.

{
  // We choose a multiplicity (1,2,3,...) and then a channel
  G4double RandNumber = G4UniformRand();
 
  if (RandNumber < _WCompoundNucleus) { 
    
    G4StatMFChannel * aChannel = new G4StatMFChannel;
    aChannel->CreateFragment(theFragment.GetA_asInt(),theFragment.GetZ_asInt());
    return aChannel;
    
  } else {
    
    G4double AccumWeight = _WCompoundNucleus;
    std::vector<G4StatMFMicroManager*>::iterator it;
    for (it = _ThePartitionManagerVector.begin(); it != _ThePartitionManagerVector.end(); ++it) {
      AccumWeight += (*it)->GetProbability();
      if (RandNumber < AccumWeight) {
    return (*it)->ChooseChannel(theFragment.GetA_asInt(),theFragment.GetZ_asInt(),__MeanTemperature);
      }
    }
    throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMicroCanonical::ChooseAandZ: wrong normalization!");
  }
 
  return 0; 
}

References G4VStatMFEnsemble::__MeanTemperature, _ThePartitionManagerVector, _WCompoundNucleus, G4StatMFChannel::CreateFragment(), G4UniformRand, G4Fragment::GetA_asInt(), and G4Fragment::GetZ_asInt().

Referenced by G4StatMF::BreakItUp().

GetMeanMultiplicity()

G4double G4VStatMFEnsemble::GetMeanMultiplicity ( void  ) const

inlineinherited

Definition at line 68 of file G4VStatMFEnsemble.hh.

{return __MeanMultiplicity;}

References G4VStatMFEnsemble::__MeanMultiplicity.

Referenced by G4StatMF::BreakItUp().

GetMeanTemperature()

G4double G4VStatMFEnsemble::GetMeanTemperature ( void  ) const

inlineinherited

Definition at line 70 of file G4VStatMFEnsemble.hh.

{return __MeanTemperature;}

References G4VStatMFEnsemble::__MeanTemperature.

Referenced by G4StatMF::BreakItUp().

Initialize()

void G4StatMFMicroCanonical::Initialize ( const G4Fragment &  theFragment )

private

Definition at line 57 of file G4StatMFMicroCanonical.cc.

{
  
  std::vector<G4StatMFMicroManager*>::iterator it;
 
  // Excitation Energy 
  G4double U = theFragment.GetExcitationEnergy();
 
  G4int A = theFragment.GetA_asInt();
  G4int Z = theFragment.GetZ_asInt();
  G4double x = 1.0 - 2.0*Z/G4double(A);
  G4Pow* g4calc = G4Pow::GetInstance();
    
  // Configuration temperature
  G4double TConfiguration = std::sqrt(8.0*U/G4double(A));
  
  // Free internal energy at Temperature T = 0
  __FreeInternalE0 = A*( 
            // Volume term (for T = 0)
            -G4StatMFParameters::GetE0() +  
            // Symmetry term
            G4StatMFParameters::GetGamma0()*x*x 
            ) + 
    // Surface term (for T = 0)
    G4StatMFParameters::GetBeta0()*g4calc->Z23(A) + 
    // Coulomb term 
    elm_coupling*0.6*Z*Z/(G4StatMFParameters::Getr0()*g4calc->Z13(A));
  
  // Statistical weight
  G4double W = 0.0;
  
  // Mean breakup multiplicity
  __MeanMultiplicity = 0.0;
  
  // Mean channel temperature
  __MeanTemperature = 0.0;
  
  // Mean channel entropy
  __MeanEntropy = 0.0;
  
  // Calculate entropy of compound nucleus
  G4double SCompoundNucleus = CalcEntropyOfCompoundNucleus(theFragment,TConfiguration);
  
  // Statistical weight of compound nucleus
  _WCompoundNucleus = 1.0; 
  
  W += _WCompoundNucleus;
    
  // Maximal fragment multiplicity allowed in direct simulation
  G4int MaxMult = G4StatMFMicroCanonical::MaxAllowedMultiplicity;
  if (A > 110) MaxMult -= 1;
  
  for (G4int im = 2; im <= MaxMult; im++) {
    G4StatMFMicroManager * aMicroManager = 
      new G4StatMFMicroManager(theFragment,im,__FreeInternalE0,SCompoundNucleus);
    _ThePartitionManagerVector.push_back(aMicroManager);
  }
  
  // W is the total probability
  W = std::accumulate(_ThePartitionManagerVector.begin(),
              _ThePartitionManagerVector.end(),
              W, [](const G4double& running_total,
                            G4StatMFMicroManager*& manager)
                         {
                   return running_total + manager->GetProbability();
                 } );
  
  // Normalization of statistical weights
  for (it =  _ThePartitionManagerVector.begin(); it !=  _ThePartitionManagerVector.end(); ++it) 
    {
      (*it)->Normalize(W);
    }
 
  _WCompoundNucleus /= W;
  
  __MeanMultiplicity += 1.0 * _WCompoundNucleus;
  __MeanTemperature += TConfiguration * _WCompoundNucleus;
  __MeanEntropy += SCompoundNucleus * _WCompoundNucleus;
  
  for (it =  _ThePartitionManagerVector.begin(); it !=  _ThePartitionManagerVector.end(); ++it) 
    {
      __MeanMultiplicity += (*it)->GetMeanMultiplicity();
      __MeanTemperature += (*it)->GetMeanTemperature();
      __MeanEntropy += (*it)->GetMeanEntropy();
    }
  
  return;
}

References G4VStatMFEnsemble::__FreeInternalE0, G4VStatMFEnsemble::__MeanEntropy, G4VStatMFEnsemble::__MeanMultiplicity, G4VStatMFEnsemble::__MeanTemperature, _ThePartitionManagerVector, _WCompoundNucleus, A, CalcEntropyOfCompoundNucleus(), source.hepunit::elm_coupling, G4Fragment::GetA_asInt(), G4StatMFParameters::GetBeta0(), G4StatMFParameters::GetE0(), G4Fragment::GetExcitationEnergy(), G4StatMFParameters::GetGamma0(), G4Pow::GetInstance(), G4StatMFParameters::Getr0(), G4Fragment::GetZ_asInt(), MaxAllowedMultiplicity, Z, G4Pow::Z13(), and G4Pow::Z23().

Referenced by G4StatMFMicroCanonical().

operator!=()

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

private

operator=()

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

private

operator==()

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

private

Field Documentation

__FreeInternalE0

G4double G4VStatMFEnsemble::__FreeInternalE0

protectedinherited

Definition at line 75 of file G4VStatMFEnsemble.hh.

Referenced by CalcEntropyOfCompoundNucleus(), G4StatMFMacroCanonical::CalculateTemperature(), G4StatMFMacroCanonical::Initialize(), and Initialize().

__MeanEntropy

G4double G4VStatMFEnsemble::__MeanEntropy

protectedinherited

Definition at line 81 of file G4VStatMFEnsemble.hh.

Referenced by G4StatMFMacroCanonical::CalculateTemperature(), and Initialize().

__MeanMultiplicity

G4double G4VStatMFEnsemble::__MeanMultiplicity

protectedinherited

Definition at line 84 of file G4VStatMFEnsemble.hh.

Referenced by G4StatMFMacroCanonical::CalculateTemperature(), G4StatMFMacroCanonical::ChooseA(), G4VStatMFEnsemble::GetMeanMultiplicity(), and Initialize().

__MeanTemperature

G4double G4VStatMFEnsemble::__MeanTemperature

protectedinherited

Definition at line 78 of file G4VStatMFEnsemble.hh.

Referenced by G4StatMFMacroCanonical::CalculateTemperature(), ChooseAandZ(), G4StatMFMacroCanonical::ChooseZ(), G4VStatMFEnsemble::GetMeanTemperature(), and Initialize().

_ThePartitionManagerVector

std::vector<G4StatMFMicroManager*> G4StatMFMicroCanonical::_ThePartitionManagerVector

private

Definition at line 97 of file G4StatMFMicroCanonical.hh.

Referenced by ChooseAandZ(), Initialize(), and ~G4StatMFMicroCanonical().

_WCompoundNucleus

G4double G4StatMFMicroCanonical::_WCompoundNucleus

private

Definition at line 100 of file G4StatMFMicroCanonical.hh.

Referenced by ChooseAandZ(), and Initialize().

---

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

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