From 219c60bd7bbb2fcc1d01cb2b00a1c122d40532bc Mon Sep 17 00:00:00 2001 From: Hirak Koley Date: Mon, 16 Mar 2026 14:09:31 +0530 Subject: [PATCH 1/2] Adjust MC label bins and update event selection logic --- .../Resonances/lambda1520analysisinpp.cxx | 22 +++++++++++++------ 1 file changed, 15 insertions(+), 7 deletions(-) diff --git a/PWGLF/Tasks/Resonances/lambda1520analysisinpp.cxx b/PWGLF/Tasks/Resonances/lambda1520analysisinpp.cxx index a494b545539..8b94d7b56c6 100644 --- a/PWGLF/Tasks/Resonances/lambda1520analysisinpp.cxx +++ b/PWGLF/Tasks/Resonances/lambda1520analysisinpp.cxx @@ -262,7 +262,7 @@ struct Lambda1520analysisinpp { AxisSpec axisTPCXrow{binsTPCXrows, "#Xrows_{TPC}"}; AxisSpec axisPIDQA{binsnSigma, "#sigma"}; AxisSpec axisTPCSignal{binsnTPCSignal, ""}; - AxisSpec axisMClabel{6, -1.5f, 6.5f, "MC Label"}; + AxisSpec axisMClabel{9, -1.5f, 7.5f, "MC Label"}; AxisSpec axisEtaPhi{binsEtaPhi, ""}; AxisSpec axisPhi{350, 0, 7, "#Phi"}; AxisSpec axisMultMix{configBkg.cfgMultPercentileBins, "Multiplicity Percentile"}; @@ -1071,6 +1071,8 @@ struct Lambda1520analysisinpp { histos.fill(HIST("QAevent/hEventsMC"), 2); } + colCuts.fillQA(collision); + fillHistograms(collision, tracks, tracks); } PROCESS_SWITCH(Lambda1520analysisinpp, processMCRec, "Process Event for MC Rec without partition", false); @@ -1136,20 +1138,19 @@ struct Lambda1520analysisinpp { histos.fill(HIST("QA/MC/h2GenEtaPt_afterEtaRapCut"), part.eta(), part.pt()); histos.fill(HIST("QA/MC/h2GenPhiRapidity_afterEtaRapCut"), part.phi(), part.y()); - // without any event selection - if (part.pdgCode() > 0) + if (part.pdgCode() > 0) // without any event selection histos.fill(HIST("Result/MC/Genlambda1520pt"), 0, part.pt(), centrality); else histos.fill(HIST("Result/MC/Genantilambda1520pt"), 0, part.pt(), centrality); - if (inVtx10) // INEL10 + if (inVtx10) // vtx10 { if (part.pdgCode() > 0) histos.fill(HIST("Result/MC/Genlambda1520pt"), 1, part.pt(), centrality); else histos.fill(HIST("Result/MC/Genantilambda1520pt"), 1, part.pt(), centrality); } - if (inVtx10 && isSel8) // INEL>10, vtx10 + if (inVtx10 && isSel8) // vtx10, sel8 { if (part.pdgCode() > 0) histos.fill(HIST("Result/MC/Genlambda1520pt"), 2, part.pt(), centrality); @@ -1163,20 +1164,27 @@ struct Lambda1520analysisinpp { else histos.fill(HIST("Result/MC/Genantilambda1520pt"), 3, part.pt(), centrality); } - if (isInAfterAllCuts) // after all event selection + if (inVtx10 && isTrueINELgt0) // vtx10, INEL>0 { if (part.pdgCode() > 0) histos.fill(HIST("Result/MC/Genlambda1520pt"), 4, part.pt(), centrality); else histos.fill(HIST("Result/MC/Genantilambda1520pt"), 4, part.pt(), centrality); } - if (isInAfterAllCuts && isTrueINELgt0) // after all event selection + if (isInAfterAllCuts) // after all event selection { if (part.pdgCode() > 0) histos.fill(HIST("Result/MC/Genlambda1520pt"), 5, part.pt(), centrality); else histos.fill(HIST("Result/MC/Genantilambda1520pt"), 5, part.pt(), centrality); } + if (isInAfterAllCuts && isTrueINELgt0) // after all event selection && INEL>0 + { + if (part.pdgCode() > 0) + histos.fill(HIST("Result/MC/Genlambda1520pt"), 6, part.pt(), centrality); + else + histos.fill(HIST("Result/MC/Genantilambda1520pt"), 6, part.pt(), centrality); + } } // QA for Trigger efficiency From e8e012686fa5df4fb69d5a09c01f79a8c6dd3c04 Mon Sep 17 00:00:00 2001 From: Hirak Koley Date: Mon, 16 Mar 2026 23:45:03 +0530 Subject: [PATCH 2/2] Adjust MC Label axis range and add Gen histogram --- PWGLF/Tasks/Resonances/lambda1520analysisinpp.cxx | 12 ++++++++++-- 1 file changed, 10 insertions(+), 2 deletions(-) diff --git a/PWGLF/Tasks/Resonances/lambda1520analysisinpp.cxx b/PWGLF/Tasks/Resonances/lambda1520analysisinpp.cxx index 8b94d7b56c6..000957d1fdf 100644 --- a/PWGLF/Tasks/Resonances/lambda1520analysisinpp.cxx +++ b/PWGLF/Tasks/Resonances/lambda1520analysisinpp.cxx @@ -262,7 +262,7 @@ struct Lambda1520analysisinpp { AxisSpec axisTPCXrow{binsTPCXrows, "#Xrows_{TPC}"}; AxisSpec axisPIDQA{binsnSigma, "#sigma"}; AxisSpec axisTPCSignal{binsnTPCSignal, ""}; - AxisSpec axisMClabel{9, -1.5f, 7.5f, "MC Label"}; + AxisSpec axisMClabel{9, -1.f, 8.f, "MC Label"}; AxisSpec axisEtaPhi{binsEtaPhi, ""}; AxisSpec axisPhi{350, 0, 7, "#Phi"}; AxisSpec axisMultMix{configBkg.cfgMultPercentileBins, "Multiplicity Percentile"}; @@ -276,7 +276,7 @@ struct Lambda1520analysisinpp { histos.add("QAevent/hPairsCounterSameE", "total valid no. of pairs sameE", HistType::kTH1F, {{1, 0.5f, 1.5f}}); histos.add("QAevent/hnTrksSameE", "n tracks per event SameE", HistType::kTH1F, {{1000, 0.0, 1000.0}}); } - // Test on Mixed event + // Gen on Mixed event if (doprocessME) { // Histograms for Mixed Event Pool characteristics @@ -401,6 +401,7 @@ struct Lambda1520analysisinpp { // MC QA histos.add("Event/hMCEventIndices", "hMCEventIndices", kTH2D, {axisMult, idxMCAxis}); if (doprocessMCGen) { + histos.add("QA/Gen", "Gen histogram", kTH1D, {{10, 0, 10, "index"}}); histos.add("QA/MC/h2GenEtaPt_beforeanycut", " #eta-#it{p}_{T} distribution of Generated #Lambda(1520); #eta; #it{p}_{T}; Counts;", HistType::kTHnSparseF, {axisEta, axisPtQA}); histos.add("QA/MC/h2GenPhiRapidity_beforeanycut", " #phi-y distribution of Generated #Lambda(1520); #phi; y; Counts;", HistType::kTHnSparseF, {axisPhi, axisRap}); histos.add("QA/MC/h2GenEtaPt_afterEtaRapCut", " #eta-#it{p}_{T} distribution of Generated #Lambda(1520); #eta; #it{p}_{T}; Counts;", HistType::kTHnSparseF, {axisEta, axisPtQA}); @@ -1138,6 +1139,7 @@ struct Lambda1520analysisinpp { histos.fill(HIST("QA/MC/h2GenEtaPt_afterEtaRapCut"), part.eta(), part.pt()); histos.fill(HIST("QA/MC/h2GenPhiRapidity_afterEtaRapCut"), part.phi(), part.y()); + histos.fill(HIST("QA/Gen"), 1); if (part.pdgCode() > 0) // without any event selection histos.fill(HIST("Result/MC/Genlambda1520pt"), 0, part.pt(), centrality); else @@ -1145,6 +1147,7 @@ struct Lambda1520analysisinpp { if (inVtx10) // vtx10 { + histos.fill(HIST("QA/Gen"), 2); if (part.pdgCode() > 0) histos.fill(HIST("Result/MC/Genlambda1520pt"), 1, part.pt(), centrality); else @@ -1152,6 +1155,7 @@ struct Lambda1520analysisinpp { } if (inVtx10 && isSel8) // vtx10, sel8 { + histos.fill(HIST("QA/Gen"), 3); if (part.pdgCode() > 0) histos.fill(HIST("Result/MC/Genlambda1520pt"), 2, part.pt(), centrality); else @@ -1159,6 +1163,7 @@ struct Lambda1520analysisinpp { } if (inVtx10 && isTriggerTVX) // vtx10, TriggerTVX { + histos.fill(HIST("QA/Gen"), 4); if (part.pdgCode() > 0) histos.fill(HIST("Result/MC/Genlambda1520pt"), 3, part.pt(), centrality); else @@ -1166,6 +1171,7 @@ struct Lambda1520analysisinpp { } if (inVtx10 && isTrueINELgt0) // vtx10, INEL>0 { + histos.fill(HIST("QA/Gen"), 5); if (part.pdgCode() > 0) histos.fill(HIST("Result/MC/Genlambda1520pt"), 4, part.pt(), centrality); else @@ -1173,6 +1179,7 @@ struct Lambda1520analysisinpp { } if (isInAfterAllCuts) // after all event selection { + histos.fill(HIST("QA/Gen"), 6); if (part.pdgCode() > 0) histos.fill(HIST("Result/MC/Genlambda1520pt"), 5, part.pt(), centrality); else @@ -1180,6 +1187,7 @@ struct Lambda1520analysisinpp { } if (isInAfterAllCuts && isTrueINELgt0) // after all event selection && INEL>0 { + histos.fill(HIST("QA/Gen"), 7); if (part.pdgCode() > 0) histos.fill(HIST("Result/MC/Genlambda1520pt"), 6, part.pt(), centrality); else