G4VEmissionProbability Class Reference

#include <G4VEmissionProbability.hh>

Inheritance diagram for G4VEmissionProbability:

Public Member Functions
	G4VEmissionProbability (G4int Z, G4int A)
virtual	~G4VEmissionProbability ()=default
void	Initialise ()
virtual G4double	EmissionProbability (const G4Fragment &fragment, G4double anEnergy)
virtual G4double	ComputeProbability (G4double anEnergy, G4double CB)
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
void	ResetIntegrator (size_t nbin, G4double de, G4double eps)
G4double	IntegrateProbability (G4double elow, G4double ehigh, G4double CB)

Protected Attributes
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 46 of file G4VEmissionProbability.hh.

Constructor & Destructor Documentation

G4VEmissionProbability() [1/2]

G4VEmissionProbability::G4VEmissionProbability ( G4int  Z, G4int  A  )

explicit

Definition at line 41 of file G4VEmissionProbability.cc.

  : OPTxs(3), pVerbose(1), theZ(Z), theA(A), elimit(CLHEP::MeV)
{
  pNuclearLevelData = G4NuclearLevelData::GetInstance(); 
  pG4pow = G4Pow::GetInstance();
  if(A > 0) { pEvapMass = G4NucleiProperties::GetNuclearMass(theA, theZ); }
}

~G4VEmissionProbability()

virtual G4VEmissionProbability::~G4VEmissionProbability ( )

virtualdefault

G4VEmissionProbability() [2/2]

G4VEmissionProbability::G4VEmissionProbability ( const G4VEmissionProbability &  right )

delete

Member Function Documentation

ComputeProbability()

G4double G4VEmissionProbability::ComputeProbability ( G4double  anEnergy, G4double  CB  )

virtual

Reimplemented in G4GEMProbabilityVI, and G4EvaporationProbability.

Definition at line 69 of file G4VEmissionProbability.cc.

{
  return 0.0;
}

Referenced by IntegrateProbability(), and SampleEnergy().

EmissionProbability()

G4double G4VEmissionProbability::EmissionProbability ( const G4Fragment &  fragment, G4double  anEnergy  )

virtual

Reimplemented in G4FissionProbability, and G4GEMProbability.

Definition at line 64 of file G4VEmissionProbability.cc.

{
  return 0.0;
}

Referenced by G4CompetitiveFission::GetEmissionProbability().

GetA()

G4int G4VEmissionProbability::GetA ( void  ) const

inline

Definition at line 63 of file G4VEmissionProbability.hh.

{ return theA; }

GetProbability()

G4double G4VEmissionProbability::GetProbability ( ) const

inline

Definition at line 75 of file G4VEmissionProbability.hh.

{ return pProbability; };

Referenced by G4EvaporationChannel::EmittedFragment().

GetRecoilExcitation()

G4double G4VEmissionProbability::GetRecoilExcitation ( ) const

inline

Definition at line 71 of file G4VEmissionProbability.hh.

{ return fExcRes; };

GetZ()

G4int G4VEmissionProbability::GetZ ( void  ) const

inline

Definition at line 61 of file G4VEmissionProbability.hh.

{ return theZ; }

Referenced by G4TritonEvaporationProbability::CalcAlphaParam(), G4AlphaEvaporationProbability::CalcAlphaParam(), G4DeuteronEvaporationProbability::CalcAlphaParam(), G4He3EvaporationProbability::CalcAlphaParam(), and G4ProtonEvaporationProbability::CalcAlphaParam().

Initialise()

void G4VEmissionProbability::Initialise

(

)

Definition at line 49 of file G4VEmissionProbability.cc.

{
  G4DeexPrecoParameters* param = pNuclearLevelData->GetParameters();
  OPTxs = param->GetDeexModelType();
  pVerbose = param->GetVerbose();
  fFD = param->GetDiscreteExcitationFlag();
  pTolerance = param->GetMinExcitation();
}

Referenced by G4EvaporationChannel::Initialise().

IntegrateProbability()

G4double G4VEmissionProbability::IntegrateProbability ( G4double  elow, G4double  ehigh, G4double  CB  )

protected

Definition at line 74 of file G4VEmissionProbability.cc.

