From eba1c0591cb057250df101e50a4fdcebc389590a Mon Sep 17 00:00:00 2001 From: Chiara De Martin Date: Thu, 17 Apr 2025 15:53:43 +0200 Subject: [PATCH 1/5] add abs for DCAtoPV --- EventFiltering/PWGLF/strangenessFilter.cxx | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/EventFiltering/PWGLF/strangenessFilter.cxx b/EventFiltering/PWGLF/strangenessFilter.cxx index 59f7f89c14b..ffaeb0cb50b 100644 --- a/EventFiltering/PWGLF/strangenessFilter.cxx +++ b/EventFiltering/PWGLF/strangenessFilter.cxx @@ -740,9 +740,9 @@ struct strangenessFilter { QAHistosTopologicalVariables.fill(HIST("hDCAV0ToPVXi"), DCAV0ToPV); QAHistosTopologicalVariables.fill(HIST("hDCAV0DaughtersXi"), mStraHelper.cascade.v0DaughterDCA); QAHistosTopologicalVariables.fill(HIST("hDCACascDaughtersXi"), mStraHelper.cascade.cascadeDaughterDCA); - QAHistosTopologicalVariables.fill(HIST("hDCABachToPVXi"), mStraHelper.cascade.bachelorDCAxy); - QAHistosTopologicalVariables.fill(HIST("hDCAPosToPVXi"), mStraHelper.cascade.positiveDCAxy); - QAHistosTopologicalVariables.fill(HIST("hDCANegToPVXi"), mStraHelper.cascade.negativeDCAxy); + QAHistosTopologicalVariables.fill(HIST("hDCABachToPVXi"), std::fabs(mStraHelper.cascade.bachelorDCAxy)); + QAHistosTopologicalVariables.fill(HIST("hDCAPosToPVXi"), std::fabs(mStraHelper.cascade.positiveDCAxy)); + QAHistosTopologicalVariables.fill(HIST("hDCANegToPVXi"), std::fabs(mStraHelper.cascade.negativeDCAxy)); QAHistosTopologicalVariables.fill(HIST("hInvMassLambdaXi"), LambdaMass); if (doextraQA) { From 0e5b719d4280e62cebf772ec228056bbcf3d9be6 Mon Sep 17 00:00:00 2001 From: Chiara De Martin Date: Thu, 17 Apr 2025 18:10:36 +0200 Subject: [PATCH 2/5] setting magnetic field in the strangeness helper --- EventFiltering/PWGLF/strangenessFilter.cxx | 1 + 1 file changed, 1 insertion(+) diff --git a/EventFiltering/PWGLF/strangenessFilter.cxx b/EventFiltering/PWGLF/strangenessFilter.cxx index ffaeb0cb50b..68999bf1939 100644 --- a/EventFiltering/PWGLF/strangenessFilter.cxx +++ b/EventFiltering/PWGLF/strangenessFilter.cxx @@ -419,6 +419,7 @@ struct strangenessFilter { mDCAFitter.setMinParamChange(minParamChange); mDCAFitter.setMinRelChi2Change(minRelChi2Change); mDCAFitter.setUseAbsDCA(useAbsDCA); + mStraHelper.fitter.setBz(mBz); } if (!mStraHelper.lut) { /// done only once ccdb->setURL(ccdbUrl); From 263725697acb7f589305772ae81d9288259fc871 Mon Sep 17 00:00:00 2001 From: Chiara De Martin Date: Thu, 17 Apr 2025 18:25:44 +0200 Subject: [PATCH 3/5] fill missing bins in histo QA --- EventFiltering/PWGLF/strangenessFilter.cxx | 30 ++++++++++++++++++++-- 1 file changed, 28 insertions(+), 2 deletions(-) diff --git a/EventFiltering/PWGLF/strangenessFilter.cxx b/EventFiltering/PWGLF/strangenessFilter.cxx index 68999bf1939..061bc63526c 100644 --- a/EventFiltering/PWGLF/strangenessFilter.cxx +++ b/EventFiltering/PWGLF/strangenessFilter.cxx @@ -668,13 +668,18 @@ struct strangenessFilter { } hCandidate->Fill(5.5); } - hCandidate->Fill(7.5); + hCandidate->Fill(6.5); //OLD: eta dau (selection now applied in strangeness helper) + hCandidate->Fill(7.5); //OLD: bachtopv (selection now applied in strangeness helper) // not striclty needed as selection are applied beforehand - just as QA (no change in number expected) if (Cascv0radius < v0radius) { continue; } hCandidate->Fill(8.5); + if (Casccascradius < cascradius) { + continue; + } + hCandidate->Fill(9.5); if (v0DauCPA < v0cospa) { continue; } @@ -683,7 +688,14 @@ struct strangenessFilter { continue; } hCandidate->Fill(11.5); - + if (mStraHelper.cascade.cascadeDaughterDCA > dcacascdau) { + continue; + } + hCandidate->Fill(12.5); + if (std::fabs(LambdaMass - constants::physics::MassLambda) > masslambdalimit) { + continue; + } + hCandidate->Fill(13.5); if (std::fabs(etaCasc) > eta) { continue; } @@ -696,6 +708,20 @@ struct strangenessFilter { } hCandidate->Fill(15.5); + // Fill selections QA for Xi + if (cascCPA > casccospaxi) { + hCandidate->Fill(16.5); + if (cascCPA > dcav0topv) { + hCandidate->Fill(17.5); + if (xiproperlifetime < properlifetimefactor * ctauxi) { + hCandidate->Fill(18.5); + if (std::fabs(yXi) < rapidity) { + hCandidate->Fill(19.5); + } + } + } + } + const auto deltaMassXi = useSigmaBasedMassCutXi ? getMassWindow(stfilter::species::Xi, ptCasc) : ximasswindow; const auto deltaMassOmega = useSigmaBasedMassCutOmega ? getMassWindow(stfilter::species::Omega, ptCasc) : omegamasswindow; From a650642814b94c3f2a62c85c3c48a20fdb17422e Mon Sep 17 00:00:00 2001 From: Chiara De Martin Date: Thu, 17 Apr 2025 18:34:05 +0200 Subject: [PATCH 4/5] small histo fix --- EventFiltering/PWGLF/strangenessFilter.cxx | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/EventFiltering/PWGLF/strangenessFilter.cxx b/EventFiltering/PWGLF/strangenessFilter.cxx index 061bc63526c..319f1d5ca35 100644 --- a/EventFiltering/PWGLF/strangenessFilter.cxx +++ b/EventFiltering/PWGLF/strangenessFilter.cxx @@ -225,7 +225,7 @@ struct strangenessFilter { hProcessedEvents->GetXaxis()->SetBinLabel(16, aod::filtering::OmegaXi::columnLabel()); hCandidate->GetXaxis()->SetBinLabel(1, "All"); - hCandidate->GetXaxis()->SetBinLabel(2, "Has_V0"); + hCandidate->GetXaxis()->SetBinLabel(2, "PassBuilderSel"); hCandidate->GetXaxis()->SetBinLabel(3, "DCA_meson"); hCandidate->GetXaxis()->SetBinLabel(4, "DCA_baryon"); hCandidate->GetXaxis()->SetBinLabel(5, "TPCNsigma_pion"); @@ -243,6 +243,7 @@ struct strangenessFilter { hCandidate->GetXaxis()->SetBinLabel(17, "CascCosPA"); hCandidate->GetXaxis()->SetBinLabel(18, "DCAV0ToPV"); hCandidate->GetXaxis()->SetBinLabel(19, "ProperLifeTime"); + hCandidate->GetXaxis()->SetBinLabel(20, "Rapidity"); std::vector centBinning = {0., 1., 5., 10., 20., 30., 40., 50., 70., 100.}; AxisSpec multAxisNTPV = {100, 0.0f, 100.0f, "N. tracks PV estimator"}; From 5c8df481e0d78d7534a61535aed694b38a5d5e41 Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Thu, 17 Apr 2025 16:46:26 +0000 Subject: [PATCH 5/5] Please consider the following formatting changes --- EventFiltering/PWGLF/strangenessFilter.cxx | 30 +++++++++++----------- 1 file changed, 15 insertions(+), 15 deletions(-) diff --git a/EventFiltering/PWGLF/strangenessFilter.cxx b/EventFiltering/PWGLF/strangenessFilter.cxx index 319f1d5ca35..ef311a4a221 100644 --- a/EventFiltering/PWGLF/strangenessFilter.cxx +++ b/EventFiltering/PWGLF/strangenessFilter.cxx @@ -669,8 +669,8 @@ struct strangenessFilter { } hCandidate->Fill(5.5); } - hCandidate->Fill(6.5); //OLD: eta dau (selection now applied in strangeness helper) - hCandidate->Fill(7.5); //OLD: bachtopv (selection now applied in strangeness helper) + hCandidate->Fill(6.5); // OLD: eta dau (selection now applied in strangeness helper) + hCandidate->Fill(7.5); // OLD: bachtopv (selection now applied in strangeness helper) // not striclty needed as selection are applied beforehand - just as QA (no change in number expected) if (Cascv0radius < v0radius) { @@ -680,7 +680,7 @@ struct strangenessFilter { if (Casccascradius < cascradius) { continue; } - hCandidate->Fill(9.5); + hCandidate->Fill(9.5); if (v0DauCPA < v0cospa) { continue; } @@ -690,11 +690,11 @@ struct strangenessFilter { } hCandidate->Fill(11.5); if (mStraHelper.cascade.cascadeDaughterDCA > dcacascdau) { - continue; + continue; } hCandidate->Fill(12.5); if (std::fabs(LambdaMass - constants::physics::MassLambda) > masslambdalimit) { - continue; + continue; } hCandidate->Fill(13.5); if (std::fabs(etaCasc) > eta) { @@ -711,16 +711,16 @@ struct strangenessFilter { // Fill selections QA for Xi if (cascCPA > casccospaxi) { - hCandidate->Fill(16.5); - if (cascCPA > dcav0topv) { - hCandidate->Fill(17.5); - if (xiproperlifetime < properlifetimefactor * ctauxi) { - hCandidate->Fill(18.5); - if (std::fabs(yXi) < rapidity) { - hCandidate->Fill(19.5); - } - } - } + hCandidate->Fill(16.5); + if (cascCPA > dcav0topv) { + hCandidate->Fill(17.5); + if (xiproperlifetime < properlifetimefactor * ctauxi) { + hCandidate->Fill(18.5); + if (std::fabs(yXi) < rapidity) { + hCandidate->Fill(19.5); + } + } + } } const auto deltaMassXi = useSigmaBasedMassCutXi ? getMassWindow(stfilter::species::Xi, ptCasc) : ximasswindow;