G4StatMFMicroPartition Class Reference

#include <G4StatMFMicroPartition.hh>

Collaboration diagram for G4StatMFMicroPartition:

Public Member Functions
	G4StatMFMicroPartition (G4int A, G4int Z)
	~G4StatMFMicroPartition ()
G4bool	operator== (const G4StatMFMicroPartition &right) const
G4bool	operator!= (const G4StatMFMicroPartition &right) const
G4StatMFChannel *	ChooseZ (G4int A0, G4int Z0, G4double MeanT)
G4double	GetProbability (void)
void	SetPartitionFragment (G4int anA)
void	Normalize (G4double Normalization)
G4double	CalcPartitionProbability (G4double U, G4double FreeInternalE0, G4double SCompound)
G4double	GetTemperature (void)
G4double	GetEntropy (void)

Private Member Functions
	G4StatMFMicroPartition ()
	G4StatMFMicroPartition (const G4StatMFMicroPartition &right)
G4StatMFMicroPartition &	operator= (const G4StatMFMicroPartition &right)
void	CoulombFreeEnergy (G4int anA)
G4double	CalcPartitionTemperature (G4double U, G4double FreeInternalE0)
G4double	GetPartitionEnergy (G4double T)
G4double	GetCoulombEnergy (void)
G4double	GetDegeneracyFactor (G4int A)
G4double	InvLevelDensity (G4double Af)

Private Attributes
G4int	theA
G4int	theZ
G4double	_Probability
G4double	_Temperature
G4double	_Entropy
std::vector< G4int >	_thePartition
std::vector< G4double >	_theCoulombFreeEnergy

Detailed Description

Definition at line 41 of file G4StatMFMicroPartition.hh.

Constructor & Destructor Documentation

G4StatMFMicroPartition() [1/3]

G4StatMFMicroPartition::G4StatMFMicroPartition ( G4int  A, G4int  Z  )

inline

Definition at line 45 of file G4StatMFMicroPartition.hh.

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

~G4StatMFMicroPartition()

G4StatMFMicroPartition::~G4StatMFMicroPartition ( )

inline

Definition at line 51 of file G4StatMFMicroPartition.hh.

{};

G4StatMFMicroPartition() [2/3]

G4StatMFMicroPartition::G4StatMFMicroPartition ( )

inlineprivate

Definition at line 56 of file G4StatMFMicroPartition.hh.

{};

Here is the call graph for this function:

G4StatMFMicroPartition() [3/3]

G4StatMFMicroPartition::G4StatMFMicroPartition ( const G4StatMFMicroPartition &  right )

private

Definition at line 42 of file G4StatMFMicroPartition.cc.

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

Here is the call graph for this function:

Member Function Documentation

CalcPartitionProbability()

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

Definition at line 229 of file G4StatMFMicroPartition.cc.

{   
  G4double T = CalcPartitionTemperature(U,FreeInternalE0);
  if ( T <= 0.0) return _Probability = 0.0;
  _Temperature = T;
  
  G4Pow* g4pow = 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) * g4pow->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*(g4pow->Z13(_thePartition.size())-1.0)
                    /(G4StatMFParameters::Getr0()*g4pow->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*g4pow->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);
}

Here is the call graph for this function:

Here is the caller graph for this function:

CalcPartitionTemperature()

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

private

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

Here is the call graph for this function:

Here is the caller graph for this function:

ChooseZ()

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

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

Here is the call graph for this function:

Here is the caller graph for this function:

CoulombFreeEnergy()

void G4StatMFMicroPartition::CoulombFreeEnergy ( G4int  anA )

private

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

Here is the call graph for this function:

Here is the caller graph for this function:

GetCoulombEnergy()

G4double G4StatMFMicroPartition::GetCoulombEnergy ( void  )

private

Definition at line 96 of file G4StatMFMicroPartition.cc.