{
  pProbability = 0.0;
  if(elow >= ehigh) { return pProbability; }
 
  emin = elow;
  emax = ehigh;
  eCoulomb = cb;
 
  const G4double edeltamin = 0.2*CLHEP::MeV;
  const G4double edeltamax = 2*CLHEP::MeV;
  G4double edelta = std::max(std::min(elimit, edeltamax), edeltamin);
  G4double xbin = (emax - emin)/edelta + 1.0;
  G4int ibin = xbin;
  if(ibin < 4) ibin = 4;
 
  // providing smart binning 
  G4int nbin = ibin*5;
  edelta = (emax - emin)/ibin;
 
  G4double x(emin), y(0.0);
  G4double edelmicro = edelta*0.02;
  probmax = ComputeProbability(x + edelmicro, eCoulomb);
  G4double problast = probmax;
  if(pVerbose > 1) {
    G4cout << "### G4VEmissionProbability::IntegrateProbability: "
       << "probmax=" << probmax << " Emin=" << emin
       << " Emax=" << emax << " QB=" << cb << " nbin=" << nbin 
       << G4endl;
  }
  fE1 = fE2 = fP2 = 0.0;
  G4double emax0 = emax - edelmicro; 
  G4bool endpoint = false;
  for(G4int i=0; i<nbin; ++i) {
    x += edelta;
    if(x >= emax0) { 
      x = emax0;
      endpoint = true;
    }
    y = ComputeProbability(x, eCoulomb);
    if(pVerbose > 2) { 
      G4cout << "    " << i << ".  E= " << x << "  prob= " << y
         << " Edel= " << edelta << G4endl;
    } 
    if(y >= probmax) {
      probmax = y;
    } else if(0.0 == fE1 && 2*y < probmax) {
      fE1 = x;
    }
    
    G4double del = (y + problast)*edelta*0.5;
    pProbability += del;
    // end of the loop
    if(del < accuracy*pProbability || endpoint) { break; }
    problast = y;
 
    // smart step definition
    if(del != pProbability && del > 0.8*pProbability && 
       0.7*edelta > edeltamin) { 
      edelta *= 0.7;
    } else if(del < 0.1*pProbability && 1.5*edelta < edeltamax) { 
      edelta *= 1.5;
    }
  }
  if(fE1 > emin && fE1 < emax) {
    fE2 = std::max(0.5*(fE1 + emax), emax - edelta);
    fP2 = 2*ComputeProbability(fE2, eCoulomb);
  }
 
  if(pVerbose > 1) { 
    G4cout << " Probability= " << pProbability << " probmax= " 
           << probmax << " emin=" << emin << " emax=" << emax 
       << " E1=" << fE1 << " E2=" << fE2 << G4endl; 
  }
  return pProbability;
}

Referenced by G4GEMProbabilityVI::ComputeTotalProbability(), and G4EvaporationProbability::TotalProbability().

operator!=()

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

delete

operator=()

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

delete

operator==()

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

delete

ResetIntegrator()

void G4VEmissionProbability::ResetIntegrator ( size_t  nbin, G4double  de, G4double  eps  )

protected

Definition at line 58 of file G4VEmissionProbability.cc.

{
  if(de > 0.0)  { elimit = de; }
  if(eps > 0.0) { accuracy = eps; }
}

Referenced by G4EvaporationProbability::G4EvaporationProbability(), and G4GEMProbabilityVI::G4GEMProbabilityVI().

ResetProbability()

void G4VEmissionProbability::ResetProbability ( )

inline

Definition at line 77 of file G4VEmissionProbability.hh.

{ pProbability = 0.0; };

Referenced by G4EvaporationChannel::GetEmissionProbability(), and G4GEMChannelVI::GetEmissionProbability().

SampleEnergy()

G4double G4VEmissionProbability::SampleEnergy

(

)

Definition at line 153 of file G4VEmissionProbability.cc.

