G4Pow Class Reference

#include <G4Pow.hh>

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

Static Public Member Functions
static G4Pow *	GetInstance ()

Private Member Functions
G4double	A13High (const G4double, const G4bool) const
G4double	A13Low (const G4double, const G4bool) const
	G4Pow ()
G4double	logBase (G4double x) const

Private Attributes
G4DataVector	ener
G4DataVector	fact
G4DataVector	fexp
G4DataVector	logen
G4DataVector	logfact
G4DataVector	lowa13
G4DataVector	lz
G4DataVector	lz2
const G4int	max2 = 5
G4double	maxA
G4double	maxA2
G4double	maxAexp
G4double	maxLowA
const G4double	onethird = 1.0 / 3.0
G4DataVector	pz13

Static Private Attributes
static G4Pow *	fpInstance = nullptr

Detailed Description

Definition at line 48 of file G4Pow.hh.

Constructor & Destructor Documentation

~G4Pow()

G4Pow::~G4Pow

(

)

Definition at line 116 of file G4Pow.cc.

{}

G4Pow()

G4Pow::G4Pow ( )

private

Definition at line 53 of file G4Pow.cc.

{
#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;
  const G4int numLowA  = 17;
 
  maxA    = -0.6 + maxZ;
  maxLowA = 4.0;
  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);
  lowa13.resize(numLowA, 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;
  }
 
  for(G4int i = 4; i < numLowA; ++i)
  {
    lowa13[i] = std::pow(0.25 * i, onethird);
  }
}

References ener, fact, FatalException, fexp, G4Exception(), G4Exp(), G4Log(), G4Threading::IsWorkerThread(), logen, logfact, lowa13, lz, lz2, max2, maxA, maxA2, maxAexp, maxLowA, maxZ, onethird, powN(), and pz13.

Member Function Documentation

A13()

G4double G4Pow::A13

(

G4double

A

)

const

Definition at line 120 of file G4Pow.cc.

{
  G4double res = 0.;
  if(A > 0.)
  {
    const bool invert = (A < 1.);
    const G4double a  = invert ? 1. / A : A;
    res               = (a < maxLowA) ? A13Low(a, invert) : A13High(a, invert);
  }
  return res;
}

References A, A13High(), A13Low(), and maxLowA.

Referenced by A23(), G4LFission::Atomas(), G4AntiNuclElastic::CalculateAm(), G4DiffuseElastic::CalculateAm(), G4DiffuseElasticV2::CalculateAm(), G4NuclNuclDiffuseElastic::CalculateAm(), G4DiffuseElastic::CalculateNuclearRad(), G4DiffuseElasticV2::CalculateNuclearRad(), G4NuclNuclDiffuseElastic::CalculateNuclearRad(), G4StatMFMacroCanonical::CalculateTemperature(), G4FermiMomentum::cbrt(), G4QMDMeanField::DoClusterJudgment(), G4KM_NucleonEqRhs::EvaluateRhsGivenB(), G4IonFluctuations::Factor(), G4StatMFMacroTemperature::FragsExcitEnergy(), G4QMDGroundStateNucleus::G4QMDGroundStateNucleus(), G4QMDParameters::G4QMDParameters(), G4fissionEvent::G4SmpNugDist(), G4WilsonAbrasionModel::GetAbradedNucleons(), G4SigmaMinusField::GetBarrier(), G4SigmaPlusField::GetBarrier(), G4StatMFMicroPartition::GetCoulombEnergy(), G4ComponentAntiNuclNuclearXS::GetInelasticElementCrossSection(), G4IonsShenCrossSection::GetIsoCrossSection(), G4QMDNucleus::GetNuclearMass(), G4StatMFMicroPartition::GetPartitionEnergy(), G4ComponentAntiNuclNuclearXS::GetTotalElementCrossSection(), G4WilsonRadius::GetWilsonRMSRadius(), G4NistManager::GetZ13(), G4mplIonisationModel::Initialise(), G4mplIonisationWithDeltaModel::Initialise(), MCGIDI_KalbachMann_S_a_or_b(), G4StatMFChannel::PlaceFragments(), G4DipBustGenerator::SampleCosTheta(), G4BetheHeitlerModel::SampleSecondaries(), G4PairProductionRelModel::SampleSecondaries(), G4WentzelOKandVIxSection::SampleSingleScattering(), G4QMDGroundStateNucleus::samplingMomentum(), G4ParticleHPKallbachMannSyst::SeparationEnergy(), and G4StatMFChannel::SolveEqOfMotion().

A13High()

G4double G4Pow::A13High ( const G4double  a, const  G4bool  ) const

private

Definition at line 134 of file G4Pow.cc.

{
  G4double res;
  if(a < maxA)
  {
    const G4int i    = static_cast<G4int>(a + 0.5);
    const G4double x = (a / i - 1.) * onethird;
    res              = pz13[i] * (1. + x - x * x * (1. - 1.666667 * x));
  }
  else
  {
    res = G4Exp(G4Log(a) * onethird);
  }
  res = invert ? 1. / res : res;
  return res;
}

