G4GEMProbabilityVI Class Reference

#include <G4GEMProbabilityVI.hh>

Inheritance diagram for G4GEMProbabilityVI:

Public Member Functions
	G4GEMProbabilityVI (G4int anA, G4int aZ, const G4LevelManager *)
	~G4GEMProbabilityVI () override=default
G4double	TotalProbability (const G4Fragment &, const G4double tmin, const G4double tmax, const G4double CB, const G4double exEnergy, const G4double exEvap)
G4double	ComputeProbability (G4double ekin, G4double CB) override
G4double	SampleEnergy (const G4double tmin, const G4double tmax, const G4double CB, const G4double exEnergy, const G4double exEvap)
	G4GEMProbabilityVI (const G4GEMProbabilityVI &right)=delete
const G4GEMProbabilityVI &	operator= (const G4GEMProbabilityVI &right)=delete
G4bool	operator== (const G4GEMProbabilityVI &right) const =delete
G4bool	operator!= (const G4GEMProbabilityVI &right) const =delete
Public Member Functions inherited from G4VEmissionProbability
	G4VEmissionProbability (G4int Z, G4int A)
virtual	~G4VEmissionProbability ()=default
virtual 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

Additional Inherited Members
Protected Member Functions inherited from G4VEmissionProbability
void	ResetIntegrator (size_t nbin, G4double de, G4double eps)
G4double	IntegrateProbability (G4double elow, G4double ehigh, G4double CB)
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
G4double	pWidth = 0.0

Detailed Description

Definition at line 36 of file G4GEMProbabilityVI.hh.

Constructor & Destructor Documentation

G4GEMProbabilityVI() [1/2]

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

explicit

Definition at line 42 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);
  }
}

Referenced by G4GEMProbabilityVI(), operator!=(), operator=(), and operator==().

~G4GEMProbabilityVI()

G4GEMProbabilityVI::~G4GEMProbabilityVI ( )

overridedefault

G4GEMProbabilityVI() [2/2]

G4GEMProbabilityVI::G4GEMProbabilityVI ( const G4GEMProbabilityVI & right )

delete

Member Function Documentation

ComputeProbability()

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

overridevirtual

Reimplemented from G4VEmissionProbability.

Definition at line 101 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 = 0.5; //CrossSection(0.0, excRes);
 
  //G4cout<<"### G4GEMProbabilityVI::ComputeProbability: Ekin(MeV)= "<<ekin 
  //<< " excRes(MeV)= " << excRes << " prob= " << prob << << G4endl;
  return prob;
}

operator!=()

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

delete

operator=()

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

delete

operator==()

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

delete

SampleEnergy()

G4double G4GEMProbabilityVI::SampleEnergy	(	const G4double	tmin,
		const G4double	tmax,
		const G4double	CB,
		const G4double	exEnergy,
		const G4double	exEvap )

Definition at line 119 of file G4GEMProbabilityVI.cc.

{
  G4double ekin = tmax - tmin - CB -exEnergy + exEvap;
  return ekin;
}

TotalProbability()

G4double G4GEMProbabilityVI::TotalProbability	(	const G4Fragment &	fragment,
		const G4double	tmin,
		const G4double	tmax,
		const G4double	CB,
		const G4double	exEnergy,
		const G4double	exEvap )

Definition at line 63 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);
  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);
  pProbability += tmax - tmin + exEnergy -exEvap;
  /*  
  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;
}

---

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

• G4GEMProbabilityVI.hh
• G4GEMProbabilityVI.cc