From 94c8475e8ec8544a85ad08c812f70d83862cbd3e Mon Sep 17 00:00:00 2001 From: rspijkers Date: Thu, 3 Jun 2021 22:10:05 +0200 Subject: [PATCH 1/2] add Ct, Y, and chi2PCA to jpsi task --- Analysis/Tasks/PWGHF/taskJpsi.cxx | 14 +++++++++++++- 1 file changed, 13 insertions(+), 1 deletion(-) diff --git a/Analysis/Tasks/PWGHF/taskJpsi.cxx b/Analysis/Tasks/PWGHF/taskJpsi.cxx index 5600c5600899d..ae92ee37e503f 100644 --- a/Analysis/Tasks/PWGHF/taskJpsi.cxx +++ b/Analysis/Tasks/PWGHF/taskJpsi.cxx @@ -114,7 +114,13 @@ struct TaskJpsiMC { { registry.add("hmassSig", "2-prong candidates (rec matched);inv. mass (e^{#plus} e^{#minus}) (GeV/#it{c}^{2});entries", {HistType::kTH2F, {{500, 0., 5.}, {(std::vector)bins, "#it{p}_{T} (GeV/#it{c})"}}}); registry.add("hmassBkg", "2-prong candidates (rec unmatched);inv. mass (e^{#plus} e^{#minus}) (GeV/#it{c}^{2});entries", {HistType::kTH2F, {{500, 0., 5.}, {(std::vector)bins, "#it{p}_{T} (GeV/#it{c})"}}}); - } + registry.add("hChi2PCASig", "3-prong candidates (rec. matched);chi2 PCA (cm);entries", {HistType::kTH2F, {{500, 0., 0.01}, {(std::vector)bins, "#it{p}_{T} (GeV/#it{c})"}}}); + registry.add("hChi2PCABg", "3-prong candidates (rec. unmatched);chi2 PCA (cm);entries", {HistType::kTH2F, {{500, 0., 0.01}, {(std::vector)bins, "#it{p}_{T} (GeV/#it{c})"}}}); + registry.add("hCtSig", "3-prong candidates (rec. matched);proper lifetime X(3872) * #it{c} (cm);entries", {HistType::kTH2F, {{100, 0., 0.1}, {(std::vector)bins, "#it{p}_{T} (GeV/#it{c})"}}}); + registry.add("hCtBg", "3-prong candidates (rec. unmatched);proper lifetime X(3872) * #it{c} (cm);entries", {HistType::kTH2F, {{100, 0., 0.1}, {(std::vector)bins, "#it{p}_{T} (GeV/#it{c})"}}}); + registry.add("hYSig", "3-prong candidates (rec. matched);candidate rapidity;entries", {HistType::kTH2F, {{100, -2., 2.}, {(std::vector)bins, "#it{p}_{T} (GeV/#it{c})"}}}); + registry.add("hYBg", "3-prong candidates (rec. unmatched);candidate rapidity;entries", {HistType::kTH2F, {{100, -2., 2.}, {(std::vector)bins, "#it{p}_{T} (GeV/#it{c})"}}}); +} Filter filterSelectCandidates = (aod::hf_selcandidate_jpsi::isSelJpsiToEE >= d_selectionFlagJpsi); @@ -139,11 +145,17 @@ struct TaskJpsiMC { registry.fill(HIST("hCPARecSig"), candidate.cpa()); registry.fill(HIST("hEtaRecSig"), candidate.eta()); registry.fill(HIST("hmassSig"), InvMassJpsiToEE(candidate), candidate.pt()); + registry.fill(HIST("hChi2PCASig"), candidate.chi2PCA(), candidate.pt()); + registry.fill(HIST("hCtSig"), CtJpsi(candidate), candidate.pt()); + registry.fill(HIST("hYSig"), YJpsi(candidate), candidate.pt()); } else { registry.fill(HIST("hPtRecBg"), candidate.pt()); registry.fill(HIST("hCPARecBg"), candidate.cpa()); registry.fill(HIST("hEtaRecBg"), candidate.eta()); registry.fill(HIST("hmassBkg"), InvMassJpsiToEE(candidate), candidate.pt()); + registry.fill(HIST("hChi2PCABg"), candidate.chi2PCA(), candidate.pt()); + registry.fill(HIST("hCtBg"), CtJpsi(candidate), candidate.pt()); + registry.fill(HIST("hYBg"), YJpsi(candidate), candidate.pt()); } } // MC gen. From fd7b37c38187deff2db765deaa54e5b42517395e Mon Sep 17 00:00:00 2001 From: rspijkers Date: Thu, 3 Jun 2021 22:25:41 +0200 Subject: [PATCH 2/2] fix clang --- Analysis/Tasks/PWGHF/taskJpsi.cxx | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/Analysis/Tasks/PWGHF/taskJpsi.cxx b/Analysis/Tasks/PWGHF/taskJpsi.cxx index ae92ee37e503f..e50a9151fd77f 100644 --- a/Analysis/Tasks/PWGHF/taskJpsi.cxx +++ b/Analysis/Tasks/PWGHF/taskJpsi.cxx @@ -120,7 +120,7 @@ struct TaskJpsiMC { registry.add("hCtBg", "3-prong candidates (rec. unmatched);proper lifetime X(3872) * #it{c} (cm);entries", {HistType::kTH2F, {{100, 0., 0.1}, {(std::vector)bins, "#it{p}_{T} (GeV/#it{c})"}}}); registry.add("hYSig", "3-prong candidates (rec. matched);candidate rapidity;entries", {HistType::kTH2F, {{100, -2., 2.}, {(std::vector)bins, "#it{p}_{T} (GeV/#it{c})"}}}); registry.add("hYBg", "3-prong candidates (rec. unmatched);candidate rapidity;entries", {HistType::kTH2F, {{100, -2., 2.}, {(std::vector)bins, "#it{p}_{T} (GeV/#it{c})"}}}); -} + } Filter filterSelectCandidates = (aod::hf_selcandidate_jpsi::isSelJpsiToEE >= d_selectionFlagJpsi);