Skip to content

Commit f877b4b

Browse files
alicja-ppalicja-pp
authored
[PWGCF] FemtoUniverse: add invariant mass vs pT and centrality histograms (#10259)
1 parent d0bb0b8 commit f877b4b

File tree

2 files changed

+21
-9
lines changed

2 files changed

+21
-9
lines changed

‎PWGCF/FemtoUniverse/Core/FemtoUniverseParticleHisto.h‎

Lines changed: 13 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -107,7 +107,9 @@ class FemtoUniverseParticleHisto
107107
mHistogramRegistry->add((folderName + folderSuffix + "/hDecayVtxY").c_str(), "; #it{Vtx}_{y} (cm)); Entries", kTH1F, {{2000, 0, 200}});
108108
mHistogramRegistry->add((folderName + folderSuffix + "/hDecayVtxZ").c_str(), "; #it{Vtx}_{z} (cm); Entries", kTH1F, {{2000, 0, 200}});
109109
mHistogramRegistry->add((folderName + folderSuffix + "/hInvMassLambda").c_str(), "; M_{#Lambda}; Entries", kTH1F, {{2000, 1.f, 3.f}});
110+
mHistogramRegistry->add((folderName + folderSuffix + "/hInvMassLambdaVsPt").c_str(), "; #it{p}_{T} (GeV/#it{c}); M_{#Lambda}; Entries", kTH2F, {{240, 0, 6}, {2000, 1.f, 3.f}});
110111
mHistogramRegistry->add((folderName + folderSuffix + "/hInvMassAntiLambda").c_str(), "; M_{#bar{#Lambda}}; Entries", kTH1F, {{2000, 1.f, 3.f}});
112+
mHistogramRegistry->add((folderName + folderSuffix + "/hInvMassAntiLambdaVsPt").c_str(), "; #it{p}_{T} (GeV/#it{c}); M_{#Lambda}; Entries", kTH2F, {{240, 0, 6}, {2000, 1.f, 3.f}});
111113
mHistogramRegistry->add((folderName + folderSuffix + "/hInvMassLambdaAntiLambda").c_str(), "; M_{#Lambda}; M_{#bar{#Lambda}}", kTH2F, {{2000, 1.f, 3.f}, {2000, 1.f, 3.f}});
112114
} else if constexpr (mParticleType == o2::aod::femtouniverseparticle::ParticleType::kCascade) {
113115
mHistogramRegistry->add((folderName + folderSuffix + "/hDaughDCA").c_str(), "; DCA^{daugh} (cm); Entries", kTH1F, {{1000, 0, 10}});
@@ -282,7 +284,9 @@ class FemtoUniverseParticleHisto
282284
mHistogramRegistry->fill(histFolder + HIST(o2::aod::femtouniverse_mc_particle::MCTypeName[mc]) + HIST("/hDecayVtxY"), part.decayVtxY());
283285
mHistogramRegistry->fill(histFolder + HIST(o2::aod::femtouniverse_mc_particle::MCTypeName[mc]) + HIST("/hDecayVtxZ"), part.decayVtxZ());
284286
mHistogramRegistry->fill(histFolder + HIST(o2::aod::femtouniverse_mc_particle::MCTypeName[mc]) + HIST("/hInvMassLambda"), part.mLambda());
287+
mHistogramRegistry->fill(histFolder + HIST(o2::aod::femtouniverse_mc_particle::MCTypeName[mc]) + HIST("/hInvMassLambdaVsPt"), part.pt(), part.mLambda());
285288
mHistogramRegistry->fill(histFolder + HIST(o2::aod::femtouniverse_mc_particle::MCTypeName[mc]) + HIST("/hInvMassAntiLambda"), part.mAntiLambda());
289+
mHistogramRegistry->fill(histFolder + HIST(o2::aod::femtouniverse_mc_particle::MCTypeName[mc]) + HIST("/hInvMassAntiLambdaVsPt"), part.pt(), part.mAntiLambda());
286290
mHistogramRegistry->fill(histFolder + HIST(o2::aod::femtouniverse_mc_particle::MCTypeName[mc]) + HIST("/hInvMassLambdaAntiLambda"), part.mLambda(), part.mAntiLambda());
287291
} else if constexpr (mParticleType == o2::aod::femtouniverseparticle::ParticleType::kCascade) {
288292
mHistogramRegistry->fill(histFolder + HIST(o2::aod::femtouniverse_mc_particle::MCTypeName[mc]) + HIST("/hDaughDCA"), part.daughDCA());
@@ -367,26 +371,26 @@ class FemtoUniverseParticleHisto
367371
if (mHistogramRegistry) {
368372
if constexpr (mParticleType == o2::aod::femtouniverseparticle::ParticleType::kTrack) {
369373
if (confPDG == mConfPDGCodePart[0]) {
370-
PDGbin = 0;
374+
binPDG = 0;
371375
} else if (confPDG == mConfPDGCodePart[1]) {
372-
PDGbin = 1;
376+
binPDG = 1;
373377
} else if (confPDG == mConfPDGCodePart[2]) {
374-
PDGbin = 2;
378+
binPDG = 2; // o2-linter: disable=pdg/explicit-code
375379
} else {
376-
PDGbin = 3;
380+
binPDG = 3; // o2-linter: disable=pdg/explicit-code
377381
}
378382
if (std::abs(pdgcode) == 211) {
379383
mHistogramRegistry->fill(histFolder + HIST("_MC/hMisidentification"),
380-
PDGbin, 0, part.pt());
384+
binPDG, 0, part.pt());
381385
} else if (std::abs(pdgcode) == 321) {
382386
mHistogramRegistry->fill(histFolder + HIST("_MC/hMisidentification"),
383-
PDGbin, 1, part.pt());
387+
binPDG, 1, part.pt());
384388
} else if (std::abs(pdgcode) == 2212) {
385389
mHistogramRegistry->fill(histFolder + HIST("_MC/hMisidentification"),
386-
PDGbin, 2, part.pt());
390+
binPDG, 2, part.pt());
387391
} else {
388392
mHistogramRegistry->fill(histFolder + HIST("_MC/hMisidentification"),
389-
PDGbin, 3, part.pt());
393+
binPDG, 3, part.pt());
390394
}
391395
}
392396
} else {
@@ -458,7 +462,7 @@ class FemtoUniverseParticleHisto
458462
static constexpr std::string_view mFolderSuffix[5] = {"", "_one", "_two", "_pos", "_neg"}; ///< Suffix for the folder name in case of analyses of pairs of the same kind (T-T, V-V, C-C) // o2-linter: disable=name/constexpr-constant
459463
int mConfPDGCodePart[4] = {211, 321, 2212, 9999}; ///< PDG code as per analysis
460464
int mPDG = 0; ///< PDG code of the selected particle
461-
int PDGbin = 0;
465+
int binPDG = 0;
462466
};
463467
} // namespace o2::analysis::femto_universe
464468

