MuonResidualsFromTrack.h

Go to the documentation of this file.

#ifndef Alignment_MuonAlignmentAlgorithms_MuonResidualsFromTrack_H
#define Alignment_MuonAlignmentAlgorithms_MuonResidualsFromTrack_H

#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Utilities/interface/Exception.h"
#include "Alignment/CommonAlignment/interface/AlignableNavigator.h"
#include "Geometry/CommonDetUnit/interface/GlobalTrackingGeometry.h"
#include "TrackingTools/PatternTools/interface/Trajectory.h"
#include "TrackingTools/TrackFitters/interface/TrajectoryStateCombiner.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/MuonDetId/interface/DTChamberId.h"
#include "DataFormats/MuonDetId/interface/DTSuperLayerId.h"
#include "DataFormats/MuonDetId/interface/CSCDetId.h"
#include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/TrackReco/interface/Track.h"

#include "TMatrixDSym.h"
#include "TMatrixD.h"

#include <vector>
#include <map>

#include "Alignment/MuonAlignmentAlgorithms/interface/MuonChamberResidual.h"

class MuonResidualsFromTrack
{
public:
  // residuals from global muon trajectories
  MuonResidualsFromTrack(edm::ESHandle<GlobalTrackingGeometry> globalGeometry,
                         const Trajectory *traj,
                         const reco::Track* trk,
                         AlignableNavigator *navigator, double maxResidual);

  // residuals from tracker muons
  MuonResidualsFromTrack(edm::ESHandle<GlobalTrackingGeometry> globalGeometry,
                         const reco::Muon *mu,
                         AlignableNavigator *navigator,
                         double maxResidual);

  ~MuonResidualsFromTrack();
  
  void clear();

  const reco::Track *getTrack() { return track; }

  int trackerNumHits() const { return m_tracker_numHits; }

  double trackerChi2() const { return m_tracker_chi2; }
  double trackerRedChi2() const;
  double normalizedChi2() const;

  bool contains_TIDTEC() const { return m_contains_TIDTEC; }

  const std::vector<DetId> chamberIds() const { return m_chamberIds; }

  MuonChamberResidual *chamberResidual(DetId chamberId, int type);
  
  TMatrixDSym covMatrix(DetId chamberId);
  TMatrixDSym corrMatrix(DetId chamberId);
  TMatrixD choleskyCorrMatrix(DetId chamberId);

private:
  TrajectoryStateCombiner m_tsoscomb;

  int m_tracker_numHits;
  double m_tracker_chi2;
  bool m_contains_TIDTEC;

  std::vector<DetId> m_chamberIds;
  std::map<DetId,MuonChamberResidual*> m_dt13, m_dt2, m_csc;
  std::map<DetId,TMatrixDSym> m_trkCovMatrix;

  void addTrkCovMatrix(DetId, TrajectoryStateOnSurface &);

  // pointer to its track
  const reco::Track *track;

  // track muon
  const reco::Muon *muon;
};

#endif // Alignment_MuonAlignmentAlgorithms_MuonResidualsFromTrack_H