G4E1SingleProbability1 Class Reference

#include <G4E1SingleProbability1.hh>

Inheritance diagram for G4E1SingleProbability1:

Public Member Functions
	G4E1SingleProbability1 ()
virtual	~G4E1SingleProbability1 ()
G4double	EmissionProbability (const G4Fragment &frag, G4double excite)
G4double	EmissionProbDensity (const G4Fragment &frag, G4double ePhoton)
Public Member Functions inherited from G4VEmissionProbability
	G4VEmissionProbability ()
virtual	~G4VEmissionProbability ()
virtual G4double	EmissionProbability (const G4Fragment &fragment, const G4double anEnergy)=0
void	SetOPTxs (G4int opt)
void	UseSICB (G4bool use)

Additional Inherited Members
Protected Attributes inherited from G4VEmissionProbability
G4int	OPTxs
G4bool	useSICB
G4Pow *	fG4pow
G4PairingCorrection *	fPairCorr
G4EvaporationLevelDensityParameter *	theEvapLDPptr

Detailed Description

Definition at line 38 of file G4E1SingleProbability1.hh.

Constructor & Destructor Documentation

G4E1SingleProbability1()

G4E1SingleProbability1::G4E1SingleProbability1

(

)

Definition at line 41 of file G4E1SingleProbability1.cc.

{}

~G4E1SingleProbability1()

G4E1SingleProbability1::~G4E1SingleProbability1 ( )

virtual

Definition at line 44 of file G4E1SingleProbability1.cc.

{}

Member Function Documentation

EmissionProbability()

G4double G4E1SingleProbability1::EmissionProbability ( const G4Fragment &  frag, G4double  excite  )

virtual

Implements G4VEmissionProbability.

Definition at line 121 of file G4E1SingleProbability1.cc.

{
 
  // From nuclear fragment properties and the excitation energy, calculate
  // the probability for photon evaporation down to the level
  // Uexcite-exciteE.
  // fragment = nuclear fragment BEFORE de-excitation
 
  G4double theProb = 0.0;
 
  G4double ScaleFactor = 1.0;     // playing with scale factors
 
  const G4double Uexcite = frag.GetExcitationEnergy();
  G4double Uafter = Uexcite - exciteE;
 
  G4double normC = 3.0;
 
  const G4double upperLim = Uexcite;
  const G4double lowerLim = Uafter;
  const G4int numIters = 25;
 
  // Need to integrate EmissionProbDensity from lowerLim to upperLim 
  // and multiply by normC
 
  G4double integ = normC *
           EmissionIntegration(frag,exciteE,lowerLim,upperLim,numIters);
 
  if(integ > 0.0) theProb = integ;
 
  return theProb * ScaleFactor;
 
}

EmissionProbDensity()

G4double G4E1SingleProbability1::EmissionProbDensity	(	const G4Fragment &	frag,
		G4double	ePhoton
	)

Definition at line 50 of file G4E1SingleProbability1.cc.

{
 
  // Calculate the probability density here
 
  // From nuclear fragment properties and the excitation energy, calculate
  // the probability density for photon evaporation from U to U - exciteE
  // (U = nucleus excitation energy, exciteE = total evaporated photon
  // energy).
  // fragment = nuclear fragment BEFORE de-excitation
 
  G4double theProb = 0.0;
 
  G4int Afrag = frag.GetA_asInt();
  G4int Zfrag = frag.GetZ_asInt();
  G4double Uexcite = frag.GetExcitationEnergy();
 
  if( (Uexcite-exciteE) < 0.0 || exciteE < 0 || Uexcite <= 0) return theProb;
 
  // Need a level density parameter.
  // For now, just use the constant approximation (not reliable near magic
  // nuclei).
 
  G4ConstantLevelDensityParameter a;
  G4double aLevelDensityParam = a.LevelDensityParameter(Afrag,Zfrag,Uexcite);
 
  G4double levelDensBef = std::exp(2.0*std::sqrt(aLevelDensityParam*Uexcite));
  G4double levelDensAft = std::exp(2.0*std::sqrt(aLevelDensityParam*(Uexcite-exciteE)));
 
  // Now form the probability density
 
  // Define constants for the photoabsorption cross-section (the reverse
  // process of our de-excitation)
 
  G4double sigma0 = 2.5 * Afrag * millibarn;  // millibarns
 
  G4double Egdp = (40.3 / G4Pow::GetInstance()->powZ(Afrag,0.2) )*MeV;
  G4double GammaR = 0.30 * Egdp;
 
  const G4double normC = 1.0 / ((pi * hbarc)*(pi * hbarc));
 
  // CD
  //cout<<"  PROB TESTS "<<G4endl;
  //cout<<" hbarc = "<<hbarc<<G4endl;
  //cout<<" pi = "<<pi<<G4endl;
  //cout<<" Uexcite, exciteE = "<<Uexcite<<"  "<<exciteE<<G4endl;
  //cout<<" Uexcite, exciteE = "<<Uexcite*MeV<<"  "<<exciteE*MeV<<G4endl;
  //cout<<" lev density param = "<<aLevelDensityParam<<G4endl;
  //cout<<" level densities = "<<levelDensBef<<"  "<<levelDensAft<<G4endl;
  //cout<<" sigma0 = "<<sigma0<<G4endl;
  //cout<<" Egdp, GammaR = "<<Egdp<<"  "<<GammaR<<G4endl;
  //cout<<" normC = "<<normC<<G4endl;
 
  G4double numerator = sigma0 * exciteE*exciteE * GammaR*GammaR;
  G4double denominator = (exciteE*exciteE - Egdp*Egdp)*
           (exciteE*exciteE - Egdp*Egdp) + GammaR*GammaR*exciteE*exciteE;
 
  G4double sigmaAbs = numerator/denominator; 
 
  theProb = normC * sigmaAbs * exciteE*exciteE *
            levelDensAft/levelDensBef;
 
  // CD
  //cout<<" sigmaAbs = "<<sigmaAbs<<G4endl;
  //cout<<" Probability = "<<theProb<<G4endl;
 
  return theProb;
 
}

---

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

• G4E1SingleProbability1.hh
• G4E1SingleProbability1.cc