#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::G4BetaFermiFunction	(	G4int const	fA,
		G4int const	fZ
	)

inline

Definition at line 39 of file G4BetaFermiFunction.hh.

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

◆ ~G4BetaFermiFunction()

G4BetaFermiFunction::~G4BetaFermiFunction ( )

inline

Definition at line 46 of file G4BetaFermiFunction.hh.

49 {};

Member Function Documentation

◆ GetFF()

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;
}

◆ GetFFN()

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: