G4E1Probability Class Reference

#include <G4E1Probability.hh>

Inheritance diagram for G4E1Probability:

Collaboration diagram for G4E1Probability:

Public Member Functions
	G4E1Probability ()
virtual	~G4E1Probability ()
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)

Private Member Functions
	G4E1Probability (const G4E1Probability &right)
const G4E1Probability &	operator= (const G4E1Probability &right)
G4bool	operator== (const G4E1Probability &right) const
G4bool	operator!= (const G4E1Probability &right) const
G4double	EmissionIntegration (const G4Fragment &frag, G4double lowLim, G4double upLim)

Private Attributes
G4Pow *	fG4pow
G4double	theLevelDensityPerNucleon
G4double	normC
G4double	aLevelDensityParam
G4double	sigma0
G4double	Egdp
G4double	GammaR
G4int	Afrag

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

Detailed Description

Definition at line 49 of file G4E1Probability.hh.

Constructor & Destructor Documentation

G4E1Probability() [1/2]

G4E1Probability::G4E1Probability

(

)

Definition at line 50 of file G4E1Probability.cc.

                                :G4VEmissionProbability()
{
  G4double x = CLHEP::pi*CLHEP::hbarc;
  normC = 1.0 / (x*x);
  theLevelDensityPerNucleon = 0.125/MeV;
  aLevelDensityParam = sigma0 = Egdp = GammaR = 0.0;
  Afrag = 0;
  fG4pow = G4Pow::GetInstance(); 
}

Here is the call graph for this function:

~G4E1Probability()

G4E1Probability::~G4E1Probability ( )

virtual

Definition at line 60 of file G4E1Probability.cc.

{}

G4E1Probability() [2/2]

G4E1Probability::G4E1Probability ( const G4E1Probability &  right )

private

Member Function Documentation

EmissionIntegration()

G4double G4E1Probability::EmissionIntegration ( const G4Fragment &  frag, G4double  lowLim, G4double  upLim  )

private

Definition at line 140 of file G4E1Probability.cc.

{
  // Simple integration
  // VI replace by direct integration over 100 point

  static const G4int numIters = 100;
  G4double Step = (upLim-lowLim)/G4double(numIters);

  G4double res = 0.0;
  G4double x = lowLim - 0.5*Step;

  for(G4int i = 0; i < numIters; ++i) {
    x += Step;
    res += EmissionProbDensity(frag, x);
  }

  if(res > 0.0) { res *= Step; }
  else { res = 0.0; }

  return res;

}

Here is the call graph for this function:

Here is the caller graph for this function:

EmissionProbability()

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

virtual

Implements G4VEmissionProbability.

Definition at line 119 of file G4E1Probability.cc.

{
  // From nuclear fragment properties and the excitation energy, calculate
  // the probability for photon evaporation down to last ground level.
  // fragment = nuclear fragment BEFORE de-excitation

  G4double upperLim = gammaE;
  G4double lowerLim = 0.0; 

  //G4cout << "G4E1Probability::EmissionProbability:  Emin= " << lowerLim
  //     << " Emax= " << upperLim << G4endl;
  if( upperLim - lowerLim <=  tolerance) { return 0.0; } 

  // Need to integrate EmissionProbDensity from lowerLim to upperLim 
  // and multiply by factor 3 (?!)

  G4double integ = EmissionIntegration(frag,lowerLim,upperLim);
  return integ;
}

Here is the call graph for this function:

EmissionProbDensity()

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

Definition at line 64 of file G4E1Probability.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 - gammaE
  // (U = nucleus excitation energy, gammaE = total evaporated photon
  // energy). Fragment = nuclear fragment BEFORE de-excitation

  G4double theProb = 0.0;

  if(gammaE > tolerance) {

    G4double Uexcite = frag.GetExcitationEnergy();
    G4double U = std::max(0.0, Uexcite - gammaE);

    // Need a level density parameter.
    // For now, just use the constant approximation (not reliable near magic
    // nuclei) - is equivalent to G4ConstantLevelDensityParameter class

    if(Afrag != frag.GetA_asInt()) {
      Afrag = frag.GetA_asInt();
      aLevelDensityParam = Afrag*theLevelDensityPerNucleon;

      // Define constants for the photoabsorption cross-section (the reverse
      // process of our de-excitation)
      sigma0 = 2.5 * Afrag * millibarn;
      Egdp   = (40.3 / fG4pow->powZ(Afrag,0.2) )*MeV;
      GammaR = 0.30 * Egdp;
    }
    // VI reduce number of calls to exp 
    G4double levelDens = 
      G4Exp(2*(std::sqrt(aLevelDensityParam*U)-
           std::sqrt(aLevelDensityParam*Uexcite)));
 
    //G4cout<<" Uexcite, gammaE = "<<Uexcite<<"  "<<gammaE<<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;

    // VI implementation 18.05.2010
    G4double gammaE2 = gammaE*gammaE;
    G4double gammaR2 = gammaE2*GammaR*GammaR;
    G4double egdp2   = gammaE2 - Egdp*Egdp;
    G4double sigmaAbs = sigma0*gammaR2/(egdp2*egdp2 + gammaR2); 
    theProb = normC * sigmaAbs * gammaE2 * levelDens;

    //G4cout<<" sigmaAbs = "<<sigmaAbs <<" Probability = "<<theProb<<G4endl;
  }
  return theProb;
}

Here is the call graph for this function:

Here is the caller graph for this function:

operator!=()

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

private

operator=()

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

private

operator==()

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

private

Member Data Documentation

Afrag

G4int G4E1Probability::Afrag

private

Definition at line 81 of file G4E1Probability.hh.

aLevelDensityParam

G4double G4E1Probability::aLevelDensityParam

private

Definition at line 77 of file G4E1Probability.hh.

Egdp

G4double G4E1Probability::Egdp

private

Definition at line 79 of file G4E1Probability.hh.

fG4pow

G4Pow* G4E1Probability::fG4pow

private

Definition at line 74 of file G4E1Probability.hh.

GammaR

G4double G4E1Probability::GammaR

private

Definition at line 80 of file G4E1Probability.hh.

normC

G4double G4E1Probability::normC

private

Definition at line 76 of file G4E1Probability.hh.

sigma0

G4double G4E1Probability::sigma0

private

Definition at line 78 of file G4E1Probability.hh.

theLevelDensityPerNucleon

G4double G4E1Probability::theLevelDensityPerNucleon

private

Definition at line 75 of file G4E1Probability.hh.

---

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

• G4E1Probability.hh
• G4E1Probability.cc