O2-4592 - adding deltaT and vertex and TPC times

Author: miranov25

Detectors/TPC/workflow/src/TPCTimeSeriesSpec.cxx

@@ -206,14 +206,14 @@ class TPCTimeSeries : public Task
       indicesITSTPC[tracksITSTPC[i].getRefTPC().getIndex()] = {i, idxVtx};
     }
 
-    std::vector<std::tuple<int, float, float,o2::track::TrackLTIntegral,float>> idxTPCTrackToTOFCluster; // store for each tpc track index the index to the TOF cluster
+    std::vector<std::tuple<int, float, float,o2::track::TrackLTIntegral,float,float>> idxTPCTrackToTOFCluster; // store for each tpc track index the index to the TOF cluster
 
     // get matches to TOF in case skimmed data is produced
     if (mUnbinnedWriter) {
       //   getLTIntegralOut(), ///< L,TOF integral calculated during the propagation
       //  getSignal()  mSignal = 0.0;              ///< TOF time in ps
       o2::track::TrackLTIntegral defLT;
-      idxTPCTrackToTOFCluster = std::vector<std::tuple<int, float, float,o2::track::TrackLTIntegral,float>>(tracksTPC.size(), {-1, -999, -999,defLT,0});
+      idxTPCTrackToTOFCluster = std::vector<std::tuple<int, float, float,o2::track::TrackLTIntegral,float,float>>(tracksTPC.size(), {-1, -999, -999,defLT,0,0});
       const std::vector<gsl::span<const o2::dataformats::MatchInfoTOF>> tofMatches{recoData.getTPCTOFMatches(), recoData.getTPCTRDTOFMatches(), recoData.getITSTPCTOFMatches(), recoData.getITSTPCTRDTOFMatches()};
 
       // loop over ITS-TPC-TRD-TOF and ITS-TPC-TOF tracks an store for each ITS-TPC track the TOF track index
@@ -223,7 +223,8 @@ class TPCTimeSeries : public Task
           if (refTPC.isIndexSet()) {
             o2::track::TrackLTIntegral ltIntegral=tpctofmatch.getLTIntegralOut();
             float signal=tpctofmatch.getSignal();
-            idxTPCTrackToTOFCluster[refTPC] = {tpctofmatch.getIdxTOFCl(), tpctofmatch.getDXatTOF(), tpctofmatch.getDZatTOF(),ltIntegral,signal};
+            float deltaT=tpctofmatch.getDeltaT();
+            idxTPCTrackToTOFCluster[refTPC] = {tpctofmatch.getIdxTOFCl(), tpctofmatch.getDXatTOF(), tpctofmatch.getDZatTOF(),ltIntegral,signal,deltaT};
           }
         }
       }
@@ -1054,7 +1055,7 @@ class TPCTimeSeries : public Task
     return isGoodTrack;
   }
 
-  void fillDCA(const gsl::span<const TrackTPC> tracksTPC, const gsl::span<const o2::dataformats::TrackTPCITS> tracksITSTPC, const gsl::span<const o2::dataformats::PrimaryVertex> vertices, const int iTrk, const int iThread, const std::unordered_map<unsigned int, std::array<int, 2>>& indicesITSTPC, const gsl::span<const o2::its::TrackITS> tracksITS, const std::vector<std::tuple<int, float, float,o2::track::TrackLTIntegral,float>>& idxTPCTrackToTOFCluster, const gsl::span<const o2::tof::Cluster> tofClusters)
+  void fillDCA(const gsl::span<const TrackTPC> tracksTPC, const gsl::span<const o2::dataformats::TrackTPCITS> tracksITSTPC, const gsl::span<const o2::dataformats::PrimaryVertex> vertices, const int iTrk, const int iThread, const std::unordered_map<unsigned int, std::array<int, 2>>& indicesITSTPC, const gsl::span<const o2::its::TrackITS> tracksITS, const std::vector<std::tuple<int, float, float,o2::track::TrackLTIntegral,float,float>>& idxTPCTrackToTOFCluster, const gsl::span<const o2::tof::Cluster> tofClusters)
   {
     const auto& trackFull = tracksTPC[iTrk];
     const bool isGoodTrack = checkTrack(trackFull);
@@ -1359,7 +1360,8 @@ class TPCTimeSeries : public Task
             }
           }
         }
-
+        float vertexTime=vertex.getTimeStamp().getTimeStamp();
+        float trackTime0=trackFull.getTime0();
         *mStreamer[iThread] << "treeTimeSeries"
                             // DCAs
                             << "triggerMask=" << triggerMask
@@ -1439,6 +1441,9 @@ class TPCTimeSeries : public Task
                             << "mDZatTOF=" << std::get<2>(idxTPCTrackToTOFCluster[iTrk])
                             << "mTOFLength=" << std::get<3>(idxTPCTrackToTOFCluster[iTrk])
                             << "mTOFSignal=" << std::get<4>(idxTPCTrackToTOFCluster[iTrk])
+                            << "mDeltaTTOFTPC=" << std::get<5>(idxTPCTrackToTOFCluster[iTrk])    /// delta T- TPC TOF
+                            << "vertexTime" << vertexTime                                        /// time stamp assigned to the vertex
+                            << "trackTime0" << trackTime0                                        /// time stamp assigned to the track
                             // TPC delta param
                             << "deltaTPCParamInOutTgl=" << deltaTPCParamInOutTgl
                             << "deltaTPCParamInOutQPt=" << deltaTPCParamInOutQPt