#include <G4VEmissionProbability.hh>

Inheritance diagram for G4VEmissionProbability:

Public Member Functions
	G4VEmissionProbability (G4int Z, G4int A)

virtual	~G4VEmissionProbability ()

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 Z, G4int A, G4double rmass, G4double fmass)

G4double	GetRecoilExcitation () const

void	SetEvapExcitation (G4double exc)

G4double	GetProbability () const

void	ResetProbability ()

G4double	SampleEnergy ()

Protected Member Functions
void	ResetIntegrator (size_t nbin, G4double de, G4double eps)

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

Protected Attributes
G4int	OPTxs

G4int	pVerbose

G4int	theZ

G4int	theA

G4int	resZ

G4int	resA

G4double	pMass

G4double	pEvapMass

G4double	pResMass

G4double	pProbability

G4NuclearLevelData *	pNuclearLevelData

G4Pow *	pG4pow

Detailed Description

Definition at line 47 of file G4VEmissionProbability.hh.

Constructor & Destructor Documentation

◆ G4VEmissionProbability()

G4VEmissionProbability::G4VEmissionProbability	(	G4int	Z,
		G4int	A
	)

explicit

Definition at line 39 of file G4VEmissionProbability.cc.

  : OPTxs(3),pVerbose(1),theZ(Z),theA(A),resZ(0),resA(0),
    pMass(0.0),pEvapMass(0.0),pResMass(0.0),fExc(0.0),fExcRes(0.0),
    elimit(CLHEP::MeV),accuracy(0.02),fFD(false)
{
  pNuclearLevelData = G4NuclearLevelData::GetInstance(); 
  pG4pow = G4Pow::GetInstance();
  if(A > 0) { pEvapMass = G4NucleiProperties::GetNuclearMass(theA, theZ); }
  // G4cout << "G4VEvaporationProbability: Z= " << theZ << " A= " << theA 
  // << " M(GeV)= " << pEvapMass/1000. << G4endl;
  length = nbin = 0;
  emin = emax = eCoulomb = pProbability = probmax = 0.0;
}

◆ ~G4VEmissionProbability()

G4VEmissionProbability::~G4VEmissionProbability ( )

virtual

Definition at line 53 of file G4VEmissionProbability.cc.

54{}

Member Function Documentation

◆ ComputeProbability()

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

virtual

Reimplemented in G4GEMProbabilityVI, and G4EvaporationProbability.

Definition at line 76 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 71 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.

63{ return theA; }

◆ GetProbability()

G4double G4VEmissionProbability::GetProbability ( ) const

inline

Definition at line 75 of file G4VEmissionProbability.hh.

75{ return pProbability; };

◆ GetRecoilExcitation()

G4double G4VEmissionProbability::GetRecoilExcitation ( ) const

inline

Definition at line 71 of file G4VEmissionProbability.hh.

71{ return fExcRes; };

◆ GetZ()

G4int G4VEmissionProbability::GetZ ( void ) const

inline

Definition at line 61 of file G4VEmissionProbability.hh.

61{ return theZ; }

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

◆ Initialise()

void G4VEmissionProbability::Initialise ( )

Definition at line 56 of file G4VEmissionProbability.cc.

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

Referenced by G4EvaporationChannel::Initialise().

◆ IntegrateProbability()

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

protected

Definition at line 81 of file G4VEmissionProbability.cc.

