G4StatMFMacroCanonical Class Reference

#include <G4StatMFMacroCanonical.hh>

Inheritance diagram for G4StatMFMacroCanonical:

Public Member Functions
	G4StatMFMacroCanonical ()
	~G4StatMFMacroCanonical () override
void	Initialise (const G4Fragment &theFragment) override
G4StatMFChannel *	ChooseAandZ (const G4Fragment &theFragment) override
	G4StatMFMacroCanonical (const G4StatMFMacroCanonical &)=delete
G4StatMFMacroCanonical &	operator= (const G4StatMFMacroCanonical &right)=delete
G4bool	operator== (const G4StatMFMacroCanonical &right) const =delete
G4bool	operator!= (const G4StatMFMacroCanonical &right) const =delete
Public Member Functions inherited from G4VStatMFEnsemble
	G4VStatMFEnsemble ()=default
virtual	~G4VStatMFEnsemble ()=default
G4double	GetMeanMultiplicity () const
G4double	GetMeanTemperature () const
	G4VStatMFEnsemble (const G4VStatMFEnsemble &right)=delete
G4VStatMFEnsemble &	operator= (const G4VStatMFEnsemble &right)=delete
G4bool	operator== (const G4VStatMFEnsemble &right) const =delete
G4bool	operator!= (const G4VStatMFEnsemble &right) const =delete

Additional Inherited Members
Protected Attributes inherited from G4VStatMFEnsemble
G4double	pFreeInternalE0 {0.0}
G4double	pMeanTemperature {0.0}
G4double	pMeanEntropy {0.0}
G4double	pMeanMultiplicity {0.0}

Detailed Description

Definition at line 49 of file G4StatMFMacroCanonical.hh.

Constructor & Destructor Documentation

G4StatMFMacroCanonical() [1/2]

G4StatMFMacroCanonical::G4StatMFMacroCanonical

(

)

Definition at line 43 of file G4StatMFMacroCanonical.cc.

{
  theTemp = new G4StatMFMacroTemperature();
  fClusters.reserve(220);
  fAcumMultiplicity.reserve(220);
 
  // Get memory for clusters
  fClusters.push_back(new G4StatMFMacroNucleon);
  fClusters.push_back(new G4StatMFMacroBiNucleon);
  fClusters.push_back(new G4StatMFMacroTriNucleon);
  fClusters.push_back(new G4StatMFMacroTetraNucleon);
}

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

~G4StatMFMacroCanonical()

G4StatMFMacroCanonical::~G4StatMFMacroCanonical ( )

override

Definition at line 57 of file G4StatMFMacroCanonical.cc.

{
  // garbage collection
  if (!fClusters.empty()) {
    for(auto const & p : fClusters) { delete p; }
  }
  delete theTemp;
}

G4StatMFMacroCanonical() [2/2]

G4StatMFMacroCanonical::G4StatMFMacroCanonical ( const G4StatMFMacroCanonical & )

delete

Member Function Documentation

ChooseAandZ()

G4StatMFChannel * G4StatMFMacroCanonical::ChooseAandZ ( const G4Fragment & theFragment )

overridevirtual

Implements G4VStatMFEnsemble.

Definition at line 110 of file G4StatMFMacroCanonical.cc.

{
  G4int A = theFragment.GetA_asInt();
  G4int Z = theFragment.GetZ_asInt();
  
  std::vector<G4int> ANumbers(A);
  
  G4double Multiplicity = ChooseA(A, ANumbers);
  
  std::vector<G4int> FragmentsA;
  
  G4int i = 0;  
  for (i = 0; i < A; i++) 
    {
      for (G4int j = 0; j < ANumbers[i]; j++) FragmentsA.push_back(i+1);
    }
  
  // Sort fragments in decreasing order
  G4int im = 0;
  for (G4int j = 0; j < Multiplicity; j++) 
    {
      G4int FragmentsAMax = 0;
      im = j;
      for (i = j; i < Multiplicity; ++i) 
    {
      if (FragmentsA[i] <= FragmentsAMax) { continue; }
      else 
        {
          im = i;
          FragmentsAMax = FragmentsA[im];
        }
    }   
      if (im != j) 
    {
      FragmentsA[im] = FragmentsA[j];
      FragmentsA[j]  = FragmentsAMax;
    }
    }
  return ChooseZ(Z,FragmentsA);
}

Initialise()

void G4StatMFMacroCanonical::Initialise ( const G4Fragment & theFragment )

overridevirtual

Implements G4VStatMFEnsemble.

Definition at line 67 of file G4StatMFMacroCanonical.cc.

{
  const G4double R0 = G4StatMFParameters::Getr0();
  const G4double E0 = G4StatMFParameters::GetE0();
  const G4double Beta0 = G4StatMFParameters::GetBeta0();
  const G4double Gamma0 = G4StatMFParameters::GetGamma0();
 
  G4int A = theFragment.GetA_asInt();
  G4int Z = theFragment.GetZ_asInt();
  G4double x = 1.0 - (2*Z)/G4double(A);
  G4Pow* g4calc = G4Pow::GetInstance();
  G4double a13 = g4calc->Z13(A);
  
  // Free Internal energy at T = 0
  pFreeInternalE0 = A*(Gamma0*x*x - E0) // Symmetry term & Volume term (for T = 0)
    + Beta0*a13*a13 +   // Surface term (for T = 0)
    0.6*CLHEP::elm_coupling*(Z*Z)/(R0*a13); // Coulomb term 
 
  G4int n = (G4int)fClusters.size();
  for (G4int i = n; i < A; ++i) {
    fClusters.push_back(new G4StatMFMacroMultiNucleon(i+1)); // Size 5 ... A
  }
  // Excitation Energy
  G4double U = theFragment.GetExcitationEnergy();
  
  // Fragment Multiplicity
  G4double FragMult = std::max((1.0+2.31*(U/CLHEP::MeV - 3.5*A))/100.0, 2.0);
 
  // Parameter Kappa
  fKappa = (1.0 + CLHEP::elm_coupling*(g4calc->A13(FragMult) - 1.0)/(R0*a13));
  fKappa = fKappa*fKappa*fKappa - 1.0;
    
  theTemp->Initialise(A, Z, U, pFreeInternalE0, fKappa, &fClusters);
  
  pMeanTemperature = theTemp->CalcTemperature();
  fChemPotentialNu = theTemp->GetChemicalPotentialNu();
  fChemPotentialMu = theTemp->GetChemicalPotentialMu();
  pMeanMultiplicity = theTemp->GetMeanMultiplicity();
  pMeanEntropy = theTemp->GetEntropy();
}

operator!=()

G4bool G4StatMFMacroCanonical::operator!= ( const G4StatMFMacroCanonical & right ) const

delete

operator=()

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

delete

operator==()

G4bool G4StatMFMacroCanonical::operator== ( const G4StatMFMacroCanonical & right ) const

delete

---

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

• G4StatMFMacroCanonical.hh
• G4StatMFMacroCanonical.cc