#include <G4GEMProbabilityVI.hh>

Inheritance diagram for G4GEMProbabilityVI:

Public Member Functions
G4double	ComputeProbability (G4double ekin, G4double CB) override

G4double	ComputeTotalProbability (const G4Fragment &, G4double CB)

virtual G4double	EmissionProbability (const G4Fragment &fragment, G4double anEnergy)

	G4GEMProbabilityVI (G4int anA, G4int aZ, const G4LevelManager *)

G4int	GetA (void) const

G4double	GetProbability () const

G4double	GetRecoilExcitation () const

G4int	GetZ (void) const

void	Initialise ()

void	ResetProbability ()

G4double	SampleEnergy ()

G4Fragment *	SampleEvaporationFragment ()

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

void	SetEvapExcitation (G4double exc)

	~G4GEMProbabilityVI () final

Protected Member Functions
G4double	IntegrateProbability (G4double elow, G4double ehigh, G4double CB)

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

Protected Attributes
G4int	OPTxs

G4double	pEvapMass

G4Pow *	pG4pow

G4double	pMass

G4NuclearLevelData *	pNuclearLevelData

G4double	pProbability

G4double	pResMass

G4int	pVerbose

G4int	resA

G4int	resZ

G4int	theA

G4int	theZ

Private Member Functions
G4double	FindRecoilExcitation (G4double e)

	G4GEMProbabilityVI (const G4GEMProbabilityVI &right)

G4double	I0 (G4double t)

G4double	I1 (G4double t, G4double tx)

G4double	I2 (G4double s0, G4double sx)

G4double	I3 (G4double s0, G4double sx)

G4double	Integrated2DProbability ()

G4bool	operator!= (const G4GEMProbabilityVI &right) const

const G4GEMProbabilityVI &	operator= (const G4GEMProbabilityVI &right)

G4bool	operator== (const G4GEMProbabilityVI &right) const

G4double	ProbabilityDistributionFunction (G4double exc, G4double resExc)

G4Fragment *	Sample2DDistribution ()

Private Attributes
G4double	a0

G4double	a1

G4double	A13

G4double	accuracy

G4double	alphaP

G4double	bCoulomb

G4double	betaP

G4double	coeff

G4double	delta0

G4double	delta1

G4double	eCoulomb

G4double	elimit

G4double	emax

G4double	emin

G4double	fExc

G4double	fExcRes

G4bool	fFD

G4int	fragA

G4int	fragZ

G4double	Gamma

G4bool	isExcited

size_t	length

const G4LevelManager *	lManager

size_t	nbin

G4double	pcoeff

G4double	probmax

G4double	resA13

G4double	U

G4double	Umax

Detailed Description

Definition at line 38 of file G4GEMProbabilityVI.hh.

Constructor & Destructor Documentation

◆ G4GEMProbabilityVI() [1/2]

G4GEMProbabilityVI::G4GEMProbabilityVI	(	G4int	anA,
		G4int	aZ,
		const G4LevelManager *	p
	)

explicit

Definition at line 69 of file G4GEMProbabilityVI.cc.

  : G4VEmissionProbability(aZ, anA), lManager(p)
{
  fragA  = fragZ = 0;
  resA13 = U = delta0 = delta1 = a0 = a1 = probmax = alphaP = betaP = 0.0;
  Umax = bCoulomb = 0.0;
  Gamma = 1.0;
  pcoeff = Gamma*pEvapMass*CLHEP::millibarn
    /((CLHEP::pi*CLHEP::hbarc)*(CLHEP::pi*CLHEP::hbarc)); 
  coeff = CLHEP::fermi*CLHEP::fermi/(CLHEP::pi*CLHEP::hbarc*CLHEP::hbarc);
 
  isExcited = (!lManager || 0.0 == lManager->MaxLevelEnergy()) ? false : true;
  A13 = pG4pow->Z13(theA);
 
  if(0 == aZ) {
    ResetIntegrator(30, 0.25*CLHEP::MeV, 0.02);
  } else {
    ResetIntegrator(30, 0.5*CLHEP::MeV, 0.03);
  }
}

References a0, a1, A13, alphaP, bCoulomb, betaP, coeff, delta0, delta1, CLHEP::fermi, fragA, fragZ, Gamma, CLHEP::hbarc, isExcited, lManager, G4LevelManager::MaxLevelEnergy(), CLHEP::MeV, CLHEP::millibarn, pcoeff, G4VEmissionProbability::pEvapMass, G4VEmissionProbability::pG4pow, CLHEP::pi, probmax, resA13, G4VEmissionProbability::ResetIntegrator(), G4VEmissionProbability::theA, U, Umax, and G4Pow::Z13().

