added flag for INEL > 0 selection

Author: Sawan Sawan

PWGLF/Tasks/Resonances/kstarqa.cxx

@@ -482,9 +482,9 @@ struct Kstarqa {
   Filter acceptanceFilter = (nabs(aod::track::eta) < cfgCutEta && nabs(aod::track::pt) > cfgCutPT);
   Filter fDCAcutFilter = (nabs(aod::track::dcaXY) < cfgCutDCAxy) && (nabs(aod::track::dcaZ) < cfgCutDCAz);
 
-  using EventCandidates = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::FT0Mults, aod::MultZeqs, aod::CentFT0Ms, aod::CentFT0As, aod::CentFT0Cs>>;
+  using EventCandidates = soa::Filtered<soa::Join<aod::Collisions, aod::EvSels, aod::FT0Mults, aod::MultZeqs, aod::CentFT0Ms, aod::CentFT0As, aod::CentFT0Cs, aod::PVMults>>;
   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, aod::FT0Mults>;
+  using EventCandidatesMC = soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels, aod::CentFT0Ms, aod::CentFT0As, aod::CentFT0Cs, aod::FT0Mults, aod::PVMults>;
 
   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>>;
 
@@ -867,7 +867,7 @@ struct Kstarqa {
       if (!collision.sel8()) {
         continue;
       }
-      if (isINELgt0 && !mcCollision.isInelGt0()) {
+      if (isINELgt0 && !collision.isInelGt0()) {
         continue;
       }
       if (ispileupGoodvtxCut && (!collision.selection_bit(aod::evsel::kNoSameBunchPileup) || !collision.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV))) {
@@ -931,41 +931,41 @@ struct Kstarqa {
   }
   PROCESS_SWITCH(Kstarqa, processGen, "Process Generated", false);
 
-  void processEvtLossSigLossMC(aod::McCollisions::iterator const&, aod::McParticles const& mcParticles, const soa::SmallGroups<EventCandidatesMC>& recCollisions)
-  {
-
-    bool isSel = false;
-    // auto multiplicity1 = -999.;
-    for (const auto& RecCollision : recCollisions) {
-      if (!selectionEvent(RecCollision))
-        continue;
-
-      // if (cSelectMultEstimator == 0) {
-      //   multiplicity1 = RecCollision.centFT0M();
-      // } else if (cSelectMultEstimator == 1) {
-      //   multiplicity1 = RecCollision.centFT0A();
-      // } else if (cSelectMultEstimator == 2) {
-      //   multiplicity1 = RecCollision.centFT0C();
-      // } else {
-      //   multiplicity1 = RecCollision.centFT0M();
-      // }
+  // void processEvtLossSigLossMC(aod::McCollisions::iterator const&, aod::McParticles const& mcParticles, const soa::SmallGroups<EventCandidatesMC>& recCollisions)
+  // {
 
-      isSel = true;
-    }
+  //   bool isSel = false;
+  //   // auto multiplicity1 = -999.;
+  //   for (const auto& RecCollision : recCollisions) {
+  //     if (!selectionEvent(RecCollision))
+  //       continue;
+
+  //     // if (cSelectMultEstimator == 0) {
+  //     //   multiplicity1 = RecCollision.centFT0M();
+  //     // } else if (cSelectMultEstimator == 1) {
+  //     //   multiplicity1 = RecCollision.centFT0A();
+  //     // } else if (cSelectMultEstimator == 2) {
+  //     //   multiplicity1 = RecCollision.centFT0C();
+  //     // } else {
+  //     //   multiplicity1 = RecCollision.centFT0M();
+  //     // }
+
+  //     isSel = true;
+  //   }
 
-    // Generated MC
-    for (const auto& mcPart : mcParticles) {
-      if (std::abs(mcPart.y()) >= 0.5 || std::abs(mcPart.pdgCode()) != 313)
-        continue;
+  //   // Generated MC
+  //   for (const auto& mcPart : mcParticles) {
+  //     if (std::abs(mcPart.y()) >= 0.5 || std::abs(mcPart.pdgCode()) != 313)
+  //       continue;
 
-      // signal loss estimation
-      hInvMass.fill(HIST("kstargenBeforeEvtSel"), mcPart.pt());
-      if (isSel) {
-        hInvMass.fill(HIST("kstargenAfterEvtSel"), mcPart.pt());
-      }
-    } // end loop on gen particles
-  }
-  PROCESS_SWITCH(Kstarqa, processEvtLossSigLossMC, "Process Signal Loss, Event Loss", false);
+  //     // signal loss estimation
+  //     hInvMass.fill(HIST("kstargenBeforeEvtSel"), mcPart.pt());
+  //     if (isSel) {
+  //       hInvMass.fill(HIST("kstargenAfterEvtSel"), mcPart.pt());
+  //     }
+  //   } // end loop on gen particles
+  // }
+  // PROCESS_SWITCH(Kstarqa, processEvtLossSigLossMC, "Process Signal Loss, Event Loss", false);
 
   void processRec(EventCandidatesMC::iterator const& collision, TrackCandidatesMC const& tracks, aod::McParticles const&, aod::McCollisions const& /*mcCollisions*/)
   {