[PWGCF] FemtoUniverse: Add shuffling in pairs from closed-pair level (#14676)

Co-authored-by: ALICE Action Bot <alibuild@cern.ch>

Author: prchakra

PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackTrackSpherHarMultKtExtended.cxx

@@ -585,40 +585,40 @@ struct femtoUniversePairTaskTrackTrackSpherHarMultKtExtended {
         if (kT < firstRealElement || kT > lastElement)
           continue;
 
+        double rand;
+        rand = randgen->Rndm();
+
+        auto part1 = p1;
+        auto part2 = p2;
+
+        if (rand > 0.5) {
+          part1 = p2;
+          part2 = p1;
+        }
+
         if (twotracksconfigs.confIsCPR.value) {
           if (twotracksconfigs.confCPRFracMax.value) {
-            if (pairCloseRejection.isClosePairAtITS(p1, p2, magFieldTesla, femto_universe_container::EventType::same)) {
+            if (pairCloseRejection.isClosePairAtITS(part1, part2, magFieldTesla, femto_universe_container::EventType::same)) {
               continue;
             }
           } else {
             if (twotracksconfigs.confIsCPRkT) {
-              if (pairCloseRejection.isClosePairkT(p1, p2, femto_universe_container::EventType::same, kT, twotracksconfigs.confCPRDphiAvgOrDist)) {
+              if (pairCloseRejection.isClosePairkT(part1, part2, femto_universe_container::EventType::same, kT, twotracksconfigs.confCPRDphiAvgOrDist)) {
                 continue;
               }
             } else {
-              if (pairCloseRejection.isClosePairFrac(p1, p2, magFieldTesla, femto_universe_container::EventType::same, twotracksconfigs.confCPRDphiAvgOrDist, twotracksconfigs.confCPRDistMax, twotracksconfigs.confCPRFracMax, twotracksconfigs.confCPRDphiAvgOrDist)) {
+              if (pairCloseRejection.isClosePairFrac(part1, part2, magFieldTesla, femto_universe_container::EventType::same, twotracksconfigs.confCPRDphiAvgOrDist, twotracksconfigs.confCPRDistMax, twotracksconfigs.confCPRFracMax, twotracksconfigs.confCPRDphiAvgOrDist)) {
                 continue;
               }
             }
           }
         }
 
         // track cleaning
-        if (!pairCleaner.isCleanPair(p1, p2, parts)) {
+        if (!pairCleaner.isCleanPair(part1, part2, parts)) {
           continue;
         }
 
-        double rand;
-        rand = randgen->Rndm();
-
-        auto part1 = p1;
-        auto part2 = p2;
-
-        if (rand > 0.5) {
-          part1 = p2;
-          part2 = p1;
-        }
-
         std::vector<double> f3d;
         double kv;
         float outsideref = 0.0;
@@ -943,35 +943,34 @@ struct femtoUniversePairTaskTrackTrackSpherHarMultKtExtended {
 
       if (kT < firstRealElement || kT > lastElement)
         continue;
+      double rand;
+      rand = randgen->Rndm();
+      auto part1 = p1;
+      auto part2 = p2;
+
+      if (rand > 0.5) {
+        part1 = p2;
+        part2 = p1;
+      }
 
       if (twotracksconfigs.confIsCPR.value) {
         if (twotracksconfigs.confCPRFracMax.value) {
-          if (pairCloseRejection.isClosePairAtITS(p1, p2, magFieldTesla, femto_universe_container::EventType::mixed)) {
+          if (pairCloseRejection.isClosePairAtITS(part1, part2, magFieldTesla, femto_universe_container::EventType::mixed)) {
             continue;
           }
         } else {
           if (twotracksconfigs.confIsCPRkT) {
-            if (pairCloseRejection.isClosePairkT(p1, p2, femto_universe_container::EventType::mixed, kT, twotracksconfigs.confCPRDphiAvgOrDist)) {
+            if (pairCloseRejection.isClosePairkT(part1, part2, femto_universe_container::EventType::mixed, kT, twotracksconfigs.confCPRDphiAvgOrDist)) {
               continue;
             }
           } else {
-            if (pairCloseRejection.isClosePairFrac(p1, p2, magFieldTesla, femto_universe_container::EventType::mixed, twotracksconfigs.confCPRDphiAvgOrDist, twotracksconfigs.confCPRDistMax, twotracksconfigs.confCPRFracMax, twotracksconfigs.confCPRDphiAvgOrDist)) {
+            if (pairCloseRejection.isClosePairFrac(part1, part2, magFieldTesla, femto_universe_container::EventType::mixed, twotracksconfigs.confCPRDphiAvgOrDist, twotracksconfigs.confCPRDistMax, twotracksconfigs.confCPRFracMax, twotracksconfigs.confCPRDphiAvgOrDist)) {
               continue;
             }
           }
         }
       }
 
-      double rand;
-      rand = randgen->Rndm();
-      auto part1 = p1;
-      auto part2 = p2;
-
-      if (rand > 0.5) {
-        part1 = p2;
-        part2 = p1;
-      }
-
       std::vector<double> f3d;
       double kv;
       float outsideref = 0.0;