From af1e51c499d6ea8f7ff19a37eedb5a690c64259e Mon Sep 17 00:00:00 2001 From: sawan Date: Wed, 21 Oct 2026 11:13:56 +0530 Subject: [PATCH 1/2] modified event loss for INEL --- PWGLF/Tasks/Resonances/kstarqa.cxx | 47 +++++++++++++++++++++++++----- 1 file changed, 39 insertions(+), 8 deletions(-) diff --git a/PWGLF/Tasks/Resonances/kstarqa.cxx b/PWGLF/Tasks/Resonances/kstarqa.cxx index c2e2a545c8e..cf169d3cbae 100644 --- a/PWGLF/Tasks/Resonances/kstarqa.cxx +++ b/PWGLF/Tasks/Resonances/kstarqa.cxx @@ -232,11 +232,11 @@ struct Kstarqa { std::vector eveCutLabels = { "All Events", - "|Vz| < 10", - std::string("Sel8") + check(configGp.isSel8.value), std::string("kIsTriggerTVX") + check(configGp.isTriggerTVX.value), + "|Vz| < 10", std::string("NoTimeFrameBorder") + check(configGp.isNoTimeFrameBorder.value), std::string("NoITSROFrameBorder") + check(configGp.isNoITSROFrameBorder.value), + std::string("Sel8") + check(configGp.isSel8.value), std::string("rctChecker") + check(rctCut.requireRCTFlagChecker.value), std::string("IsINELgt0") + check(configGp.isINELgt0.value), std::string("IsGoodITSLayersAll") + check(configGp.isGoodITSLayersAll.value), @@ -355,6 +355,7 @@ struct Kstarqa { hInvMass.add("MCcorrections/MultiplicityRec", "Multiplicity in generated MC with at least 1 reconstruction", kTH1F, {multiplicityAxis}); hInvMass.add("MCcorrections/hImpactParameterGen", "Impact parameter in generated MC", kTH1F, {impactParAxis}); hInvMass.add("MCcorrections/MultiplicityGen", "Multiplicity in generated MC", kTH1F, {multiplicityAxis}); + hInvMass.add("MCcorrections/MultiplicityGenINEL", "Multiplicity in generated MC", kTH1F, {multiplicityAxis}); hInvMass.add("MCcorrections/hImpactParametervsMultiplicity", "Impact parameter vs multiplicity in reconstructed MC", kTH2F, {{impactParAxis}, {multiplicityAxis}}); } @@ -473,27 +474,27 @@ struct Kstarqa { if (fillHist) rEventSelection.fill(HIST("hEventCut"), 0); - if (std::abs(collision.posZ()) > configGp.cutzvertex) + if (configGp.isTriggerTVX && !collision.selection_bit(aod::evsel::kIsTriggerTVX)) return false; if (fillHist) rEventSelection.fill(HIST("hEventCut"), 1); - if (configGp.isSel8 && !collision.sel8()) + if (std::abs(collision.posZ()) > configGp.cutzvertex) return false; if (fillHist) rEventSelection.fill(HIST("hEventCut"), 2); - if (configGp.isTriggerTVX && !collision.selection_bit(aod::evsel::kIsTriggerTVX)) + if (configGp.isNoTimeFrameBorder && !collision.selection_bit(aod::evsel::kNoTimeFrameBorder)) return false; if (fillHist) rEventSelection.fill(HIST("hEventCut"), 3); - if (configGp.isNoTimeFrameBorder && !collision.selection_bit(aod::evsel::kNoTimeFrameBorder)) + if (configGp.isNoITSROFrameBorder && !collision.selection_bit(aod::evsel::kNoITSROFrameBorder)) return false; if (fillHist) rEventSelection.fill(HIST("hEventCut"), 4); - if (configGp.isNoITSROFrameBorder && !collision.selection_bit(aod::evsel::kNoITSROFrameBorder)) + if (configGp.isSel8 && !collision.sel8()) return false; if (fillHist) rEventSelection.fill(HIST("hEventCut"), 5); @@ -1911,6 +1912,7 @@ struct Kstarqa { { // auto multiplicityRec = -1; bool isSelectedEvent = false; + bool isSelectedEventINEL = false; auto multiplicity1 = -999.; if (configGp.isINELgt0Gen && !mcCollision.isInelGt0()) { @@ -1942,6 +1944,35 @@ struct Kstarqa { isSelectedEvent = true; } + for (const auto& RecCollision : recCollisions) { + if (!RecCollision.has_mcCollision()) + continue; + + if (configGp.isTriggerTVX && !RecCollision.selection_bit(aod::evsel::kIsTriggerTVX)) + return false; + + if (std::abs(RecCollision.posZ()) > configGp.cutzvertex) + return false; + + if (configGp.isNoTimeFrameBorder && !RecCollision.selection_bit(aod::evsel::kNoTimeFrameBorder)) + return false; + + if (configGp.isNoITSROFrameBorder && !RecCollision.selection_bit(aod::evsel::kNoITSROFrameBorder)) + return false; + + if (configGp.isSel8 && !RecCollision.sel8()) + return false; + + if (rctCut.requireRCTFlagChecker && !rctChecker(RecCollision)) + return false; + + isSelectedEventINEL = true; + } + + if (isSelectedEventINEL) { + hInvMass.fill(HIST("MCcorrections/MultiplicityGenINEL"), multiplicityGen); + } + bool isINELgt0true = false; if (pwglf::isINELgtNmc(mcParticles, 0, pdgDB)) { @@ -3200,7 +3231,7 @@ struct Kstarqa { genDaughter1 = ROOT::Math::PxPyPzMVector(mctrack1.px(), mctrack1.py(), mctrack1.pz(), massKa); genDaughter2 = ROOT::Math::PxPyPzMVector(mctrack2.px(), mctrack2.py(), mctrack2.pz(), massKa); - mother = daughter1 + daughter2; // Phi meson + mother = daughter1 + daughter2; // Phi meson genMother = genDaughter1 + genDaughter2; // Phi meson recMass = mother.M(); From da2f9fc1d2c2acf5062d5f8ab217f4a8c82f6567 Mon Sep 17 00:00:00 2001 From: sawan Date: Wed, 21 Oct 2026 11:26:00 +0530 Subject: [PATCH 2/2] fixed errors --- PWGLF/Tasks/Resonances/kstarqa.cxx | 22 +++++++++++----------- 1 file changed, 11 insertions(+), 11 deletions(-) diff --git a/PWGLF/Tasks/Resonances/kstarqa.cxx b/PWGLF/Tasks/Resonances/kstarqa.cxx index cf169d3cbae..21e7db35a50 100644 --- a/PWGLF/Tasks/Resonances/kstarqa.cxx +++ b/PWGLF/Tasks/Resonances/kstarqa.cxx @@ -355,7 +355,7 @@ struct Kstarqa { hInvMass.add("MCcorrections/MultiplicityRec", "Multiplicity in generated MC with at least 1 reconstruction", kTH1F, {multiplicityAxis}); hInvMass.add("MCcorrections/hImpactParameterGen", "Impact parameter in generated MC", kTH1F, {impactParAxis}); hInvMass.add("MCcorrections/MultiplicityGen", "Multiplicity in generated MC", kTH1F, {multiplicityAxis}); - hInvMass.add("MCcorrections/MultiplicityGenINEL", "Multiplicity in generated MC", kTH1F, {multiplicityAxis}); + hInvMass.add("MCcorrections/MultiplicityRecINEL", "Multiplicity in generated MC", kTH1F, {multiplicityAxis}); hInvMass.add("MCcorrections/hImpactParametervsMultiplicity", "Impact parameter vs multiplicity in reconstructed MC", kTH2F, {{impactParAxis}, {multiplicityAxis}}); } @@ -1949,30 +1949,26 @@ struct Kstarqa { continue; if (configGp.isTriggerTVX && !RecCollision.selection_bit(aod::evsel::kIsTriggerTVX)) - return false; + continue; if (std::abs(RecCollision.posZ()) > configGp.cutzvertex) - return false; + continue; if (configGp.isNoTimeFrameBorder && !RecCollision.selection_bit(aod::evsel::kNoTimeFrameBorder)) - return false; + continue; if (configGp.isNoITSROFrameBorder && !RecCollision.selection_bit(aod::evsel::kNoITSROFrameBorder)) - return false; + continue; if (configGp.isSel8 && !RecCollision.sel8()) - return false; + continue; if (rctCut.requireRCTFlagChecker && !rctChecker(RecCollision)) - return false; + continue; isSelectedEventINEL = true; } - if (isSelectedEventINEL) { - hInvMass.fill(HIST("MCcorrections/MultiplicityGenINEL"), multiplicityGen); - } - bool isINELgt0true = false; if (pwglf::isINELgtNmc(mcParticles, 0, pdgDB)) { @@ -1988,6 +1984,10 @@ struct Kstarqa { hInvMass.fill(HIST("MCcorrections/hImpactParameterGen"), impactPar); hInvMass.fill(HIST("MCcorrections/MultiplicityGen"), multiplicityGen); + if (isSelectedEventINEL) { + hInvMass.fill(HIST("MCcorrections/MultiplicityRecINEL"), multiplicityGen); + } + // Event loss if (isSelectedEvent) { hInvMass.fill(HIST("MCcorrections/hImpactParameterRec"), impactPar);