G4GEMProbabilityVI Class Reference

#include <G4GEMProbabilityVI.hh>

Inheritance diagram for G4GEMProbabilityVI:

Public Member Functions
	G4GEMProbabilityVI (G4int anA, G4int aZ, const G4LevelManager *)
	~G4GEMProbabilityVI () final
G4double	ComputeTotalProbability (const G4Fragment &, G4double CB)
G4double	ComputeProbability (G4double ekin, G4double CB) override
G4Fragment *	SampleEvaporationFragment ()
Public Member Functions inherited from G4VEmissionProbability
	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

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

Detailed Description

Definition at line 38 of file G4GEMProbabilityVI.hh.

Constructor & Destructor Documentation

G4GEMProbabilityVI()

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

~G4GEMProbabilityVI()

G4GEMProbabilityVI::~G4GEMProbabilityVI ( )

final

Definition at line 90 of file G4GEMProbabilityVI.cc.

{}

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

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

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

Referenced by G4GEMChannelVI::EmittedFragment().

---

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

• G4GEMProbabilityVI.hh
• G4GEMProbabilityVI.cc