JSIRENE::JPythia Struct Reference

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

JSIRENE::JPythia::JPythia ( )

inline

Default constructor.

Definition at line 25 of file JPythia.hh.

    {}

Member Function Documentation

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;
      }
    }

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