G4BetaFermiFunction Class Reference

#include <G4BetaFermiFunction.hh>

Public Member Functions
	G4BetaFermiFunction (G4int const fA, G4int const fZ)
	~G4BetaFermiFunction ()
G4double	GetFF (const G4double E)
G4double	GetFFN (const G4double E0)

Detailed Description

Definition at line 31 of file G4BetaFermiFunction.hh.

Constructor & Destructor Documentation

G4BetaFermiFunction::G4BetaFermiFunction ( G4int const  fA, G4int const  fZ  )

inline

Definition at line 39 of file G4BetaFermiFunction.hh.

                                                      :
    A(fA), Z(fZ)
  {};

G4BetaFermiFunction::~G4BetaFermiFunction ( )

inline

Definition at line 46 of file G4BetaFermiFunction.hh.

  {};

Member Function Documentation

G4double G4BetaFermiFunction::GetFF

(

const G4double

E

)

Definition at line 36 of file G4BetaFermiFunction.cc.

{
  G4double A1, A2;
  G4double P, U, S, Y;
  G4double F2;
  G4double E = E0+1.;  
  P=std::sqrt(E*E-1.0) ;
  U=Z/137.0;
  S=std::sqrt(1.0-U*U) - 1.;
  Y = 2*PI*U*E/P;
  A1 = U*U*E*E + P*P/4.;
  A2 = std::fabs(Y/(1-std::exp(-Y)));
  F2 = std::pow(A1,S) * A2; 
  return F2;
}

G4double G4BetaFermiFunction::GetFFN

(

const G4double

E0

)

Definition at line 57 of file G4BetaFermiFunction.cc.

{

  G4double A1, A2;
  G4double P, U, S, Y;
  G4double F2,E;
  G4double EE = E0/100.;
  U=Z/137.0;
  S=std::sqrt(1.0-U*U) - 1.;
  G4double F1 = 1E-10;
  for (G4int i = 1; i<=100 ; i++) {
    E = G4double(i)*EE + 1.;
    P=std::sqrt(E*E-1.0) ;
    Y = 2*PI*U*E/P;
    A1 = U*U*E*E + P*P/4.;
    A2 = std::fabs(Y/(1-std::exp(-Y)));
    F2 = std::pow(A1,S) * A2; 
    if (F2 > F1) F1 = F2;
  }        
  return F1;
}

---

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

• G4BetaFermiFunction.hh
• G4BetaFermiFunction.cc