RdChannelBandstopFilter.h

Go to the documentation of this file.

#ifndef _RdChannelBandstopFilter_RdChannelBandstopFilter_h_
#define _RdChannelBandstopFilter_RdChannelBandstopFilter_h_

#include <fwk/VModule.h>

#include <revt/Channel.h>
#include <rdet/Channel.h>

#include <utl/AugerUnits.h>


namespace RdChannelBandstopFilter {

  class RdChannelBandstopFilter : public fwk::VModule {

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

  private:
    //xml settings
    std::vector<double> fLowerBoundsNS;
    std::vector<double> fUpperBoundsNS;
    std::vector<double> fLowerBoundsEW;
    std::vector<double> fUpperBoundsEW;

    bool fUseOnlineBandstopFilter = false;
    std::string fOnlineBandstopFilterSpectrumBaseline;
    double fOnlineBandstopFilterDeviationFactor = 0;

    // needed to assign channels with vectors of frequencies to cut
    const double fTolerance = 2 * utl::deg;
    const double fNorthSouth = 1.51006;
    const double fEastWest = -0.0607375;

    bool CloseTo(const double a, const double b, const double tolerance);
    void CutNoiseFrequencies(revt::Channel& channel, const rdet::Channel& detChannel);

    REGISTER_MODULE("RdChannelBandstopFilter", RdChannelBandstopFilter);
  };
}

#endif