{
  static const G4double fact = 1.05;
  static const G4double alim = 0.05;
  static const G4double blim = 20.;
  probmax *= fact;
 
  // two regions with flat and exponential majorant 
  G4double del = emax - emin;
  G4double p1 = 1.0;
  G4double p2 = 0.0;
  G4double a0 = 0.0;
  G4double a1 = 1.0;
  G4double x;
  if(fE1 > 0.0 && fP2 > 0.0 && fP2 < 0.5*probmax) {
    a0 = G4Log(probmax/fP2)/(fE2 - fE1);
    del= fE1 - emin;
    p1 = del;
    x = a0*(emax - fE1);
    if(x < blim) {
      a1 = (x > alim) ? 1.0 - G4Exp(-x) : x*(1.0 - 0.5*x);
    }
    p2 = a1/a0;
    p1 /= (p1 + p2);
    p2 = 1.0 - p1;
  }  
 
  if(pVerbose > 1) {
    G4cout << "### G4VEmissionProbability::SampleEnergy: " 
       << " Emin= " << emin << " Emax= " << emax 
           << "/n    E1=" << fE1 << " p1=" << p1 
       << " probmax=" << probmax << " P2=" << fP2 << G4endl;
  }
 
  CLHEP::HepRandomEngine* rndm = G4Random::getTheEngine();
  const G4int nmax = 1000;
  G4double ekin, g, gmax;
  G4int n = 0;
  do {
    ++n;
    G4double q = rndm->flat();
    if(q <= p1) {
      gmax = probmax;
      ekin = del*q/p1 + emin;
    } else {
      ekin = fE1 - G4Log(1.0 - (q - p1)*a1/p2)/a0;
      x = a0*(ekin - fE1);
      gmax = fP2;
      if(x < blim) {
    gmax = probmax*((x > alim) ? G4Exp(-x) : 1.0 - x*(1.0 - 0.5*x));
      }
    }
    g = ComputeProbability(ekin, eCoulomb);
    if(pVerbose > 2) {
      G4cout << "    " << n
         << ". prob= " << g << " probmax= " << probmax
         << " Ekin= " << ekin << G4endl;
    }
    if((g > gmax || n > nmax) && pVerbose > 1) {
      G4cout << "### G4VEmissionProbability::SampleEnergy for Z= " << theZ 
             << " A= " << theA << " Eex(MeV)=" << fExc << " p1=" << p1
             << "\n    Warning n= " << n
         << " prob/gmax=" << g/gmax 
         << " prob=" << g << " gmax=" << gmax << " probmax=" << probmax 
         << "\n    Ekin= " << ekin << " Emin= " << emin
         << " Emax= " << emax << G4endl;
    }
  } while(gmax*rndm->flat() > g && n < nmax);
  G4double enew = FindRecoilExcitation(ekin);
  if(pVerbose > 1) {
    G4cout << "### SampleEnergy: Efinal= " 
       << enew << " E=" << ekin << "  Eexc=" << fExcRes << G4endl;
  }
  return enew;
}

Referenced by G4EvaporationChannel::EmittedFragment(), and G4GEMProbabilityVI::SampleEvaporationFragment().

SetDecayKinematics()

void G4VEmissionProbability::SetDecayKinematics ( G4int  rZ, G4int  rA, G4double  rmass, G4double  fmass  )

inline

Definition at line 133 of file G4VEmissionProbability.hh.

{
  resZ = rZ;
  resA = rA;
  pMass = fmass;
  pResMass = rmass;
}

Referenced by G4EvaporationChannel::GetEmissionProbability(), and G4GEMChannelVI::GetEmissionProbability().

SetEvapExcitation()

void G4VEmissionProbability::SetEvapExcitation ( G4double  exc )

inline

Definition at line 73 of file G4VEmissionProbability.hh.

{ fExc = exc; };

Member Data Documentation

OPTxs

G4int G4VEmissionProbability::OPTxs

protected

Definition at line 98 of file G4VEmissionProbability.hh.

Referenced by G4EvaporationProbability::CrossSection(), Initialise(), and G4EvaporationProbability::TotalProbability().

pEvapMass

G4double G4VEmissionProbability::pEvapMass = 0.0

protected

Definition at line 106 of file G4VEmissionProbability.hh.

Referenced by G4GEMProbabilityVI::ComputeProbability(), G4EvaporationProbability::ComputeProbability(), G4GEMProbabilityVI::ComputeTotalProbability(), G4EvaporationProbability::G4EvaporationProbability(), G4GEMProbabilityVI::G4GEMProbabilityVI(), G4VEmissionProbability(), G4GEMProbabilityVI::SampleEvaporationFragment(), and G4EvaporationProbability::TotalProbability().

pG4pow

G4Pow* G4VEmissionProbability::pG4pow

protected

Definition at line 96 of file G4VEmissionProbability.hh.

