Functions
G4double	toDegrees (G4double radians)
G4int	heaviside (G4int n)
G4double	pow13 (G4double x)
G4double	powMinus13 (G4double x)
G4double	pow23 (G4double x)
G4double	aSinH (G4double x)
template<typename T>
G4int	sign (const T t)
template<typename T>
T	max (const T t1, const T t2)
	brief Return the largest of the two arguments
template<typename T>
T	min (const T t1, const T t2)
	brief Return the smallest of the two arguments
G4double	gaussianCDF (const G4double x)
	Cumulative distribution function for Gaussian.
G4double	gaussianCDF (const G4double x, const G4double x0, const G4double sigma)
	Generic cumulative distribution function for Gaussian.
G4double	inverseGaussianCDF (const G4double x)
	Inverse cumulative distribution function for Gaussian.
G4double	arcSin (const G4double x)
	Calculates arcsin with some tolerance on illegal arguments.
G4double	arcCos (const G4double x)
	Calculates arccos with some tolerance on illegal arguments.

Variables
const G4double	pi = 3.14159265358979323846264338328
const G4double	twoPi = 2.0 * pi
const G4double	tenPi = 10.0 * pi
const G4double	piOverTwo = 0.5 * pi
const G4double	oneOverSqrtTwo = 1./std::sqrt((G4double)2.)
const G4double	oneOverSqrtThree = 1./std::sqrt((G4double)3.)
const G4double	oneThird = 1./3.
const G4double	twoThirds = 2./3.
const G4double	sqrtFiveThirds = std::sqrt(5./3.)
const G4double	sqrtThreeFifths = std::sqrt(3./5.)

Function Documentation

◆ arcCos()

G4double G4INCL::Math::arcCos ( const G4double x )

Calculates arccos with some tolerance on illegal arguments.

Definition at line 103 of file G4INCLGlobals.cc.

                                      {
// assert(x>-1.000001 && x<1.000001);
      return ((x > 1.) ? 0. : ((x<-1.) ? pi : std::acos(x)));
    }

Referenced by G4INCL::CoulombNonRelativistic::distortOut(), G4INCL::EventInfo::fillInverseKinematics(), G4INCL::SurfaceAvatar::getTransmissionProbability(), and G4INCL::EventInfo::remnantToParticle().

◆ arcSin()

G4double G4INCL::Math::arcSin ( const G4double x )

Calculates arcsin with some tolerance on illegal arguments.

Definition at line 98 of file G4INCLGlobals.cc.

                                      {
// assert(x>-1.000001 && x<1.000001);
      return ((x > 1.) ? 0. : ((x<-1.) ? pi : std::asin(x)));
    }

◆ aSinH()

G4double G4INCL::Math::aSinH ( G4double x )

inline

Definition at line 101 of file G4INCLGlobals.hh.

                                      {
      return std::log(x + std::sqrt(x*x+1.));
    }

◆ gaussianCDF() [1/2]

G4double G4INCL::Math::gaussianCDF ( const G4double x )

Cumulative distribution function for Gaussian.

A public-domain approximation taken from Abramowitz and Stegun. Applies to a Gaussian with mean=0 and sigma=1.

Parameters

x	a Gaussian variable

Definition at line 74 of file G4INCLGlobals.cc.

    {
      // Save the sign of x
      const G4double sgn = sign(x);
      const G4double z = std::fabs(x) * oneOverSqrtTwo;
 
      // A&S formula 7.1.26
      G4double t = 1.0/(1.0 + gcdfp*z);
      G4double y = 1.0 - (((((gcdfa5*t + gcdfa4)*t) + gcdfa3)*t + gcdfa2)*t + gcdfa1)*t*std::exp(-z*z);
 
      return 0.5*(1.0 + sgn*y);
    }

Referenced by G4INCL::Random::correlatedUniform(), and gaussianCDF().

◆ gaussianCDF() [2/2]

G4double G4INCL::Math::gaussianCDF	(	const G4double	x,
		const G4double	x0,
		const G4double	sigma )

Generic cumulative distribution function for Gaussian.

A public-domain approximation taken from Abramowitz and Stegun. Applies to a generic Gaussian.

Parameters

x	a Gaussian variable
x0	mean of the Gaussian
sigma	standard deviation of the Gaussian

Definition at line 87 of file G4INCLGlobals.cc.

                                                                                    {
      return gaussianCDF((x-x0)/sigma);
    }

◆ heaviside()

