JAAnet/event_selector.cc

Go to the documentation of this file.

#include "km3net-dataformat/offline/Evt.hh"
#include "km3net-dataformat/definitions/fitparameters.hh"
#include "km3net-dataformat/tools/reconstruction.hh"
#include "JAAnet/JAAnetToolkit.hh"
#include "JMath/JMathToolkit.hh"

/**
 * Event selection.
 *
 * \param  trk           track
 * \param  evt           event
 * \return               true if accepted; else false
 */
extern "C" bool accept(const Trk& trk, const Evt& evt)
{
  using namespace std;
  using namespace JPP;
  /*
  if (has_muon(evt)) {

    Trk muon;

    for (const auto& t1 : evt.mc_trks) {
      if (is_muon(t1)) {
        if (t1.E > muon.E) {
          muon = t1;
        }
      }
    }

    return (is_muon(muon) && getAngle(getDirection(trk), getDirection(muon)) > 5.0);
  }
  */
  return true;
}