From 1faf7892e35fac2e3104089d6253b31fa3013503 Mon Sep 17 00:00:00 2001
From: David Dobrigkeit Chinellato <david.dobrigkeit.chinellato.cern.ch>
Date: Tue, 11 Jul 2023 09:04:06 -0300
Subject: [PATCH 1/2] Propagator QA for matLUTinner tuning

---
 Common/Tasks/propagatorQa.cxx | 32 ++++++++++++++++++++++++++++++++
 1 file changed, 32 insertions(+)

diff --git a/Common/Tasks/propagatorQa.cxx b/Common/Tasks/propagatorQa.cxx
index be6cf5dee2d..29aa15160d5 100644
--- a/Common/Tasks/propagatorQa.cxx
+++ b/Common/Tasks/propagatorQa.cxx
@@ -438,6 +438,38 @@ struct propagatorQa {
     }
   }
   PROCESS_SWITCH(propagatorQa, processData, "process data", false);
+
+    void processMatLUTTest(aod::Collision const& collision, soa::Join<aod::Tracks, aod::TracksCov, aod::TracksExtra, aod::TracksDCA> const& tracks, soa::Join<aod::TracksIU, aod::TracksCovIU, aod::TracksExtra> const& tracksIU, aod::BCsWithTimestamps const&)
+  {
+    /* check the previous run number */
+    auto bc = collision.bc_as<aod::BCsWithTimestamps>();
+    initCCDB(bc);
+    gpu::gpustd::array<float, 2> dcaInfo;
+
+    for (auto& trackIU : tracksIU) {
+      if (trackIU.tpcNClsFound() < minTPCClustersRequired)
+        continue; // skip if not enough TPC clusters
+
+      if (trackIU.trackType() != aod::track::TrackIU && trackIU.x() > maxXtoConsider)
+        continue; // skip if not track IU or if beyong the max X to be considered 
+
+      o2::track::TrackParCov trackParCov = getTrackParCov(trackIU);
+
+      dcaInfo[0] = 999;
+      dcaInfo[1] = 999;
+
+      // Recalculate the propagation with this instance of the matLUT
+      o2::base::Propagator::Instance()->propagateToDCABxByBz({collision.posX(), collision.posY(), collision.posZ()}, trackParCov, maxPropagStep, matCorr, &dcaInfo);
+
+      auto track = tracks.iteratorAt(trackIU.globalIndex());
+
+      histos.fill(HIST("hDeltaTanLambdaVsPt"), track.tgl(), track.tgl() - trackParCov.getTgl());
+      histos.fill(HIST("hDeltaPtVsPt"), track.pt(), track.pt() - trackParCov.getPt());
+      histos.fill(HIST("hDeltaDCAs"), track.dcaXY() - dcaInfo[0]);
+      histos.fill(HIST("hDeltaDCAsVsPt"), track.pt(), track.dcaXY() - dcaInfo[0]);
+    }
+  }
+  PROCESS_SWITCH(propagatorQa, processMatLUTTest, "process mat lut test", false);
 };
 
 WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)

From 478f069e3bee54ea9081fb1770fafd3f487fada2 Mon Sep 17 00:00:00 2001
From: ALICE Builder <alibuild@users.noreply.github.com>
Date: Tue, 11 Jul 2023 15:11:01 +0200
Subject: [PATCH 2/2] Please consider the following formatting changes (#152)

---
 Common/Tasks/propagatorQa.cxx | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/Common/Tasks/propagatorQa.cxx b/Common/Tasks/propagatorQa.cxx
index 29aa15160d5..d2753157caa 100644
--- a/Common/Tasks/propagatorQa.cxx
+++ b/Common/Tasks/propagatorQa.cxx
@@ -439,7 +439,7 @@ struct propagatorQa {
   }
   PROCESS_SWITCH(propagatorQa, processData, "process data", false);
 
-    void processMatLUTTest(aod::Collision const& collision, soa::Join<aod::Tracks, aod::TracksCov, aod::TracksExtra, aod::TracksDCA> const& tracks, soa::Join<aod::TracksIU, aod::TracksCovIU, aod::TracksExtra> const& tracksIU, aod::BCsWithTimestamps const&)
+  void processMatLUTTest(aod::Collision const& collision, soa::Join<aod::Tracks, aod::TracksCov, aod::TracksExtra, aod::TracksDCA> const& tracks, soa::Join<aod::TracksIU, aod::TracksCovIU, aod::TracksExtra> const& tracksIU, aod::BCsWithTimestamps const&)
   {
     /* check the previous run number */
     auto bc = collision.bc_as<aod::BCsWithTimestamps>();
@@ -451,7 +451,7 @@ struct propagatorQa {
         continue; // skip if not enough TPC clusters
 
       if (trackIU.trackType() != aod::track::TrackIU && trackIU.x() > maxXtoConsider)
-        continue; // skip if not track IU or if beyong the max X to be considered 
+        continue; // skip if not track IU or if beyong the max X to be considered
 
       o2::track::TrackParCov trackParCov = getTrackParCov(trackIU);