References G4Exp(), G4Log(), maxA, onethird, and pz13.

Referenced by A13().

A13Low()

G4double G4Pow::A13Low ( const G4double  a, const G4bool  invert  ) const

private

Definition at line 153 of file G4Pow.cc.

{
  G4double res;
  const G4int i    = static_cast<G4int>(4. * (a + 0.125));
  const G4double y = 0.25 * i;
  const G4double x = (a / y - 1.) * onethird;
  res              = lowa13[i] * (1. + x - x * x * (1. - 1.666667 * x));
  res              = invert ? 1. / res : res;
  return res;
}

References lowa13, and onethird.

Referenced by A13().

A23()

G4double G4Pow::A23 ( G4double  A ) const

inline

Definition at line 131 of file G4Pow.hh.

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

References A, and A13().

Referenced by G4LFission::Atomas(), G4AtimaEnergyLossModel::Bethek_dedx_e(), G4DiffuseElastic::CalculateNuclearRad(), G4DiffuseElasticV2::CalculateNuclearRad(), G4RKFieldIntegrator::CalculateTotalEnergy(), G4IonisParamMat::ComputeIonParameters(), G4QMDMeanField::DoClusterJudgment(), G4ionEffectiveCharge::EffectiveCharge(), G4KM_NucleonEqRhs::G4KM_NucleonEqRhs(), G4NuclearAbrasionGeometry::GetExcitationEnergyOfProjectile(), G4NuclearAbrasionGeometry::GetExcitationEnergyOfTarget(), G4QMDNucleus::GetNuclearMass(), G4ParticleHPJENDLHEData::getXSfromThisIsotope(), G4ComponentBarNucleonNucleusXsc::LoadData(), and G4ParticleHPKallbachMannSyst::SeparationEnergy().

expA()

G4double G4Pow::expA ( G4double  A ) const

inline

Definition at line 203 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;
}

References A, fexp, G4Exp(), maxAexp, and onethird.

Referenced by G4IonFluctuations::Factor(), powA(), and powZ().

factorial()

G4double G4Pow::factorial ( G4int  Z ) const

inline

Definition at line 235 of file G4Pow.hh.

{ return fact[Z]; }

References fact, and Z.

Referenced by G4LegendrePolynomial::EvalAssocLegendrePoly().

GetInstance()

G4Pow * G4Pow::GetInstance ( )

static

Definition at line 41 of file G4Pow.cc.

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

References fpInstance.