◆ ~G4GEMProbabilityVI()

G4GEMProbabilityVI::~G4GEMProbabilityVI ( )

final

Definition at line 90 of file G4GEMProbabilityVI.cc.

91{}

◆ G4GEMProbabilityVI() [2/2]

G4GEMProbabilityVI::G4GEMProbabilityVI ( const G4GEMProbabilityVI & right )

private

Member Function Documentation

◆ ComputeProbability()

G4double G4GEMProbabilityVI::ComputeProbability	(	G4double	ekin,
		G4double	CB
	)

overridevirtual

Reimplemented from G4VEmissionProbability.

Definition at line 151 of file G4GEMProbabilityVI.cc.

{ 
  // abnormal case - should never happens
  if(pMass < pEvapMass + pResMass) { return 0.0; }
    
  const G4double m02   = pMass*pMass;
  const G4double m12   = pEvapMass*pEvapMass;
  const G4double mres  = std::sqrt(m02 + m12 - 2.*pMass*(pEvapMass + ekin));
 
  G4double excRes = std::max(mres - pResMass, 0.0);
  a1 = pNuclearLevelData->GetLevelDensity(resZ,resA,excRes);
  G4double prob = ProbabilityDistributionFunction(0.0, excRes);
 
  //G4cout<<"### G4GEMProbabilityVI::ComputeProbability: Ekin(MeV)= "<<ekin 
  //<< " excRes(MeV)= " << excRes << " prob= " << prob << << G4endl;
  return prob;
}

References a1, G4NuclearLevelData::GetLevelDensity(), G4INCL::Math::max(), G4VEmissionProbability::pEvapMass, G4VEmissionProbability::pMass, G4VEmissionProbability::pNuclearLevelData, G4VEmissionProbability::pResMass, ProbabilityDistributionFunction(), G4VEmissionProbability::resA, and G4VEmissionProbability::resZ.

◆ ComputeTotalProbability()

G4double G4GEMProbabilityVI::ComputeTotalProbability	(	const G4Fragment &	fragment,
		G4double	CB
	)

Definition at line 93 of file G4GEMProbabilityVI.cc.

{
  fragA = fragment.GetA_asInt();
  fragZ = fragment.GetZ_asInt();
 
  bCoulomb = CB;
  U = fragment.GetExcitationEnergy();
  delta0 = pNuclearLevelData->GetPairingCorrection(fragZ,fragA);
  delta1 = pNuclearLevelData->GetPairingCorrection(resZ,resA);
  Umax = pMass - pEvapMass - pResMass - CB;
  if(0.0 >= Umax) { return 0.0; }
 
  resA13 = pG4pow->Z13(resA);
  a0 = pNuclearLevelData->GetLevelDensity(fragZ,fragA,U);
 
  G4double C = 0.0;
  G4int Z2 = theZ*theZ;
  G4int Z3 = Z2*theZ;
  G4int Z4 = Z2*Z2;
 
  if(resA >= 50) {
    C = -0.10/(G4double)theA;
  } else if(resZ > 20) {
    C = (0.123482-0.00534691*theZ-0.0000610624*Z2+5.93719*1e-7*Z3+
     1.95687*1e-8*Z4)/(G4double)theA;
  }
  if(0 == theZ) {
    alphaP = 0.76+1.93/resA13;
    betaP = (1.66/(resA13*resA13)-0.05)*CLHEP::MeV/alphaP;
  } else {
    alphaP = 1.0 + C;
    betaP = - bCoulomb;
  }
  if(isExcited) {
    pProbability = Integrated2DProbability();
 
  } else {
    const G4double twoMass = pMass + pMass;
    const G4double evapMass2 = pEvapMass*pEvapMass;
    G4double ekinmax = 
     ((pMass-pResMass)*(pMass+pResMass) + evapMass2)/twoMass - pEvapMass;
    G4double ekinmin = 
      std::max((CB*(twoMass - CB) + evapMass2)/twoMass - pEvapMass,0.0);
    if(ekinmax <= ekinmin) { return 0.0; }
    pProbability = IntegrateProbability(ekinmin, ekinmax, CB);
  }
  /*  
  G4cout << "G4GEMProbabilityVI: Z= " << theZ << " A= " << theA 
     << " resZ= " << resZ << " resA= " << resA 
     << " fragZ= " << fragZ << " fragA= " << fragA 
         << " prob= " << pProbability 
     << "\n   U= " << U << " Umax= " << Umax << " d0= " << delta0 
         << " a0= " << a0 << G4endl;
  */
  return pProbability;
}

