#include <G4StatMFMicroPartition.hh>

Public Member Functions
G4double	CalcPartitionProbability (G4double U, G4double FreeInternalE0, G4double SCompound)

G4StatMFChannel *	ChooseZ (G4int A0, G4int Z0, G4double MeanT)

	G4StatMFMicroPartition (G4int A, G4int Z)

G4double	GetEntropy (void)

G4double	GetProbability (void)

G4double	GetTemperature (void)

void	Normalize (G4double Normalization)

G4bool	operator!= (const G4StatMFMicroPartition &right) const

G4bool	operator== (const G4StatMFMicroPartition &right) const

void	SetPartitionFragment (G4int anA)

	~G4StatMFMicroPartition ()

Private Member Functions
G4double	CalcPartitionTemperature (G4double U, G4double FreeInternalE0)

void	CoulombFreeEnergy (G4int anA)

	G4StatMFMicroPartition ()

	G4StatMFMicroPartition (const G4StatMFMicroPartition &right)

G4double	GetCoulombEnergy (void)

G4double	GetDegeneracyFactor (G4int A)

G4double	GetPartitionEnergy (G4double T)

G4double	InvLevelDensity (G4double Af)

G4StatMFMicroPartition &	operator= (const G4StatMFMicroPartition &right)

Private Attributes
G4double	_Entropy

G4double	_Probability

G4double	_Temperature

std::vector< G4double >	_theCoulombFreeEnergy

std::vector< G4int >	_thePartition

G4int	theA

G4int	theZ

Detailed Description

Definition at line 40 of file G4StatMFMicroPartition.hh.

Constructor & Destructor Documentation

◆ G4StatMFMicroPartition() [1/3]

G4StatMFMicroPartition::G4StatMFMicroPartition	(	G4int	A,
		G4int	Z
	)

inline

Definition at line 44 of file G4StatMFMicroPartition.hh.

                                             :
    theA(A), theZ(Z), _Probability(0.0), _Temperature(0.0), 
    _Entropy(0.0) {};

◆ ~G4StatMFMicroPartition()

G4StatMFMicroPartition::~G4StatMFMicroPartition ( )

inline

Definition at line 50 of file G4StatMFMicroPartition.hh.

50{};

◆ G4StatMFMicroPartition() [2/3]

G4StatMFMicroPartition::G4StatMFMicroPartition ( )

inlineprivate

Definition at line 55 of file G4StatMFMicroPartition.hh.

55{};

◆ G4StatMFMicroPartition() [3/3]

G4StatMFMicroPartition::G4StatMFMicroPartition ( const G4StatMFMicroPartition & right )

private

Definition at line 41 of file G4StatMFMicroPartition.cc.

{
  throw G4HadronicException(__FILE__, __LINE__, "G4StatMFMicroPartition::copy_constructor meant to not be accessible");
}

Member Function Documentation

◆ CalcPartitionProbability()

G4double G4StatMFMicroPartition::CalcPartitionProbability	(	G4double	U,
		G4double	FreeInternalE0,
		G4double	SCompound
	)

Definition at line 228 of file G4StatMFMicroPartition.cc.

