MuonKinkFinder.h

Go to the documentation of this file.

#ifndef RecoMuon_MuonIdentification_MuonKinkFinder_h
#define RecoMuon_MuonIdentification_MuonKinkFinder_h

#include "DataFormats/MuonReco/interface/MuonQuality.h"
#include "TrackingTools/TrackRefitter/interface/TrackTransformer.h"

class MuonKinkFinder {
    public:
        MuonKinkFinder(const edm::ParameterSet &iConfig);
        ~MuonKinkFinder() ;

        // set event setup
        void init(const edm::EventSetup &iSetup) ;

        // fill data, return false if refit failed or too few hits
        bool fillTrkKink(reco::MuonQuality & quality, const Trajectory &trajectory) const ;

        // fill data, return false if refit failed or too few hits
        bool fillTrkKink(reco::MuonQuality & quality, const reco::Track &track) const ;

    private:
        bool diagonalOnly_;
        bool usePosition_;

        TrackTransformer refitter_;

        // compute chi2 between track states
        double getChi2(const TrajectoryStateOnSurface &start, const TrajectoryStateOnSurface &other) const ;

        // possibly crop matrix or set to zero off-diagonal elements, then invert
        void cropAndInvert(AlgebraicSymMatrix55 &cov) const ;
};
#endif