{
  pProbability = 0.0;
  if(elow >= ehigh) { return pProbability; }
 
  emin = elow;
  emax = ehigh;
  eCoulomb = cb;
 
  G4double edelta = elimit;
  nbin = (size_t)((emax - emin)/edelta) + 1;
  const G4double edeltamin = 0.2*CLHEP::MeV;
  const G4double edeltamax = 2*CLHEP::MeV;
  if(nbin < 4) { 
    nbin = 4;
    edelta = (emax - emin)/(G4double)nbin;
  } else if(nbin > length) {
    nbin = length;
  }
 
  G4double x(emin), del, y; 
  G4double edelmicro= edelta*0.02;
  probmax = ComputeProbability(x + edelmicro, eCoulomb);
  G4double problast = probmax;
  if(pVerbose > 2) {
    G4cout << "### G4VEmissionProbability::IntegrateProbability: " 
       << " Emax= " << emax << " QB= " << cb << " nbin= " << nbin 
       << G4endl;
    G4cout << "    0.  E= " << emin << "  prob= " << probmax << G4endl;
  }
  for(size_t i=1; i<=nbin; ++i) {
    x += edelta;
    if(x > emax) { 
      edelta += (emax - x);
      x = emax; 
    }
    G4bool endpoint = (std::abs(x - emax) < edelmicro) ? true : false;
    G4double xx = endpoint ? x - edelmicro : x;
    y = ComputeProbability(xx, eCoulomb);
    if(pVerbose > 2) { 
      G4cout << "    " << i << ".  E= " << x << "  prob= " << y 
         << " Edel= " << edelta << G4endl;
    } 
    probmax = std::max(probmax, y);
    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(pVerbose > 1) { 
    G4cout << " Probability= " << pProbability << " probmax= " 
           << probmax << G4endl; 
  }
  return pProbability;
}

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

◆ ResetIntegrator()

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

protected

Definition at line 64 of file G4VEmissionProbability.cc.

{
  if(nbins > 0) { length = nbins; }
  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.

77{ pProbability = 0.0; };

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

◆ SampleEnergy()

G4double G4VEmissionProbability::SampleEnergy ( )

Definition at line 149 of file G4VEmissionProbability.cc.

{
  static const G4double fact = 1.05;
  probmax *= fact;
 
  if(pVerbose > 1) {
    G4cout << "### G4VEmissionProbability::SampleEnergy: " 
       << " Emin= " << emin << " Emax= " << emax 
       << " probmax= " << probmax << G4endl;
  }
 
  CLHEP::HepRandomEngine* rndm = G4Random::getTheEngine();
  const G4int nmax = 100;
  G4double del = emax - emin;
  G4double ekin, g;
  G4int n = 0;
  do {
    ekin = del*rndm->flat() + emin; 
    ++n;
    g = ComputeProbability(ekin, eCoulomb);
    if(pVerbose > 2) {
      G4cout << "    " << n
         << ". prob= " << g << " probmax= " << probmax
         << " Ekin= " << ekin << G4endl;
    }
    if((g > probmax || n > nmax) && pVerbose > 1) {
      G4cout << "### G4VEmissionProbability::SampleEnergy for Z= " << theZ 
             << " A= " << theA 
             << "\n    Warning n= " << n
         << " prob/probmax= " << g/probmax 
         << " prob= " << g << " probmax= " << probmax 
         << "\n    Ekin= " << ekin << " Emin= " << emin
         << " Emax= " << emax << G4endl;
    }
  } while(probmax*rndm->flat() > g && n < nmax);
  return (fFD) ? FindRecoilExcitation(ekin) : ekin;
}

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

◆ SetDecayKinematics()

void G4VEmissionProbability::SetDecayKinematics	(	G4int	Z,
		G4int	A,
		G4double	rmass,
		G4double	fmass
	)

inline

Definition at line 131 of file G4VEmissionProbability.hh.

{
  resZ = Z;
  resA = A;
  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.

73{ fExc = exc; };

Member Data Documentation

◆ OPTxs

G4int G4VEmissionProbability::OPTxs

protected

Definition at line 89 of file G4VEmissionProbability.hh.

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

◆ pEvapMass

G4double G4VEmissionProbability::pEvapMass

protected

Definition at line 97 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 102 of file G4VEmissionProbability.hh.

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

◆ pMass

G4double G4VEmissionProbability::pMass

protected

Definition at line 96 of file G4VEmissionProbability.hh.

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

◆ pNuclearLevelData

G4NuclearLevelData* G4VEmissionProbability::pNuclearLevelData

protected

Definition at line 101 of file G4VEmissionProbability.hh.

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

◆ pProbability

G4double G4VEmissionProbability::pProbability

protected

Definition at line 99 of file G4VEmissionProbability.hh.

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

◆ pResMass

G4double G4VEmissionProbability::pResMass

protected

Definition at line 98 of file G4VEmissionProbability.hh.

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

◆ pVerbose

G4int G4VEmissionProbability::pVerbose

protected

Definition at line 90 of file G4VEmissionProbability.hh.

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

◆ resA

G4int G4VEmissionProbability::resA

protected

Definition at line 94 of file G4VEmissionProbability.hh.

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

◆ resZ

G4int G4VEmissionProbability::resZ

protected

Definition at line 93 of file G4VEmissionProbability.hh.

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

◆ theA

G4int G4VEmissionProbability::theA

protected

Definition at line 92 of file G4VEmissionProbability.hh.

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

◆ theZ

G4int G4VEmissionProbability::theZ

protected

Definition at line 91 of file G4VEmissionProbability.hh.

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

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()

◆ ~G4VEmissionProbability()

Member Function Documentation

◆ ComputeProbability()

◆ EmissionProbability()

◆ GetA()

◆ GetProbability()

◆ GetRecoilExcitation()

◆ GetZ()

◆ Initialise()

◆ IntegrateProbability()

◆ ResetIntegrator()

◆ ResetProbability()

◆ SampleEnergy()

◆ SetDecayKinematics()

◆ SetEvapExcitation()

Member Data Documentation

◆ OPTxs

◆ pEvapMass

◆ pG4pow

◆ pMass

◆ pNuclearLevelData

◆ pProbability

◆ pResMass

◆ pVerbose

◆ resA

◆ resZ

◆ theA

◆ theZ