CalibratedSpacePoint.h

Go to the documentation of this file.

/*
   Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
#ifndef MUONSPACEPOINT_CALIBSPACEPOINT_H
#define MUONSPACEPOINT_CALIBSPACEPOINT_H
 
#include <MuonSpacePoint/SpacePoint.h>
 
namespace MuonR4{
    class CalibratedSpacePoint {
        public:
            using Cov_t = SpacePoint::Cov_t;
            using CovIdx = SpacePoint::CovIdx;
            enum class State : std::uint8_t {
                Valid = 1,
                FailedCalib,
                Outlier,
                Duplicate
            };

            static std::string toString(const State s);
            CalibratedSpacePoint(const SpacePoint* uncalibSpacePoint,
                                 Amg::Vector3D&& posInChamber,
                                 State st = State::Valid);
 
            ~CalibratedSpacePoint() = default;
            /*** @brief: Position of the space point inside the chamber */
            const Amg::Vector3D& localPosition() const;
            /*** @brief: Returns the direction parallel to the primary channel, i.e. the strip or the wire */
            const Amg::Vector3D& sensorDirection() const;
            /*** @brief: Returns the vector pointing to the adjacent channel in the chamber */
            const Amg::Vector3D& toNextSensor() const;
            const Amg::Vector3D& planeNormal() const;
            double time() const;
            bool isStraw() const;
            bool hasTime() const;
            bool measuresLoc0() const;
            bool measuresLoc1() const;
            double driftRadius() const;
            const Cov_t& covariance() const; 

            void setCovariance(const Cov_t& cov);
            void setDriftRadius(const double r);
            
            const SpacePoint* spacePoint() const;
            xAOD::UncalibMeasType type() const;
            void setTimeMeasurement(double t);
            bool measuresPhi() const;
            bool measuresEta() const;
            State fitState() const;
            void setFitState(State st);
            unsigned dimension() const;
            void setBeamDirection(Amg::Vector3D&& beamDir);
 
            friend std::ostream& operator<<(std::ostream& ostr, const CalibratedSpacePoint& sp) {
                    sp.print(ostr);
                    return ostr;
            }
        private:
            void print(std::ostream& ostr) const;
            Amg::Vector3D m_posInChamber{Amg::Vector3D::Zero()};
            Cov_t m_cov{Acts::filledArray<double, 3>(0.)};
            double m_driftRadius{0.};
            double m_time{0.};
            const SpacePoint* m_parent{nullptr};
            std::shared_ptr<Amg::Vector3D> m_beamLine{};
            State m_state{State::Valid};
            bool m_measuresTime{type() == xAOD::UncalibMeasType::MdtDriftCircleType};
        
        };
        static_assert(Acts::Experimental::CompositeSpacePoint<CalibratedSpacePoint>);
 
 
}
 
#endif