diff --git a/PWGUD/DataModel/TauThreeProngEventTables.h b/PWGUD/DataModel/TauThreeProngEventTables.h index bd7da95f14c..e01daca155a 100644 --- a/PWGUD/DataModel/TauThreeProngEventTables.h +++ b/PWGUD/DataModel/TauThreeProngEventTables.h @@ -50,6 +50,8 @@ DECLARE_SOA_COLUMN(ZvtxFT0vsPv, zvtxFT0vsPv, int8_t); DECLARE_SOA_COLUMN(VtxITSTPC, vtxITSTPC, int8_t); DECLARE_SOA_COLUMN(ZdcAenergy, zdcAenergy, float); DECLARE_SOA_COLUMN(ZdcCenergy, zdcCenergy, float); +DECLARE_SOA_COLUMN(ZdcAtime, zdcAtime, float); +DECLARE_SOA_COLUMN(ZdcCtime, zdcCtime, float); // DECLARE_SOA_COLUMN(Qtot, qtot, int8_t); // FIT info DECLARE_SOA_COLUMN(TotalFT0AmplitudeA, totalFT0AmplitudeA, float); @@ -113,6 +115,7 @@ DECLARE_SOA_TABLE(DataTauFourTracks, "AOD", "TAUFOURTRACK", tautree::Trs, tautree::Trofs, tautree::Hmpr, tautree::Tfb, tautree::ItsRofb, tautree::Sbp, tautree::ZvtxFT0vsPv, tautree::VtxITSTPC, tautree::ZdcAenergy, tautree::ZdcCenergy, + tautree::ZdcAtime, tautree::ZdcCtime, // tautree::Qtot, tautree::TotalFT0AmplitudeA, tautree::TotalFT0AmplitudeC, tautree::TotalFV0AmplitudeA, // tautree::TimeFT0A, tautree::TimeFT0C, tautree::TimeFV0A, @@ -135,6 +138,7 @@ DECLARE_SOA_TABLE(TrueTauFourTracks, "AOD", "TRUETAU", tautree::Trs, tautree::Trofs, tautree::Hmpr, tautree::Tfb, tautree::ItsRofb, tautree::Sbp, tautree::ZvtxFT0vsPv, tautree::VtxITSTPC, tautree::ZdcAenergy, tautree::ZdcCenergy, + tautree::ZdcAtime, tautree::ZdcCtime, // tautree::Qtot, tautree::TotalFT0AmplitudeA, tautree::TotalFT0AmplitudeC, tautree::TotalFV0AmplitudeA, // tautree::TimeFT0A, tautree::TimeFT0C, tautree::TimeFV0A, diff --git a/PWGUD/TableProducer/tauThreeProngEventTableProducer.cxx b/PWGUD/TableProducer/tauThreeProngEventTableProducer.cxx index 87878b7b364..c90139faac4 100644 --- a/PWGUD/TableProducer/tauThreeProngEventTableProducer.cxx +++ b/PWGUD/TableProducer/tauThreeProngEventTableProducer.cxx @@ -63,10 +63,11 @@ using namespace o2::constants::physics; enum MyRecoProblem { NO_PROBLEM = 0, // no problem - MANY_RECO = 1, // more than 1 reconstructed collision + MANY_RECO = 1, // more than 1 reconstructed collision, event not rejected TOO_MANY_DAUGHTERS = 2, // more than 6 daughters from 2 taus TWO_TRACKS = 3, // more than 1 associated track to MC particle (tau daughter) - NO_TRACK = 4 // No associated track to MC particle (tau daughter) + NO_TRACK = 4, // No associated track to MC particle (tau daughter) + MIXED_TRACKS = 5 // 4 (or 6) tracks reconstructed but from 2 collisions }; enum MyParticle { @@ -718,6 +719,8 @@ struct TauThreeProngEventTableProducer { energyZNA = -1.; if (energyZNC < 0) energyZNC = -1.; + float timeZNA = dgcand.timeZNA(); + float timeZNC = dgcand.timeZNC(); int nTofTrk = 0; int nEtaIn15 = 0; @@ -872,6 +875,7 @@ struct TauThreeProngEventTableProducer { dgcand.trs(), dgcand.trofs(), dgcand.hmpr(), // to test it dgcand.tfb(), dgcand.itsROFb(), dgcand.sbp(), dgcand.zVtxFT0vPV(), dgcand.vtxITSTPC(), energyZNA, energyZNC, + timeZNA, timeZNC, // qtot, <<-------- comment out dgcand.totalFT0AmplitudeA(), dgcand.totalFT0AmplitudeC(), dgcand.totalFV0AmplitudeA(), // dgcand.timeFT0A(), dgcand.timeFT0C(), dgcand.timeFV0A(), @@ -1202,6 +1206,8 @@ struct TauThreeProngEventTableProducer { // zdc information - there i sno information in MC float energyZNA = -999.; float energyZNC = -999.; + float timeZNA = -999.; + float timeZNC = -999.; float amplitudesFIT[3] = {-999., -999., -999.}; // FT0A, FT0C, FV0 // float timesFIT[3] = {-999., -999., -999.}; // FT0A, FT0C, FV0 @@ -1577,6 +1583,7 @@ struct TauThreeProngEventTableProducer { bcSels[0], bcSels[1], bcSels[2], // to test it bcSels[3], bcSels[4], bcSels[5], bcSels[6], bcSels[7], energyZNA, energyZNC, + timeZNA, timeZNC, // qtot, <<-------- comment out amplitudesFIT[0], amplitudesFIT[1], amplitudesFIT[2], // timesFIT[0], timesFIT[1], timesFIT[2],