From c1c70b3abad78c1fdb5026a8a2db639f61848a72 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Nicol=C3=B2=20Jacazio?= Date: Sun, 19 Sep 2021 17:03:55 +0200 Subject: [PATCH] Add TrackTime, TOFChi2 and TrackTimeRes --- examples/aod/createO2tables.C | 6 ++++++ examples/aod/createO2tables.h | 16 ++++++++++------ 2 files changed, 16 insertions(+), 6 deletions(-) diff --git a/examples/aod/createO2tables.C b/examples/aod/createO2tables.C index 77f9a2b..64a5428 100644 --- a/examples/aod/createO2tables.C +++ b/examples/aod/createO2tables.C @@ -419,14 +419,20 @@ int createO2tables(const char* inputFile = "delphes.root", } aod_track.fLength = track->L * 0.1; // [cm] + aod_track.fTOFChi2 = 1.f; // Negative if TOF is not available aod_track.fTOFSignal = track->TOuter * 1.e12; // [ps] + aod_track.fTrackTime = track->TOuter * 1.e9; // [ns] + aod_track.fTrackTimeRes = 200 * 1.e9; // [ns] aod_track.fTOFExpMom = track->P * 0.029979246; // if primary push to TOF tracks if (fabs(aod_track.fY) < 3. * aod_track.fSigmaY && fabs(aod_track.fZ) < 3. * aod_track.fSigmaZ) tof_tracks.push_back(track); } else { + aod_track.fTOFChi2 = -1.f; aod_track.fLength = -999.f; aod_track.fTOFSignal = -999.f; + aod_track.fTrackTime = -999.f; + aod_track.fTrackTimeRes = 2000 * 1.e9; aod_track.fTOFExpMom = -999.f; } diff --git a/examples/aod/createO2tables.h b/examples/aod/createO2tables.h index 34c2a43..d163eb4 100644 --- a/examples/aod/createO2tables.h +++ b/examples/aod/createO2tables.h @@ -282,11 +282,13 @@ struct { Float_t fTOFChi2 = -999.f; /// chi2 TOF match (?) // PID - Float_t fTPCSignal = -999.f; /// dE/dX TPC - Float_t fTRDSignal = -999.f; /// dE/dX TRD - Float_t fTOFSignal = -999.f; /// TOFsignal - Float_t fLength = -999.f; /// Int.Lenght @ TOF - Float_t fTOFExpMom = -999.f; /// TOF Expected momentum based on the expected time of pions + Float_t fTPCSignal = -999.f; /// dE/dX TPC + Float_t fTRDSignal = -999.f; /// dE/dX TRD + Float_t fTOFSignal = -999.f; /// TOFsignal + Float_t fTrackTime = -999.f; /// fTrackTime + Float_t fTrackTimeRes = -999.f; /// fTrackTimeRes + Float_t fLength = -999.f; /// Int.Lenght @ TOF + Float_t fTOFExpMom = -999.f; /// TOF Expected momentum based on the expected time of pions // Track extrapolation to EMCAL surface Float_t fTrackEtaEMCAL = -999.f; /// Track eta at the EMCAL surface @@ -349,6 +351,8 @@ void MakeTreeO2trackExtra() tTracksExtra->Branch("fTPCSignal", &aod_track.fTPCSignal, "fTPCSignal/F"); tTracksExtra->Branch("fTRDSignal", &aod_track.fTRDSignal, "fTRDSignal/F"); tTracksExtra->Branch("fTOFSignal", &aod_track.fTOFSignal, "fTOFSignal/F"); + tTracksExtra->Branch("fTrackTime", &aod_track.fTrackTime, "fTrackTime/F"); + tTracksExtra->Branch("fTrackTimeRes", &aod_track.fTrackTimeRes, "fTrackTimeRes/F"); tTracksExtra->Branch("fLength", &aod_track.fLength, "fLength/F"); tTracksExtra->Branch("fTOFExpMom", &aod_track.fTOFExpMom, "fTOFExpMom/F"); tTracksExtra->Branch("fTrackEtaEMCAL", &aod_track.fTrackEtaEMCAL, "fTrackEtaEMCAL/F"); @@ -373,7 +377,7 @@ struct { Float_t fRICHNsigmaPi = -999.f; /// Nsigma for Pi Float_t fRICHNsigmaKa = -999.f; /// Nsigma for Ka Float_t fRICHNsigmaPr = -999.f; /// Nsigma for Pr -} rich, frich; //! structure to keep RICH info +} rich, frich; //! structure to keep RICH info void MakeTreeO2rich() {