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

Constructor & Destructor Documentation

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

inline

Definition at line 51 of file G4StatMFMicroPartition.hh.

{};

Member Function Documentation

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

Definition at line 225 of file G4StatMFMicroPartition.cc.

References G4StatMFParameters::DBetaDT(), python.hepunit::elm_coupling, python.hepunit::fermi, G4StatMFParameters::Getr0(), G4INCL::Math::max(), and python.hepunit::pi.

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

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

Definition at line 315 of file G4StatMFMicroPartition.cc.

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

{
  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 = Z0 - SumZ;
    } 
  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;
}

G4double G4StatMFMicroPartition::GetEntropy ( void  )

inline

Definition at line 93 of file G4StatMFMicroPartition.hh.

    {
        return _Entropy;
    }

G4double G4StatMFMicroPartition::GetProbability ( void  )

inline

Definition at line 72 of file G4StatMFMicroPartition.hh.

    { return _Probability; }

G4double G4StatMFMicroPartition::GetTemperature ( void  )

inline

Definition at line 88 of file G4StatMFMicroPartition.hh.

    {
        return _Temperature;
    }

void G4StatMFMicroPartition::Normalize ( G4double  Normalization )

inline

Definition at line 81 of file G4StatMFMicroPartition.hh.

    { _Probability /= Normalization; }

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

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

void G4StatMFMicroPartition::SetPartitionFragment ( G4int  anA )

inline

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