TPC timeseries: Adding pt from ITS-TPC track to tree (#12229)

* TPC timeseries: Adding pt from ITS-TPC track to tree * Adding chi2 matching and qPt
AliceO2Group · Nov 9, 2023 · d37c907 · d37c907
1 parent 694c71e
commit d37c907
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/Detectors/TPC/workflow/src/TPCTimeSeriesSpec.cxx b/Detectors/TPC/workflow/src/TPCTimeSeriesSpec.cxx
@@ -1040,6 +1040,8 @@ class TPCTimeSeries : public Task
           const auto& trkOrig = tracksTPC[iTrk];
           const bool isNearestVtx = (idxITSTPC.back() == -1); // is nearest vertex in case no vertex was found
           const float mx_ITS = hasITSTPC ? tracksITSTPC[idxITSTPC.front()].getX() : -1;
+          const float pt_ITS = hasITSTPC ? tracksITSTPC[idxITSTPC.front()].getQ2Pt() : -1;
+          const float chi2match_ITSTPC = hasITSTPC ? tracksITSTPC[idxITSTPC.front()].getChi2Match() : -1;
           const int nClITS = idxITSCheck ? tracksITS[idxITSTrack].getNClusters() : -1;
           const int chi2ITS = idxITSCheck ? tracksITS[idxITSTrack].getChi2() : -1;
           int typeSide = 2; // A- and C-Side cluster
@@ -1069,6 +1071,7 @@ class TPCTimeSeries : public Task
                               << "isNearestVertex=" << isNearestVtx
                               // tpc track properties
                               << "pt=" << trkOrig.getPt()
+                              << "qpt_ITSTPC=" << pt_ITS
                               << "tpc_timebin=" << trkOrig.getTime0()
                               << "qpt=" << trkOrig.getParam(4)
                               << "ncl=" << trkOrig.getNClusters()
@@ -1086,6 +1089,7 @@ class TPCTimeSeries : public Task
                               << "mX_ITS=" << mx_ITS
                               << "nClITS=" << nClITS
                               << "chi2ITS=" << chi2ITS
+                              << "chi2match_ITSTPC=" << chi2match_ITSTPC
                               // meta
                               << "mult=" << mNTracksWindow[iTrk]
                               << "time_window_mult=" << mTimeWindowMUS