ShowerRRecData.h

Go to the documentation of this file.

#ifndef _evt_ShowerRRecData_h_
#define _evt_ShowerRRecData_h_

#include <fwk/RunController.h>

#include <utl/Point.h>
#include <utl/Vector.h>
#include <utl/LameShadowPtr.h>
#include <utl/ParameterStorage.h>

#include <det/Detector.h>

#include <evt/ShowerRRecDataQuantities.h>
#include <evt/Event.h>

#include <vector>

namespace evt {

  class ShowerRecData;


  class ShowerRRecData {

  public:
    // enums for reference core/position
    enum ReferenceCorePosition {
      eUseSDReferenceCore,
      eUseFDReferenceCore,
      eUseMCReferenceCore,
      eUseRDReferenceCore,
      eUseExternalReferenceCore,
      eUseRDCoordinateOrigin
    };

    enum ReferenceAxis {
      eUseSDReferenceAxis,
      eUseFDReferenceAxis,
      eUseMCReferenceAxis,
      eUseRDReferenceAxis,
      eUseExternalReferenceAxis
    };

    // Parameters Storage
    typedef revt::ShowerRRecDataQuantities Parameter;

    double GetParameter(const Parameter i) const
    { return fShowerQuantities.GetParameter(i, CurrentModule()); }

    bool GetParameterLockStatus(const Parameter i) const
    { return fShowerQuantities.GetParameterLockStatus(i, CurrentModule()); }

    double GetParameterError(const Parameter i) const
    { return std::sqrt(fShowerQuantities.GetParameterCovariance(i, i, CurrentModule())); }

    bool GetParameterErrorLockStatus(const Parameter i) const
    { return fShowerQuantities.GetParameterCovarianceLockStatus(i, i, CurrentModule()); }

    double GetParameterCovariance(const Parameter i1, const Parameter i2) const
    { return fShowerQuantities.GetParameterCovariance(i1, i2, CurrentModule()); }

    bool GetParameterCovarianceLockStatus(const Parameter i1, const Parameter i2) const
    { return fShowerQuantities.GetParameterCovarianceLockStatus(i1, i2, CurrentModule()); }

    bool HasParameter(const Parameter i) const
    { return fShowerQuantities.HasParameter(i); }

    bool HasParameterCovariance(const Parameter i1, const Parameter i2) const
    { return fShowerQuantities.HasParameterCovariance(i1, i2); }

    bool HasParameterError(const Parameter i) const
    { return fShowerQuantities.HasParameterCovariance(i, i); }

    void SetParameter(Parameter i, double value, bool lock = true)
    { fShowerQuantities.SetParameter(i, value, lock, CurrentModule()); }

    void SetParameterError(Parameter i, double value, bool lock = true)
    { fShowerQuantities.SetParameterCovariance(i, i, value * value, lock, CurrentModule()); }

    void SetParameterCovariance(Parameter i1, Parameter i2, double value, bool lock = true)
    { fShowerQuantities.SetParameterCovariance(i1, i2, value, lock, CurrentModule()); }

    void DeleteParameter(const Parameter i)
    { fShowerQuantities.DeleteParameter(i, CurrentModule()); }

    void DeleteParameterCovariance(const Parameter i1, const Parameter i2)
    { fShowerQuantities.DeleteParameterCovariance(i1, i2, CurrentModule()); }

    std::vector<Parameter> GetEnumVector() const
    { return fShowerQuantities.GetEnumVector(); }

    std::vector<std::pair<Parameter, Parameter>> GetCovarianceEnumVector() const
    { return fShowerQuantities.GetCovarianceEnumVector(); }

    double GetRadius() const;

    double GetRadiusError() const;

    double GetAzimuth() const;

    double GetAzimuthError() const;

    double GetZenith() const;

    double GetZenithError() const;

    double GetAzimuthPreFit() const;

    double GetAzimuthPreFitError() const;

    double GetZenithPreFit() const;

    double GetZenithPreFitError() const;

    utl::Point GetCoordinateOrigin() const;

    bool HasCorePosition() const;

    utl::Point GetCorePosition() const;

    utl::Vector GetCoreError() const;

    double GetCoreCorrelationXY() const;

    bool HasAxis() const;

    utl::Vector GetAxis() const;

    utl::Vector GetMagneticFieldVector() const;

    ReferenceCorePosition GetReferenceCoreFlag() const { return fReferenceCoreFlag; }

    ReferenceAxis GetReferenceAxisFlag() const { return fReferenceAxisFlag; }

    void SetReferenceCoreFlag(const ReferenceCorePosition refCoreFlag) { fReferenceCoreFlag = refCoreFlag; }

    void SetReferenceAxisFlag(const ReferenceAxis refAxisFlag) { fReferenceAxisFlag = refAxisFlag; }

    bool HasReferenceCorePosition(const Event& event) const;

    bool HasReferenceAxis(const Event& event) const;

    utl::Point GetReferenceCorePosition(const Event& event) const;

    utl::Vector GetReferenceCoreError(const Event& event) const;

    double GetReferenceCoreErrorCorrelationXY(const Event& event) const;

    utl::Vector GetReferenceAxis(const Event& event) const;

    unsigned int GetWhichEye() const { return fWhichEye; }

    void SetWhichEye(const unsigned int eye) { fWhichEye = eye; }

  private:
    ShowerRRecData() = default;

    static const std::string& CurrentModule() { return fwk::RunController::GetInstance().GetCurrentModule(); }

    utl::ParameterStorage<Parameter> fShowerQuantities{"Shower:ShowerRRecData"};

    ReferenceCorePosition fReferenceCoreFlag;
    ReferenceAxis fReferenceAxisFlag;

    unsigned int fWhichEye = 0;

    friend class evt::ShowerRecData;
    friend class utl::LameShadowPtr<ShowerRRecData>;

  };

}


#endif