Referenced by G4ParticleHPKallbachMannSyst::A(), G4FTFAnnihilation::Annihilate(), G4WilsonAbrasionModel::ApplyYourself(), G4LENDElastic::ApplyYourself(), G4QMDReaction::ApplyYourself(), G4ParticleHPElasticFS::ApplyYourself(), G4LFission::Atomas(), G4FissionBarrier::BarashenkovFissionBarrier(), G4CascadeFinalStateAlgorithm::BetaKopylov(), G4HadPhaseSpaceKopylov::BetaKopylov(), G4WilsonAblationModel::BreakItUp(), G4UPiNuclearCrossSection::BuildPhysicsTable(), G4ParticleHPNBodyPhaseSpace::C(), G4QMDMeanField::Cal2BodyQuantities(), G4StatMFMacroMultiplicity::CalcChemicalPotentialMu(), G4StatMFMacroChemicalPotential::CalcChemicalPotentialNu(), G4QuasiElRatios::CalcElTot(), G4StatMFMacroBiNucleon::CalcEnergy(), G4StatMFMacroMultiNucleon::CalcEnergy(), G4StatMFMacroTetraNucleon::CalcEnergy(), G4StatMFMacroTriNucleon::CalcEnergy(), G4StatMFMacroMultiNucleon::CalcEntropy(), G4StatMFMicroCanonical::CalcEntropyOfCompoundNucleus(), G4StatMFFragment::CalcExcitationEnergy(), G4StatMFMicroCanonical::CalcFreeInternalEnergy(), G4StatMFMacroBiNucleon::CalcMeanMultiplicity(), G4StatMFMacroMultiNucleon::CalcMeanMultiplicity(), G4StatMFMacroTetraNucleon::CalcMeanMultiplicity(), G4StatMFMacroTriNucleon::CalcMeanMultiplicity(), G4StatMFMicroPartition::CalcPartitionProbability(), G4QuasiElRatios::CalcQF2IN_Ratio(), G4AntiNuclElastic::CalculateAm(), G4DiffuseElastic::CalculateAm(), G4DiffuseElasticV2::CalculateAm(), G4NuclNuclDiffuseElastic::CalculateAm(), G4MuonVDNuclearModel::CalculateEMVertex(), G4DiffuseElastic::CalculateNuclearRad(), G4DiffuseElasticV2::CalculateNuclearRad(), G4NuclNuclDiffuseElastic::CalculateNuclearRad(), G4StatMFMacroCanonical::CalculateTemperature(), G4RKFieldIntegrator::CalculateTotalEnergy(), G4StatMFMacroMultiNucleon::CalcZARatio(), G4QMDNucleus::CalEnergyAndAngularMomentumInCM(), G4QMDCollision::CalFinalStateOfTheBinaryCollisionJQMD(), G4QMDMeanField::CalGraduate(), G4FermiMomentum::cbrt(), G4Reggeons::Chi_pomeron(), G4Reggeons::Chi_reggeon(), G4DiffractiveExcitation::ChooseP(), G4QGSDiffractiveExcitation::ChooseP(), G4SingleDiffractiveExcitation::ChooseX(), G4StatMFMacroCanonical::ChooseZ(), G4Clebsch::ClebschGordanCoeff(), G4UrbanAdjointMscModel::ComputeCrossSectionPerAtom(), G4UrbanMscModel::ComputeCrossSectionPerAtom(), G4GoudsmitSaundersonMscModel::ComputeGeomPathLength(), G4IonisParamMat::ComputeIonParameters(), G4KalbachCrossSection::ComputePowerParameter(), G4GoudsmitSaundersonMscModel::ComputeTrueStepLength(), G4RToEConvForElectron::ComputeValue(), G4RToEConvForGamma::ComputeValue(), G4RToEConvForPositron::ComputeValue(), G4Fancy3DNucleus::CoulombBarrier(), G4StatMFMicroPartition::CoulombFreeEnergy(), G4DiffractiveExcitation::CreateStrings(), G4XAqmElastic::CrossSection(), G4XAqmTotal::CrossSection(), G4XPDGElastic::CrossSection(), G4XPDGTotal::CrossSection(), G4ChipsAntiBaryonInelasticXS::CrossSectionFormula(), G4QMDMeanField::DoClusterJudgment(), G4LENDCrossSection::DumpPhysicsTable(), G4ParticleHPCaptureData::DumpPhysicsTable(), G4ParticleHPElasticData::DumpPhysicsTable(), G4ParticleHPFissionData::DumpPhysicsTable(), G4ParticleHPInelasticData::DumpPhysicsTable(), G4LegendrePolynomial::EvalAssocLegendrePoly(), G4KM_NucleonEqRhs::EvaluateRhsGivenB(), G4FPYSamplingOps::EvaluateWattConstants(), G4QuarkExchange::ExciteParticipants(), G4NuclearAbrasionGeometry::F(), G4VCrossSectionSource::FcrossX(), G4InuclSpecialFunctions::FermiEnergy(), G4JTPolynomialSolver::FindRoots(), G4CompetitiveFission::FissionKineticEnergy(), G4StatMFMacroTemperature::FragsExcitEnergy(), G4AtimaEnergyLossModel::G4AtimaEnergyLossModel(), G4AtimaFluctuations::G4AtimaFluctuations(), G4BGGPionElasticXS::G4BGGPionElasticXS(), G4BGGPionInelasticXS::G4BGGPionInelasticXS(), G4InuclSpecialFunctions::G4cbrt(), G4EmCorrections::G4EmCorrections(), G4ExcitationHandler::G4ExcitationHandler(), G4FermiPhaseSpaceDecay::G4FermiPhaseSpaceDecay(), G4GEMChannel::G4GEMChannel(), G4GEMProbability::G4GEMProbability(), G4Generator2BS::G4Generator2BS(), G4FissionProductYieldDist::G4GetFission(), G4HadronNucleonXsc::G4HadronNucleonXsc(), G4ICRU49NuclearStoppingModel::G4ICRU49NuclearStoppingModel(), G4IonCoulombCrossSection::G4IonCoulombCrossSection(), G4ionEffectiveCharge::G4ionEffectiveCharge(), G4IonFluctuations::G4IonFluctuations(), G4IonisParamElm::G4IonisParamElm(), G4IonisParamMat::G4IonisParamMat(), G4KM_NucleonEqRhs::G4KM_NucleonEqRhs(), G4LevelManager::G4LevelManager(), G4LevelReader::G4LevelReader(), G4LindhardSorensenData::G4LindhardSorensenData(), G4NistManager::G4NistManager(), G4NuclearFermiDensity::G4NuclearFermiDensity(), G4NuclearLevelData::G4NuclearLevelData(), G4NuclearShellModelDensity::G4NuclearShellModelDensity(), G4PreCompoundEmission::G4PreCompoundEmission(), G4QMDGroundStateNucleus::G4QMDGroundStateNucleus(), G4QMDParameters::G4QMDParameters(), G4FPYSamplingOps::G4SampleWatt(), G4ScreeningMottCrossSection::G4ScreeningMottCrossSection(), G4fissionEvent::G4SmpGEng(), G4fissionEvent::G4SmpNEngCf252(), G4fissionEvent::G4SmpNuDistDataPu239(), G4fissionEvent::G4SmpNuDistDataPu239_241(), G4fissionEvent::G4SmpNuDistDataU232_234_236_238(), G4fissionEvent::G4SmpNuDistDataU233_235(), G4fissionEvent::G4SmpNuDistDataU235(), G4fissionEvent::G4SmpNuDistDataU238(), G4fissionEvent::G4SmpNugDist(), G4VEmissionProbability::G4VEmissionProbability(), G4VPreCompoundFragment::G4VPreCompoundFragment(), G4WentzelOKandVIxSection::G4WentzelOKandVIxSection(), G4ParticleHPMadlandNixSpectrum::Gamma25(), G4QMDParticipant::Get4Momentum(), G4WilsonAbrasionModel::GetAbradedNucleons(), G4AntiProtonField::GetBarrier(), G4KaonMinusField::GetBarrier(), G4KaonPlusField::GetBarrier(), G4PionMinusField::GetBarrier(), G4PionPlusField::GetBarrier(), G4ProtonField::GetBarrier(), G4SigmaMinusField::GetBarrier(), G4SigmaPlusField::GetBarrier(), G4ParticleHPInterpolator::GetBinIntegral(), G4EMDissociationSpectrum::GetClosestApproach(), G4AntiNuclElastic::GetcosTeta1(), G4StatMFMicroPartition::GetCoulombEnergy(), G4EMDissociationCrossSection::GetCrossSectionForProjectile(), G4ChargeExchangeProcess::GetElementCrossSection(), G4StatMFFragment::GetEnergy(), G4ChipsAntiBaryonElasticXS::GetExchangeT(), G4ChipsHyperonElasticXS::GetExchangeT(), G4ChipsPionMinusElasticXS::GetExchangeT(), G4ChipsPionPlusElasticXS::GetExchangeT(), G4NuclearAbrasionGeometry::GetExcitationEnergyOfProjectile(), G4NuclearAbrasionGeometry::GetExcitationEnergyOfTarget(), G4EMDissociationSpectrum::GetGeneralE2Spectrum(), G4NuclNuclDiffuseElastic::GetHadronNucleonXscNS(), G4IonsShenCrossSection::GetIsoCrossSection(), G4LundStringFragmentation::GetLightConeZ(), G4QGSMFragmentation::GetLightConeZ(), G4InuclParamMomDst::GetMomentum(), G4QMDNucleus::GetNuclearMass(), G4StatMFMicroPartition::GetPartitionEnergy(), G4QMDMeanField::GetPotential(), G4ChipsAntiBaryonElasticXS::GetPTables(), G4ChipsHyperonElasticXS::GetPTables(), G4ChipsPionMinusElasticXS::GetPTables(), G4ChipsPionPlusElasticXS::GetPTables(), G4ChipsAntiBaryonElasticXS::GetTabValues(), G4ChipsHyperonElasticXS::GetTabValues(), G4ChipsPionMinusElasticXS::GetTabValues(), G4ChipsPionPlusElasticXS::GetTabValues(), G4QMDMeanField::GetTotalPotential(), G4ParticleHPInterpolator::GetWeightedBinIntegral(), G4WilsonRadius::GetWilsonRMSRadius(), G4ParticleHPJENDLHEData::getXSfromThisIsotope(), G4ParticleHPMadlandNixSpectrum::GIntegral(), G4Bessel::I0(), G4Bessel::I1(), G4FTFParameters::InitForInteraction(), G4mplIonisationModel::Initialise(), G4mplIonisationWithDeltaModel::Initialise(), G4PhotonEvaporation::InitialiseGRData(), G4StatMFMacroCanonical::Initialize(), G4StatMFMicroCanonical::Initialize(), G4ParticleHPVector::Integrate(), G4PolarizationTransition::LnFactorial(), G4ComponentBarNucleonNucleusXsc::LoadData(), G4ParticleHPMadlandNixSpectrum::Madland(), MCGIDI_energy_NBodyPhaseSpacePDF_callback(), MCGIDI_KalbachMann_S_a_or_b(), MCGIDI_sampling_interpolationValues(), MCGIDI_sampling_sampleX_from_pdfsOfXGivenW(), G4NuclearAbrasionGeometry::P(), G4QMDGroundStateNucleus::packNucleons(), G4StatMFChannel::PlaceFragments(), G4ParticleHPNBodyPhaseSpace::Prob(), ptwXY_f_integrate(), ptwXY_LogLogToLinLin(), ptwXY_pow_callback(), G4LEHadronProtonElastic::RandCosThetaDipPen(), G4InuclSpecialFunctions::randomInuclPowers(), G4DipBustGenerator::SampleCosTheta(), G4AntiNuclElastic::SampleInvariantT(), G4HadronElastic::SampleInvariantT(), G4LivermorePolarizedComptonModel::SampleSecondaries(), G4ChargeExchange::SampleT(), G4QGSMSplitableHadron::SampleX(), G4QGSParticipants::SampleX(), G4QMDGroundStateNucleus::samplingMomentum(), G4ParticleHPKallbachMannSyst::SeparationEnergy(), G4FissionProductYieldDist::SetNubar(), G4StatMFChannel::SolveEqOfMotion(), G4LundStringFragmentation::Splitup(), stirf(), G4ChipsPionPlusInelasticXS::ThresholdMomentum(), G4ChipsProtonInelasticXS::ThresholdMomentum(), G4Clebsch::TriangleCoeff(), G4Clebsch::Wigner6J(), G4Clebsch::WignerLittleD(), and G4BigBanger::xProbability().

