Correction rapidity + compactification of the skimming function

Author: Valerio Di Bella

PWGHF/HFC/TableProducer/correlatorDplusDplusReduced.cxx

@@ -163,14 +163,14 @@ struct HfCorrelatorDplusDplusReduced {
     if (candidate.pt() < ptCutSkimming[0] || candidate.pt() > ptCutSkimming[ptCutSkimming.size() - 1]) {
       return false;
     }
-    for (size_t i = 1; i < ptCutSkimming.size(); i++) {
-      if (candidate.pt() <= ptCutSkimming[i]) {
-        if (hfHelper.invMassDplusToPiKPi(candidate) < massMinimum[i - 1] ||
-            hfHelper.invMassDplusToPiKPi(candidate) > massMaximum[i - 1] ||
-            candidate.cpa() < cosThetaminimum[i - 1] ||
-            candidate.decayLength() < decayLengthminimum[i - 1] ||
-            abs(candidate.nSigTofKa1()) > nSigmaTofmaximum[i - 1] ||
-            abs(candidate.nSigTpcKa1()) > nSigmaTpcmaximum[i - 1]) {
+    for (size_t i = 0; i < ptCutSkimming.size() - 1; i++) {
+      if (candidate.pt() <= ptCutSkimming[i + 1]) {
+        if (hfHelper.invMassDplusToPiKPi(candidate) < massMinimum[i] ||
+            hfHelper.invMassDplusToPiKPi(candidate) > massMaximum[i] ||
+            candidate.cpa() < cosThetaminimum[i] ||
+            candidate.decayLength() < decayLengthminimum[i] ||
+            abs(candidate.nSigTofKa1()) > nSigmaTofmaximum[i] ||
+            abs(candidate.nSigTpcKa1()) > nSigmaTpcmaximum[i]) {
           return false;
         }
         return true;

PWGHF/HFC/Tasks/taskCorrelationDplusDplusReduced.cxx

@@ -37,7 +37,7 @@ using namespace o2::framework::expressions;
 
 struct HfTaskCorrelationDplusDplusReduced {
   Configurable<int> selectionFlagDplus{"selectionFlagDplus", 1, "Selection Flag for Dplus"};
-  Configurable<double> selectionCutRapidity{"selectionCutRapidity", 1, "Selection the cut of rapidity"};
+  Configurable<double> rapidityMax{"rapidityMax", 1, "Selection the cut of rapidity"};
   Configurable<bool> applyMl{"applyMl", false, "Flag to apply ML selections"};
   Configurable<bool> applySkimming{"applySkimming", false, "Flag to apply Skimming selections"};
 
@@ -237,7 +237,7 @@ struct HfTaskCorrelationDplusDplusReduced {
 
     for (const auto& part1 : localMcParticles) {
       for (auto part2 = part1 + 1; part2 != localMcParticles.end(); ++part2) {
-        if (part1.eta() < selectionCutRapidity && part2.eta() < selectionCutRapidity) {
+        if (std::abs(part1.eta()) < rapidityMax && std::abs(part2.eta()) < rapidityMax) {
           registry.fill(HIST("hDltPhiMcGen"), part2.phi() - part1.phi());
         }
       }