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