SdCalibrator.h

Go to the documentation of this file.

#ifndef _SdCalibratorOG_SdCalibrator_h_
#define _SdCalibratorOG_SdCalibrator_h_

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


namespace sevt {
  class PMTRecData;
}

namespace SdCalibratorOG {

  class SdCalibrator : 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 ApplyTimeCorrection(sevt::Station& station);
    void CalculatePeakAndCharge(sevt::Station& station);

    //bool ApplyStartTimeCorrection(sevt::Station& station);
    void SumPMTComponents(sevt::Station& station) const;

    // new offline calibration
    bool CopySmallPMTCalibData(sevt::Station& station);
    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;
    bool MakeComponentVEMTraces(sevt::PMT& pmt) const;
    int  BuildSignals(sevt::Station& station) const;
    bool BuildSignals(sevt::PMT& pmt, const unsigned int traceLength, const int saturationValue) const;
    bool MergeSignals(sevt::Station& station) const;
    bool SelectSignal(sevt::Station& station) const;
    void TreatSaturated(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); }

    // XML config
    double fPMTSummationCutoff = 0;

    typedef std::pair<double, double> PairD;
    PairD fPeakFitRange{0, 0};
    double fPeakFitChi2Accept = 0;
    double fPeakOnlineToVEMFactor = 0;
    double fPeakHistogramToVEMFactor = 0;
    double fPeakOnlineToMIPFactor = 0;
    double fPeakHistogramToMIPFactor = 0;
    double fChargeFitShoulderHeadRatio = 0;
    double fChargeFitChi2Accept = 0;
    double fChargeOnlineToVEMFactor = 0;
    double fChargeHistogramToVEMFactor = 0;
    double fChargeOnlineToMIPFactor = 0;
    double fChargeHistogramToMIPFactor = 0;
    PairD fShapeFitRangeBefore12{0, 0};
    PairD fShapeFitRangeSince12{0, 0};
    PairD fRiseTimeFractions{0, 0};
    PairD fFallTimeFractions{0, 0};
    std::set<int> fTestStations;
    double fFindSignalThreshold = 0;
    int fFindSignalBinsAboveThreshold = 0;
    bool fTreatHGLGEqualInSignalSearch = false;
    bool fApplyBackwardFlatPieceCheck = false;
    bool fDecreaseLGFlatPieceTolerance = false;
    bool fIncludeWaterCherenkovDetectorInScintillatorStartStopDetermination = false;
    int fBinsBeforeStartForSPMT = 0;

    bool fToldYaPeak = false;
    bool fToldYaCharge = false;
    bool fToldYaShape = false;

    bool fIsUUB = false;

    const utl::TimeRange kCommsCrisis{utl::TimeStamp{928378815}, utl::TimeStamp{942375615}};

    REGISTER_MODULE("SdCalibratorOG", SdCalibrator);

  };

}


#endif