From 7c180d63fbb2fad3e766938fb213762f25f57b18 Mon Sep 17 00:00:00 2001 From: sawan <124118453+sawankumawat@users.noreply.github.com> Date: Sat, 18 Jan 2025 01:12:42 +0530 Subject: [PATCH] [PWGLF] Included all resonances for generated and rec MC (#9358) --- PWGLF/Tasks/Resonances/higherMassResonances.cxx | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/PWGLF/Tasks/Resonances/higherMassResonances.cxx b/PWGLF/Tasks/Resonances/higherMassResonances.cxx index f6de828ef46..ac1d54bf404 100644 --- a/PWGLF/Tasks/Resonances/higherMassResonances.cxx +++ b/PWGLF/Tasks/Resonances/higherMassResonances.cxx @@ -119,6 +119,8 @@ struct HigherMassResonances { // Configurable for MC Configurable allGenCollisions{"allGenCollisions", true, "To fill all generated collisions for the signal loss calculations"}; Configurable cTVXEvsel{"cTVXEvsel", true, "Triggger selection"}; + Configurable selectMCparticles{"selectMCparticles", 1, "0: f0(1710), 1: f2(1525), 2: a2(1320), 3: f0(1370), 4: f0(1500)"}; + std::vector pdgCodes = {10331, 335, 115, 10221, 9030221}; // output THnSparses Configurable activateTHnSparseCosThStarHelicity{"activateTHnSparseCosThStarHelicity", false, "Activate the THnSparse with cosThStar w.r.t. helicity axis"}; @@ -1014,7 +1016,7 @@ struct HigherMassResonances { // std::cout << "PDG code mother " << mcParticle.pdgCode() << std::endl; // } // counter++; - if (std::abs(mcParticle.pdgCode()) != 10331) // f2(1525), f0(1710) + if (std::abs(mcParticle.pdgCode()) != pdgCodes[selectMCparticles]) // f2(1525), f0(1710) { continue; } @@ -1181,7 +1183,7 @@ struct HigherMassResonances { // if (counter2 < 1e4) // std::cout << "Mother2 pdg code: " << motpdgs2 << " p_{T} " << mothertrack2.pt() << "Global index " << mothertrack1.globalIndex() << " event " << eventCounter << std::endl; - if (mothertrack1.pdgCode() != 10331) { + if (mothertrack1.pdgCode() != pdgCodes[selectMCparticles]) { continue; } hMChists.fill(HIST("events_checkrec"), 15.5);