G4Pow Class Reference

#include <G4Pow.hh>

Collaboration diagram for G4Pow:

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 ()

Private Member Functions
	G4Pow ()
G4double	logBase (G4double x) const

Private Attributes
const G4double	onethird
const G4int	max2
G4double	maxA
G4double	maxA2
G4double	maxAexp
G4DataVector	ener
G4DataVector	logen
G4DataVector	pz13
G4DataVector	lz
G4DataVector	lz2
G4DataVector	fexp
G4DataVector	fact
G4DataVector	logfact

Static Private Attributes
static G4Pow *	fpInstance = 0

Detailed Description

Definition at line 56 of file G4Pow.hh.

Constructor & Destructor Documentation

~G4Pow()

G4Pow::~G4Pow

(

)

Definition at line 123 of file G4Pow.cc.

{}

G4Pow()

G4Pow::G4Pow ( )

private

Definition at line 67 of file G4Pow.cc.

  : onethird(1.0/3.0), max2(5)
{  
#ifdef G4MULTITHREADED
  if(G4Threading::IsWorkerThread())
  { 
    G4Exception ("G4Pow::G4Pow()", "InvalidSetup", FatalException, 
                 "Attempt to instantiate G4Pow in worker thread!");
  }
#endif
  const G4int maxZ = 512; 
  const G4int maxZfact = 170; 

  maxA    = -0.6 + maxZ;
  maxA2   = 1.25 + max2*0.2;
  maxAexp = -0.76+ maxZfact*0.5;

  ener.resize(max2+1,1.0);
  logen.resize(max2+1,0.0);
  lz2.resize(max2+1,0.0);
  pz13.resize(maxZ,0.0);
  lz.resize(maxZ,0.0);
  fexp.resize(maxZfact,0.0);
  fact.resize(maxZfact,0.0);
  logfact.resize(maxZ,0.0);

  G4double f = 1.0;
  G4double logf = 0.0;
  fact[0] = 1.0;
  fexp[0] = 1.0;

  for(G4int i=1; i<=max2; ++i)
  {
    ener[i] = powN(500.,i); 
    logen[i]= G4Log(ener[i]); 
    lz2[i]  = G4Log(1.0 + i*0.2);
  }

  for(G4int i=1; i<maxZ; ++i)
  {
    G4double x  = G4double(i);
    pz13[i] = std::pow(x,onethird);
    lz[i]   = G4Log(x);
    if(i < maxZfact)
    { 
      f *= x; 
      fact[i] = f;
      fexp[i] = G4Exp(0.5*i);
    }
    logf += lz[i];
    logfact[i] = logf;
  }
}

Here is the call graph for this function:

Member Function Documentation

A13()

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:

A23()

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:

expA()

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:

factorial()

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:

GetInstance()

G4Pow * G4Pow::GetInstance ( void  )

static

Definition at line 55 of file G4Pow.cc.

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

log10A()

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:

Here is the caller graph for this function:

log10Z()

G4double G4Pow::log10Z ( G4int  Z ) const

inline

Definition at line 225 of file G4Pow.hh.

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

logA()

G4double G4Pow::logA ( G4double  A ) const

inline

Definition at line 194 of file G4Pow.hh.

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

Here is the call graph for this function:

logBase()

G4double G4Pow::logBase ( G4double  x ) const

inlineprivate

Definition at line 171 of file G4Pow.hh.

{
  G4double res;
  if(a <= maxA2) 
  {
    G4int i = G4int(max2*(a - 1) + 0.5);
    if(i > max2) { i = max2; }
    G4double x = a/(G4double(i)/max2 + 1) - 1;
    res = lz2[i] + x*(1.0 - (0.5 - onethird*x)*x);
  }
  else if(a <= maxA)
  {
    G4int i = G4int(a + 0.5);
    G4double x = a/G4double(i) - 1;
    res = lz[i] + x*(1.0 - (0.5 - onethird*x)*x);
  }
  else
  {
    res = G4Log(a);
  }
  return res;
}

Here is the call graph for this function:

Here is the caller graph for this function:

logfactorial()

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:

logX()

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:

logZ()

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:

powA()

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:

powN()

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:

powZ()

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:

Z13()

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:

Z23()

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:

Member Data Documentation

ener

G4DataVector G4Pow::ener

private

Definition at line 115 of file G4Pow.hh.

fact

G4DataVector G4Pow::fact

private

Definition at line 121 of file G4Pow.hh.

fexp

G4DataVector G4Pow::fexp

private

Definition at line 120 of file G4Pow.hh.

fpInstance

G4Pow * G4Pow::fpInstance = 0

staticprivate

Definition at line 106 of file G4Pow.hh.

logen

G4DataVector G4Pow::logen

private

Definition at line 116 of file G4Pow.hh.

logfact

G4DataVector G4Pow::logfact

private

Definition at line 122 of file G4Pow.hh.

lz

G4DataVector G4Pow::lz

private

Definition at line 118 of file G4Pow.hh.

lz2

G4DataVector G4Pow::lz2

private

Definition at line 119 of file G4Pow.hh.

max2

const G4int G4Pow::max2

private

Definition at line 109 of file G4Pow.hh.

maxA

G4double G4Pow::maxA

private

Definition at line 111 of file G4Pow.hh.

maxA2

G4double G4Pow::maxA2

private

Definition at line 112 of file G4Pow.hh.

maxAexp

G4double G4Pow::maxAexp

private

Definition at line 113 of file G4Pow.hh.

onethird

const G4double G4Pow::onethird

private

Definition at line 108 of file G4Pow.hh.

pz13

G4DataVector G4Pow::pz13

private

Definition at line 117 of file G4Pow.hh.

---

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

• G4Pow.hh
• G4Pow.cc