utl::InvisibleEnergy Namespace Reference

Enumerations
enum	ECompositionModel { eMixedComposition, eProton, eIron, eData }
enum	EInteractionModel { eQGSJET01, eSIBYLL21, eDATA }

Functions
double	EnergyLimits (const double energy)
double	Factor (const double Eem, const EInteractionModel iMod, const ECompositionModel iCompo, const double CosTheta)
	invisible energy factor, finv=Etot/Eem, given Eem. CosTheta only needed when using data driven estimation More...
double	FactorDerivative (const double Eem, const EInteractionModel iMod, const ECompositionModel iCompo, const double CosTheta)
	derivative of invisible energy factor dfinv/dEem given Eem. CosTheta only needed when using data driven estimation More...
double	FactorVariance (const double eCal, const double eTot)
const double *	FitParameters (const ECompositionModel composition)
double	ModelFactor (const EInteractionModel interaction, const ECompositionModel composition)

Variables
constexpr double	kData [] = { 0.179, 0.947, 1.95, 4, 0.846, -0.265, 0.489, -0.441 }
constexpr double	kQ2Sib [] = { 1.016, 1.012 }
constexpr double	kQGSJetIron [] = { 0.977, 0.109, 0.130 }
constexpr double	kQGSJetMixed [] = { 0.967, 0.078, 0.140 }
constexpr double	kQGSJetProton [] = { 0.958, 0.048, 0.162 }

Enumeration Type Documentation

enum utl::InvisibleEnergy::ECompositionModel

Enumerator
eMixedComposition
eProton
eIron
eData

Definition at line 206 of file PhysicalFunctions.h.

enum utl::InvisibleEnergy::EInteractionModel

Enumerator
eQGSJET01
eSIBYLL21
eDATA

Definition at line 199 of file PhysicalFunctions.h.

Function Documentation

double utl::InvisibleEnergy::EnergyLimits

(

const double

energy

)

Definition at line 537 of file PhysicalFunctions.cc.

References galactic::EeV, and max.

Referenced by Factor(), and FactorDerivative().

double utl::InvisibleEnergy::Factor	(	const double	energyEM,
		const EInteractionModel	intMod,
		const ECompositionModel	compMod,
		const double	cosTheta
	)

invisible energy factor, finv=Etot/Eem, given Eem. CosTheta only needed when using data driven estimation

returns the the invisible energy factor from H.M.J. Barbosa et al., Astropart. Phys. 22 (2004) 159-166 values are for 45 degree

input: electromagnetic energy

Definition at line 588 of file PhysicalFunctions.cc.

References RdGeoCeLDFFitter::a, RdGeoCeLDFFitter::b, degree, eDATA, galactic::EeV, EnergyLimits(), FitParameters(), ModelFactor(), G4StationSimulatorOG::p, and std::pow().

Referenced by FdProfileReconstructorKG::FdProfileReconstructor::CalculateTotalErrors(), FdEnergyFinderOG::FdEnergyFinder::FindTotalEnergy(), UnivFitterKG::WCDFitFunction::operator()(), and FdEnergyDepositFinderKG::FdEnergyDepositFinder::PropagateUncertainties().

double utl::InvisibleEnergy::FactorDerivative	(	const double	energyEM,
		const EInteractionModel	intMod,
		const ECompositionModel	compMod,
		const double	cosTheta
	)

derivative of invisible energy factor dfinv/dEem given Eem. CosTheta only needed when using data driven estimation

   \returns the derivative of the invisible energy factor from
   H.M.J. Barbosa et al., Astropart. Phys. 22 (2004) 159-166
   values are for 45 degree

Add the Invisible Energy model from GAPs GAP2012-124 & GAP2013-026

   input: electromagnetic energy
   output: df/dE

return (p[3] - 1) * p[1] * pow(p[5], p[2]) * pow(energyLim/EeV/p[5], p[3]-1) /p[5] * p[4] / energyLim; PRD Paper

Definition at line 642 of file PhysicalFunctions.cc.

References RdGeoCeLDFFitter::a, RdGeoCeLDFFitter::b, degree, eDATA, galactic::EeV, EnergyLimits(), FitParameters(), ModelFactor(), G4StationSimulatorOG::p, std::pow(), and utl::Sqr().

Referenced by FdProfileReconstructorKG::FdProfileReconstructor::CalculateTotalErrors(), and FdEnergyDepositFinderKG::FdEnergyDepositFinder::PropagateUncertainties().

double utl::InvisibleEnergy::FactorVariance	(	const double	eCal,
		const double	eTot
	)

Returns

statistical variance of invisible energy due to shower to shower fluctuations

derived from fit to table from T. Pierog for QGSJET/proton (note: iron has smaller fluctuations–>upper limit)

Returns

statistical variance of invisible energy due to shower to shower fluctuations and to the mass variance

Definition at line 710 of file PhysicalFunctions.cc.

References max, std::pow(), utl::Sqr(), and sqrt().

Referenced by FdProfileReconstructorKG::FdProfileReconstructor::CalculateTotalErrors(), and FdEnergyDepositFinderKG::FdEnergyDepositFinder::PropagateUncertainties().

const double * utl::InvisibleEnergy::FitParameters

(

const ECompositionModel

composition

)

Definition at line 567 of file PhysicalFunctions.cc.

References eData, RdREASSimPreparator::eIron, eMixedComposition, RdREASSimPreparator::eProton, ERROR, kData, kQGSJetIron, kQGSJetMixed, and kQGSJetProton.

Referenced by Factor(), FactorDerivative(), and UnivFitterKG::WCDFitFunction::operator()().

double utl::InvisibleEnergy::ModelFactor	(	const EInteractionModel	interaction,
		const ECompositionModel	composition
	)

Definition at line 546 of file PhysicalFunctions.cc.

References eDATA, eData, RdREASSimPreparator::eIron, eMixedComposition, RdREASSimPreparator::eProton, eQGSJET01, ERROR, eSIBYLL21, and kQ2Sib.

Referenced by Factor(), FactorDerivative(), and UnivFitterKG::WCDFitFunction::operator()().

Variable Documentation

constexpr double utl::InvisibleEnergy::kData[] = { 0.179, 0.947, 1.95, 4, 0.846, -0.265, 0.489, -0.441 }

Definition at line 533 of file PhysicalFunctions.cc.

Referenced by FitParameters().

constexpr double utl::InvisibleEnergy::kQ2Sib[] = { 1.016, 1.012 }

Definition at line 521 of file PhysicalFunctions.cc.

Referenced by ModelFactor().

constexpr double utl::InvisibleEnergy::kQGSJetIron[] = { 0.977, 0.109, 0.130 }

Definition at line 518 of file PhysicalFunctions.cc.

Referenced by FitParameters().

constexpr double utl::InvisibleEnergy::kQGSJetMixed[] = { 0.967, 0.078, 0.140 }

Definition at line 519 of file PhysicalFunctions.cc.

Referenced by FitParameters().

constexpr double utl::InvisibleEnergy::kQGSJetProton[] = { 0.958, 0.048, 0.162 }

Definition at line 517 of file PhysicalFunctions.cc.

Referenced by FitParameters().