add more monitoring histograms to shortened track resolution

Author: mmusich

DQM/TrackingMonitorSource/plugins/ShortenedTrackResolution.cc

@@ -39,7 +39,15 @@ class ShortenedTrackResolution : public DQMEDAnalyzer {
   const edm::EDGetTokenT<std::vector<reco::Track>> tracksToken_;
   const std::vector<edm::EDGetTokenT<std::vector<reco::Track>>> tracksRerecoToken_;
 
+  std::vector<MonitorElement *> histsPtRatioAll_;
+  std::vector<MonitorElement *> histsPtDiffAll_;
+  std::vector<MonitorElement *> histsEtaDiffAll_;
+  std::vector<MonitorElement *> histsPhiDiffAll_;
+  std::vector<MonitorElement *> histsPtRatioVsDeltaRAll_;
+  std::vector<MonitorElement *> histsDeltaPtOverPtAll_;
   std::vector<MonitorElement *> histsPtAll_;
+  std::vector<MonitorElement *> histsNhitsAll_;
+  std::vector<MonitorElement *> histsDeltaRAll_;
 };
 
 // -----------------------------
@@ -58,25 +66,103 @@ ShortenedTrackResolution::ShortenedTrackResolution(const edm::ParameterSet &ps)
       tracksToken_(consumes<std::vector<reco::Track>>(tracksTag_)),
       tracksRerecoToken_(edm::vector_transform(
           tracksRerecoTag_, [this](edm::InputTag const &tag) { return consumes<std::vector<reco::Track>>(tag); })) {
+  histsPtRatioAll_.clear();
+  histsPtDiffAll_.clear();
+  histsEtaDiffAll_.clear();
+  histsPhiDiffAll_.clear();
+  histsPtRatioVsDeltaRAll_.clear();
+  histsDeltaPtOverPtAll_.clear();
   histsPtAll_.clear();
+  histsNhitsAll_.clear();
+  histsDeltaRAll_.clear();
+
+  const size_t n = hitsRemain_.size();
+  histsPtRatioAll_.reserve(n);
+  histsPtDiffAll_.reserve(n);
+  histsEtaDiffAll_.reserve(n);
+  histsPhiDiffAll_.reserve(n);
+  histsPtRatioVsDeltaRAll_.reserve(n);
+  histsDeltaPtOverPtAll_.reserve(n);
+  histsPtAll_.reserve(n);
+  histsNhitsAll_.reserve(n);
+  histsDeltaRAll_.reserve(n);
 }
 
 //__________________________________________________________________________________
 void ShortenedTrackResolution::bookHistograms(DQMStore::IBooker &iBook,
                                               edm::Run const &iRun,
                                               edm::EventSetup const &iSetup) {
-  std::string currentFolder = folderName_ + "/";
+  auto book1D = [&](const std::string &name, const std::string &title, int bins, double min, double max) {
+    return iBook.book1D(name.c_str(), title.c_str(), bins, min, max);
+  };
+
+  auto book2D = [&](const std::string &name,
+                    const std::string &title,
+                    int binsX,
+                    double minX,
+                    double maxX,
+                    int binsY,
+                    double minY,
+                    double maxY) {
+    return iBook.book2D(name.c_str(), title.c_str(), binsX, minX, maxX, binsY, minY, maxY);
+  };
+
+  std::string currentFolder = folderName_ + "/Resolutions";
   iBook.setCurrentFolder(currentFolder);
 
-  for (int i = 0; i < int(hitsRemain_.size()); ++i) {
-    histsPtAll_.push_back(iBook.book1D(
-        fmt::sprintf("trackPtRatio_%s", hitsRemain_[i]).c_str(),
-        fmt::sprintf("Short Track p_{T} / Full Track p_{T} - %s layers;p_{T}^{short}/p_{T}^{full};n. tracks",
-                     hitsRemain_[i])
-            .c_str(),
-        101,
-        -0.05,
-        2.05));
+  for (const auto &label : hitsRemain_) {
+    std::string name, title;
+
+    name = fmt::sprintf("trackPtRatio_%s", label);
+    title =
+        fmt::sprintf("Short Track p_{T} / Full Track p_{T} - %s layers;p_{T}^{short}/p_{T}^{full};n. tracks", label);
+    histsPtRatioAll_.push_back(book1D(name, title, 100, 0.5, 1.5));
+
+    name = fmt::sprintf("trackPtDiff_%s", label);
+    title = fmt::sprintf(
+        "Short Track p_{T} - Full Track p_{T} - %s layers;p_{T}^{short} - p_{T}^{full} [GeV];n. tracks", label);
+    histsPtDiffAll_.push_back(book1D(name, title, 100, -10., 10.));
+
+    name = fmt::sprintf("trackEtaDiff_%s", label);
+    title = fmt::sprintf("Short Track #eta - Full Track #eta - %s layers;#eta^{short} - #eta^{full};n. tracks", label);
+    histsEtaDiffAll_.push_back(book1D(name, title, 100, -0.001, 0.001));
+
+    name = fmt::sprintf("trackPhiDiff_%s", label);
+    title = fmt::sprintf("Short Track #phi - Full Track #phi - %s layers;#phi^{short} - #phi^{full};n. tracks", label);
+    histsPhiDiffAll_.push_back(book1D(name, title, 100, -0.001, 0.001));
+
+    name = fmt::sprintf("trackPtRatioVsDeltaR_%s", label);
+    title = fmt::sprintf(
+        "Short Track p_{T} / Full Track p_{T} - %s layers vs "
+        "#DeltaR;#DeltaR(short,full);p_{T}^{short}/p_{T}^{full} [GeV];n. tracks",
+        label);
+    histsPtRatioVsDeltaRAll_.push_back(book2D(name, title, 100, 0., 0.01, 101, -0.05, 2.05));
+
+    name = fmt::sprintf("trackDeltaPtOverPt_%s", label);
+    title = fmt::sprintf(
+        "Short Track p_{T} - Full Track p_{T} / Full Track p_{T} - %s layers;"
+        "p_{T}^{short} - p_{T}^{full} / p^{full}_{T};n. tracks",
+        label);
+    histsDeltaPtOverPtAll_.push_back(book1D(name, title, 101, -5., 5.));
+  }
+
+  currentFolder = folderName_ + "/TrackProperties";
+  iBook.setCurrentFolder(currentFolder);
+
+  for (const auto &label : hitsRemain_) {
+    std::string name, title;
+
+    name = fmt::sprintf("trackPt_%s", label);
+    title = fmt::sprintf("Short Track p_{T} - %s layers;p_{T}^{short} [GeV];n. tracks", label);
+    histsPtAll_.push_back(book1D(name, title, 100, 0., 100.));
+
+    name = fmt::sprintf("trackNhits_%s", label);
+    title = fmt::sprintf("Short Track n. hits - %s layers;n. hits per track;n. tracks", label);
+    histsNhitsAll_.push_back(book1D(name, title, 20, -0.5, 19.5));
+
+    name = fmt::sprintf("trackDeltaR_%s", label);
+    title = fmt::sprintf("Short Track / Full Track #DeltaR - %s layers;#DeltaR(short,full);n. tracks", label);
+    histsDeltaRAll_.push_back(book1D(name, title, 100, 0., 0.005));
   }
 }
 
@@ -110,7 +196,15 @@ void ShortenedTrackResolution::analyze(edm::Event const &iEvent, edm::EventSetup
         if (deltaR < maxDr_) {
           if (track_rereco.pt() >= minTracksPt_ && track_rereco.pt() <= maxTracksPt_ &&
               std::abs(track_rereco.eta()) >= minTracksEta_ && std::abs(track_rereco.eta()) <= maxTracksEta_) {
-            histsPtAll_[i]->Fill(1.0 * track_rereco.pt() / track.pt());
+            histsPtRatioAll_[i]->Fill(1.0 * track_rereco.pt() / track.pt());
+            histsPtDiffAll_[i]->Fill(track_rereco.pt() - track.pt());
+            histsDeltaPtOverPtAll_[i]->Fill((track_rereco.pt() - track.pt()) / track.pt());
+            histsEtaDiffAll_[i]->Fill(track_rereco.eta() - track.eta());
+            histsPhiDiffAll_[i]->Fill(track_rereco.phi() - track.phi());
+            histsPtRatioVsDeltaRAll_[i]->Fill(deltaR, track_rereco.pt() / track.pt());
+            histsPtAll_[i]->Fill(track_rereco.pt());
+            histsNhitsAll_[i]->Fill(track_rereco.numberOfValidHits());
+            histsDeltaRAll_[i]->Fill(deltaR);
           }
         }
       }