[PWGJE] EMCal: event selection for emcTmMonitor, track matching optio… (#18606)

Author: mhemmer-cern

PWGJE/TableProducer/emcalCorrectionTask.cxx

@@ -133,7 +133,6 @@ struct EmcalCorrectionTask {
   Configurable<bool> useWeightExotic{"useWeightExotic", false, "States if weights should be used for exotic cell cut"};
   Configurable<bool> isMC{"isMC", false, "States if run over MC"};
   Configurable<bool> applyCellTimeCorrection{"applyCellTimeCorrection", true, "apply a correction to the cell time for data and MC: Shift both average cell times to 0 and smear MC time distribution to fit data better. For MC requires isMC to be true"};
-  Configurable<float> trackMinPt{"trackMinPt", 0.3, "Minimum pT for tracks to perform track matching, to reduce computing time. Tracks below a certain pT will be loopers anyway."};
   Configurable<bool> fillQA{"fillQA", false, "Switch to turn on QA histograms."};
   Configurable<bool> useCCDBAlignment{"useCCDBAlignment", false, "EXPERTS ONLY! Switch to use the alignment object stored in CCDB instead of using the default alignment from the global geometry object."};
   Configurable<bool> applyTempCalib{"applyTempCalib", false, "Switch to turn on Temperature calibration."};
@@ -147,6 +146,12 @@ struct EmcalCorrectionTask {
   Configurable<bool> storePerDFInfo{"storePerDFInfo", false, "store addition information per DF."};
   ConfigurableAxis thConfigAxisClusters{"thConfigAxisClusters", {1000, 0.5f, 1000.5f}, ""};
   ConfigurableAxis thConfigAxisCells{"thConfigAxisCells", {1000, 0.5f, 1000.5f}, ""};
+  // Track matching config
+  struct : ConfigurableGroup {
+    std::string prefix = "trackmatching";
+    Configurable<float> trackMinPt{"trackMinPt", 0.3, "Minimum pT (GeV) for tracks to perform track matching, to reduce computing time. Tracks below a certain pT will be loopers anyway."};
+    Configurable<bool> useAllTracks{"useAllTracks", false, "Flag to use all tracks instead of just global tracks."};
+  } trackmatching;
   // cross talk emulation configs
   EmcCrossTalkConf emcCrossTalkConf;
 
@@ -1494,7 +1499,7 @@ struct EmcalCorrectionTask {
       if (trackEtaEmcal < TrackNotOnEMCal || trackPhiEmcal < TrackNotOnEMCal) {
         continue;
       }
-      if (trackMinPt > 0 && track.pt() < trackMinPt) {
+      if (trackmatching.trackMinPt > 0 && track.pt() < trackmatching.trackMinPt) {
         continue;
       }
       trackPhi.emplace_back(RecoDecay::constrainAngle(trackPhiEmcal));
@@ -1508,15 +1513,15 @@ struct EmcalCorrectionTask {
   void fillTrackInfo(Tracks const& tracks, std::vector<float>& trackPhi, std::vector<float>& trackEta, std::vector<int64_t>& trackGlobalIndex)
   {
     for (const auto& track : tracks) {
-      if (!track.isGlobalTrack()) { // only global tracks
+      if (!trackmatching.useAllTracks && !track.isGlobalTrack()) { // only global tracks depending on trackmatching.useAllTracks flag
         continue;
       }
       // Tracks that do not point to the EMCal/DCal/PHOS get default values of -999
       // This way we can cut out tracks that do not point to the EMCal+DCal
       if (track.trackEtaEmcal() < TrackNotOnEMCal || track.trackPhiEmcal() < TrackNotOnEMCal) {
         continue;
       }
-      if (trackMinPt > 0 && track.pt() < trackMinPt) {
+      if (trackmatching.trackMinPt > 0 && track.pt() < trackmatching.trackMinPt) {
         continue;
       }
       trackPhi.emplace_back(RecoDecay::constrainAngle(track.trackPhiEmcal()));