From 29319b990d0313e72657a1fbac8aa3ef62a485dd Mon Sep 17 00:00:00 2001 From: Sawan Sawan Date: Wed, 31 Dec 2025 11:28:10 +0530 Subject: [PATCH 1/6] added sel8 configurable, to remove it during MC checks --- PWGLF/Tasks/Resonances/kstarqa.cxx | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/PWGLF/Tasks/Resonances/kstarqa.cxx b/PWGLF/Tasks/Resonances/kstarqa.cxx index 7d530010757..18fb2df2477 100644 --- a/PWGLF/Tasks/Resonances/kstarqa.cxx +++ b/PWGLF/Tasks/Resonances/kstarqa.cxx @@ -78,6 +78,7 @@ struct Kstarqa { struct : ConfigurableGroup { // Configurables for event selections Configurable isINELgt0{"isINELgt0", true, "INEL>0 selection"}; + Configurable isSel8{"isSel8", false, "Event selection sel8"}; Configurable isTriggerTVX{"isTriggerTVX", false, "TriggerTVX"}; Configurable isGoodZvtxFT0vsPV{"isGoodZvtxFT0vsPV", false, "IsGoodZvtxFT0vsPV"}; Configurable isApplyOccCut{"isApplyOccCut", true, "Apply occupancy cut"}; @@ -446,7 +447,7 @@ struct Kstarqa { if (fillHist) rEventSelection.fill(HIST("hEventCut"), 1); - if (!collision.sel8()) + if (selectionConfig.isSel8 && !collision.sel8()) return false; if (fillHist) rEventSelection.fill(HIST("hEventCut"), 2); From a300a14922c5e5b40ec2eabe0fc82cf5fa62d453 Mon Sep 17 00:00:00 2001 From: Sawan Sawan Date: Wed, 31 Dec 2025 12:43:29 +0530 Subject: [PATCH 2/6] different event/signal loss method --- PWGLF/Tasks/Resonances/kstarqa.cxx | 19 ++++++++++++++++++- 1 file changed, 18 insertions(+), 1 deletion(-) diff --git a/PWGLF/Tasks/Resonances/kstarqa.cxx b/PWGLF/Tasks/Resonances/kstarqa.cxx index 18fb2df2477..0849722b6f5 100644 --- a/PWGLF/Tasks/Resonances/kstarqa.cxx +++ b/PWGLF/Tasks/Resonances/kstarqa.cxx @@ -378,8 +378,13 @@ struct Kstarqa { hGenEvents->GetYaxis()->SetBinLabel(2, "Generated events with Mc collision V_{z} cut"); hGenEvents->GetYaxis()->SetBinLabel(3, "Generated events with Mc INEL>0"); hGenEvents->GetYaxis()->SetBinLabel(4, "Generated events with at least one reconstructed event"); + hInvMass.add("CorrFactors/h2dGenKstar", "Centrality vs p_{T}", kTH2D, {{101, 0.0f, 101.0f}, ptAxis}); hInvMass.add("CorrFactors/h3dGenKstarVsMultMCVsMultiplicity", "MC centrality vs centrality vs p_{T}", kTH3D, {axisNch, {101, 0.0f, 101.0f}, ptAxis}); + hInvMass.add("CorrFactors/hSignalLossDenominator", "Kstar generated before event selection", kTH2F, {{ptAxis}, {multiplicityAxis}}); + hInvMass.add("CorrFactors/hSignalLossNumerator", "Kstar generated after event selection", kTH2F, {{ptAxis}, {multiplicityAxis}}); + hInvMass.add("CorrFactors/MultiplicityRec", "Multiplicity in generated MC with at least 1 reconstruction", kTH1F, {multiplicityAxis}); + hInvMass.add("CorrFactors/MultiplicityGen", "Multiplicity in generated MC", kTH1F, {multiplicityAxis}); } rEventSelection.add("tracksCheckData", "No. of events in the data", kTH1I, {{10, 0, 10}}); @@ -1673,12 +1678,13 @@ struct Kstarqa { auto multiplicityRec = -1; bool isSelectedEvent = false; auto multiplicity1 = -999.; + for (const auto& RecCollision : recCollisions) { if (!RecCollision.has_mcCollision()) continue; if (!selectionEvent(RecCollision, false)) // don't fill event cut histogram continue; - // multiplicity1 = RecCollision.centFT0M(); + const auto& mcCollisionRec = RecCollision.mcCollision_as(); multiplicityRec = mcCollisionRec.centFT0M(); @@ -1752,6 +1758,8 @@ struct Kstarqa { hInvMass.fill(HIST("CorrFactors/hGenEvents"), mcCollision.multMCNParticlesEta08(), 2.5); float multiplicity = 100.5f; + bool isSelectedEvent = false; + for (auto const& collision : collisions) { if (!selectionEvent(collision, false)) // don't fill event cut histogram continue; @@ -1767,10 +1775,15 @@ struct Kstarqa { } else { multiplicity = collision.centFT0M(); // default } + isSelectedEvent = true; } hInvMass.fill(HIST("CorrFactors/hMultiplicityVsMultMC"), multiplicity, mcCollision.multMCNParticlesEta08()); hInvMass.fill(HIST("CorrFactors/hNrecInGen"), collisions.size()); + hInvMass.fill(HIST("CorrFactors/MultiplicityGen"), multiplicity); + if (isSelectedEvent) { + hInvMass.fill(HIST("CorrFactors/MultiplicityRec"), multiplicity); + } for (const auto& mcParticle : mcParticles) { @@ -1816,6 +1829,10 @@ struct Kstarqa { hInvMass.fill(HIST("CorrFactors/h2dGenKstar"), multiplicity, mother.Pt()); hInvMass.fill(HIST("CorrFactors/h3dGenKstarVsMultMCVsMultiplicity"), mcCollision.multMCNParticlesEta08(), multiplicity, mother.Pt()); + hInvMass.fill(HIST("CorrFactors/hSignalLossDenominator"), mother.pt(), multiplicity); + if (isSelectedEvent) { + hInvMass.fill(HIST("CorrFactors/hSignalLossNumerator"), mother.pt(), multiplicity); + } } } } From 3e660b72a82a477055c49151a5257a8443d84e53 Mon Sep 17 00:00:00 2001 From: Sawan Sawan Date: Wed, 31 Dec 2025 14:42:24 +0530 Subject: [PATCH 3/6] added another method for event/signal loss corrections --- PWGLF/Tasks/Resonances/kstarqa.cxx | 41 +++++++++++++++++++++--------- 1 file changed, 29 insertions(+), 12 deletions(-) diff --git a/PWGLF/Tasks/Resonances/kstarqa.cxx b/PWGLF/Tasks/Resonances/kstarqa.cxx index 0849722b6f5..590c246c6bb 100644 --- a/PWGLF/Tasks/Resonances/kstarqa.cxx +++ b/PWGLF/Tasks/Resonances/kstarqa.cxx @@ -347,6 +347,10 @@ struct Kstarqa { hInvMass.add("hAllKstarGenCollisisons1Rec", "All generated Kstar in events with at least one rec event in rapidity in 0.5", kTH2F, {{multiplicityAxis}, {ptAxis}}); hInvMass.add("hAllRecCollisions", "All reconstructed events", kTH1F, {multiplicityAxis}); hInvMass.add("hAllRecCollisionsCalib", "All reconstructed events", kTH1F, {multiplicityAxis}); + hInvMass.add("sigEvLossFromGenRec/MultiplicityGen", "Multiplicity in generated MC", kTH1F, {multiplicityAxis}); + hInvMass.add("sigEvLossFromGenRec/MultiplicityRec", "Multiplicity in generated MC with at least 1 reconstruction", kTH1F, {multiplicityAxis}); + hInvMass.add("sigEvLossFromGenRec/hSignalLossDenominator", "Kstar generated before event selection", kTH2F, {{ptAxis}, {multiplicityAxis}}); + hInvMass.add("sigEvLossFromGenRec/hSignalLossNumerator", "Kstar generated after event selection", kTH2F, {{ptAxis}, {multiplicityAxis}}); if (doprocessEvtLossSigLossMC || doprocessEvtLossSigLossMCPhi) { hInvMass.add("MCcorrections/hSignalLossDenominator", "Kstar generated before event selection", kTH2F, {{ptAxis}, {multiplicityAxis}}); @@ -1611,15 +1615,23 @@ struct Kstarqa { } } - const auto evtReconstructedAndSelected = std::find(selectedEvents.begin(), selectedEvents.end(), mcCollision.globalIndex()) != selectedEvents.end(); - hInvMass.fill(HIST("hAllGenCollisions"), multiplicity); - if (!evtReconstructedAndSelected) { // Check that the event is reconstructed and that the reconstructed events pass the selection + if (std::abs(mcCollision.posZ()) < selectionConfig.cutzvertex) { return; } + + const auto evtReconstructedAndSelected = std::find(selectedEvents.begin(), selectedEvents.end(), mcCollision.globalIndex()) != selectedEvents.end(); + hInvMass.fill(HIST("hAllGenCollisions"), multiplicity); + // if (!evtReconstructedAndSelected) { // Check that the event is reconstructed and that the reconstructed events pass the selection + // return; + // } double genMultiplicity = mcCollision.centFT0M(); - hInvMass.fill(HIST("h1GenMult2"), genMultiplicity); - hInvMass.fill(HIST("hAllGenCollisions1Rec"), multiplicity); - rEventSelection.fill(HIST("eventsCheckGen"), 3.5); + hInvMass.fill(HIST("sigEvLossFromGenRec/MultiplicityGen"), genMultiplicity); + if (evtReconstructedAndSelected) { + hInvMass.fill(HIST("sigEvLossFromGenRec/MultiplicityRec"), genMultiplicity); + hInvMass.fill(HIST("h1GenMult2"), genMultiplicity); + hInvMass.fill(HIST("hAllGenCollisions1Rec"), multiplicity); + rEventSelection.fill(HIST("eventsCheckGen"), 3.5); + } for (const auto& mcParticle : mcParticles) { @@ -1637,7 +1649,8 @@ struct Kstarqa { if (std::abs(mcParticle.pdgCode()) != o2::constants::physics::kK0Star892) { continue; } - hInvMass.fill(HIST("hAllKstarGenCollisisons1Rec"), multiplicity, mcParticle.pt()); + if (evtReconstructedAndSelected) + hInvMass.fill(HIST("hAllKstarGenCollisisons1Rec"), multiplicity, mcParticle.pt()); auto kDaughters = mcParticle.daughters_as(); if (kDaughters.size() != selectionConfig.noOfDaughters) { @@ -1662,11 +1675,15 @@ struct Kstarqa { } if (passkaon && passpion) { mother = daughter1 + daughter2; // Kstar meson - hInvMass.fill(HIST("hk892GenpT"), mcParticle.pt(), multiplicity); - hInvMass.fill(HIST("hk892GenpT2"), mother.Pt(), multiplicity); - hInvMass.fill(HIST("hk892GenpTCalib1"), mcParticle.pt(), genMultiplicity); - hInvMass.fill(HIST("hk892GenpTCalib2"), mother.Pt(), genMultiplicity); - hInvMass.fill(HIST("h1genmass"), mother.M()); + hInvMass.fill(HIST("sigEvLossFromGenRec/hSignalLossDenominator"), mother.Pt(), genMultiplicity); + if (evtReconstructedAndSelected) { + hInvMass.fill(HIST("hk892GenpT"), mcParticle.pt(), multiplicity); + hInvMass.fill(HIST("hk892GenpT2"), mother.Pt(), multiplicity); + hInvMass.fill(HIST("hk892GenpTCalib1"), mcParticle.pt(), genMultiplicity); + hInvMass.fill(HIST("hk892GenpTCalib2"), mother.Pt(), genMultiplicity); + hInvMass.fill(HIST("h1genmass"), mother.M()); + hInvMass.fill(HIST("sigEvLossFromGenRec/hSignalLossNumerator"), mother.Pt(), genMultiplicity); + } } } } From f5981ec4eb191c2d0cb93cc55686c78576b7695e Mon Sep 17 00:00:00 2001 From: Sawan Sawan Date: Thu, 1 Jan 2026 23:03:39 +0530 Subject: [PATCH 4/6] minor fix to code --- PWGLF/Tasks/Resonances/kstarqa.cxx | 23 ++++++++++++++--------- 1 file changed, 14 insertions(+), 9 deletions(-) diff --git a/PWGLF/Tasks/Resonances/kstarqa.cxx b/PWGLF/Tasks/Resonances/kstarqa.cxx index 590c246c6bb..57a33bb072f 100644 --- a/PWGLF/Tasks/Resonances/kstarqa.cxx +++ b/PWGLF/Tasks/Resonances/kstarqa.cxx @@ -385,10 +385,10 @@ struct Kstarqa { hInvMass.add("CorrFactors/h2dGenKstar", "Centrality vs p_{T}", kTH2D, {{101, 0.0f, 101.0f}, ptAxis}); hInvMass.add("CorrFactors/h3dGenKstarVsMultMCVsMultiplicity", "MC centrality vs centrality vs p_{T}", kTH3D, {axisNch, {101, 0.0f, 101.0f}, ptAxis}); - hInvMass.add("CorrFactors/hSignalLossDenominator", "Kstar generated before event selection", kTH2F, {{ptAxis}, {multiplicityAxis}}); - hInvMass.add("CorrFactors/hSignalLossNumerator", "Kstar generated after event selection", kTH2F, {{ptAxis}, {multiplicityAxis}}); - hInvMass.add("CorrFactors/MultiplicityRec", "Multiplicity in generated MC with at least 1 reconstruction", kTH1F, {multiplicityAxis}); - hInvMass.add("CorrFactors/MultiplicityGen", "Multiplicity in generated MC", kTH1F, {multiplicityAxis}); + hInvMass.add("CorrFactors/hSignalLossDenominator3D", "Kstar generated before event selection", kTH3F, {{ptAxis}, {multiplicityAxis}, axisNch}); + hInvMass.add("CorrFactors/hSignalLossNumerator3D", "Kstar generated after event selection", kTH3F, {{ptAxis}, {multiplicityAxis}, axisNch}); + hInvMass.add("CorrFactors/MultiplicityRec2D", "Multiplicity in generated MC with at least 1 reconstruction", kTH2F, {{multiplicityAxis}, axisNch}); + hInvMass.add("CorrFactors/MultiplicityGen2D", "Multiplicity in generated MC", kTH2F, {{multiplicityAxis}, axisNch}); } rEventSelection.add("tracksCheckData", "No. of events in the data", kTH1I, {{10, 0, 10}}); @@ -1615,7 +1615,7 @@ struct Kstarqa { } } - if (std::abs(mcCollision.posZ()) < selectionConfig.cutzvertex) { + if (std::abs(mcCollision.posZ()) > selectionConfig.cutzvertex) { return; } @@ -1778,6 +1778,8 @@ struct Kstarqa { bool isSelectedEvent = false; for (auto const& collision : collisions) { + if (!collision.has_mcCollision()) + continue; if (!selectionEvent(collision, false)) // don't fill event cut histogram continue; @@ -1795,11 +1797,14 @@ struct Kstarqa { isSelectedEvent = true; } + // auto multiplicityGen = -1; + // multiplicityGen = mcCollision.centFT0M(); + hInvMass.fill(HIST("CorrFactors/hMultiplicityVsMultMC"), multiplicity, mcCollision.multMCNParticlesEta08()); hInvMass.fill(HIST("CorrFactors/hNrecInGen"), collisions.size()); - hInvMass.fill(HIST("CorrFactors/MultiplicityGen"), multiplicity); + hInvMass.fill(HIST("CorrFactors/MultiplicityGen2D"), multiplicity, mcCollision.multMCNParticlesEta08()); if (isSelectedEvent) { - hInvMass.fill(HIST("CorrFactors/MultiplicityRec"), multiplicity); + hInvMass.fill(HIST("CorrFactors/MultiplicityRec2D"), multiplicity, mcCollision.multMCNParticlesEta08()); } for (const auto& mcParticle : mcParticles) { @@ -1846,9 +1851,9 @@ struct Kstarqa { hInvMass.fill(HIST("CorrFactors/h2dGenKstar"), multiplicity, mother.Pt()); hInvMass.fill(HIST("CorrFactors/h3dGenKstarVsMultMCVsMultiplicity"), mcCollision.multMCNParticlesEta08(), multiplicity, mother.Pt()); - hInvMass.fill(HIST("CorrFactors/hSignalLossDenominator"), mother.pt(), multiplicity); + hInvMass.fill(HIST("CorrFactors/hSignalLossDenominator3D"), mother.pt(), multiplicity, mcCollision.multMCNParticlesEta08()); if (isSelectedEvent) { - hInvMass.fill(HIST("CorrFactors/hSignalLossNumerator"), mother.pt(), multiplicity); + hInvMass.fill(HIST("CorrFactors/hSignalLossNumerator3D"), mother.pt(), multiplicity, mcCollision.multMCNParticlesEta08()); } } } From face16c3a6803c7a6b9f9b8aec95abd76c6a1194 Mon Sep 17 00:00:00 2001 From: Sawan Sawan Date: Fri, 2 Jan 2026 19:26:41 +0530 Subject: [PATCH 5/6] corrected event loss histograms --- PWGLF/Tasks/Resonances/kstarqa.cxx | 40 ++++++++++++++---------------- 1 file changed, 19 insertions(+), 21 deletions(-) diff --git a/PWGLF/Tasks/Resonances/kstarqa.cxx b/PWGLF/Tasks/Resonances/kstarqa.cxx index 57a33bb072f..a0887e9a1a7 100644 --- a/PWGLF/Tasks/Resonances/kstarqa.cxx +++ b/PWGLF/Tasks/Resonances/kstarqa.cxx @@ -357,7 +357,6 @@ struct Kstarqa { hInvMass.add("MCcorrections/hSignalLossNumerator", "Kstar generated after event selection", kTH2F, {{ptAxis}, {multiplicityAxis}}); hInvMass.add("MCcorrections/hImpactParameterRec", "Impact parameter in reconstructed MC", kTH1F, {impactParAxis}); hInvMass.add("MCcorrections/MultiplicityRec", "Multiplicity in generated MC with at least 1 reconstruction", kTH1F, {multiplicityAxis}); - hInvMass.add("MCcorrections/MultiplicityRec2", "Multiplicity in reconstructed MC", 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/hImpactParametervsMultiplicity", "Impact parameter vs multiplicity in reconstructed MC", kTH2F, {{impactParAxis}, {multiplicityAxis}}); @@ -385,10 +384,10 @@ struct Kstarqa { hInvMass.add("CorrFactors/h2dGenKstar", "Centrality vs p_{T}", kTH2D, {{101, 0.0f, 101.0f}, ptAxis}); hInvMass.add("CorrFactors/h3dGenKstarVsMultMCVsMultiplicity", "MC centrality vs centrality vs p_{T}", kTH3D, {axisNch, {101, 0.0f, 101.0f}, ptAxis}); - hInvMass.add("CorrFactors/hSignalLossDenominator3D", "Kstar generated before event selection", kTH3F, {{ptAxis}, {multiplicityAxis}, axisNch}); - hInvMass.add("CorrFactors/hSignalLossNumerator3D", "Kstar generated after event selection", kTH3F, {{ptAxis}, {multiplicityAxis}, axisNch}); - hInvMass.add("CorrFactors/MultiplicityRec2D", "Multiplicity in generated MC with at least 1 reconstruction", kTH2F, {{multiplicityAxis}, axisNch}); - hInvMass.add("CorrFactors/MultiplicityGen2D", "Multiplicity in generated MC", kTH2F, {{multiplicityAxis}, axisNch}); + hInvMass.add("CorrFactors/hSignalLossDenominator", "Kstar generated before event selection", kTH2F, {{ptAxis}, {axisNch}}); + hInvMass.add("CorrFactors/hSignalLossNumerator", "Kstar generated after event selection", kTH2F, {{ptAxis}, {axisNch}}); + hInvMass.add("CorrFactors/MultiplicityRec", "Multiplicity in generated MC with at least 1 reconstruction", kTH1F, {axisNch}); + hInvMass.add("CorrFactors/MultiplicityGen", "Multiplicity in generated MC", kTH1F, {axisNch}); } rEventSelection.add("tracksCheckData", "No. of events in the data", kTH1I, {{10, 0, 10}}); @@ -1743,7 +1742,6 @@ struct Kstarqa { if (isSelectedEvent) { hInvMass.fill(HIST("MCcorrections/hImpactParameterRec"), impactPar); hInvMass.fill(HIST("MCcorrections/MultiplicityRec"), multiplicityGen); - hInvMass.fill(HIST("MCcorrections/MultiplicityRec2"), multiplicityRec); hInvMass.fill(HIST("MCcorrections/hImpactParametervsMultiplicity"), impactPar, multiplicity1); } @@ -1762,19 +1760,21 @@ struct Kstarqa { void processEvSigLossFactors(McCollisionMults::iterator const& mcCollision, soa::SmallGroups const& collisions, LabeledTracks const&, aod::McParticles const& mcParticles) { - hInvMass.fill(HIST("CorrFactors/hGenEvents"), mcCollision.multMCNParticlesEta08(), 0.5); + auto multiplicityNch = -1; + multiplicityNch = mcCollision.multMCNParticlesEta05(); + hInvMass.fill(HIST("CorrFactors/hGenEvents"), multiplicityNch, 0.5); - if (std::abs(mcCollision.posZ()) > selectionConfig.cutzvertex) + if (selectionConfig.checkVzEvSigLoss && std::abs(mcCollision.posZ()) > selectionConfig.cutzvertex) return; - hInvMass.fill(HIST("CorrFactors/hGenEvents"), mcCollision.multMCNParticlesEta08(), 1.5); + hInvMass.fill(HIST("CorrFactors/hGenEvents"), multiplicityNch, 1.5); if (selectionConfig.isINELgt0 && !mcCollision.isInelGt0()) { return; } - hInvMass.fill(HIST("CorrFactors/hGenEvents"), mcCollision.multMCNParticlesEta08(), 2.5); + hInvMass.fill(HIST("CorrFactors/hGenEvents"), multiplicityNch, 2.5); - float multiplicity = 100.5f; + float multiplicity = -1.0; bool isSelectedEvent = false; for (auto const& collision : collisions) { @@ -1800,11 +1800,11 @@ struct Kstarqa { // auto multiplicityGen = -1; // multiplicityGen = mcCollision.centFT0M(); - hInvMass.fill(HIST("CorrFactors/hMultiplicityVsMultMC"), multiplicity, mcCollision.multMCNParticlesEta08()); + hInvMass.fill(HIST("CorrFactors/hMultiplicityVsMultMC"), multiplicity, multiplicityNch); hInvMass.fill(HIST("CorrFactors/hNrecInGen"), collisions.size()); - hInvMass.fill(HIST("CorrFactors/MultiplicityGen2D"), multiplicity, mcCollision.multMCNParticlesEta08()); + hInvMass.fill(HIST("CorrFactors/MultiplicityGen"), multiplicityNch); if (isSelectedEvent) { - hInvMass.fill(HIST("CorrFactors/MultiplicityRec2D"), multiplicity, mcCollision.multMCNParticlesEta08()); + hInvMass.fill(HIST("CorrFactors/MultiplicityRec"), multiplicityNch); } for (const auto& mcParticle : mcParticles) { @@ -1850,10 +1850,10 @@ struct Kstarqa { mother = daughter1 + daughter2; // Kstar meson hInvMass.fill(HIST("CorrFactors/h2dGenKstar"), multiplicity, mother.Pt()); - hInvMass.fill(HIST("CorrFactors/h3dGenKstarVsMultMCVsMultiplicity"), mcCollision.multMCNParticlesEta08(), multiplicity, mother.Pt()); - hInvMass.fill(HIST("CorrFactors/hSignalLossDenominator3D"), mother.pt(), multiplicity, mcCollision.multMCNParticlesEta08()); + hInvMass.fill(HIST("CorrFactors/h3dGenKstarVsMultMCVsMultiplicity"), multiplicityNch, multiplicity, mother.Pt()); + hInvMass.fill(HIST("CorrFactors/hSignalLossDenominator"), mother.pt(), multiplicityNch); if (isSelectedEvent) { - hInvMass.fill(HIST("CorrFactors/hSignalLossNumerator3D"), mother.pt(), multiplicity, mcCollision.multMCNParticlesEta08()); + hInvMass.fill(HIST("CorrFactors/hSignalLossNumerator"), mother.pt(), multiplicityNch); } } } @@ -1862,7 +1862,7 @@ struct Kstarqa { if (collisions.size() == 0) return; - hInvMass.fill(HIST("CorrFactors/hGenEvents"), mcCollision.multMCNParticlesEta08(), 3.5); + hInvMass.fill(HIST("CorrFactors/hGenEvents"), multiplicityNch, 3.5); } PROCESS_SWITCH(Kstarqa, processEvSigLossFactors, "Process Event and Signal loss", false); @@ -2171,7 +2171,7 @@ struct Kstarqa { hInvMass.fill(HIST("hAllRecCollisions"), multiplicity); hInvMass.fill(HIST("hAllRecCollisionsCalib"), multiplicityRec); - if (!selectionEvent(collision, false)) { // don't fill event cut histogram + if (!selectionEvent(collision, true)) { // don't fill event cut histogram return; } @@ -3079,7 +3079,6 @@ struct Kstarqa { } else { multiplicity1 = RecCollision.centFT0M(); // default } - isSelectedEvent = true; } @@ -3087,7 +3086,6 @@ struct Kstarqa { if (isSelectedEvent) { hInvMass.fill(HIST("MCcorrections/hImpactParameterRec"), impactPar); hInvMass.fill(HIST("MCcorrections/MultiplicityRec"), multiplicityGen); - hInvMass.fill(HIST("MCcorrections/MultiplicityRec2"), multiplicityRec); hInvMass.fill(HIST("MCcorrections/hImpactParametervsMultiplicity"), impactPar, multiplicity1); } From 0b6982e33eaa7226b3762f77f080c8643c59e32c Mon Sep 17 00:00:00 2001 From: Sawan Sawan Date: Fri, 2 Jan 2026 21:34:36 +0530 Subject: [PATCH 6/6] removed unused variable --- PWGLF/Tasks/Resonances/kstarqa.cxx | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/PWGLF/Tasks/Resonances/kstarqa.cxx b/PWGLF/Tasks/Resonances/kstarqa.cxx index a0887e9a1a7..ffe82bef9b4 100644 --- a/PWGLF/Tasks/Resonances/kstarqa.cxx +++ b/PWGLF/Tasks/Resonances/kstarqa.cxx @@ -1691,7 +1691,7 @@ struct Kstarqa { void processEvtLossSigLossMC(EventMCGenerated::iterator const& mcCollision, aod::McParticles const& mcParticles, const soa::SmallGroups& recCollisions) // void processEvtLossSigLossMC(aod::McCollisions::iterator const& mcCollision, aod::McParticles const& mcParticles, const soa::SmallGroups& recCollisions) { - auto multiplicityRec = -1; + // auto multiplicityRec = -1; bool isSelectedEvent = false; auto multiplicity1 = -999.; @@ -1701,8 +1701,8 @@ struct Kstarqa { if (!selectionEvent(RecCollision, false)) // don't fill event cut histogram continue; - const auto& mcCollisionRec = RecCollision.mcCollision_as(); - multiplicityRec = mcCollisionRec.centFT0M(); + // const auto& mcCollisionRec = RecCollision.mcCollision_as(); + // multiplicityRec = mcCollisionRec.centFT0M(); if (cSelectMultEstimator == kFT0M) { multiplicity1 = RecCollision.centFT0M(); @@ -3051,7 +3051,7 @@ struct Kstarqa { } auto impactPar = mcCollision.impactParameter(); - auto multiplicityRec = -1; + // auto multiplicityRec = -1; auto multiplicityGen = -1; multiplicityGen = mcCollision.centFT0M(); hInvMass.fill(HIST("MCcorrections/hImpactParameterGen"), impactPar); @@ -3065,8 +3065,8 @@ struct Kstarqa { if (!selectionEvent(RecCollision, false)) // don't fill event cut histogram continue; // multiplicity1 = RecCollision.centFT0M(); - const auto& mcCollisionRec = RecCollision.mcCollision_as(); - multiplicityRec = mcCollisionRec.centFT0M(); + // const auto& mcCollisionRec = RecCollision.mcCollision_as(); + // multiplicityRec = mcCollisionRec.centFT0M(); if (cSelectMultEstimator == kFT0M) { multiplicity1 = RecCollision.centFT0M();