TelescopeRecData.h

Go to the documentation of this file.

#ifndef _fevt_TelescopeRecData_h_
#define _fevt_TelescopeRecData_h_

#include <fevt/Eye.h>

#include <utl/TimeStamp.h>
#include <utl/AxialVector.h>
#include <utl/MultiTabulatedFunctionErrors.h>
#include <utl/ShadowPtr.h>
#include <utl/LameShadowPtr.h>

#include <boost/iterator/indirect_iterator.hpp>

#include <vector>
#include <list>


namespace fevt {

  class Pixel;


  class TelescopeRecData {

  public:
    double GetZeta() const { return fZeta; }

    const utl::AxialVector& GetSDP() const { return fSDP; }
    double GetSDPThetaError() const { return fSDPThetaError; }
    double GetSDPPhiError() const { return fSDPPhiError; }
    double GetSDPCorrThetaPhi() const { return fSDPThetaPhiCorrelation; }
    double GetSDPFitChiSquare() const { return fSDPFitChiSquare; }
    unsigned int GetSDPFitNDof() const { return fSDPFitNDof; }

    double GetTZero() const { return fTZero; }
    double GetTZeroError() const { return fTZeroError; }
    double GetChiZero() const { return fChiZero; }
    double GetChiZeroError() const { return fChiZeroError; }
    double GetRp() const { return fRp; }
    double GetChi0TZeroCorrelation() const { return fChi0TZeroCorrelation; }
    double GetRpChi0Correlation() const { return fRpChi0Correlation; }
    double GetRpTZeroCorrelation() const { return fRpTZeroCorrelation; }
    double GetRpError() const { return fRpError; }
    double GetTimeFitChiSquare() const { return fTimeFitChiSquare; }
    unsigned int GetTimeFitNDof() const { return fTimeFitNDof; }

    double GetNorthEastCorrelation() const { return fNorthEastCorrelation; }
    double GetNorthThetaCorrelation() const { return fNorthThetaCorrelation; }
    double GetNorthPhiCorrelation() const { return fNorthPhiCorrelation; }
    double GetNorthTCoreCorrelation() const { return fNorthTCoreCorrelation; }
    double GetEastThetaCorrelation() const { return fEastThetaCorrelation; }
    double GetEastPhiCorrelation() const { return fEastPhiCorrelation; }
    double GetEastTCoreCorrelation() const { return fEastTCoreCorrelation; }
    double GetThetaPhiCorrelation() const { return fThetaPhiCorrelation; }
    double GetThetaTCoreCorrelation() const { return fThetaTCoreCorrelation; }
    double GetPhiTCoreCorrelation() const { return fPhiTCoreCorrelation; }
    double GetAxisFitChiSquare() const { return fAxisFitChiSquare; }
    unsigned int GetAxisFitNDof() const { return fAxisFitNDof; }

    unsigned int GetEyeId() const { return fEyeId; }
    unsigned int GetTelescopeId() const { return fTelescopeId; }

    void SetZeta(const double zeta) { fZeta = zeta; }

    void SetSDP(const utl::AxialVector& vec) { fSDP = vec; }

    void SetSDPThetaError(const double sdpThetaError)
    { fSDPThetaError = sdpThetaError; }
    void SetSDPPhiError(const double sdpPhiError)
    { fSDPPhiError = sdpPhiError; }
    void SetSDPCorrThetaPhi (double sdpCorrThetaPhi)
    { fSDPThetaPhiCorrelation=sdpCorrThetaPhi; }
    void SetSDPFitChiSquare(const double sdpChi2, const unsigned int ndof)
    { fSDPFitChiSquare = sdpChi2; fSDPFitNDof = ndof; }

    void SetTimeFitChiSquare(const double tfitChi2, const unsigned int ndof)
    { fTimeFitChiSquare = tfitChi2; fTimeFitNDof = ndof; }

    void SetTZero(const double tzero, const double error)
    { fTZero = tzero; fTZeroError = error; }
    void SetChiZero(const double chiZero, const double error)
    { fChiZero = chiZero; fChiZeroError = error; }
    void SetRp(const double rp, const double error)
    {  fRp = rp; fRpError = error; }
    void SetTimeFitCorrelations(double rRpT0, double rRpChi0, double rChi0T0);

    void SetAxisFitCorrelations(double northEast,
                                double northTheta,
                                double northPhi,
                                double northTCore,
                                double eastTheta,
                                double eastPhi,
                                double eastTCore,
                                double thetaPhi,
                                double thetaTCore,
                                double phiTCore)
    {
      fNorthEastCorrelation = northEast;
      fNorthThetaCorrelation = northTheta;
      fNorthPhiCorrelation = northPhi;
      fNorthTCoreCorrelation = northTCore;
      fEastThetaCorrelation = eastTheta;
      fEastPhiCorrelation = eastPhi;
      fEastTCoreCorrelation = eastTCore;
      fThetaPhiCorrelation = thetaPhi;
      fThetaTCoreCorrelation = thetaTCore;
      fPhiTCoreCorrelation = phiTCore;
    }
    void SetAxisFitChiSquare(const double globalChi2, const unsigned int ndof)
    { fAxisFitChiSquare = globalChi2; fAxisFitNDof = ndof; }


