SdTraceCalibrator.h

Go to the documentation of this file.

#ifndef _SdTraceCalibratorOG_SdTraceCalibrator_h_
#define _SdTraceCalibratorOG_SdTraceCalibrator_h_

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


namespace SdTraceCalibratorOG {

  class SdTraceCalibrator : 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:
    int BuildSignals(sevt::Station& station) const;
    bool BuildSignals(sevt::PMT& pmt, const unsigned int traceLength, const int saturationValue) const;
    bool MakeComponentVEMTraces(sevt::PMT& pmt) const;
    void SumPMTComponents(sevt::Station& station) const;
    bool MergeSignals(sevt::Station& station) const;
    bool SelectSignal(sevt::Station& station) const;
    void ComputeShapeRiseFallPeak(sevt::PMTRecData& pmtRecData,
                                  const double binTiming,
                                  const unsigned int startBin,
                                  const unsigned int startIntegration,
                                  const unsigned int endIntegration,
                                  const double traceIntegral) const;

    bool IsTestStation(const int id) const { return utl::Is(id).In(fTestStations); }

    // from XML config file
    double fPMTSummationCutoff = 0;
    bool fTreatHGLGEqualInSignalSearch = false;
    double fFindSignalThreshold = 0;
    int fFindSignalBinsAboveThreshold = 0;
    bool fIncludeWaterCherenkovDetectorInScintillatorStartStopDetermination = false;
    int fBinsBeforeStartForSPMT = 0;

    typedef std::pair<double, double> PairD;
    PairD fRiseTimeFractions{0, 0};
    PairD fFallTimeFractions{0, 0};

    std::set<int> fTestStations;

    bool fIsUUB = false;

    REGISTER_MODULE("SdTraceCalibratorOG", SdTraceCalibrator);

  };

}


#endif