From f35031c474c4afe2d3844efb57d744b70d21088b Mon Sep 17 00:00:00 2001 From: Francesco Mazzaschi Date: Fri, 4 Oct 2024 15:03:31 +0200 Subject: [PATCH 1/2] Move abs to std::abs for hypernuclei --- PWGLF/TableProducer/Nuspex/hyperRecoTask.cxx | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/PWGLF/TableProducer/Nuspex/hyperRecoTask.cxx b/PWGLF/TableProducer/Nuspex/hyperRecoTask.cxx index a42d65a4296..0068d7a506e 100644 --- a/PWGLF/TableProducer/Nuspex/hyperRecoTask.cxx +++ b/PWGLF/TableProducer/Nuspex/hyperRecoTask.cxx @@ -490,7 +490,7 @@ struct hyperRecoTask { o2::base::Propagator::Instance()->propagateToDCABxByBz({collision.posX(), collision.posY(), collision.posZ()}, piTrackCov, 2.f, fitter.getMatCorrType(), &dcaInfo); hypCand.piDCAXY = dcaInfo[0]; - if (abs(hypCand.piDCAXY) < dcaToPvPion || abs(hypCand.he3DCAXY) < dcaToPvHe) { + if (std::abs(hypCand.piDCAXY) < dcaToPvPion || std::abs(hypCand.he3DCAXY) < dcaToPvHe) { return; } @@ -597,7 +597,7 @@ struct hyperRecoTask { for (auto& piMother : mcTrackPi.mothers_as()) { if (heMother.globalIndex() != piMother.globalIndex()) continue; - if (abs(mcTrackHe.pdgCode()) != heDauPdg || abs(mcTrackPi.pdgCode()) != 211) + if (std::abs(mcTrackHe.pdgCode()) != heDauPdg || std::abs(mcTrackPi.pdgCode()) != 211) continue; if (std::abs(heMother.pdgCode()) != hyperPdg) continue; From c7ab430e774fe474ec72f79c40df4e77d7aa0058 Mon Sep 17 00:00:00 2001 From: Francesco Mazzaschi Date: Fri, 4 Oct 2024 15:07:28 +0200 Subject: [PATCH 2/2] Add kinks --- PWGLF/TableProducer/Nuspex/hyperKinkRecoTask.cxx | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/PWGLF/TableProducer/Nuspex/hyperKinkRecoTask.cxx b/PWGLF/TableProducer/Nuspex/hyperKinkRecoTask.cxx index 658ab180d67..8e46e3326c1 100644 --- a/PWGLF/TableProducer/Nuspex/hyperKinkRecoTask.cxx +++ b/PWGLF/TableProducer/Nuspex/hyperKinkRecoTask.cxx @@ -332,7 +332,7 @@ struct hyperKinkRecoTask { bool isGoodTPCCand = false; if (candidate.itsNClsInnerBarrel() == 0 && candidate.itsNCls() < 4 && candidate.tpcNClsCrossedRows() >= 70 && candidate.tpcChi2NCl() < 4.f && - candidate.tpcNClsCrossedRows() > 0.8 * candidate.tpcNClsFindable() && candidate.tpcNClsFound() > 80 && abs(nSigmaTrit) < nSigmaTPCCutTrit) { + candidate.tpcNClsCrossedRows() > 0.8 * candidate.tpcNClsFindable() && candidate.tpcNClsFound() > 80 && std::abs(nSigmaTrit) < nSigmaTPCCutTrit) { isGoodTPCCand = true; } @@ -340,7 +340,7 @@ struct hyperKinkRecoTask { return false; } - if (candidate.hasTOF() && abs(nSigmaTOFTrit) > nSigmaTOFCutTrit) { + if (candidate.hasTOF() && std::abs(nSigmaTOFTrit) > nSigmaTOFCutTrit) { return false; } @@ -394,7 +394,7 @@ struct hyperKinkRecoTask { gpu::gpustd::array dcaInfoHyp; o2::base::Propagator::Instance()->propagateToDCABxByBz({primaryVertex.getX(), primaryVertex.getY(), primaryVertex.getZ()}, trackParCovHyperPV, 2.f, static_cast(cfgMaterialCorrection.value), &dcaInfoHyp); - if (abs(dcaInfoHyp[0]) > maxDCAHypToPV) { + if (std::abs(dcaInfoHyp[0]) > maxDCAHypToPV) { continue; } @@ -412,7 +412,7 @@ struct hyperKinkRecoTask { // propagate to PV gpu::gpustd::array dcaInfoTrit; o2::base::Propagator::Instance()->propagateToDCABxByBz({primaryVertex.getX(), primaryVertex.getY(), primaryVertex.getZ()}, trackParCovTrit, 2.f, static_cast(cfgMaterialCorrection.value), &dcaInfoTrit); - if (abs(dcaInfoTrit[0]) < minDCATritToPV) { + if (std::abs(dcaInfoTrit[0]) < minDCATritToPV) { continue; } @@ -492,13 +492,13 @@ struct hyperKinkRecoTask { auto mcTrackHyper = mcLabHyper.mcParticle_as(); auto mcTrackTrit = mcLabTrit.mcParticle_as(); - if (abs(mcTrackHyper.pdgCode()) != hyperPdg || abs(mcTrackTrit.pdgCode()) != tritDauPdg) { + if (std::abs(mcTrackHyper.pdgCode()) != hyperPdg || std::abs(mcTrackTrit.pdgCode()) != tritDauPdg) { continue; } auto tritIdx = mcTrackTrit.globalIndex(); kinkCand.isSignal = false; for (auto& dauMCTracks : mcTrackHyper.daughters_as()) { - if (abs(dauMCTracks.pdgCode()) == tritDauPdg) { + if (std::abs(dauMCTracks.pdgCode()) == tritDauPdg) { if (dauMCTracks.globalIndex() == tritIdx) { kinkCand.isSignal = true; break;