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

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

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

91{}

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:

Public Member Functions

Additional Inherited Members