Auxiliary class to determine EM-equivalent energy as a function of PDG particle code and energy. More...

#include <JPythia.hh>

Public Member Functions
	JPythia ()
	Default constructor. More...

double	operator() (const int type, const double Es) const
	Get EM-equivalent energy. More...

Static Public Member Functions
static double	getEnergy (const double Es)
	Get equivalent EM-energy for given pion energy. More...

Detailed Description

Auxiliary class to determine EM-equivalent energy as a function of PDG particle code and energy.

Definition at line 21 of file JPythia.hh.

Constructor & Destructor Documentation

◆ JPythia()

JSIRENE::JPythia::JPythia ( )

inline

Default constructor.

Definition at line 25 of file JPythia.hh.

26 {}

Member Function Documentation

◆ operator()()

double JSIRENE::JPythia::operator()	(	const int	type,
		const double	Es
	)		const

inline

Get EM-equivalent energy.

Parameters

type	particle type [PDG]
Es	particle energy [GeV]

Returns: EM-equivalent energy [GeV]

Definition at line 36 of file JPythia.hh.

     {
       using namespace JAANET;
  
       switch (type) {
         
       case TRACK_TYPE_PHOTON:
       case TRACK_TYPE_NEUTRAL_PION:
       case TRACK_TYPE_NEUTRAL_ANTIPION:
       case TRACK_TYPE_ELECTRON:
       case TRACK_TYPE_ANTIELECTRON:
         return Es;
         
       case TRACK_TYPE_CHARGED_PION_PLUS:
       case TRACK_TYPE_CHARGED_PION_MINUS:
         return getEnergy(Es);
  
       default:
         return 0.0;
       }
     }

◆ getEnergy()

static double JSIRENE::JPythia::getEnergy ( const double Es )

inlinestatic

Get equivalent EM-energy for given pion energy.

Reference: Mona Dentler, "Investigation of the One-Particle Approximation in the ANTARES simulation package KM3",
Bachelorarbeit Erlangen Centre for Astroparticle Physics,
Friedrich-Alexander-Universität, Erlangen-Nürnberg.

Parameters

Es	particle energy [GeV]

Returns: EM-equivalent energy [GeV]

Definition at line 70 of file JPythia.hh.

     {
       static const double a     =  72.425;
       static const double b     = -49.417;
       static const double c     =   5.858;
       static const double d     = 207.252;
       static const double e     = 132.784;
       static const double f     = -10.277;
       static const double g     = -19.441;
       static const double h     =  58.598;
       static const double i     =  53.161;
       static const double kref  =   2.698;
       
       static const double u     = (a - f) / kref;
  
       const double x = log10(Es);
       const double y = (e + x*(d + x*(c + x*(b + x*(a + x*u)))))  /  (i + x*(h + x*(g + x*(f + x*u))));
  
       return pow(10.0, y - kref);
     }

The documentation for this struct was generated from the following file:

JPythia.hh