{
  G4Pow* g4pow = G4Pow::GetInstance();
  G4double  CoulombFactor = 1.0/g4pow->A13(1.0+G4StatMFParameters::GetKappaCoulomb());  
            
  G4double CoulombEnergy = elm_coupling*0.6*theZ*theZ*CoulombFactor/
    (G4StatMFParameters::Getr0()*g4pow->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]*g4pow->Z23(_thePartition[i])/
      G4StatMFParameters::Getr0();
        
  return CoulombEnergy;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetDegeneracyFactor()

G4double G4StatMFMicroPartition::GetDegeneracyFactor ( G4int  A )

private

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

Here is the caller graph for this function:

GetEntropy()

G4double G4StatMFMicroPartition::GetEntropy ( void  )

inline

Definition at line 93 of file G4StatMFMicroPartition.hh.

    {
        return _Entropy;
    }

Here is the call graph for this function:

Here is the caller graph for this function:

GetPartitionEnergy()

G4double G4StatMFMicroPartition::GetPartitionEnergy ( G4double  T )

private

Definition at line 113 of file G4StatMFMicroPartition.cc.

{
  G4Pow* g4pow = G4Pow::GetInstance();
  G4double  CoulombFactor = 1.0/g4pow->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))*
            g4pow->Z23(_thePartition[i]) +
            
            // Coulomb term 
            _theCoulombFreeEnergy[i];
    }
    }
    
  PartitionEnergy += elm_coupling*0.6*theZ*theZ*CoulombFactor/
    (G4StatMFParameters::Getr0()*g4pow->Z13(theA))
    + 1.5*T*(_thePartition.size()-1);
  
  return PartitionEnergy;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetProbability()

G4double G4StatMFMicroPartition::GetProbability ( void  )

inline

Definition at line 72 of file G4StatMFMicroPartition.hh.

    { return _Probability; }

GetTemperature()

G4double G4StatMFMicroPartition::GetTemperature ( void  )

inline

Definition at line 88 of file G4StatMFMicroPartition.hh.

    {
        return _Temperature;
    }

Here is the caller graph for this function:

InvLevelDensity()

G4double G4StatMFMicroPartition::InvLevelDensity ( G4double  Af )

inlineprivate

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

Here is the call graph for this function:

Here is the caller graph for this function:

Normalize()

void G4StatMFMicroPartition::Normalize ( G4double  Normalization )

inline

Definition at line 81 of file G4StatMFMicroPartition.hh.

    { _Probability /= Normalization; }

Here is the call graph for this function:

operator!=()

G4bool G4StatMFMicroPartition::operator!=

(

const G4StatMFMicroPartition &

right

)

const

Definition at line 64 of file G4StatMFMicroPartition.cc.

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

Here is the caller graph for this function:

operator=()

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

private

Definition at line 50 of file G4StatMFMicroPartition.cc.

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

Here is the caller graph for this function:

operator==()

G4bool G4StatMFMicroPartition::operator==

(

const G4StatMFMicroPartition &

right

)

const

Definition at line 57 of file G4StatMFMicroPartition.cc.

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

Here is the caller graph for this function:

SetPartitionFragment()

void G4StatMFMicroPartition::SetPartitionFragment ( G4int  anA )

inline

Definition at line 75 of file G4StatMFMicroPartition.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

_Entropy

G4double G4StatMFMicroPartition::_Entropy

private

Definition at line 132 of file G4StatMFMicroPartition.hh.

_Probability

G4double G4StatMFMicroPartition::_Probability

private

Definition at line 126 of file G4StatMFMicroPartition.hh.

_Temperature

G4double G4StatMFMicroPartition::_Temperature

private

Definition at line 129 of file G4StatMFMicroPartition.hh.

_theCoulombFreeEnergy

std::vector<G4double> G4StatMFMicroPartition::_theCoulombFreeEnergy

private

Definition at line 137 of file G4StatMFMicroPartition.hh.

_thePartition

std::vector<G4int> G4StatMFMicroPartition::_thePartition

private

Definition at line 135 of file G4StatMFMicroPartition.hh.

theA

G4int G4StatMFMicroPartition::theA

private

Definition at line 122 of file G4StatMFMicroPartition.hh.

theZ

G4int G4StatMFMicroPartition::theZ

private

Definition at line 123 of file G4StatMFMicroPartition.hh.

---

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

• G4StatMFMicroPartition.hh
• G4StatMFMicroPartition.cc