    utl::TabulatedFunctionErrors& GetLightFlux(const FdConstants::LightSource source = fevt::FdConstants::eTotal)
    { return fLightProfile->GetTabulatedFunctionErrors(source); }
    const utl::TabulatedFunctionErrors& GetLightFlux(const FdConstants::LightSource source = fevt::FdConstants::eTotal) const
    { return fLightProfile->GetTabulatedFunctionErrors(source); }

    void MakeLightFlux(const FdConstants::LightSource source = FdConstants::eTotal)
    { fLightProfile->AddTabulatedFunctionErrors(source); }

    bool HasLightFlux(const FdConstants::LightSource source = FdConstants::eTotal) const
    { return fLightProfile->HasLabel(source); }

    typedef utl::MultiTabulatedFunctionErrors::Iterator MultiLightFluxIterator;
    typedef utl::MultiTabulatedFunctionErrors::ConstIterator ConstMultiLightFluxIterator;

    MultiLightFluxIterator LightFluxesBegin()
    { return fLightProfile->TabulatedFunctionsErrorsBegin(); }
    ConstMultiLightFluxIterator LightFluxesBegin() const
    { return fLightProfile->TabulatedFunctionsErrorsBegin(); }

    MultiLightFluxIterator LightFluxesEnd()
    { return fLightProfile->TabulatedFunctionsErrorsEnd(); }
    ConstMultiLightFluxIterator LightFluxesEnd() const
    { return fLightProfile->TabulatedFunctionsErrorsEnd(); }

    utl::TimeStamp GetPhotonsStartTime() const { return fTraceStartTime; }

    void SetPhotonsStartTime(const utl::TimeStamp& ts)
    { fTraceStartTime = ts; }


    const std::vector<std::vector<unsigned int>>& GetPixelsInZetaOverTime() const
    { return fPixelsInZeta; }
    std::vector<std::vector<unsigned int>>& GetPixelsInZetaOverTime()
    { return fPixelsInZeta; }

    utl::MultiTabulatedFunctionErrors& GetLightCollectionEfficiency()
    { return *fLightCollectionEfficiency; }
    const utl::MultiTabulatedFunctionErrors& GetLightCollectionEfficiency() const
    { return *fLightCollectionEfficiency; }
    void MakeLightCollectionEfficiency();
    bool HasLightCollectionEfficiency() const { return bool(fLightCollectionEfficiency); }

    bool IsSpotFarFromBorder(const double timeRelToEyeTriggerTime) const;
    void SetSpotFarFromBorderTimeRanges(const std::list<std::pair<double, double>>& timeRanges)
    { fSpotFarFromBorderTimes = timeRanges; }
    const std::list<std::pair<double, double>>& GetSpotFarFromBorderTimeRanges() const
    { return fSpotFarFromBorderTimes; }
    std::list<std::pair<double, double>>& GetSpotFarFromBorderTimeRanges()
    { return fSpotFarFromBorderTimes; }

  private:
    TelescopeRecData(const unsigned int eyeId, const unsigned int telId)
      : fEyeId(eyeId), fTelescopeId(telId) { }

    ~TelescopeRecData() { }

    void CorrectIds(const unsigned int eyeId, const unsigned int telescopeId)
    { fEyeId = eyeId; fTelescopeId = telescopeId; }

    unsigned int fEyeId = 0;
    unsigned int fTelescopeId = 0;

    double fZeta = 0;

    utl::AxialVector fSDP;
    double fSDPThetaError = 0;
    double fSDPPhiError = 0;
    double fSDPThetaPhiCorrelation = 0;
    double fSDPFitChiSquare = 0;
    unsigned int fSDPFitNDof = 0;

    double fTZero = 0;
    double fTZeroError = 0;
    double fChiZero = 0;
    double fChiZeroError = 0;
    double fRp = 0;
    double fRpError = 0;
    double fRpTZeroCorrelation = 0;
    double fRpChi0Correlation = 0;
    double fChi0TZeroCorrelation = 0;
    double fTimeFitChiSquare = 0;
    unsigned int fTimeFitNDof = 0;

    double fNorthEastCorrelation = 0;
    double fNorthThetaCorrelation = 0;
    double fNorthPhiCorrelation = 0;
    double fNorthTCoreCorrelation = 0;
    double fEastThetaCorrelation = 0;
    double fEastPhiCorrelation = 0;
    double fEastTCoreCorrelation = 0;
    double fThetaPhiCorrelation = 0;
    double fThetaTCoreCorrelation = 0;
    double fPhiTCoreCorrelation = 0;
    double fAxisFitChiSquare = 0;
    unsigned int fAxisFitNDof = 0;

    utl::InitializedShadowPtr<utl::MultiTabulatedFunctionErrors> fLightProfile;

    std::list<std::pair<double, double>> fSpotFarFromBorderTimes;
    utl::TimeStamp fTraceStartTime;
    std::vector<std::vector<unsigned int>> fPixelsInZeta;
    utl::ShadowPtr<utl::MultiTabulatedFunctionErrors> fLightCollectionEfficiency;

    friend class Telescope;
    friend class utl::ShadowPtr<TelescopeRecData>;
  };

}


#endif