SmallPMTCalibData.h

Go to the documentation of this file.

#ifndef _sevt_SmallPMTCalibData_h_
#define _sevt_SmallPMTCalibData_h_

#include <utl/ShadowPtr_fwd.h>


namespace sevt {

  class SmallPMTData;


  class SmallPMTCalibData {

  public:
    bool IsTubeOk() const { return fIsTubeOk; }
    // Get version of the small-shower inter-calibration algorithm
    int GetVersion() const { return fVersion; }

    // Calibration factor - Q_LPMT(VEM) = beta * Q_SPMT(FADC counts)
    // where Q_LPMT is the average over the unmasked LPMT signals
    double GetBeta() const { return fBeta; }
    double GetBetaError() const { return fBetaError; }
    // Calibration factors and errors obtained wrt each single LPMT
    double GetBeta(const unsigned int lpmtId) const { return Checked(fBetaLPMTs, lpmtId); }
    double GetBetaError(const unsigned int lpmtId) const { return Checked(fBetaErrorLPMTs, lpmtId); }

    // Chi^2 obtained from the SPMT calibration procedure
    // using the average over the unmasked LPMT signals
    double GetChi2() const { return fChi2; }
    // Same wrt each LPMT
    double GetChi2(const unsigned int lpmtId) const { return Checked(fChi2LPMTs, lpmtId); }

    // Correction factor (usually = 1)
    // For instance, correction wrt temperature
    // Also, negative if the calibration is unsuccessful or problematic
    double GetCorrectionFactor() const { return fCorrectionFactor; }
    // Same wrt each LPMT
    double GetCorrectionFactor(const unsigned int lpmtId) const { return Checked(fCorrectionFactorLPMTs, lpmtId); }

    void SetIsTubeOk(const bool ok) { fIsTubeOk = ok; }
    void SetVersion(const int vers) { fVersion = vers; }

    void SetBeta(const double beta) { fBeta = beta; }
    void SetBetaError(const double betaerr) { fBetaError = betaerr; }
    void SetChi2(const double chi2) { fChi2 = chi2; }
    void SetCorrectionFactor(const double fact) { fCorrectionFactor = fact; }
    // Set values for each LPMT
    void SetBeta(const double beta, const unsigned int lpmtId) { Checked(fBetaLPMTs, lpmtId) = beta; }
    void SetBetaError(const double betaerr, const unsigned int lpmtId) { Checked(fBetaErrorLPMTs, lpmtId) = betaerr; }
    void SetChi2(const double chi2, const unsigned int lpmtId) { Checked(fChi2LPMTs, lpmtId) = chi2; }
    void SetCorrectionFactor(const double fact, const unsigned int lpmtId) { Checked(fCorrectionFactorLPMTs, lpmtId) = fact; }

  private:
    SmallPMTCalibData() = default;
    ~SmallPMTCalibData() = default;

    template<typename T, std::size_t n>
    static T& Checked(T (&)[n], const unsigned int i);

    bool fIsTubeOk = false;

    int fVersion = 0;
    double fBeta = 0;
    double fBetaError = 0;
    double fChi2 = 0;
    double fCorrectionFactor = -1;

    double fBetaLPMTs[3] = { 0 };
    double fBetaErrorLPMTs[3] = { 0 };
    double fChi2LPMTs[3] = { 0 };
    double fCorrectionFactorLPMTs[3] = { -1, -1, -1 };

    friend class SmallPMTData;
    friend class utl::ShadowPtr<SmallPMTCalibData>;

  };

}


#endif