#include <G4VEmissionProbability.hh>

Inheritance diagram for G4VEmissionProbability:

Public Member Functions
	G4VEmissionProbability (G4int Z, G4int A)

virtual	~G4VEmissionProbability ()=default

virtual 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

G4double	pWidth = 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.

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

Referenced by G4EvaporationProbability::G4EvaporationProbability(), G4FissionProbability::G4FissionProbability(), G4GEMProbability::G4GEMProbability(), G4GEMProbabilityVI::G4GEMProbabilityVI(), G4VEmissionProbability(), operator!=(), operator=(), and operator==().

◆ ~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 G4EvaporationProbability, and G4GEMProbabilityVI.

Definition at line 71 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 66 of file G4VEmissionProbability.cc.

{
  return 0.0;
}

◆ GetA()

G4int G4VEmissionProbability::GetA ( void ) const

inline

Definition at line 63 of file G4VEmissionProbability.hh.

63{ return theA; }

◆ GetProbability()

G4double G4VEmissionProbability::GetProbability ( ) const

inline

Definition at line 80 of file G4VEmissionProbability.hh.

80{ return pProbability; };

G4VEmissionProbability::pProbability

G4double pProbability

Definition G4VEmissionProbability.hh:113

◆ GetRecoilExcitation()

G4double G4VEmissionProbability::GetRecoilExcitation ( ) const

inline

Definition at line 76 of file G4VEmissionProbability.hh.

76{ return fExcRes; };

◆ GetZ()

G4int G4VEmissionProbability::GetZ ( void ) const

inline

Definition at line 61 of file G4VEmissionProbability.hh.

61{ return theZ; }

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

◆ Initialise()

void G4VEmissionProbability::Initialise ( )

virtual

Definition at line 51 of file G4VEmissionProbability.cc.

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

◆ IntegrateProbability()

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

protected

Definition at line 76 of file G4VEmissionProbability.cc.

{
  pProbability = 0.0;
  if(elow >= ehigh) { return pProbability; }
 
  emin = elow;
  emax = ehigh;
  eCoulomb = cb;
 
  const G4double edeltamin = 0.1*CLHEP::MeV;
  const G4double edeltamax = 2*CLHEP::MeV;
  G4double edelta = std::min(std::min(elimit, edeltamax), edeltamin);
  G4double xbin = (emax - emin)/edelta + 1.0;
  G4int ibin = std::max((G4int)xbin, 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 G4EvaporationProbability::TotalProbability(), and G4GEMProbabilityVI::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 60 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 82 of file G4VEmissionProbability.hh.

82{ pProbability = 0.0; };

◆ SampleEnergy()

G4double G4VEmissionProbability::SampleEnergy ( )

Definition at line 154 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, gg, gmax;
  G4int n = 0;
  do {
    ++n;
    G4double q = rndm->flat();
    if (p2 == 0.0) {
      gmax = probmax;
      ekin = del*q + emin;
    } else 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));
      }
    }
    gg = ComputeProbability(ekin, eCoulomb);
    if(pVerbose > 2) {
      G4cout << "    " << n
         << ". prob= " << gg << " probmax= " << probmax
         << " Ekin= " << ekin << G4endl;
    }
    if((gg > gmax || n > nmax) && pVerbose > 1) {
      G4cout << "### G4VEmissionProbability::SampleEnergy for Z= " << theZ 
             << " A= " << theA << " Eex(MeV)=" << fExc << " p1=" << p1
             << "\n    Warning n= " << n
         << " prob/gmax=" << gg/gmax 
         << " prob=" << gg << " gmax=" << gmax << " probmax=" << probmax 
         << "\n    Ekin= " << ekin << " Emin= " << emin
         << " Emax= " << emax << G4endl;
    }
  } while(gmax*rndm->flat() > gg && n < nmax);
  G4double enew = FindRecoilExcitation(ekin);
  if(pVerbose > 1) {
    G4cout << "### SampleEnergy: Efinal= " 
       << enew << " E=" << ekin << "  Eexc=" << fExcRes << G4endl;
  }
  return enew;
}

◆ SetDecayKinematics()

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

inline

Definition at line 68 of file G4VEmissionProbability.hh.

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

◆ SetEvapExcitation()

void G4VEmissionProbability::SetEvapExcitation ( G4double exc )

inline

Definition at line 78 of file G4VEmissionProbability.hh.

78{ fExc = exc; };

Member Data Documentation

◆ OPTxs

G4int G4VEmissionProbability::OPTxs

protected

Definition at line 103 of file G4VEmissionProbability.hh.

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

◆ pEvapMass

G4double G4VEmissionProbability::pEvapMass = 0.0

protected

Definition at line 111 of file G4VEmissionProbability.hh.

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

◆ pG4pow

G4Pow* G4VEmissionProbability::pG4pow

protected

Definition at line 101 of file G4VEmissionProbability.hh.

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

◆ pMass

G4double G4VEmissionProbability::pMass = 0.0

protected

Definition at line 110 of file G4VEmissionProbability.hh.

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

◆ pNuclearLevelData

G4NuclearLevelData* G4VEmissionProbability::pNuclearLevelData

protected

Definition at line 100 of file G4VEmissionProbability.hh.

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

◆ pProbability

G4double G4VEmissionProbability::pProbability = 0.0

protected

Definition at line 113 of file G4VEmissionProbability.hh.

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

◆ pResMass

G4double G4VEmissionProbability::pResMass = 0.0

protected

Definition at line 112 of file G4VEmissionProbability.hh.

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

◆ pTolerance

G4double G4VEmissionProbability::pTolerance = 0.0

protected

Definition at line 114 of file G4VEmissionProbability.hh.

Referenced by Initialise().

◆ pVerbose

G4int G4VEmissionProbability::pVerbose

protected

Definition at line 104 of file G4VEmissionProbability.hh.

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

◆ pWidth

G4double G4VEmissionProbability::pWidth = 0.0

protected

Definition at line 115 of file G4VEmissionProbability.hh.

Referenced by Initialise().

◆ resA

G4int G4VEmissionProbability::resA = 0

protected

Definition at line 108 of file G4VEmissionProbability.hh.

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

◆ resZ

G4int G4VEmissionProbability::resZ = 0

protected

Definition at line 107 of file G4VEmissionProbability.hh.

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

◆ theA

G4int G4VEmissionProbability::theA

protected

Definition at line 106 of file G4VEmissionProbability.hh.

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

◆ theZ

G4int G4VEmissionProbability::theZ

protected

Definition at line 105 of file G4VEmissionProbability.hh.

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

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

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4VEmissionProbability() [1/2]

◆ ~G4VEmissionProbability()

◆ G4VEmissionProbability() [2/2]

Member Function Documentation

◆ ComputeProbability()

◆ EmissionProbability()

◆ GetA()

◆ GetProbability()

◆ GetRecoilExcitation()

◆ GetZ()

◆ Initialise()

◆ IntegrateProbability()

◆ operator!=()

◆ operator=()

◆ operator==()

◆ ResetIntegrator()

◆ ResetProbability()

◆ SampleEnergy()

◆ SetDecayKinematics()

◆ SetEvapExcitation()

Member Data Documentation

◆ OPTxs

◆ pEvapMass

◆ pG4pow

◆ pMass

◆ pNuclearLevelData

◆ pProbability

◆ pResMass

◆ pTolerance

◆ pVerbose

◆ pWidth

◆ resA

◆ resZ

◆ theA

◆ theZ