Low Pt Track Selection (#5594)

* Low Pt Track Selection

* Please consider the following formatting changes

---------

Co-authored-by: ALICE Action Bot <alibuild@cern.ch>

Author: fkrizek

EventFiltering/PWGJE/jetFilter.cxx

@@ -44,15 +44,17 @@ using namespace o2::framework::expressions;
 struct jetFilter {
   enum { kJetChLowPt = 0,
          kJetChHighPt = 1,
-         kTrackPt = 2,
-         kTriggerObjects = 3
+         kTrackLowPt = 2,
+         kTrackHighPt = 3,
+         kTriggerObjects = 4
   };
 
   enum { kBinAllEvents = 0,
          kBinJetChLowPt = 1,
          kBinJetChHighPt = 2,
-         kBinTrackPt = 3,
-         kBins = 4
+         kBinTrackLowPt = 3,
+         kBinTrackHighPt = 4,
+         kBins = 5
   };
 
   Produces<aod::JetFilters> tags;
@@ -62,7 +64,8 @@ struct jetFilter {
 
   Configurable<std::string> trackSelections{"trackSelections", "globalTracks", "set track selections"};
   Configurable<float> trackPtMin{"trackPtMin", 0.15, "min track pT in filter"};
-  Configurable<float> trackPtTriggerThreshold{"trackPtThreshold", 25.0, "pT threshold for track trigger"};
+  Configurable<float> trackLowPtTriggerThreshold{"trackLowPtThreshold", 10.0, "low pT track trigger threshold"};
+  Configurable<float> trackHighPtTriggerThreshold{"trackHighPtThreshold", 25.0, "high pT track trigger threshold"};
   Configurable<float> jetPtLowThreshold{"jetPtLowThreshold", 30.0, "threshold for low pT jet trigger"};
   Configurable<float> jetPtHighThreshold{"jetPtHighThreshold", 50.0, "threshold for high pT jet trigger"};
 
@@ -107,7 +110,11 @@ struct jetFilter {
                 {{150, 0., +150., "charged jet #it{p}_{T} (GeV/#it{c})"},
                  {60, 0, TMath::TwoPi(), "#varphi"}});
 
-    spectra.add("ptphiTrackSelected_trackpttrigger", "#it{p}_{T} of selected high pT tracks vs #varphi", HistType::kTH2F,
+    spectra.add("ptphiTrackSelected_lowtrackpttrigger", "#it{p}_{T} of selected low pT tracks vs #varphi", HistType::kTH2F,
+                {{150, 0., +150., "track #it{p}_{T} (GeV/#it{c})"},
+                 {60, 0, TMath::TwoPi(), "#varphi"}});
+
+    spectra.add("ptphiTrackSelected_hightrackpttrigger", "#it{p}_{T} of selected high pT tracks vs #varphi", HistType::kTH2F,
                 {{150, 0., +150., "track #it{p}_{T} (GeV/#it{c})"},
                  {60, 0, TMath::TwoPi(), "#varphi"}});
 
@@ -119,7 +126,11 @@ struct jetFilter {
                 {{150, 0., +150., "charged jet #it{p}_{T} (GeV/#it{c})"},
                  {40, -1.0, 1.0, "#eta"}});
 
-    spectra.add("ptetaTrackSelected_trackpttrigger", "#it{p}_{T} of selected high pT tracks vs #eta", HistType::kTH2F,
+    spectra.add("ptetaTrackSelected_lowtrackpttrigger", "#it{p}_{T} of selected low pT tracks vs #eta", HistType::kTH2F,
+                {{150, 0., +150., "track #it{p}_{T} (GeV/#it{c})"},
+                 {40, -1.0, 1.0, "#eta"}});
+
+    spectra.add("ptetaTrackSelected_hightrackpttrigger", "#it{p}_{T} of selected high pT tracks vs #eta", HistType::kTH2F,
                 {{150, 0., +150., "track #it{p}_{T} (GeV/#it{c})"},
                  {40, -1.0, 1.0, "#eta"}});
 
@@ -134,12 +145,18 @@ struct jetFilter {
       jetRFidVolume.push_back(cfgEtaTPC - jetRadiiAxis.binEdges[ir]);
       jetIntR.push_back(TMath::Nint(100 * jetRadiiAxis.binEdges[ir]));
     }
+    spectra.add("hLeadingTrackPt", "Leading track pT in |#eta| < 0.9;",
+                {HistType::kTH1F, {axisTrackPt}});
+
     spectra.add("hEtaVsPtTracksInclusive", "#eta of tracks |#eta| < 0.9;",
                 {HistType::kTH2F, {axisTrackPt, axisEta}});
 
     spectra.add("hPhiVsPtTracksInclusive", "#varphi of tracks |#eta| < 0.9;",
                 {HistType::kTH2F, {axisTrackPt, axisPhi}});
 
+    spectra.add("hLeadingAKTJetR06Pt", "#it{p}_{T} of AKT R=0.6 charged jets in |#eta| < 0.9;",
+                {HistType::kTH1F, {axisJetPt}});
+
     spectra.add("hPtAKTJetsInclusive", "#it{p}_{T} of AKT charged jets in |#eta| < 0.9 - #it{R};",
                 {HistType::kTH2F, {jetRadiiAxis, axisJetPt}});
 
@@ -157,8 +174,9 @@ struct jetFilter {
 
     hProcessedEvents->GetXaxis()->SetBinLabel(kBinAllEvents + 1, "Processed events");
     hProcessedEvents->GetXaxis()->SetBinLabel(kBinJetChLowPt + 1, o2::aod::filtering::JetChLowPt::columnLabel());
-    hProcessedEvents->GetXaxis()->SetBinLabel(kBinJetChHighPt + 1, o2::aod::filtering::JetChLowPt::columnLabel());
-    hProcessedEvents->GetXaxis()->SetBinLabel(kBinTrackPt + 1, o2::aod::filtering::TrackHighPt::columnLabel());
+    hProcessedEvents->GetXaxis()->SetBinLabel(kBinJetChHighPt + 1, o2::aod::filtering::JetChHighPt::columnLabel());
+    hProcessedEvents->GetXaxis()->SetBinLabel(kBinTrackLowPt + 1, o2::aod::filtering::TrackLowPt::columnLabel());
+    hProcessedEvents->GetXaxis()->SetBinLabel(kBinTrackHighPt + 1, o2::aod::filtering::TrackHighPt::columnLabel());
   }
 
   template <bool withRho, typename T, typename U, typename D>
@@ -168,7 +186,7 @@ struct jetFilter {
     // collision process loop
     bool keepEvent[kTriggerObjects]{false};
     if (!collision.selection_bit(aod::evsel::kNoTimeFrameBorder)) {
-      tags(keepEvent[kJetChLowPt], keepEvent[kJetChHighPt], keepEvent[kTrackPt]);
+      tags(keepEvent[kJetChLowPt], keepEvent[kJetChHighPt], keepEvent[kTrackLowPt], keepEvent[kTrackHighPt]);
       return;
     }
 
@@ -204,6 +222,7 @@ struct jetFilter {
     for (const auto& jet : jets) { // jets are ordered by pT
       if (jet.r() != triggerJetR)
         continue;
+      spectra.fill(HIST("hLeadingAKTJetR06Pt"), jet.pt());
 
       if (jet.pt() >= jetPtLowThreshold) {
         spectra.fill(HIST("ptphiJetChSelected_lowptjettrigger"), jet.pt(), jet.phi()); // charged jet pT vs phi
@@ -234,20 +253,29 @@ struct jetFilter {
         leadingTrackEta = track.eta();
       }
     }
-    if (leadingTrackPt > trackPtTriggerThreshold) {
-      keepEvent[kTrackPt] = true;
-      spectra.fill(HIST("ptphiTrackSelected_trackpttrigger"), leadingTrackPt, leadingTrackPhi);
-      spectra.fill(HIST("ptetaTrackSelected_trackpttrigger"), leadingTrackPt, leadingTrackEta);
+    spectra.fill(HIST("hLeadingTrackPt"), leadingTrackPt);
+
+    if (leadingTrackPt > trackLowPtTriggerThreshold) {
+      keepEvent[kTrackLowPt] = true;
+      spectra.fill(HIST("ptphiTrackSelected_lowtrackpttrigger"), leadingTrackPt, leadingTrackPhi);
+      spectra.fill(HIST("ptetaTrackSelected_lowtrackpttrigger"), leadingTrackPt, leadingTrackEta);
+    }
+    if (leadingTrackPt > trackHighPtTriggerThreshold) {
+      keepEvent[kTrackHighPt] = true;
+      spectra.fill(HIST("ptphiTrackSelected_hightrackpttrigger"), leadingTrackPt, leadingTrackPhi);
+      spectra.fill(HIST("ptetaTrackSelected_hightrackpttrigger"), leadingTrackPt, leadingTrackEta);
     }
 
     if (keepEvent[kJetChLowPt])
       hProcessedEvents->Fill(static_cast<float>(kBinJetChLowPt) + 0.1f);
     if (keepEvent[kJetChHighPt])
       hProcessedEvents->Fill(static_cast<float>(kBinJetChHighPt) + 0.1f);
-    if (keepEvent[kTrackPt])
-      hProcessedEvents->Fill(static_cast<float>(kBinTrackPt) + 0.1f);
+    if (keepEvent[kTrackLowPt])
+      hProcessedEvents->Fill(static_cast<float>(kBinTrackLowPt) + 0.1f);
+    if (keepEvent[kTrackHighPt])
+      hProcessedEvents->Fill(static_cast<float>(kBinTrackHighPt) + 0.1f);
 
-    tags(keepEvent[kJetChLowPt], keepEvent[kJetChHighPt], keepEvent[kTrackPt]);
+    tags(keepEvent[kJetChLowPt], keepEvent[kJetChHighPt], keepEvent[kTrackLowPt], keepEvent[kTrackHighPt]);
   }
 
   void processWithoutRho(soa::Join<JetCollisions, aod::EvSels>::iterator const& collision, o2::aod::ChargedJets const& jets, soa::Filtered<JetTracks> const& tracks)

EventFiltering/filterTables.h

@@ -66,6 +66,7 @@ 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
+DECLARE_SOA_COLUMN(TrackLowPt, hasTrackLowPt, bool);   //! low-pT track
 DECLARE_SOA_COLUMN(TrackHighPt, hasTrackHighPt, bool); //! high-pT track
 
 // hf-jets
@@ -179,6 +180,7 @@ using CfFilter = CFFilters::iterator;
 DECLARE_SOA_TABLE(JetFilters, "AOD", "JetFilters", //!
                   filtering::JetChLowPt,
                   filtering::JetChHighPt,
+                  filtering::TrackLowPt,
                   filtering::TrackHighPt);
 
 using JetFilter = JetFilters::iterator;

PWGJE/Core/JetDerivedDataUtilities.h

@@ -90,7 +90,8 @@ enum JTrigSelCh {
   noChargedTigger = 0,
   chargedLow = 1,
   chargedHigh = 2,
-  trackPt = 3
+  trackLowPt = 3,
+  trackHighPt = 4
 };
 
 template <typename T>
@@ -110,9 +111,13 @@ int initialiseChargedTriggerSelection(std::string triggerSelection)
   if (triggerSelection == "chargedHigh") {
     return JTrigSelCh::chargedHigh;
   }
-  if (triggerSelection == "trackPt") {
-    return JTrigSelCh::trackPt;
+  if (triggerSelection == "trackLowPt") {
+    return JTrigSelCh::trackLowPt;
   }
+  if (triggerSelection == "trackHighPt") {
+    return JTrigSelCh::trackHighPt;
+  }
+
   return -1;
 }
 
@@ -127,9 +132,13 @@ uint8_t setChargedTriggerSelectionBit(T const& collision)
   if (collision.hasJetChHighPt()) {
     SETBIT(bit, JTrigSelCh::chargedHigh);
   }
+  if (collision.hasTrackLowPt()) {
+    SETBIT(bit, JTrigSelCh::trackLowPt);
+  }
   if (collision.hasTrackHighPt()) {
-    SETBIT(bit, JTrigSelCh::trackPt);
+    SETBIT(bit, JTrigSelCh::trackHighPt);
   }
+
   return bit;
 }
 

PWGJE/Tasks/ChJetTriggerQATask.cxx

@@ -73,7 +73,8 @@ struct ChJetTriggerQATask {
 
   Configurable<bool> bLowPtTrigger{"bLowPtTrigger", false, "charged jet low pT trigger selection"};
   Configurable<bool> bHighPtTrigger{"bHighPtTrigger", false, "charged jet high pT trigger selection"};
-  Configurable<bool> bTrackPtTrigger{"bTrackPtTrigger", false, "track pT trigger selection"};
+  Configurable<bool> bTrackLowPtTrigger{"bTrackLowPtTrigger", false, "track low pT trigger selection"};
+  Configurable<bool> bTrackHighPtTrigger{"bTrackHighPtTrigger", false, "track high pT trigger selection"};
 
   Configurable<bool> bAddSupplementHistosToOutput{"bAddAdditionalHistosToOutput", false, "add supplementary histos to the output"};
 
@@ -144,12 +145,19 @@ struct ChJetTriggerQATask {
       return;
     }
 
-    if ((bLowPtTrigger && jetderiveddatautilities::selectChargedTrigger(collision, jetderiveddatautilities::JTrigSelCh::chargedLow)) || (bHighPtTrigger && jetderiveddatautilities::selectChargedTrigger(collision, jetderiveddatautilities::JTrigSelCh::chargedHigh)) || (bTrackPtTrigger && jetderiveddatautilities::selectChargedTrigger(collision, jetderiveddatautilities::JTrigSelCh::trackPt)) || ((!bLowPtTrigger) && (!bHighPtTrigger) && (!bTrackPtTrigger))) {
+    bool bLowPtJet = (bLowPtTrigger && jetderiveddatautilities::selectChargedTrigger(collision, jetderiveddatautilities::JTrigSelCh::chargedLow));
+    bool bHighPtJet = (bHighPtTrigger && jetderiveddatautilities::selectChargedTrigger(collision, jetderiveddatautilities::JTrigSelCh::chargedHigh));
+    bool bLowPtTrack = (bTrackLowPtTrigger && jetderiveddatautilities::selectChargedTrigger(collision, jetderiveddatautilities::JTrigSelCh::trackLowPt));
+    bool bHighPtTrack = (bTrackHighPtTrigger && jetderiveddatautilities::selectChargedTrigger(collision, jetderiveddatautilities::JTrigSelCh::trackHighPt));
+    bool bMinimumBias = ((!bLowPtTrigger) && (!bHighPtTrigger) && (!bTrackLowPtTrigger) && (!bTrackHighPtTrigger));
+
+    if (bLowPtJet || bHighPtJet || bLowPtTrack || bHighPtTrack || bMinimumBias) {
       // bLowPtTrigger=1  and bHighPtTrigger=0 --> fill histos with low trigger only
       // bLowPtTrigger=0  and bHighPtTrigger=1 --> fill histos with high trigger only
       // bLowPtTrigger=1  and bHighPtTrigger=1 --> fill histos with mixture of low and high trigger
-      // bTrackPtTrigger=1 --> fill histos for high pt track trigger
-      // bLowPtTrigger=0  and bHighPtTrigger=0 and bTrackPtTrigger=0 --> fill histos with minimum bias ie. ignore trigger decision
+      // bTrackLowPtTrigger=1 --> fill histos for low pt track trigger
+      // bTrackHighPtTrigger=1 --> fill histos for high pt track trigger
+      // bLowPtTrigger=0 and bHighPtTrigger=0 and bTrackLowPtTrigger=0 and bTrackHighPtTrigger=0 --> fill histos with minimum bias ie. ignore trigger decision
 
       float leadingJetPt = -1.0;
       float leadingJetEta = -2.0;

PWGJE/Tasks/triggerCorrelations.cxx

@@ -55,7 +55,7 @@ struct TriggerCorrelationsTask {
 
   void init(o2::framework::InitContext&)
   {
-    std::vector<std::string> trigSelChLabels = {"chargedLow", "chargedHigh", "trackPt"};
+    std::vector<std::string> trigSelChLabels = {"chargedLow", "chargedHigh", "trackLowPt", "trackHighPt"};
     std::vector<std::string> trigSelFullLabels = {"fullHigh", "fullLow", "neutralHigh", "neutralLow", "gammaVeryHighEMCAL", "gammaHighEMCAL", "gammaLowEMCAL", "gammaVeryLowEMCAL", "gammaVeryHighDCAL", "gammaHighDCAL", "gammaLowDCAL", "gammaVeryLowDCAL"};
     std::vector<std::string> trigSelChHFLabels = {"chargedD0Low", "chargedD0High", "chargedLcLow", "chargedLcHigh"};
     nChTrigs = trigSelChLabels.size();