OuterTrackerMonitorTrackingParticles.h

Go to the documentation of this file.

#ifndef SiOuterTrackerV_OuterTrackerMonitorTrackingParticles_h
#define SiOuterTrackerV_OuterTrackerMonitorTrackingParticles_h

#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include "DataFormats/Common/interface/DetSetVector.h"
#include "DataFormats/Common/interface/DetSetVectorNew.h"
#include "DataFormats/L1TrackTrigger/interface/TTTypes.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/EDGetToken.h"
#include "SimDataFormats/TrackingAnalysis/interface/TrackingParticle.h"
#include "SimTracker/TrackTriggerAssociation/interface/TTClusterAssociationMap.h"
#include "SimTracker/TrackTriggerAssociation/interface/TTStubAssociationMap.h"
#include "SimTracker/TrackTriggerAssociation/interface/TTTrackAssociationMap.h"
#include <memory>
#include <string>
#include <vector>

class DQMStore;

class OuterTrackerMonitorTrackingParticles : public DQMEDAnalyzer {
public:
  explicit OuterTrackerMonitorTrackingParticles(const edm::ParameterSet &);
  ~OuterTrackerMonitorTrackingParticles() override;
  void analyze(const edm::Event &, const edm::EventSetup &) override;
  void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;
  int Layer(const float R_, const float Z_) const;

  // Tracking particle distributions
  MonitorElement *trackParts_Eta = nullptr;
  MonitorElement *trackParts_Phi = nullptr;
  MonitorElement *trackParts_Pt = nullptr;

  // Plots for correctly matched tracks
  MonitorElement *Track_MatchedChi2 = nullptr;     // Chi2 for only tracks correctly matched to truth level
  MonitorElement *Track_MatchedChi2Red = nullptr;  // Chi2/dof for only tracks correctly matched to truth level

  // pT and eta for efficiency plots
  MonitorElement *tp_pt = nullptr;             // denominator
  MonitorElement *tp_pt_zoom = nullptr;        // denominator
  MonitorElement *tp_eta = nullptr;            // denominator
  MonitorElement *tp_d0 = nullptr;             // denominator
  MonitorElement *tp_VtxR = nullptr;           // denominator (also known as vxy)
  MonitorElement *tp_VtxZ = nullptr;           // denominator
  MonitorElement *match_tp_pt = nullptr;       // numerator
  MonitorElement *match_tp_pt_zoom = nullptr;  // numerator
  MonitorElement *match_tp_eta = nullptr;      // numerator
  MonitorElement *match_tp_d0 = nullptr;       // numerator
  MonitorElement *match_tp_VtxR = nullptr;     // numerator (also known as vxy)
  MonitorElement *match_tp_VtxZ = nullptr;     // numerator

