diff --git a/PWGLF/DataModel/Vtx3BodyTables.h b/PWGLF/DataModel/Vtx3BodyTables.h index d9070875778..e0714108b5f 100644 --- a/PWGLF/DataModel/Vtx3BodyTables.h +++ b/PWGLF/DataModel/Vtx3BodyTables.h @@ -337,11 +337,6 @@ using Vtx3BodyDatasCovsIndexed = soa::Join, NUCLEI_THREE_BODY_RECO_COLUMNS, - nuclei3body::McLabelHe3, nuclei3body::McLabelDaughter1, nuclei3body::McLabelDaughter2, nuclei3body::McPdgHe3, nuclei3body::McPdgDaughter1, nuclei3body::McPdgDaughter2, - nuclei3body::McMatchStatus, nuclei3body::McMotherLabel, nuclei3body::McMotherPdg, + nuclei3body::McMatchStatus, nuclei3body::McMotherPdg, nuclei3body::GenMotherPt, nuclei3body::GenMotherEta, nuclei3body::GenMotherPhi, nuclei3body::GenDecayLength); diff --git a/PWGLF/TableProducer/Nuspex/nucleiThreeBodyBuilder.cxx b/PWGLF/TableProducer/Nuspex/nucleiThreeBodyBuilder.cxx index a3ce4fb2d1a..154747598cc 100644 --- a/PWGLF/TableProducer/Nuspex/nucleiThreeBodyBuilder.cxx +++ b/PWGLF/TableProducer/Nuspex/nucleiThreeBodyBuilder.cxx @@ -195,8 +195,6 @@ struct NucleiThreeBodyBuilder { std::vector preselectedTracks; struct RecoCandidate { - std::array trackIds{InvalidLabel, InvalidLabel, InvalidLabel}; - int collisionId = InvalidLabel; std::array signs{}; float pt = 0.f; float eta = 0.f; @@ -218,10 +216,8 @@ struct NucleiThreeBodyBuilder { }; struct MCInfo { - std::array labels{InvalidLabel, InvalidLabel, InvalidLabel}; std::array pdgs{}; uint8_t matchStatus = MissingDaughterLabel; - int motherLabel = InvalidLabel; int motherPdg = 0; float motherPt = InvalidFloat; float motherEta = InvalidFloat; @@ -472,10 +468,6 @@ struct NucleiThreeBodyBuilder { } registry.fill(HIST("triplets"), 1.); - out.trackIds = {static_cast(he3Track.globalIndex()), - static_cast(daughter1Track.globalIndex()), - static_cast(daughter2Track.globalIndex())}; - out.collisionId = collision.globalIndex(); out.signs = {static_cast(he3Track.sign() > 0 ? 1 : -1), static_cast(daughter1Track.sign() > 0 ? 1 : -1), static_cast(daughter2Track.sign() > 0 ? 1 : -1)}; @@ -546,8 +538,7 @@ struct NucleiThreeBodyBuilder { void fillDataTable(RecoCandidate const& c) { - outputData(c.trackIds[0], c.trackIds[1], c.trackIds[2], c.collisionId, - c.signs[0], c.signs[1], c.signs[2], + outputData(c.signs[0], c.signs[1], c.signs[2], c.pt, c.eta, c.phi, c.daughterMomentum[0], c.daughterMomentum[1], c.daughterMomentum[2], c.chi2, c.dcaDaughters, c.cosPA, c.decayLength, @@ -575,9 +566,6 @@ struct NucleiThreeBodyBuilder { const auto he3Particle = he3Track.template mcParticle_as(); const auto daughter1Particle = daughter1Track.template mcParticle_as(); const auto daughter2Particle = daughter2Track.template mcParticle_as(); - info.labels = {static_cast(he3Particle.globalIndex()), - static_cast(daughter1Particle.globalIndex()), - static_cast(daughter2Particle.globalIndex())}; info.pdgs = {he3Particle.pdgCode(), daughter1Particle.pdgCode(), daughter2Particle.pdgCode()}; info.matchStatus = NoCommonImmediateMother; @@ -601,7 +589,6 @@ struct NucleiThreeBodyBuilder { const auto selectedMother = mcParticles.rawIteratorAt(commonMotherLabel); info.matchStatus = CommonImmediateMother; - info.motherLabel = commonMotherLabel; info.motherPdg = selectedMother.pdgCode(); const int he3Sign = he3Track.sign() > 0 ? 1 : -1; const bool expectedDaughters = info.pdgs[0] == he3Sign * He3Pdg && @@ -627,8 +614,7 @@ struct NucleiThreeBodyBuilder { void fillMCTable(RecoCandidate const& c, MCInfo const& m) { - outputMC(c.trackIds[0], c.trackIds[1], c.trackIds[2], c.collisionId, - c.signs[0], c.signs[1], c.signs[2], + outputMC(c.signs[0], c.signs[1], c.signs[2], c.pt, c.eta, c.phi, c.daughterMomentum[0], c.daughterMomentum[1], c.daughterMomentum[2], c.chi2, c.dcaDaughters, c.cosPA, c.decayLength, @@ -643,9 +629,8 @@ struct NucleiThreeBodyBuilder { c.tpcNCls[0], c.tpcNCls[1], c.tpcNCls[2], c.tpcCrossedRows[0], c.tpcCrossedRows[1], c.tpcCrossedRows[2], c.itsNCls[0], c.itsNCls[1], c.itsNCls[2], - m.labels[0], m.labels[1], m.labels[2], m.pdgs[0], m.pdgs[1], m.pdgs[2], - m.matchStatus, m.motherLabel, m.motherPdg, + m.matchStatus, m.motherPdg, m.motherPt, m.motherEta, m.motherPhi, m.decayLength); }