Referenced by G4GEMChannelVI::GetEmissionProbability().

◆ EmissionProbability()

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

virtualinherited

Reimplemented in G4FissionProbability, and G4GEMProbability.

Definition at line 71 of file G4VEmissionProbability.cc.

{
  return 0.0;
}

Referenced by G4CompetitiveFission::GetEmissionProbability().

◆ FindRecoilExcitation()

G4double G4VEmissionProbability::FindRecoilExcitation ( G4double e )

privateinherited

Definition at line 187 of file G4VEmissionProbability.cc.

{
  fExcRes = 0.0;
  G4double mass = pEvapMass + fExc;
  // abnormal case - should never happens
  if(pMass < mass + pResMass) { return 0.0; }
    
  G4double m02   = pMass*pMass;
  G4double m12   = mass*mass;
  G4double m22   = pResMass*pResMass;
  G4double mres  = std::sqrt(m02 + m12 - 2.*pMass*(mass + e));
 
  fExcRes = mres - pResMass;
  const G4double tolerance = 0.1*CLHEP::keV;
 
  if(pVerbose > 1) {
    G4cout << "### G4VEmissionProbability::FindRecoilExcitation for resZ= " 
           << resZ << " resA= " << resA 
           << " evaporated Z= " << theZ << " A= " << theA
       << " Ekin= " << e << " Eexc= " << fExcRes << G4endl;
  }
 
  // residual nucleus is in the ground state
  if(fExcRes < tolerance) {
    fExcRes = 0.0;
    //G4cout<<"Ground state Ekin= "<< 0.5*(m02 + m12 - m22)/pMass - mass<<G4endl;
    return std::max(0.5*(m02 + m12 - m22)/pMass - mass,0.0);
  }
  // select final state excitation
  auto lManager = pNuclearLevelData->GetLevelManager(resZ, resA);
  if(!lManager) { return e; }
 
  //G4cout<<"ExcMax= "<< lManager->MaxLevelEnergy()<<" CB= "<<eCoulomb<<G4endl;
  // levels are not known
  if(fExcRes > lManager->MaxLevelEnergy() + tolerance) { return e; }
 
  // find level
  auto idx = lManager->NearestLevelIndex(fExcRes);
  //G4cout << "idx= " << idx << " Exc= " << fExcRes 
  //     << " Elevel= " << lManager->LevelEnergy(idx) << G4endl;
  for(; idx > 0; --idx) {
    fExcRes = lManager->LevelEnergy(idx);
    // excited level
    if(pMass > mass + pResMass + fExcRes && lManager->FloatingLevel(idx) == 0) { 
      G4double massR = pResMass + fExcRes;
      G4double mr2 = massR*massR;
      //G4cout << "Result idx= " << idx << " Eexc= " << fExcRes
      //     << " Ekin= " << 0.5*(m02 + m12 - mr2)/pMass - mass << G4endl;
      return std::max(0.5*(m02 + m12 - mr2)/pMass - mass,0.0);
    }
  }
  // ground level
  fExcRes = 0.0;
  //G4cout << "Ground state Ekin= " << 0.5*(m02 + m12 - m22)/pMass - mass << G4endl;
  return std::max(0.5*(m02 + m12 - m22)/pMass - mass,0.0);
}

References G4VEmissionProbability::fExc, G4VEmissionProbability::fExcRes, G4cout, G4endl, G4NuclearLevelData::GetLevelManager(), CLHEP::keV, G4INCL::Math::max(), G4LevelManager::NearestLevelIndex(), G4VEmissionProbability::pEvapMass, G4VEmissionProbability::pMass, G4VEmissionProbability::pNuclearLevelData, G4VEmissionProbability::pResMass, G4VEmissionProbability::pVerbose, G4VEmissionProbability::resA, G4VEmissionProbability::resZ, G4VEmissionProbability::theA, and G4VEmissionProbability::theZ.

Referenced by G4VEmissionProbability::SampleEnergy().

◆ GetA()

G4int G4VEmissionProbability::GetA ( void ) const

inlineinherited

Definition at line 63 of file G4VEmissionProbability.hh.

63{ return theA; }

References G4VEmissionProbability::theA.

◆ GetProbability()

G4double G4VEmissionProbability::GetProbability ( ) const

inlineinherited

Definition at line 75 of file G4VEmissionProbability.hh.

75{ return pProbability; };

References G4VEmissionProbability::pProbability.

◆ GetRecoilExcitation()

G4double G4VEmissionProbability::GetRecoilExcitation ( ) const

