event_selector.cc

Go to the documentation of this file.

#include "km3net-dataformat/offline/Evt.hh"
#include "km3net-dataformat/definitions/fitparameters.hh"
#include "JReconstruction/JEvt.hh"
#include "JReconstruction/JEvtToolkit.hh"
#include "JAAnet/JAAnetToolkit.hh"
#include "JMath/JMathToolkit.hh"

using JRECONSTRUCTION::JEvt;

/**
 * Event selection.
 *
 * \param  in            input event
 * \param  event         pointer to Monte Carlo event
 * \return               true if accepted; else false
 */
extern "C" bool accept(const JEvt& in, const Evt* event)
{
  using namespace JPP;

  if (!in.empty() && in[0].getQ() >= 100) {

    if (event != NULL) {

      Trk muon;

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

      return (is_muon(muon) && getAngle(getDirection(in[0]), getDirection(muon)) > 5.0);
    }
  }

  return true;
  //return false;
}