Skip to content

Commit 26d8ca6

Browse files
authored
[PWGCF] FU D0 task: adding y histogram and non-prompt cut for reco D0s (#18074)
1 parent cb2bcbc commit 26d8ca6

1 file changed

Lines changed: 3 additions & 1 deletion

File tree

‎PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackD0.cxx‎

Lines changed: 3 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -314,6 +314,7 @@ struct FemtoUniversePairTaskTrackD0 {
314314
{"hEtaD0D0bar", ";#eta ;counts", {HistType::kTH1F, {{200, -1., 1.}}}},
315315
{"hEtaD0", ";#eta ;counts", {HistType::kTH1F, {{200, -1., 1.}}}},
316316
{"hEtaD0bar", ";#eta ;counts", {HistType::kTH1F, {{200, -1., 1.}}}},
317+
{"hYD0D0bar", ";y ;counts", {HistType::kTH1F, {{200, -1., 1.}}}},
317318
{"hPtDaughters", ";#it{p}_{T} (GeV/#it{c});counts", {HistType::kTH1F, {{300, 0., 12.}}}},
318319
{"hSignDaughters", ";sign ;counts", {HistType::kTH1F, {{10, -2.5, 2.5}}}}}};
319320

@@ -675,6 +676,7 @@ struct FemtoUniversePairTaskTrackD0 {
675676
registry.fill(HIST("hPtD0D0bar"), d0d0bar.pt());
676677
registry.fill(HIST("hPhiD0D0bar"), d0d0bar.phi());
677678
registry.fill(HIST("hEtaD0D0bar"), d0d0bar.eta());
679+
registry.fill(HIST("hYD0D0bar"), d0d0bar.transRadius());
678680
// BDT score classes
679681
registry.fill(HIST("DebugBdt/hBdtScore0"), d0d0bar.decayVtxX());
680682
registry.fill(HIST("DebugBdt/hBdtScore1"), d0d0bar.decayVtxY());
@@ -1388,7 +1390,7 @@ struct FemtoUniversePairTaskTrackD0 {
13881390
{
13891391
for (auto const& part : recoParts) {
13901392
// filling the histograms for identified hadrons
1391-
if ((part.partType() == aod::femtouniverseparticle::ParticleType::kD0) && (part.pt() > ConfDmesons.confMinPtD0D0barReco) && (part.pt() < ConfDmesons.confMaxPtD0D0barReco)) {
1393+
if ((part.partType() == aod::femtouniverseparticle::ParticleType::kD0) && (part.pt() > ConfDmesons.confMinPtD0D0barReco) && (part.pt() < ConfDmesons.confMaxPtD0D0barReco) && (part.decayVtxZ() < ConfMlProb.confMlProbNonPromptMax)) {
13921394
// getting the efficiency value
13931395
if (doEfficiencyCorr) {
13941396
weight = efficiencyCalculator.getWeight(ParticleNo::TWO, part.pt());

0 commit comments

Comments
 (0)