G4StatMFMacroCanonical Class Reference

#include <G4StatMFMacroCanonical.hh>

Inheritance diagram for G4StatMFMacroCanonical:

Data Structures
struct	DeleteFragment

Public Member Functions
G4StatMFChannel *	ChooseAandZ (const G4Fragment &theFragment)
	G4StatMFMacroCanonical (G4Fragment const &theFragment)
G4double	GetMeanMultiplicity (void) const
G4double	GetMeanTemperature (void) const
	~G4StatMFMacroCanonical ()

Protected Attributes
G4double	__FreeInternalE0
G4double	__MeanEntropy
G4double	__MeanMultiplicity
G4double	__MeanTemperature

Private Member Functions
void	CalculateTemperature (const G4Fragment &theFragment)
G4double	ChooseA (G4int A, std::vector< G4int > &ANumbers)
G4StatMFChannel *	ChooseZ (G4int &Z, std::vector< G4int > &FragmentsA)
	G4StatMFMacroCanonical ()
	G4StatMFMacroCanonical (const G4StatMFMacroCanonical &)
void	Initialize (const G4Fragment &theFragment)
G4bool	operator!= (const G4StatMFMacroCanonical &right) const
G4StatMFMacroCanonical &	operator= (const G4StatMFMacroCanonical &right)
G4bool	operator== (const G4StatMFMacroCanonical &right) const

Private Attributes
G4double	_ChemPotentialMu
G4double	_ChemPotentialNu
G4double	_Kappa
std::vector< G4VStatMFMacroCluster * >	_theClusters

Detailed Description

Definition at line 49 of file G4StatMFMacroCanonical.hh.

Constructor & Destructor Documentation

G4StatMFMacroCanonical() [1/3]

G4StatMFMacroCanonical::G4StatMFMacroCanonical

(

G4Fragment const &

theFragment

)

Definition at line 43 of file G4StatMFMacroCanonical.cc.

{
 
  // Get memory for clusters
  _theClusters.push_back(new G4StatMFMacroNucleon);              // Size 1
  _theClusters.push_back(new G4StatMFMacroBiNucleon);            // Size 2
  _theClusters.push_back(new G4StatMFMacroTriNucleon);           // Size 3
  _theClusters.push_back(new G4StatMFMacroTetraNucleon);         // Size 4
  for (G4int i = 4; i < theFragment.GetA_asInt(); i++)   
    _theClusters.push_back(new G4StatMFMacroMultiNucleon(i+1)); // Size 5 ... A
  
  // Perform class initialization
  Initialize(theFragment);
    
}

References _theClusters, G4Fragment::GetA_asInt(), and Initialize().

~G4StatMFMacroCanonical()

G4StatMFMacroCanonical::~G4StatMFMacroCanonical

(

)

Definition at line 60 of file G4StatMFMacroCanonical.cc.

{
  // garbage collection
  if (!_theClusters.empty()) 
    {
      std::for_each(_theClusters.begin(),_theClusters.end(),DeleteFragment());
    }
}

References _theClusters.

G4StatMFMacroCanonical() [2/3]

G4StatMFMacroCanonical::G4StatMFMacroCanonical ( )

inlineprivate

Definition at line 61 of file G4StatMFMacroCanonical.hh.

{};

G4StatMFMacroCanonical() [3/3]

G4StatMFMacroCanonical::G4StatMFMacroCanonical ( const G4StatMFMacroCanonical &  )

inlineprivate

Definition at line 65 of file G4StatMFMacroCanonical.hh.

: G4VStatMFEnsemble() {};

Member Function Documentation

CalculateTemperature()

void G4StatMFMacroCanonical::CalculateTemperature ( const G4Fragment &  theFragment )

private

Definition at line 89 of file G4StatMFMacroCanonical.cc.

