MdEventSelector.h

Go to the documentation of this file.

#ifndef _MdEventSelectorAG_MdEventSelector_h_
#define _MdEventSelectorAG_MdEventSelector_h_

#include <fwk/VModule.h>
//
#include <mevt/MEvent.h>
#include <mevt/Counter.h>
#include <mevt/Module.h>
#include <sevt/SEvent.h>

#include <mdet/MDetector.h>
#include <mdet/Counter.h>


namespace evt {
  class Event;
}

namespace MdEventSelectorAG {

  class MdEventSelector : public fwk::VModule {

  public:
    MdEventSelector();
    virtual ~MdEventSelector();

    fwk::VModule::ResultFlag Init();
    fwk::VModule::ResultFlag Run(evt::Event& event);
    fwk::VModule::ResultFlag Finish();

  private:
    const mdet::MDetector* theMDetector = nullptr;

    unsigned int fMinNumberOfCounters = 0;
    double fSilentRange = 0;
    bool fSelectEventByTheta = false;
    double fMaxTheta = 0;

    bool fRejectNonT5Events = false;
    int fT5NumberOfActiveStations = 0;

    bool fRejectSaturatedEvents = false;
    bool fRejectDenseStations = false;

    // parameters to filter time outliers
    bool fRejectTimeOutlierCounters = false;
    int fTimeOutlierMuons = 0;
    double fTimeOutlierDistance = 0;
    double fTimeOutlierConstant = 0;
    double fTimeOutlierSlope = 0;

    // selection parameters
    double fFirstCrownMinDistance = 0;
    double fFirstCrownMaxDistance = 0;

    bool HasEnoughDetectors(const mevt::MEvent& me) const;

    void RejectFarSilents(mevt::MEvent& me) const;

    const mevt::Counter* FindHottestCounter(const mevt::MEvent& me) const;
    const mevt::Counter* FindHottestCounterUnsaturated(const mevt::MEvent& me) const;
    const mevt::Counter* FindHottestCounterSaturated(const mevt::MEvent& me) const;

    bool IsHighTheta(const evt::Event& e) const;

    void RejectCloseSimDetectors(mevt::MEvent& me) const;

    void RejectDenseDetectors(mevt::MEvent& me) const;

    void RejectTimeOutliers(mevt::MEvent& me , const sevt::SEvent& se) const;

    bool IsMEventT5(const mevt::MEvent& me) const;

    void SetT5(evt::Event& e, const bool t5) const;

    // Check if the counter t50 is out of time with the SD start time
    bool IsTimeOutlier(const double md2sdTimeInterval, const double distance) const;

    REGISTER_MODULE("MdEventSelectorAG", MdEventSelector);

  };

}


#endif