d8/d6e/OuterTrackerMonitorTrackingParticles_8h_source.html

#ifndef SiOuterTrackerV_OuterTrackerMonitorTrackingParticles_h

#define SiOuterTrackerV_OuterTrackerMonitorTrackingParticles_h


#include "DQMServices/Core/interface/DQMEDAnalyzer.h"

#include "DQMServices/Core/interface/DQMStore.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 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;


  // 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<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;

  int TP_minNLayersStub;

  double TP_minPt;

  double TP_maxPt;

  double TP_maxEta;

  double TP_maxVtxZ;

  int TP_select_eventid;

  std::string topFolderName_;

};

#endif