Add time info to clusters

Author: Marcello Di Costanzo

Detectors/Upgrades/ALICE3/IOTOF/DataFormatsIOTOF/include/DataFormatsIOTOF/Cluster.h

@@ -30,6 +30,7 @@ struct Cluster {
   float xCoord = 0.f;
   float yCoord = 0.f;
   float zCoord = 0.f;
+  float time = 0.f;
 
   std::string asString() const;
 

Detectors/Upgrades/ALICE3/IOTOF/macros/CheckClustersIOTOF.C

@@ -84,7 +84,7 @@ void CheckClustersIOTOF(std::string digiFilePath = "tf3digits.root", std::string
   o2::dataformats::ConstMCTruthContainer<o2::MCCompLabel> labels;
   clsPlabelsArr->copyandflatten(labels);
 
-  auto clsTuple = new TNtuple("clsTuple", "clsTuple", "chip_id:x:y:z:subdet_id:row:col");
+  auto clsTuple = new TNtuple("clsTuple", "clsTuple", "chip_id:x:y:z:subdet_id:row:col:time");
   clsTuple->SetDirectory(nullptr);
 
   TH1F* histXCoordCls = new TH1F("histXCoordCls", "histXCoordCls", 8000, -100, 100);
@@ -144,14 +144,15 @@ void CheckClustersIOTOF(std::string digiFilePath = "tf3digits.root", std::string
 
       o2::math_utils::Point3D<float> localClsCoords(x, y, z); // local Digit
       const auto globalClsCoords = tofGeo->getMatrixL2G(chipID)(localClsCoords); // convert to global
-      std::cout << "Cluster " << iCls << ": chipID = " << chipID << ", X=" << globalClsCoords.x() << ", Y=" << globalClsCoords.y() << ", Z=" << globalClsCoords.z() << std::endl;
+      std::cout << "Cluster " << iCls << ": chipID = " << chipID << ", X=" << globalClsCoords.x() << ", Y=" << globalClsCoords.y() << ", Z=" << globalClsCoords.z() << ", time=" << (*clsArray)[iCls].time << std::endl;
       clsTuple->Fill((*clsArray)[iCls].chipID,
                      globalClsCoords.x(),
                      globalClsCoords.y(),
                      globalClsCoords.z(),
                      (*clsArray)[iCls].subDetID,
                      (*clsArray)[iCls].row,
-                     (*clsArray)[iCls].col);
+                     (*clsArray)[iCls].col,
+                     (*clsArray)[iCls].time);
       histXCoordCls->Fill(globalClsCoords.x());
       histYCoordCls->Fill(globalClsCoords.y());
       histZCoordCls->Fill(globalClsCoords.z());

Detectors/Upgrades/ALICE3/IOTOF/reconstruction/src/Clusterer.cxx

@@ -30,7 +30,7 @@ o2::math_utils::Point3D<float> Clusterer::getClusterGlobalCoordinates(const Clus
 
   Segmentation::Instance()->detectorToLocal(cluster.row, cluster.col, coords, cluster.subDetID);
   LOG(info) << "[Clusterer] Cluster local coordinates: x=" << coords.x() << ", y=" << coords.y() << ", z=" << coords.z();
-  GeometryTGeo::Instance()->getMatrixL2G(cluster.subDetID)(coords);
+  GeometryTGeo::Instance()->getMatrixL2G(cluster.chipID)(coords);
 
   LOG(info) << "[Clusterer] Cluster global coordinates: x=" << coords.x() << ", y=" << coords.y() << ", z=" << coords.z();
   return coords;
@@ -183,10 +183,11 @@ void Clusterer::ClustererThread::finishChipSingleHitFast(gsl::span<const Digit>
                                                          ClusterTruth* labelsClusPtr,
                                                          GeometryTGeo* geom)
 {
-  const auto& d = digits[digitIdx];
-  const uint16_t chipID = d.getChipIndex();
-  const uint16_t row = d.getRow();
-  const uint16_t col = d.getColumn();
+  const auto& digit = digits[digitIdx];
+  const uint16_t chipID = digit.getChipIndex();
+  const uint16_t row    = digit.getRow();
+  const uint16_t col    = digit.getColumn();
+  const float time      = digit.getTime();
 
   if (labelsClusPtr) {
     int nlab = 0;
@@ -216,7 +217,7 @@ void Clusterer::ClustererThread::finishChipSingleHitFast(gsl::span<const Digit>
   cluster.xCoord = localClsCoords.x();
   cluster.yCoord = localClsCoords.y();
   cluster.zCoord = localClsCoords.z();
-
+  cluster.time = time;
   clusters.emplace_back(cluster);
 }