G4FissionParameters Class Reference

#include <G4FissionParameters.hh>

Public Member Functions
	G4FissionParameters (G4int A, G4int Z, G4double ExEnergy, G4double FissionBarrier)
	~G4FissionParameters ()
G4double	GetA1 (void) const
G4double	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 41 of file G4FissionParameters.hh.

Constructor & Destructor Documentation

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

Definition at line 43 of file G4FissionParameters.cc.

References G4FissionParameters(), G4INCL::Math::max(), and python.hepunit::MeV.

Referenced by G4FissionParameters().

{
  G4double U = ExEnergy; 
  
  As = A/2.0;
    
  if (A <= 235) Sigma2 = 5.6;  // MeV
  else Sigma2 = 5.6 + 0.096*(A-235); // MeV
    
  Sigma1 = 0.5*Sigma2; // MeV
    
  SigmaS = std::exp(0.00553*U/MeV + 2.1386); // MeV
    
  //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;
  //
    
  G4double FasymAsym = 2.0*std::exp(-((A2-As)*(A2-As))/(2.0*Sigma2*Sigma2)) + 
    std::exp(-((A1-As)*(A1-As))/(2.0*Sigma1*Sigma1));
    
  G4double FsymA1A2 = std::exp(-((As-(A1+A2)/2.0)*(As-(A1+A2)/2.0))
                   /(2.0*SigmaS*SigmaS));
    
    
  G4double wa = 0.0;
  w = 0.0;
  if (Z >= 90) {         // Z >= 90
    if (U <= 16.25) wa = std::exp(0.5385*U/MeV-9.9564);  // U <= 16.25 MeV
    else wa = std::exp(0.09197*U/MeV-2.7003);            // U  > 16.25 MeV
  } else if (Z == 89) {  // Z == 89
    wa = std::exp(0.09197*U-1.0808);
  } else if (Z >= 82) {  //  82 <= Z <= 88
    G4double X = FissionBarrier - 7.5*MeV;
    if (X < 0.0) X = 0.0;
    wa = std::exp(0.09197*(U-X)/MeV-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 *= std::exp(0.3*(227-A));
  }
    
}

G4FissionParameters::~G4FissionParameters

(

)

Definition at line 97 of file G4FissionParameters.cc.

{}

Member Function Documentation

G4double G4FissionParameters::GetA1 ( void  ) const

inline

Definition at line 62 of file G4FissionParameters.hh.

{ return A1; }

G4double G4FissionParameters::GetA2 ( void  ) const

inline

Definition at line 63 of file G4FissionParameters.hh.

{ return A2; }

G4double G4FissionParameters::GetAs ( void  ) const

inline

Definition at line 65 of file G4FissionParameters.hh.

{ return As; }

G4double G4FissionParameters::GetSigma1 ( void  ) const

inline

Definition at line 66 of file G4FissionParameters.hh.

{ return Sigma1; }

G4double G4FissionParameters::GetSigma2 ( void  ) const

inline

Definition at line 67 of file G4FissionParameters.hh.

{ return Sigma2; }

G4double G4FissionParameters::GetSigmaS ( void  ) const

inline

Definition at line 68 of file G4FissionParameters.hh.

{ return SigmaS; }

G4double G4FissionParameters::GetW ( void  ) const

inline

Definition at line 69 of file G4FissionParameters.hh.

{ return w; }

---

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

• G4FissionParameters.hh
• G4FissionParameters.cc