inlineinherited

Definition at line 71 of file G4VEmissionProbability.hh.

71{ return fExcRes; };

References G4VEmissionProbability::fExcRes.

◆ GetZ()

G4int G4VEmissionProbability::GetZ ( void ) const

inlineinherited

Definition at line 61 of file G4VEmissionProbability.hh.

61{ return theZ; }

References G4VEmissionProbability::theZ.

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

◆ I0()

G4double G4GEMProbabilityVI::I0 ( G4double t )

private

Definition at line 262 of file G4GEMProbabilityVI.cc.

{
  return G4Exp(t) - 1.0;
}

References G4Exp().

Referenced by ProbabilityDistributionFunction().

◆ I1()

G4double G4GEMProbabilityVI::I1	(	G4double	t,
		G4double	tx
	)

private

Definition at line 267 of file G4GEMProbabilityVI.cc.

{
  return (t - tx + 1.0)*G4Exp(tx) - t - 1.0;
}

References G4Exp().

Referenced by ProbabilityDistributionFunction().

◆ I2()

G4double G4GEMProbabilityVI::I2	(	G4double	s0,
		G4double	sx
	)

private

Definition at line 272 of file G4GEMProbabilityVI.cc.

{
  G4double S = 1.0/std::sqrt(s0);
  G4double Sx = 1.0/std::sqrt(sx);
  
  G4double p1 = S*S*S*( 1.0 + S*S*( 1.5 + 3.75*S*S) );
  G4double p2 = Sx*Sx*Sx*( 1.0 + Sx*Sx*( 1.5 + 3.75*Sx*Sx) )*G4Exp(sx-s0);
  
  return p1-p2;
}

References G4Exp(), S(), and G4InuclParticleNames::s0.

Referenced by ProbabilityDistributionFunction().

◆ I3()

G4double G4GEMProbabilityVI::I3	(	G4double	s0,
		G4double	sx
	)

private

Definition at line 283 of file G4GEMProbabilityVI.cc.

{
  G4double s2 = s0*s0;
  G4double sx2 = sx*sx;
  G4double S = 1.0/std::sqrt(s0);
  G4double S2 = S*S;
  G4double Sx = 1.0/std::sqrt(sx);
  G4double Sx2 = Sx*Sx;
  
  G4double p1 = S *(2.0 + S2 *( 4.0 + S2 *( 13.5 + S2 *( 60.0 + S2 * 325.125 ))));
  G4double p2 = Sx*Sx2 *((s2-sx2) + Sx2 *((1.5*s2+0.5*sx2) 
          + Sx2 *((3.75*s2+0.25*sx2) + Sx2 *((12.875*s2+0.625*sx2) 
          + Sx2 *((59.0625*s2+0.9375*sx2) + Sx2 *(324.8*s2+3.28*sx2))))));
  p2 *= G4Exp(sx-s0);
  return p1-p2; 
}

References G4Exp(), S(), and G4InuclParticleNames::s0.

Referenced by ProbabilityDistributionFunction().

◆ Initialise()

void G4VEmissionProbability::Initialise ( )

inherited

Definition at line 56 of file G4VEmissionProbability.cc.

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

References G4VEmissionProbability::fFD, G4DeexPrecoParameters::GetDeexModelType(), G4DeexPrecoParameters::GetDiscreteExcitationFlag(), G4NuclearLevelData::GetParameters(), G4DeexPrecoParameters::GetVerbose(), G4VEmissionProbability::OPTxs, G4VEmissionProbability::pNuclearLevelData, and G4VEmissionProbability::pVerbose.

Referenced by G4EvaporationChannel::Initialise().

◆ Integrated2DProbability()

G4double G4GEMProbabilityVI::Integrated2DProbability ( )

private

Definition at line 179 of file G4GEMProbabilityVI.cc.

{
  return 0.0;
}

Referenced by ComputeTotalProbability().

◆ IntegrateProbability()

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

protectedinherited

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;
}

References G4VEmissionProbability::accuracy, G4VEmissionProbability::ComputeProbability(), G4VEmissionProbability::eCoulomb, G4VEmissionProbability::elimit, G4VEmissionProbability::emax, G4VEmissionProbability::emin, G4cout, G4endl, G4VEmissionProbability::length, G4INCL::Math::max(), CLHEP::MeV, G4VEmissionProbability::nbin, G4VEmissionProbability::pProbability, G4VEmissionProbability::probmax, and G4VEmissionProbability::pVerbose.

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

◆ operator!=()

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

private

◆ operator=()

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

private