G4int G4INCL::Math::heaviside ( G4int n )

inline

Definition at line 84 of file G4INCLGlobals.hh.

                                    {
      if(n < 0) return 0;
      else return 1;
    }

Referenced by G4INCL::Nucleus::insertParticle().

◆ inverseGaussianCDF()

G4double G4INCL::Math::inverseGaussianCDF ( const G4double x )

Inverse cumulative distribution function for Gaussian.

A public-domain approximation taken from Abramowitz and Stegun. Applies to a Gaussian with mean=0 and sigma=1.

Parameters

x	a uniform variate

Returns: a Gaussian variate

Definition at line 91 of file G4INCLGlobals.cc.

                                            {
      if (x < 0.5)
        return -inverseGaussianCDFRational( std::sqrt(-2.0*std::log(x)) );
      else
        return inverseGaussianCDFRational( std::sqrt(-2.0*std::log(1.-x)) );
    }

◆ max()

template<typename T>

T G4INCL::Math::max	(	const T	t1,
		const T	t2 )

inline

brief Return the largest of the two arguments

Definition at line 113 of file G4INCLGlobals.hh.

                                                               {
      return t1 > t2 ? t1 : t2;
    }

Referenced by G4INCL::ParticleTable::getLargestNuclearRadius().

◆ min()

template<typename T>

T G4INCL::Math::min	(	const T	t1,
		const T	t2 )

inline

brief Return the smallest of the two arguments

Definition at line 118 of file G4INCLGlobals.hh.

                                                               {
      return t1 < t2 ? t1 : t2;
    }

◆ pow13()

G4double G4INCL::Math::pow13 ( G4double x )

inline

Definition at line 89 of file G4INCLGlobals.hh.

                                      {
      return std::pow(x, oneThird);
    }

Referenced by G4INCL::NuclearDensityFactory::createRPCorrelationTable(), G4INCL::NuclearPotential::INuclearPotential::INuclearPotential(), and G4INCL::Random::sphereVector().

◆ pow23()

G4double G4INCL::Math::pow23 ( G4double x )

inline

Definition at line 97 of file G4INCLGlobals.hh.

                                      {
      return std::pow(x, twoThirds);
    }

◆ powMinus13()

G4double G4INCL::Math::powMinus13 ( G4double x )

inline

Definition at line 93 of file G4INCLGlobals.hh.

                                           {
      return std::pow(x, -oneThird);
    }

◆ sign()

template<typename T>

G4int G4INCL::Math::sign ( const T t )

inline

A simple sign function that allows us to port fortran code to c++ more easily.

Definition at line 108 of file G4INCLGlobals.hh.

                                                       {
      return t > 0 ? 1: t < 0 ? -1 : 0;
    }

Referenced by G4INCL::DeltaProductionChannel::fillFinalState(), G4INCL::ElasticChannel::fillFinalState(), G4INCL::NNToMissingStrangenessChannel::fillFinalState(), G4INCL::NpiToMissingStrangenessChannel::fillFinalState(), gaussianCDF(), and G4INCL::RootFinder::solve().

◆ toDegrees()

G4double G4INCL::Math::toDegrees ( G4double radians )

inline

Definition at line 80 of file G4INCLGlobals.hh.

                                                {
      return radians * (180.0 / pi);
    }

Referenced by G4INCL::Nucleus::fillEventInfo(), G4INCL::EventInfo::fillInverseKinematics(), and G4INCL::EventInfo::remnantToParticle().

Variable Documentation

◆ oneOverSqrtThree

const G4double G4INCL::Math::oneOverSqrtThree = 1./std::sqrt((G4double)3.)

Definition at line 74 of file G4INCLGlobals.hh.

Referenced by G4INCL::NuclearDensityFactory::createPCDFTable(), G4INCL::NuclearDensityFactory::createRCDFTable(), G4INCL::NuclearDensityFactory::createRPCorrelationTable(), G4INCL::NbarAtrestEntryChannel::densityN(), G4INCL::PbarAtrestEntryChannel::densityN(), G4INCL::NbarAtrestEntryChannel::densityP(), G4INCL::PbarAtrestEntryChannel::densityP(), and G4INCL::Random::gaussVector().

◆ oneOverSqrtTwo

const G4double G4INCL::Math::oneOverSqrtTwo = 1./std::sqrt((G4double)2.)

Definition at line 73 of file G4INCLGlobals.hh.

Referenced by gaussianCDF().