log10A()

G4double G4Pow::log10A ( G4double  A ) const

inline

Definition at line 201 of file G4Pow.hh.

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

References A, logX(), and lz.

log10Z()

G4double G4Pow::log10Z ( G4int  Z ) const

inline

Definition at line 199 of file G4Pow.hh.

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

References lz, and Z.

logA()

G4double G4Pow::logA ( G4double  A ) const

inline

Definition at line 165 of file G4Pow.hh.

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

References A, and logBase().

logBase()

G4double G4Pow::logBase ( G4double  x ) const

inlineprivate

Definition at line 139 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;
}

References G4Log(), lz, lz2, max2, maxA, maxA2, and onethird.

Referenced by logA(), and logX().

logfactorial()

G4double G4Pow::logfactorial ( G4int  Z ) const

inline

Definition at line 237 of file G4Pow.hh.

{ return logfact[Z]; }

References logfact, and Z.

Referenced by G4Clebsch::ClebschGordanCoeff(), G4LegendrePolynomial::EvalAssocLegendrePoly(), G4PolarizationTransition::LnFactorial(), G4PreCompoundEmission::rho(), G4Clebsch::TriangleCoeff(), G4Clebsch::Wigner6J(), and G4Clebsch::WignerLittleD().

logX()

G4double G4Pow::logX ( G4double  x ) const

