ITS3: add TopoDict print

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

Author: f3sch

Detectors/Upgrades/ITS3/base/include/ITS3Base/SpecsV2.h

@@ -135,7 +135,6 @@ constexpr std::array<double, nLayers> radiiInner{radii[0] - silicon::thicknessIn
 constexpr std::array<double, nLayers> radiiOuter{radii[0] + silicon::thicknessOut, radii[1] + silicon::thicknessOut, radii[2] + silicon::thicknessOut};              // outer silicon radius
 constexpr std::array<double, nLayers> radiiMiddle{(radiiInner[0] + radiiOuter[0]) / 2., (radiiInner[1] + radiiOuter[1]) / 2., (radiiInner[2] + radiiOuter[2]) / 2.}; // middle silicon radius
 constexpr double nominalYShift{-metalstack::thickness / 2.};                                                                                                         // shift to position in silicion volume to the chip volume (silicon+metalstack)
-/*constexpr double nominalYShift{0}; // shift to position in silicion volume to the chip volume (silicon+metalstack)*/
 
 // extra information of pixels and their response functions
 namespace pixelarray::pixels

Detectors/Upgrades/ITS3/macros/test/CheckClustersITS3.C

@@ -42,7 +42,7 @@
 void CheckClustersITS3(const std::string& clusfile = "o2clus_its.root",
                        const std::string& hitfile = "o2sim_HitsIT3.root",
                        const std::string& inputGeom = "",
-                       std::string dictfile = "../ccdb/IT3/Calib/ClusterDictionary/snapshot.root",
+                       std::string dictfile = "./ccdb/IT3/Calib/ClusterDictionary/snapshot.root",
                        bool batch = false)
 {
   gROOT->SetBatch(batch);
@@ -63,7 +63,8 @@ void CheckClustersITS3(const std::string& clusfile = "o2clus_its.root",
   std::vector<HitVec*> hitVecPool;
   std::vector<MC2HITS_map> mc2hitVec;
 
-  ULong_t cPattValid{0}, cPattInvalid{0}, cLabelInvalid{0}, cNoMC{0};
+  ULong_t cPattValidIB{0}, cPattInvalidIB{0}, cLabelInvalidIB{0}, cNoMCIB{0};
+  ULong_t cPattValidOB{0}, cPattInvalidOB{0}, cLabelInvalidOB{0}, cNoMCOB{0};
 
   TFile fout("CheckClusters.root", "recreate");
   TNtuple nt("ntc", "cluster ntuple", "ev:lab:hlx:hlz:hgx:hgz:tx:tz:cgx:cgy:cgz:clx:cly:clz:dx:dy:dz:ex:ez:patid:rof:npx:id:eta:row:col:lay");
@@ -103,6 +104,7 @@ void CheckClustersITS3(const std::string& clusfile = "o2clus_its.root",
   } else {
     LOG(info) << "Running without dictionary !";
   }
+  dict.print();
 
   // ROFrecords
   std::vector<ROFRec> rofRecVec, *rofRecVecP = &rofRecVec;
@@ -179,22 +181,22 @@ void CheckClustersITS3(const std::string& clusfile = "o2clus_its.root",
         o2::itsmft::ClusterPattern patt(pattIt);
         locC = dict.getClusterCoordinates(cluster, patt, false);
         LOGP(debug, "I am invalid and I am on chip {}", chipID);
-        ++cPattInvalid;
+        (isIB) ? ++cPattInvalidIB : ++cPattInvalidOB;
         continue;
       } else {
         locC = dict.getClusterCoordinates(cluster);
         errX = dict.getErrX(pattID, isIB);
         errZ = dict.getErrZ(pattID, isIB);
         npix = dict.getNpixels(pattID, isIB);
-        ++cPattValid;
+        (isIB) ? ++cPattValidIB : ++cPattValidOB;
       }
 
       // Transformation to the local --> global
       auto gloC = gman->getMatrixL2G(chipID)(locC);
       const auto& lab = (clusLabArr->getLabels(clEntry))[0];
 
       if (!lab.isValid()) {
-        ++cLabelInvalid;
+        (isIB) ? ++cLabelInvalidIB : ++cLabelInvalidOB;
         continue;
       }
 
@@ -206,7 +208,7 @@ void CheckClustersITS3(const std::string& clusfile = "o2clus_its.root",
       auto hitEntry = mc2hit.find(key);
       if (hitEntry == mc2hit.end()) {
         LOG(debug) << "Failed to find MC hit entry for Tr" << trID << " chipID" << chipID;
-        ++cNoMC;
+        (isIB) ? ++cNoMCIB : ++cNoMCOB;
         continue;
       }
       const auto& hit = (*hitArray)[hitEntry->second];
@@ -261,8 +263,10 @@ void CheckClustersITS3(const std::string& clusfile = "o2clus_its.root",
     }
   }
 
-  LOGP(info, "There were {} valid PatternIDs and {} ({:.1f}%) invalid ones", cPattValid, cPattInvalid, ((float)cPattInvalid / (float)(cPattInvalid + cPattValid)) * 100);
-  LOGP(info, "There were {} invalid Labels and {} with No MC Hit information ", cLabelInvalid, cNoMC);
+  LOGP(info, "IB {} valid PatternIDs and {} ({:.1f}%) invalid ones", cPattValidIB, cPattInvalidIB, ((float)cPattInvalidIB / (float)(cPattInvalidIB + cPattValidIB)) * 100);
+  LOGP(info, "IB {} invalid Labels and {} with No MC Hit information ", cLabelInvalidIB, cNoMCIB);
+  LOGP(info, "OB {} valid PatternIDs and {} ({:.1f}%) invalid ones", cPattValidOB, cPattInvalidOB, ((float)cPattInvalidOB / (float)(cPattInvalidOB + cPattValidOB)) * 100);
+  LOGP(info, "OB {} invalid Labels and {} with No MC Hit information ", cLabelInvalidOB, cNoMCOB);
 
   auto canvCgXCgY = new TCanvas("canvCgXCgY", "", 1600, 1600);
   canvCgXCgY->Divide(2, 2);

Detectors/Upgrades/ITS3/reconstruction/include/ITS3Reconstruction/TopologyDictionary.h

@@ -59,6 +59,7 @@ class TopologyDictionary
   /// Reads the dictionary from a binary file
   void readBinaryFile(const std::string& fileName);
   void readFromFile(const std::string& fileName);
+  void print() const noexcept;
 
   /// Returns the x position of the COG for the n_th element
   [[nodiscard]] float getXCOG(int n, bool IB) const

Detectors/Upgrades/ITS3/reconstruction/src/TopologyDictionary.cxx

@@ -40,6 +40,15 @@ TopologyDictionary::TopologyDictionary(const std::string& fileName)
   readFromFile(fileName);
 }
 
+void TopologyDictionary::print() const noexcept
+{
+  LOG(info) << "ITS3 TopologyDictionary";
+  LOG(info) << "InnerBarrel";
+  mDataIB.print();
+  LOG(info) << "OuterBarrel";
+  mDataOB.print();
+}
+
 void TopologyDictionary::reset() noexcept
 {
   mDataIB.mSmallTopologiesLUT.fill(-1);