{
  // Excitation Energy
  G4double U = theFragment.GetExcitationEnergy();
  
  G4int A = theFragment.GetA_asInt();
  G4int Z = theFragment.GetZ_asInt();
  
  // Fragment Multiplicity
  G4double FragMult = std::max((1.0+(2.31/MeV)*(U/A - 3.5*MeV))*A/100.0, 2.0);
 
  // Parameter Kappa
  G4Pow* g4calc = G4Pow::GetInstance();
  _Kappa = (1.0+elm_coupling*(g4calc->A13(FragMult)-1)/
        (G4StatMFParameters::Getr0()*g4calc->Z13(A)));
  _Kappa = _Kappa*_Kappa*_Kappa - 1.0;
    
  G4StatMFMacroTemperature * theTemp = new  
    G4StatMFMacroTemperature(A,Z,U,__FreeInternalE0,_Kappa,&_theClusters);
  
  __MeanTemperature = theTemp->CalcTemperature();
  _ChemPotentialNu = theTemp->GetChemicalPotentialNu();
  _ChemPotentialMu = theTemp->GetChemicalPotentialMu();
  __MeanMultiplicity = theTemp->GetMeanMultiplicity();
  __MeanEntropy = theTemp->GetEntropy();
    
  delete theTemp;           
  
  return;
}

References G4VStatMFEnsemble::__FreeInternalE0, G4VStatMFEnsemble::__MeanEntropy, G4VStatMFEnsemble::__MeanMultiplicity, G4VStatMFEnsemble::__MeanTemperature, _ChemPotentialMu, _ChemPotentialNu, _Kappa, _theClusters, A, G4Pow::A13(), G4StatMFMacroTemperature::CalcTemperature(), source.hepunit::elm_coupling, G4Fragment::GetA_asInt(), G4StatMFMacroTemperature::GetChemicalPotentialMu(), G4StatMFMacroTemperature::GetChemicalPotentialNu(), G4StatMFMacroTemperature::GetEntropy(), G4Fragment::GetExcitationEnergy(), G4Pow::GetInstance(), G4StatMFMacroTemperature::GetMeanMultiplicity(), G4StatMFParameters::Getr0(), G4Fragment::GetZ_asInt(), G4INCL::Math::max(), MeV, Z, and G4Pow::Z13().

Referenced by Initialize().

ChooseA()

G4double G4StatMFMacroCanonical::ChooseA ( G4int  A, std::vector< G4int > &  ANumbers  )

private

Definition at line 164 of file G4StatMFMacroCanonical.cc.

{
  G4double multiplicity = 0.0;
  G4int i;
    
  std::vector<G4double> AcumMultiplicity;
  AcumMultiplicity.reserve(A);
 
  AcumMultiplicity.push_back((*(_theClusters.begin()))->GetMeanMultiplicity());
  for (std::vector<G4VStatMFMacroCluster*>::iterator it = _theClusters.begin()+1;
       it != _theClusters.end(); ++it)
    {
      AcumMultiplicity.push_back((*it)->GetMeanMultiplicity()+AcumMultiplicity.back());
    }
  
  G4int CheckA;
  do {
    CheckA = -1;
    G4int SumA = 0;
    G4int ThisOne = 0;
    multiplicity = 0.0;
    for (i = 0; i < A; i++) ANumbers[i] = 0;
    do {
      G4double RandNumber = G4UniformRand()*__MeanMultiplicity;
      for (i = 0; i < A; i++) {
    if (RandNumber < AcumMultiplicity[i]) {
      ThisOne = i;
      break;
    }
      }
      multiplicity++;
      ANumbers[ThisOne] = ANumbers[ThisOne]+1;
      SumA += ThisOne+1;
      CheckA = A - SumA;
      
      // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
    } while (CheckA > 0);
    
    // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
  } while (CheckA < 0 || std::abs(__MeanMultiplicity - multiplicity) > std::sqrt(__MeanMultiplicity) + 0.5);
  
  return multiplicity;
}

References G4VStatMFEnsemble::__MeanMultiplicity, _theClusters, A, and G4UniformRand.

Referenced by ChooseAandZ().

ChooseAandZ()

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

virtual

Implements G4VStatMFEnsemble.

Definition at line 122 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);
}

References A, ChooseA(), ChooseZ(), G4Fragment::GetA_asInt(), G4Fragment::GetZ_asInt(), and Z.

