G4StatMFMicroPartition Class Reference

#include <G4StatMFMicroPartition.hh>

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

Detailed Description

Definition at line 41 of file G4StatMFMicroPartition.hh.

Constructor & Destructor Documentation

G4StatMFMicroPartition()

G4StatMFMicroPartition::G4StatMFMicroPartition ( const G4int  A, const G4double  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.

{};

Member Function Documentation

CalcPartitionProbability()

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

Definition at line 230 of file G4StatMFMicroPartition.cc.

{   
  G4double T = CalcPartitionTemperature(U,FreeInternalE0);
  if ( T <= 0.0) return _Probability = 0.0;
  _Temperature = T;
  
  
  // 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(static_cast<G4int>(_thePartition[i]));
      ProbA32 *= static_cast<G4double>(_thePartition[i])*
        std::sqrt(static_cast<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)
            * std::pow(static_cast<G4double>(_thePartition[i]),2.0/3.0);
        } 
    }
    
  // 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*(std::pow(static_cast<G4double>(_thePartition.size()),1./3.)-1.0)
                    /(G4StatMFParameters::Getr0()*std::pow(static_cast<G4double>(theA),1./3.)));
  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, std::log(ProbA32/Fact) +
                                     (_thePartition.size()-1.0)*std::log(FreeVolume/ThermalWaveLenght3) +
                                     1.5*(_thePartition.size()-1.0) - (3./2.)*std::log(theA));
  
  PartitionEntropy += std::log(ProbDegeneracy) + TranslationalS;
  _Entropy = PartitionEntropy;
    
  // And finally compute probability of fragment configuration
  G4double exponent = PartitionEntropy-SCompound;
  if (exponent > 700.0) exponent = 700.0;
  return _Probability = std::exp(exponent);
}

ChooseZ()

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

Definition at line 323 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));
        } 
          while (Zf < 0 || Zf > Af);
          FragmentsZ.push_back(Zf);
          SumZ += Zf;
    }
      ZBalance = static_cast<G4int>(Z0) - SumZ;
    } 
  while (std::abs(ZBalance) > 1.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 96 of file G4StatMFMicroPartition.hh.

    {
        return _Entropy;
    }

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 91 of file G4StatMFMicroPartition.hh.

    {
        return _Temperature;
    }

Normalize()

void G4StatMFMicroPartition::Normalize ( const G4double  Normalization )

inline

Definition at line 83 of file G4StatMFMicroPartition.hh.

    { _Probability /= Normalization; }

operator!=()

G4bool G4StatMFMicroPartition::operator!=

(

const G4StatMFMicroPartition &

right

)

const

Definition at line 60 of file G4StatMFMicroPartition.cc.

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

operator==()

G4bool G4StatMFMicroPartition::operator==

(

const G4StatMFMicroPartition &

right

)

const

Definition at line 53 of file G4StatMFMicroPartition.cc.

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

SetPartitionFragment()

void G4StatMFMicroPartition::SetPartitionFragment ( const G4int  anA )

inline

Definition at line 77 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