Add dumping of deltas on file in CreateDictionaries

Author: mconcas

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

@@ -32,7 +32,8 @@
 #include <unordered_map>
 #endif
 
-void CreateDictionaries(std::string clusfile = "o2clus_its.root",
+void CreateDictionaries(bool saveDeltas = false,
+                        std::string clusfile = "o2clus_its.root",
                         std::string hitfile = "o2sim_HitsITS.root",
                         std::string collContextfile = "collisioncontext.root",
                         std::string inputGeom = "",
@@ -56,6 +57,14 @@ void CreateDictionaries(std::string clusfile = "o2clus_its.root",
   std::unordered_map<int, int> hadronicMCMap;            // mapping from MC event entry to hadronic event ID
   std::vector<HitVec*> hitVecPool;
   std::vector<MC2HITS_map> mc2hitVec;
+
+  TFile* fout = nullptr;
+  TNtuple* nt = nullptr;
+  if (saveDeltas) {
+    fout = TFile::Open("CreateDictionaries.root", "recreate");
+    nt = new TNtuple("nt", "hashes ntuple", "hash:dx:dz");
+  }
+
   const o2::steer::DigitizationContext* digContext = nullptr;
   TStopwatch sw;
   sw.Start();
@@ -214,6 +223,9 @@ void CreateDictionaries(std::string clusfile = "o2clus_its.root",
               const auto locC = o2::itsmft::TopologyDictionary::getClusterCoordinates(cluster, pattern, false);
               dX = locH.X() - locC.X();
               dZ = locH.Z() - locC.Z();
+              if (saveDeltas) {
+                nt->Fill(topology.getHash(), dX, dZ);
+              }
               if (checkOutliers > 0.) {
                 if (std::abs(dX) > topology.getRowSpan() * o2::itsmft::SegmentationAlpide::PitchRow * checkOutliers ||
                     std::abs(dZ) > topology.getColumnSpan() * o2::itsmft::SegmentationAlpide::PitchCol * checkOutliers) { // ignore outlier
@@ -299,4 +311,8 @@ void CreateDictionaries(std::string clusfile = "o2clus_its.root",
     sw.Stop();
     sw.Print();
   }
+  if (saveDeltas) {
+    fout->cd();
+    nt->Write();
+  }
 }