Referenced by G4StatMF::BreakItUp().

ChooseZ()

G4StatMFChannel * G4StatMFMacroCanonical::ChooseZ ( G4int &  Z, std::vector< G4int > &  FragmentsA  )

private

Definition at line 209 of file G4StatMFMacroCanonical.cc.

{
  G4Pow* g4calc = G4Pow::GetInstance();
  std::vector<G4int> FragmentsZ;
  
  G4int DeltaZ = 0;
  G4double CP =  G4StatMFParameters::GetCoulomb();
  G4int multiplicity = FragmentsA.size();
  
  do {
    FragmentsZ.clear();
    G4int SumZ = 0;
    for (G4int i = 0; i < multiplicity; i++) 
      {
    G4int A = FragmentsA[i];
    if (A <= 1) 
      {
        G4double RandNumber = G4UniformRand();
        if (RandNumber < (*_theClusters.begin())->GetZARatio()) 
          {
        FragmentsZ.push_back(1);
        SumZ += FragmentsZ[i];
          } 
        else FragmentsZ.push_back(0);
      } 
    else 
      {
        G4double RandZ;
        G4double CC = 8.0*G4StatMFParameters::GetGamma0()
          + 2*CP*g4calc->Z23(FragmentsA[i]);
        G4double ZMean;
        if (FragmentsA[i] > 1 && FragmentsA[i] < 5) { ZMean = 0.5*FragmentsA[i]; }
        else {
          ZMean = FragmentsA[i]*(4.0*G4StatMFParameters::GetGamma0()
                     + _ChemPotentialNu)/CC; 
        }
        G4double ZDispersion = std::sqrt(FragmentsA[i]*__MeanTemperature/CC);
        G4int z;
        do 
          {
        RandZ = G4RandGauss::shoot(ZMean,ZDispersion);
        z = G4lrint(RandZ+0.5);
        // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
          } while (z < 0 || z > A);
        FragmentsZ.push_back(z);
        SumZ += z;
      }
      }
    DeltaZ = Z - SumZ;
  // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
  } while (std::abs(DeltaZ) > 1);
    
  // DeltaZ can be 0, 1 or -1
  G4int idx = 0;
  if (DeltaZ < 0.0)
    {
      while (FragmentsZ[idx] < 1) { ++idx; }
    }
  FragmentsZ[idx] += DeltaZ;
  
  G4StatMFChannel * theChannel = new G4StatMFChannel;
  for (G4int i = multiplicity-1; i >= 0; i--) 
    {
      theChannel->CreateFragment(FragmentsA[i],FragmentsZ[i]);
    }
 
  return theChannel;
}

References G4VStatMFEnsemble::__MeanTemperature, _ChemPotentialNu, _theClusters, A, anonymous_namespace{paraMaker.cc}::CP, G4StatMFChannel::CreateFragment(), G4lrint(), G4UniformRand, G4StatMFParameters::GetCoulomb(), G4StatMFParameters::GetGamma0(), G4Pow::GetInstance(), G4INCL::DeJongSpin::shoot(), Z, and G4Pow::Z23().

Referenced by ChooseAandZ().

GetMeanMultiplicity()

G4double G4VStatMFEnsemble::GetMeanMultiplicity ( void  ) const

inlineinherited

Definition at line 68 of file G4VStatMFEnsemble.hh.

{return __MeanMultiplicity;}

References G4VStatMFEnsemble::__MeanMultiplicity.

Referenced by G4StatMF::BreakItUp().

GetMeanTemperature()

G4double G4VStatMFEnsemble::GetMeanTemperature ( void  ) const

inlineinherited

Definition at line 70 of file G4VStatMFEnsemble.hh.

{return __MeanTemperature;}

References G4VStatMFEnsemble::__MeanTemperature.

Referenced by G4StatMF::BreakItUp().

Initialize()

void G4StatMFMacroCanonical::Initialize ( const G4Fragment &  theFragment )

private

Definition at line 70 of file G4StatMFMacroCanonical.cc.

