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

Referenced by operator!=(), and operator==().

~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();
  G4int n = (G4int)_thePartition.size();
  
  // Factorial of fragment multiplicity
  G4double Fact = g4calc->factorial(n);
    
  G4double ProbDegeneracy = 1.0;
  G4double ProbA32 = 1.0;   
  G4double PartitionEntropy = 0.0;
  G4double db = G4StatMFParameters::DBetaDT(T);
    
  for (G4int i = 0; i < n; ++i) {
    G4int par = _thePartition[i];
    ProbDegeneracy *= GetDegeneracyFactor(par);
    ProbA32 *= _thePartition[i]*std::sqrt((G4double)par);
 
    // interaction entropy for alpha
    if (par == 4) {
      PartitionEntropy += 2.0 * T * par/InvLevelDensity(par);
    }
    // interaction entropy for Af > 4
    else if (par > 4) {
      PartitionEntropy += 2.0 * T * par/InvLevelDensity(par) - db * g4calc->Z23(par);
    } 
  }
    
  // 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 r0 = G4StatMFParameters::Getr0();
  G4double kappa = 1. + elm_coupling*(g4calc->Z13(n) - 1.0)/(r0*g4calc->Z13(theA));
  G4double V0 = (4./3.)*pi*theA*r0*r0*r0;
  G4double FreeVolume = (kappa*kappa*kappa - 1.0)*V0;
  G4double TranslationalS = G4Log(ProbA32/Fact)
    + (n - 1)*G4Log(FreeVolume/ThermalWaveLenght3)
    + 1.5*(n - 1) - 1.5*g4calc->logZ(theA);
  TranslationalS = std::max(TranslationalS, 0.0);
  
  PartitionEntropy += G4Log(ProbDegeneracy) + TranslationalS;
  _Entropy = PartitionEntropy;
    
  // And finally compute probability of fragment configuration
  G4double exponent = std::min(PartitionEntropy - SCompound, 200.);
  return _Probability = G4Exp(exponent);
}

ChooseZ()

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

Definition at line 297 of file G4StatMFMicroPartition.cc.

{
  std::vector<G4int> FragmentsZ;
  G4int n = (G4int)_thePartition.size();
  
  G4int ZBalance = 0;
  do {
    G4double CC = G4StatMFParameters::GetGamma0()*8.0;
    G4int SumZ = 0;
    for (G4int i = 0; i < n; ++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 (G4int i = 0; i < n; ++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