MuCCaFitter.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/

// MuCCaFitter.h
//   Header file for class MuCCaFitter
// (c) ATLAS Detector software
//
 
#ifndef MuCCaFITTER_H
#define MuCCaFITTER_H
 
#include <vector>
 
#include "MdtCalibInterfaces/IMdtSegmentFitter.h"
#include "MuonCalibEventBase/MuonCalibSegment.h"
 
namespace MuonCalib {

    class MuCCaFitter : public IMdtSegmentFitter {
    public:
        bool fit(MuonCalibSegment& seg) const;

        bool fit(MuonCalibSegment& seg, HitSelection selection) const;

        void printLevel(int level);
 
    private:
        bool m_debug;

        //  double getY( const Amg::Vector3D& p ) const { return p.y(); }
        //  double getX( const Amg::Vector3D& p ) const { return p.x(); }
        //  double getZ( const Amg::Vector3D& p ) const { return p.z(); }
        //  Amg::Vector3D getVec( double x, double y, double z) const { return Amg::Vector3D( x, y, z ); }
        double getY(const Amg::Vector3D& p) const { return p.z(); }
        double getX(const Amg::Vector3D& p) const { return p.x(); }
        double getZ(const Amg::Vector3D& p) const { return p.y(); }
        Amg::Vector3D getVec(double x, double y, double z) const { return Amg::Vector3D(x, z, y); }
    };

    class MuCCaFitterImplementation {
    public:
        void Computelin(double x1, double y1, double r1, double x2, double y2, double r2);
        void Computelinparnew(double x1, double y1, double r1, double x2, double y2, double r2);
        void computehitsfromcircles(double x0, double y0, double r0, double x1, double y1, double r1, double a, double b);
        //  void Computehitpsemes(int nhit,double *xcirc, double *ycirc,
        //                        double *rcirc, double a, double b);
        void Computehitpsemes(int nhit, const std::vector<double>& xcirc, const std::vector<double>& ycirc, const std::vector<double>& rcirc, double a,
                              double b);
        //  void Computeparam3(int number_of_hits,double x[100],double y[100],
        //                     double r[100],double sr[100]);
        void Computeparam3(int number_of_hits, const std::vector<double>& x, const std::vector<double>& y, const std::vector<double>& r, const std::vector<double>& sr);
        double get_a();
        double get_b();
        double get_da();
        double get_db();
        double get_corrab();
        double get_chi2f();
 
    private:
        double m_angularcoefficient[4];
        double m_bfpar[4];
        double m_xout0;
        double m_yout0;
        double m_xout1;
        double m_yout1;
        double m_xpoint[100];
        double m_ypoint[100];
        double m_aoutn;   
        double m_sig2a;   
        double m_bout;    
        double m_sig2b;   
        double m_corrab;  
        double m_chi2f;   
    };
 
}  // namespace MuonCalib
#endif