{
  
  G4int A = theFragment.GetA_asInt();
  G4int Z = theFragment.GetZ_asInt();
  G4double x = 1.0 - 2.0*Z/G4double(A);
  G4Pow* g4calc = G4Pow::GetInstance();
  
  // Free Internal energy at T = 0
  __FreeInternalE0 = A*( -G4StatMFParameters::GetE0() +     // Volume term (for T = 0)
             G4StatMFParameters::GetGamma0()*x*x) // Symmetry term
    + G4StatMFParameters::GetBeta0()*g4calc->Z23(A) +   // Surface term (for T = 0)
    0.6*elm_coupling*Z*Z/(G4StatMFParameters::Getr0()*     // Coulomb term 
              g4calc->Z13(A));
  
  CalculateTemperature(theFragment);
  return;
}

References G4VStatMFEnsemble::__FreeInternalE0, A, CalculateTemperature(), source.hepunit::elm_coupling, G4Fragment::GetA_asInt(), G4StatMFParameters::GetBeta0(), G4StatMFParameters::GetE0(), G4StatMFParameters::GetGamma0(), G4Pow::GetInstance(), G4StatMFParameters::Getr0(), G4Fragment::GetZ_asInt(), Z, G4Pow::Z13(), and G4Pow::Z23().

Referenced by G4StatMFMacroCanonical().

operator!=()

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

private

operator=()

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

private

operator==()

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

private

Field Documentation

__FreeInternalE0

G4double G4VStatMFEnsemble::__FreeInternalE0

protectedinherited

Definition at line 75 of file G4VStatMFEnsemble.hh.

Referenced by G4StatMFMicroCanonical::CalcEntropyOfCompoundNucleus(), CalculateTemperature(), Initialize(), and G4StatMFMicroCanonical::Initialize().

__MeanEntropy

G4double G4VStatMFEnsemble::__MeanEntropy

protectedinherited

Definition at line 81 of file G4VStatMFEnsemble.hh.

Referenced by CalculateTemperature(), and G4StatMFMicroCanonical::Initialize().

__MeanMultiplicity

G4double G4VStatMFEnsemble::__MeanMultiplicity

protectedinherited

Definition at line 84 of file G4VStatMFEnsemble.hh.

Referenced by CalculateTemperature(), ChooseA(), G4VStatMFEnsemble::GetMeanMultiplicity(), and G4StatMFMicroCanonical::Initialize().

__MeanTemperature

G4double G4VStatMFEnsemble::__MeanTemperature

protectedinherited

Definition at line 78 of file G4VStatMFEnsemble.hh.

Referenced by CalculateTemperature(), G4StatMFMicroCanonical::ChooseAandZ(), ChooseZ(), G4VStatMFEnsemble::GetMeanTemperature(), and G4StatMFMicroCanonical::Initialize().

_ChemPotentialMu

G4double G4StatMFMacroCanonical::_ChemPotentialMu

private

Definition at line 98 of file G4StatMFMacroCanonical.hh.

Referenced by CalculateTemperature().

_ChemPotentialNu

G4double G4StatMFMacroCanonical::_ChemPotentialNu

private

Definition at line 101 of file G4StatMFMacroCanonical.hh.

Referenced by CalculateTemperature(), and ChooseZ().

_Kappa

G4double G4StatMFMacroCanonical::_Kappa

private

Definition at line 105 of file G4StatMFMacroCanonical.hh.

Referenced by CalculateTemperature().

_theClusters

std::vector<G4VStatMFMacroCluster*> G4StatMFMacroCanonical::_theClusters

private

Definition at line 108 of file G4StatMFMacroCanonical.hh.

Referenced by CalculateTemperature(), ChooseA(), ChooseZ(), G4StatMFMacroCanonical(), and ~G4StatMFMacroCanonical().

---

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

• source/processes/hadronic/models/de_excitation/multifragmentation/include/G4StatMFMacroCanonical.hh
• source/processes/hadronic/models/de_excitation/multifragmentation/src/G4StatMFMacroCanonical.cc