G4EvaporationProbability Class Reference

#include <G4EvaporationProbability.hh>

Inheritance diagram for G4EvaporationProbability:

Public Member Functions
	G4EvaporationProbability (G4int anA, G4int aZ, G4double aGamma)
	~G4EvaporationProbability () override
virtual G4double	TotalProbability (const G4Fragment &fragment, G4double minKinEnergy, G4double maxKinEnergy, G4double CB, G4double exEnergy)
G4double	ComputeProbability (G4double K, G4double CB) override
G4double	CrossSection (G4double K, G4double CB)
G4double	RecentXS () const
	G4EvaporationProbability (const G4EvaporationProbability &right)=delete
const G4EvaporationProbability &	operator= (const G4EvaporationProbability &right)=delete
G4bool	operator== (const G4EvaporationProbability &right) const =delete
G4bool	operator!= (const G4EvaporationProbability &right) const =delete
Public Member Functions inherited from G4VEmissionProbability
	G4VEmissionProbability (G4int Z, G4int A)
	~G4VEmissionProbability () override=default
G4double	ProbabilityDensityFunction (G4double energy) override
virtual void	Initialise ()
virtual G4double	EmissionProbability (const G4Fragment &fragment, G4double anEnergy)
G4int	GetZ (void) const
G4int	GetA (void) const
void	SetDecayKinematics (G4int rZ, G4int rA, G4double rmass, G4double fmass)
G4double	GetRecoilExcitation () const
void	SetEvapExcitation (G4double exc)
G4double	GetProbability () const
void	ResetProbability ()
G4double	SampleEnergy ()
	G4VEmissionProbability (const G4VEmissionProbability &right)=delete
const G4VEmissionProbability &	operator= (const G4VEmissionProbability &right)=delete
G4bool	operator== (const G4VEmissionProbability &right) const =delete
G4bool	operator!= (const G4VEmissionProbability &right) const =delete

Protected Member Functions
virtual G4double	CalcAlphaParam (const G4Fragment &fragment)
virtual G4double	CalcBetaParam (const G4Fragment &fragment)
Protected Member Functions inherited from G4VEmissionProbability
void	ResetIntegrator (G4double de, G4double eps)
G4double	IntegrateProbability (G4double elow, G4double ehigh, G4double CB)

Additional Inherited Members
Protected Attributes inherited from G4VEmissionProbability
G4NuclearLevelData *	pNuclearLevelData
G4Pow *	pG4pow
G4int	OPTxs
G4int	pVerbose
G4int	theZ
G4int	theA
G4int	resZ = 0
G4int	resA = 0
G4double	pMass = 0.0
G4double	pEvapMass = 0.0
G4double	pResMass = 0.0
G4double	pProbability = 0.0
G4double	pTolerance = 0.0
G4double	pWidth = 0.0

Detailed Description

Definition at line 42 of file G4EvaporationProbability.hh.

Constructor & Destructor Documentation

G4EvaporationProbability() [1/2]

G4EvaporationProbability::G4EvaporationProbability ( G4int anA, G4int aZ, G4double aGamma )

explicit

Definition at line 71 of file G4EvaporationProbability.cc.

  : G4VEmissionProbability(aZ, anA), fGamma(aGamma)
{
  resA13 = lastA = muu = freeU = a0 = a1 = delta0 = delta1 = 0.0;
  pcoeff = fGamma*pEvapMass*CLHEP::millibarn
    /((CLHEP::pi*CLHEP::hbarc)*(CLHEP::pi*CLHEP::hbarc)); 
 
  if (1 == theZ && 1 == theA) { index = 1; }
  else if (1 == theZ && 2 == theA) { index = 2; }
  else if (1 == theZ && 3 == theA) { index = 3; }
  else if (2 == theZ && 3 == theA) { index = 4; }
  else if (2 == theZ && 4 == theA) { index = 5; }
 
  if (OPTxs == 1) {
    const G4ParticleDefinition* part = nullptr;
    if (index == 1) { part = G4Proton::Proton(); }
    else if (index == 2) { part = G4Deuteron::Deuteron(); }
    else if (index == 3) { part = G4Triton::Triton(); }
    else if (index == 4) { part = G4He3::He3(); }
    else if (index == 5) { part = G4Alpha::Alpha(); }
    else { part = G4Neutron::Neutron(); }
    fXSection = new G4InterfaceToXS(part, index);
  }
  
  if (0 == aZ) {
    ResetIntegrator(0.15*CLHEP::MeV, 0.01);
  } else {
    ResetIntegrator(0.20*CLHEP::MeV, 0.01);
  }
}