◆ operator==()

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

private

◆ ProbabilityDistributionFunction()

G4double G4GEMProbabilityVI::ProbabilityDistributionFunction	(	G4double	exc,
		G4double	resExc
	)

private

Definition at line 184 of file G4GEMProbabilityVI.cc.

{
  G4double Ux = (2.5 + 150.0/G4double(resA))*CLHEP::MeV;
  G4double Ex = Ux + delta1;
  G4double T  = 1.0/(std::sqrt(a0/Ux) - 1.5/Ux);
  G4double E0 = Ex - T*(G4Log(T) - G4Log(a0)*0.25 
    - 1.25*G4Log(Ux) + 2.0*std::sqrt(a0*Ux));
 
  G4double UxCN = (2.5 + 150.0/(G4double)theA)*CLHEP::MeV;
  G4double ExCN = UxCN + delta0;
  G4double TCN  = 1.0/(std::sqrt(a0/UxCN) - 1.5/UxCN);
 
  G4double mass1 = pEvapMass + exc;
  G4double mass2 = pResMass + resExc;
 
  G4double maxKinEnergy = std::max(0.5*((pMass - mass2)*(pMass + mass2) 
                    + mass1*mass1)/pMass - mass1, 0.0);
 
  G4double Width = 0.0;
  G4double t = maxKinEnergy/T;
  if ( maxKinEnergy < Ex ) {
    Width = (I1(t,t)*T + (betaP+bCoulomb)*I0(t))/G4Exp(E0/T);
 
  } else {
 
    G4double tx = Ex/T;
    G4double s0 = 2.0*std::sqrt(a0*(maxKinEnergy-delta0));
    G4double sx = 2.0*std::sqrt(a0*(Ex-delta0));
 
    // VI: protection against FPE exception
    s0 = std::min(s0, 350.);
 
    G4double expE0T = G4Exp(E0/T);
    G4double exps0  = G4Exp(s0);
    const G4double sqrt2 = std::sqrt(2.0);
 
    Width = I1(t,tx)*T/expE0T + I3(s0,sx)*exps0/(sqrt2*a0);
 
    if (0 == theZ) {
      Width += (betaP+bCoulomb)*(I0(tx)/expE0T + 2.0*sqrt2*I2(s0,sx)*exps0);
    }
  }
  Width *= alphaP*pMass;
 
  //JMQ 190709 fix on Rb and  geometrical cross sections according to 
  //           Furihata's paper (JAERI-Data/Code 2001-105, p6)
  G4double Rb = 0.0;
  if (theA > 4) {
    Rb = 1.12*(resA13 + A13) - 0.86*((resA13 + A13)/(resA13*A13))+2.85;
  } else if (theA > 1) {
    Rb=1.5*(resA13 + A13);
  } else {
    Rb = 1.5*resA13;
  }
 
  G4double ild;
  if (exc < ExCN ) {
    G4double E0CN = ExCN - TCN*(G4Log(TCN) - 0.25*G4Log(a0) 
                - 1.25*G4Log(UxCN) 
                + 2.0*std::sqrt(a0*UxCN));
    ild = G4Exp((exc-E0CN)/TCN)/TCN;
  } else {
    G4double x  = exc - delta0;
    G4double x1 = std::sqrt(a0*x);
    ild = G4Exp(2*x1)/(x*std::sqrt(x1));
  }
 
  Width *= (Rb*Rb/ild); 
  return Width;
}

References a0, A13, alphaP, bCoulomb, betaP, delta0, delta1, G4Exp(), G4Log(), I0(), I1(), I2(), I3(), G4INCL::Math::max(), CLHEP::MeV, G4INCL::Math::min(), G4VEmissionProbability::pEvapMass, G4VEmissionProbability::pMass, G4VEmissionProbability::pResMass, G4VEmissionProbability::resA, resA13, G4InuclParticleNames::s0, G4VEmissionProbability::theA, and G4VEmissionProbability::theZ.

Referenced by ComputeProbability().

◆ ResetIntegrator()

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

protectedinherited

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; }
}

References G4VEmissionProbability::accuracy, G4VEmissionProbability::elimit, eps, and G4VEmissionProbability::length.

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

◆ ResetProbability()

void G4VEmissionProbability::ResetProbability ( )

inlineinherited

Definition at line 77 of file G4VEmissionProbability.hh.

77{ pProbability = 0.0; };

References G4VEmissionProbability::pProbability.

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

◆ Sample2DDistribution()

G4Fragment * G4GEMProbabilityVI::Sample2DDistribution ( )

private

