G4EvaporationProbability Class Reference

#include <G4EvaporationProbability.hh>

Inheritance diagram for G4EvaporationProbability:

Public Member Functions
	G4EvaporationProbability (G4int anA, G4int aZ, G4double aGamma)
	~G4EvaporationProbability () override=default
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)
	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)
virtual	~G4VEmissionProbability ()=default
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 (size_t nbin, 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

Detailed Description

Definition at line 41 of file G4EvaporationProbability.hh.

Constructor & Destructor Documentation

G4EvaporationProbability() [1/2]

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

explicit

Definition at line 57 of file G4EvaporationProbability.cc.

  : G4VEmissionProbability(aZ, anA), fGamma(aGamma)
{
  resA13 = lastA = muu = freeU = a0 = delta1 = 0.0;
  pcoeff = fGamma*pEvapMass*CLHEP::millibarn
    /((CLHEP::pi*CLHEP::hbarc)*(CLHEP::pi*CLHEP::hbarc)); 
 
  if(0 == theZ)      { index = 0; }
  else if(1 == theZ) { index = theA; }
  else               { index = theA + 1; }
  if(0 == aZ) {
    ResetIntegrator(30, 0.15*CLHEP::MeV, 0.02);
  } else {
    ResetIntegrator(30, 0.25*CLHEP::MeV, 0.03);
  }
}

~G4EvaporationProbability()

G4EvaporationProbability::~G4EvaporationProbability ( )

overridedefault

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 75 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 80 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 142 of file G4EvaporationProbability.cc.

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

Referenced by G4EvaporationChannel::ComputeProbability().

CrossSection()

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

Definition at line 162 of file G4EvaporationProbability.cc.

{
  // compute power once
  if(resA != lastA) {
    lastA = resA;
    if(0 < index) 
      muu = G4KalbachCrossSection::ComputePowerParameter(resA, index);
  }
  G4double res = 0.0;
  if(OPTxs <= 2) { 
    res = G4ChatterjeeCrossSection::ComputeCrossSection(K, CB, resA13, muu,
                            index, theZ, resA); 
  } else {
    // added barrier penetration factor
    G4double elim = 0.6*CB;
    if (K > elim) {
      res = G4KalbachCrossSection::ComputeCrossSection(K, CB, resA13, muu,
                                   index, theZ, theA, resA);
      //res *= (1.0 - elim/K);
    }
  }
  return res;
}

Referenced by G4EvaporationChannel::ComputeInverseXSection(), and 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

TotalProbability()

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

virtual

Definition at line 85 of file G4EvaporationProbability.cc.

{
  G4int fragA = fragment.GetA_asInt();
  G4int fragZ = fragment.GetZ_asInt();
  G4double U = fragment.GetExcitationEnergy(); 
  a0 = pNuclearLevelData->GetLevelDensity(fragZ,fragA,U);
  freeU = exEnergy;
  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    
    G4double 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;
}

Referenced by G4EvaporationChannel::GetEmissionProbability().

---

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

• G4EvaporationProbability.hh
• G4EvaporationProbability.cc