G4StatMFMicroPartition Class Reference

#include <G4StatMFMicroPartition.hh>

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)

Detailed Description

Definition at line 40 of file G4StatMFMicroPartition.hh.

Constructor & Destructor Documentation

G4StatMFMicroPartition()

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.

{};

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((G4int)_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);
}

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;
}

GetEntropy()

G4double G4StatMFMicroPartition::GetEntropy ( void )

inline

Definition at line 92 of file G4StatMFMicroPartition.hh.

    {
        return _Entropy;
    }

GetProbability()

G4double G4StatMFMicroPartition::GetProbability ( void )

inline

Definition at line 71 of file G4StatMFMicroPartition.hh.

    { return _Probability; }

GetTemperature()

G4double G4StatMFMicroPartition::GetTemperature ( void )

inline

Definition at line 87 of file G4StatMFMicroPartition.hh.

    {
        return _Temperature;
    }

Normalize()

void G4StatMFMicroPartition::Normalize ( G4double Normalization )

inline

Definition at line 80 of file G4StatMFMicroPartition.hh.

    { _Probability /= Normalization; }

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

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);
    }

---

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

• G4StatMFMicroPartition.hh
• G4StatMFMicroPartition.cc