inline

Definition at line 170 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;
}

References ener, G4Log(), logBase(), logen, and maxA.

Referenced by log10A(), and powA().

logZ()

G4double G4Pow::logZ ( G4int  Z ) const

inline

Definition at line 137 of file G4Pow.hh.

{ return lz[Z]; }

References lz, and Z.

Referenced by G4StatMFMicroPartition::CalcPartitionProbability(), G4MuonVDNuclearModel::CalculateEMVertex(), G4RToEConvForElectron::ComputeValue(), G4RToEConvForGamma::ComputeValue(), G4RToEConvForPositron::ComputeValue(), G4ScreeningMottCrossSection::FormFactor2ExpHof(), G4ScreeningMottCrossSection::FormFactor2Gauss(), G4GEMProbability::G4GEMProbability(), G4IonisParamElm::G4IonisParamElm(), G4IonisParamMat::G4IonisParamMat(), G4LevelReader::G4LevelReader(), and G4NistManager::GetLOGZ().

powA()

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

inline

Definition at line 230 of file G4Pow.hh.

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

References A, expA(), and logX().

Referenced by G4FTFAnnihilation::Annihilate(), G4WilsonAbrasionModel::ApplyYourself(), G4LENDElastic::ApplyYourself(), G4ParticleHPElasticFS::ApplyYourself(), G4LFission::Atomas(), G4AtimaEnergyLossModel::Bethek_dedx_e(), G4WilsonAblationModel::BreakItUp(), G4UPiNuclearCrossSection::BuildPhysicsTable(), G4ParticleHPNBodyPhaseSpace::C(), G4QuasiElRatios::CalcElTot(), G4QuasiElRatios::CalcQF2IN_Ratio(), G4DiffuseElastic::CalculateNuclearRad(), G4DiffuseElasticV2::CalculateNuclearRad(), G4RKFieldIntegrator::CalculateTotalEnergy(), G4QMDMeanField::CalGraduate(), G4Reggeons::Chi_pomeron(), G4Reggeons::Chi_reggeon(), G4DiffractiveExcitation::ChooseP(), G4QGSDiffractiveExcitation::ChooseP(), G4SingleDiffractiveExcitation::ChooseX(), G4GoudsmitSaundersonMscModel::ComputeGeomPathLength(), G4GoudsmitSaundersonMscModel::ComputeTrueStepLength(), G4XAqmElastic::CrossSection(), G4XPDGElastic::CrossSection(), G4XPDGTotal::CrossSection(), G4ChipsAntiBaryonInelasticXS::CrossSectionFormula(), G4AtimaEnergyLossModel::dedx_n(), G4NonEquilibriumEvaporator::deExcite(), G4AtimaFluctuations::Dispersion(), G4LENDCrossSection::DumpPhysicsTable(), G4ParticleHPCaptureData::DumpPhysicsTable(), G4ParticleHPElasticData::DumpPhysicsTable(), G4ParticleHPFissionData::DumpPhysicsTable(), G4ParticleHPInelasticData::DumpPhysicsTable(), G4FPYSamplingOps::EvaluateWattConstants(), G4QuarkExchange::ExciteParticipants(), G4NuclearAbrasionGeometry::F(), G4IonFluctuations::Factor(), G4VCrossSectionSource::FcrossX(), G4QMDParameters::G4QMDParameters(), G4fissionEvent::G4SmpGEng(), G4fissionEvent::G4SmpNEngCf252(), G4fissionEvent::G4SmpNuDistDataU233_235(), G4fissionEvent::G4SmpNugDist(), G4ParticleHPMadlandNixSpectrum::Gamma25(), G4ParticleHPInterpolator::GetBinIntegral(), G4EMDissociationSpectrum::GetClosestApproach(), G4EMDissociationCrossSection::GetCrossSectionForProjectile(), G4ChargeExchangeProcess::GetElementCrossSection(), G4ChipsAntiBaryonElasticXS::GetExchangeT(), G4ChipsHyperonElasticXS::GetExchangeT(), G4ChipsPionMinusElasticXS::GetExchangeT(), G4ChipsPionPlusElasticXS::GetExchangeT(), G4EMDissociationSpectrum::GetGeneralE2Spectrum(), G4NuclNuclDiffuseElastic::GetHadronNucleonXscNS(), G4ComponentAntiNuclNuclearXS::GetInelasticElementCrossSection(), G4LundStringFragmentation::GetLightConeZ(), G4QGSMFragmentation::GetLightConeZ(), G4QMDMeanField::GetPotential(), G4ChipsAntiBaryonElasticXS::GetPTables(), G4ChipsHyperonElasticXS::GetPTables(), G4ChipsPionMinusElasticXS::GetPTables(), G4ChipsPionPlusElasticXS::GetPTables(), G4ChipsAntiBaryonElasticXS::GetTabValues(), G4ChipsHyperonElasticXS::GetTabValues(), G4ChipsPionMinusElasticXS::GetTabValues(), G4ChipsPionPlusElasticXS::GetTabValues(), G4ComponentAntiNuclNuclearXS::GetTotalElementCrossSection(), G4QMDMeanField::GetTotalPotential(), G4ParticleHPInterpolator::GetWeightedBinIntegral(), G4ParticleHPMadlandNixSpectrum::GIntegral(), G4FTFParameters::InitForInteraction(), G4ParticleHPVector::Integrate(), G4ParticleHPMadlandNixSpectrum::Madland(), MCGIDI_energy_NBodyPhaseSpacePDF_callback(), MCGIDI_sampling_interpolationValues(), MCGIDI_sampling_sampleX_from_pdfsOfXGivenW(), G4QMDGroundStateNucleus::packNucleons(), G4ParticleHPNBodyPhaseSpace::Prob(), ptwXY_f_integrate(), ptwXY_LogLogToLinLin(), ptwXY_pow_callback(), G4LEHadronProtonElastic::RandCosThetaDipPen(), G4QGSMSplitableHadron::SampleX(), G4QGSParticipants::SampleX(), G4ParticleHPKallbachMannSyst::SeparationEnergy(), G4FissionProductYieldDist::SetNubar(), G4AtimaEnergyLossModel::sezi_dedx_e(), G4AtimaEnergyLossModel::sezi_p_se(), G4LundStringFragmentation::Splitup(), stirf(), G4ChipsPionPlusInelasticXS::ThresholdMomentum(), and G4ChipsProtonInelasticXS::ThresholdMomentum().

