From c5eded1ff68e986184f1c562aed782c219dc728a Mon Sep 17 00:00:00 2001 From: yuanzhe Date: Fri, 27 Jun 2025 16:02:56 +0200 Subject: [PATCH 1/2] Update datamodel of hyphe4s analysis --- PWGLF/DataModel/LFHyperhelium4sigmaTables.h | 15 +++- .../Nuspex/hyperhelium4sigmaRecoTask.cxx | 69 ++++++++++++++----- 2 files changed, 66 insertions(+), 18 deletions(-) diff --git a/PWGLF/DataModel/LFHyperhelium4sigmaTables.h b/PWGLF/DataModel/LFHyperhelium4sigmaTables.h index 5d50a007f66..4e6d28085a3 100644 --- a/PWGLF/DataModel/LFHyperhelium4sigmaTables.h +++ b/PWGLF/DataModel/LFHyperhelium4sigmaTables.h @@ -30,6 +30,9 @@ DECLARE_SOA_COLUMN(ZPrimVtx, zPrimVtx, float); // Pri DECLARE_SOA_COLUMN(XDecVtx, xDecVtx, float); // Decay vertex of the candidate (x direction) DECLARE_SOA_COLUMN(YDecVtx, yDecVtx, float); // Decay vertex of the candidate (y direction) DECLARE_SOA_COLUMN(ZDecVtx, zDecVtx, float); // Decay vertex of the candidate (z direction) +DECLARE_SOA_COLUMN(XMoth, xMoth, float); // X of the mother track at the radii of ITS layer which has the outermost update +DECLARE_SOA_COLUMN(YMoth, yMoth, float); // Y of the mother track at the radii of ITS layer which has the outermost update +DECLARE_SOA_COLUMN(ZMoth, zMoth, float); // Z of the mother track at the radii of ITS layer which has the outermost update DECLARE_SOA_COLUMN(PxMoth, pxMoth, float); //! Px of the mother track at the decay vertex DECLARE_SOA_COLUMN(PyMoth, pyMoth, float); //! Py of the mother track at the decay vertex DECLARE_SOA_COLUMN(PzMoth, pzMoth, float); //! Pz of the mother track at the decay vertex @@ -56,6 +59,9 @@ DECLARE_SOA_COLUMN(TrueZDecVtx, trueZDecVtx, float); // true z decay DECLARE_SOA_COLUMN(GenPxMoth, genPxMoth, float); // Generated px of the mother track DECLARE_SOA_COLUMN(GenPyMoth, genPyMoth, float); // Generated py of the mother track DECLARE_SOA_COLUMN(GenPzMoth, genPzMoth, float); // Generated pz of the mother track +DECLARE_SOA_COLUMN(TruePxMoth, truePxMoth, float); // true px of the mother track at the decay vertex +DECLARE_SOA_COLUMN(TruePyMoth, truePyMoth, float); // true py of the mother track at the decay vertex +DECLARE_SOA_COLUMN(TruePzMoth, truePzMoth, float); // true pz of the mother track at the decay vertex DECLARE_SOA_COLUMN(GenPxAlpha, genPxAlpha, float); // true px of the daughter alpha track DECLARE_SOA_COLUMN(GenPyAlpha, genPyAlpha, float); // true py of the daughter alpha track DECLARE_SOA_COLUMN(GenPzAlpha, genPzAlpha, float); // true pz of the daughter alpha track @@ -68,7 +74,9 @@ DECLARE_SOA_TABLE(He4S2BCands, "AOD", "HE4S2BCANDS", o2::soa::Index<>, he4scand::XPrimVtx, he4scand::YPrimVtx, he4scand::ZPrimVtx, he4scand::XDecVtx, he4scand::YDecVtx, he4scand::ZDecVtx, - he4scand::IsMatter, he4scand::PxMoth, he4scand::PyMoth, he4scand::PzMoth, + he4scand::IsMatter, + he4scand::XMoth, he4scand::YMoth, he4scand::ZMoth, + he4scand::PxMoth, he4scand::PyMoth, he4scand::PzMoth, he4scand::PxAlpha, he4scand::PyAlpha, he4scand::PzAlpha, he4scand::DcaMothPv, he4scand::DcaAlphaPv, he4scand::DcaKinkTopo, he4scand::ItsChi2Moth, he4scand::ItsClusterSizesMoth, he4scand::ItsClusterSizesAlpha, @@ -78,7 +86,9 @@ DECLARE_SOA_TABLE(MCHe4S2BCands, "AOD", "MCHE4S2BCANDS", o2::soa::Index<>, he4scand::XPrimVtx, he4scand::YPrimVtx, he4scand::ZPrimVtx, he4scand::XDecVtx, he4scand::YDecVtx, he4scand::ZDecVtx, - he4scand::IsMatter, he4scand::PxMoth, he4scand::PyMoth, he4scand::PzMoth, + he4scand::IsMatter, + he4scand::XMoth, he4scand::YMoth, he4scand::ZMoth, + he4scand::PxMoth, he4scand::PyMoth, he4scand::PzMoth, he4scand::PxAlpha, he4scand::PyAlpha, he4scand::PzAlpha, he4scand::DcaMothPv, he4scand::DcaAlphaPv, he4scand::DcaKinkTopo, he4scand::ItsChi2Moth, he4scand::ItsClusterSizesMoth, he4scand::ItsClusterSizesAlpha, @@ -86,6 +96,7 @@ DECLARE_SOA_TABLE(MCHe4S2BCands, "AOD", "MCHE4S2BCANDS", he4scand::IsSignal, he4scand::IsSignalReco, he4scand::IsCollReco, he4scand::IsSurvEvSelection, he4scand::TrueXDecVtx, he4scand::TrueYDecVtx, he4scand::TrueZDecVtx, he4scand::GenPxMoth, he4scand::GenPyMoth, he4scand::GenPzMoth, + he4scand::TruePxMoth, he4scand::TruePyMoth, he4scand::TruePzMoth, he4scand::GenPxAlpha, he4scand::GenPyAlpha, he4scand::GenPzAlpha, he4scand::IsMothReco, he4scand::RecoPtMoth, he4scand::RecoPzMoth); diff --git a/PWGLF/TableProducer/Nuspex/hyperhelium4sigmaRecoTask.cxx b/PWGLF/TableProducer/Nuspex/hyperhelium4sigmaRecoTask.cxx index 3e35c99e2df..2c48ea924f1 100644 --- a/PWGLF/TableProducer/Nuspex/hyperhelium4sigmaRecoTask.cxx +++ b/PWGLF/TableProducer/Nuspex/hyperhelium4sigmaRecoTask.cxx @@ -68,7 +68,7 @@ enum DaughterType { namespace { -// constexpr std::array LayerRadii{2.33959f, 3.14076f, 3.91924f, 19.6213f, 24.5597f, 34.388f, 39.3329f}; +constexpr std::array LayerRadii{2.33959f, 3.14076f, 3.91924f, 19.6213f, 24.5597f, 34.388f, 39.3329f}; constexpr int kITSLayers = 7; constexpr int kITSInnerBarrelLayers = 3; // constexpr int kITSOuterBarrelLayers = 4; @@ -231,6 +231,7 @@ struct Hyphe4sCandidate { std::array primVtx = {0.0f, 0.0f, 0.0f}; std::array decVtx = {0.0f, 0.0f, 0.0f}; + std::array lastPosMoth = {0.0f, 0.0f, 0.0f}; // last position of mother track at the radii of ITS layer which has the outermost update std::array momMoth = {0.0f, 0.0f, 0.0f}; std::array momDaug = {0.0f, 0.0f, 0.0f}; @@ -251,7 +252,8 @@ struct Hyphe4sCandidate { bool isSurvEvSelection = false; std::array trueDecVtx = {0.0f, 0.0f, 0.0f}; - std::array gMomMoth = {0.0f, 0.0f, 0.0f}; + std::array gMomMoth = {0.0f, 0.0f, 0.0f}; // generated mother momentum + std::array trueMomMoth = {0.0f, 0.0f, 0.0f}; // true mother momentum at decay vertex std::array gMomDau = {0.0f, 0.0f, 0.0f}; bool isMothReco = false; @@ -322,8 +324,8 @@ struct Hyperhelium4sigmaRecoTask { registry.add("h2TrueMotherDiffPyVsRecSVR", ";Reconstruced SV R (cm);#Delta #it{p}_{T} (GeV/#it{c});", HistType::kTH2F, {radiusAxis, diffPyAxis}); registry.add("h2TrueMotherDiffPtVsRecSVR", ";Reconstruced SV R (cm);#Delta #it{p}_{T} (GeV/#it{c});", HistType::kTH2F, {radiusAxis, diffPtAxis}); registry.add("h2TrueMotherDiffPzVsRecSVR", ";Reconstruced SV R (cm);#Delta #it{p}_{z} (GeV/#it{c});", HistType::kTH2F, {radiusAxis, diffPzAxis}); - registry.add("h2TrueMotherDiffTglVsRecSVR", ";Reconstruced SV R (cm);#Delta tan#lambda;", HistType::kTH2F, {radiusAxis, {200, -1.f, 1.f}}); - registry.add("h2TrueMotherDiffEtaVsRecSVR", ";Reconstruced SV R (cm);#Delta #eta;", HistType::kTH2F, {radiusAxis, {200, -1.f, 1.f}}); + registry.add("h2TrueMotherDiffTglVsRecSVR", ";Reconstruced SV R (cm);#Delta tan#lambda;", HistType::kTH2F, {radiusAxis, {200, -0.1, 0.1}}); + registry.add("h2TrueMotherDiffEtaVsRecSVR", ";Reconstruced SV R (cm);#Delta #eta;", HistType::kTH2F, {radiusAxis, {200, -0.1, 0.1}}); registry.add("hDiffDauPx", ";#Delta p_{x} (GeV/#it{c}); ", HistType::kTH1D, {{200, -10, 10}}); registry.add("hDiffDauPy", ";#Delta p_{y} (GeV/#it{c}); ", HistType::kTH1D, {{200, -10, 10}}); registry.add("hDiffDauPz", ";#Delta p_{z} (GeV/#it{c}); ", HistType::kTH1D, {{200, -10, 10}}); @@ -394,6 +396,22 @@ struct Hyperhelium4sigmaRecoTask { hyphe4sCand.itsClusterSizeDau = trackDau.itsClusterSizes(); hyphe4sCand.nSigmaTPCDau = trackDau.tpcNSigmaAl(); hyphe4sCand.nSigmaITSDau = itsResponse.nSigmaITS(trackDau); + + int lastLayerMoth = 0; + for (int i = 6; i >= 0; i--) { + if (trackMoth.itsClusterMap() & (1 << i)) { + lastLayerMoth = i; + break; + } + } + auto trackparMother = getTrackParCov(trackMoth); + o2::base::Propagator::Instance()->PropagateToXBxByBz(trackparMother, LayerRadii[lastLayerMoth]); + std::array vecLab{0.f}; + if (trackparMother.getPosDirGlo(vecLab)) { + hyphe4sCand.lastPosMoth[0] = vecLab[0]; + hyphe4sCand.lastPosMoth[1] = vecLab[1]; + hyphe4sCand.lastPosMoth[2] = vecLab[2]; + } } template @@ -407,7 +425,7 @@ struct Hyperhelium4sigmaRecoTask { } template - void fillCandidateMCInfo(Hyphe4sCandidate& hyphe4sCand, TMCParticle const& mcMothTrack, TMCParticle const& mcDauTrack) + void fillCandidateMCInfo(Hyphe4sCandidate& hyphe4sCand, TMCParticle const& mcMothTrack, TMCParticle const& mcDauTrack, TMCParticle const& mcNeutDauTrack) { hyphe4sCand.trueDecVtx[0] = mcDauTrack.vx(); hyphe4sCand.trueDecVtx[1] = mcDauTrack.vy(); @@ -415,6 +433,9 @@ struct Hyperhelium4sigmaRecoTask { hyphe4sCand.gMomMoth[0] = mcMothTrack.px(); hyphe4sCand.gMomMoth[1] = mcMothTrack.py(); hyphe4sCand.gMomMoth[2] = mcMothTrack.pz(); + hyphe4sCand.trueMomMoth[0] = mcDauTrack.px() + mcNeutDauTrack.px(); + hyphe4sCand.trueMomMoth[1] = mcDauTrack.py() + mcNeutDauTrack.py(); + hyphe4sCand.trueMomMoth[2] = mcDauTrack.pz() + mcNeutDauTrack.pz(); hyphe4sCand.gMomDau[0] = mcDauTrack.px(); hyphe4sCand.gMomDau[1] = mcDauTrack.py(); hyphe4sCand.gMomDau[2] = mcDauTrack.pz(); @@ -448,6 +469,8 @@ struct Hyperhelium4sigmaRecoTask { registry.fill(HIST("h2MassHyperhelium4sigmaPt"), kinkCand.mothSign() * kinkCand.ptMoth(), invMass); registry.fill(HIST("h2NSigmaAlPt"), kinkCand.mothSign() * kinkCand.ptDaug(), dauTrack.tpcNSigmaAl()); + auto bc = collision.bc_as(); + initCCDB(bc); auto motherTrack = kinkCand.trackMoth_as(); Hyphe4sCandidate hyphe4sCand; fillCandidate(hyphe4sCand, collision, kinkCand, motherTrack, dauTrack); @@ -455,7 +478,9 @@ struct Hyperhelium4sigmaRecoTask { outputDataTable( hyphe4sCand.primVtx[0], hyphe4sCand.primVtx[1], hyphe4sCand.primVtx[2], hyphe4sCand.decVtx[0], hyphe4sCand.decVtx[1], hyphe4sCand.decVtx[2], - hyphe4sCand.isMatter, hyphe4sCand.momMoth[0], hyphe4sCand.momMoth[1], hyphe4sCand.momMoth[2], + hyphe4sCand.isMatter, + hyphe4sCand.lastPosMoth[0], hyphe4sCand.lastPosMoth[1], hyphe4sCand.lastPosMoth[2], + hyphe4sCand.momMoth[0], hyphe4sCand.momMoth[1], hyphe4sCand.momMoth[2], hyphe4sCand.momDaug[0], hyphe4sCand.momDaug[1], hyphe4sCand.momDaug[2], hyphe4sCand.dcaXYMothPv, hyphe4sCand.dcaXYDauPv, hyphe4sCand.dcaKinkTopo, hyphe4sCand.chi2ITSMoth, hyphe4sCand.itsClusterSizeMoth, hyphe4sCand.itsClusterSizeDau, @@ -524,6 +549,8 @@ struct Hyperhelium4sigmaRecoTask { registry.fill(HIST("h2MassHyperhelium4sigmaPt"), kinkCand.mothSign() * kinkCand.ptMoth(), invMass); registry.fill(HIST("h2NSigmaAlPt"), kinkCand.mothSign() * kinkCand.ptDaug(), dauTrack.tpcNSigmaAl()); + auto bc = collision.bc_as(); + initCCDB(bc); Hyphe4sCandidate hyphe4sCand; fillCandidate(hyphe4sCand, collision, kinkCand, motherTrack, dauTrack); @@ -531,6 +558,7 @@ struct Hyperhelium4sigmaRecoTask { if (isTrueSignal) { auto mcMotherTrack = motherTrack.mcParticle_as(); auto mcDauTrack = dauTrack.mcParticle_as(); + auto mcNeutTrack = particlesMC.rawIteratorAt(dauIDList[2]); float posDecVtx[3] = {kinkCand.xDecVtx() + collision.posX(), kinkCand.yDecVtx() + collision.posY(), kinkCand.zDecVtx() + collision.posZ()}; float recSVR = std::sqrt(posDecVtx[0] * posDecVtx[0] + posDecVtx[1] * posDecVtx[1]); registry.fill(HIST("hDiffSVx"), posDecVtx[0] - mcDauTrack.vx()); @@ -551,11 +579,9 @@ struct Hyperhelium4sigmaRecoTask { hyphe4sCand.isSignalReco = true; hyphe4sCand.isCollReco = true; hyphe4sCand.isSurvEvSelection = true; - fillCandidateMCInfo(hyphe4sCand, mcMotherTrack, mcDauTrack); + fillCandidateMCInfo(hyphe4sCand, mcMotherTrack, mcDauTrack, mcNeutTrack); mcHe4sIndices.push_back(mcMotherTrack.globalIndex()); - auto bc = collision.bc_as(); - initCCDB(bc); std::array dcaInfo; auto mcMotherTrackPar = getTrackParFromMC(mcMotherTrack); o2::base::Propagator::Instance()->propagateToDCABxByBz({posDecVtx[0], posDecVtx[1], posDecVtx[2]}, mcMotherTrackPar, 2.f, matCorr, &dcaInfo); @@ -591,7 +617,9 @@ struct Hyperhelium4sigmaRecoTask { outputMCTable( hyphe4sCand.primVtx[0], hyphe4sCand.primVtx[1], hyphe4sCand.primVtx[2], hyphe4sCand.decVtx[0], hyphe4sCand.decVtx[1], hyphe4sCand.decVtx[2], - hyphe4sCand.isMatter, hyphe4sCand.momMoth[0], hyphe4sCand.momMoth[1], hyphe4sCand.momMoth[2], + hyphe4sCand.isMatter, + hyphe4sCand.lastPosMoth[0], hyphe4sCand.lastPosMoth[1], hyphe4sCand.lastPosMoth[2], + hyphe4sCand.momMoth[0], hyphe4sCand.momMoth[1], hyphe4sCand.momMoth[2], hyphe4sCand.momDaug[0], hyphe4sCand.momDaug[1], hyphe4sCand.momDaug[2], hyphe4sCand.dcaXYMothPv, hyphe4sCand.dcaXYDauPv, hyphe4sCand.dcaKinkTopo, hyphe4sCand.chi2ITSMoth, hyphe4sCand.itsClusterSizeMoth, hyphe4sCand.itsClusterSizeDau, @@ -599,6 +627,7 @@ struct Hyperhelium4sigmaRecoTask { hyphe4sCand.isSignal, hyphe4sCand.isSignalReco, hyphe4sCand.isCollReco, hyphe4sCand.isSurvEvSelection, hyphe4sCand.trueDecVtx[0], hyphe4sCand.trueDecVtx[1], hyphe4sCand.trueDecVtx[2], hyphe4sCand.gMomMoth[0], hyphe4sCand.gMomMoth[1], hyphe4sCand.gMomMoth[2], + hyphe4sCand.trueMomMoth[0], hyphe4sCand.trueMomMoth[1], hyphe4sCand.trueMomMoth[2], hyphe4sCand.gMomDau[0], hyphe4sCand.gMomDau[1], hyphe4sCand.gMomDau[2], hyphe4sCand.isMothReco, hyphe4sCand.ptMoth, hyphe4sCand.pzMoth); } @@ -615,7 +644,8 @@ struct Hyperhelium4sigmaRecoTask { Hyphe4sCandidate hyphe4sCand; auto mcDauTrack = particlesMC.rawIteratorAt(dauIDList[0]); - fillCandidateMCInfo(hyphe4sCand, mcparticle, mcDauTrack); + auto mcNeutTrack = particlesMC.rawIteratorAt(dauIDList[2]); + fillCandidateMCInfo(hyphe4sCand, mcparticle, mcDauTrack, mcNeutTrack); if (mcPartIndices[mcparticle.globalIndex()] != -1) { auto mothTrack = tracks.rawIteratorAt(mcPartIndices[mcparticle.globalIndex()]); @@ -625,7 +655,9 @@ struct Hyperhelium4sigmaRecoTask { outputMCTable( -1, -1, -1, -1, -1, -1, - -1, -1, -1, -1, + -1, + -1, -1, -1, + -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, @@ -633,6 +665,7 @@ struct Hyperhelium4sigmaRecoTask { true, false, isReconstructedMCCollisions[mcparticle.mcCollisionId()], isSelectedMCCollisions[mcparticle.mcCollisionId()], hyphe4sCand.trueDecVtx[0], hyphe4sCand.trueDecVtx[1], hyphe4sCand.trueDecVtx[2], hyphe4sCand.gMomMoth[0], hyphe4sCand.gMomMoth[1], hyphe4sCand.gMomMoth[2], + hyphe4sCand.trueMomMoth[0], hyphe4sCand.trueMomMoth[1], hyphe4sCand.trueMomMoth[2], hyphe4sCand.gMomDau[0], hyphe4sCand.gMomDau[1], hyphe4sCand.gMomDau[2], hyphe4sCand.isMothReco, hyphe4sCand.ptMoth, hyphe4sCand.pzMoth); } @@ -731,6 +764,7 @@ struct Hyperhelium4sigmaQa { hDauTPCNSigma[kDauTriton] = recoQAHist.add("hDauTritonTPCNSigma", "", HistType::kTH2F, {rigidityAxis, nsigmaAxis}); hDauTPCNSigma[kDauProton] = recoQAHist.add("hDauProtonTPCNSigma", "", HistType::kTH2F, {rigidityAxis, nsigmaAxis}); hDauTPCNSigma[kDauChargedPion] = recoQAHist.add("hDauPionTPCNSigma", "", HistType::kTH2F, {rigidityAxis, nsigmaAxis}); + recoQAHist.add("hDauAlphaITSNSigmaCheck", "", HistType::kTH2F, {rigidityAxis, itsnsigmaAxis}); hRecoMotherCounter = recoQAHist.add("hRecoMotherCounter", "", HistType::kTH1F, {{9, 0.f, 9.f}}); hRecoDauAlphaCounter = recoQAHist.add("hRecoDauAlphaCounter", "", HistType::kTH1F, {{9, 0.f, 9.f}}); @@ -750,8 +784,8 @@ struct Hyperhelium4sigmaQa { recoQAHist.add("hMotherITSCls", "", HistType::kTH1F, {{8, 0.f, 8.f}}); recoQAHist.add("hDauAlphaIsPVContributer", "", HistType::kTH1F, {{2, 0.f, 2.f}}); recoQAHist.add("hDauAlphaITSCls", "", HistType::kTH1F, {{8, 0.f, 8.f}}); - recoQAHist.add("hDauAlphaPVsITSNSigma", "", HistType::kTH2F, {rigidityAxis, itsnsigmaAxis}); - recoQAHist.add("h2BCandDauAlphaPVsITSNSigma", "", HistType::kTH2F, {rigidityAxis, nsigmaAxis}); + recoQAHist.add("hDauAlphaITSNSigma", "", HistType::kTH2F, {rigidityAxis, itsnsigmaAxis}); + recoQAHist.add("hReco2BDauAlphaPVsITSNSigma", "", HistType::kTH2F, {rigidityAxis, itsnsigmaAxis}); recoQAHist.add("hReco2BCandidateCount", "", HistType::kTH1F, {{4, 0.f, 4.f}}); } } @@ -950,6 +984,9 @@ struct Hyperhelium4sigmaQa { if (track.hasTPC()) { hDauTPCNSigma[type]->Fill(track.p() * track.sign(), tpcNSigma); } + if (type == kDauAlpha && track.itsNCls() > kITSLayers - 2) { + recoQAHist.fill(HIST("hDauAlphaITSNSigmaCheck"), track.p() * track.sign(), itsResponse.nSigmaITS(track)); + } } } } @@ -1000,7 +1037,7 @@ struct Hyperhelium4sigmaQa { float itsNSigma = itsResponse.nSigmaITS(daughterTrack); if (daughterTrack.hasITS()) { - recoQAHist.fill(HIST("hDauAlphaPVsITSNSigma"), daughterTrack.sign() * daughterTrack.p(), itsNSigma); + recoQAHist.fill(HIST("hDauAlphaITSNSigma"), daughterTrack.sign() * daughterTrack.p(), itsNSigma); recoQAHist.fill(HIST("hDauAlphaITSCls"), daughterTrack.itsNCls()); } @@ -1013,7 +1050,7 @@ struct Hyperhelium4sigmaQa { if (isMoth && isDaug) { recoQAHist.fill(HIST("hReco2BCandidateCount"), 3.5); - recoQAHist.fill(HIST("h2BCandDauAlphaPVsITSNSigma"), daughterTrack.sign() * daughterTrack.p(), itsNSigma); + recoQAHist.fill(HIST("hReco2BDauAlphaPVsITSNSigma"), daughterTrack.sign() * daughterTrack.p(), itsNSigma); } } } From 462f1a6c93c042ba352db06dd7ae27dabaf06cf9 Mon Sep 17 00:00:00 2001 From: yuanzhe Date: Fri, 27 Jun 2025 16:10:14 +0200 Subject: [PATCH 2/2] Fix processData --- PWGLF/TableProducer/Nuspex/hyperhelium4sigmaRecoTask.cxx | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/PWGLF/TableProducer/Nuspex/hyperhelium4sigmaRecoTask.cxx b/PWGLF/TableProducer/Nuspex/hyperhelium4sigmaRecoTask.cxx index 2c48ea924f1..92859fe0ae7 100644 --- a/PWGLF/TableProducer/Nuspex/hyperhelium4sigmaRecoTask.cxx +++ b/PWGLF/TableProducer/Nuspex/hyperhelium4sigmaRecoTask.cxx @@ -441,7 +441,7 @@ struct Hyperhelium4sigmaRecoTask { hyphe4sCand.gMomDau[2] = mcDauTrack.pz(); } - void processData(CollisionsFull const& collisions, aod::KinkCands const& KinkCands, FullTracksExtIU const&) + void processData(CollisionsFull const& collisions, aod::KinkCands const& KinkCands, FullTracksExtIU const&, aod::BCs const&) { for (const auto& collision : collisions) { registry.fill(HIST("hEventCounter"), 0);