Definition at line 256 of file G4GEMProbabilityVI.cc.

{
  G4Fragment* aFragment = nullptr;
  return aFragment;
}

Referenced by SampleEvaporationFragment().

◆ SampleEnergy()

G4double G4VEmissionProbability::SampleEnergy ( )

inherited

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;
}

References G4VEmissionProbability::ComputeProbability(), G4VEmissionProbability::eCoulomb, G4VEmissionProbability::emax, G4VEmissionProbability::emin, G4VEmissionProbability::fFD, G4VEmissionProbability::FindRecoilExcitation(), CLHEP::HepRandomEngine::flat(), g, G4cout, G4endl, CLHEP::detail::n, G4VEmissionProbability::probmax, G4VEmissionProbability::pVerbose, G4VEmissionProbability::theA, and G4VEmissionProbability::theZ.

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

◆ SampleEvaporationFragment()

G4Fragment * G4GEMProbabilityVI::SampleEvaporationFragment ( )

Definition at line 169 of file G4GEMProbabilityVI.cc.

{
  if(isExcited) { return Sample2DDistribution(); } 
  G4double ekin = SampleEnergy();
  G4LorentzVector lv(std::sqrt(ekin*(ekin + 2.0*pEvapMass))
                     *G4RandomDirection(), ekin + pEvapMass);
  G4Fragment* evFragment = new G4Fragment(theA, theZ, lv);
  return evFragment;
}

References G4RandomDirection(), isExcited, G4VEmissionProbability::pEvapMass, Sample2DDistribution(), G4VEmissionProbability::SampleEnergy(), G4VEmissionProbability::theA, and G4VEmissionProbability::theZ.

Referenced by G4GEMChannelVI::EmittedFragment().

◆ SetDecayKinematics()

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

inlineinherited

Definition at line 131 of file G4VEmissionProbability.hh.

{
  resZ = Z;
  resA = A;
  pMass = fmass;
  pResMass = rmass;
}

References A, G4VEmissionProbability::pMass, G4VEmissionProbability::pResMass, G4VEmissionProbability::resA, G4VEmissionProbability::resZ, and Z.

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

◆ SetEvapExcitation()

void G4VEmissionProbability::SetEvapExcitation ( G4double exc )

inlineinherited

Definition at line 73 of file G4VEmissionProbability.hh.

73{ fExc = exc; };

References G4VEmissionProbability::fExc.

Field Documentation

◆ a0

G4double G4GEMProbabilityVI::a0

private

Definition at line 81 of file G4GEMProbabilityVI.hh.

Referenced by ComputeTotalProbability(), G4GEMProbabilityVI(), and ProbabilityDistributionFunction().

◆ a1

G4double G4GEMProbabilityVI::a1

private

Definition at line 81 of file G4GEMProbabilityVI.hh.

Referenced by ComputeProbability(), and G4GEMProbabilityVI().

◆ A13

G4double G4GEMProbabilityVI::A13

private

Definition at line 83 of file G4GEMProbabilityVI.hh.

Referenced by G4GEMProbabilityVI(), and ProbabilityDistributionFunction().

◆ accuracy

G4double G4VEmissionProbability::accuracy

privateinherited

Definition at line 124 of file G4VEmissionProbability.hh.

Referenced by G4VEmissionProbability::IntegrateProbability(), and G4VEmissionProbability::ResetIntegrator().

◆ alphaP

G4double G4GEMProbabilityVI::alphaP

private

Definition at line 82 of file G4GEMProbabilityVI.hh.

Referenced by ComputeTotalProbability(), G4GEMProbabilityVI(), and ProbabilityDistributionFunction().

◆ bCoulomb

G4double G4GEMProbabilityVI::bCoulomb

private

Definition at line 79 of file G4GEMProbabilityVI.hh.

Referenced by ComputeTotalProbability(), G4GEMProbabilityVI(), and ProbabilityDistributionFunction().

◆ betaP

G4double G4GEMProbabilityVI::betaP

private

Definition at line 82 of file G4GEMProbabilityVI.hh.

Referenced by ComputeTotalProbability(), G4GEMProbabilityVI(), and ProbabilityDistributionFunction().

◆ coeff

G4double G4GEMProbabilityVI::coeff

private

Definition at line 89 of file G4GEMProbabilityVI.hh.

Referenced by G4GEMProbabilityVI().

◆ delta0

G4double G4GEMProbabilityVI::delta0

private

Definition at line 81 of file G4GEMProbabilityVI.hh.

Referenced by ComputeTotalProbability(), G4GEMProbabilityVI(), and ProbabilityDistributionFunction().