Referenced by G4AlphaEvaporationProbability::G4AlphaEvaporationProbability(), G4DeuteronEvaporationProbability::G4DeuteronEvaporationProbability(), G4EvaporationProbability(), G4He3EvaporationProbability::G4He3EvaporationProbability(), G4NeutronEvaporationProbability::G4NeutronEvaporationProbability(), G4ProtonEvaporationProbability::G4ProtonEvaporationProbability(), G4TritonEvaporationProbability::G4TritonEvaporationProbability(), operator!=(), operator=(), and operator==().

~G4EvaporationProbability()

G4EvaporationProbability::~G4EvaporationProbability ( )

override

Definition at line 103 of file G4EvaporationProbability.cc.

{
  delete fXSection;
}

G4EvaporationProbability() [2/2]

G4EvaporationProbability::G4EvaporationProbability ( const G4EvaporationProbability & right )

delete

Member Function Documentation

CalcAlphaParam()

G4double G4EvaporationProbability::CalcAlphaParam ( const G4Fragment & fragment )

protectedvirtual

Reimplemented in G4AlphaEvaporationProbability, G4DeuteronEvaporationProbability, G4He3EvaporationProbability, G4NeutronEvaporationProbability, G4ProtonEvaporationProbability, and G4TritonEvaporationProbability.

Definition at line 108 of file G4EvaporationProbability.cc.

{
  return 1.0;
}

Referenced by TotalProbability().

CalcBetaParam()

G4double G4EvaporationProbability::CalcBetaParam ( const G4Fragment & fragment )

protectedvirtual

Reimplemented in G4AlphaEvaporationProbability, G4DeuteronEvaporationProbability, G4He3EvaporationProbability, G4NeutronEvaporationProbability, G4ProtonEvaporationProbability, and G4TritonEvaporationProbability.

Definition at line 113 of file G4EvaporationProbability.cc.

{
  return 1.0;
}

Referenced by TotalProbability().

ComputeProbability()

G4double G4EvaporationProbability::ComputeProbability ( G4double K, G4double CB )

overridevirtual

Reimplemented from G4VEmissionProbability.

Definition at line 175 of file G4EvaporationProbability.cc.

{
  G4double K = std::max(kinE, kmin);
  // abnormal case - should never happens
  if(pMass < pEvapMass + pResMass + K) { return 0.0; }
    
  G4double K1 = pMass - pEvapMass - K;
  G4double mres = std::sqrt(K1*K1 - K*(2*pEvapMass + K));
 
  G4double excRes = mres - pResMass;
  if (excRes < 0.0) { return 0.0; }
  G4double K2 = 0.5*(pMass + pEvapMass + mres)*(pMass - pEvapMass - mres)/mres;
  G4double xs = CrossSection(K2, CB);
  if (xs <= 0.0) { return 0.0; }
 
  a1 = pNuclearLevelData->GetLevelDensity(resZ, resA, excRes);
  G4double E0 = std::max(freeU - delta0, 0.0);
  G4double E1 = std::max(excRes - delta1, 0.0);
  G4double prob = pcoeff*G4Exp(2.0*(std::sqrt(a1*E1) - std::sqrt(a0*E0)))*K*xs;
  return prob;
}

CrossSection()

G4double G4EvaporationProbability::CrossSection	(	G4double	K,
		G4double	CB )

Definition at line 198 of file G4EvaporationProbability.cc.

