inel check added

Author: Sawan Sawan

PWGLF/Tasks/Resonances/kstarqa.cxx

@@ -78,6 +78,7 @@ struct Kstarqa {
   HistogramRegistry hOthers{"hOthers", {}, OutputObjHandlingPolicy::AnalysisObject, true, true};
 
   // Confugrable for QA histograms
+  Configurable<bool> isINELgt0{"isINELgt0", true, "INEL>0 selection"};
   Configurable<bool> calcLikeSign{"calcLikeSign", true, "Calculate Like Sign"};
   Configurable<bool> calcRotational{"calcRotational", false, "Calculate Rotational"};
   Configurable<bool> cQAplots{"cQAplots", true, "cQAplots"};
@@ -236,6 +237,9 @@ struct Kstarqa {
   template <typename Coll>
   bool selectionEvent(const Coll& collision)
   {
+    if (isINELgt0 && !collision.isInelGt0()) {
+      return false;
+    }
     if (std::abs(collision.posZ()) > cutzvertex)
       return false;
     if (!collision.sel8())
@@ -480,7 +484,7 @@ struct Kstarqa {
 
   using EventCandidates = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::FT0Mults, aod::MultZeqs, aod::CentFT0Ms, aod::CentFT0As, aod::CentFT0Cs>>;
   using TrackCandidates = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::TrackSelection, aod::pidTPCFullPi, aod::pidTOFFullPi, aod::pidTPCFullKa, aod::pidTOFFullKa, aod::pidTOFbeta>>;
-  using EventCandidatesMC = soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels, aod::CentFT0Ms, aod::CentFT0As, aod::CentFT0Cs>;
+  using EventCandidatesMC = soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels, aod::CentFT0Ms, aod::CentFT0As, aod::CentFT0Cs, aod::FT0Mults>;
 
   using TrackCandidatesMC = soa::Filtered<soa::Join<aod::Tracks, aod::TracksExtra, aod::TracksDCA, aod::TrackSelection, aod::pidTPCFullKa, aod::pidTOFFullKa, aod::pidTPCFullPi, aod::pidTOFFullPi, aod::McTrackLabels, aod::pidTOFbeta>>;
 
@@ -863,6 +867,9 @@ struct Kstarqa {
       if (!collision.sel8()) {
         continue;
       }
+      if (isINELgt0 && !mcCollision.isInelGt0()) {
+        continue;
+      }
       if (ispileupGoodvtxCut && (!collision.selection_bit(aod::evsel::kNoSameBunchPileup) || !collision.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV))) {
         continue;
       }
@@ -970,6 +977,10 @@ struct Kstarqa {
     }
     rEventSelection.fill(HIST("events_checkrec"), 1.5);
 
+    if (isINELgt0 && !collision.isInelGt0()) {
+      return;
+    }
+
     // if (std::abs(collision.mcCollision().posZ()) > cutzvertex || !collision.sel8()) {
     if (std::abs(collision.mcCollision().posZ()) > cutzvertex) {
       return;