Geant4  10.02.p03
G4FissionParameters Class Reference

#include <G4FissionParameters.hh>

Collaboration diagram for G4FissionParameters:

Public Member Functions

 G4FissionParameters ()
 
 ~G4FissionParameters ()
 
void DefineParameters (G4int A, G4int Z, G4double ExEnergy, G4double FissionBarrier)
 
G4int GetA1 (void) const
 
G4int GetA2 (void) const
 
G4double GetAs (void) const
 
G4double GetSigma1 (void) const
 
G4double GetSigma2 (void) const
 
G4double GetSigmaS (void) const
 
G4double GetW (void) const
 

Private Member Functions

 G4FissionParameters (const G4FissionParameters &right)
 
const G4FissionParametersoperator= (const G4FissionParameters &right)
 
G4bool operator== (const G4FissionParameters &right) const
 
G4bool operator!= (const G4FissionParameters &right) const
 

Private Attributes

G4int A1
 
G4int A2
 
G4double A3
 
G4double As
 
G4double Sigma1
 
G4double Sigma2
 
G4double SigmaS
 
G4double w
 

Detailed Description

Definition at line 38 of file G4FissionParameters.hh.

Constructor & Destructor Documentation

◆ G4FissionParameters() [1/2]

G4FissionParameters::G4FissionParameters ( )

Definition at line 39 of file G4FissionParameters.cc.

Here is the caller graph for this function:

◆ ~G4FissionParameters()

G4FissionParameters::~G4FissionParameters ( )

Definition at line 44 of file G4FissionParameters.cc.

45 {}

◆ G4FissionParameters() [2/2]

G4FissionParameters::G4FissionParameters ( const G4FissionParameters right)
private

Member Function Documentation

◆ DefineParameters()

void G4FissionParameters::DefineParameters ( G4int  A,
G4int  Z,
G4double  ExEnergy,
G4double  FissionBarrier 
)

Definition at line 47 of file G4FissionParameters.cc.

49 {
50  // to avoid usage of units
51  G4double U = ExEnergy/CLHEP::MeV;
52 
53  As = A*0.5;
54 
55  if (A <= 235) { Sigma2 = 5.6; }
56  else { Sigma2 = 5.6 + 0.096*(A-235); }
57 
58  Sigma1 = 0.5*Sigma2;
59 
60  //JMQ 310509
61  // if (SigmaS > 20.0) SigmaS = 20.0;
62  // SigmaS*=1.3;
63  //JMQ 301009: retuning (after CEM transition prob.have been chosen as default)
64  SigmaS = 0.8*G4Exp(0.00553*U + 2.1386);
65 
66  G4double x1 = (A1-As)/Sigma1;
67  G4double x2 = (A2-As)/Sigma2;
68  G4double FasymAsym = 2*G4Exp(-0.5*x2*x2) + G4Exp(-0.5*x1*x1);
69 
70  G4double x3 = (As-A3)/SigmaS;
71  G4double FsymA1A2 = G4Exp(-0.5*x3*x3);
72 
73  G4double wa = 0.0;
74  w = 0.0;
75  if (Z >= 90) {
76  if (U <= 16.25) { wa = G4Exp(0.5385*U-9.9564); }
77  else { wa = G4Exp(0.09197*U-2.7003); }
78  } else if (Z == 89) {
79  wa = G4Exp(0.09197*U-1.0808);
80  } else if (Z >= 82) {
81  G4double X = std::max(0.0, FissionBarrier/CLHEP::MeV - 7.5);
82  wa = G4Exp(0.09197*(U-X) - 1.0808);
83  } else { // Z < 82
84  w = 1001.0;
85  }
86 
87  if (w == 0.0) {
88  G4double w1 = std::max(1.03*wa - FasymAsym, 0.0001);
89  G4double w2 = std::max(1.0 - FsymA1A2*wa, 0.0001);
90 
91  w = w1/w2;
92 
93  if (82 <= Z && Z < 89 && A < 227) { w *= G4Exp(0.3*(227-A)); }
94  }
95 }
Double_t x2[nxs]
Float_t X
double A(double temperature)
Float_t Z
Double_t x1[nxs]
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
Definition: G4Exp.hh:183
double G4double
Definition: G4Types.hh:76
static const double MeV
Here is the call graph for this function:
Here is the caller graph for this function:

◆ GetA1()

G4int G4FissionParameters::GetA1 ( void  ) const
inline

Definition at line 51 of file G4FissionParameters.hh.

51 { return A1; }
Here is the caller graph for this function:

◆ GetA2()

G4int G4FissionParameters::GetA2 ( void  ) const
inline

Definition at line 52 of file G4FissionParameters.hh.

52 { return A2; }
Here is the caller graph for this function:

◆ GetAs()

G4double G4FissionParameters::GetAs ( void  ) const
inline

Definition at line 54 of file G4FissionParameters.hh.

54 { return As; }
Here is the caller graph for this function:

◆ GetSigma1()

G4double G4FissionParameters::GetSigma1 ( void  ) const
inline

Definition at line 55 of file G4FissionParameters.hh.

55 { return Sigma1; }
Here is the caller graph for this function:

◆ GetSigma2()

G4double G4FissionParameters::GetSigma2 ( void  ) const
inline

Definition at line 56 of file G4FissionParameters.hh.

56 { return Sigma2; }
Here is the caller graph for this function:

◆ GetSigmaS()

G4double G4FissionParameters::GetSigmaS ( void  ) const
inline

Definition at line 57 of file G4FissionParameters.hh.

57 { return SigmaS; }
Here is the caller graph for this function:

◆ GetW()

G4double G4FissionParameters::GetW ( void  ) const
inline

Definition at line 58 of file G4FissionParameters.hh.

58 { return w; }
Here is the call graph for this function:
Here is the caller graph for this function:

◆ operator!=()

G4bool G4FissionParameters::operator!= ( const G4FissionParameters right) const
private
Here is the caller graph for this function:

◆ operator=()

const G4FissionParameters& G4FissionParameters::operator= ( const G4FissionParameters right)
private
Here is the caller graph for this function:

◆ operator==()

G4bool G4FissionParameters::operator== ( const G4FissionParameters right) const
private
Here is the caller graph for this function:

Member Data Documentation

◆ A1

G4int G4FissionParameters::A1
private

Definition at line 69 of file G4FissionParameters.hh.

◆ A2

G4int G4FissionParameters::A2
private

Definition at line 70 of file G4FissionParameters.hh.

◆ A3

G4double G4FissionParameters::A3
private

Definition at line 71 of file G4FissionParameters.hh.

◆ As

G4double G4FissionParameters::As
private

Definition at line 74 of file G4FissionParameters.hh.

◆ Sigma1

G4double G4FissionParameters::Sigma1
private

Definition at line 78 of file G4FissionParameters.hh.

◆ Sigma2

G4double G4FissionParameters::Sigma2
private

Definition at line 79 of file G4FissionParameters.hh.

◆ SigmaS

G4double G4FissionParameters::SigmaS
private

Definition at line 82 of file G4FissionParameters.hh.

◆ w

G4double G4FissionParameters::w
private

Definition at line 86 of file G4FissionParameters.hh.


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