◆ oneThird

const G4double G4INCL::Math::oneThird = 1./3.

Definition at line 75 of file G4INCLGlobals.hh.

Referenced by pow13(), and powMinus13().

◆ pi

const G4double G4INCL::Math::pi = 3.14159265358979323846264338328

Definition at line 69 of file G4INCLGlobals.hh.

Referenced by arcCos(), arcSin(), G4INCL::NNToNSKpiChannel::fillFinalState(), G4INCL::NpiToLK2piChannel::fillFinalState(), G4INCL::NpiToLKpiChannel::fillFinalState(), G4INCL::NpiToNKKbChannel::fillFinalState(), G4INCL::NpiToSK2piChannel::fillFinalState(), G4INCL::NpiToSKpiChannel::fillFinalState(), G4INCL::PauliStandard::getBlockingProbability(), G4INCL::CrossSectionsAntiparticles::NNbarToAnnihilation(), G4INCL::INCL::prepareReaction(), and toDegrees().

◆ piOverTwo

const G4double G4INCL::Math::piOverTwo = 0.5 * pi

Definition at line 72 of file G4INCLGlobals.hh.

Referenced by G4INCL::CoulombNonRelativistic::distortOut().

◆ sqrtFiveThirds

const G4double G4INCL::Math::sqrtFiveThirds = std::sqrt(5./3.)

Definition at line 77 of file G4INCLGlobals.hh.

Referenced by G4INCL::ParticleTable::getFermiMomentumConstantLight().

◆ sqrtThreeFifths

const G4double G4INCL::Math::sqrtThreeFifths = std::sqrt(3./5.)

Definition at line 78 of file G4INCLGlobals.hh.

Referenced by G4INCL::ParticleTable::getMomentumRMS().

◆ tenPi

const G4double G4INCL::Math::tenPi = 10.0 * pi

Definition at line 71 of file G4INCLGlobals.hh.

Referenced by G4INCL::StandardPropagationModel::generateBinaryCollisionAvatar(), G4INCL::BinaryCollisionAvatar::getChannel(), G4INCL::CrossSections::interactionDistanceKbarN(), G4INCL::CrossSections::interactionDistanceKN(), G4INCL::CrossSections::interactionDistanceNbarN(), G4INCL::CrossSections::interactionDistancenbarN(), G4INCL::CrossSections::interactionDistanceNN(), G4INCL::CrossSections::interactionDistancePiN(), G4INCL::CrossSections::interactionDistanceYN(), and G4INCL::INCL::prepareReaction().

◆ twoPi

const G4double G4INCL::Math::twoPi = 2.0 * pi

Definition at line 70 of file G4INCLGlobals.hh.

Referenced by G4INCL::DeuteronDensity::derivDensityR(), G4INCL::DeltaProductionChannel::fillFinalState(), G4INCL::ElasticChannel::fillFinalState(), G4INCL::Random::gaussWithMemory(), G4INCL::NbarAtrestEntryChannel::getAnnihilationPosition(), G4INCL::PbarAtrestEntryChannel::getAnnihilationPosition(), G4INCL::PauliStandard::getBlockingProbability(), G4INCL::NKbElasticChannel::KaonMomentum(), G4INCL::NKbToLpiChannel::KaonMomentum(), G4INCL::NKbToNKbChannel::KaonMomentum(), G4INCL::NKbToSpiChannel::KaonMomentum(), G4INCL::NpiToLKChannel::KaonMomentum(), G4INCL::NpiToSKChannel::KaonMomentum(), and G4INCL::Random::normVector().

◆ twoThirds

const G4double G4INCL::Math::twoThirds = 2./3.

Definition at line 76 of file G4INCLGlobals.hh.

Referenced by G4INCL::InteractionAvatar::postInteraction(), and pow23().

Functions

Variables

Function Documentation

◆ arcCos()

◆ arcSin()

◆ aSinH()

◆ gaussianCDF() [1/2]

◆ gaussianCDF() [2/2]

◆ heaviside()

◆ inverseGaussianCDF()

◆ max()

◆ min()

◆ pow13()

◆ pow23()

◆ powMinus13()

◆ sign()

◆ toDegrees()

Variable Documentation

◆ oneOverSqrtThree

◆ oneOverSqrtTwo

◆ oneThird

◆ pi

◆ piOverTwo

◆ sqrtFiveThirds

◆ sqrtThreeFifths

◆ tenPi

◆ twoPi

◆ twoThirds