Geant4  10.03.p01
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
G4FissionProbability Class Reference

#include <G4FissionProbability.hh>

Inheritance diagram for G4FissionProbability:
Collaboration diagram for G4FissionProbability:

Public Member Functions

 G4FissionProbability ()
 
 ~G4FissionProbability ()
 
G4double EmissionProbability (const G4Fragment &fragment, G4double MaximalKineticEnergy)
 
void SetEvaporationLevelDensityParameter (G4VLevelDensityParameter *aLevelDensity)
 
void SetFissionLevelDensityParameter (G4VLevelDensityParameter *aLevelDensity)
 
- Public Member Functions inherited from G4VEmissionProbability
 G4VEmissionProbability ()
 
virtual ~G4VEmissionProbability ()
 
void Initialise ()
 
void SetOPTxs (G4int opt)
 
void UseSICB (G4bool use)
 

Additional Inherited Members

- Protected Attributes inherited from G4VEmissionProbability
G4int OPTxs
 
G4bool useSICB
 
G4double LevelDensity
 
G4PowfG4pow
 
G4PairingCorrectionfPairCorr
 

Detailed Description

Definition at line 42 of file G4FissionProbability.hh.

Constructor & Destructor Documentation

G4FissionProbability::G4FissionProbability ( )

Definition at line 43 of file G4FissionProbability.cc.

G4FissionProbability::~G4FissionProbability ( )

Definition at line 50 of file G4FissionProbability.cc.

51 {
52  if (ownEvapLDP) delete theEvapLDP;
53  if (ownFissLDP) delete theFissLDP;
54 }

Member Function Documentation

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

Implements G4VEmissionProbability.

Definition at line 58 of file G4FissionProbability.cc.

61 {
62  if (MaximalKineticEnergy <= 0.0) return 0.0;
63  G4int A = fragment.GetA_asInt();
64  G4int Z = fragment.GetZ_asInt();
65  G4double U = fragment.GetExcitationEnergy();
66 
67  G4double Ucompound = U -
69 
70  G4double Ufission = U -
72 
73  G4double SystemEntropy =
74  2.0*std::sqrt(theEvapLDP->LevelDensityParameter(A,Z,Ucompound)*Ucompound);
75 
76  G4double afission = theFissLDP->LevelDensityParameter(A,Z,Ufission);
77 
78  G4double Cf = 2.0*std::sqrt(afission*MaximalKineticEnergy);
79 
80  // G4double Q1 = 1.0 + (Cf - 1.0)*G4Exp(Cf);
81  // G4double Q2 = 4.0*pi*afission*G4Exp(SystemEntropy);
82 
83  // G4double probability = Q1/Q2;
84 
85  G4double Exp1 = 0.0;
86  if (SystemEntropy <= 160.0) { Exp1 = G4Exp(-SystemEntropy); }
87  // @@@@@@@@@@@@@@@@@ hpw changed max to min - cannot notify vicente now
88  G4double Exp2 = G4Exp( std::min(300.0,Cf-SystemEntropy) );
89 
90  // JMQ 14/02/09 BUG fixed in fission probability (missing parenthesis
91  // at denominator)
92  //AH fix from Vincente: G4double probability =
93  // (Exp1 + (1.0-Cf)*Exp2) / 4.0*pi*afission;
94  // G4double probability = (Exp1 + (Cf-1.0)*Exp2) / 4.0*pi*afission;
95  G4double probability = (Exp1 + (Cf-1.0)*Exp2) / (4.0*pi*afission);
96 
97  return probability;
98 }
int G4int
Definition: G4Types.hh:78
double A(double temperature)
G4int GetA_asInt() const
Definition: G4Fragment.hh:266
G4double GetPairingCorrection(G4int A, G4int Z) const
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
Definition: G4Exp.hh:183
static G4PairingCorrection * GetInstance()
G4double GetFissionPairingCorrection(G4int A, G4int Z) const
virtual G4double LevelDensityParameter(G4int A, G4int Z, G4double U) const =0
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4int GetZ_asInt() const
Definition: G4Fragment.hh:271
static constexpr double pi
Definition: G4SIunits.hh:75
double G4double
Definition: G4Types.hh:76
G4double GetExcitationEnergy() const
Definition: G4Fragment.hh:283

Here is the call graph for this function:

void G4FissionProbability::SetEvaporationLevelDensityParameter ( G4VLevelDensityParameter aLevelDensity)
inline

Definition at line 61 of file G4FissionProbability.hh.

62  {
63  if (ownEvapLDP) delete theEvapLDP;
64  theEvapLDP = aLevelDensity;
65  ownEvapLDP = false;
66  }
void G4FissionProbability::SetFissionLevelDensityParameter ( G4VLevelDensityParameter aLevelDensity)
inline

Definition at line 68 of file G4FissionProbability.hh.

69  {
70  if (ownFissLDP) delete theFissLDP;
71  theFissLDP = aLevelDensity;
72  ownFissLDP = false;
73  }

Here is the caller graph for this function:


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