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