◆ delta1

G4double G4GEMProbabilityVI::delta1

private

Definition at line 81 of file G4GEMProbabilityVI.hh.

Referenced by ComputeTotalProbability(), G4GEMProbabilityVI(), and ProbabilityDistributionFunction().

◆ eCoulomb

G4double G4VEmissionProbability::eCoulomb

privateinherited

Definition at line 123 of file G4VEmissionProbability.hh.

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

◆ elimit

G4double G4VEmissionProbability::elimit

privateinherited

Definition at line 122 of file G4VEmissionProbability.hh.

Referenced by G4VEmissionProbability::IntegrateProbability(), and G4VEmissionProbability::ResetIntegrator().

◆ emax

G4double G4VEmissionProbability::emax

privateinherited

Definition at line 121 of file G4VEmissionProbability.hh.

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

◆ emin

G4double G4VEmissionProbability::emin

privateinherited

Definition at line 120 of file G4VEmissionProbability.hh.

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

◆ fExc

G4double G4VEmissionProbability::fExc

privateinherited

Definition at line 117 of file G4VEmissionProbability.hh.

Referenced by G4VEmissionProbability::FindRecoilExcitation(), and G4VEmissionProbability::SetEvapExcitation().

◆ fExcRes

G4double G4VEmissionProbability::fExcRes

privateinherited

Definition at line 118 of file G4VEmissionProbability.hh.

Referenced by G4VEmissionProbability::FindRecoilExcitation(), and G4VEmissionProbability::GetRecoilExcitation().

◆ fFD

G4bool G4VEmissionProbability::fFD

privateinherited

Definition at line 127 of file G4VEmissionProbability.hh.

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

◆ fragA

G4int G4GEMProbabilityVI::fragA

private

Definition at line 76 of file G4GEMProbabilityVI.hh.

Referenced by ComputeTotalProbability(), and G4GEMProbabilityVI().

◆ fragZ

G4int G4GEMProbabilityVI::fragZ

private

Definition at line 77 of file G4GEMProbabilityVI.hh.

Referenced by ComputeTotalProbability(), and G4GEMProbabilityVI().

◆ Gamma

G4double G4GEMProbabilityVI::Gamma

private

Definition at line 88 of file G4GEMProbabilityVI.hh.

Referenced by G4GEMProbabilityVI().

◆ isExcited

G4bool G4GEMProbabilityVI::isExcited

private

Definition at line 94 of file G4GEMProbabilityVI.hh.

Referenced by ComputeTotalProbability(), G4GEMProbabilityVI(), and SampleEvaporationFragment().

◆ length

size_t G4VEmissionProbability::length

privateinherited

Definition at line 114 of file G4VEmissionProbability.hh.

Referenced by G4VEmissionProbability::G4VEmissionProbability(), G4VEmissionProbability::IntegrateProbability(), and G4VEmissionProbability::ResetIntegrator().

◆ lManager

const G4LevelManager* G4GEMProbabilityVI::lManager

private

Definition at line 74 of file G4GEMProbabilityVI.hh.

Referenced by G4GEMProbabilityVI().

◆ nbin

size_t G4VEmissionProbability::nbin

privateinherited

Definition at line 115 of file G4VEmissionProbability.hh.

Referenced by G4VEmissionProbability::G4VEmissionProbability(), and G4VEmissionProbability::IntegrateProbability().

◆ OPTxs

G4int G4VEmissionProbability::OPTxs

protectedinherited

Definition at line 89 of file G4VEmissionProbability.hh.

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

◆ pcoeff

G4double G4GEMProbabilityVI::pcoeff

private

Definition at line 90 of file G4GEMProbabilityVI.hh.

Referenced by G4GEMProbabilityVI().

◆ pEvapMass

G4double G4VEmissionProbability::pEvapMass

protectedinherited

Definition at line 97 of file G4VEmissionProbability.hh.

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

◆ pG4pow

G4Pow* G4VEmissionProbability::pG4pow

protectedinherited

Definition at line 102 of file G4VEmissionProbability.hh.

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

◆ pMass

G4double G4VEmissionProbability::pMass

protectedinherited

Definition at line 96 of file G4VEmissionProbability.hh.

Referenced by ComputeProbability(), G4EvaporationProbability::ComputeProbability(), ComputeTotalProbability(), G4VEmissionProbability::FindRecoilExcitation(), ProbabilityDistributionFunction(), and G4VEmissionProbability::SetDecayKinematics().

◆ pNuclearLevelData

G4NuclearLevelData* G4VEmissionProbability::pNuclearLevelData

