Skip to content

Commit e0481d5

Browse files
ercolessiercolessi
authored
[PWGLF] V0 analysis: add feeddown for Lambdas (#12971)
1 parent ef83bb6 commit e0481d5

File tree

1 file changed

+25
-26
lines changed

1 file changed

+25
-26
lines changed

‎PWGLF/Tasks/Strangeness/v0postprocessing.cxx‎

Lines changed: 25 additions & 26 deletions
Original file line numberDiff line numberDiff line change
@@ -13,14 +13,16 @@
1313
/// \author Francesca Ercolessi (francesca.ercolessi@cern.ch)
1414
/// \since
1515

16-
#include "Framework/runDataProcessing.h"
17-
#include "Framework/AnalysisTask.h"
18-
#include "Common/DataModel/TrackSelectionTables.h"
1916
#include "PWGLF/DataModel/LFStrangenessTables.h"
2017
#include "PWGLF/DataModel/v0qaanalysis.h"
18+
2119
#include "Common/DataModel/EventSelection.h"
2220
#include "Common/DataModel/PIDResponse.h"
21+
#include "Common/DataModel/TrackSelectionTables.h"
22+
2323
#include "CommonConstants/PhysicsConstants.h"
24+
#include "Framework/AnalysisTask.h"
25+
#include "Framework/runDataProcessing.h"
2426

2527
using namespace o2;
2628
using namespace o2::framework;
@@ -122,11 +124,11 @@ struct v0postprocessing {
122124
registry.add("hMassLambda_MC", "hMassLambda", {HistType::kTH1F, {{200, 1.016f, 1.216f}}});
123125
registry.add("hMassVsPtLambdaVsCentFT0M_MC", ";p_{T} [GeV/c];M_{p^{+}#pi^{-}} [GeV/c^{2}]", {HistType::kTH3F, {{250, 0.0f, 25.0f}, {100, 0.f, 100.f}, {200, 1.016f, 1.216f}}});
124126
registry.add("hMassAntiLambda_MC", "hMassAntiLambda", {HistType::kTH1F, {{200, 1.016f, 1.216f}}});
125-
registry.add("hMassVsPtLambdaVsMotherPt_DoubleCharged_MC", ";p_{T} [GeV/c] (V0);p_{T}^{gen} [GeV/c] (Xi);M_{p^{-}#pi^{+}} [GeV/c^{2}]", {HistType::kTH3F, {{250, 0.0f, 25.0f}, {250, 0.0f, 25.0f}, {200, 1.016f, 1.216f}}});
126-
registry.add("hMassVsPtLambdaVsMotherPt_MCRatio_MC", ";p_{T} [GeV/c] (V0);p_{T}^{gen} [GeV/c] (Xi);M_{p^{-}#pi^{+}} [GeV/c^{2}]", {HistType::kTH3F, {{250, 0.0f, 25.0f}, {250, 0.0f, 25.0f}, {200, 1.016f, 1.216f}}});
127+
registry.add("hFDVsPtLambdaVsMotherPt_DoubleCharged_MC", ";p_{T} [GeV/c] (V0);p_{T}^{gen} [GeV/c] (#Xi^{-}); percentile", {HistType::kTH3F, {{250, 0.0f, 25.0f}, {250, 0.0f, 25.0f}, {100, 0.f, 100.f}}});
128+
registry.add("hFDVsPtLambdaVsMotherPt_MCRatio_MC", ";p_{T} [GeV/c] (V0);p_{T}^{gen} [GeV/c] (#Xi^{-/0}); percentile", {HistType::kTH3F, {{250, 0.0f, 25.0f}, {250, 0.0f, 25.0f}, {100, 0.f, 100.f}}});
127129
registry.add("hMassVsPtAntiLambdaVsCentFT0M_MC", ";p_{T} [GeV/c];M_{p^{-}#pi^{+}} [GeV/c^{2}]", {HistType::kTH3F, {{250, 0.0f, 25.0f}, {100, 0.f, 100.f}, {200, 1.016f, 1.216f}}});
128-
registry.add("hMassVsPtAntiLambdaVsMotherPt_DoubleCharged_MC", ";p_{T} [GeV/c] (V0);p_{T}^{gen} [GeV/c] (Xi);M_{p^{-}#pi^{+}} [GeV/c^{2}]", {HistType::kTH3F, {{250, 0.0f, 25.0f}, {250, 0.0f, 25.0f}, {200, 1.016f, 1.216f}}});
129-
registry.add("hMassVsPtAntiLambdaVsMotherPt_MCRatio_MC", ";p_{T} [GeV/c] (V0);p_{T}^{gen} [GeV/c] (Xi);M_{p^{-}#pi^{+}} [GeV/c^{2}]", {HistType::kTH3F, {{250, 0.0f, 25.0f}, {250, 0.0f, 25.0f}, {200, 1.016f, 1.216f}}});
130+
registry.add("hFDVsPtAntiLambdaVsMotherPt_DoubleCharged_MC", ";p_{T} [GeV/c] (V0);p_{T}^{gen} [GeV/c] (#bar{#Xi}^{+});percentile", {HistType::kTH3F, {{250, 0.0f, 25.0f}, {250, 0.0f, 25.0f}, {100, 0.f, 100.f}}});
131+
registry.add("hFDVsPtAntiLambdaVsMotherPt_MCRatio_MC", ";p_{T} [GeV/c] (V0);p_{T}^{gen} [GeV/c] (#bar{#Xi}^{+/0});percentile", {HistType::kTH3F, {{250, 0.0f, 25.0f}, {250, 0.0f, 25.0f}, {100, 0.f, 100.f}}});
130132
}
131133

132134
if (doQA) {
@@ -398,7 +400,7 @@ struct v0postprocessing {
398400
if (candidate.v0cospa() > cospaK0s &&
399401
std::abs(candidate.rapk0short()) < rap &&
400402
candidate.ctauk0short() < ctauK0s &&
401-
std::abs(candidate.massk0short() - o2::constants::physics::MassK0Short) < 0.075 &&
403+
std::abs(candidate.massk0short() - o2::constants::physics::MassK0Short) < 0.1 &&
402404
std::abs(candidate.masslambda() - o2::constants::physics::MassLambda0) > v0rejK0s &&
403405
std::abs(candidate.ntpcsigmanegpi()) <= ntpcsigma &&
404406
std::abs(candidate.ntpcsigmapospi()) <= ntpcsigma &&
@@ -448,20 +450,17 @@ struct v0postprocessing {
448450
registry.fill(HIST("hMassVsPtLambda"), candidate.v0pt(), candidate.masslambda());
449451
registry.fill(HIST("hMassVsPtLambdaVsCentFT0M"), candidate.v0pt(), candidate.multft0m(), candidate.masslambda());
450452

451-
if (isMC) {
452-
453-
if (candidate.pdgcode() == 3122 && candidate.isdaulambda()) {
454-
455-
if (candidate.isphysprimary() == 1) {
456-
registry.fill(HIST("hMassLambda_MC"), candidate.masslambda());
457-
registry.fill(HIST("hMassVsPtLambdaVsCentFT0M_MC"), candidate.v0pt(), candidate.multft0m(), candidate.masslambda());
458-
}
453+
if (isMC && candidate.pdgcode() == 3122 && candidate.isdaulambda()) {
459454

455+
if (candidate.isphysprimary() == 1) {
456+
registry.fill(HIST("hMassLambda_MC"), candidate.masslambda());
457+
registry.fill(HIST("hMassVsPtLambdaVsCentFT0M_MC"), candidate.v0pt(), candidate.multft0m(), candidate.masslambda());
458+
} else if (std::abs(candidate.masslambda() - o2::constants::physics::MassLambda0) < 0.01) {
460459
if (candidate.pdgcodemother() == 3312) {
461-
registry.fill(HIST("hMassVsPtLambdaVsMotherPt_DoubleCharged_MC"), candidate.v0pt(), candidate.v0motherpt(), candidate.masslambda());
460+
registry.fill(HIST("hFDVsPtLambdaVsMotherPt_DoubleCharged_MC"), candidate.v0pt(), candidate.v0motherpt(), candidate.multft0m());
462461
}
463462
if (candidate.pdgcodemother() == 3312 || candidate.pdgcodemother() == 3322) {
464-
registry.fill(HIST("hMassVsPtLambdaVsMotherPt_MCRatio_MC"), candidate.v0pt(), candidate.v0motherpt(), candidate.masslambda());
463+
registry.fill(HIST("hFDVsPtLambdaVsMotherPt_MCRatio_MC"), candidate.v0pt(), candidate.v0motherpt(), candidate.multft0m());
465464
}
466465
}
467466
}
@@ -489,18 +488,18 @@ struct v0postprocessing {
489488
registry.fill(HIST("hMassVsPtAntiLambda"), candidate.v0pt(), candidate.massantilambda());
490489
registry.fill(HIST("hMassVsPtAntiLambdaVsCentFT0M"), candidate.v0pt(), candidate.multft0m(), candidate.massantilambda());
491490

492-
if (candidate.pdgcode() == -3122 && candidate.isdauantilambda()) {
491+
if (isMC && candidate.pdgcode() == -3122 && candidate.isdauantilambda()) {
493492

494493
if (candidate.isphysprimary() == 1) {
495494
registry.fill(HIST("hMassAntiLambda_MC"), candidate.massantilambda());
496495
registry.fill(HIST("hMassVsPtAntiLambdaVsCentFT0M_MC"), candidate.v0pt(), candidate.v0motherpt(), candidate.massantilambda());
497-
}
498-
499-
if (candidate.pdgcodemother() == -3312) {
500-
registry.fill(HIST("hMassVsPtAntiLambdaVsMotherPt_DoubleCharged_MC"), candidate.v0pt(), candidate.v0motherpt(), candidate.massantilambda());
501-
}
502-
if (candidate.pdgcodemother() == -3312 || candidate.pdgcodemother() == -3322) {
503-
registry.fill(HIST("hMassVsPtAntiLambdaVsMotherPt_MCRatio_MC"), candidate.v0pt(), candidate.v0motherpt(), candidate.massantilambda());
496+
} else if (std::abs(candidate.massantilambda() - o2::constants::physics::MassLambda0) < 0.01) {
497+
if (candidate.pdgcodemother() == -3312) {
498+
registry.fill(HIST("hFDVsPtAntiLambdaVsMotherPt_DoubleCharged_MC"), candidate.v0pt(), candidate.v0motherpt(), candidate.multft0m());
499+
}
500+
if (candidate.pdgcodemother() == -3312 || candidate.pdgcodemother() == -3322) {
501+
registry.fill(HIST("hFDVsPtAntiLambdaVsMotherPt_MCRatio_MC"), candidate.v0pt(), candidate.v0motherpt(), candidate.multft0m());
502+
}
504503
}
505504
}
506505

0 commit comments

Comments
 (0)