diff --git a/EventFiltering/PWGJE/jetFilter.cxx b/EventFiltering/PWGJE/jetFilter.cxx index cd7e219e371..acb47bdf1f1 100644 --- a/EventFiltering/PWGJE/jetFilter.cxx +++ b/EventFiltering/PWGJE/jetFilter.cxx @@ -43,23 +43,17 @@ using namespace o2::framework::expressions; static const std::vector highPtObjectsNames{"JetChHighPt"}; struct jetFilter { - enum { kJetChHighPt = 0, + enum { kJetChLowPt = 0, + kJetChHighPt = 1, kHighPtObjects }; - // event selection cuts - Configurable selectionJetChHighPt{ - "selectionJetChHighPt", 33., - "Minimum charged jet pT trigger threshold"}; // we want to keep all events - // having a charged jet with - // pT above this - Produces tags; Configurable cfgJetR{"cfgJetR", 0.6, "jet resolution parameter"}; // jet cone radius - Configurable cfgJetPtMin{ - "cfgJetPtMin", 0.1, - "minimum jet pT constituent cut"}; // minimum jet constituent pT + + Configurable jetPtLowThreshold{"jetPtLowThreshold", 30.0, "threshold for low pT jet trigger"}; + Configurable jetPtHighThreshold{"jetPtHighThreshold", 50.0, "threshold for high pT jet trigger"}; HistogramRegistry spectra{ "spectra", @@ -73,12 +67,20 @@ struct jetFilter { spectra.add("fCollZpos", "collision z position", HistType::kTH1F, {{200, -20., +20., "#it{z}_{vtx} position (cm)"}}); - spectra.add("ptphiJetChSelected", - "pT of selected high pT charged jets vs phi", HistType::kTH2F, + + spectra.add("ptphiJetChSelected_lowptjettrigger", "pT of selected low pT charged jet trigger vs phi", HistType::kTH2F, + {{150, 0., +150., "charged jet #it{p}_{T} (GeV/#it{c})"}, + {60, 0, TMath::TwoPi()}}); + + spectra.add("ptphiJetChSelected_highptjettrigger", "pT of selected high pT charged jet trigger vs phi", HistType::kTH2F, {{150, 0., +150., "charged jet #it{p}_{T} (GeV/#it{c})"}, {60, 0, TMath::TwoPi()}}); - spectra.add("ptetaJetChSelected", - "pT of selected high pT charged jets vs eta", HistType::kTH2F, + + spectra.add("ptetaJetChSelected_lowptjettrigger", "pT of selected low pT charged jet trigger vs eta", HistType::kTH2F, + {{150, 0., +150., "charged jet #it{p}_{T} (GeV/#it{c})"}, + {40, -1.0, 1.0}}); + + spectra.add("ptetaJetChSelected_highptjettrigger", "pT of selected high pT charged jet trigger vs eta", HistType::kTH2F, {{150, 0., +150., "charged jet #it{p}_{T} (GeV/#it{c})"}, {40, -1.0, 1.0}}); @@ -96,24 +98,25 @@ struct jetFilter { Filter jetRadiusSelection = o2::aod::jet::r == nround(cfgJetR.node() * 100.0f); using filteredJets = o2::soa::Filtered; - void - process(aod::JCollision const& collision, filteredJets const& jets) + void process(aod::JCollision const& collision, filteredJets const& jets) { // collision process loop bool keepEvent[kHighPtObjects]{false}; - // spectra.fill(HIST("fCollZpos"), collision.posZ()); - // Check whether there is a high pT charged jet - for (const auto& jet : jets) { - if (jet.pt() >= selectionJetChHighPt) { - spectra.fill(HIST("ptphiJetChSelected"), jet.pt(), - jet.phi()); // charged jet pT vs phi - spectra.fill(HIST("ptetaJetChSelected"), jet.pt(), - jet.eta()); // charged jet pT vs eta + for (const auto& jet : jets) { // jets are ordered by pT + + if (jet.pt() >= jetPtLowThreshold) { + spectra.fill(HIST("ptphiJetChSelected_lowptjettrigger"), jet.pt(), jet.phi()); // charged jet pT vs phi + spectra.fill(HIST("ptetaJetChSelected_lowptjettrigger"), jet.pt(), jet.eta()); // charged jet pT vs eta + keepEvent[kJetChLowPt] = true; + } + if (jet.pt() >= jetPtHighThreshold) { + spectra.fill(HIST("ptphiJetChSelected_highptjettrigger"), jet.pt(), jet.phi()); // charged jet pT vs phi + spectra.fill(HIST("ptetaJetChSelected_highptjettrigger"), jet.pt(), jet.eta()); // charged jet pT vs eta keepEvent[kJetChHighPt] = true; - break; } + break; // only looks at the highest pT jet in the event } for (int iDecision{0}; iDecision < kHighPtObjects; ++iDecision) { @@ -121,7 +124,7 @@ struct jetFilter { spectra.fill(HIST("fProcessedEvents"), iDecision); } } - tags(keepEvent[kJetChHighPt]); + tags(keepEvent[kJetChLowPt], keepEvent[kJetChHighPt]); } }; diff --git a/EventFiltering/filterTables.h b/EventFiltering/filterTables.h index 1df6da30140..3ab56d13088 100644 --- a/EventFiltering/filterTables.h +++ b/EventFiltering/filterTables.h @@ -60,6 +60,7 @@ DECLARE_SOA_COLUMN(PLL, hasPLL, bool); //! has p-L-L triplet DECLARE_SOA_COLUMN(LLL, hasLLL, bool); //! has L-L-L tripletD // jets +DECLARE_SOA_COLUMN(JetChLowPt, hasJetChLowPt, bool); //! low-pT charged jet DECLARE_SOA_COLUMN(JetChHighPt, hasJetChHighPt, bool); //! high-pT charged jet // full jets @@ -160,6 +161,7 @@ using CfFilter = CFFilters::iterator; // jets DECLARE_SOA_TABLE(JetFilters, "AOD", "JetFilters", //! + filtering::JetChLowPt, filtering::JetChHighPt); using JetFilter = JetFilters::iterator;