CurvedLine.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
#ifndef MuonCalib_CurvedLineH
#define MuonCalib_CurvedLineH
 
//::::::::::::::::::::::
//:: CLASS CurvedLine ::
//::::::::::::::::::::::
 
 
//::::::::::::::::::
//:: HEADER FILES ::
//::::::::::::::::::
 
#include <vector>
 
#include "EventPrimitives/EventPrimitives.h"
#include "GeoPrimitives/GeoPrimitives.h"
#include "MdtCalibFitters/MTStraightLine.h"
#include "MuonCalibEventBase/MuonCalibSegment.h"
#include "MuonCalibMath/Legendre_polynomial.h"
 
namespace MuonCalib {
 
    class CurvedLine {
    public:
        using MdtHitVec = MuonCalibSegment::MdtHitVec;
        using MdtHitPtr = MuonCalibSegment::MdtHitPtr;
 
        // Constructors //
        CurvedLine();
 
        CurvedLine(std::vector<Amg::Vector3D> &points);
 
        CurvedLine(std::vector<Amg::Vector3D> &points, std::vector<Amg::Vector3D> &x_and_y_errors);
 
        // Methods //
        // get-methods //
        Amg::Vector3D getPointOnLine(const double &loc_z) const;
 
        MTStraightLine getTangent(const double &loc_z) const;
 
        void setChi2(double chi2);
        double chi2() const;
 
        void setNumberOfTrackHits(unsigned int n_hits);
        unsigned int numberOfTrackHits() const;
        double chi2PerDegreesOfFreedom() const;
 
        void setUsedHits(const MdtHitVec &hits);
        const MdtHitVec &trackHits() const;
 
    private:
        // parameters of the curved line //
        const Legendre_polynomial *m_Legendre = nullptr;  // pointer to the Legendre polynomial
                                                // describing the curved line
        Amg::VectorX m_coeff_xz;                // coefficients of the straight line in the local
                                                // xz plane
        Amg::VectorX m_coeff_yz;                // coefficients of the curved line in the local
                                                // yz plane
 
        double m_chi2{-1};
        unsigned int m_numTrkHits{0};
        MdtHitVec m_used_hits;
 
        // private methods //
        void init(std::vector<Amg::Vector3D> &points, std::vector<Amg::Vector3D> &x_and_y_errors);
        // initialization routine
    };
 
}  // namespace MuonCalib
 
#endif