G4Pow Class Reference

#include <G4Pow.hh>

Public Member Functions
	~G4Pow ()
G4double	Z13 (G4int Z) const
G4double	A13 (G4double A) const
G4double	Z23 (G4int Z) const
G4double	A23 (G4double A) const
G4double	logZ (G4int Z) const
G4double	logA (G4double A) const
G4double	logX (G4double x) const
G4double	log10Z (G4int Z) const
G4double	log10A (G4double A) const
G4double	expA (G4double A) const
G4double	powZ (G4int Z, G4double y) const
G4double	powA (G4double A, G4double y) const
G4double	powN (G4double x, G4int n) const
G4double	factorial (G4int Z) const
G4double	logfactorial (G4int Z) const

Static Public Member Functions
static G4Pow *	GetInstance ()

Detailed Description

Definition at line 56 of file G4Pow.hh.

Constructor & Destructor Documentation

G4Pow::~G4Pow

(

)

Definition at line 123 of file G4Pow.cc.

{}

Member Function Documentation

G4double G4Pow::A13 ( G4double  A ) const

inline

Definition at line 132 of file G4Pow.hh.

{
  G4double res = 0.0;
  if(A > 0.0) 
  {
    G4double a = (1.0 <= A) ? A : 1.0/A;
    if(1.0 > A) { a = 1.0/A; }
    if(a <= maxA)
    {
      G4int i = G4int(a + 0.5);
      G4double x = (a/G4double(i) - 1.0)*onethird;
      res = pz13[i]*(1.0 + x - x*x*(1.0 - 1.66666666*x));
      if(1.0 > A) { res = 1.0/res; }
    }
    else
    {
      res = std::pow(A, onethird); 
    }
  }
  return res;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4Pow::A23 ( G4double  A ) const

inline

Definition at line 160 of file G4Pow.hh.

{
  G4double x = A13(A);
  return x*x;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4Pow::expA ( G4double  A ) const

inline

Definition at line 235 of file G4Pow.hh.

{
  G4double res;
  G4double a = (0.0 <= A) ? A : -A;

  if(a <= maxAexp)
  {
    G4int i = G4int(2*a + 0.5);
    G4double x = a - i*0.5;
    res = fexp[i]*(1.0 + x*(1.0 + 0.5*(1.0 + onethird*x)*x));
  }
  else
  {
    res = G4Exp(a);
  }
  if(0.0 > A) { res = 1.0/res; }
  return res;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4Pow::factorial ( G4int  Z ) const

inline

Definition at line 264 of file G4Pow.hh.

{
  return fact[Z];
}

Here is the caller graph for this function:

G4Pow * G4Pow::GetInstance ( void  )

static

Definition at line 55 of file G4Pow.cc.

{
  if (fpInstance == 0)
  {
    static G4Pow geant4pow;
    fpInstance = &geant4pow;
  }
  return fpInstance;
}

G4double G4Pow::log10A ( G4double  A ) const

inline

Definition at line 230 of file G4Pow.hh.

{
  return logX(A)/lz[10];
}

Here is the call graph for this function:

G4double G4Pow::log10Z ( G4int  Z ) const

inline

Definition at line 225 of file G4Pow.hh.

{
  return lz[Z]/lz[10];
}

G4double G4Pow::logA ( G4double  A ) const

inline

Definition at line 194 of file G4Pow.hh.

{
  return (1.0 <= A ? logBase(A) : -logBase(1./A));
}

G4double G4Pow::logfactorial ( G4int  Z ) const

inline

Definition at line 269 of file G4Pow.hh.

{
  return logfact[Z];
}

Here is the caller graph for this function:

G4double G4Pow::logX ( G4double  x ) const

inline

Definition at line 199 of file G4Pow.hh.

{
  G4double res = 0.0;
  G4double a = (1.0 <= x) ? x : 1.0/x;

  if(a <= maxA) 
  {
    res = logBase(a);
  }
  else if(a <= ener[2])
  {
    res = logen[1] + logBase(a/ener[1]);
  }
  else if(a <= ener[3])
  {
    res = logen[2] + logBase(a/ener[2]);
  }
  else
  {
    res = G4Log(a);
  }

  if(1.0 > x) { res = -res; }
  return res;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4Pow::logZ ( G4int  Z ) const

inline

Definition at line 166 of file G4Pow.hh.

{
  return lz[Z];
}

Here is the caller graph for this function:

G4double G4Pow::powA ( G4double  A, G4double  y  ) const

inline

Definition at line 259 of file G4Pow.hh.

{
  return (0.0 == A ? 0.0 : expA(y*logX(A)));
}

Here is the call graph for this function:

G4double G4Pow::powN	(	G4double	x,
		G4int	n
	)		const

Definition at line 128 of file G4Pow.cc.

{
  if(0.0 == x)        { return 0.0; }
  if(std::abs(n) > 8) { return std::pow(x, G4double(n)); }
  G4double res = 1.0;
  if(n >= 0) { for(G4int i=0; i<n; ++i) { res *= x; } }
  else if(n < 0)
  {
    G4double y = 1.0/x;
    G4int nn = -n;
    for(G4int i=0; i<nn; ++i) { res *= y; }
  }
  return res;
}

Here is the caller graph for this function:

G4double G4Pow::powZ ( G4int  Z, G4double  y  ) const

inline

Definition at line 254 of file G4Pow.hh.

{
  return expA(y*lz[Z]);
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4Pow::Z13 ( G4int  Z ) const

inline

Definition at line 127 of file G4Pow.hh.

{
  return pz13[Z];
}

Here is the caller graph for this function:

G4double G4Pow::Z23 ( G4int  Z ) const

inline

Definition at line 154 of file G4Pow.hh.

{
  G4double x = Z13(Z);
  return x*x;
}

Here is the call graph for this function:

Here is the caller graph for this function:

---

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

• source/geant4.10.03.p03/source/global/management/include/G4Pow.hh
• source/geant4.10.03.p03/source/global/management/src/G4Pow.cc