CalibratedSpacePoint.h

Go to the documentation of this file.

/*
   Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
#ifndef MUONSPACEPOINT_CALIBSPACEPOINT_H
#define MUONSPACEPOINT_CALIBSPACEPOINT_H
 
#include <MuonSpacePoint/SpacePoint.h>
#include <variant>
 
namespace MuonR4{
    class CalibratedSpacePoint {
        public:
            
            enum class State : uint8_t {
                Valid = 0,
                FailedCalib = 1,
                Outlier = 2,
            };
            
            CalibratedSpacePoint(const SpacePoint* uncalibSpacePoint,
                                 Amg::Vector3D&& posInChamber,
                                 Amg::Vector3D&& dirInChamber,
                                 State st = State::Valid);
 
            ~CalibratedSpacePoint() = default;
            const Amg::Vector3D& positionInChamber() const;
            const Amg::Vector3D& directionInChamber() const;
            double driftRadius() const;
            void setDriftRadius(const double r);
            using Covariance_t = std::variant<AmgSymMatrix(2), AmgSymMatrix(3)>;
            const Covariance_t& covariance() const;
            template <unsigned k>
                void setCovariance(const AmgSymMatrix(k)& cov){
                    static_assert(k==2 || k==3, "Covariance dimension needs to be 2 or 3");
                    m_cov = cov;
                }
            const SpacePoint* spacePoint() const;
            xAOD::UncalibMeasType type() const;
            double time() const;
            void setTimeMeasurement(double t);
            bool measuresTime() const;
            bool measuresPhi() const;
            bool measuresEta() const;
            State fitState() const;
            void setFitState(State st);
            unsigned dimension() const;
        private:
            const SpacePoint* m_parent{nullptr};
            Amg::Vector3D m_posInChamber{Amg::Vector3D::Zero()};
            Amg::Vector3D m_dirInChamber{Amg::Vector3D::Zero()};
            
            double m_driftRadius{0.};
            Covariance_t m_cov{};
 
            double m_time{0.};
            bool m_measuresTime{type() == xAOD::UncalibMeasType::MdtDriftCircleType};
            State m_state{State::Valid};
    };
 
}
 
#endif