{   
  G4double T = CalcPartitionTemperature(U,FreeInternalE0);
  if ( T <= 0.0) return _Probability = 0.0;
  _Temperature = T;
  
  G4Pow* g4calc = G4Pow::GetInstance();
  
  // Factorial of fragment multiplicity
  G4double Fact = 1.0;
  unsigned int i;
  for (i = 0; i < _thePartition.size() - 1; i++) 
    {
      G4double f = 1.0;
      for (unsigned int ii = i+1; i< _thePartition.size(); i++) 
        {
          if (_thePartition[i] == _thePartition[ii]) f++;
        }
      Fact *= f;
  }
    
  G4double ProbDegeneracy = 1.0;
  G4double ProbA32 = 1.0;   
    
  for (i = 0; i < _thePartition.size(); i++) 
    {
      ProbDegeneracy *= GetDegeneracyFactor(_thePartition[i]);
      ProbA32 *= _thePartition[i]*std::sqrt((G4double)_thePartition[i]);
    }
    
  // Compute entropy
  G4double PartitionEntropy = 0.0;
  for (i = 0; i < _thePartition.size(); i++) 
    {
      // interaction entropy for alpha
      if (_thePartition[i] == 4) 
        {
          PartitionEntropy += 
            2.0*T*_thePartition[i]/InvLevelDensity(_thePartition[i]);
        }
      // interaction entropy for Af > 4
      else if (_thePartition[i] > 4) 
        {
          PartitionEntropy += 
            2.0*T*_thePartition[i]/InvLevelDensity(_thePartition[i])
            - G4StatMFParameters::DBetaDT(T) * g4calc->Z23(_thePartition[i]);
        } 
    }
    
  // Thermal Wave Lenght = std::sqrt(2 pi hbar^2 / nucleon_mass T)
  G4double ThermalWaveLenght3 = 16.15*fermi/std::sqrt(T);
  ThermalWaveLenght3 = ThermalWaveLenght3*ThermalWaveLenght3*ThermalWaveLenght3;
  
  // Translational Entropy
  G4double kappa = 1. + elm_coupling*(g4calc->Z13(_thePartition.size())-1.0)
                    /(G4StatMFParameters::Getr0()*g4calc->Z13(theA));
  kappa = kappa*kappa*kappa;
  kappa -= 1.;
  G4double V0 = (4./3.)*pi*theA*G4StatMFParameters::Getr0()*G4StatMFParameters::Getr0()*
    G4StatMFParameters::Getr0();
  G4double FreeVolume = kappa*V0;
  G4double TranslationalS = std::max(0.0, G4Log(ProbA32/Fact) +
      (_thePartition.size()-1.0)*G4Log(FreeVolume/ThermalWaveLenght3) +
      1.5*(_thePartition.size()-1.0) - 1.5*g4calc->logZ(theA));
  
  PartitionEntropy += G4Log(ProbDegeneracy) + TranslationalS;
  _Entropy = PartitionEntropy;
    
  // And finally compute probability of fragment configuration
  G4double exponent = PartitionEntropy-SCompound;
  if (exponent > 300.0) exponent = 300.0;
  return _Probability = G4Exp(exponent);
}

References _Entropy, _Probability, _Temperature, _thePartition, CalcPartitionTemperature(), G4StatMFParameters::DBetaDT(), source.hepunit::elm_coupling, fermi, G4Exp(), G4Log(), GetDegeneracyFactor(), G4Pow::GetInstance(), G4StatMFParameters::Getr0(), InvLevelDensity(), G4Pow::logZ(), G4INCL::Math::max(), pi, theA, G4Pow::Z13(), and G4Pow::Z23().

Referenced by G4StatMFMicroManager::Initialize().

◆ CalcPartitionTemperature()

G4double G4StatMFMicroPartition::CalcPartitionTemperature	(	G4double	U,
		G4double	FreeInternalE0
	)

private

Definition at line 173 of file G4StatMFMicroPartition.cc.

{
  G4double PartitionEnergy = GetPartitionEnergy(0.0);
  
  // If this happens, T = 0 MeV, which means that probability for this
  // partition will be 0
  if (std::fabs(U + FreeInternalE0 - PartitionEnergy) < 0.003) return -1.0;
    
  // Calculate temperature by midpoint method
    
  // Bracketing the solution
  G4double Ta = 0.001;
  G4double Tb = std::max(std::sqrt(8.0*U/theA),0.0012*MeV);
  G4double Tmid = 0.0;
  
  G4double Da = (U + FreeInternalE0 - GetPartitionEnergy(Ta))/U;
  G4double Db = (U + FreeInternalE0 - GetPartitionEnergy(Tb))/U;
  
  G4int maxit = 0;
  // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
  while (Da*Db > 0.0 && maxit < 1000) 
    {
      ++maxit;
      Tb += 0.5*Tb;     
      Db = (U + FreeInternalE0 - GetPartitionEnergy(Tb))/U;
    }
  
  G4double eps = 1.0e-14*std::abs(Ta-Tb);
  
  for (G4int i = 0; i < 1000; i++) 
    {
      Tmid = (Ta+Tb)/2.0;
      if (std::fabs(Ta-Tb) <= eps) return Tmid;
      G4double Dmid = (U + FreeInternalE0 - GetPartitionEnergy(Tmid))/U;
      if (std::fabs(Dmid) < 0.003) return Tmid;
      if (Da*Dmid < 0.0) 
        {
          Tb = Tmid;
          Db = Dmid;
        } 
      else 
        {
          Ta = Tmid;
          Da = Dmid;
        } 
    }
  // if we arrive here the temperature could not be calculated
  G4cout << "G4StatMFMicroPartition::CalcPartitionTemperature: I can't calculate the temperature"  
         << G4endl;
  // and set probability to 0 returning T < 0
  return -1.0;
  
}

