Update

Author: njacazio

Detectors/Upgrades/ALICE3/GlobalReconstruction/reconstruction/CMakeLists.txt

@@ -10,7 +10,10 @@
 # or submit itself to any jurisdiction.
 
 if(Acts_FOUND)
-  set(actsTarget Acts::Core)
+  set(actsTarget
+        Acts::Core
+        O2::ACTSInterface
+        O2::ALICE3ACTS)
 endif()
 
 o2_add_library(ALICE3GlobalReconstruction

Detectors/Upgrades/ALICE3/GlobalReconstruction/reconstruction/include/ALICE3GlobalReconstruction/TrackerACTS.h

@@ -19,11 +19,11 @@
 #ifndef ALICE3_GLOBALRECONSTRUCTION_INCLUDE_TRACKERACTS_H_
 #define ALICE3_GLOBALRECONSTRUCTION_INCLUDE_TRACKERACTS_H_
 
-#include "Acts/Definitions/Units.hpp"
-#include "Framework/Logger.h"
-
-#include "ITStracking/TimeFrame.h"
-#include "TH2F.h"
+#include <Acts/Definitions/Units.hpp>
+#include <Framework/Logger.h>
+#include <DataFormatsITSMFT/CompCluster.h>
+#include <ITStracking/TimeFrame.h>
+#include <TH2F.h>
 
 namespace o2::trk
 {
@@ -62,6 +62,8 @@ struct SpacePoint {
   int layer{-1};
   int clusterId{-1};
   int rof{-1};
+  int sensorId{-1};
+  // const Acts::Surface* surface{nullptr};
 
   // Derived quantities
   float r() const { return std::hypot(x, y); }
@@ -99,6 +101,11 @@ class TrackerACTS
   /// Adopt a TimeFrame for processing
   void adoptTimeFrame(o2::its::TimeFrame<nLayers>& tf);
 
+  void adoptClusters(gsl::span<const o2::itsmft::CompClusterExt> clusters)
+  {
+    mInputClusters = clusters;
+  }
+
   /// Main tracking entry point: convert clusters to tracks
   void clustersToTracks();
 
@@ -155,6 +162,7 @@ class TrackerACTS
   float mBz{0.5f}; ///< Magnetic field in Tesla
   unsigned int mTimeFrameCounter{0};
   double mTotalTime{0.};
+  gsl::span<const o2::itsmft::CompClusterExt> mInputClusters;
 
   // Tracking states for logging
   enum State {

Detectors/Upgrades/ALICE3/GlobalReconstruction/reconstruction/src/TrackerACTS.cxx

@@ -18,6 +18,10 @@
 
 #include "ALICE3GlobalReconstruction/TrackerACTS.h"
 
+#include <TRKBase/GeometryTGeo.h>
+#include <MathUtils/Utils.h>
+#include <Acts/Surfaces/Surface.hpp>
+#include "ACTSInterface/TrackingGeometryManager.h"
 #include <Acts/EventData/Seed.hpp>
 #include <Acts/EventData/SpacePointContainer.hpp>
 #include <Acts/Seeding/BinnedGroup.hpp>
@@ -273,7 +277,7 @@ bool TrackerACTS<nLayers>::estimateTrackParams(const SeedACTS& seed, o2::its::Tr
   const float x = p[Acts::eFreePos0];
   const float y = p[Acts::eFreePos1];
   const float z = p[Acts::eFreePos2];
-  const float px = p[Acts::eFreeDir0] * pMag;
+  const float px = p[Acts::eFreeDir0] * 1;
   const float py = p[Acts::eFreeDir1] * pMag;
   const float pz = p[Acts::eFreeDir2] * pMag;
   const int charge = (p[Acts::eFreeQOverP] >= 0) ? 1 : -1;
@@ -395,6 +399,19 @@ void TrackerACTS<nLayers>::clustersToTracks()
   //                             StateNames[mCurState]);
   // }
 
+  auto& actsGeometryManager = o2::acts::TrackingGeometryManager::instance();
+
+  const auto& trackingGeometry = actsGeometryManager.get();
+  const auto& geometryContext = actsGeometryManager.getNominalContext();
+
+  auto* trkGeometry = o2::trk::GeometryTGeo::Instance();
+
+  trkGeometry->fillMatrixCache(o2::math_utils::bit2Mask(o2::math_utils::TransformType::L2G));
+
+  const auto& surfaceIndex = actsGeometryManager.getIndex(*trkGeometry, 1.e-3, [trkGeometry](int sensorID) {
+    return std::string(trkGeometry->getMatrixPath(sensorID).Data());
+  });
+
   // MC labeling
   if (mTimeFrame->hasMCinformation()) {
     computeTracksMClabels();