RdEventPreSelector.h

Go to the documentation of this file.

#ifndef _RdEventPreSelector_RdEventPreSelector_h_
#define _RdEventPreSelector_RdEventPreSelector_h_

#include <fwk/VModule.h>

#include <utl/AugerUnits.h>
#include <utl/TimeStamp.h>

#include <set>
#include <vector>


namespace revt {
  class EventTrigger;
}


namespace RdEventPreSelector {

  class RdEventPreSelector : public fwk::VModule {

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

  private:
    int fMinEventId = 0;

    int fMaxEventId = 0;

    bool fUseTriggerInformation = false;
    bool fUseSelfTriggeredEvent = false;
    bool fUseExternallyTriggeredEvent = false;
    bool fUseCalibrationTriggeredEvent = false;
    bool fUseScintillatorTriggeredEvent = false;
    bool fUseSDTriggeredEvent = false;
    bool fUseGUITriggeredEvent = false;
    bool fUseFDTriggeredEvent = false;
    bool fUseHEATTriggeredEvent = false;
    bool fUseAERAletTriggeredEvent = false;
    bool fUseAIRPLANETriggeredEvent = false;
    bool fUsePeriodicTriggeredEvent = false;
    bool fUsePassThroughTriggeredEvent = false;

    std::set<int> fSelectedEventIds;

    std::set<int> fDeselectedEventIds;

    std::set<utl::TimeStamp> fSelectedEventDates;

    std::set<std::pair<int, int> > fSelectedRunEventIds;

    std::set<std::pair<int, int> > fDeselectedRunEventIds;

    int fMinNumberOfStations = 0;

    int fMaxNumberOfStations = 0;

    std::vector<int> fAllParticipatingStationIds;

    std::vector<int> fAtLeastOneParticipatingStationIds;

    std::vector<int> fNotParticipatingStationIds;

    unsigned int fNRejectedMinNumberOfStations = 0;

    unsigned int fNRejectedStationRejectionList = 0;

    unsigned int fNRejectedEventRejectionList = 0;

    unsigned int fNRejectedTriggerConditions = 0;

    unsigned int fNumberOfEvents = 0;

    unsigned int fNRejected = 0;

    unsigned int fNAccepted = 0;

    utl::TimeStamp fAllowedTimeBegin;
    utl::TimeStamp fAllowedTimeStop;

    bool fBreakLoopAtMaxDateTime = true;

    // min ADC signal to noise ratio
    double fMinSignalToNoise = 0;

    // min number of station with pulse
    double fMinNumberOfStationsWithPulse = 0;

    // max number of stations available in event with most stations still rejected by MinNumberOfStations cut
    double fMaxNumberOfAvailableStationsInLargestEvents = 0;

    unsigned int fNRejectedMinNumberOfStationsWithPulse = 0;

    // minimum allowed SD angle
    double fMinSDZenith = 0 * utl::deg;

    // maximum allowed SD angle
    double fMaxSDZenith = 90 * utl::deg;

    // minimum allowed SD energy
    double fMinSDEnergy = 0;

    std::vector<unsigned long> fEventGPSSeconds;
    double fAllowedTimeAroundEventTime = 10 * utl::minute;

    bool fFDInput = false;

    unsigned int fWhichEye = 6;

    double fMaxFDZenith = 90 * utl::deg;

    bool CheckTrigger(const revt::EventTrigger& evtrig);

    // counts the number of stations with pulse
    unsigned int CountStationWithPulse(evt::Event& event) const;

    bool CloseTo(const double a, const double b, const double tolerance) const;

    REGISTER_MODULE("RdEventPreSelector", RdEventPreSelector);
  };

}


#endif