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

(

)

Definition at line 41 of file G4E1SingleProbability1.cc.

{}

G4E1SingleProbability1::~G4E1SingleProbability1 ( )

virtual

Definition at line 44 of file G4E1SingleProbability1.cc.

{}

Member Function Documentation

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

virtual

Implements G4VEmissionProbability.

Definition at line 121 of file G4E1SingleProbability1.cc.

References G4Fragment::GetExcitationEnergy().

{

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

  G4double Uexcite = frag.GetExcitationEnergy();
  G4double Uafter = Uexcite - exciteE;

  static const G4double normC = 3.0;

  G4double upperLim = Uexcite;
  G4double lowerLim = Uafter;
  static 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;

}

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

Definition at line 50 of file G4E1SingleProbability1.cc.

References test::a, G4Fragment::GetA_asInt(), G4Fragment::GetExcitationEnergy(), G4Pow::GetInstance(), G4Fragment::GetZ_asInt(), python.hepunit::hbarc, G4ConstantLevelDensityParameter::LevelDensityParameter(), python.hepunit::MeV, python.hepunit::millibarn, python.hepunit::pi, and G4Pow::powZ().

{

  // 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;
 
  static 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