SdBaselineFinderKG.h

Go to the documentation of this file.

#ifndef _SdBaselineFinderKG_SdBaselineFinderKG_h_
#define _SdBaselineFinderKG_SdBaselineFinderKG_h_

#include <fwk/VModule.h>
#include <sevt/Station.h>
#include <utl/Is.h>


namespace SdBaselineFinderKG {

  class SdBaselineFinderKG : public fwk::VModule {

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

  private:
    void MakeFlatBaseline(sevt::PMT& pmt, const sdet::PMTConstants::PMTGain gain) const;
    bool ComputeBaselines(sevt::Station& station) const;
    bool ComputeBaseline(const sevt::Station& station, sevt::PMT& pmt, const sdet::PMTConstants::PMTGain gain) const;

    std::vector<double> ComputeBaselineEstimates(const bool isUUB,
                                                 const utl::TraceI& trace, const int start, const int baselineLength) const;

    // from XML config
    int fNumberOfIterations = 0;

    typedef std::pair<int, int> PairI;
    PairI fBaselineEstimationLengths{0, 0};
    PairI fUbBaselineEstimationLengths{0, 0};
    PairI fUubBaselineEstimationLengths{0, 0};

    typedef std::pair<double, double> PairD;
    PairD fNegativeSignalFactors{0, 0};
    PairD fUbNegativeSignalFactors{0, 0};
    PairD fUubNegativeSignalFactors{0, 0};

    PairD fUubDecayConstant{0, 0};
    std::set<int> fRecalculateLSTriggerForTestStations;

    PairD fUpperSigma{0, 0};
    PairD fUbUpperSigma{0, 0};
    PairD fUubUpperSigma{0, 0};

    double fLowerSigma = 0;
    double fUbLowerSigma = 0;
    double fUubLowerSigma = 0;

    double fTruncationSigma = 0;
    double fUbTruncationSigma = 0;
    double fUubTruncationSigma = 0;

    int fMinimumPeakSize = 0;
    int fUbMinimumPeakSize = 0;
    int fUubMinimumPeakSize = 0;

    REGISTER_MODULE("SdBaselineFinderKG", SdBaselineFinderKG);

  };

}


#endif