Referenced by G4NeutronEvaporationProbability::CalcAlphaParam(), G4NeutronEvaporationProbability::CalcBetaParam(), G4GEMProbabilityVI::ComputeTotalProbability(), G4GEMProbabilityVI::G4GEMProbabilityVI(), G4VEmissionProbability(), and G4EvaporationProbability::TotalProbability().

pMass

G4double G4VEmissionProbability::pMass = 0.0

protected

Definition at line 105 of file G4VEmissionProbability.hh.

Referenced by G4GEMProbabilityVI::ComputeProbability(), G4EvaporationProbability::ComputeProbability(), G4GEMProbabilityVI::ComputeTotalProbability(), and SetDecayKinematics().

pNuclearLevelData

G4NuclearLevelData* G4VEmissionProbability::pNuclearLevelData

protected

Definition at line 95 of file G4VEmissionProbability.hh.

Referenced by G4GEMProbabilityVI::ComputeProbability(), G4EvaporationProbability::ComputeProbability(), G4GEMProbabilityVI::ComputeTotalProbability(), G4VEmissionProbability(), Initialise(), and G4EvaporationProbability::TotalProbability().

pProbability

G4double G4VEmissionProbability::pProbability = 0.0

protected

Definition at line 108 of file G4VEmissionProbability.hh.

Referenced by G4GEMProbabilityVI::ComputeTotalProbability(), GetProbability(), IntegrateProbability(), ResetProbability(), and G4EvaporationProbability::TotalProbability().

pResMass

G4double G4VEmissionProbability::pResMass = 0.0

protected

Definition at line 107 of file G4VEmissionProbability.hh.

Referenced by G4GEMProbabilityVI::ComputeProbability(), G4EvaporationProbability::ComputeProbability(), G4GEMProbabilityVI::ComputeTotalProbability(), and SetDecayKinematics().

pTolerance

G4double G4VEmissionProbability::pTolerance = 0.0

protected

Definition at line 109 of file G4VEmissionProbability.hh.

Referenced by Initialise().

pVerbose

G4int G4VEmissionProbability::pVerbose

protected

Definition at line 99 of file G4VEmissionProbability.hh.

Referenced by Initialise(), IntegrateProbability(), and SampleEnergy().

resA

G4int G4VEmissionProbability::resA = 0

protected

Definition at line 103 of file G4VEmissionProbability.hh.

Referenced by G4GEMProbabilityVI::ComputeProbability(), G4EvaporationProbability::ComputeProbability(), G4GEMProbabilityVI::ComputeTotalProbability(), G4EvaporationProbability::CrossSection(), SetDecayKinematics(), and G4EvaporationProbability::TotalProbability().

resZ

G4int G4VEmissionProbability::resZ = 0

protected

Definition at line 102 of file G4VEmissionProbability.hh.

Referenced by G4GEMProbabilityVI::ComputeProbability(), G4EvaporationProbability::ComputeProbability(), G4GEMProbabilityVI::ComputeTotalProbability(), SetDecayKinematics(), and G4EvaporationProbability::TotalProbability().

theA

G4int G4VEmissionProbability::theA

protected

Definition at line 101 of file G4VEmissionProbability.hh.

Referenced by G4GEMProbability::CalcAlphaParam(), G4GEMProbability::CalcBetaParam(), G4GEMProbabilityVI::ComputeTotalProbability(), G4EvaporationProbability::CrossSection(), G4GEMProbability::Dump(), G4EvaporationProbability::G4EvaporationProbability(), G4GEMProbabilityVI::G4GEMProbabilityVI(), G4VEmissionProbability(), GetA(), G4GEMProbability::GetCoulombBarrier(), SampleEnergy(), and G4GEMProbabilityVI::SampleEvaporationFragment().

theZ

G4int G4VEmissionProbability::theZ

protected

Definition at line 100 of file G4VEmissionProbability.hh.

Referenced by G4GEMProbability::CalcAlphaParam(), G4GEMProbability::CalcBetaParam(), G4GEMProbabilityVI::ComputeTotalProbability(), G4EvaporationProbability::CrossSection(), G4GEMProbability::Dump(), G4EvaporationProbability::G4EvaporationProbability(), G4VEmissionProbability(), G4GEMProbability::GetCoulombBarrier(), GetZ(), SampleEnergy(), and G4GEMProbabilityVI::SampleEvaporationFragment().

---

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

• G4VEmissionProbability.hh
• G4VEmissionProbability.cc