[PWGMM] dndeta-hi: add SD/DD/ND classification and fix sliceBy on filtered tables

- Add kDD, kSD, kND enum values for diffractive event classification
- Join McCollisions with HepMCXSections to access processId() in processMCCounting
- Classify events by Pythia8 process code: 101=ND, 103/104=SD, 105/106=DD
- Replace sliceByCached with sliceBy+Preslice for soa::Filtered tables (FiTracks, FiLTracks)
  to fix runtime error: Disabled cache /fIndexCollisions is requested

Author: hyungjun-lee

PWGMM/Mult/Tasks/dndeta-hi.cxx

@@ -75,6 +75,9 @@ enum {
   kDATA = 1,
   kINEL,
   kINELg0,
+  kDD,
+  kSD,
+  kND,
   kECend
 };
 enum {
@@ -173,9 +176,11 @@ static constexpr TrackSelectionFlags::flagtype trackSelectionTPC =
 static constexpr TrackSelectionFlags::flagtype trackSelectionDCA =
   TrackSelectionFlags::kDCAz | TrackSelectionFlags::kDCAxy;
 
+
 struct MultiplicityCounter {
   SliceCache cache;
-  Preslice<aod::McParticles> perMCCol = aod::mcparticle::mcCollisionId;
+  Preslice<Particles> perMCCol = aod::mcparticle::mcCollisionId;
+  Preslice<FiTracks> perCol = aod::track::collisionId;
 
   Service<o2::framework::O2DatabasePDG> pdg;
 
@@ -318,7 +323,7 @@ struct MultiplicityCounter {
       Bool_1d btrigc(kTrigend, false);
       registry.fill(HIST("Selection"), 1.);
       auto z = collision.posZ();
-      auto pertracks = tSample3->sliceByCached(aod::track::collisionId, collision.globalIndex(), cache);
+      auto pertracks = tSample3->sliceBy(perCol, collision.globalIndex());
       auto Ntrk = 0;
 
       // if (collision.selection_bit(aod::evsel::kIsTriggerTVX) && collision.selection_bit(aod::evsel::kNoTimeFrameBorder) && collision.selection_bit(aod::evsel::kNoITSROFrameBorder)) {
@@ -388,13 +393,23 @@ struct MultiplicityCounter {
 
   PresliceUnsorted<soa::Join<MyCollisions, aod::McCollisionLabels>> perMcCol = o2::aod::mccollisionlabel::mcCollisionId;
   Preslice<aod::McParticles> perMCColparticles = aod::mcparticle::mcCollisionId;
+  Preslice<FiLTracks> perColFiLTracks = aod::track::collisionId;
+  using MCex = soa::Join<aod::McCollisions, aod::HepMCXSections>;
   void processMCCounting(
-    aod::McCollisions const& mcCollisions, soa::Join<MyCollisionsCent, aod::McCollisionLabels> const& collisions, Particles const& mcParticles,
+    MCex const& mcCollisions, soa::Join<MyCollisionsCent, aod::McCollisionLabels> const& collisions, Particles const& mcParticles,
     FiLTracks const& tracks)
   {
     for (auto& mcCollision : mcCollisions) {
       Bool_1d bevtc(kECend, false);
       bevtc[kINEL] = true;
+      auto procId = mcCollision.processId();
+      if (procId == 101) {
+        bevtc[kND] = true;
+      } else if (procId == 103 || procId == 104) {
+        bevtc[kSD] = true;
+      } else if (procId == 105 || procId == 106) {
+        bevtc[kDD] = true;
+      }
       registry.fill(HIST("Selection"), 1.);
 
       auto mcz = mcCollision.posZ();
@@ -565,7 +580,7 @@ struct MultiplicityCounter {
         }
 
         auto Ntrk_rec = 0;
-        auto trackspart = tracks.sliceByCached(aod::track::collisionId, collision.globalIndex(), cache);
+        auto trackspart = tracks.sliceBy(perColFiLTracks, collision.globalIndex());
         for (auto& track : trackspart) {
           if (std::abs(track.eta()) < 1) {
             Ntrk_rec++;
@@ -677,7 +692,7 @@ struct MultiplicityCounter {
     }
 
     auto mcCollision = collision.mcCollision();
-    auto particlesPerCol = particles.sliceByCached(aod::mcparticle::mcCollisionId, mcCollision.globalIndex(), cache);
+    auto particlesPerCol = particles.sliceBy(perMCCol, mcCollision.globalIndex());
   }
 
   PROCESS_SWITCH(MultiplicityCounter, processTrackEfficiencyGeneral, "MC Count tracks", false);