Optimize binning for vertex time and z distributions

cms-sw · Oct 22, 2024 · 397bbc9 · 397bbc9
1 parent 442eda9
commit 397bbc9
Showing 1 changed file with 8 additions and 4 deletions.
diff --git a/Validation/MtdValidation/plugins/Primary4DVertexValidation.cc b/Validation/MtdValidation/plugins/Primary4DVertexValidation.cc
@@ -373,6 +373,7 @@ class Primary4DVertexValidation : public DQMEDAnalyzer {
   MonitorElement* meSimVerNumber_;
   MonitorElement* meSimPVZ_;
   MonitorElement* meSimPVT_;
+  MonitorElement* meSimPVTvsZ_;
 
   MonitorElement* meVtxTrackMult_;
   MonitorElement* meVtxTrackW_;
@@ -698,10 +699,12 @@ void Primary4DVertexValidation::bookHistograms(DQMStore::IBooker& ibook,
                    20);
   meRecVerNumber_ = ibook.book1D("RecVerNumber", "RECO Vertex Number: Number of vertices", 50, 0, 250);
   meSimVerNumber_ = ibook.book1D("SimVerNumber", "SIM Vertex Number: Number of vertices", 50, 0, 250);
-  meRecPVZ_ = ibook.book1D("recPVZ", "#Rec vertices/mm", 400, -20., 20.);
-  meRecPVT_ = ibook.book1D("recPVT", "#Rec vertices/10 ps", 200, -1., 1.);
-  meSimPVZ_ = ibook.book1D("simPVZ", "#Sim vertices/mm", 400, -20., 20.);
-  meSimPVT_ = ibook.book1D("simPVT", "#Sim vertices/10 ps", 200, -1., 1.);
+  meRecPVZ_ = ibook.book1D("recPVZ", "#Rec vertices/10 mm", 30, -15., 15.);
+  meRecPVT_ = ibook.book1D("recPVT", "#Rec vertices/50 ps", 30, -0.75, 0.75);
+  meSimPVZ_ = ibook.book1D("simPVZ", "#Sim vertices/10 mm", 30, -15., 15.);
+  meSimPVT_ = ibook.book1D("simPVT", "#Sim vertices/50 ps", 30, -0.75, 0.75);
+  meSimPVTvsZ_ = ibook.bookProfile("simPVTvsZ", "PV Time vs Z", 30, -15., 15., 30, -0.75, 0.75);
+
 
   meVtxTrackMult_ = ibook.book1D("VtxTrackMult", "Log10(Vertex track multiplicity)", 80, 0.5, 2.5);
   meVtxTrackW_ = ibook.book1D("VtxTrackW", "Vertex track weight (all)", 50, 0., 1.);
@@ -2641,6 +2644,7 @@ void Primary4DVertexValidation::analyze(const edm::Event& iEvent, const edm::Eve
   for (unsigned int is = 0; is < simpv.size(); is++) {
     meSimPVZ_->Fill(simpv.at(is).z);
     meSimPVT_->Fill(simpv.at(is).t * simUnit_);
+    meSimPVTvsZ_->Fill(simpv.at(is).z, simpv.at(is).t * simUnit_);
     if (is == 0 && optionalPlots_) {
       meSimPosInSimOrigCollection_->Fill(simpv.at(is).OriginalIndex);
     }