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