[PWGLF] Add new ITS based event selection for resonance flow analysis and add globaltrack selection criteria for f1 reconstruction (#9410)

Author: skundu692

PWGLF/TableProducer/Resonances/f1protonreducedtable.cxx

@@ -248,6 +248,15 @@ struct f1protonreducedtable {
     return true;
   }
 
+  template <typename T>
+  bool selectionGlobalTrack(const T& candidate)
+  {
+    if (!(candidate.isGlobalTrack() && candidate.isPVContributor())) {
+      return false;
+    }
+    return true;
+  }
+
   template <typename T>
   double updatePID(T const& track, double bgScaling, std::vector<double> BB)
   {
@@ -586,8 +595,12 @@ struct f1protonreducedtable {
       if (zorroSelected) {
         hProcessedEvents->Fill(1.5);
         for (auto& track : tracks) {
-          if (!isSelectedTrack(track))
+          if (!isSelectedTrack(track)) {
             continue;
+          }
+          if (!selectionGlobalTrack(track)) {
+            continue;
+          }
           qaRegistry.fill(HIST("hDCAxy"), track.dcaXY());
           qaRegistry.fill(HIST("hDCAz"), track.dcaZ());
           qaRegistry.fill(HIST("hEta"), track.eta());

PWGLF/Tasks/Resonances/highmasslambda.cxx

@@ -77,6 +77,7 @@ struct highmasslambda {
   // fill output
   Configurable<bool> fillQA{"fillQA", false, "fillQA"};
   Configurable<bool> useSP{"useSP", false, "useSP"};
+  // Configurable<bool> additionalEvselITS{"additionalEvselITS", true, "additionalEvselITS"};
   Configurable<bool> useSignDCAV0{"useSignDCAV0", true, "useSignDCAV0"};
   Configurable<bool> fillDefault{"fillDefault", false, "fill Occupancy"};
   Configurable<int> cfgOccupancyCut{"cfgOccupancyCut", 2500, "Occupancy cut"};
@@ -560,7 +561,9 @@ struct highmasslambda {
     if (!collision.sel8() || !collision.triggereventep() || !collision.selection_bit(aod::evsel::kNoTimeFrameBorder) || !collision.selection_bit(aod::evsel::kNoITSROFrameBorder) || !collision.selection_bit(aod::evsel::kNoSameBunchPileup) || !collision.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV) || !collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard)) {
       return;
     }
-
+    if (!collision.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) {
+      return;
+    }
     auto centrality = collision.centFT0C();
     auto multTPC = collision.multNTracksPV();
     histos.fill(HIST("hFTOCvsTPCNoCut"), centrality, multTPC);
@@ -765,6 +768,12 @@ struct highmasslambda {
       if (collision1.bcId() == collision2.bcId()) {
         continue;
       }
+      if (!collision1.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) {
+        continue;
+      }
+      if (!collision2.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) {
+        continue;
+      }
       auto centrality = collision1.centFT0C();
       auto psiFT0C = collision1.psiFT0C();
       auto QFT0C = collision1.qFT0C();

PWGLF/Tasks/Resonances/phipbpb.cxx

@@ -80,6 +80,7 @@ struct phipbpb {
   Configurable<int> cfgCutOccupancy{"cfgCutOccupancy", 3000, "Occupancy cut"};
   // track
   Configurable<bool> additionalEvsel{"additionalEvsel", false, "Additional event selcection"};
+  Configurable<bool> additionalEvselITS{"additionalEvselITS", true, "Additional event selcection for ITS"};
   Configurable<bool> removefaketrak{"removefaketrack", true, "Remove fake track from momentum difference"};
   Configurable<float> ConfFakeKaonCut{"ConfFakeKaonCut", 0.1, "Cut based on track from momentum difference"};
   Configurable<bool> useGlobalTrack{"useGlobalTrack", true, "use Global track"};
@@ -387,6 +388,9 @@ struct phipbpb {
     if (additionalEvsel && !eventSelected(collision, centrality)) {
       return;
     }
+    if (additionalEvselITS && !collision.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) {
+      return;
+    }
     histos.fill(HIST("hFTOCvsTPCSelected"), centrality, multTPC);
     histos.fill(HIST("hPsiFT0C"), centrality, occupancy, psiFT0C);
     histos.fill(HIST("hPsiFT0A"), centrality, occupancy, psiFT0A);
@@ -525,7 +529,12 @@ struct phipbpb {
       if (additionalEvsel && !eventSelected(collision2, centrality2)) {
         continue;
       }
-
+      if (additionalEvselITS && !collision1.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) {
+        continue;
+      }
+      if (additionalEvselITS && !collision2.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) {
+        continue;
+      }
       for (auto& [track1, track2] : o2::soa::combinations(o2::soa::CombinationsFullIndexPolicy(tracks1, tracks2))) {
         if (track1.sign() * track2.sign() > 0) {
           continue;