‎PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackV0Extended.cxx‎

Lines changed: 8 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -219,6 +219,10 @@ struct FemtoUniversePairTaskTrackV0Extended {
219219
posChildHistos.init(&qaRegistry, confChildTempFitVarpTBins, confChildTempFitVarBins, false, 0, true);
220220
negChildHistos.init(&qaRegistry, confChildTempFitVarpTBins, confChildTempFitVarBins, false, 0, true);
221221

222+
qaRegistry.add("V0Type1/hInvMassLambdaVsCent", "; Centrality; M_{#Lambda}; Entries", kTH2F, {confMultBins, {2000, 1.f, 3.f}});
223+
qaRegistry.add("V0Type2/hInvMassLambdaVsCent", "; Centrality; M_{#Lambda}; Entries", kTH2F, {confMultBins, {2000, 1.f, 3.f}});
224+
qaRegistry.add("V0Type1/hInvMassAntiLambdaVsCent", "; Centrality; M_{#Lambda}; Entries", kTH2F, {confMultBins, {2000, 1.f, 3.f}});
225+
qaRegistry.add("V0Type2/hInvMassAntiLambdaVsCent", "; Centrality; M_{#Lambda}; Entries", kTH2F, {confMultBins, {2000, 1.f, 3.f}});
222226
trackHistoV0Type1.init(&qaRegistry, confV0TempFitVarpTBins, confV0TempFitVarBins, confIsMC, confV0PDGCodePartTwo, true, "V0Type1");
223227
posChildV0Type1.init(&qaRegistry, confChildTempFitVarpTBins, confChildTempFitVarBins, false, 0, true, "posChildV0Type1");
224228
negChildV0Type1.init(&qaRegistry, confChildTempFitVarpTBins, confChildTempFitVarBins, false, 0, true, "negChildV0Type1");
@@ -410,12 +414,16 @@ struct FemtoUniversePairTaskTrackV0Extended {
410414
trackHistoV0Type1.fillQABase<false, true>(part, HIST("V0Type1"));
411415
posChildV0Type1.fillQABase<false, true>(posChild, HIST("posChildV0Type1"));
412416
negChildV0Type1.fillQABase<false, true>(negChild, HIST("negChildV0Type1"));
417+
qaRegistry.fill(HIST("V0Type1/hInvMassLambdaVsCent"), multCol, part.mLambda());
418+
qaRegistry.fill(HIST("V0Type1/hInvMassAntiLambdaVsCent"), multCol, part.mAntiLambda());
413419
}
414420
/// Check daughters of second V0 particle
415421
if (isParticleTPC(posChild, V0ChildTable[confV0Type2][0]) && isParticleTPC(negChild, V0ChildTable[confV0Type2][1])) {
416422
trackHistoV0Type2.fillQABase<false, true>(part, HIST("V0Type2"));
417423
posChildV0Type2.fillQABase<false, true>(posChild, HIST("posChildV0Type2"));
418424
negChildV0Type2.fillQABase<false, true>(negChild, HIST("negChildV0Type2"));
425+
qaRegistry.fill(HIST("V0Type2/hInvMassLambdaVsCent"), multCol, part.mLambda());
426+
qaRegistry.fill(HIST("V0Type2/hInvMassAntiLambdaVsCent"), multCol, part.mAntiLambda());
419427
}
420428
}
421429

0 commit comments

Comments
 (0)