#include <G4Pow.hh>

Public Member Functions
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.

Member Function Documentation

G4double G4Pow::A13 ( G4double A ) const

inline

Definition at line 134 of file G4Pow.hh.

References test::a, and test::x.

Referenced by A23(), G4AntiNuclElastic::CalculateAm(), G4ionEffectiveCharge::EffectiveCharge(), G4TripathiLightCrossSection::GetElementCrossSection(), G4IonsShenCrossSection::GetIsoCrossSection(), G4NistManager::GetZ13(), and G4BetheHeitlerModel::SampleSecondaries().

 {
   G4double res;
   G4double a = 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;
 }

G4double G4Pow::A23 ( G4double A ) const

inline

Definition at line 159 of file G4Pow.hh.

References A13(), and test::x.

Referenced by G4ionEffectiveCharge::EffectiveCharge().

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

G4double G4Pow::expA ( G4double A ) const

inline

Definition at line 238 of file G4Pow.hh.

References G4Exp(), and test::x.

Referenced by powA(), and powZ().

 {
   G4double res;
   G4double a = A;
   if(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;
 }

G4double G4Pow::factorial ( G4int Z ) const

inline

Definition at line 268 of file G4Pow.hh.

Referenced by G4PomeronCrossSection::GetCutPomeronProbability().

 {
   return fact[Z];
 }

G4Pow * G4Pow::GetInstance ( void )

static

Definition at line 53 of file G4Pow.cc.

Referenced by G4HadPhaseSpaceKopylov::BetaKopylov(), G4CascadeFinalStateAlgorithm::BetaKopylov(), G4GEMCoulombBarrier::CalcCompoundRadius(), G4GEMCoulombBarrierHE::CalcCompoundRadius(), G4AntiNuclElastic::CalculateAm(), G4UrbanMscModel::ComputeCrossSectionPerAtom(), G4E1SingleProbability1::EmissionProbDensity(), G4BetheHeitlerModel::G4BetheHeitlerModel(), G4BGGPionInelasticXS::G4BGGPionInelasticXS(), G4ContinuumGammaTransition::G4ContinuumGammaTransition(), G4E1Probability::G4E1Probability(), G4ExcitationHandler::G4ExcitationHandler(), G4FermiConfigurationList::G4FermiConfigurationList(), G4FermiPhaseSpaceDecay::G4FermiPhaseSpaceDecay(), G4GEMChannel::G4GEMChannel(), G4GEMProbability::G4GEMProbability(), G4Generator2BS::G4Generator2BS(), G4HadronCrossSections::G4HadronCrossSections(), G4ICRU49NuclearStoppingModel::G4ICRU49NuclearStoppingModel(), G4ionEffectiveCharge::G4ionEffectiveCharge(), G4IonFluctuations::G4IonFluctuations(), G4IonisParamElm::G4IonisParamElm(), G4IonisParamMat::G4IonisParamMat(), G4PreCompoundEmission::G4PreCompoundEmission(), G4PreCompoundTransitions::G4PreCompoundTransitions(), G4VEmissionProbability::G4VEmissionProbability(), G4VPreCompoundFragment::G4VPreCompoundFragment(), G4WentzelOKandVIxSection::G4WentzelOKandVIxSection(), G4WentzelVIRelModel::G4WentzelVIRelModel(), G4WentzelVIRelXSection::G4WentzelVIRelXSection(), G4AntiNuclElastic::GetcosTeta1(), G4CoulombBarrier::GetCoulombBarrier(), G4NeutronInelasticCrossSection::GetCrossSection(), G4PomeronCrossSection::GetCutPomeronProbability(), G4IonsSihverCrossSection::GetElementCrossSection(), G4TripathiLightCrossSection::GetElementCrossSection(), G4IonsShenCrossSection::GetIsoCrossSection(), G4InuclParamMomDst::GetMomentum(), G4InuclSpecialFunctions::randomInuclPowers(), G4AntiNuclElastic::SampleInvariantT(), G4HadronElastic::SampleInvariantT(), and G4DiscreteGammaTransition::SelectGamma().

 {
   if (fpInstance == 0)
   {
     fpInstance = new G4Pow;
   }
   return fpInstance;
 }

G4double G4Pow::log10A ( G4double A ) const

inline

Definition at line 233 of file G4Pow.hh.

References logX().

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

G4double G4Pow::log10Z ( G4int Z ) const

inline

Definition at line 228 of file G4Pow.hh.

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

G4double G4Pow::logA ( G4double A ) const

inline

Definition at line 193 of file G4Pow.hh.

Referenced by G4NistManager::GetLOGA().

 {
   G4double res;
   if(1.0 <= A) { res = logBase(A); }
   else         { res = -logBase(1./A); }
   return res;
 }

G4double G4Pow::logfactorial ( G4int Z ) const

inline

Definition at line 273 of file G4Pow.hh.

 {
   return logfact[Z];
 }

G4double G4Pow::logX ( G4double x ) const

inline

Definition at line 201 of file G4Pow.hh.

References G4Log(), and test::x.

Referenced by log10A(), and powA().

 {
   G4double res = 0.0;
   G4double a = x;
   if(1.0 > x) { a = 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;
 }

G4double G4Pow::logZ ( G4int Z ) const

inline

Definition at line 165 of file G4Pow.hh.

Referenced by G4GEMProbability::G4GEMProbability(), G4IonisParamElm::G4IonisParamElm(), G4IonisParamMat::G4IonisParamMat(), G4NeutronInelasticCrossSection::GetCrossSection(), G4NistManager::GetLOGZ(), G4PreCompoundProton::GetOpt2(), and G4DiscreteGammaTransition::SelectGamma().

 {
   return lz[Z];
 }

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

inline

Definition at line 263 of file G4Pow.hh.

References expA(), and logX().

Referenced by G4NonEquilibriumEvaporator::deExcite().

 {
   return expA(y*logX(A));
 }

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

Definition at line 125 of file G4Pow.cc.

References n, G4InuclParticleNames::nn, and test::x.

Referenced by G4HadPhaseSpaceKopylov::BetaKopylov(), G4CascadeFinalStateAlgorithm::BetaKopylov(), G4PreCompoundTransitions::CalculateProbability(), G4NonEquilibriumEvaporator::deExcite(), G4FermiConfigurationList::G4FermiConfigurationList(), G4InuclParamMomDst::GetMomentum(), G4PreCompoundNucleon::ProbabilityDistributionFunction(), G4PreCompoundIon::ProbabilityDistributionFunction(), and G4InuclSpecialFunctions::randomInuclPowers().

 {
   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) && (x != 0.0))
   {
     G4double y = 1.0/x;
     G4int nn = -n;
     for(G4int i=0; i<nn; ++i) { res *= y; }
   }
   return res;
 }

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

