drof macro

Signed-off-by: Felix Schlepper <felix.schlepper@cern.ch>

Author: f3sch

Detectors/ITSMFT/ITS/macros/test/CheckDROF.C

@@ -707,39 +707,39 @@ void CheckDROF(bool plot = false, bool write = false, const std::string& tracfil
       }
     }
 
-    std::vector<VertexInfo> vertexInfo;
-    std::cout << "** Creating vertices/particles correspondance ... " << std::flush;
-    for (int frame = 0; frame < recTree->GetEntriesFast(); frame++) { // Vertices frames
-      if (!recTree->GetEvent(frame)) {
-        continue;
-      }
-      int contLabIdx{0}; // contributor labels are stored as flat vector
-      for (size_t iRecord{0}; iRecord < recVerROFArr->size(); ++iRecord) {
-        auto& rec = recVerROFArr->at(iRecord);
-        auto verStartIdx = rec.getFirstEntry(), verSize = rec.getNEntries();
-        for (int iVertex{rec.getFirstEntry()}; iVertex < verStartIdx + verSize; ++iVertex) {
-          auto& info = vertexInfo.emplace_back();
-          info.vertex = recVerArr->at(iVertex);
-          info.mainLabel = recVerLabelsArr->at(contLabIdx);
-          info.purity = recVerPurityArr->at(contLabIdx);
-          info.event = info.mainLabel.getEventID();
-          ++contLabIdx;
-          if (info.mainLabel.isSet()) {
-            const auto& ir = irs[info.event];
-            // LOGP(info, "iROF={} {} to {}", iRecord, info.mainLabel.asString(), ir.asString());
-            if (!ir.isDummy()) {
-              info.bcInROF = (ir.toLong() - roFrameBiasInBC) % roFrameLengthInBC;
-              info.rofId = (ir.toLong() - roFrameBiasInBC) / roFrameLengthInBC;
+    //////////////////////
+    // DROF Vertices
+    if constexpr (false) {
+      std::vector<VertexInfo> vertexInfo;
+      std::cout << "** Creating vertices/particles correspondance ... " << std::flush;
+      for (int frame = 0; frame < recTree->GetEntriesFast(); frame++) { // Vertices frames
+        if (!recTree->GetEvent(frame)) {
+          continue;
+        }
+        int contLabIdx{0}; // contributor labels are stored as flat vector
+        for (size_t iRecord{0}; iRecord < recVerROFArr->size(); ++iRecord) {
+          auto& rec = recVerROFArr->at(iRecord);
+          auto verStartIdx = rec.getFirstEntry(), verSize = rec.getNEntries();
+          for (int iVertex{rec.getFirstEntry()}; iVertex < verStartIdx + verSize; ++iVertex) {
+            auto& info = vertexInfo.emplace_back();
+            info.vertex = recVerArr->at(iVertex);
+            info.mainLabel = recVerLabelsArr->at(contLabIdx);
+            info.purity = recVerPurityArr->at(contLabIdx);
+            info.event = info.mainLabel.getEventID();
+            ++contLabIdx;
+            if (info.mainLabel.isSet()) {
+              const auto& ir = irs[info.event];
+              // LOGP(info, "iROF={} {} to {}", iRecord, info.mainLabel.asString(), ir.asString());
+              if (!ir.isDummy()) {
+                info.bcInROF = (ir.toLong() - roFrameBiasInBC) % roFrameLengthInBC;
+                info.rofId = (ir.toLong() - roFrameBiasInBC) / roFrameLengthInBC;
+              }
             }
           }
         }
       }
-    }
-    std::cout << "done." << std::endl;
+      std::cout << "done." << std::endl;
 
-    //////////////////////
-    // DROF Vertices
-    {
       auto hMCVtxZ = new TH1F("hMCVtxZ", "MC Vertex;Z", 50, -16, 16);
       auto hReVtxZ = new TH1F("hRecoVtxZ", "Reco Vertex;Z", 50, -16, 16);
 
@@ -891,7 +891,7 @@ void CheckDROF(bool plot = false, bool write = false, const std::string& tracfil
       fOut->Close();
     }
   } else {
-    auto fWO = TFile::Open("checkDROF_wo.root");
+    auto fWO = TFile::Open("checkDROF_w_new.root");
     auto f = TFile::Open("checkDROF_w.root");
     plotHistos(fWO, f, "");
     plotHistos(fWO, f, "Val");