{
  G4double K = std::max(kine, kmin);
  // compute power once
  if (OPTxs > 1 && 0 < index && resA != lastA) {
    lastA = resA;
    muu = G4KalbachCrossSection::ComputePowerParameter(resA, index);
  }
  // In the case of OPTxs = 0 this method is not called
  if (OPTxs == 1) {
    G4int Z = std::min(resZ, ZMAXNUCLEARDATA);
    if (0 == index) {
      G4double e1 = lowEnergyLimitMeV[Z];
      if (e1 == 0.0) { e1 = lim; }
      K = std::max(K, e1);
    } else {
      if (K < 0.5*CB) {
    recentXS = 0.0;
    return recentXS;
      }
      K = std::max(K, 2*CB);
    }
    G4double corr = G4DeexPrecoUtility::CorrectionFactor(index, theZ, resA13, CB, kine);
    recentXS = corr*fXSection->GetElementCrossSection(K, Z)/CLHEP::millibarn;
 
  } else if (OPTxs == 2) { 
    recentXS = G4ChatterjeeCrossSection::ComputeCrossSection(K, CB, resA13, muu,
                                                             index, theZ, resA); 
  } else if (OPTxs == 3) {
    recentXS = G4KalbachCrossSection::ComputeCrossSection(K, CB, resA13, muu,
                                                          index, theZ, theA, resA);
  }
  return recentXS;
}

Referenced by ComputeProbability().

operator!=()

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

delete

operator=()

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

delete

operator==()

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

delete

RecentXS()

G4double G4EvaporationProbability::RecentXS ( ) const

inline

Definition at line 62 of file G4EvaporationProbability.hh.

{ return recentXS; };

TotalProbability()

G4double G4EvaporationProbability::TotalProbability ( const G4Fragment & fragment, G4double minKinEnergy, G4double maxKinEnergy, G4double CB, G4double exEnergy )

virtual

Definition at line 118 of file G4EvaporationProbability.cc.

{
  G4int fragA = fragment.GetA_asInt();
  G4int fragZ = fragment.GetZ_asInt(); 
  freeU = exEnergy;
  a0 = pNuclearLevelData->GetLevelDensity(fragZ, fragA, freeU);
  delta0 = pNuclearLevelData->GetPairingCorrection(fragZ, fragA);
  delta1 = pNuclearLevelData->GetPairingCorrection(resZ, resA);
  resA13 = pG4pow->Z13(resA);
  /*
  G4cout << "G4EvaporationProbability: Z= " << theZ << " A= " << theA 
     << " resZ= " << resZ << " resA= " << resA 
     << " fragZ= " << fragZ << " fragA= " << fragA 
     << "\n   freeU= " << freeU  
     << " a0= " << a0 << " OPT= " << OPTxs << " emin= " 
         << minEnergy << " emax= " << maxEnergy 
     << " CB= " << CB << G4endl;
  */
  if (OPTxs==0) {
 
    G4double SystemEntropy = 2.0*std::sqrt(a0*freeU);
    const G4double RN2 = 2.25*CLHEP::fermi*CLHEP::fermi
      /(CLHEP::twopi*CLHEP::hbar_Planck*hbar_Planck);
 
    G4double Alpha = CalcAlphaParam(fragment);
    G4double Beta = CalcBetaParam(fragment);
 
    // to be checked where to use a0, where - a1    
    a1 = pNuclearLevelData->GetLevelDensity(resZ,resA,freeU);
    G4double GlobalFactor = fGamma*Alpha*pEvapMass*RN2*resA13*resA13/(a1*a1);
    
    G4double maxea = maxEnergy*a1;
    G4double Term1 = Beta*a1 - 1.5 + maxea;
    G4double Term2 = (2.0*Beta*a1-3.0)*std::sqrt(maxea) + 2*maxea;
    
    G4double ExpTerm1 = (SystemEntropy <= explim) ? G4Exp(-SystemEntropy) : 0.0;
    
    G4double ExpTerm2 = 2.*std::sqrt(maxea) - SystemEntropy;
    ExpTerm2 = std::min(ExpTerm2, explim);
    ExpTerm2 = G4Exp(ExpTerm2);
    
    pProbability = GlobalFactor*(Term1*ExpTerm1 + Term2*ExpTerm2);
             
  } else {
    // if Coulomb barrier cutoff is superimposed for all cross sections 
    // then the limit is the Coulomb Barrier
    pProbability = IntegrateProbability(minEnergy, maxEnergy, CB);
  }
  /*
  G4cout << "TotalProbability:  Emin=" << minEnergy << " Emax= " << maxEnergy 
     << " CB= " << CB << " prob=" << pProbability << G4endl;
  */
  return pProbability;
}

---

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

• G4EvaporationProbability.hh
• G4EvaporationProbability.cc