From c7363656ae31dad5378823b5d39c966f2c5bb6c5 Mon Sep 17 00:00:00 2001 From: FDUEnrich <23210190047@m.fudan.edu.cn> Date: Thu, 13 Feb 2025 14:21:07 +0800 Subject: [PATCH 1/4] Update treeCreatorOmegacToOmegaPi.cxx --- PWGHF/TableProducer/treeCreatorOmegacToOmegaPi.cxx | 9 ++++++++- 1 file changed, 8 insertions(+), 1 deletion(-) diff --git a/PWGHF/TableProducer/treeCreatorOmegacToOmegaPi.cxx b/PWGHF/TableProducer/treeCreatorOmegacToOmegaPi.cxx index 11f7d20a63e..1706ce6075c 100644 --- a/PWGHF/TableProducer/treeCreatorOmegacToOmegaPi.cxx +++ b/PWGHF/TableProducer/treeCreatorOmegacToOmegaPi.cxx @@ -207,6 +207,7 @@ struct HfTreeCreatorOmegac0ToOmegaPi { Produces rowEv; Configurable zPvCut{"zPvCut", 10., "Cut on absolute value of primary vertex z coordinate"}; + Configurable keepOnlyMcSignal{"keepOnlyMcSignal", true, "Fill MC tree only with signal candidates"}; using MyTrackTable = soa::Join; using MyEventTable = soa::Join; @@ -425,7 +426,13 @@ struct HfTreeCreatorOmegac0ToOmegaPi { // Filling candidate properties rowCandidateLite.reserve(candidates.size()); for (const auto& candidate : candidates) { - fillKfCandidate(candidate, candidate.flagMcMatchRec(), candidate.originRec(), candidate.collisionMatched()); + if (keepOnlyMcSignal){ + if (candidate.originRec() != 0) { + fillKfCandidate(candidate, candidate.flagMcMatchRec(), candidate.originRec(), candidate.collisionMatched()); + } + } else { + fillKfCandidate(candidate, candidate.flagMcMatchRec(), candidate.originRec(), candidate.collisionMatched()); + } } } PROCESS_SWITCH(HfTreeCreatorOmegac0ToOmegaPi, processKFMcFull, "Process KF MC", false); From 708a7eafd3f51f0095806c43459fba23e48e90a9 Mon Sep 17 00:00:00 2001 From: FDUEnrich <23210190047@m.fudan.edu.cn> Date: Thu, 13 Feb 2025 14:21:49 +0800 Subject: [PATCH 2/4] Update candidateSelectorOmegac0ToOmegaPi.cxx --- .../candidateSelectorOmegac0ToOmegaPi.cxx | 11 +++++++++++ 1 file changed, 11 insertions(+) diff --git a/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx b/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx index 97e98b26d87..ee351abd1d7 100644 --- a/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx +++ b/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx @@ -156,6 +156,7 @@ struct HfCandidateSelectorToOmegaPi { HistogramRegistry registry{"registry"}; // for QA of selections OutputObj hInvMassCharmBaryon{TH1D("hInvMassCharmBaryon", "Charm baryon invariant mass;inv mass;entries", 500, 2.3, 3.1)}; + OutputObj PtCharmBaryon{TH1D("PtCharmBaryon", "Charm baryon transverse momentum before sel;Pt;entries", 8000, 0., 80)}; void init(InitContext const&) { @@ -192,6 +193,7 @@ struct HfCandidateSelectorToOmegaPi { registry.add("hStatusCheck", "Check consecutive selections status;status;entries", {HistType::kTH1D, {{12, 0., 12.}}}); // for QA of the selections (bin 0 -> candidates that did not pass the selection, bin 1 -> candidates that passed the selection) + registry.add("hSelPtOmegac", "hSelPtOmegac;status;entries", {HistType::kTH1D, {axisSel}}); registry.add("hSelSignDec", "hSelSignDec;status;entries", {HistType::kTH1D, {axisSel}}); registry.add("hSelEtaPosV0Dau", "hSelEtaPosV0Dau;status;entries", {HistType::kTH1D, {axisSel}}); registry.add("hSelEtaNegV0Dau", "hSelEtaNegV0Dau;status;entries", {HistType::kTH1D, {axisSel}}); @@ -425,6 +427,15 @@ struct HfCandidateSelectorToOmegaPi { } } + //Omegac Pt selection + PtCharmBaryon->Fill(std::abs(candidate.kfptOmegac())); + if (std::abs(candidate.kfptOmegac()) < ptCandMin || std::abs(candidate.kfptOmegac()) > ptCandMax) { + resultSelections = false; + registry.fill(HIST("hSelPtOmegac"), 0); + } else { + registry.fill(HIST("hSelPtOmegac"), 1); + } + // v0&Casc&Omegac ldl selection if ((candidate.v0ldl() < v0LdlMin) || (candidate.cascldl() < cascLdlMin) || (candidate.omegacldl() > omegacLdlMax)) { resultSelections = false; From 65bbf6093f0a984056f3c4ccfd773f769617efb3 Mon Sep 17 00:00:00 2001 From: FDUEnrich <23210190047@m.fudan.edu.cn> Date: Thu, 13 Feb 2025 14:24:49 +0800 Subject: [PATCH 3/4] Update candidateSelectorOmegac0ToOmegaPi.cxx --- PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx b/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx index ee351abd1d7..757e58d514d 100644 --- a/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx +++ b/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx @@ -156,7 +156,7 @@ struct HfCandidateSelectorToOmegaPi { HistogramRegistry registry{"registry"}; // for QA of selections OutputObj hInvMassCharmBaryon{TH1D("hInvMassCharmBaryon", "Charm baryon invariant mass;inv mass;entries", 500, 2.3, 3.1)}; - OutputObj PtCharmBaryon{TH1D("PtCharmBaryon", "Charm baryon transverse momentum before sel;Pt;entries", 8000, 0., 80)}; + OutputObj hPtCharmBaryon{TH1D("hPtCharmBaryon", "Charm baryon transverse momentum before sel;Pt;entries", 8000, 0., 80)}; void init(InitContext const&) { @@ -428,7 +428,7 @@ struct HfCandidateSelectorToOmegaPi { } //Omegac Pt selection - PtCharmBaryon->Fill(std::abs(candidate.kfptOmegac())); + hPtCharmBaryon->Fill(std::abs(candidate.kfptOmegac())); if (std::abs(candidate.kfptOmegac()) < ptCandMin || std::abs(candidate.kfptOmegac()) > ptCandMax) { resultSelections = false; registry.fill(HIST("hSelPtOmegac"), 0); From 6552b4208694cbd1c0a7ea45beef93bf4f6c910d Mon Sep 17 00:00:00 2001 From: ALICE Action Bot Date: Thu, 13 Feb 2025 06:56:17 +0000 Subject: [PATCH 4/4] Please consider the following formatting changes --- PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx | 2 +- PWGHF/TableProducer/treeCreatorOmegacToOmegaPi.cxx | 4 ++-- 2 files changed, 3 insertions(+), 3 deletions(-) diff --git a/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx b/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx index 757e58d514d..117dee5430c 100644 --- a/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx +++ b/PWGHF/TableProducer/candidateSelectorOmegac0ToOmegaPi.cxx @@ -427,7 +427,7 @@ struct HfCandidateSelectorToOmegaPi { } } - //Omegac Pt selection + // Omegac Pt selection hPtCharmBaryon->Fill(std::abs(candidate.kfptOmegac())); if (std::abs(candidate.kfptOmegac()) < ptCandMin || std::abs(candidate.kfptOmegac()) > ptCandMax) { resultSelections = false; diff --git a/PWGHF/TableProducer/treeCreatorOmegacToOmegaPi.cxx b/PWGHF/TableProducer/treeCreatorOmegacToOmegaPi.cxx index 1706ce6075c..0b0c3a2a4e4 100644 --- a/PWGHF/TableProducer/treeCreatorOmegacToOmegaPi.cxx +++ b/PWGHF/TableProducer/treeCreatorOmegacToOmegaPi.cxx @@ -426,9 +426,9 @@ struct HfTreeCreatorOmegac0ToOmegaPi { // Filling candidate properties rowCandidateLite.reserve(candidates.size()); for (const auto& candidate : candidates) { - if (keepOnlyMcSignal){ + if (keepOnlyMcSignal) { if (candidate.originRec() != 0) { - fillKfCandidate(candidate, candidate.flagMcMatchRec(), candidate.originRec(), candidate.collisionMatched()); + fillKfCandidate(candidate, candidate.flagMcMatchRec(), candidate.originRec(), candidate.collisionMatched()); } } else { fillKfCandidate(candidate, candidate.flagMcMatchRec(), candidate.originRec(), candidate.collisionMatched());