  // 1D intermediate resolution plots (pT and eta)
  MonitorElement *res_eta = nullptr;    // for all eta and pT
  MonitorElement *res_pt = nullptr;     // for all eta and pT
  MonitorElement *res_ptRel = nullptr;  // for all eta and pT (delta(pT)/pT)
  MonitorElement *respt_eta0to0p7_pt2to3 = nullptr;
  MonitorElement *respt_eta0p7to1_pt2to3 = nullptr;
  MonitorElement *respt_eta1to1p2_pt2to3 = nullptr;
  MonitorElement *respt_eta1p2to1p6_pt2to3 = nullptr;
  MonitorElement *respt_eta1p6to2_pt2to3 = nullptr;
  MonitorElement *respt_eta2to2p4_pt2to3 = nullptr;
  MonitorElement *respt_eta0to0p7_pt3to8 = nullptr;
  MonitorElement *respt_eta0p7to1_pt3to8 = nullptr;
  MonitorElement *respt_eta1to1p2_pt3to8 = nullptr;
  MonitorElement *respt_eta1p2to1p6_pt3to8 = nullptr;
  MonitorElement *respt_eta1p6to2_pt3to8 = nullptr;
  MonitorElement *respt_eta2to2p4_pt3to8 = nullptr;
  MonitorElement *respt_eta0to0p7_pt8toInf = nullptr;
  MonitorElement *respt_eta0p7to1_pt8toInf = nullptr;
  MonitorElement *respt_eta1to1p2_pt8toInf = nullptr;
  MonitorElement *respt_eta1p2to1p6_pt8toInf = nullptr;
  MonitorElement *respt_eta1p6to2_pt8toInf = nullptr;
  MonitorElement *respt_eta2to2p4_pt8toInf = nullptr;
  MonitorElement *reseta_eta0to0p7 = nullptr;
  MonitorElement *reseta_eta0p7to1 = nullptr;
  MonitorElement *reseta_eta1to1p2 = nullptr;
  MonitorElement *reseta_eta1p2to1p6 = nullptr;
  MonitorElement *reseta_eta1p6to2 = nullptr;
  MonitorElement *reseta_eta2to2p4 = nullptr;
  MonitorElement *resphi_eta0to0p7 = nullptr;
  MonitorElement *resphi_eta0p7to1 = nullptr;
  MonitorElement *resphi_eta1to1p2 = nullptr;
  MonitorElement *resphi_eta1p2to1p6 = nullptr;
  MonitorElement *resphi_eta1p6to2 = nullptr;
  MonitorElement *resphi_eta2to2p4 = nullptr;
  MonitorElement *resVtxZ_eta0to0p7 = nullptr;
  MonitorElement *resVtxZ_eta0p7to1 = nullptr;
  MonitorElement *resVtxZ_eta1to1p2 = nullptr;
  MonitorElement *resVtxZ_eta1p2to1p6 = nullptr;
  MonitorElement *resVtxZ_eta1p6to2 = nullptr;
  MonitorElement *resVtxZ_eta2to2p4 = nullptr;

  // For d0
  MonitorElement *resd0_eta0to0p7 = nullptr;
  MonitorElement *resd0_eta0p7to1 = nullptr;
  MonitorElement *resd0_eta1to1p2 = nullptr;
  MonitorElement *resd0_eta1p2to1p6 = nullptr;
  MonitorElement *resd0_eta1p6to2 = nullptr;
  MonitorElement *resd0_eta2to2p4 = nullptr;

private:
  edm::ParameterSet conf_;
  edm::EDGetTokenT<std::vector<TrackingParticle>> trackingParticleToken_;
  edm::EDGetTokenT<edmNew::DetSetVector<TTStub<Ref_Phase2TrackerDigi_>>> ttStubToken_;
  edm::EDGetTokenT<std::vector<TTTrack<Ref_Phase2TrackerDigi_>>> ttTrackToken_;
  edm::EDGetTokenT<TTClusterAssociationMap<Ref_Phase2TrackerDigi_>>
      ttClusterMCTruthToken_;  // MC truth association map for clusters
  edm::EDGetTokenT<TTStubAssociationMap<Ref_Phase2TrackerDigi_>>
      ttStubMCTruthToken_;  // MC truth association map for stubs
  edm::EDGetTokenT<TTTrackAssociationMap<Ref_Phase2TrackerDigi_>>
      ttTrackMCTruthToken_;  // MC truth association map for tracks
  int L1Tk_nPar;
  int L1Tk_minNStub;
  double L1Tk_maxChi2;
  double L1Tk_maxChi2dof;
  int TP_minNStub;
  double TP_minPt;
  double TP_maxPt;
  double TP_maxEta;
  double TP_maxVtxZ;
  int TP_select_eventid;
  std::string topFolderName_;

  // The following adds a variable that is the number of layers hit for each
  // tracking particle
  struct TpStruct {
    int TpId;
    std::vector<bool> layer;
    int Nlayers() {  // Counts how many layers are set to "true" for their hit status
      int layers = 0;
      for (unsigned l = 0; l < layer.size(); ++l)
        if (layer[l])
          ++layers;
      return layers;
    }
  };
};
#endif