References eps, G4cout, G4endl, GetPartitionEnergy(), G4INCL::Math::max(), MeV, and theA.

Referenced by CalcPartitionProbability().

◆ ChooseZ()

G4StatMFChannel * G4StatMFMicroPartition::ChooseZ	(	G4int	A0,
		G4int	Z0,
		G4double	MeanT
	)

Definition at line 317 of file G4StatMFMicroPartition.cc.

{
  std::vector<G4int> FragmentsZ;
  
  G4int ZBalance = 0;
  do 
    {
      G4double CC = G4StatMFParameters::GetGamma0()*8.0;
      G4int SumZ = 0;
      for (unsigned int i = 0; i < _thePartition.size(); i++) 
        {
          G4double ZMean;
          G4double Af = _thePartition[i];
          if (Af > 1.5 && Af < 4.5) ZMean = 0.5*Af;
          else ZMean = Af*Z0/A0;
          G4double ZDispersion = std::sqrt(Af * MeanT/CC);
          G4int Zf;
          do 
            {
              Zf = static_cast<G4int>(G4RandGauss::shoot(ZMean,ZDispersion));
        } 
          // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
          while (Zf < 0 || Zf > Af);
          FragmentsZ.push_back(Zf);
          SumZ += Zf;
    }
      ZBalance = Z0 - SumZ;
    } 
  // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
  while (std::abs(ZBalance) > 1);
  FragmentsZ[0] += ZBalance;
    
  G4StatMFChannel * theChannel = new G4StatMFChannel;
  for (unsigned int i = 0; i < _thePartition.size(); i++)
    {
      theChannel->CreateFragment(_thePartition[i],FragmentsZ[i]);
    }
 
  return theChannel;
}

References _thePartition, G4StatMFChannel::CreateFragment(), G4StatMFParameters::GetGamma0(), and G4INCL::DeJongSpin::shoot().

◆ CoulombFreeEnergy()

void G4StatMFMicroPartition::CoulombFreeEnergy ( G4int anA )

private

Definition at line 69 of file G4StatMFMicroPartition.cc.

{
  // This Z independent factor in the Coulomb free energy 
  G4double  CoulombConstFactor = G4StatMFParameters::GetCoulomb();
 
  // We use the aproximation Z_f ~ Z/A * A_f
 
  G4double ZA = G4double(theZ)/G4double(theA);
                                        
  if (anA == 0 || anA == 1) 
    {
      _theCoulombFreeEnergy.push_back(CoulombConstFactor*ZA*ZA);
    } 
  else if (anA == 2 || anA == 3 || anA == 4) 
    {
      // Z/A ~ 1/2
      _theCoulombFreeEnergy.push_back(CoulombConstFactor*0.5
                      *anA*G4Pow::GetInstance()->Z23(anA));
    } 
  else  // anA > 4
    {
      _theCoulombFreeEnergy.push_back(CoulombConstFactor*ZA*ZA  
                      *anA*G4Pow::GetInstance()->Z23(anA));
    }
}

References _theCoulombFreeEnergy, G4StatMFParameters::GetCoulomb(), G4Pow::GetInstance(), theA, and theZ.

Referenced by SetPartitionFragment().

◆ GetCoulombEnergy()

G4double G4StatMFMicroPartition::GetCoulombEnergy ( void )

private

Definition at line 95 of file G4StatMFMicroPartition.cc.

{
  G4Pow* g4calc = G4Pow::GetInstance();
  G4double  CoulombFactor = 1.0/g4calc->A13(1.0+G4StatMFParameters::GetKappaCoulomb()); 
            
  G4double CoulombEnergy = elm_coupling*0.6*theZ*theZ*CoulombFactor/
    (G4StatMFParameters::Getr0()*g4calc->Z13(theA));
    
  G4double ZA = G4double(theZ)/G4double(theA);
  for (unsigned int i = 0; i < _thePartition.size(); i++) 
    CoulombEnergy += _theCoulombFreeEnergy[i] - elm_coupling*0.6*
      ZA*ZA*_thePartition[i]*g4calc->Z23(_thePartition[i])/
      G4StatMFParameters::Getr0();
        
  return CoulombEnergy;
}

References _theCoulombFreeEnergy, _thePartition, G4Pow::A13(), source.hepunit::elm_coupling, G4Pow::GetInstance(), G4StatMFParameters::GetKappaCoulomb(), G4StatMFParameters::Getr0(), theA, theZ, G4Pow::Z13(), and G4Pow::Z23().

