G4FissionProbability Class Reference

#include <G4FissionProbability.hh>

Inheritance diagram for G4FissionProbability:

Public Member Functions
	G4FissionProbability ()
	~G4FissionProbability ()
G4double	EmissionProbability (const G4Fragment &fragment, G4double MaximalKineticEnergy)
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 43 of file G4FissionProbability.hh.

Constructor & Destructor Documentation

G4FissionProbability()

G4FissionProbability::G4FissionProbability

(

)

Definition at line 40 of file G4FissionProbability.cc.

{}

~G4FissionProbability()

G4FissionProbability::~G4FissionProbability

(

)

Definition at line 43 of file G4FissionProbability.cc.

{}

Member Function Documentation

EmissionProbability()

G4double G4FissionProbability::EmissionProbability ( const G4Fragment &  fragment, G4double  MaximalKineticEnergy  )

virtual

Implements G4VEmissionProbability.

Definition at line 48 of file G4FissionProbability.cc.

{
  if (MaximalKineticEnergy <= 0.0) return 0.0;
  G4int A = fragment.GetA_asInt();
  G4int Z = fragment.GetZ_asInt();
  G4double U = fragment.GetExcitationEnergy();
  
  G4double Ucompound = U - 
    G4PairingCorrection::GetInstance()->GetPairingCorrection(A,Z);
 
  G4double Ufission = U - 
    G4PairingCorrection::GetInstance()->GetFissionPairingCorrection(A,Z);
  
  G4double SystemEntropy = 
    2.0*std::sqrt(theEvapLDP.LevelDensityParameter(A,Z,Ucompound)*Ucompound);
    
  G4double afission = theFissLDP.LevelDensityParameter(A,Z,Ufission);
 
  G4double Cf = 2.0*std::sqrt(afission*MaximalKineticEnergy);
 
  //    G4double Q1 = 1.0 + (Cf - 1.0)*std::exp(Cf);
  //    G4double Q2 = 4.0*pi*afission*std::exp(SystemEntropy);
  
  //    G4double probability = Q1/Q2;
   
  G4double Exp1 = 0.0;
  if (SystemEntropy <= 160.0) { Exp1 = std::exp(-SystemEntropy); }
  // @@@@@@@@@@@@@@@@@ hpw changed max to min - cannot notify vicente now
  G4double Exp2 = std::exp( std::min(700.0,Cf-SystemEntropy) ); 
 
  // JMQ 14/02/09 BUG fixed in fission probability (missing parenthesis 
  // at denominator)
  //AH fix from Vincente: G4double probability = 
  //        (Exp1 + (1.0-Cf)*Exp2) / 4.0*pi*afission;
  //    G4double probability = (Exp1 + (Cf-1.0)*Exp2) / 4.0*pi*afission;
  G4double probability = (Exp1 + (Cf-1.0)*Exp2) / (4.0*pi*afission);
 
  return probability;
}

---

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

• G4FissionProbability.hh
• G4FissionProbability.cc