powN()

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

Definition at line 166 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;
}

References CLHEP::detail::n, and G4InuclParticleNames::nn.

Referenced by G4ParticleHPKallbachMannSyst::A(), G4QMDReaction::ApplyYourself(), G4CascadeFinalStateAlgorithm::BetaKopylov(), G4HadPhaseSpaceKopylov::BetaKopylov(), G4FermiPhaseSpaceDecay::BetaKopylov(), G4ParticleHPNBodyPhaseSpace::C(), G4QMDMeanField::Cal2BodyQuantities(), G4QMDCollision::CalFinalStateOfTheBinaryCollisionJQMD(), G4XAqmTotal::CrossSection(), G4NonEquilibriumEvaporator::deExcite(), G4DensityEffectCalculator::DEll(), G4DensityEffectCalculator::DeltaOnceSolved(), G4DensityEffectCalculator::DFRho(), G4DensityEffectCalculator::Ell(), G4NuclearAbrasionGeometry::F(), G4DensityEffectCalculator::FermiDeltaCalculation(), G4JTPolynomialSolver::FindRoots(), G4ScreeningMottCrossSection::FormFactor2UniformHelm(), G4DensityEffectCalculator::FRho(), G4Pow(), G4fissionEvent::G4SmpNEngCf252(), G4fissionEvent::G4SmpNuDistDataPu239(), G4fissionEvent::G4SmpNuDistDataPu239_241(), G4fissionEvent::G4SmpNuDistDataU232_234_236_238(), G4fissionEvent::G4SmpNuDistDataU233_235(), G4fissionEvent::G4SmpNuDistDataU235(), G4fissionEvent::G4SmpNuDistDataU238(), G4ComponentAntiNuclNuclearXS::GetAntiHadronNucleonElCrSc(), G4ComponentAntiNuclNuclearXS::GetAntiHadronNucleonTotCrSc(), G4InuclParamMomDst::GetMomentum(), G4HadronNucleonXsc::HadronNucleonXscNS(), G4Bessel::I0(), G4Bessel::I1(), G4HETCFragment::IntegrateEmissionProbability(), G4NuclearAbrasionGeometry::P(), G4PreCompoundIon::ProbabilityDistributionFunction(), G4PreCompoundNucleon::ProbabilityDistributionFunction(), G4NuclearRadii::RadiusECS(), G4InuclSpecialFunctions::randomInuclPowers(), G4QGSMSplitableHadron::SampleX(), G4FissionProductYieldDist::SetNubar(), and G4BigBanger::xProbability().