inline

Definition at line 258 of file G4Pow.hh.

References expA().

Referenced by G4E1SingleProbability1::EmissionProbDensity(), G4E1Probability::EmissionProbDensity(), G4ContinuumGammaTransition::G4ContinuumGammaTransition(), G4PreCompoundProton::GetOpt1(), G4PreCompoundHe3::GetOpt12(), G4PreCompoundAlpha::GetOpt12(), G4PreCompoundTriton::GetOpt12(), G4PreCompoundDeuteron::GetOpt12(), G4PreCompoundProton::GetOpt3(), G4PreCompoundHe3::GetOpt34(), G4PreCompoundTriton::GetOpt34(), G4PreCompoundAlpha::GetOpt34(), G4PreCompoundDeuteron::GetOpt34(), and G4HadronElastic::SampleInvariantT().

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

G4double G4Pow::Z13 ( G4int Z ) const

inline

Definition at line 129 of file G4Pow.hh.

Referenced by G4GEMProbability::CalcAlphaParam(), G4GEMCoulombBarrier::CalcCompoundRadius(), G4GEMCoulombBarrierHE::CalcCompoundRadius(), G4PreCompoundAlpha::CrossSection(), G4PreCompoundDeuteron::CrossSection(), G4PreCompoundHe3::CrossSection(), G4PreCompoundNeutron::CrossSection(), G4PreCompoundProton::CrossSection(), G4PreCompoundTriton::CrossSection(), G4FermiConfigurationList::G4FermiConfigurationList(), G4IonisParamElm::G4IonisParamElm(), G4WentzelOKandVIxSection::G4WentzelOKandVIxSection(), G4WentzelVIRelXSection::G4WentzelVIRelXSection(), G4HETCNeutron::GetAlpha(), G4CoulombBarrier::GetCoulombBarrier(), G4NeutronInelasticCrossSection::GetCrossSection(), G4IonsSihverCrossSection::GetElementCrossSection(), G4TripathiLightCrossSection::GetElementCrossSection(), G4HadronCrossSections::GetFissionCrossSection(), G4IonsShenCrossSection::GetIsoCrossSection(), G4NistManager::GetZ13(), G4VPreCompoundFragment::Initialize(), G4Generator2BS::SampleDirection(), G4AntiNuclElastic::SampleInvariantT(), G4HadronElastic::SampleInvariantT(), and Z23().

 {
   return pz13[Z];
 }

G4double G4Pow::Z23 ( G4int Z ) const

inline

Definition at line 153 of file G4Pow.hh.

References test::x, and Z13().

Referenced by G4GEMProbability::CalcBetaParam(), G4UrbanMscModel::ComputeCrossSectionPerAtom(), G4HETCNeutron::GetBeta(), G4AntiNuclElastic::GetcosTeta1(), G4HadronElastic::SampleInvariantT(), and G4WentzelOKandVIxSection::SetupTarget().

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

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

Public Member Functions

Static Public Member Functions

Detailed Description

Member Function Documentation