G4FissionParameters Class Reference

#include <G4FissionParameters.hh>

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

Detailed Description

Definition at line 38 of file G4FissionParameters.hh.

Constructor & Destructor Documentation

G4FissionParameters::G4FissionParameters

(

)

Definition at line 39 of file G4FissionParameters.cc.

  : A1(134),A2(141),A3((A1 + A2)*0.5),As(0.0),Sigma1(0.0),Sigma2(0.0),
    SigmaS(0.),w(0.0)
{}

G4FissionParameters::~G4FissionParameters

(

)

Definition at line 44 of file G4FissionParameters.cc.

{}

Member Function Documentation

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

Definition at line 47 of file G4FissionParameters.cc.

{
  // to avoid usage of units
  G4double U = ExEnergy/CLHEP::MeV; 

  As = A*0.5;
    
  if (A <= 235) { Sigma2 = 5.6; }  
  else          { Sigma2 = 5.6 + 0.096*(A-235); }

  Sigma1 = 0.5*Sigma2;     
    
  //JMQ 310509 
  //    if (SigmaS > 20.0) SigmaS = 20.0;
  //   SigmaS*=1.3;
  //JMQ 301009: retuning (after CEM transition prob.have been chosen as default)
  SigmaS = 0.8*G4Exp(0.00553*U + 2.1386); 
    
  G4double x1 = (A1-As)/Sigma1;
  G4double x2 = (A2-As)/Sigma2;
  G4double FasymAsym = 2*G4Exp(-0.5*x2*x2) + G4Exp(-0.5*x1*x1);
    
  G4double x3 = (As-A3)/SigmaS;
  G4double FsymA1A2 = G4Exp(-0.5*x3*x3);
    
  G4double wa = 0.0;
  w = 0.0;
  if (Z >= 90) {  
    if (U <= 16.25) { wa = G4Exp(0.5385*U-9.9564); }  
    else            { wa = G4Exp(0.09197*U-2.7003); }            
  } else if (Z == 89) {  
    wa = G4Exp(0.09197*U-1.0808);
  } else if (Z >= 82) {  
    G4double X = std::max(0.0, FissionBarrier/CLHEP::MeV - 7.5);
    wa = G4Exp(0.09197*(U-X) - 1.0808);
  } else {               // Z < 82
    w = 1001.0;
  }
    
  if (w == 0.0) {
    G4double w1 = std::max(1.03*wa - FasymAsym, 0.0001);
    G4double w2 = std::max(1.0 - FsymA1A2*wa,   0.0001);
    
    w = w1/w2;
      
    if (82 <= Z && Z < 89 && A < 227) { w *= G4Exp(0.3*(227-A)); }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4int G4FissionParameters::GetA1 ( void  ) const

inline

Definition at line 51 of file G4FissionParameters.hh.

{ return A1; }

G4int G4FissionParameters::GetA2 ( void  ) const

inline

Definition at line 52 of file G4FissionParameters.hh.

{ return A2; }

G4double G4FissionParameters::GetAs ( void  ) const

inline

Definition at line 54 of file G4FissionParameters.hh.

{ return As; }

G4double G4FissionParameters::GetSigma1 ( void  ) const

inline

Definition at line 55 of file G4FissionParameters.hh.

{ return Sigma1; }

G4double G4FissionParameters::GetSigma2 ( void  ) const

inline

Definition at line 56 of file G4FissionParameters.hh.

{ return Sigma2; }

G4double G4FissionParameters::GetSigmaS ( void  ) const

inline

Definition at line 57 of file G4FissionParameters.hh.

{ return SigmaS; }

G4double G4FissionParameters::GetW ( void  ) const

inline

Definition at line 58 of file G4FissionParameters.hh.

{ return w; }

---

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

• source/geant4.10.03.p03/source/processes/hadronic/models/de_excitation/fission/include/G4FissionParameters.hh
• source/geant4.10.03.p03/source/processes/hadronic/models/de_excitation/fission/src/G4FissionParameters.cc