QuasianalyticLineReconstruction.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef QuasianalyticLineReconstructionH
#define QuasianalyticLineReconstructionH
 
//:::::::::::::::::::::::::::::::::::::::::::
//:: CLASS QuasianalyticLineReconstruction ::
//:::::::::::::::::::::::::::::::::::::::::::
 
 
#include <vector>
 
#include "CxxUtils/checker_macros.h"
#include "MdtCalibFitters/DCSLFitter.h"
#include "MdtCalibFitters/IndexSet.h"
#include "MdtCalibFitters/MTStraightLine.h"
#include "MdtCalibInterfaces/IMdtPatRecFitter.h"
#include "MuonCalibEventBase/MuonCalibSegment.h"
 
namespace MuonCalib {
 
    class QuasianalyticLineReconstruction : public IMdtPatRecFitter {
    public:
        using MdtHitPtr = MuonCalibSegment::MdtHitPtr;
        using MdtHitVec = MuonCalibSegment::MdtHitVec;
        // Constructors //
        QuasianalyticLineReconstruction() { init(); }
 
        QuasianalyticLineReconstruction(const double  r_road_width) { init(r_road_width); }
 
        // Methods //
        // get-methods //
        double roadWidth() const;
 
        // set-method //
        void setRoadWidth(const double  r_road_width);
        void setTimeOut(const double  time_out);
        void setMaxRadius(const double  maxR);
        // methods required by the base class "IMdtSegmentFitter" //
        bool fit(MuonCalibSegment& r_segment) const;
        bool fit(MuonCalibSegment& r_segment, HitSelection r_selection) const;
        bool fit(MuonCalibSegment& r_segment, HitSelection r_selection, MTStraightLine& final_track) const;
        void printLevel(int /*level*/){};
 
    private:
        // internal co-ordinate definition //
        //                      x3
        //                      ^
        //     o o o o o o      |
        //  ... o o o o o ...   o--> x2
        //     o o o o o o     x1
        //
 
        // parameters for the adjustment of the track reconstruction //
        double m_r_max;       // maximum radius
        double m_road_width;  // road width for pattern recognition
        double m_time_out;    // time-out for track finding
 
        // chi^2 refitter //
        DCSLFitter m_nfitter;  // NIKHEF straight line reconstruction
 
        // internal storage vectors //
        std::vector<MTStraightLine> m_tangent;
        // vector of tangents (for 1 track)
        std::vector<int> m_nb_cand_hits;          // number of hits on track candidate
        std::vector<MTStraightLine> m_candidate;  // track candidates
 
        // initialization methods //
        void init();
        // default initialization:  road width = 0.5 CLHEP::mm
        void init(const double  r_road_width);
        // initialization with user-defined road width
 
        // auxiliary methods //
        MTStraightLine tangent(const Amg::Vector3D& r_w1, const double  r_r1, const double  r_sigma12, const Amg::Vector3D& r_w2,
                               const double  r_r2, const double  r_sigma22, const int& r_case) const;
        // method for the calculation of tangents with errors;
        // r_w1: wire position for the first hit,
        // r_r1: drift radius of the first hit,
        // r_sigma12: sigma(r_r1)^2,
        // r_w2: wire position for the second hit,
        // r_r2: drift radius of the second hit,
        // r_sigma22: sigma(r_r2)^2,
        // r_case = 1, 2, 3, 4: select one of the four cases of a tangent
 
        MTStraightLine track_candidate(const IndexSet& r_index_set, const int& r_k_cand, const int& r_l_cand, const int& r_cand_case,
                                       const std::vector<Amg::Vector3D>& r_w, const std::vector<double>& r_r, const std::vector<double>& r_sigma2,
                                       double& r_chi2) const;
        // method for the calculatation of a straight line from tangents;
        // r_nb_selected_hits: number of selected hits,
        // r_k_cand: index of the 1st wire of the candidate tangent
        //           (starting from 0),
        // r_l_cand: index of the 2nd wire of the candidate tangent
        //           (starting from 0),
        // r_cand_case: configuration of the tangent to be used (1, 2, 3, or 4)
        // r_w: wire positions,
        // r_r: drift radii,
        // r_sigma2: sigma(r)^2,
        // r_chi2: chi^2 for the reconstructed track
    };
 
}  // namespace MuonCalib
 
#endif