km3-opa_efrac.hh

Go to the documentation of this file.

#ifndef __KM3_OPA_EFRAC__
#define __KM3_OPA_EFRAC__
 
#include "JAAnet/JPDB.hh"
 
 
/**
 * \author mdejong
 */
 
extern "C"
{
  float opa_efrac_(int& ipart, float& ekin);
  float opa_weight_high_e_(float& ekin);
}
 
namespace JSIRENE {}
namespace JPP { using namespace JSIRENE; }
 
namespace JSIRENE {
 
  using JAANET::JPDB;
  
  /**
   * Get equivalent EM-energy for given pion energy.
   *
   * \param   type    particle type   [PDG]
   * \param   E       particle energy [GeV]
   * \return          EM-equivalent energy [GeV]
   */
  inline double pythia(const int type, const double E)
  {
    int   ipart = JPDB::getInstance().getPDG(type).geant;
    float ekin  = (float) E;
 
    return opa_efrac_(ipart, ekin) * E;
  }
}
 
#endif