◆ GetDegeneracyFactor()

G4double G4StatMFMicroPartition::GetDegeneracyFactor ( G4int A )

private

Definition at line 304 of file G4StatMFMicroPartition.cc.

{
  // Degeneracy factors are statistical factors
  // DegeneracyFactor for nucleon is (2S_n + 1)(2I_n + 1) = 4
  G4double DegFactor = 0;
  if (A > 4) DegFactor = 1.0;
  else if (A == 1) DegFactor = 4.0;     // nucleon
  else if (A == 2) DegFactor = 3.0;     // Deuteron
  else if (A == 3) DegFactor = 4.0;     // Triton + He3
  else if (A == 4) DegFactor = 1.0;     // alpha
  return DegFactor;
}

References A.

Referenced by CalcPartitionProbability().

◆ GetEntropy()

G4double G4StatMFMicroPartition::GetEntropy ( void )

inline

Definition at line 92 of file G4StatMFMicroPartition.hh.

    {
        return _Entropy;
    }

References _Entropy.

Referenced by G4StatMFMicroManager::Initialize().

◆ GetPartitionEnergy()

G4double G4StatMFMicroPartition::GetPartitionEnergy ( G4double T )

private

Definition at line 112 of file G4StatMFMicroPartition.cc.

{
  G4Pow* g4calc = G4Pow::GetInstance();
  G4double  CoulombFactor = 1.0/g4calc->A13(1.0+G4StatMFParameters::GetKappaCoulomb()); 
  
  G4double PartitionEnergy = 0.0;
  
  // We use the aprox that Z_f ~ Z/A * A_f
  for (unsigned int i = 0; i < _thePartition.size(); i++) 
    {
      if (_thePartition[i] == 0 || _thePartition[i] == 1) 
        {   
          PartitionEnergy += _theCoulombFreeEnergy[i];
        }
      else if (_thePartition[i] == 2) 
        {       
          PartitionEnergy +=    
            -2.796 // Binding Energy of deuteron ??????
            + _theCoulombFreeEnergy[i];     
    }
      else if (_thePartition[i] == 3) 
        {   
          PartitionEnergy +=    
            -9.224 // Binding Energy of trtion/He3 ??????
            + _theCoulombFreeEnergy[i];     
    } 
      else if (_thePartition[i] == 4) 
        {   
          PartitionEnergy +=
            -30.11 // Binding Energy of ALPHA ??????
            + _theCoulombFreeEnergy[i] 
            + 4.*T*T/InvLevelDensity(4.);
    } 
      else 
        {
          PartitionEnergy +=
            //Volume term                       
            (- G4StatMFParameters::GetE0() + 
             T*T/InvLevelDensity(_thePartition[i]))
            *_thePartition[i] + 
            
            // Symmetry term
            G4StatMFParameters::GetGamma0()*
            (1.0-2.0*theZ/theA)*(1.0-2.0*theZ/theA)*_thePartition[i] +  
            
            // Surface term
            (G4StatMFParameters::Beta(T) - T*G4StatMFParameters::DBetaDT(T))*
            g4calc->Z23(_thePartition[i]) +
            
            // Coulomb term 
            _theCoulombFreeEnergy[i];
    }
    }
    
  PartitionEnergy += elm_coupling*0.6*theZ*theZ*CoulombFactor/
    (G4StatMFParameters::Getr0()*g4calc->Z13(theA))
    + 1.5*T*(_thePartition.size()-1);
  
  return PartitionEnergy;
}

References _theCoulombFreeEnergy, _thePartition, G4Pow::A13(), G4StatMFParameters::Beta(), G4StatMFParameters::DBetaDT(), source.hepunit::elm_coupling, G4StatMFParameters::GetE0(), G4StatMFParameters::GetGamma0(), G4Pow::GetInstance(), G4StatMFParameters::GetKappaCoulomb(), G4StatMFParameters::Getr0(), InvLevelDensity(), theA, theZ, G4Pow::Z13(), and G4Pow::Z23().

Referenced by CalcPartitionTemperature().

◆ GetProbability()

G4double G4StatMFMicroPartition::GetProbability ( void )

inline

Definition at line 71 of file G4StatMFMicroPartition.hh.

72 { return _Probability; }

References _Probability.

◆ GetTemperature()

G4double G4StatMFMicroPartition::GetTemperature ( void )

inline