powZ()

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

inline

Definition at line 225 of file G4Pow.hh.

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

References expA(), lz, and Z.

Referenced by G4WilsonAblationModel::BreakItUp(), G4KalbachCrossSection::ComputePowerParameter(), G4DensityEffectCalculator::DeltaOnceSolved(), G4DensityEffectCalculator::Ell(), G4NistManager::GetA27(), G4ICRU49NuclearStoppingModel::InitialiseArray(), G4PhotonEvaporation::InitialiseGRData(), G4NuclearRadii::Radius(), G4NuclearRadii::RadiusRMS(), G4HadronElastic::SampleInvariantT(), G4ChargeExchange::SampleT(), and G4IonCoulombCrossSection::SetScreenRSquare().

Z13()

G4double G4Pow::Z13 ( G4int  Z ) const

inline

Definition at line 123 of file G4Pow.hh.

{ return pz13[Z]; }

References pz13, and Z.

Referenced by G4GEMProbability::CalcAlphaParam(), G4NeutronEvaporationProbability::CalcAlphaParam(), G4StatMFMacroMultiplicity::CalcChemicalPotentialMu(), G4GEMCoulombBarrier::CalcCompoundRadius(), G4StatMFMicroCanonical::CalcEntropyOfCompoundNucleus(), G4StatMFMicroCanonical::CalcFreeInternalEnergy(), G4StatMFMicroPartition::CalcPartitionProbability(), G4GEMProbability::CalcProbability(), G4StatMFMacroCanonical::CalculateTemperature(), G4GEMProbabilityVI::ComputeTotalProbability(), G4Fancy3DNucleus::CoulombBarrier(), G4ionEffectiveCharge::EffectiveCharge(), G4CompetitiveFission::FissionKineticEnergy(), G4ScreeningMottCrossSection::FormFactor2UniformHelm(), G4StatMFMacroTemperature::FragsExcitEnergy(), G4InuclSpecialFunctions::G4cbrt(), G4GEMCoulombBarrier::G4GEMCoulombBarrier(), G4GEMProbabilityVI::G4GEMProbabilityVI(), G4IonisParamElm::G4IonisParamElm(), G4LevelManager::G4LevelManager(), G4NuclearFermiDensity::G4NuclearFermiDensity(), G4AntiProtonField::GetBarrier(), G4KaonMinusField::GetBarrier(), G4KaonPlusField::GetBarrier(), G4PionMinusField::GetBarrier(), G4PionPlusField::GetBarrier(), G4ProtonField::GetBarrier(), G4StatMFMicroPartition::GetCoulombEnergy(), G4IonsShenCrossSection::GetIsoCrossSection(), G4NuclearLevelData::GetLevelDensity(), G4StatMFMicroPartition::GetPartitionEnergy(), G4NistManager::GetZ13(), G4WentzelOKandVIxSection::InitialiseA(), G4VPreCompoundFragment::Initialize(), G4StatMFMacroCanonical::Initialize(), G4StatMFMicroCanonical::Initialize(), G4ICRU49NuclearStoppingModel::NuclearStoppingPower(), G4StatMFChannel::PlaceFragments(), G4NuclearRadii::Radius(), G4NuclearRadii::RadiusCB(), G4NuclearRadii::RadiusECS(), G4NuclearRadii::RadiusHNGG(), G4NuclearRadii::RadiusKNGG(), G4NuclearRadii::RadiusNNGG(), G4Generator2BS::SampleDirection(), G4AntiNuclElastic::SampleInvariantT(), G4HadronElastic::SampleInvariantT(), G4GEMChannel::SampleKineticEnergy(), G4IonCoulombCrossSection::SetScreenRSquare(), G4ScreeningMottCrossSection::SetupKinematic(), G4StatMFChannel::SolveEqOfMotion(), G4EvaporationProbability::TotalProbability(), and Z23().

