G4StatMFMicroCanonical Class Reference

#include <G4StatMFMicroCanonical.hh>

Inheritance diagram for G4StatMFMicroCanonical:

Public Member Functions
	G4StatMFMicroCanonical ()
	~G4StatMFMicroCanonical () override
void	Initialise (const G4Fragment &theFragment) override
G4StatMFChannel *	ChooseAandZ (const G4Fragment &theFragment) override
G4double	Function (G4double T)
	G4StatMFMicroCanonical (const G4StatMFMicroCanonical &right)=delete
G4StatMFMicroCanonical &	operator= (const G4StatMFMicroCanonical &right)=delete
G4bool	operator== (const G4StatMFMicroCanonical &right) const =delete
G4bool	operator!= (const G4StatMFMicroCanonical &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 G4StatMFMicroCanonical.hh.

Constructor & Destructor Documentation

G4StatMFMicroCanonical() [1/2]

G4StatMFMicroCanonical::G4StatMFMicroCanonical

(

)

Definition at line 45 of file G4StatMFMicroCanonical.cc.

{
  fSolver = new G4FunctionSolver<G4StatMFMicroCanonical>(this, 100, 5.e-4);
  fSolver->SetIntervalLimits(t1, t2);
  fPartitionManagerVector.reserve(fMaxMultiplicity);
  g4calc = G4Pow::GetInstance();
}

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

~G4StatMFMicroCanonical()

G4StatMFMicroCanonical::~G4StatMFMicroCanonical ( )

override

Definition at line 54 of file G4StatMFMicroCanonical.cc.

{
  delete fSolver;
  if (!fPartitionManagerVector.empty()) {
    for (auto const & p : fPartitionManagerVector) { delete p; }
  }
}

G4StatMFMicroCanonical() [2/2]

G4StatMFMicroCanonical::G4StatMFMicroCanonical ( const G4StatMFMicroCanonical & right )

delete

Member Function Documentation

ChooseAandZ()

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

overridevirtual

Implements G4VStatMFEnsemble.

Definition at line 150 of file G4StatMFMicroCanonical.cc.

{
  // Choice of fragment atomic numbers and charges 
  // We choose a multiplicity (1,2,3,...) and then a channel
  G4int AA = theFragment.GetA_asInt();
  G4int ZZ = theFragment.GetZ_asInt();
  
  if (G4UniformRand() < fWCompoundNucleus) { 
    
    G4StatMFChannel * aChannel = new G4StatMFChannel;
    aChannel->CreateFragment(AA, ZZ);
    return aChannel;
    
  } else {
    G4double rand = G4UniformRand();
    G4double AccumWeight = fWCompoundNucleus;
    for (auto & ptr : fPartitionManagerVector) {
      AccumWeight += ptr->GetProbability();
      if (rand <= AccumWeight) {
    return ptr->ChooseChannel(A, Z, pMeanTemperature);
      }
    }
  }
  return nullptr;
}

Function()

G4double G4StatMFMicroCanonical::Function ( G4double T )

inline

Definition at line 63 of file G4StatMFMicroCanonical.hh.

  { return (fExEnergy + pFreeInternalE0 - CalcFreeInternalEnergy(T)); }

Initialise()

void G4StatMFMicroCanonical::Initialise ( const G4Fragment & theFragment )

overridevirtual

Implements G4VStatMFEnsemble.

Definition at line 62 of file G4StatMFMicroCanonical.cc.

{
  fPartitionManagerVector.clear();
  // Excitation Energy 
  fExEnergy = theFragment.GetExcitationEnergy();
 
  A = theFragment.GetA_asInt();
  Z = theFragment.GetZ_asInt();
  A13 = g4calc->Z13(A);
  
  fInvLevelDensity = G4StatMFParameters::GetEpsilon0()*(1.0 + 3.0/G4double(A-1));
  
  fSymmetryTerm = G4StatMFParameters::GetGamma0()*(A - 2*Z)*(A - 2*Z)/(G4double)A;
 
  fCoulombTerm = elm_coupling*0.6*Z*Z/(G4StatMFParameters::Getr0()*A13);
 
  // Configuration temperature
  G4double TConf = std::sqrt(8.0*fExEnergy/(G4double)A);
  
  // Free internal energy at Temperature T = 0 (SurfaceTerm at T = 0)
  pFreeInternalE0 = -G4StatMFParameters::GetE0()*A + fSymmetryTerm  
    + G4StatMFParameters::GetBeta0()*A13*A13 + fCoulombTerm;
  
  //G4cout << "Tconf=" <<  TConf << " freeE=" << pFreeInternalE0 << G4endl;
    
  // Mean breakup multiplicity
  pMeanMultiplicity = 0.0;
  
  // Mean channel temperature
  pMeanTemperature = 0.0;
  
  // Mean channel entropy
  pMeanEntropy = 0.0;
  
  // Calculate entropy of compound nucleus
  G4double SCompoundNucleus = CalcEntropyOfCompoundNucleus(TConf);
  
  // Statistical weight of compound nucleus
  fWCompoundNucleus = 1.0; 
  
  // Statistical weight
  G4double W = fWCompoundNucleus;
  // Maximal fragment multiplicity allowed in direct simulation  
 
  for (G4int im = 2; im <= fMaxMultiplicity; ++im) {
    auto ptr = new G4StatMFMicroManager(theFragment, im, pFreeInternalE0, SCompoundNucleus);
    fPartitionManagerVector.push_back(ptr);
    W += ptr->GetProbability();
  }
  
  // Normalization of statistical weights
  for (auto & ptr : fPartitionManagerVector) {
    ptr->Normalize(W);
    pMeanMultiplicity += ptr->GetMeanMultiplicity();
    pMeanTemperature += ptr->GetMeanTemperature();
    pMeanEntropy += ptr->GetMeanEntropy();
  }
 
  fWCompoundNucleus /= W;
  
  pMeanMultiplicity += fWCompoundNucleus;
  pMeanTemperature += TConf * fWCompoundNucleus;
  pMeanEntropy += SCompoundNucleus * fWCompoundNucleus;
}

operator!=()

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

delete

operator=()

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

delete

operator==()

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

delete

---

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

• G4StatMFMicroCanonical.hh
• G4StatMFMicroCanonical.cc