diff --git a/PWGCF/EbyEFluctuations/Tasks/radialFlowDecorr.cxx b/PWGCF/EbyEFluctuations/Tasks/radialFlowDecorr.cxx
index bc7c5078594..36fd6cd2254 100644
--- a/PWGCF/EbyEFluctuations/Tasks/radialFlowDecorr.cxx
+++ b/PWGCF/EbyEFluctuations/Tasks/radialFlowDecorr.cxx
@@ -550,10 +550,10 @@ struct RadialFlowDecorr {
     if (cfgApplyLinPupCut) {
       if (trksize > (cfgLinPupParam0 + cfgLinPupParam1 * multPV))
         return false;
-        histos.fill(HIST("hEvtCount"), 7.5);
-     if (trksize < (cfgLinPupParam2 + cfgLinPupParam3 * multPV))
+      histos.fill(HIST("hEvtCount"), 7.5);
+      if (trksize < (cfgLinPupParam2 + cfgLinPupParam3 * multPV))
         return false;
-        histos.fill(HIST("hEvtCount"), 8.5);
+      histos.fill(HIST("hEvtCount"), 8.5);
     }
     return true;
   }
@@ -800,7 +800,7 @@ struct RadialFlowDecorr {
     aod::CentFT0Cs, aod::CentFT0Ms, aod::CentFDDMs, aod::CentFV0As,
     aod::CentNGlobals, aod::McCollisionLabels>;
 
-    PresliceUnsorted<MyRun3MCCollisions> colPerMcCollision = aod::mccollisionlabel::mcCollisionId;
+  PresliceUnsorted<MyRun3MCCollisions> colPerMcCollision = aod::mccollisionlabel::mcCollisionId;
 
   void declareCommonQA()
   {
@@ -810,8 +810,8 @@ struct RadialFlowDecorr {
     histos.add("Hist2D_globalTracks_PVTracks", ";N_{global};N_{PV}", kTH2F, {{nChAxis}, {nChAxis}});
     histos.add("Hist2D_cent_nch", ";N_{PV};cent (%)", kTH2F, {{nChAxis}, {centAxis1Per}});
 
-    histos.add("Hist2D_globalTracks_cent", "cent (%);N_{global}", kTH2F, {{centAxis1Per},{nChAxis}});
-    histos.add("Hist2D_PVTracks_cent", "cent (%);N_{PV}", kTH2F, {{centAxis1Per},{nChAxis}});
+    histos.add("Hist2D_globalTracks_cent", "cent (%);N_{global}", kTH2F, {{centAxis1Per}, {nChAxis}});
+    histos.add("Hist2D_PVTracks_cent", "cent (%);N_{PV}", kTH2F, {{centAxis1Per}, {nChAxis}});
 
     histos.add("hP", ";p (GeV/c)", kTH1F, {{KNbinsPt, KPMin, KPMax}});
     histos.add("hPt", ";p_{T} (GeV/c)", kTH1F, {{KNbinsPt, KPtMin, KPtMax}});
@@ -1436,375 +1436,375 @@ struct RadialFlowDecorr {
 
   void processMCGetMeanNsig(MyRun3MCCollisions::iterator const& mcCollision, FilteredTCs const& mcTracks)
   {
-        histos.fill(HIST("hVtxZ"), mcCollision.posZ());
-        if (!mcCollision.has_mcCollision() || !isEventSelected(mcCollision))
-          return;
-        float cent = getCentrality(mcCollision);
-        if (cent > KCentMax)
-          return;
-        float multPV = mcCollision.multNTracksPV();
-
-        histos.fill(HIST("hVtxZ_after_sel"), mcCollision.posZ());
-        histos.fill(HIST("hCentrality"), cent);
-        histos.fill(HIST("Hist2D_globalTracks_PVTracks"), multPV, mcTracks.size());
-        histos.fill(HIST("Hist2D_cent_nch"), mcTracks.size(), cent);
-        histos.fill(HIST("Hist2D_globalTracks_cent"), cent, mcTracks.size());
-        histos.fill(HIST("Hist2D_PVTracks_cent"), cent, multPV);
+    histos.fill(HIST("hVtxZ"), mcCollision.posZ());
+    if (!mcCollision.has_mcCollision() || !isEventSelected(mcCollision))
+      return;
+    float cent = getCentrality(mcCollision);
+    if (cent > KCentMax)
+      return;
+    float multPV = mcCollision.multNTracksPV();
 
-        for (const auto& track : mcTracks) {
-          if (!isTrackSelected(track))
-            continue;
-          fillNSigmaBefCut(track, cent);
-        }
-      }
+    histos.fill(HIST("hVtxZ_after_sel"), mcCollision.posZ());
+    histos.fill(HIST("hCentrality"), cent);
+    histos.fill(HIST("Hist2D_globalTracks_PVTracks"), multPV, mcTracks.size());
+    histos.fill(HIST("Hist2D_cent_nch"), mcTracks.size(), cent);
+    histos.fill(HIST("Hist2D_globalTracks_cent"), cent, mcTracks.size());
+    histos.fill(HIST("Hist2D_PVTracks_cent"), cent, multPV);
+
+    for (const auto& track : mcTracks) {
+      if (!isTrackSelected(track))
+        continue;
+      fillNSigmaBefCut(track, cent);
+    }
+  }
   PROCESS_SWITCH(RadialFlowDecorr, processMCGetMeanNsig, "process MC to calculate Mean values of nSig Plots", cfgRunMCGetNSig);
 
   void processGetEffHists(MyRun3MCCollisions::iterator const& mcCollision, FilteredTCs const& mcTracks, aod::McParticles const& mcParticles)
   {
-        histos.fill(HIST("hVtxZ"), mcCollision.posZ());
-        if (!mcCollision.has_mcCollision() || !isEventSelected(mcCollision))
-          return;
-        float cent = getCentrality(mcCollision);
-        if (cent > KCentMax)
-          return;
-        float multPV = mcCollision.multNTracksPV();
-        float vz = mcCollision.posZ();
-        if (!isPassAddPileup(multPV, mcTracks.size(), cent))
-          return;
-        histos.fill(HIST("hVtxZ_after_sel"), mcCollision.posZ());
-        histos.fill(HIST("hCentrality"), cent);
-
-        histos.fill(HIST("Hist2D_globalTracks_PVTracks"), multPV, mcTracks.size());
-        histos.fill(HIST("Hist2D_cent_nch"), mcTracks.size(), cent);
-        histos.fill(HIST("Hist2D_globalTracks_cent"), cent, mcTracks.size());
-        histos.fill(HIST("Hist2D_PVTracks_cent"), cent, multPV);
-
-        for (const auto& particle : mcParticles) {
-          if (!isParticleSelected(particle) || !particle.isPhysicalPrimary())
-            continue;
+    histos.fill(HIST("hVtxZ"), mcCollision.posZ());
+    if (!mcCollision.has_mcCollision() || !isEventSelected(mcCollision))
+      return;
+    float cent = getCentrality(mcCollision);
+    if (cent > KCentMax)
+      return;
+    float multPV = mcCollision.multNTracksPV();
+    float vz = mcCollision.posZ();
+    if (!isPassAddPileup(multPV, mcTracks.size(), cent))
+      return;
+    histos.fill(HIST("hVtxZ_after_sel"), mcCollision.posZ());
+    histos.fill(HIST("hCentrality"), cent);
 
-          const int pdg = particle.pdgCode();
-          const int absPdg = std::abs(pdg);
-          float pt = particle.pt(), eta = particle.eta();
+    histos.fill(HIST("Hist2D_globalTracks_PVTracks"), multPV, mcTracks.size());
+    histos.fill(HIST("Hist2D_cent_nch"), mcTracks.size(), cent);
+    histos.fill(HIST("Hist2D_globalTracks_cent"), cent, mcTracks.size());
+    histos.fill(HIST("Hist2D_PVTracks_cent"), cent, multPV);
 
-          bool isSpecies[KNsp] = {
-            true,              // kInclusiveIdx
-            pdg == -KPiPlus,   // kPiMinusIdx
-            pdg == KPiPlus,    // kPiPlusIdx
-            absPdg == KPiPlus, // kPiAllIdx
-            pdg == -KKPlus,    // kKaMinusIdx
-            pdg == KKPlus,     // kKaPlusIdx
-            absPdg == KKPlus,  // kKaAllIdx
-            pdg == -KProton,   // kAntiPrIdx
-            pdg == KProton,    // kPrIdx
-            absPdg == KProton  // kPrAllIdx
-          };
+    for (const auto& particle : mcParticles) {
+      if (!isParticleSelected(particle) || !particle.isPhysicalPrimary())
+        continue;
 
-          histos.fill(HIST("h3_AllPrimary"), multPV, pt, eta);
-          if (isSpecies[kPiMinusIdx])
-            histos.fill(HIST("h3_AllPrimary_PiMinus"), multPV, pt, eta);
-          else if (isSpecies[kPiPlusIdx])
-            histos.fill(HIST("h3_AllPrimary_PiPlus"), multPV, pt, eta);
-          if (isSpecies[kPiAllIdx])
-            histos.fill(HIST("h3_AllPrimary_PiAll"), multPV, pt, eta);
-
-          if (isSpecies[kKaMinusIdx])
-            histos.fill(HIST("h3_AllPrimary_KaMinus"), multPV, pt, eta);
-          else if (isSpecies[kKaPlusIdx])
-            histos.fill(HIST("h3_AllPrimary_KaPlus"), multPV, pt, eta);
-          if (isSpecies[kKaAllIdx])
-            histos.fill(HIST("h3_AllPrimary_KaAll"), multPV, pt, eta);
-
-          if (isSpecies[kAntiPrIdx])
-            histos.fill(HIST("h3_AllPrimary_AntiPr"), multPV, pt, eta);
-          else if (isSpecies[kPrIdx])
-            histos.fill(HIST("h3_AllPrimary_Pr"), multPV, pt, eta);
-          if (isSpecies[kPrAllIdx])
-            histos.fill(HIST("h3_AllPrimary_PrAll"), multPV, pt, eta);
-        }
+      const int pdg = particle.pdgCode();
+      const int absPdg = std::abs(pdg);
+      float pt = particle.pt(), eta = particle.eta();
 
-        for (const auto& track : mcTracks) {
-          if (!isTrackSelected(track))
-            continue;
+      bool isSpecies[KNsp] = {
+        true,              // kInclusiveIdx
+        pdg == -KPiPlus,   // kPiMinusIdx
+        pdg == KPiPlus,    // kPiPlusIdx
+        absPdg == KPiPlus, // kPiAllIdx
+        pdg == -KKPlus,    // kKaMinusIdx
+        pdg == KKPlus,     // kKaPlusIdx
+        absPdg == KKPlus,  // kKaAllIdx
+        pdg == -KProton,   // kAntiPrIdx
+        pdg == KProton,    // kPrIdx
+        absPdg == KProton  // kPrAllIdx
+      };
 
-          float pt = track.pt(), eta = track.eta();
-          auto sign = track.sign();
-          fillNSigmaBefCut(track, cent);
+      histos.fill(HIST("h3_AllPrimary"), multPV, pt, eta);
+      if (isSpecies[kPiMinusIdx])
+        histos.fill(HIST("h3_AllPrimary_PiMinus"), multPV, pt, eta);
+      else if (isSpecies[kPiPlusIdx])
+        histos.fill(HIST("h3_AllPrimary_PiPlus"), multPV, pt, eta);
+      if (isSpecies[kPiAllIdx])
+        histos.fill(HIST("h3_AllPrimary_PiAll"), multPV, pt, eta);
+
+      if (isSpecies[kKaMinusIdx])
+        histos.fill(HIST("h3_AllPrimary_KaMinus"), multPV, pt, eta);
+      else if (isSpecies[kKaPlusIdx])
+        histos.fill(HIST("h3_AllPrimary_KaPlus"), multPV, pt, eta);
+      if (isSpecies[kKaAllIdx])
+        histos.fill(HIST("h3_AllPrimary_KaAll"), multPV, pt, eta);
+
+      if (isSpecies[kAntiPrIdx])
+        histos.fill(HIST("h3_AllPrimary_AntiPr"), multPV, pt, eta);
+      else if (isSpecies[kPrIdx])
+        histos.fill(HIST("h3_AllPrimary_Pr"), multPV, pt, eta);
+      if (isSpecies[kPrAllIdx])
+        histos.fill(HIST("h3_AllPrimary_PrAll"), multPV, pt, eta);
+    }
+
+    for (const auto& track : mcTracks) {
+      if (!isTrackSelected(track))
+        continue;
 
-          int id = identifyTrack(track, cent);
-          bool isPi = (id == KPidPionOne);
-          bool isKa = (id == KPidKaonTwo);
-          bool isPr = (id == KPidProtonThree);
-          bool isSpecies[KNsp] = {
-            true,
-            isPi && sign < 0, isPi && sign > 0, isPi,
-            isKa && sign < 0, isKa && sign > 0, isKa,
-            isPr && sign < 0, isPr && sign > 0, isPr};
+      float pt = track.pt(), eta = track.eta();
+      auto sign = track.sign();
+      fillNSigmaBefCut(track, cent);
+
+      int id = identifyTrack(track, cent);
+      bool isPi = (id == KPidPionOne);
+      bool isKa = (id == KPidKaonTwo);
+      bool isPr = (id == KPidProtonThree);
+      bool isSpecies[KNsp] = {
+        true,
+        isPi && sign < 0, isPi && sign > 0, isPi,
+        isKa && sign < 0, isKa && sign > 0, isKa,
+        isPr && sign < 0, isPr && sign > 0, isPr};
 
-          fillNSigmaAftCut(track, cent, isSpecies);
+      fillNSigmaAftCut(track, cent, isSpecies);
 
-          for (int isp = 0; isp < KNsp; ++isp) {
-            if (!isSpecies[isp])
-              continue;
+      for (int isp = 0; isp < KNsp; ++isp) {
+        if (!isSpecies[isp])
+          continue;
 
-            if (isp == kInclusiveIdx) {
-              histos.fill(HIST("h3_AllReco"), multPV, pt, eta);
-              if (track.has_mcParticle()) {
-                auto mcP = track.mcParticle();
-                if (mcP.isPhysicalPrimary()) {
-                  histos.fill(HIST("ptResolution"), mcP.pt(), (pt - mcP.pt()) / mcP.pt());
-                  histos.fill(HIST("etaResolution"), mcP.eta(), eta - mcP.eta());
-                  histos.fill(HIST("etaTruthReco"), mcP.eta(), eta);
-                  histos.fill(HIST("vzResolution"), mcP.vz(), (vz - mcP.vz()) / mcP.vz());
-                  histos.fill(HIST("TruthTracKVz"), mcP.vz(), vz);
-                  histos.fill(HIST("h3_RecoMatchedToPrimary"), multPV, mcP.pt(), mcP.eta());
-                } else {
-                  histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary"), multPV, pt, eta);
-                }
-              } else {
-                histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake"), multPV, pt, eta);
-              }
-            } else if (isp == kPiMinusIdx) {
-              histos.fill(HIST("h3_AllReco_PiMinus"), multPV, pt, eta);
-              if (track.has_mcParticle()) {
-                auto mcP = track.mcParticle();
-                if (mcP.isPhysicalPrimary()) {
-                  if (mcP.pdgCode() == -KPiPlus) {
-                    histos.fill(HIST("h3_RecoMatchedToPrimary_PiMinus"), multPV, mcP.pt(), mcP.eta());
-                  } else { // Misidentified
-                    histos.fill(HIST("h3_RecoMatchedToPrimary_MisID_PiMinus"), multPV, pt, eta);
-                  }
-                } else {
-                  histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary_PiMinus"), multPV, pt, eta);
-                }
-              } else { // No MC
-                histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake_PiMinus"), multPV, pt, eta);
+        if (isp == kInclusiveIdx) {
+          histos.fill(HIST("h3_AllReco"), multPV, pt, eta);
+          if (track.has_mcParticle()) {
+            auto mcP = track.mcParticle();
+            if (mcP.isPhysicalPrimary()) {
+              histos.fill(HIST("ptResolution"), mcP.pt(), (pt - mcP.pt()) / mcP.pt());
+              histos.fill(HIST("etaResolution"), mcP.eta(), eta - mcP.eta());
+              histos.fill(HIST("etaTruthReco"), mcP.eta(), eta);
+              histos.fill(HIST("vzResolution"), mcP.vz(), (vz - mcP.vz()) / mcP.vz());
+              histos.fill(HIST("TruthTracKVz"), mcP.vz(), vz);
+              histos.fill(HIST("h3_RecoMatchedToPrimary"), multPV, mcP.pt(), mcP.eta());
+            } else {
+              histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary"), multPV, pt, eta);
+            }
+          } else {
+            histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake"), multPV, pt, eta);
+          }
+        } else if (isp == kPiMinusIdx) {
+          histos.fill(HIST("h3_AllReco_PiMinus"), multPV, pt, eta);
+          if (track.has_mcParticle()) {
+            auto mcP = track.mcParticle();
+            if (mcP.isPhysicalPrimary()) {
+              if (mcP.pdgCode() == -KPiPlus) {
+                histos.fill(HIST("h3_RecoMatchedToPrimary_PiMinus"), multPV, mcP.pt(), mcP.eta());
+              } else { // Misidentified
+                histos.fill(HIST("h3_RecoMatchedToPrimary_MisID_PiMinus"), multPV, pt, eta);
               }
-            } else if (isp == kPiPlusIdx) {
-              histos.fill(HIST("h3_AllReco_PiPlus"), multPV, pt, eta);
-              if (track.has_mcParticle()) {
-                auto mcP = track.mcParticle();
-                if (mcP.isPhysicalPrimary()) {
-                  if (mcP.pdgCode() == KPiPlus) {
-                    histos.fill(HIST("h3_RecoMatchedToPrimary_PiPlus"), multPV, mcP.pt(), mcP.eta());
-                  } else { // Misidentified
-                    histos.fill(HIST("h3_RecoMatchedToPrimary_MisID_PiPlus"), multPV, pt, eta);
-                  }
-                } else {
-                  histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary_PiPlus"), multPV, pt, eta);
-                }
-              } else {
-                histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake_PiPlus"), multPV, pt, eta);
+            } else {
+              histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary_PiMinus"), multPV, pt, eta);
+            }
+          } else { // No MC
+            histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake_PiMinus"), multPV, pt, eta);
+          }
+        } else if (isp == kPiPlusIdx) {
+          histos.fill(HIST("h3_AllReco_PiPlus"), multPV, pt, eta);
+          if (track.has_mcParticle()) {
+            auto mcP = track.mcParticle();
+            if (mcP.isPhysicalPrimary()) {
+              if (mcP.pdgCode() == KPiPlus) {
+                histos.fill(HIST("h3_RecoMatchedToPrimary_PiPlus"), multPV, mcP.pt(), mcP.eta());
+              } else { // Misidentified
+                histos.fill(HIST("h3_RecoMatchedToPrimary_MisID_PiPlus"), multPV, pt, eta);
               }
-            } else if (isp == kPiAllIdx) {
-              histos.fill(HIST("h3_AllReco_PiAll"), multPV, pt, eta);
-              if (track.has_mcParticle()) {
-                auto mcP = track.mcParticle();
-                if (mcP.isPhysicalPrimary()) {
-                  if (std::abs(mcP.pdgCode()) == KPiPlus) {
-                    histos.fill(HIST("h3_RecoMatchedToPrimary_PiAll"), multPV, mcP.pt(), mcP.eta());
-                  } else { // Misidentified
-                    histos.fill(HIST("h3_RecoMatchedToPrimary_MisID_PiAll"), multPV, pt, eta);
-                  }
-                } else {
-                  histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary_PiAll"), multPV, pt, eta);
-                }
-              } else {
-                histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake_PiAll"), multPV, pt, eta);
+            } else {
+              histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary_PiPlus"), multPV, pt, eta);
+            }
+          } else {
+            histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake_PiPlus"), multPV, pt, eta);
+          }
+        } else if (isp == kPiAllIdx) {
+          histos.fill(HIST("h3_AllReco_PiAll"), multPV, pt, eta);
+          if (track.has_mcParticle()) {
+            auto mcP = track.mcParticle();
+            if (mcP.isPhysicalPrimary()) {
+              if (std::abs(mcP.pdgCode()) == KPiPlus) {
+                histos.fill(HIST("h3_RecoMatchedToPrimary_PiAll"), multPV, mcP.pt(), mcP.eta());
+              } else { // Misidentified
+                histos.fill(HIST("h3_RecoMatchedToPrimary_MisID_PiAll"), multPV, pt, eta);
               }
-            } else if (isp == kKaMinusIdx) {
-              histos.fill(HIST("h3_AllReco_KaMinus"), multPV, pt, eta);
-              if (track.has_mcParticle()) {
-                auto mcP = track.mcParticle();
-                if (mcP.isPhysicalPrimary()) {
-                  if (mcP.pdgCode() == -KKPlus) {
-                    histos.fill(HIST("h3_RecoMatchedToPrimary_KaMinus"), multPV, mcP.pt(), mcP.eta());
-                  } else { // Misidentified
-                    histos.fill(HIST("h3_RecoMatchedToPrimary_MisID_KaMinus"), multPV, pt, eta);
-                  }
-                } else {
-                  histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary_KaMinus"), multPV, pt, eta);
-                }
-              } else {
-                histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake_KaMinus"), multPV, pt, eta);
+            } else {
+              histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary_PiAll"), multPV, pt, eta);
+            }
+          } else {
+            histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake_PiAll"), multPV, pt, eta);
+          }
+        } else if (isp == kKaMinusIdx) {
+          histos.fill(HIST("h3_AllReco_KaMinus"), multPV, pt, eta);
+          if (track.has_mcParticle()) {
+            auto mcP = track.mcParticle();
+            if (mcP.isPhysicalPrimary()) {
+              if (mcP.pdgCode() == -KKPlus) {
+                histos.fill(HIST("h3_RecoMatchedToPrimary_KaMinus"), multPV, mcP.pt(), mcP.eta());
+              } else { // Misidentified
+                histos.fill(HIST("h3_RecoMatchedToPrimary_MisID_KaMinus"), multPV, pt, eta);
               }
-            } else if (isp == kKaPlusIdx) {
-              histos.fill(HIST("h3_AllReco_KaPlus"), multPV, pt, eta);
-              if (track.has_mcParticle()) {
-                auto mcP = track.mcParticle();
-                if (mcP.isPhysicalPrimary()) {
-                  if (mcP.pdgCode() == KKPlus) {
-                    histos.fill(HIST("h3_RecoMatchedToPrimary_KaPlus"), multPV, mcP.pt(), mcP.eta());
-                  } else { // Misidentified
-                    histos.fill(HIST("h3_RecoMatchedToPrimary_MisID_KaPlus"), multPV, pt, eta);
-                  }
-                } else {
-                  histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary_KaPlus"), multPV, pt, eta);
-                }
-              } else {
-                histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake_KaPlus"), multPV, pt, eta);
+            } else {
+              histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary_KaMinus"), multPV, pt, eta);
+            }
+          } else {
+            histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake_KaMinus"), multPV, pt, eta);
+          }
+        } else if (isp == kKaPlusIdx) {
+          histos.fill(HIST("h3_AllReco_KaPlus"), multPV, pt, eta);
+          if (track.has_mcParticle()) {
+            auto mcP = track.mcParticle();
+            if (mcP.isPhysicalPrimary()) {
+              if (mcP.pdgCode() == KKPlus) {
+                histos.fill(HIST("h3_RecoMatchedToPrimary_KaPlus"), multPV, mcP.pt(), mcP.eta());
+              } else { // Misidentified
+                histos.fill(HIST("h3_RecoMatchedToPrimary_MisID_KaPlus"), multPV, pt, eta);
               }
-            } else if (isp == kKaAllIdx) {
-              histos.fill(HIST("h3_AllReco_KaAll"), multPV, pt, eta);
-              if (track.has_mcParticle()) {
-                auto mcP = track.mcParticle();
-                if (mcP.isPhysicalPrimary()) {
-                  if (std::abs(mcP.pdgCode()) == KKPlus) {
-                    histos.fill(HIST("h3_RecoMatchedToPrimary_KaAll"), multPV, mcP.pt(), mcP.eta());
-                  } else { // Misidentified
-                    histos.fill(HIST("h3_RecoMatchedToPrimary_MisID_KaAll"), multPV, pt, eta);
-                  }
-                } else {
-                  histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary_KaAll"), multPV, pt, eta);
-                }
-              } else {
-                histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake_KaAll"), multPV, pt, eta);
+            } else {
+              histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary_KaPlus"), multPV, pt, eta);
+            }
+          } else {
+            histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake_KaPlus"), multPV, pt, eta);
+          }
+        } else if (isp == kKaAllIdx) {
+          histos.fill(HIST("h3_AllReco_KaAll"), multPV, pt, eta);
+          if (track.has_mcParticle()) {
+            auto mcP = track.mcParticle();
+            if (mcP.isPhysicalPrimary()) {
+              if (std::abs(mcP.pdgCode()) == KKPlus) {
+                histos.fill(HIST("h3_RecoMatchedToPrimary_KaAll"), multPV, mcP.pt(), mcP.eta());
+              } else { // Misidentified
+                histos.fill(HIST("h3_RecoMatchedToPrimary_MisID_KaAll"), multPV, pt, eta);
               }
-            } else if (isp == kAntiPrIdx) {
-              histos.fill(HIST("h3_AllReco_AntiPr"), multPV, pt, eta);
-              if (track.has_mcParticle()) {
-                auto mcP = track.mcParticle();
-                if (mcP.isPhysicalPrimary()) {
-                  if (mcP.pdgCode() == -KProton) {
-                    histos.fill(HIST("h3_RecoMatchedToPrimary_AntiPr"), multPV, mcP.pt(), mcP.eta());
-                  } else { // Misidentified
-                    histos.fill(HIST("h3_RecoMatchedToPrimary_MisID_AntiPr"), multPV, pt, eta);
-                  }
-                } else {
-                  histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary_AntiPr"), multPV, pt, eta);
-                }
-              } else {
-                histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake_AntiPr"), multPV, pt, eta);
+            } else {
+              histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary_KaAll"), multPV, pt, eta);
+            }
+          } else {
+            histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake_KaAll"), multPV, pt, eta);
+          }
+        } else if (isp == kAntiPrIdx) {
+          histos.fill(HIST("h3_AllReco_AntiPr"), multPV, pt, eta);
+          if (track.has_mcParticle()) {
+            auto mcP = track.mcParticle();
+            if (mcP.isPhysicalPrimary()) {
+              if (mcP.pdgCode() == -KProton) {
+                histos.fill(HIST("h3_RecoMatchedToPrimary_AntiPr"), multPV, mcP.pt(), mcP.eta());
+              } else { // Misidentified
+                histos.fill(HIST("h3_RecoMatchedToPrimary_MisID_AntiPr"), multPV, pt, eta);
               }
-            } else if (isp == kPrIdx) {
-              histos.fill(HIST("h3_AllReco_Pr"), multPV, pt, eta);
-              if (track.has_mcParticle()) {
-                auto mcP = track.mcParticle();
-                if (mcP.isPhysicalPrimary()) {
-                  if (mcP.pdgCode() == KProton) {
-                    histos.fill(HIST("h3_RecoMatchedToPrimary_Pr"), multPV, mcP.pt(), mcP.eta());
-                  } else { // Misidentified
-                    histos.fill(HIST("h3_RecoMatchedToPrimary_MisID_Pr"), multPV, pt, eta);
-                  }
-                } else {
-                  histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary_Pr"), multPV, pt, eta);
-                }
-              } else {
-                histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake_Pr"), multPV, pt, eta);
+            } else {
+              histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary_AntiPr"), multPV, pt, eta);
+            }
+          } else {
+            histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake_AntiPr"), multPV, pt, eta);
+          }
+        } else if (isp == kPrIdx) {
+          histos.fill(HIST("h3_AllReco_Pr"), multPV, pt, eta);
+          if (track.has_mcParticle()) {
+            auto mcP = track.mcParticle();
+            if (mcP.isPhysicalPrimary()) {
+              if (mcP.pdgCode() == KProton) {
+                histos.fill(HIST("h3_RecoMatchedToPrimary_Pr"), multPV, mcP.pt(), mcP.eta());
+              } else { // Misidentified
+                histos.fill(HIST("h3_RecoMatchedToPrimary_MisID_Pr"), multPV, pt, eta);
               }
-            } else if (isp == kPrAllIdx) {
-              histos.fill(HIST("h3_AllReco_PrAll"), multPV, pt, eta);
-              if (track.has_mcParticle()) {
-                auto mcP = track.mcParticle();
-                if (mcP.isPhysicalPrimary()) {
-                  if (std::abs(mcP.pdgCode()) == KProton) {
-                    histos.fill(HIST("h3_RecoMatchedToPrimary_PrAll"), multPV, mcP.pt(), mcP.eta());
-                  } else { // Misidentified
-                    histos.fill(HIST("h3_RecoMatchedToPrimary_MisID_PrAll"), multPV, pt, eta);
-                  }
-                } else {
-                  histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary_PrAll"), multPV, pt, eta);
-                }
-              } else {
-                histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake_PrAll"), multPV, pt, eta);
+            } else {
+              histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary_Pr"), multPV, pt, eta);
+            }
+          } else {
+            histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake_Pr"), multPV, pt, eta);
+          }
+        } else if (isp == kPrAllIdx) {
+          histos.fill(HIST("h3_AllReco_PrAll"), multPV, pt, eta);
+          if (track.has_mcParticle()) {
+            auto mcP = track.mcParticle();
+            if (mcP.isPhysicalPrimary()) {
+              if (std::abs(mcP.pdgCode()) == KProton) {
+                histos.fill(HIST("h3_RecoMatchedToPrimary_PrAll"), multPV, mcP.pt(), mcP.eta());
+              } else { // Misidentified
+                histos.fill(HIST("h3_RecoMatchedToPrimary_MisID_PrAll"), multPV, pt, eta);
               }
+            } else {
+              histos.fill(HIST("h3_RecoUnMatchedToPrimary_Secondary_PrAll"), multPV, pt, eta);
             }
+          } else {
+            histos.fill(HIST("h3_RecoUnMatchedToPrimary_Fake_PrAll"), multPV, pt, eta);
           }
         }
+      }
+    }
   }
   PROCESS_SWITCH(RadialFlowDecorr, processGetEffHists, "process MC to calculate EffWeights", cfgRunGetEff);
 
   void processMCFlat(MyRun3MCCollisions::iterator const& mcCollision, FilteredTCs const& mcTracks)
   {
-        histos.fill(HIST("hVtxZ"), mcCollision.posZ());
-        if (!mcCollision.has_mcCollision() || !isEventSelected(mcCollision))
-          return;
-
-        float cent = getCentrality(mcCollision);
-        if (cent > KCentMax)
-          return;
-
-        float multPV = mcCollision.multNTracksPV();
-        float vz = mcCollision.posZ();
-
-        if (!isPassAddPileup(multPV, mcTracks.size(), cent))
-          return;
-        histos.fill(HIST("hVtxZ_after_sel"), mcCollision.posZ());
-        histos.fill(HIST("hCentrality"), cent);
-
-        histos.fill(HIST("Hist2D_globalTracks_PVTracks"), mcCollision.multNTracksPV(), mcTracks.size());
-        histos.fill(HIST("Hist2D_cent_nch"), mcTracks.size(), cent);
-        histos.fill(HIST("Hist2D_globalTracks_cent"), cent, mcTracks.size());
-        histos.fill(HIST("Hist2D_PVTracks_cent"), cent, multPV);
-        for (const auto& track : mcTracks) {
-          if (!isTrackSelected(track))
-            continue;
+    histos.fill(HIST("hVtxZ"), mcCollision.posZ());
+    if (!mcCollision.has_mcCollision() || !isEventSelected(mcCollision))
+      return;
 
-          float pt = track.pt(), eta = track.eta(), phi = track.phi();
-          auto sign = track.sign();
-          int id = identifyTrack(track, cent);
-          bool isPi = (id == KPidPionOne);
-          bool isKa = (id == KPidKaonTwo);
-          bool isPr = (id == KPidProtonThree);
-          bool isSpecies[KNsp] = {
-            true,
-            isPi && sign < 0, isPi && sign > 0, isPi,
-            isKa && sign < 0, isKa && sign > 0, isKa,
-            isPr && sign < 0, isPr && sign > 0, isPr};
+    float cent = getCentrality(mcCollision);
+    if (cent > KCentMax)
+      return;
 
-          for (int isp = 0; isp < KNsp; ++isp) {
-            if (!isSpecies[isp])
-              continue;
+    float multPV = mcCollision.multNTracksPV();
+    float vz = mcCollision.posZ();
 
-            float eff = getEfficiency(multPV, pt, eta, static_cast<PIDIdx>(isp), 0, cfgEff);
-            float fake = getEfficiency(multPV, pt, eta, static_cast<PIDIdx>(isp), 1, cfgEff);
-            float w = (eff > KFloatEpsilon) ? (1.0f - fake) / eff : 0.0f;
+    if (!isPassAddPileup(multPV, mcTracks.size(), cent))
+      return;
+    histos.fill(HIST("hVtxZ_after_sel"), mcCollision.posZ());
+    histos.fill(HIST("hCentrality"), cent);
 
-            if (std::isfinite(w) && w > 0.f) {
-              if (isp == kInclusiveIdx) {
-                histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd"), vz, sign, pt, eta, phi, w);
-                histos.fill(HIST("MCReco/hEtaPhiReco"), vz, sign, pt, eta, phi, 1.0);
-                histos.fill(HIST("MCReco/hEtaPhiRecoWtd"), vz, sign, pt, eta, phi, w);
-              } else if (isp == kPiMinusIdx) {
-                histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd_PiMinus"), vz, sign, pt, eta, phi, w);
-                histos.fill(HIST("MCReco/hEtaPhiReco_PiMinus"), vz, sign, pt, eta, phi, 1.0);
-                histos.fill(HIST("MCReco/hEtaPhiRecoWtd_PiMinus"), vz, sign, pt, eta, phi, w);
-              } else if (isp == kPiPlusIdx) {
-                histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd_PiPlus"), vz, sign, pt, eta, phi, w);
-                histos.fill(HIST("MCReco/hEtaPhiReco_PiPlus"), vz, sign, pt, eta, phi, 1.0);
-                histos.fill(HIST("MCReco/hEtaPhiRecoWtd_PiPlus"), vz, sign, pt, eta, phi, w);
-              } else if (isp == kPiAllIdx) {
-                histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd_PiAll"), vz, sign, pt, eta, phi, w);
-                histos.fill(HIST("MCReco/hEtaPhiReco_PiAll"), vz, sign, pt, eta, phi, 1.0);
-                histos.fill(HIST("MCReco/hEtaPhiRecoWtd_PiAll"), vz, sign, pt, eta, phi, w);
-              } else if (isp == kKaMinusIdx) {
-                histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd_KaMinus"), vz, sign, pt, eta, phi, w);
-                histos.fill(HIST("MCReco/hEtaPhiReco_KaMinus"), vz, sign, pt, eta, phi, 1.0);
-                histos.fill(HIST("MCReco/hEtaPhiRecoWtd_KaMinus"), vz, sign, pt, eta, phi, w);
-              } else if (isp == kKaPlusIdx) {
-                histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd_KaPlus"), vz, sign, pt, eta, phi, w);
-                histos.fill(HIST("MCReco/hEtaPhiReco_KaPlus"), vz, sign, pt, eta, phi, 1.0);
-                histos.fill(HIST("MCReco/hEtaPhiRecoWtd_KaPlus"), vz, sign, pt, eta, phi, w);
-              } else if (isp == kKaAllIdx) {
-                histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd_KaAll"), vz, sign, pt, eta, phi, w);
-                histos.fill(HIST("MCReco/hEtaPhiReco_KaAll"), vz, sign, pt, eta, phi, 1.0);
-                histos.fill(HIST("MCReco/hEtaPhiRecoWtd_KaAll"), vz, sign, pt, eta, phi, w);
-              } else if (isp == kAntiPrIdx) {
-                histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd_AntiPr"), vz, sign, pt, eta, phi, w);
-                histos.fill(HIST("MCReco/hEtaPhiReco_AntiPr"), vz, sign, pt, eta, phi, 1.0);
-                histos.fill(HIST("MCReco/hEtaPhiRecoWtd_AntiPr"), vz, sign, pt, eta, phi, w);
-              } else if (isp == kPrIdx) {
-                histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd_Pr"), vz, sign, pt, eta, phi, w);
-                histos.fill(HIST("MCReco/hEtaPhiReco_Pr"), vz, sign, pt, eta, phi, 1.0);
-                histos.fill(HIST("MCReco/hEtaPhiRecoWtd_Pr"), vz, sign, pt, eta, phi, w);
-              } else if (isp == kPrAllIdx) {
-                histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd_PrAll"), vz, sign, pt, eta, phi, w);
-                histos.fill(HIST("MCReco/hEtaPhiReco_PrAll"), vz, sign, pt, eta, phi, 1.0);
-                histos.fill(HIST("MCReco/hEtaPhiRecoWtd_PrAll"), vz, sign, pt, eta, phi, w);
-              }
-            }
+    histos.fill(HIST("Hist2D_globalTracks_PVTracks"), mcCollision.multNTracksPV(), mcTracks.size());
+    histos.fill(HIST("Hist2D_cent_nch"), mcTracks.size(), cent);
+    histos.fill(HIST("Hist2D_globalTracks_cent"), cent, mcTracks.size());
+    histos.fill(HIST("Hist2D_PVTracks_cent"), cent, multPV);
+    for (const auto& track : mcTracks) {
+      if (!isTrackSelected(track))
+        continue;
+
+      float pt = track.pt(), eta = track.eta(), phi = track.phi();
+      auto sign = track.sign();
+      int id = identifyTrack(track, cent);
+      bool isPi = (id == KPidPionOne);
+      bool isKa = (id == KPidKaonTwo);
+      bool isPr = (id == KPidProtonThree);
+      bool isSpecies[KNsp] = {
+        true,
+        isPi && sign < 0, isPi && sign > 0, isPi,
+        isKa && sign < 0, isKa && sign > 0, isKa,
+        isPr && sign < 0, isPr && sign > 0, isPr};
+
+      for (int isp = 0; isp < KNsp; ++isp) {
+        if (!isSpecies[isp])
+          continue;
+
+        float eff = getEfficiency(multPV, pt, eta, static_cast<PIDIdx>(isp), 0, cfgEff);
+        float fake = getEfficiency(multPV, pt, eta, static_cast<PIDIdx>(isp), 1, cfgEff);
+        float w = (eff > KFloatEpsilon) ? (1.0f - fake) / eff : 0.0f;
+
+        if (std::isfinite(w) && w > 0.f) {
+          if (isp == kInclusiveIdx) {
+            histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd"), vz, sign, pt, eta, phi, w);
+            histos.fill(HIST("MCReco/hEtaPhiReco"), vz, sign, pt, eta, phi, 1.0);
+            histos.fill(HIST("MCReco/hEtaPhiRecoWtd"), vz, sign, pt, eta, phi, w);
+          } else if (isp == kPiMinusIdx) {
+            histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd_PiMinus"), vz, sign, pt, eta, phi, w);
+            histos.fill(HIST("MCReco/hEtaPhiReco_PiMinus"), vz, sign, pt, eta, phi, 1.0);
+            histos.fill(HIST("MCReco/hEtaPhiRecoWtd_PiMinus"), vz, sign, pt, eta, phi, w);
+          } else if (isp == kPiPlusIdx) {
+            histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd_PiPlus"), vz, sign, pt, eta, phi, w);
+            histos.fill(HIST("MCReco/hEtaPhiReco_PiPlus"), vz, sign, pt, eta, phi, 1.0);
+            histos.fill(HIST("MCReco/hEtaPhiRecoWtd_PiPlus"), vz, sign, pt, eta, phi, w);
+          } else if (isp == kPiAllIdx) {
+            histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd_PiAll"), vz, sign, pt, eta, phi, w);
+            histos.fill(HIST("MCReco/hEtaPhiReco_PiAll"), vz, sign, pt, eta, phi, 1.0);
+            histos.fill(HIST("MCReco/hEtaPhiRecoWtd_PiAll"), vz, sign, pt, eta, phi, w);
+          } else if (isp == kKaMinusIdx) {
+            histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd_KaMinus"), vz, sign, pt, eta, phi, w);
+            histos.fill(HIST("MCReco/hEtaPhiReco_KaMinus"), vz, sign, pt, eta, phi, 1.0);
+            histos.fill(HIST("MCReco/hEtaPhiRecoWtd_KaMinus"), vz, sign, pt, eta, phi, w);
+          } else if (isp == kKaPlusIdx) {
+            histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd_KaPlus"), vz, sign, pt, eta, phi, w);
+            histos.fill(HIST("MCReco/hEtaPhiReco_KaPlus"), vz, sign, pt, eta, phi, 1.0);
+            histos.fill(HIST("MCReco/hEtaPhiRecoWtd_KaPlus"), vz, sign, pt, eta, phi, w);
+          } else if (isp == kKaAllIdx) {
+            histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd_KaAll"), vz, sign, pt, eta, phi, w);
+            histos.fill(HIST("MCReco/hEtaPhiReco_KaAll"), vz, sign, pt, eta, phi, 1.0);
+            histos.fill(HIST("MCReco/hEtaPhiRecoWtd_KaAll"), vz, sign, pt, eta, phi, w);
+          } else if (isp == kAntiPrIdx) {
+            histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd_AntiPr"), vz, sign, pt, eta, phi, w);
+            histos.fill(HIST("MCReco/hEtaPhiReco_AntiPr"), vz, sign, pt, eta, phi, 1.0);
+            histos.fill(HIST("MCReco/hEtaPhiRecoWtd_AntiPr"), vz, sign, pt, eta, phi, w);
+          } else if (isp == kPrIdx) {
+            histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd_Pr"), vz, sign, pt, eta, phi, w);
+            histos.fill(HIST("MCReco/hEtaPhiReco_Pr"), vz, sign, pt, eta, phi, 1.0);
+            histos.fill(HIST("MCReco/hEtaPhiRecoWtd_Pr"), vz, sign, pt, eta, phi, w);
+          } else if (isp == kPrAllIdx) {
+            histos.fill(HIST("MCReco/hEtaPhiRecoEffWtd_PrAll"), vz, sign, pt, eta, phi, w);
+            histos.fill(HIST("MCReco/hEtaPhiReco_PrAll"), vz, sign, pt, eta, phi, 1.0);
+            histos.fill(HIST("MCReco/hEtaPhiRecoWtd_PrAll"), vz, sign, pt, eta, phi, w);
           }
         }
+      }
+    }
   }
   PROCESS_SWITCH(RadialFlowDecorr, processMCFlat, "process MC to calculate FlatWeights", cfgRunGetMCFlat);
 
@@ -1813,66 +1813,66 @@ struct RadialFlowDecorr {
     double sumWiTruth[KNsp][KNEta]{}, sumWiptiTruth[KNsp][KNEta]{};
     double sumWiReco[KNsp][KNEta]{}, sumWiptiReco[KNsp][KNEta]{};
     double sumWiRecoEffCorr[KNsp][KNEta]{}, sumWiptiRecoEffCorr[KNsp][KNEta]{};
-        histos.fill(HIST("hVtxZ"), mcCollision.posZ());
-        if (!mcCollision.has_mcCollision() || !isEventSelected(mcCollision))
-          return;
-        float cent = getCentrality(mcCollision);
-        if (cent > KCentMax)
-          return;
-        float multPV = mcCollision.multNTracksPV();
-        float vz = mcCollision.posZ();
-        if (!isPassAddPileup(multPV, mcTracks.size(), cent))
-          return;
-
-        histos.fill(HIST("hVtxZ_after_sel"), mcCollision.posZ());
-        histos.fill(HIST("hCentrality"), cent);
-
-        histos.fill(HIST("Hist2D_globalTracks_PVTracks"), mcCollision.multNTracksPV(), mcTracks.size());
-        histos.fill(HIST("Hist2D_cent_nch"), mcTracks.size(), cent);
-        histos.fill(HIST("Hist2D_globalTracks_cent"), cent, mcTracks.size());
-        histos.fill(HIST("Hist2D_PVTracks_cent"), cent, multPV);
-
-        memset(sumWiTruth, 0, sizeof(sumWiTruth));
-        memset(sumWiptiTruth, 0, sizeof(sumWiptiTruth));
-        memset(sumWiReco, 0, sizeof(sumWiReco));
-        memset(sumWiptiReco, 0, sizeof(sumWiptiReco));
-        memset(sumWiRecoEffCorr, 0, sizeof(sumWiRecoEffCorr));
-        memset(sumWiptiRecoEffCorr, 0, sizeof(sumWiptiRecoEffCorr));
-
-        for (const auto& particle : mcParticles) {
-          if (!isParticleSelected(particle) || !particle.isPhysicalPrimary())
-            continue;
-          float pt = particle.pt(), eta = particle.eta();
-          if (pt <= cfgPtMin || pt > cfgPtMax)
-            continue;
-          int pdgCode = particle.pdgCode();
-          int absPdg = std::abs(pdgCode);
+    histos.fill(HIST("hVtxZ"), mcCollision.posZ());
+    if (!mcCollision.has_mcCollision() || !isEventSelected(mcCollision))
+      return;
+    float cent = getCentrality(mcCollision);
+    if (cent > KCentMax)
+      return;
+    float multPV = mcCollision.multNTracksPV();
+    float vz = mcCollision.posZ();
+    if (!isPassAddPileup(multPV, mcTracks.size(), cent))
+      return;
 
-          bool isSpecies[KNsp] = {
-            true,                // kInclusiveIdx
-            pdgCode == -KPiPlus, // kPiMinusIdx
-            pdgCode == KPiPlus,  // kPiPlusIdx
-            absPdg == KPiPlus,   // kPiAllIdx
-            pdgCode == -KKPlus,  // kKaMinusIdx
-            pdgCode == KKPlus,   // kKaPlusIdx
-            absPdg == KKPlus,    // kKaAllIdx
-            pdgCode == -KProton, // kAntiPrIdx
-            pdgCode == KProton,  // kPrIdx
-            absPdg == KProton    // kPrAllIdx
-          };
+    histos.fill(HIST("hVtxZ_after_sel"), mcCollision.posZ());
+    histos.fill(HIST("hCentrality"), cent);
 
-          for (int ieta = 0; ieta < KNEta; ++ieta) {
-            if (eta <= etaLw[ieta] || eta > etaUp[ieta])
-              continue;
+    histos.fill(HIST("Hist2D_globalTracks_PVTracks"), mcCollision.multNTracksPV(), mcTracks.size());
+    histos.fill(HIST("Hist2D_cent_nch"), mcTracks.size(), cent);
+    histos.fill(HIST("Hist2D_globalTracks_cent"), cent, mcTracks.size());
+    histos.fill(HIST("Hist2D_PVTracks_cent"), cent, multPV);
 
-            for (int isp = 0; isp < KNsp; ++isp) {
-              if (isSpecies[isp]) {
-                sumWiTruth[isp][ieta]++;
-                sumWiptiTruth[isp][ieta] += pt;
-              }
-            }
+    memset(sumWiTruth, 0, sizeof(sumWiTruth));
+    memset(sumWiptiTruth, 0, sizeof(sumWiptiTruth));
+    memset(sumWiReco, 0, sizeof(sumWiReco));
+    memset(sumWiptiReco, 0, sizeof(sumWiptiReco));
+    memset(sumWiRecoEffCorr, 0, sizeof(sumWiRecoEffCorr));
+    memset(sumWiptiRecoEffCorr, 0, sizeof(sumWiptiRecoEffCorr));
+
+    for (const auto& particle : mcParticles) {
+      if (!isParticleSelected(particle) || !particle.isPhysicalPrimary())
+        continue;
+      float pt = particle.pt(), eta = particle.eta();
+      if (pt <= cfgPtMin || pt > cfgPtMax)
+        continue;
+      int pdgCode = particle.pdgCode();
+      int absPdg = std::abs(pdgCode);
+
+      bool isSpecies[KNsp] = {
+        true,                // kInclusiveIdx
+        pdgCode == -KPiPlus, // kPiMinusIdx
+        pdgCode == KPiPlus,  // kPiPlusIdx
+        absPdg == KPiPlus,   // kPiAllIdx
+        pdgCode == -KKPlus,  // kKaMinusIdx
+        pdgCode == KKPlus,   // kKaPlusIdx
+        absPdg == KKPlus,    // kKaAllIdx
+        pdgCode == -KProton, // kAntiPrIdx
+        pdgCode == KProton,  // kPrIdx
+        absPdg == KProton    // kPrAllIdx
+      };
+
+      for (int ieta = 0; ieta < KNEta; ++ieta) {
+        if (eta <= etaLw[ieta] || eta > etaUp[ieta])
+          continue;
+
+        for (int isp = 0; isp < KNsp; ++isp) {
+          if (isSpecies[isp]) {
+            sumWiTruth[isp][ieta]++;
+            sumWiptiTruth[isp][ieta] += pt;
           }
         }
+      }
+    }
 
         for (int isp = 0; isp < KNsp; ++isp) {
           histos.fill(HIST("MCGen/Prof_Cent_Nsp_Nchrec"), cent, isp, sumWiTruth[isp][0]);
@@ -2152,182 +2152,179 @@ struct RadialFlowDecorr {
     double p1kBarTru[KNsp][KNEta]{}, p1kBarReco[KNsp][KNEta]{}, p1kBarRecoEffCor[KNsp][KNEta]{};
     double p1kBarTruMult[KNsp][KNEta]{}, p1kBarRecoMult[KNsp][KNEta]{}, p1kBarRecoEffCorMult[KNsp][KNEta]{};
 
-        if (!mcCollision.has_mcCollision() || !isEventSelected(mcCollision))
-          return;
-        float cent = getCentrality(mcCollision);
-        if (cent > KCentMax)
-          return;
-        float multPV = mcCollision.multNTracksPV();
-        float vz = mcCollision.posZ();
-        if (!isPassAddPileup(multPV, mcTracks.size(), cent))
-          return;
-
-        histos.fill(HIST("hVtxZ_after_sel"), mcCollision.posZ());
-        histos.fill(HIST("hCentrality"), cent);
-
-        histos.fill(HIST("Hist2D_globalTracks_PVTracks"), multPV, mcTracks.size());
-        histos.fill(HIST("Hist2D_cent_nch"), mcTracks.size(), cent);
-        histos.fill(HIST("Hist2D_globalTracks_cent"), cent, mcTracks.size());
-        histos.fill(HIST("Hist2D_PVTracks_cent"), cent, multPV);
-
-        memset(sumPmwkTru, 0, sizeof(sumPmwkTru));
-        memset(sumWkTru, 0, sizeof(sumWkTru));
-        memset(sumPmwkReco, 0, sizeof(sumPmwkReco));
-        memset(sumWkReco, 0, sizeof(sumWkReco));
-        memset(sumPmwkRecoEffCor, 0, sizeof(sumPmwkRecoEffCor));
-        memset(sumWkRecoEffCor, 0, sizeof(sumWkRecoEffCor));
-
-        memset(meanTru, 0, sizeof(meanTru));
-        memset(c2Tru, 0, sizeof(c2Tru));
-        memset(meanReco, 0, sizeof(meanReco));
-        memset(c2Reco, 0, sizeof(c2Reco));
-        memset(meanRecoEffCor, 0, sizeof(meanRecoEffCor));
-        memset(c2RecoEffCor, 0, sizeof(c2RecoEffCor));
-
-        memset(meanTruMult, 0, sizeof(meanTruMult));
-        memset(meanRecoMult, 0, sizeof(meanRecoMult));
-        memset(meanRecoEffCorMult, 0, sizeof(meanRecoEffCorMult));
-
-        memset(p1kBarTru, 0, sizeof(p1kBarTru));
-        memset(p1kBarReco, 0, sizeof(p1kBarReco));
-        memset(p1kBarRecoEffCor, 0, sizeof(p1kBarRecoEffCor));
-
-        memset(p1kBarTruMult, 0, sizeof(p1kBarTruMult));
-        memset(p1kBarRecoMult, 0, sizeof(p1kBarRecoMult));
-        memset(p1kBarRecoEffCorMult, 0, sizeof(p1kBarRecoEffCorMult));
-
-        double p1kBarFt0A = 0.0, p1kBarFt0C = 0.0;
-
+    if (!mcCollision.has_mcCollision() || !isEventSelected(mcCollision))
+      return;
+    float cent = getCentrality(mcCollision);
+    if (cent > KCentMax)
+      return;
+    float multPV = mcCollision.multNTracksPV();
+    float vz = mcCollision.posZ();
+    if (!isPassAddPileup(multPV, mcTracks.size(), cent))
+      return;
 
+    histos.fill(HIST("hVtxZ_after_sel"), mcCollision.posZ());
+    histos.fill(HIST("hCentrality"), cent);
 
-        for (const auto& particle : mcParticles) {
-          if (!isParticleSelected(particle) || !particle.isPhysicalPrimary())
-            continue;
+    histos.fill(HIST("Hist2D_globalTracks_PVTracks"), multPV, mcTracks.size());
+    histos.fill(HIST("Hist2D_cent_nch"), mcTracks.size(), cent);
+    histos.fill(HIST("Hist2D_globalTracks_cent"), cent, mcTracks.size());
+    histos.fill(HIST("Hist2D_PVTracks_cent"), cent, multPV);
+
+    memset(sumPmwkTru, 0, sizeof(sumPmwkTru));
+    memset(sumWkTru, 0, sizeof(sumWkTru));
+    memset(sumPmwkReco, 0, sizeof(sumPmwkReco));
+    memset(sumWkReco, 0, sizeof(sumWkReco));
+    memset(sumPmwkRecoEffCor, 0, sizeof(sumPmwkRecoEffCor));
+    memset(sumWkRecoEffCor, 0, sizeof(sumWkRecoEffCor));
+
+    memset(meanTru, 0, sizeof(meanTru));
+    memset(c2Tru, 0, sizeof(c2Tru));
+    memset(meanReco, 0, sizeof(meanReco));
+    memset(c2Reco, 0, sizeof(c2Reco));
+    memset(meanRecoEffCor, 0, sizeof(meanRecoEffCor));
+    memset(c2RecoEffCor, 0, sizeof(c2RecoEffCor));
+
+    memset(meanTruMult, 0, sizeof(meanTruMult));
+    memset(meanRecoMult, 0, sizeof(meanRecoMult));
+    memset(meanRecoEffCorMult, 0, sizeof(meanRecoEffCorMult));
+
+    memset(p1kBarTru, 0, sizeof(p1kBarTru));
+    memset(p1kBarReco, 0, sizeof(p1kBarReco));
+    memset(p1kBarRecoEffCor, 0, sizeof(p1kBarRecoEffCor));
+
+    memset(p1kBarTruMult, 0, sizeof(p1kBarTruMult));
+    memset(p1kBarRecoMult, 0, sizeof(p1kBarRecoMult));
+    memset(p1kBarRecoEffCorMult, 0, sizeof(p1kBarRecoEffCorMult));
+
+    double p1kBarFt0A = 0.0, p1kBarFt0C = 0.0;
+
+    for (const auto& particle : mcParticles) {
+      if (!isParticleSelected(particle) || !particle.isPhysicalPrimary())
+        continue;
 
-          float pt = particle.pt();
-          if (pt <= cfgPtMin || pt > cfgPtMax)
-            continue;
-          float eta = particle.eta();
-          int pdgCode = particle.pdgCode();
-          int absPdg = std::abs(pdgCode);
+      float pt = particle.pt();
+      if (pt <= cfgPtMin || pt > cfgPtMax)
+        continue;
+      float eta = particle.eta();
+      int pdgCode = particle.pdgCode();
+      int absPdg = std::abs(pdgCode);
 
-          bool isSpecies[KNsp] = {
-            true,                // kInclusiveIdx
-            pdgCode == -KPiPlus, // kPiMinusIdx
-            pdgCode == KPiPlus,  // kPiPlusIdx
-            absPdg == KPiPlus,   // kPiAllIdx
-            pdgCode == -KKPlus,  // kKaMinusIdx
-            pdgCode == KKPlus,   // kKaPlusIdx
-            absPdg == KKPlus,    // kKaAllIdx
-            pdgCode == -KProton, // kAntiPrIdx
-            pdgCode == KProton,  // kPrIdx
-            absPdg == KProton    // kPrAllIdx
-          };
+      bool isSpecies[KNsp] = {
+        true,                // kInclusiveIdx
+        pdgCode == -KPiPlus, // kPiMinusIdx
+        pdgCode == KPiPlus,  // kPiPlusIdx
+        absPdg == KPiPlus,   // kPiAllIdx
+        pdgCode == -KKPlus,  // kKaMinusIdx
+        pdgCode == KKPlus,   // kKaPlusIdx
+        absPdg == KKPlus,    // kKaAllIdx
+        pdgCode == -KProton, // kAntiPrIdx
+        pdgCode == KProton,  // kPrIdx
+        absPdg == KProton    // kPrAllIdx
+      };
 
-          for (int ieta = 0; ieta < KNEta; ++ieta) {
-            if (eta <= etaLw[ieta] || eta > etaUp[ieta])
-              continue;
-            for (int isp = 0; isp < KNsp; ++isp) {
-              if (isSpecies[isp]) {
-                for (int k = 0; k < KIntK; ++k) {
-                  for (int m = 0; m < KIntM; ++m) {
-                    sumPmwkTru[isp][ieta][m][k] += std::pow(pt, m);
-                  }
-                  sumWkTru[isp][ieta][k]++;
-                }
+      for (int ieta = 0; ieta < KNEta; ++ieta) {
+        if (eta <= etaLw[ieta] || eta > etaUp[ieta])
+          continue;
+        for (int isp = 0; isp < KNsp; ++isp) {
+          if (isSpecies[isp]) {
+            for (int k = 0; k < KIntK; ++k) {
+              for (int m = 0; m < KIntM; ++m) {
+                sumPmwkTru[isp][ieta][m][k] += std::pow(pt, m);
               }
+              sumWkTru[isp][ieta][k]++;
             }
           }
-        } // end truth loop
-
-
-        for (const auto& track : mcTracks) {
-          if (!isTrackSelected(track))
-            continue;
+        }
+      }
+    } // end truth loop
 
-          float pt = track.pt();
-          if (pt <= cfgPtMin || pt > cfgPtMax)
-            continue;
-          float eta = track.eta();
-          float phi = track.phi();
-          auto sign = track.sign();
-          int id = identifyTrack(track, cent);
-          bool isPi = (id == KPidPionOne);
-          bool isKa = (id == KPidKaonTwo);
-          bool isPr = (id == KPidProtonThree);
-          bool isSpecies[KNsp] = {
-            true,
-            isPi && sign < 0, isPi && sign > 0, isPi,
-            isKa && sign < 0, isKa && sign > 0, isKa,
-            isPr && sign < 0, isPr && sign > 0, isPr};
+    for (const auto& track : mcTracks) {
+      if (!isTrackSelected(track))
+        continue;
 
-          for (int isp = 0; isp < KNsp; ++isp) {
-            if (!isSpecies[isp])
-              continue;
-            float eff = getEfficiency(multPV, pt, eta, static_cast<PIDIdx>(isp), 0, cfgEff);
-            float fake = getEfficiency(multPV, pt, eta, static_cast<PIDIdx>(isp), 1, cfgEff);
-            float flatW = getFlatteningWeight(vz, sign, pt, eta, phi, static_cast<PIDIdx>(isp), cfgFlat);
-            float w = flatW * (1.0 - fake) / eff;
+      float pt = track.pt();
+      if (pt <= cfgPtMin || pt > cfgPtMax)
+        continue;
+      float eta = track.eta();
+      float phi = track.phi();
+      auto sign = track.sign();
+      int id = identifyTrack(track, cent);
+      bool isPi = (id == KPidPionOne);
+      bool isKa = (id == KPidKaonTwo);
+      bool isPr = (id == KPidProtonThree);
+      bool isSpecies[KNsp] = {
+        true,
+        isPi && sign < 0, isPi && sign > 0, isPi,
+        isKa && sign < 0, isKa && sign > 0, isKa,
+        isPr && sign < 0, isPr && sign > 0, isPr};
 
-            if (!std::isfinite(w) || w <= 0.f || eff <= KFloatEpsilon)
-              continue;
+      for (int isp = 0; isp < KNsp; ++isp) {
+        if (!isSpecies[isp])
+          continue;
+        float eff = getEfficiency(multPV, pt, eta, static_cast<PIDIdx>(isp), 0, cfgEff);
+        float fake = getEfficiency(multPV, pt, eta, static_cast<PIDIdx>(isp), 1, cfgEff);
+        float flatW = getFlatteningWeight(vz, sign, pt, eta, phi, static_cast<PIDIdx>(isp), cfgFlat);
+        float w = flatW * (1.0 - fake) / eff;
 
-            for (int ieta = 0; ieta < KNEta; ++ieta) {
-              if (eta <= etaLw[ieta] || eta > etaUp[ieta])
-                continue;
-              for (int k = 0; k < KIntK; ++k) {
-                for (int m = 0; m < KIntM; ++m) {
-                  sumPmwkReco[isp][ieta][m][k] += std::pow(1.0, k) * std::pow(pt, m);
-                  sumPmwkRecoEffCor[isp][ieta][m][k] += std::pow(w, k) * std::pow(pt, m);
-                }
-                sumWkReco[isp][ieta][k] += std::pow(1.0, k);
-                sumWkRecoEffCor[isp][ieta][k] += std::pow(w, k);
-              }
-            }
+        if (!std::isfinite(w) || w <= 0.f || eff <= KFloatEpsilon)
+          continue;
 
-            if (isp == kInclusiveIdx) {
-              histos.fill(HIST("hEtaPhiReco"), vz, sign, pt, eta, phi);
-              histos.fill(HIST("hEtaPhiRecoWtd"), vz, sign, pt, eta, phi, w);
-              histos.fill(HIST("hEtaPhiRecoEffWtd"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
-            } else if (isp == kPiMinusIdx) {
-              histos.fill(HIST("hEtaPhiReco_PiMinus"), vz, sign, pt, eta, phi);
-              histos.fill(HIST("hEtaPhiRecoWtd_PiMinus"), vz, sign, pt, eta, phi, w);
-              histos.fill(HIST("hEtaPhiRecoEffWtd_PiMinus"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
-            } else if (isp == kPiPlusIdx) {
-              histos.fill(HIST("hEtaPhiReco_PiPlus"), vz, sign, pt, eta, phi);
-              histos.fill(HIST("hEtaPhiRecoWtd_PiPlus"), vz, sign, pt, eta, phi, w);
-              histos.fill(HIST("hEtaPhiRecoEffWtd_PiPlus"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
-            } else if (isp == kPiAllIdx) {
-              histos.fill(HIST("hEtaPhiReco_PiAll"), vz, sign, pt, eta, phi);
-              histos.fill(HIST("hEtaPhiRecoWtd_PiAll"), vz, sign, pt, eta, phi, w);
-              histos.fill(HIST("hEtaPhiRecoEffWtd_PiAll"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
-            } else if (isp == kKaMinusIdx) {
-              histos.fill(HIST("hEtaPhiReco_KaMinus"), vz, sign, pt, eta, phi);
-              histos.fill(HIST("hEtaPhiRecoWtd_KaMinus"), vz, sign, pt, eta, phi, w);
-              histos.fill(HIST("hEtaPhiRecoEffWtd_KaMinus"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
-            } else if (isp == kKaPlusIdx) {
-              histos.fill(HIST("hEtaPhiReco_KaPlus"), vz, sign, pt, eta, phi);
-              histos.fill(HIST("hEtaPhiRecoWtd_KaPlus"), vz, sign, pt, eta, phi, w);
-              histos.fill(HIST("hEtaPhiRecoEffWtd_KaPlus"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
-            } else if (isp == kKaAllIdx) {
-              histos.fill(HIST("hEtaPhiReco_KaAll"), vz, sign, pt, eta, phi);
-              histos.fill(HIST("hEtaPhiRecoWtd_KaAll"), vz, sign, pt, eta, phi, w);
-              histos.fill(HIST("hEtaPhiRecoEffWtd_KaAll"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
-            } else if (isp == kPrIdx) {
-              histos.fill(HIST("hEtaPhiReco_Pr"), vz, sign, pt, eta, phi);
-              histos.fill(HIST("hEtaPhiRecoWtd_Pr"), vz, sign, pt, eta, phi, w);
-              histos.fill(HIST("hEtaPhiRecoEffWtd_Pr"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
-            } else if (isp == kAntiPrIdx) {
-              histos.fill(HIST("hEtaPhiReco_AntiPr"), vz, sign, pt, eta, phi);
-              histos.fill(HIST("hEtaPhiRecoWtd_AntiPr"), vz, sign, pt, eta, phi, w);
-              histos.fill(HIST("hEtaPhiRecoEffWtd_AntiPr"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
-            } else if (isp == kPrAllIdx) {
-              histos.fill(HIST("hEtaPhiReco_PrAll"), vz, sign, pt, eta, phi);
-              histos.fill(HIST("hEtaPhiRecoWtd_PrAll"), vz, sign, pt, eta, phi, w);
-              histos.fill(HIST("hEtaPhiRecoEffWtd_PrAll"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
+        for (int ieta = 0; ieta < KNEta; ++ieta) {
+          if (eta <= etaLw[ieta] || eta > etaUp[ieta])
+            continue;
+          for (int k = 0; k < KIntK; ++k) {
+            for (int m = 0; m < KIntM; ++m) {
+              sumPmwkReco[isp][ieta][m][k] += std::pow(1.0, k) * std::pow(pt, m);
+              sumPmwkRecoEffCor[isp][ieta][m][k] += std::pow(w, k) * std::pow(pt, m);
             }
+            sumWkReco[isp][ieta][k] += std::pow(1.0, k);
+            sumWkRecoEffCor[isp][ieta][k] += std::pow(w, k);
           }
-        } // trkslice
+        }
+
+        if (isp == kInclusiveIdx) {
+          histos.fill(HIST("hEtaPhiReco"), vz, sign, pt, eta, phi);
+          histos.fill(HIST("hEtaPhiRecoWtd"), vz, sign, pt, eta, phi, w);
+          histos.fill(HIST("hEtaPhiRecoEffWtd"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
+        } else if (isp == kPiMinusIdx) {
+          histos.fill(HIST("hEtaPhiReco_PiMinus"), vz, sign, pt, eta, phi);
+          histos.fill(HIST("hEtaPhiRecoWtd_PiMinus"), vz, sign, pt, eta, phi, w);
+          histos.fill(HIST("hEtaPhiRecoEffWtd_PiMinus"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
+        } else if (isp == kPiPlusIdx) {
+          histos.fill(HIST("hEtaPhiReco_PiPlus"), vz, sign, pt, eta, phi);
+          histos.fill(HIST("hEtaPhiRecoWtd_PiPlus"), vz, sign, pt, eta, phi, w);
+          histos.fill(HIST("hEtaPhiRecoEffWtd_PiPlus"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
+        } else if (isp == kPiAllIdx) {
+          histos.fill(HIST("hEtaPhiReco_PiAll"), vz, sign, pt, eta, phi);
+          histos.fill(HIST("hEtaPhiRecoWtd_PiAll"), vz, sign, pt, eta, phi, w);
+          histos.fill(HIST("hEtaPhiRecoEffWtd_PiAll"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
+        } else if (isp == kKaMinusIdx) {
+          histos.fill(HIST("hEtaPhiReco_KaMinus"), vz, sign, pt, eta, phi);
+          histos.fill(HIST("hEtaPhiRecoWtd_KaMinus"), vz, sign, pt, eta, phi, w);
+          histos.fill(HIST("hEtaPhiRecoEffWtd_KaMinus"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
+        } else if (isp == kKaPlusIdx) {
+          histos.fill(HIST("hEtaPhiReco_KaPlus"), vz, sign, pt, eta, phi);
+          histos.fill(HIST("hEtaPhiRecoWtd_KaPlus"), vz, sign, pt, eta, phi, w);
+          histos.fill(HIST("hEtaPhiRecoEffWtd_KaPlus"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
+        } else if (isp == kKaAllIdx) {
+          histos.fill(HIST("hEtaPhiReco_KaAll"), vz, sign, pt, eta, phi);
+          histos.fill(HIST("hEtaPhiRecoWtd_KaAll"), vz, sign, pt, eta, phi, w);
+          histos.fill(HIST("hEtaPhiRecoEffWtd_KaAll"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
+        } else if (isp == kPrIdx) {
+          histos.fill(HIST("hEtaPhiReco_Pr"), vz, sign, pt, eta, phi);
+          histos.fill(HIST("hEtaPhiRecoWtd_Pr"), vz, sign, pt, eta, phi, w);
+          histos.fill(HIST("hEtaPhiRecoEffWtd_Pr"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
+        } else if (isp == kAntiPrIdx) {
+          histos.fill(HIST("hEtaPhiReco_AntiPr"), vz, sign, pt, eta, phi);
+          histos.fill(HIST("hEtaPhiRecoWtd_AntiPr"), vz, sign, pt, eta, phi, w);
+          histos.fill(HIST("hEtaPhiRecoEffWtd_AntiPr"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
+        } else if (isp == kPrAllIdx) {
+          histos.fill(HIST("hEtaPhiReco_PrAll"), vz, sign, pt, eta, phi);
+          histos.fill(HIST("hEtaPhiRecoWtd_PrAll"), vz, sign, pt, eta, phi, w);
+          histos.fill(HIST("hEtaPhiRecoEffWtd_PrAll"), vz, sign, pt, eta, phi, (1.0 - fake) / eff);
+        }
+      }
+    } // trkslice
 
         for (int ieta = 0; ieta < KNEta; ++ieta) {
           const int ibx = state.pmeanTruNchEtabinSpbinStep2->GetXaxis()->FindBin(mcCollision.multNTracksPV());