TrackEfficiencyMonitor.h

Go to the documentation of this file.

#ifndef TrackEfficiencyMonitor_H
#define TrackEfficiencyMonitor_H
// -*- C++ -*-
//
// Package:    TrackEfficiencyMonitor
// Class:      TrackEfficiencyMonitor
//
//  Original Author:  Jeremy Andrea
// Insertion in DQM:  Anne-Catherine Le Bihan
//          Created:  Thu 28 22:45:30 CEST 2008
 
#include <memory>
#include <fstream>
#include "FWCore/Utilities/interface/EDGetToken.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "DQMServices/Core/interface/DQMStore.h"
 
#include <DQMServices/Core/interface/DQMEDAnalyzer.h>
 
#include "RecoMuon/GlobalTrackingTools/interface/DirectTrackerNavigation.h"
#include "RecoMuon/TrackingTools/interface/MuonServiceProxy.h"
#include "RecoTracker/MeasurementDet/interface/MeasurementTracker.h"
#include "TrackingTools/TransientTrack/interface/TransientTrackBuilder.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/MuonReco/interface/MuonSelectors.h"
 
namespace reco {
  class TransientTrack;
}
class NavigationSchool;
 
class TrackEfficiencyMonitor : public DQMEDAnalyzer {
public:
  typedef reco::Track Track;
  typedef reco::TrackCollection TrackCollection;
  explicit TrackEfficiencyMonitor(const edm::ParameterSet&);
  ~TrackEfficiencyMonitor() override;
  void analyze(const edm::Event&, const edm::EventSetup&) override;
 
  void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;
 
  enum SemiCylinder { Up, Down };
  std::pair<TrajectoryStateOnSurface, const DetLayer*> findNextLayer(
      TrajectoryStateOnSurface startTSOS, const std::vector<const DetLayer*>& trackCompatibleLayers, bool isUpMuon);
  SemiCylinder checkSemiCylinder(const Track&);
  void testTrackerTracks(edm::Handle<TrackCollection> tkTracks,
                         edm::Handle<TrackCollection> staTracks,
                         const NavigationSchool& navigationSchool);
  void testSTATracks(edm::Handle<TrackCollection> tkTracks, edm::Handle<TrackCollection> staTracks);
  bool trackerAcceptance(TrajectoryStateOnSurface theTSOS, double theRadius, double theMaxZ);
  int compatibleLayers(const NavigationSchool& navigationSchool, TrajectoryStateOnSurface theTSOS);
 
private:
  // ----------member data ---------------------------
 
  std::string histname;  //for naming the histograms
 
  DQMStore* dqmStore_;
  edm::ParameterSet conf_;
 
  double theRadius_;
  double theMaxZ_;
  bool isBFieldOff_;
  bool trackEfficiency_;  //1 if one wants to measure the tracking efficiency
                          //0 if one wants to measure the muon reco efficiency
 
  edm::InputTag theTKTracksLabel_;
  edm::InputTag theSTATracksLabel_;
  edm::EDGetTokenT<reco::TrackCollection> theTKTracksToken_;
  edm::EDGetTokenT<reco::TrackCollection> theSTATracksToken_;
 
  int failedToPropagate;
  int nCompatibleLayers;
  bool findDetLayer;
 
  MonitorElement* muonX;
  MonitorElement* muonY;
  MonitorElement* muonZ;
  MonitorElement* muonEta;
  MonitorElement* muonPhi;
  MonitorElement* muonD0;
  MonitorElement* muonCompatibleLayers;
 
  MonitorElement* trackX;
  MonitorElement* trackY;
  MonitorElement* trackZ;
  MonitorElement* trackEta;
  MonitorElement* trackPhi;
  MonitorElement* trackD0;
  MonitorElement* trackCompatibleLayers;
 
  MonitorElement* deltaX;
  MonitorElement* deltaY;
  MonitorElement* signDeltaX;
  MonitorElement* signDeltaY;
  MonitorElement* GlobalMuonPtEtaPhiLowPt;
  MonitorElement* StandaloneMuonPtEtaPhiLowPt;
  MonitorElement* GlobalMuonPtEtaPhiHighPt;
  MonitorElement* StandaloneMuonPtEtaPhiHighPt;
  const DirectTrackerNavigation* theNavigation;
  MuonServiceProxy* theMuonServiceProxy;
  edm::EDGetTokenT<edm::View<reco::Muon> > muonToken_;
  edm::ESHandle<GeometricSearchTracker> theGeometricSearchTracker;
  edm::ESHandle<Propagator> thePropagator;
  edm::ESHandle<Propagator> thePropagatorCyl;
  edm::ESHandle<TransientTrackBuilder> theTTrackBuilder;
  edm::ESHandle<GeometricSearchTracker> theTracker;
  edm::ESHandle<MagneticField> bField;
 
  edm::ESHandle<MeasurementTracker> measurementTrackerHandle;
};
#endif