From 15b1573f210cfdb63c64c9b3eca65a1d3b5f8117 Mon Sep 17 00:00:00 2001 From: ZFederica Date: Fri, 16 Aug 2024 17:19:45 +0200 Subject: [PATCH 01/10] Change track table for LF daughters --- .../candidateCreatorXic0Omegac0.cxx | 117 +++++++++++------- .../candidateSelectorOmegac0ToOmegaKa.cxx | 10 +- .../candidateSelectorOmegac0ToOmegaPi.cxx | 10 +- .../TableProducer/candidateSelectorToXiPi.cxx | 10 +- 4 files changed, 93 insertions(+), 54 deletions(-) diff --git a/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx b/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx index d68fbdb7486..229ff962366 100644 --- a/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx +++ b/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx @@ -115,6 +115,7 @@ struct HfCandidateCreatorXic0Omegac0 { using MyCascTable = soa::Join; // to use strangeness tracking, use aod::TraCascDatas instead of aod::CascDatas using CascadesLinked = soa::Join; using MyV0Table = soa::Join; + using MyLFTracksWCov = soa::Join; using MyKfTracks = soa::Join; using MyKfCascTable = soa::Join; @@ -256,6 +257,7 @@ struct HfCandidateCreatorXic0Omegac0 { void runXic0Omegac0Creator(Coll const&, aod::BCsWithTimestamps const& /*bcWithTimeStamps*/, TracksWCovDca const&, + MyLFTracksWCov const&, MyCascTable const&, CascadesLinked const&, aod::HfCascLf2Prongs const& candidates, Hist& hInvMassCharmBaryon, @@ -317,9 +319,9 @@ struct HfCandidateCreatorXic0Omegac0 { } auto casc = cascAodElement.cascData_as(); hCascadesCounter->Fill(1); - auto trackCascDauCharged = casc.bachelor_as(); // pion <- xi track - auto trackV0Dau0 = casc.posTrack_as(); // V0 positive daughter track - auto trackV0Dau1 = casc.negTrack_as(); // V0 negative daughter track + auto trackCascDauCharged = casc.bachelor_as(); // pion <- xi track + auto trackV0Dau0 = casc.posTrack_as(); // V0 positive daughter track + auto trackV0Dau1 = casc.negTrack_as(); // V0 negative daughter track //-------------------------- V0 info--------------------------- // pseudorapidity @@ -398,20 +400,27 @@ struct HfCandidateCreatorXic0Omegac0 { // pseudorapidity float pseudorapCharmBachelor = trackCharmBachelor.eta(); - // DCAxy (computed with propagateToDCABxByBz method) - float dcaxyV0Dau0 = trackV0Dau0.dcaXY(); - float dcaxyV0Dau1 = trackV0Dau1.dcaXY(); - float dcaxyCascBachelor = trackCascDauCharged.dcaXY(); - - // DCAz (computed with propagateToDCABxByBz method) - float dcazV0Dau0 = trackV0Dau0.dcaZ(); - float dcazV0Dau1 = trackV0Dau1.dcaZ(); - float dcazCascBachelor = trackCascDauCharged.dcaZ(); - // primary vertex of the collision auto primaryVertex = getPrimaryVertex(collision); // get the associated covariance matrix with auto covMatrixPV = primaryVertex.getCov(); std::array pvCoord = {collision.posX(), collision.posY(), collision.posZ()}; + // DCAxy and DCAz (computed with propagateToDCABxByBz method) + o2::dataformats::DCA impactParameterV0Dau0; + o2::dataformats::DCA impactParameterV0Dau1; + o2::dataformats::DCA impactParameterCascDauCharged; + auto trackParVarV0Dau0 = getTrackParCov(trackV0Dau0); + auto trackParVarV0Dau1 = getTrackParCov(trackV0Dau1); + auto trackParVarCascDauCharged = getTrackParCov(trackCascDauCharged); + o2::base::Propagator::Instance()->propagateToDCABxByBz(primaryVertex, trackParVarV0Dau0, 2.f, matCorr, &impactParameterV0Dau0); + o2::base::Propagator::Instance()->propagateToDCABxByBz(primaryVertex, trackParVarV0Dau1, 2.f, matCorr, &impactParameterV0Dau1); + o2::base::Propagator::Instance()->propagateToDCABxByBz(primaryVertex, trackParVarCascDauCharged, 2.f, matCorr, &impactParameterCascDauCharged); + float dcaxyV0Dau0 = impactParameterV0Dau0.getY(); + float dcaxyV0Dau1 = impactParameterV0Dau1.getY(); + float dcaxyCascBachelor = impactParameterCascDauCharged.getY(); + float dcazV0Dau0 = impactParameterV0Dau0.dcaZ(); + float dcazV0Dau1 = impactParameterV0Dau1.dcaZ(); + float dcazCascBachelor = impactParameterCascDauCharged.dcaZ(); + // impact parameters o2::dataformats::DCA impactParameterCasc; o2::dataformats::DCA impactParameterCharmBachelor; @@ -991,28 +1000,31 @@ struct HfCandidateCreatorXic0Omegac0 { void processNoCentToXiPi(soa::Join const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, TracksWCovDca const& tracks, + MyLFTracksWCov const &lfTracks, MyCascTable const& cascades, CascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) { - runXic0Omegac0Creator(collisions, bcWithTimeStamps, tracks, cascades, cascadeLinks, candidates, hInvMassCharmBaryonToXiPi, hFitterStatusToXiPi, hCandidateCounterToXiPi, hCascadesCounterToXiPi); + runXic0Omegac0Creator(collisions, bcWithTimeStamps, tracks, lfTracks, cascades, cascadeLinks, candidates, hInvMassCharmBaryonToXiPi, hFitterStatusToXiPi, hCandidateCounterToXiPi, hCascadesCounterToXiPi); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0, processNoCentToXiPi, "Run candidate creator w/o centrality selections for xi pi decay channel", true); void processNoCentToOmegaPi(soa::Join const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, TracksWCovDca const& tracks, + MyLFTracksWCov const &lfTracks, MyCascTable const& cascades, CascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) { - runXic0Omegac0Creator(collisions, bcWithTimeStamps, tracks, cascades, cascadeLinks, candidates, hInvMassCharmBaryonToOmegaPi, hFitterStatusToOmegaPi, hCandidateCounterToOmegaPi, hCascadesCounterToOmegaPi); + runXic0Omegac0Creator(collisions, bcWithTimeStamps, tracks, lfTracks, cascades, cascadeLinks, candidates, hInvMassCharmBaryonToOmegaPi, hFitterStatusToOmegaPi, hCandidateCounterToOmegaPi, hCascadesCounterToOmegaPi); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0, processNoCentToOmegaPi, "Run candidate creator w/o centrality selections for omega pi decay channel", false); void processOmegacToOmegaPiWithKFParticle(aod::Collisions const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, MyKfTracks const& tracks, + MyLFTracksWCov const &lfTracks, MyKfCascTable const& cascades, KFCascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) @@ -1024,11 +1036,12 @@ struct HfCandidateCreatorXic0Omegac0 { void processNoCentToOmegaK(soa::Join const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, TracksWCovDca const& tracks, + MyLFTracksWCov const &lfTracks, MyCascTable const& cascades, CascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) { - runXic0Omegac0Creator(collisions, bcWithTimeStamps, tracks, cascades, cascadeLinks, candidates, hInvMassCharmBaryonToOmegaK, hFitterStatusToOmegaK, hCandidateCounterToOmegaK, hCascadesCounterToOmegaK); + runXic0Omegac0Creator(collisions, bcWithTimeStamps, tracks, lfTracks, cascades, cascadeLinks, candidates, hInvMassCharmBaryonToOmegaK, hFitterStatusToOmegaK, hCandidateCounterToOmegaK, hCascadesCounterToOmegaK); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0, processNoCentToOmegaK, "Run candidate creator w/o centrality selections for omega K decay channel", false); @@ -1036,33 +1049,36 @@ struct HfCandidateCreatorXic0Omegac0 { void processCentFT0CToXiPi(soa::Join const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, TracksWCovDca const& tracks, + MyLFTracksWCov const &lfTracks, MyCascTable const& cascades, CascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) { - runXic0Omegac0Creator(collisions, bcWithTimeStamps, tracks, cascades, cascadeLinks, candidates, hInvMassCharmBaryonToXiPi, hFitterStatusToXiPi, hCandidateCounterToXiPi, hCascadesCounterToXiPi); + runXic0Omegac0Creator(collisions, bcWithTimeStamps, tracks, lfTracks, cascades, cascadeLinks, candidates, hInvMassCharmBaryonToXiPi, hFitterStatusToXiPi, hCandidateCounterToXiPi, hCascadesCounterToXiPi); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0, processCentFT0CToXiPi, "Run candidate creator w/ centrality selection on FT0C for xi pi channel", false); void processCentFT0CToOmegaPi(soa::Join const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, TracksWCovDca const& tracks, + MyLFTracksWCov const &lfTracks, MyCascTable const& cascades, CascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) { - runXic0Omegac0Creator(collisions, bcWithTimeStamps, tracks, cascades, cascadeLinks, candidates, hInvMassCharmBaryonToOmegaPi, hFitterStatusToOmegaPi, hCandidateCounterToOmegaPi, hCascadesCounterToOmegaPi); + runXic0Omegac0Creator(collisions, bcWithTimeStamps, tracks, lfTracks, cascades, cascadeLinks, candidates, hInvMassCharmBaryonToOmegaPi, hFitterStatusToOmegaPi, hCandidateCounterToOmegaPi, hCascadesCounterToOmegaPi); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0, processCentFT0CToOmegaPi, "Run candidate creator w/ centrality selection on FT0C for omega pi channel", false); void processCentFT0CToOmegaK(soa::Join const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, TracksWCovDca const& tracks, + MyLFTracksWCov const &lfTracks, MyCascTable const& cascades, CascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) { - runXic0Omegac0Creator(collisions, bcWithTimeStamps, tracks, cascades, cascadeLinks, candidates, hInvMassCharmBaryonToOmegaK, hFitterStatusToOmegaK, hCandidateCounterToOmegaK, hCascadesCounterToOmegaK); + runXic0Omegac0Creator(collisions, bcWithTimeStamps, tracks, lfTracks, cascades, cascadeLinks, candidates, hInvMassCharmBaryonToOmegaK, hFitterStatusToOmegaK, hCandidateCounterToOmegaK, hCascadesCounterToOmegaK); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0, processCentFT0CToOmegaK, "Run candidate creator w/ centrality selection on FT0C for omega K channel", false); @@ -1070,33 +1086,36 @@ struct HfCandidateCreatorXic0Omegac0 { void processCentFT0MToXiPi(soa::Join const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, TracksWCovDca const& tracks, + MyLFTracksWCov const &lfTracks, MyCascTable const& cascades, CascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) { - runXic0Omegac0Creator(collisions, bcWithTimeStamps, tracks, cascades, cascadeLinks, candidates, hInvMassCharmBaryonToXiPi, hFitterStatusToXiPi, hCandidateCounterToXiPi, hCascadesCounterToXiPi); + runXic0Omegac0Creator(collisions, bcWithTimeStamps, tracks, lfTracks, cascades, cascadeLinks, candidates, hInvMassCharmBaryonToXiPi, hFitterStatusToXiPi, hCandidateCounterToXiPi, hCascadesCounterToXiPi); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0, processCentFT0MToXiPi, "Run candidate creator w/ centrality selection on FT0M for xi pi channel", false); void processCentFT0MToOmegaPi(soa::Join const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, TracksWCovDca const& tracks, + MyLFTracksWCov const &lfTracks, MyCascTable const& cascades, CascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) { - runXic0Omegac0Creator(collisions, bcWithTimeStamps, tracks, cascades, cascadeLinks, candidates, hInvMassCharmBaryonToOmegaPi, hFitterStatusToOmegaPi, hCandidateCounterToOmegaPi, hCascadesCounterToOmegaPi); + runXic0Omegac0Creator(collisions, bcWithTimeStamps, tracks, lfTracks, cascades, cascadeLinks, candidates, hInvMassCharmBaryonToOmegaPi, hFitterStatusToOmegaPi, hCandidateCounterToOmegaPi, hCascadesCounterToOmegaPi); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0, processCentFT0MToOmegaPi, "Run candidate creator w/ centrality selection on FT0M for omega pi channel", false); void processCentFT0MToOmegaK(soa::Join const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, TracksWCovDca const& tracks, + MyLFTracksWCov const &lfTracks, MyCascTable const& cascades, CascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) { - runXic0Omegac0Creator(collisions, bcWithTimeStamps, tracks, cascades, cascadeLinks, candidates, hInvMassCharmBaryonToOmegaK, hFitterStatusToOmegaK, hCandidateCounterToOmegaK, hCascadesCounterToOmegaK); + runXic0Omegac0Creator(collisions, bcWithTimeStamps, tracks, lfTracks, cascades, cascadeLinks, candidates, hInvMassCharmBaryonToOmegaK, hFitterStatusToOmegaK, hCandidateCounterToOmegaK, hCascadesCounterToOmegaK); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0, processCentFT0MToOmegaK, "Run candidate creator w/ centrality selection on FT0M for omega K channel", false); @@ -1173,6 +1192,7 @@ struct HfCandidateCreatorXic0Omegac0Mc { // Configuration o2::framework::Configurable rejectBackground{"rejectBackground", true, "Reject particles from background events"}; + using lfTracksWMc = soa::Join; using McCollisionsNoCents = soa::Join; using McCollisionsFT0Cs = soa::Join; using McCollisionsFT0Ms = soa::Join; @@ -1233,6 +1253,7 @@ struct HfCandidateCreatorXic0Omegac0Mc { template void runXic0Omegac0Mc(TMyRecoCand const& candidates, aod::TracksWMc const&, + lfTracksWMc const&, aod::McParticles const& mcParticles, Colls const& collsWithMcLabels, aod::McCollisions const& mcCollisions, @@ -1273,14 +1294,14 @@ struct HfCandidateCreatorXic0Omegac0Mc { std::vector idxBhadMothers{}; auto arrayDaughters = std::array{candidate.template bachelorFromCharmBaryon_as(), // bachelor <- charm baryon - candidate.template bachelor_as(), // bachelor <- cascade - candidate.template posTrack_as(), // p <- lambda - candidate.template negTrack_as()}; // pi <- lambda - auto arrayDaughtersCasc = std::array{candidate.template bachelor_as(), - candidate.template posTrack_as(), - candidate.template negTrack_as()}; - auto arrayDaughtersV0 = std::array{candidate.template posTrack_as(), - candidate.template negTrack_as()}; + candidate.template bachelor_as(), // bachelor <- cascade + candidate.template posTrack_as(), // p <- lambda + candidate.template negTrack_as()}; // pi <- lambda + auto arrayDaughtersCasc = std::array{candidate.template bachelor_as(), + candidate.template posTrack_as(), + candidate.template negTrack_as()}; + auto arrayDaughtersV0 = std::array{candidate.template posTrack_as(), + candidate.template negTrack_as()}; // Check whether the particle is from background events. If so, reject it. if (rejectBackground) { @@ -1703,133 +1724,145 @@ struct HfCandidateCreatorXic0Omegac0Mc { void processMcXicToXiPi(aod::HfCandToXiPi const& candidates, aod::TracksWMc const& tracks, + lfTracksWMc const& lfTracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsNoCents const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcXicToXiPi, "Run Xic0 to xi pi MC process function - no centrality", false); void processMcXicToXiPiFT0m(aod::HfCandToXiPi const& candidates, aod::TracksWMc const& tracks, + lfTracksWMc const& lfTracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsFT0Ms const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcXicToXiPiFT0m, "Run Xic0 to xi pi MC process function - FT0M", false); void processMcXicToXiPiFT0c(aod::HfCandToXiPi const& candidates, aod::TracksWMc const& tracks, + lfTracksWMc const& lfTracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsFT0Cs const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcXicToXiPiFT0c, "Run Xic0 to xi pi MC process function - FT0C", false); void processMcOmegacToXiPi(aod::HfCandToXiPi const& candidates, aod::TracksWMc const& tracks, + lfTracksWMc const& lfTracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsNoCents const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcOmegacToXiPi, "Run Omegac0 to xi pi MC process function - no centrality", false); void processMcOmegacToXiPiFT0m(aod::HfCandToXiPi const& candidates, aod::TracksWMc const& tracks, + lfTracksWMc const& lfTracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsFT0Ms const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcOmegacToXiPiFT0m, "Run Omegac0 to xi pi MC process function - FT0M", false); void processMcOmegacToXiPiFT0c(aod::HfCandToXiPi const& candidates, aod::TracksWMc const& tracks, + lfTracksWMc const& lfTracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsFT0Cs const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcOmegacToXiPiFT0c, "Run Omegac0 to xi pi MC process function - FT0C", false); void processMcOmegacToOmegaPi(aod::HfCandToOmegaPi const& candidates, aod::TracksWMc const& tracks, + lfTracksWMc const& lfTracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsNoCents const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcOmegacToOmegaPi, "Run Omegac0 to omega pi MC process function - no centrality", false); void processMcOmegacToOmegaPiFT0m(aod::HfCandToOmegaPi const& candidates, aod::TracksWMc const& tracks, + lfTracksWMc const& lfTracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsFT0Ms const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcOmegacToOmegaPiFT0m, "Run Omegac0 to omega pi MC process function - FT0M", false); void processMcOmegacToOmegaPiFT0c(aod::HfCandToOmegaPi const& candidates, aod::TracksWMc const& tracks, + lfTracksWMc const& lfTracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsFT0Cs const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcOmegacToOmegaPiFT0c, "Run Omegac0 to omega pi MC process function - FT0C", false); void processMcOmegacToOmegaK(aod::HfCandToOmegaK const& candidates, aod::TracksWMc const& tracks, + lfTracksWMc const& lfTracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsNoCents const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcOmegacToOmegaK, "Run Omegac0 to omega K MC process function - no centrality", false); void processMcOmegacToOmegaKFT0m(aod::HfCandToOmegaK const& candidates, aod::TracksWMc const& tracks, + lfTracksWMc const& lfTracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsFT0Ms const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcOmegacToOmegaKFT0m, "Run Omegac0 to omega K MC process function - FT0M", false); void processMcOmegacToOmegaKFT0c(aod::HfCandToOmegaK const& candidates, aod::TracksWMc const& tracks, + lfTracksWMc const& lfTracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsFT0Cs const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcOmegacToOmegaKFT0c, "Run Omegac0 to omega K MC process function - FT0C", false); diff --git a/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaKa.cxx b/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaKa.cxx index 12f2dddc6f6..144e51108ee 100644 --- a/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaKa.cxx +++ b/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaKa.cxx @@ -129,6 +129,7 @@ struct HfCandidateSelectorToOmegaKa { TrackSelectorKa selectorKaon; using TracksSel = soa::Join; + using TracksSelLf = soa::Join; HistogramRegistry registry{"registry"}; // for QA of selections @@ -194,7 +195,8 @@ struct HfCandidateSelectorToOmegaKa { } void process(aod::HfCandToOmegaK const& candidates, - TracksSel const&) + TracksSel const&, + TracksSelLf const&) { // looping over charm baryon candidates @@ -202,9 +204,9 @@ struct HfCandidateSelectorToOmegaKa { bool resultSelections = true; // True if the candidate passes all the selections, False otherwise - auto trackV0PosDau = candidate.posTrack_as(); // positive V0 daughter - auto trackV0NegDau = candidate.negTrack_as(); // negative V0 daughter - auto trackKaFromCasc = candidate.bachelor_as(); // kaon <- cascade + auto trackV0PosDau = candidate.posTrack_as(); // positive V0 daughter + auto trackV0NegDau = candidate.negTrack_as(); // negative V0 daughter + auto trackKaFromCasc = candidate.bachelor_as(); // kaon <- cascade auto trackKaFromCharm = candidate.bachelorFromCharmBaryon_as(); // kaon <- charm baryon auto trackPiFromLam = trackV0NegDau; diff --git a/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx b/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx index 2ecc7afc038..65d12548f5c 100644 --- a/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx +++ b/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx @@ -129,6 +129,7 @@ struct HfCandidateSelectorToOmegaPi { TrackSelectorKa selectorKaon; using TracksSel = soa::Join; + using TracksSelLf = soa::Join; HistogramRegistry registry{"registry"}; // for QA of selections @@ -194,16 +195,17 @@ struct HfCandidateSelectorToOmegaPi { } void process(aod::HfCandToOmegaPi const& candidates, - TracksSel const&) + TracksSel const&, + TracksSelLf const&) { // looping over charm baryon candidates for (const auto& candidate : candidates) { bool resultSelections = true; // True if the candidate passes all the selections, False otherwise - auto trackV0PosDau = candidate.posTrack_as(); // positive V0 daughter - auto trackV0NegDau = candidate.negTrack_as(); // negative V0 daughter - auto trackKaFromCasc = candidate.bachelor_as(); // kaon <- cascade + auto trackV0PosDau = candidate.posTrack_as(); // positive V0 daughter + auto trackV0NegDau = candidate.negTrack_as(); // negative V0 daughter + auto trackKaFromCasc = candidate.bachelor_as(); // kaon <- cascade auto trackPiFromCharm = candidate.bachelorFromCharmBaryon_as(); // pion <- charm baryon auto trackPiFromLam = trackV0NegDau; diff --git a/PWGHF/TableProducer/candidateSelectorToXiPi.cxx b/PWGHF/TableProducer/candidateSelectorToXiPi.cxx index eff8377e44a..86b15b24ee8 100644 --- a/PWGHF/TableProducer/candidateSelectorToXiPi.cxx +++ b/PWGHF/TableProducer/candidateSelectorToXiPi.cxx @@ -118,6 +118,7 @@ struct HfCandidateSelectorToXiPi { TrackSelectorPr selectorProton; using TracksSel = soa::Join; + using TracksSelLf = soa::Join; HistogramRegistry registry{"registry"}; // for QA of selections @@ -176,7 +177,8 @@ struct HfCandidateSelectorToXiPi { } void process(aod::HfCandToXiPi const& candidates, - TracksSel const&) + TracksSel const&, + TracksSelLf const&) { double massLambdaFromPDG = o2::constants::physics::MassLambda0; @@ -187,9 +189,9 @@ struct HfCandidateSelectorToXiPi { bool resultSelections = true; // True if the candidate passes all the selections, False otherwise - auto trackV0PosDau = candidate.posTrack_as(); // positive V0 daughter - auto trackV0NegDau = candidate.negTrack_as(); // negative V0 daughter - auto trackPiFromCasc = candidate.bachelor_as(); // pion <- cascade + auto trackV0PosDau = candidate.posTrack_as(); // positive V0 daughter + auto trackV0NegDau = candidate.negTrack_as(); // negative V0 daughter + auto trackPiFromCasc = candidate.bachelor_as(); // pion <- cascade auto trackPiFromCharm = candidate.bachelorFromCharmBaryon_as(); // pion <- charm baryon auto trackPiFromLam = trackV0NegDau; From edaf4454b628ad8e08311c4aa2c7767fe9acb752 Mon Sep 17 00:00:00 2001 From: ZFederica Date: Fri, 16 Aug 2024 17:22:12 +0200 Subject: [PATCH 02/10] Remove unnecessary change --- PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx | 1 - 1 file changed, 1 deletion(-) diff --git a/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx b/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx index 229ff962366..3ccba2d70a9 100644 --- a/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx +++ b/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx @@ -1024,7 +1024,6 @@ struct HfCandidateCreatorXic0Omegac0 { void processOmegacToOmegaPiWithKFParticle(aod::Collisions const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, MyKfTracks const& tracks, - MyLFTracksWCov const &lfTracks, MyKfCascTable const& cascades, KFCascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) From a20477203b712a57ca5ac4a9e294ad6aaf298e63 Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Fri, 16 Aug 2024 15:26:07 +0000 Subject: [PATCH 03/10] Please consider the following formatting changes --- .../candidateCreatorXic0Omegac0.cxx | 24 +++++++++---------- .../candidateSelectorOmegac0ToOmegaKa.cxx | 6 ++--- .../candidateSelectorOmegac0ToOmegaPi.cxx | 6 ++--- .../TableProducer/candidateSelectorToXiPi.cxx | 6 ++--- 4 files changed, 21 insertions(+), 21 deletions(-) diff --git a/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx b/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx index 3ccba2d70a9..8e48a7faeed 100644 --- a/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx +++ b/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx @@ -1000,7 +1000,7 @@ struct HfCandidateCreatorXic0Omegac0 { void processNoCentToXiPi(soa::Join const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, TracksWCovDca const& tracks, - MyLFTracksWCov const &lfTracks, + MyLFTracksWCov const& lfTracks, MyCascTable const& cascades, CascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) @@ -1012,7 +1012,7 @@ struct HfCandidateCreatorXic0Omegac0 { void processNoCentToOmegaPi(soa::Join const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, TracksWCovDca const& tracks, - MyLFTracksWCov const &lfTracks, + MyLFTracksWCov const& lfTracks, MyCascTable const& cascades, CascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) @@ -1035,7 +1035,7 @@ struct HfCandidateCreatorXic0Omegac0 { void processNoCentToOmegaK(soa::Join const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, TracksWCovDca const& tracks, - MyLFTracksWCov const &lfTracks, + MyLFTracksWCov const& lfTracks, MyCascTable const& cascades, CascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) @@ -1048,7 +1048,7 @@ struct HfCandidateCreatorXic0Omegac0 { void processCentFT0CToXiPi(soa::Join const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, TracksWCovDca const& tracks, - MyLFTracksWCov const &lfTracks, + MyLFTracksWCov const& lfTracks, MyCascTable const& cascades, CascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) @@ -1060,7 +1060,7 @@ struct HfCandidateCreatorXic0Omegac0 { void processCentFT0CToOmegaPi(soa::Join const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, TracksWCovDca const& tracks, - MyLFTracksWCov const &lfTracks, + MyLFTracksWCov const& lfTracks, MyCascTable const& cascades, CascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) @@ -1072,7 +1072,7 @@ struct HfCandidateCreatorXic0Omegac0 { void processCentFT0CToOmegaK(soa::Join const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, TracksWCovDca const& tracks, - MyLFTracksWCov const &lfTracks, + MyLFTracksWCov const& lfTracks, MyCascTable const& cascades, CascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) @@ -1085,7 +1085,7 @@ struct HfCandidateCreatorXic0Omegac0 { void processCentFT0MToXiPi(soa::Join const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, TracksWCovDca const& tracks, - MyLFTracksWCov const &lfTracks, + MyLFTracksWCov const& lfTracks, MyCascTable const& cascades, CascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) @@ -1097,7 +1097,7 @@ struct HfCandidateCreatorXic0Omegac0 { void processCentFT0MToOmegaPi(soa::Join const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, TracksWCovDca const& tracks, - MyLFTracksWCov const &lfTracks, + MyLFTracksWCov const& lfTracks, MyCascTable const& cascades, CascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) @@ -1109,7 +1109,7 @@ struct HfCandidateCreatorXic0Omegac0 { void processCentFT0MToOmegaK(soa::Join const& collisions, aod::BCsWithTimestamps const& bcWithTimeStamps, TracksWCovDca const& tracks, - MyLFTracksWCov const &lfTracks, + MyLFTracksWCov const& lfTracks, MyCascTable const& cascades, CascadesLinked const& cascadeLinks, aod::HfCascLf2Prongs const& candidates) @@ -1293,9 +1293,9 @@ struct HfCandidateCreatorXic0Omegac0Mc { std::vector idxBhadMothers{}; auto arrayDaughters = std::array{candidate.template bachelorFromCharmBaryon_as(), // bachelor <- charm baryon - candidate.template bachelor_as(), // bachelor <- cascade - candidate.template posTrack_as(), // p <- lambda - candidate.template negTrack_as()}; // pi <- lambda + candidate.template bachelor_as(), // bachelor <- cascade + candidate.template posTrack_as(), // p <- lambda + candidate.template negTrack_as()}; // pi <- lambda auto arrayDaughtersCasc = std::array{candidate.template bachelor_as(), candidate.template posTrack_as(), candidate.template negTrack_as()}; diff --git a/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaKa.cxx b/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaKa.cxx index 144e51108ee..aba94cb1add 100644 --- a/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaKa.cxx +++ b/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaKa.cxx @@ -204,9 +204,9 @@ struct HfCandidateSelectorToOmegaKa { bool resultSelections = true; // True if the candidate passes all the selections, False otherwise - auto trackV0PosDau = candidate.posTrack_as(); // positive V0 daughter - auto trackV0NegDau = candidate.negTrack_as(); // negative V0 daughter - auto trackKaFromCasc = candidate.bachelor_as(); // kaon <- cascade + auto trackV0PosDau = candidate.posTrack_as(); // positive V0 daughter + auto trackV0NegDau = candidate.negTrack_as(); // negative V0 daughter + auto trackKaFromCasc = candidate.bachelor_as(); // kaon <- cascade auto trackKaFromCharm = candidate.bachelorFromCharmBaryon_as(); // kaon <- charm baryon auto trackPiFromLam = trackV0NegDau; diff --git a/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx b/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx index 65d12548f5c..aabb13137e8 100644 --- a/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx +++ b/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx @@ -203,9 +203,9 @@ struct HfCandidateSelectorToOmegaPi { bool resultSelections = true; // True if the candidate passes all the selections, False otherwise - auto trackV0PosDau = candidate.posTrack_as(); // positive V0 daughter - auto trackV0NegDau = candidate.negTrack_as(); // negative V0 daughter - auto trackKaFromCasc = candidate.bachelor_as(); // kaon <- cascade + auto trackV0PosDau = candidate.posTrack_as(); // positive V0 daughter + auto trackV0NegDau = candidate.negTrack_as(); // negative V0 daughter + auto trackKaFromCasc = candidate.bachelor_as(); // kaon <- cascade auto trackPiFromCharm = candidate.bachelorFromCharmBaryon_as(); // pion <- charm baryon auto trackPiFromLam = trackV0NegDau; diff --git a/PWGHF/TableProducer/candidateSelectorToXiPi.cxx b/PWGHF/TableProducer/candidateSelectorToXiPi.cxx index 86b15b24ee8..035fd7ed9dc 100644 --- a/PWGHF/TableProducer/candidateSelectorToXiPi.cxx +++ b/PWGHF/TableProducer/candidateSelectorToXiPi.cxx @@ -189,9 +189,9 @@ struct HfCandidateSelectorToXiPi { bool resultSelections = true; // True if the candidate passes all the selections, False otherwise - auto trackV0PosDau = candidate.posTrack_as(); // positive V0 daughter - auto trackV0NegDau = candidate.negTrack_as(); // negative V0 daughter - auto trackPiFromCasc = candidate.bachelor_as(); // pion <- cascade + auto trackV0PosDau = candidate.posTrack_as(); // positive V0 daughter + auto trackV0NegDau = candidate.negTrack_as(); // negative V0 daughter + auto trackPiFromCasc = candidate.bachelor_as(); // pion <- cascade auto trackPiFromCharm = candidate.bachelorFromCharmBaryon_as(); // pion <- charm baryon auto trackPiFromLam = trackV0NegDau; From 1c792e9c15883ba893a03c79d0acf8406dde9dcb Mon Sep 17 00:00:00 2001 From: ZFederica Date: Fri, 16 Aug 2024 23:07:03 +0200 Subject: [PATCH 04/10] Fix build error --- PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaKa.cxx | 2 +- PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaKa.cxx b/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaKa.cxx index aba94cb1add..6f232004033 100644 --- a/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaKa.cxx +++ b/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaKa.cxx @@ -129,7 +129,7 @@ struct HfCandidateSelectorToOmegaKa { TrackSelectorKa selectorKaon; using TracksSel = soa::Join; - using TracksSelLf = soa::Join; + using TracksSelLf = soa::Join; HistogramRegistry registry{"registry"}; // for QA of selections diff --git a/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx b/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx index aabb13137e8..fca1f240721 100644 --- a/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx +++ b/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx @@ -129,7 +129,7 @@ struct HfCandidateSelectorToOmegaPi { TrackSelectorKa selectorKaon; using TracksSel = soa::Join; - using TracksSelLf = soa::Join; + using TracksSelLf = soa::Join; HistogramRegistry registry{"registry"}; // for QA of selections From 4f88d423c11d62a2cce794067e4f94858147929b Mon Sep 17 00:00:00 2001 From: ZFederica Date: Fri, 16 Aug 2024 23:22:01 +0200 Subject: [PATCH 05/10] Retrieve eta at the correct vertex --- PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx b/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx index 8e48a7faeed..802a47c8ada 100644 --- a/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx +++ b/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx @@ -325,8 +325,8 @@ struct HfCandidateCreatorXic0Omegac0 { //-------------------------- V0 info--------------------------- // pseudorapidity - float pseudorapV0Dau0 = trackV0Dau0.eta(); - float pseudorapV0Dau1 = trackV0Dau1.eta(); + float pseudorapV0Dau0 = casc.positiveeta(); + float pseudorapV0Dau1 = casc.negativeeta(); // info from LF table std::array pVecV0 = {casc.pxlambda(), casc.pylambda(), casc.pzlambda()}; @@ -336,7 +336,7 @@ struct HfCandidateCreatorXic0Omegac0 { //-------------------reconstruct cascade track------------------ // pseudorapidity - float pseudorapCascBachelor = trackCascDauCharged.eta(); + float pseudorapCascBachelor = casc.bacheloreta(); // info from LF table std::array vertexCasc = {casc.x(), casc.y(), casc.z()}; From 7e4ce86f548a96ac3c5cc38d31f4e3453655af4a Mon Sep 17 00:00:00 2001 From: ZFederica Date: Sat, 17 Aug 2024 15:45:14 +0200 Subject: [PATCH 06/10] Fix build error --- .../candidateCreatorXic0Omegac0.cxx | 16 ++++++++-------- 1 file changed, 8 insertions(+), 8 deletions(-) diff --git a/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx b/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx index 802a47c8ada..af973cc67cb 100644 --- a/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx +++ b/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx @@ -1293,14 +1293,14 @@ struct HfCandidateCreatorXic0Omegac0Mc { std::vector idxBhadMothers{}; auto arrayDaughters = std::array{candidate.template bachelorFromCharmBaryon_as(), // bachelor <- charm baryon - candidate.template bachelor_as(), // bachelor <- cascade - candidate.template posTrack_as(), // p <- lambda - candidate.template negTrack_as()}; // pi <- lambda - auto arrayDaughtersCasc = std::array{candidate.template bachelor_as(), - candidate.template posTrack_as(), - candidate.template negTrack_as()}; - auto arrayDaughtersV0 = std::array{candidate.template posTrack_as(), - candidate.template negTrack_as()}; + candidate.template bachelor_as(), // bachelor <- cascade + candidate.template posTrack_as(), // p <- lambda + candidate.template negTrack_as()}; // pi <- lambda + auto arrayDaughtersCasc = std::array{candidate.template bachelor_as(), + candidate.template posTrack_as(), + candidate.template negTrack_as()}; + auto arrayDaughtersV0 = std::array{candidate.template posTrack_as(), + candidate.template negTrack_as()}; // Check whether the particle is from background events. If so, reject it. if (rejectBackground) { From dd7930f01cc4a67707462bd89517ae1333c8e24e Mon Sep 17 00:00:00 2001 From: ZFederica Date: Sat, 17 Aug 2024 15:46:38 +0200 Subject: [PATCH 07/10] Fix build error --- PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx b/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx index af973cc67cb..c0b53001386 100644 --- a/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx +++ b/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx @@ -417,9 +417,9 @@ struct HfCandidateCreatorXic0Omegac0 { float dcaxyV0Dau0 = impactParameterV0Dau0.getY(); float dcaxyV0Dau1 = impactParameterV0Dau1.getY(); float dcaxyCascBachelor = impactParameterCascDauCharged.getY(); - float dcazV0Dau0 = impactParameterV0Dau0.dcaZ(); - float dcazV0Dau1 = impactParameterV0Dau1.dcaZ(); - float dcazCascBachelor = impactParameterCascDauCharged.dcaZ(); + float dcazV0Dau0 = impactParameterV0Dau0.getZ(); + float dcazV0Dau1 = impactParameterV0Dau1.getZ(); + float dcazCascBachelor = impactParameterCascDauCharged.getZ(); // impact parameters o2::dataformats::DCA impactParameterCasc; From 0af39dd113727772955d99ee49216934a91a36ec Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Sat, 17 Aug 2024 13:47:06 +0000 Subject: [PATCH 08/10] Please consider the following formatting changes --- PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx b/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx index c0b53001386..59ada87e466 100644 --- a/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx +++ b/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx @@ -1293,9 +1293,9 @@ struct HfCandidateCreatorXic0Omegac0Mc { std::vector idxBhadMothers{}; auto arrayDaughters = std::array{candidate.template bachelorFromCharmBaryon_as(), // bachelor <- charm baryon - candidate.template bachelor_as(), // bachelor <- cascade - candidate.template posTrack_as(), // p <- lambda - candidate.template negTrack_as()}; // pi <- lambda + candidate.template bachelor_as(), // bachelor <- cascade + candidate.template posTrack_as(), // p <- lambda + candidate.template negTrack_as()}; // pi <- lambda auto arrayDaughtersCasc = std::array{candidate.template bachelor_as(), candidate.template posTrack_as(), candidate.template negTrack_as()}; From 89a2a957f96cd6883015f0cb91e3502dd5aecea4 Mon Sep 17 00:00:00 2001 From: ZFederica Date: Sat, 17 Aug 2024 23:46:22 +0200 Subject: [PATCH 09/10] Fix track table --- .../candidateCreatorXic0Omegac0.cxx | 83 ++++++++----------- 1 file changed, 35 insertions(+), 48 deletions(-) diff --git a/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx b/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx index 59ada87e466..ce353e4b90f 100644 --- a/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx +++ b/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx @@ -1191,7 +1191,7 @@ struct HfCandidateCreatorXic0Omegac0Mc { // Configuration o2::framework::Configurable rejectBackground{"rejectBackground", true, "Reject particles from background events"}; - using lfTracksWMc = soa::Join; + using MyTracksWMc = soa::Join; using McCollisionsNoCents = soa::Join; using McCollisionsFT0Cs = soa::Join; using McCollisionsFT0Ms = soa::Join; @@ -1251,8 +1251,7 @@ struct HfCandidateCreatorXic0Omegac0Mc { template void runXic0Omegac0Mc(TMyRecoCand const& candidates, - aod::TracksWMc const&, - lfTracksWMc const&, + MyTracksWMc const&, aod::McParticles const& mcParticles, Colls const& collsWithMcLabels, aod::McCollisions const& mcCollisions, @@ -1292,15 +1291,15 @@ struct HfCandidateCreatorXic0Omegac0Mc { collisionMatched = false; std::vector idxBhadMothers{}; - auto arrayDaughters = std::array{candidate.template bachelorFromCharmBaryon_as(), // bachelor <- charm baryon - candidate.template bachelor_as(), // bachelor <- cascade - candidate.template posTrack_as(), // p <- lambda - candidate.template negTrack_as()}; // pi <- lambda - auto arrayDaughtersCasc = std::array{candidate.template bachelor_as(), - candidate.template posTrack_as(), - candidate.template negTrack_as()}; - auto arrayDaughtersV0 = std::array{candidate.template posTrack_as(), - candidate.template negTrack_as()}; + auto arrayDaughters = std::array{candidate.template bachelorFromCharmBaryon_as(), // bachelor <- charm baryon + candidate.template bachelor_as(), // bachelor <- cascade + candidate.template posTrack_as(), // p <- lambda + candidate.template negTrack_as()}; // pi <- lambda + auto arrayDaughtersCasc = std::array{candidate.template bachelor_as(), + candidate.template posTrack_as(), + candidate.template negTrack_as()}; + auto arrayDaughtersV0 = std::array{candidate.template posTrack_as(), + candidate.template negTrack_as()}; // Check whether the particle is from background events. If so, reject it. if (rejectBackground) { @@ -1722,146 +1721,134 @@ struct HfCandidateCreatorXic0Omegac0Mc { PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processDoNoMc, "Do not run any MC process function", true); void processMcXicToXiPi(aod::HfCandToXiPi const& candidates, - aod::TracksWMc const& tracks, - lfTracksWMc const& lfTracks, + MyTracksWMc const& tracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsNoCents const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcXicToXiPi, "Run Xic0 to xi pi MC process function - no centrality", false); void processMcXicToXiPiFT0m(aod::HfCandToXiPi const& candidates, - aod::TracksWMc const& tracks, - lfTracksWMc const& lfTracks, + MyTracksWMc const& tracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsFT0Ms const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcXicToXiPiFT0m, "Run Xic0 to xi pi MC process function - FT0M", false); void processMcXicToXiPiFT0c(aod::HfCandToXiPi const& candidates, - aod::TracksWMc const& tracks, - lfTracksWMc const& lfTracks, + MyTracksWMc const& tracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsFT0Cs const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcXicToXiPiFT0c, "Run Xic0 to xi pi MC process function - FT0C", false); void processMcOmegacToXiPi(aod::HfCandToXiPi const& candidates, - aod::TracksWMc const& tracks, - lfTracksWMc const& lfTracks, + MyTracksWMc const& tracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsNoCents const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcOmegacToXiPi, "Run Omegac0 to xi pi MC process function - no centrality", false); void processMcOmegacToXiPiFT0m(aod::HfCandToXiPi const& candidates, - aod::TracksWMc const& tracks, - lfTracksWMc const& lfTracks, + MyTracksWMc const& tracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsFT0Ms const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcOmegacToXiPiFT0m, "Run Omegac0 to xi pi MC process function - FT0M", false); void processMcOmegacToXiPiFT0c(aod::HfCandToXiPi const& candidates, - aod::TracksWMc const& tracks, - lfTracksWMc const& lfTracks, + MyTracksWMc const& tracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsFT0Cs const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcOmegacToXiPiFT0c, "Run Omegac0 to xi pi MC process function - FT0C", false); void processMcOmegacToOmegaPi(aod::HfCandToOmegaPi const& candidates, - aod::TracksWMc const& tracks, - lfTracksWMc const& lfTracks, + MyTracksWMc const& tracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsNoCents const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcOmegacToOmegaPi, "Run Omegac0 to omega pi MC process function - no centrality", false); void processMcOmegacToOmegaPiFT0m(aod::HfCandToOmegaPi const& candidates, - aod::TracksWMc const& tracks, - lfTracksWMc const& lfTracks, + MyTracksWMc const& tracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsFT0Ms const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcOmegacToOmegaPiFT0m, "Run Omegac0 to omega pi MC process function - FT0M", false); void processMcOmegacToOmegaPiFT0c(aod::HfCandToOmegaPi const& candidates, - aod::TracksWMc const& tracks, - lfTracksWMc const& lfTracks, + MyTracksWMc const& tracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsFT0Cs const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcOmegacToOmegaPiFT0c, "Run Omegac0 to omega pi MC process function - FT0C", false); void processMcOmegacToOmegaK(aod::HfCandToOmegaK const& candidates, - aod::TracksWMc const& tracks, - lfTracksWMc const& lfTracks, + MyTracksWMc const& tracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsNoCents const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcOmegacToOmegaK, "Run Omegac0 to omega K MC process function - no centrality", false); void processMcOmegacToOmegaKFT0m(aod::HfCandToOmegaK const& candidates, - aod::TracksWMc const& tracks, - lfTracksWMc const& lfTracks, + MyTracksWMc const& tracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsFT0Ms const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcOmegacToOmegaKFT0m, "Run Omegac0 to omega K MC process function - FT0M", false); void processMcOmegacToOmegaKFT0c(aod::HfCandToOmegaK const& candidates, - aod::TracksWMc const& tracks, - lfTracksWMc const& lfTracks, + MyTracksWMc const& tracks, aod::McParticles const& mcParticles, aod::McCollisions const& mcColls, McCollisionsFT0Cs const& collsWithMcLabels, BCsInfo const& bcs) { - runXic0Omegac0Mc(candidates, tracks, lfTracks, mcParticles, collsWithMcLabels, mcColls, bcs); + runXic0Omegac0Mc(candidates, tracks, mcParticles, collsWithMcLabels, mcColls, bcs); } PROCESS_SWITCH(HfCandidateCreatorXic0Omegac0Mc, processMcOmegacToOmegaKFT0c, "Run Omegac0 to omega K MC process function - FT0C", false); From fdd18723ce3362aa8460f32bee95018d24f0eb9c Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Sat, 17 Aug 2024 21:46:58 +0000 Subject: [PATCH 10/10] Please consider the following formatting changes --- PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx b/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx index ce353e4b90f..961df7c6b87 100644 --- a/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx +++ b/PWGHF/TableProducer/candidateCreatorXic0Omegac0.cxx @@ -1292,9 +1292,9 @@ struct HfCandidateCreatorXic0Omegac0Mc { std::vector idxBhadMothers{}; auto arrayDaughters = std::array{candidate.template bachelorFromCharmBaryon_as(), // bachelor <- charm baryon - candidate.template bachelor_as(), // bachelor <- cascade - candidate.template posTrack_as(), // p <- lambda - candidate.template negTrack_as()}; // pi <- lambda + candidate.template bachelor_as(), // bachelor <- cascade + candidate.template posTrack_as(), // p <- lambda + candidate.template negTrack_as()}; // pi <- lambda auto arrayDaughtersCasc = std::array{candidate.template bachelor_as(), candidate.template posTrack_as(), candidate.template negTrack_as()};