Definition at line 87 of file G4StatMFMicroPartition.hh.

    {
        return _Temperature;
    }

References _Temperature.

Referenced by G4StatMFMicroManager::Initialize().

◆ InvLevelDensity()

G4double G4StatMFMicroPartition::InvLevelDensity ( G4double Af )

inlineprivate

Definition at line 110 of file G4StatMFMicroPartition.hh.

    {
        // Calculate Inverse Density Level
        // Epsilon0*(1 + 3 /(Af - 1))
        if (Af < 1.5) return 0.0;
        else return G4StatMFParameters::GetEpsilon0()*(1.0+3.0/(Af - 1.0));
    }

References G4StatMFParameters::GetEpsilon0().

Referenced by CalcPartitionProbability(), and GetPartitionEnergy().

◆ Normalize()

void G4StatMFMicroPartition::Normalize ( G4double Normalization )

inline

Definition at line 80 of file G4StatMFMicroPartition.hh.

81 { _Probability /= Normalization; }

References _Probability.

◆ operator!=()

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

Definition at line 63 of file G4StatMFMicroPartition.cc.

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

◆ operator=()

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

private

Definition at line 48 of file G4StatMFMicroPartition.cc.

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

◆ operator==()

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

Definition at line 56 of file G4StatMFMicroPartition.cc.

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

◆ SetPartitionFragment()

void G4StatMFMicroPartition::SetPartitionFragment ( G4int anA )

inline

Definition at line 74 of file G4StatMFMicroPartition.hh.

    { 
        _thePartition.push_back(anA);
        CoulombFreeEnergy(anA);
    }

References _thePartition, and CoulombFreeEnergy().

Referenced by G4StatMFMicroManager::Initialize().

Field Documentation

◆ _Entropy

G4double G4StatMFMicroPartition::_Entropy

private

Definition at line 131 of file G4StatMFMicroPartition.hh.

Referenced by CalcPartitionProbability(), and GetEntropy().

◆ _Probability

G4double G4StatMFMicroPartition::_Probability

private

Definition at line 125 of file G4StatMFMicroPartition.hh.

Referenced by CalcPartitionProbability(), GetProbability(), and Normalize().

◆ _Temperature

G4double G4StatMFMicroPartition::_Temperature

private

Definition at line 128 of file G4StatMFMicroPartition.hh.

Referenced by CalcPartitionProbability(), and GetTemperature().

◆ _theCoulombFreeEnergy

std::vector<G4double> G4StatMFMicroPartition::_theCoulombFreeEnergy

private

Definition at line 136 of file G4StatMFMicroPartition.hh.

Referenced by CoulombFreeEnergy(), GetCoulombEnergy(), and GetPartitionEnergy().

◆ _thePartition

std::vector<G4int> G4StatMFMicroPartition::_thePartition

private

Definition at line 134 of file G4StatMFMicroPartition.hh.

Referenced by CalcPartitionProbability(), ChooseZ(), GetCoulombEnergy(), GetPartitionEnergy(), and SetPartitionFragment().

◆ theA

G4int G4StatMFMicroPartition::theA

private

Definition at line 121 of file G4StatMFMicroPartition.hh.

Referenced by CalcPartitionProbability(), CalcPartitionTemperature(), CoulombFreeEnergy(), GetCoulombEnergy(), and GetPartitionEnergy().

◆ theZ

G4int G4StatMFMicroPartition::theZ

private

Definition at line 122 of file G4StatMFMicroPartition.hh.

Referenced by CoulombFreeEnergy(), GetCoulombEnergy(), and GetPartitionEnergy().

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

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

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4StatMFMicroPartition() [1/3]

◆ ~G4StatMFMicroPartition()

◆ G4StatMFMicroPartition() [2/3]

◆ G4StatMFMicroPartition() [3/3]

Member Function Documentation

◆ CalcPartitionProbability()

◆ CalcPartitionTemperature()

◆ ChooseZ()

◆ CoulombFreeEnergy()

◆ GetCoulombEnergy()

◆ GetDegeneracyFactor()

◆ GetEntropy()

◆ GetPartitionEnergy()

◆ GetProbability()

◆ GetTemperature()

◆ InvLevelDensity()

◆ Normalize()

◆ operator!=()

◆ operator=()

◆ operator==()

◆ SetPartitionFragment()

Field Documentation

◆ _Entropy

◆ _Probability

◆ _Temperature

◆ _theCoulombFreeEnergy

◆ _thePartition

◆ theA

◆ theZ