protectedinherited

Definition at line 101 of file G4VEmissionProbability.hh.

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

◆ pProbability

G4double G4VEmissionProbability::pProbability

protectedinherited

Definition at line 99 of file G4VEmissionProbability.hh.

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

◆ pResMass

G4double G4VEmissionProbability::pResMass

protectedinherited

Definition at line 98 of file G4VEmissionProbability.hh.

Referenced by ComputeProbability(), G4EvaporationProbability::ComputeProbability(), ComputeTotalProbability(), G4VEmissionProbability::FindRecoilExcitation(), ProbabilityDistributionFunction(), and G4VEmissionProbability::SetDecayKinematics().

◆ probmax

G4double G4GEMProbabilityVI::probmax

private

Definition at line 92 of file G4GEMProbabilityVI.hh.

Referenced by G4GEMProbabilityVI().

◆ pVerbose

G4int G4VEmissionProbability::pVerbose

protectedinherited

Definition at line 90 of file G4VEmissionProbability.hh.

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

◆ resA

G4int G4VEmissionProbability::resA

protectedinherited

Definition at line 94 of file G4VEmissionProbability.hh.

Referenced by ComputeProbability(), G4EvaporationProbability::ComputeProbability(), ComputeTotalProbability(), G4EvaporationProbability::CrossSection(), G4VEmissionProbability::FindRecoilExcitation(), ProbabilityDistributionFunction(), G4VEmissionProbability::SetDecayKinematics(), and G4EvaporationProbability::TotalProbability().

◆ resA13

G4double G4GEMProbabilityVI::resA13

private

Definition at line 80 of file G4GEMProbabilityVI.hh.

Referenced by ComputeTotalProbability(), G4GEMProbabilityVI(), and ProbabilityDistributionFunction().

◆ resZ

G4int G4VEmissionProbability::resZ

protectedinherited

Definition at line 93 of file G4VEmissionProbability.hh.

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

◆ theA

G4int G4VEmissionProbability::theA

protectedinherited

Definition at line 92 of file G4VEmissionProbability.hh.

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

◆ theZ

G4int G4VEmissionProbability::theZ

protectedinherited

Definition at line 91 of file G4VEmissionProbability.hh.

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

◆ U

G4double G4GEMProbabilityVI::U

private

Definition at line 81 of file G4GEMProbabilityVI.hh.

Referenced by ComputeTotalProbability(), and G4GEMProbabilityVI().

◆ Umax

G4double G4GEMProbabilityVI::Umax

private

Definition at line 83 of file G4GEMProbabilityVI.hh.

Referenced by ComputeTotalProbability(), and G4GEMProbabilityVI().

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

source/processes/hadronic/models/de_excitation/gem_evaporation/include/G4GEMProbabilityVI.hh
source/processes/hadronic/models/de_excitation/gem_evaporation/src/G4GEMProbabilityVI.cc

Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4GEMProbabilityVI() [1/2]

◆ ~G4GEMProbabilityVI()

◆ G4GEMProbabilityVI() [2/2]

Member Function Documentation

◆ ComputeProbability()

◆ ComputeTotalProbability()

◆ EmissionProbability()

◆ FindRecoilExcitation()

◆ GetA()

◆ GetProbability()

◆ GetRecoilExcitation()

◆ GetZ()

◆ I0()

◆ I1()

◆ I2()

◆ I3()

◆ Initialise()

◆ Integrated2DProbability()

◆ IntegrateProbability()

◆ operator!=()

◆ operator=()

◆ operator==()

◆ ProbabilityDistributionFunction()

◆ ResetIntegrator()

◆ ResetProbability()

◆ Sample2DDistribution()

◆ SampleEnergy()

◆ SampleEvaporationFragment()

◆ SetDecayKinematics()

◆ SetEvapExcitation()

Field Documentation

◆ a0

◆ a1

◆ A13

◆ accuracy

◆ alphaP

◆ bCoulomb

◆ betaP

◆ coeff

◆ delta0

◆ delta1

◆ eCoulomb

◆ elimit

◆ emax

◆ emin

◆ fExc

◆ fExcRes

◆ fFD

◆ fragA

◆ fragZ

◆ Gamma

◆ isExcited

◆ length

◆ lManager

◆ nbin

◆ OPTxs

◆ pcoeff

◆ pEvapMass

◆ pG4pow

◆ pMass

◆ pNuclearLevelData

◆ pProbability

◆ pResMass

◆ probmax

◆ pVerbose

◆ resA

◆ resA13

◆ resZ

◆ theA

◆ theZ

◆ U

◆ Umax