Z23()

G4double G4Pow::Z23 ( G4int  Z ) const

inline

Definition at line 125 of file G4Pow.hh.

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

References Z, and Z13().

Referenced by G4FissionBarrier::BarashenkovFissionBarrier(), G4GEMProbability::CalcBetaParam(), G4NeutronEvaporationProbability::CalcBetaParam(), G4StatMFMacroMultiplicity::CalcChemicalPotentialMu(), G4StatMFMacroChemicalPotential::CalcChemicalPotentialNu(), G4StatMFMacroBiNucleon::CalcEnergy(), G4StatMFMacroMultiNucleon::CalcEnergy(), G4StatMFMacroTetraNucleon::CalcEnergy(), G4StatMFMacroTriNucleon::CalcEnergy(), G4StatMFMacroMultiNucleon::CalcEntropy(), G4StatMFFragment::CalcExcitationEnergy(), G4StatMFMacroBiNucleon::CalcMeanMultiplicity(), G4StatMFMacroMultiNucleon::CalcMeanMultiplicity(), G4StatMFMacroTetraNucleon::CalcMeanMultiplicity(), G4StatMFMacroTriNucleon::CalcMeanMultiplicity(), G4StatMFMicroPartition::CalcPartitionProbability(), G4StatMFMacroMultiNucleon::CalcZARatio(), G4StatMFMacroCanonical::ChooseZ(), G4UrbanAdjointMscModel::ComputeCrossSectionPerAtom(), G4UrbanMscModel::ComputeCrossSectionPerAtom(), G4AtimaFluctuations::Dispersion(), G4InuclSpecialFunctions::FermiEnergy(), G4NuclearShellModelDensity::G4NuclearShellModelDensity(), G4AntiNuclElastic::GetcosTeta1(), G4StatMFMicroPartition::GetCoulombEnergy(), G4LindhardSorensenData::GetDeltaL(), G4StatMFFragment::GetEnergy(), G4StatMFMicroPartition::GetPartitionEnergy(), G4StatMFMacroCanonical::Initialize(), G4StatMFMicroCanonical::Initialize(), G4HETCFragment::IntegrateEmissionProbability(), G4HadronElastic::SampleInvariantT(), and G4WentzelOKandVIxSection::SetupTarget().

Field Documentation

ener

G4DataVector G4Pow::ener

private

Definition at line 108 of file G4Pow.hh.

Referenced by G4Pow(), and logX().

fact

G4DataVector G4Pow::fact

private

Definition at line 115 of file G4Pow.hh.

Referenced by factorial(), and G4Pow().

fexp

G4DataVector G4Pow::fexp

private

Definition at line 114 of file G4Pow.hh.

Referenced by expA(), and G4Pow().

fpInstance

G4Pow * G4Pow::fpInstance = nullptr

staticprivate

Definition at line 98 of file G4Pow.hh.

Referenced by GetInstance().

logen

G4DataVector G4Pow::logen

private

Definition at line 109 of file G4Pow.hh.

Referenced by G4Pow(), and logX().

logfact

G4DataVector G4Pow::logfact

private

Definition at line 116 of file G4Pow.hh.

Referenced by G4Pow(), and logfactorial().

lowa13

G4DataVector G4Pow::lowa13

private

Definition at line 111 of file G4Pow.hh.

Referenced by A13Low(), and G4Pow().

lz

G4DataVector G4Pow::lz

private

Definition at line 112 of file G4Pow.hh.

Referenced by G4Pow(), log10A(), log10Z(), logBase(), logZ(), and powZ().

lz2

G4DataVector G4Pow::lz2

private

Definition at line 113 of file G4Pow.hh.

Referenced by G4Pow(), and logBase().

max2

const G4int G4Pow::max2 = 5

private

Definition at line 101 of file G4Pow.hh.

Referenced by G4Pow(), and logBase().

maxA

G4double G4Pow::maxA

private

Definition at line 103 of file G4Pow.hh.

Referenced by A13High(), G4Pow(), logBase(), and logX().

maxA2

G4double G4Pow::maxA2

private

Definition at line 105 of file G4Pow.hh.

Referenced by G4Pow(), and logBase().

maxAexp

G4double G4Pow::maxAexp

private

Definition at line 106 of file G4Pow.hh.

Referenced by expA(), and G4Pow().

maxLowA

G4double G4Pow::maxLowA

private

Definition at line 104 of file G4Pow.hh.

Referenced by A13(), and G4Pow().

onethird

const G4double G4Pow::onethird = 1.0 / 3.0

private

Definition at line 100 of file G4Pow.hh.

Referenced by A13High(), A13Low(), expA(), G4Pow(), and logBase().

pz13

G4DataVector G4Pow::pz13

private

Definition at line 110 of file G4Pow.hh.

Referenced by A13High(), G4Pow(), and Z13().

---

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

• source/global/management/include/G4Pow.hh
• source/global/management/src/G4Pow.cc