Please consider the following formatting changes

Author: alibuild

PWGLF/Tasks/Nuspex/antinucleiInJets.cxx

@@ -498,43 +498,37 @@ struct AntinucleiInJets {
       4, 3, 6, 3, 7, 5, 4, 6, 5, 7,
       6, 5, 3, 5, 4, 3, 6, 6, 4, 7,
       3, 4, 3, 5, 7, 6, 6, 4, 3, 5,
-      4, 7, 3, 6, 4, 5, 6, 3, 7, 5
-    };
+      4, 7, 3, 6, 4, 5, 6, 3, 7, 5};
     static std::vector<int> minTpcNcrossedRowsSyst = {
       72, 99, 87, 66, 83, 71, 78, 74, 84, 94,
       89, 85, 66, 80, 95, 74, 76, 77, 67, 88,
       61, 100, 97, 78, 81, 91, 64, 63, 93, 69,
       87, 99, 65, 62, 96, 60, 75, 92, 82, 89,
-      98, 70, 90, 85, 73, 79, 68, 86, 97, 61
-    };
+      98, 70, 90, 85, 73, 79, 68, 86, 97, 61};
     static std::vector<double> maxChiSquareTpcSyst = {
       4.28, 4.81, 4.43, 4.02, 3.38, 3.58, 3.11, 4.17, 3.51, 4.53,
       4.90, 3.07, 3.20, 4.86, 4.62, 3.91, 3.98, 4.38, 3.66, 3.84,
       3.03, 3.14, 4.96, 4.07, 4.75, 4.32, 3.31, 3.78, 4.11, 3.23,
       3.87, 3.70, 4.99, 4.48, 4.69, 4.25, 3.93, 3.45, 4.58, 3.35,
-      3.18, 3.60, 4.21, 3.75, 4.64, 4.35, 3.26, 3.42, 4.15, 3.09
-    };
+      3.18, 3.60, 4.21, 3.75, 4.64, 4.35, 3.26, 3.42, 4.15, 3.09};
     static std::vector<double> maxChiSquareItsSyst = {
       42.84, 48.66, 39.27, 34.09, 43.73, 36.98, 30.23, 49.11, 37.67, 35.10,
       44.55, 46.79, 38.92, 40.66, 47.14, 33.46, 30.88, 41.32, 45.90, 39.68,
       31.42, 32.71, 43.17, 36.04, 49.80, 33.95, 31.89, 38.37, 48.08, 35.87,
       47.61, 44.02, 32.15, 46.21, 34.75, 40.17, 37.14, 30.55, 45.42, 42.30,
-      41.79, 33.21, 39.12, 47.98, 36.52, 31.58, 49.44, 38.02, 35.56, 43.49
-    };
+      41.79, 33.21, 39.12, 47.98, 36.52, 31.58, 49.44, 38.02, 35.56, 43.49};
     static std::vector<double> minEtaSyst = {
       -0.804, -0.844, -0.751, -0.784, -0.819, -0.823, -0.768, -0.781, -0.845, -0.787,
       -0.758, -0.828, -0.776, -0.842, -0.808, -0.763, -0.849, -0.770, -0.799, -0.754,
       -0.825, -0.847, -0.806, -0.783, -0.796, -0.835, -0.777, -0.752, -0.838, -0.813,
       -0.785, -0.802, -0.795, -0.846, -0.780, -0.829, -0.817, -0.773, -0.765, -0.789,
-      -0.800, -0.839, -0.758, -0.820, -0.833, -0.762, -0.792, -0.809, -0.827, -0.751
-    };
+      -0.800, -0.839, -0.758, -0.820, -0.833, -0.762, -0.792, -0.809, -0.827, -0.751};
     static std::vector<double> maxEtaSyst = {
       0.804, 0.844, 0.751, 0.784, 0.819, 0.823, 0.768, 0.781, 0.845, 0.787,
       0.758, 0.828, 0.776, 0.842, 0.808, 0.763, 0.849, 0.770, 0.799, 0.754,
       0.825, 0.847, 0.806, 0.783, 0.796, 0.835, 0.777, 0.752, 0.838, 0.813,
       0.785, 0.802, 0.795, 0.846, 0.780, 0.829, 0.817, 0.773, 0.765, 0.789,
-      0.800, 0.839, 0.758, 0.820, 0.833, 0.762, 0.792, 0.809, 0.827, 0.751
-    };
+      0.800, 0.839, 0.758, 0.820, 0.833, 0.762, 0.792, 0.809, 0.827, 0.751};
 
     // Track Selection
     if (requirePvContributor && !(track.isPVContributor()))
@@ -555,7 +549,7 @@ struct AntinucleiInJets {
       return false;
     if (track.pt() < minPt)
       return false;
-     
+
     return true;
   }
 
@@ -596,7 +590,7 @@ struct AntinucleiInJets {
     // Retrieve the bunch crossing information with timestamps from the collision
     auto bc = collision.template bc_as<aod::BCsWithTimestamps>();
     initCCDB(bc);
-      
+
     // If skimmed processing is enabled, apply Zorro trigger selection
     if (cfgSkimmedProcessing && !zorro.isSelected(collision.template bc_as<aod::BCsWithTimestamps>().globalBC())) {
       return;
@@ -942,7 +936,7 @@ struct AntinucleiInJets {
       fastjet::PseudoJet jetMinusBkg = backgroundSub.doRhoAreaSub(jetForSub, rhoPerp, rhoMPerp);
       if (jetMinusBkg.pt() < minJetPt)
         continue;
-        
+
       // Apply area cut if required
       double normalizedJetArea = jet.area() / (PI * rJet * rJet);
       if (applyAreaCut && normalizedJetArea > maxNormalizedJetArea)
@@ -1334,7 +1328,7 @@ struct AntinucleiInJets {
         for (const auto& particle : jetConstituents) {
           if (particle.user_index() != PDG_t::kProtonBar)
             continue;
-            
+
           if (particle.eta() < minEta || particle.eta() > maxEta)
             continue;
 
@@ -1374,7 +1368,7 @@ struct AntinucleiInJets {
           // Reject tracks that lie outside the maxConeRadius from both UE axes
           if (deltaRUe1 > maxConeRadius && deltaRUe2 > maxConeRadius)
             continue;
-            
+
           // Select antiprotons based on PDG
           if (particle.pdgCode() != PDG_t::kProtonBar)
             continue;
@@ -1405,7 +1399,7 @@ struct AntinucleiInJets {
 
       // Event counter: after event selection
       registryMC.fill(HIST("recEvents"), 1.5);
-        
+
       // Reject events near the ITS Read-Out Frame border
       if (rejectITSROFBorder && !collision.selection_bit(o2::aod::evsel::kNoITSROFrameBorder))
         continue;
@@ -1435,7 +1429,7 @@ struct AntinucleiInJets {
       if (requireIsVertexTOFmatched && !collision.selection_bit(o2::aod::evsel::kIsVertexTOFmatched))
         continue;
       registryMC.fill(HIST("recEvents"), 7.5);
-  
+
       // Loop over reconstructed tracks
       int id(-1);
       std::vector<fastjet::PseudoJet> fjParticles;
@@ -1569,7 +1563,7 @@ struct AntinucleiInJets {
           if (!track.has_mcParticle())
             continue;
           const auto mcparticle = track.mcParticle();
-            
+
           // Antiproton selection
           if (track.sign() > 0 || mcparticle.pdgCode() != PDG_t::kProtonBar)
             continue;
@@ -1580,7 +1574,7 @@ struct AntinucleiInJets {
           double pt = track.pt();
           double dcaxy = track.dcaXY();
           double dcaz = track.dcaZ();
-            
+
           // Fill DCA templates
           if (mcparticle.pdgCode() == PDG_t::kProtonBar && std::fabs(dcaz) < maxDcaz) {
             if (mcparticle.isPhysicalPrimary()) {
@@ -1621,11 +1615,11 @@ struct AntinucleiInJets {
 
           // Fill inclusive antiproton spectrum
           registryMC.fill(HIST("antiproton_incl_ue"), pt);
-          
+
           // Select physical primary antiprotons
           if (!mcparticle.isPhysicalPrimary())
             continue;
-            
+
           // Fill antiproton spectrum for physical primaries
           registryMC.fill(HIST("antiproton_prim_ue"), pt);
 
@@ -1653,39 +1647,39 @@ struct AntinucleiInJets {
 
     // Apply standard event selection
     if (!collision.sel8() || std::fabs(collision.posZ()) > zVtx)
-        return;
+      return;
 
     // Event counter: after event selection
     registryData.fill(HIST("number_of_events_data_syst"), 1.5);
 
     // Reject events near the ITS Read-Out Frame border
     if (rejectITSROFBorder && !collision.selection_bit(o2::aod::evsel::kNoITSROFrameBorder))
-        return;
+      return;
     registryData.fill(HIST("number_of_events_data_syst"), 2.5);
 
     // Reject events at the Time Frame border
     if (rejectTFBorder && !collision.selection_bit(o2::aod::evsel::kNoTimeFrameBorder))
-        return;
+      return;
     registryData.fill(HIST("number_of_events_data_syst"), 3.5);
 
     // Require at least one ITS-TPC matched track
     if (requireVtxITSTPC && !collision.selection_bit(o2::aod::evsel::kIsVertexITSTPC))
-        return;
+      return;
     registryData.fill(HIST("number_of_events_data_syst"), 4.5);
 
     // Reject events with same-bunch pileup
     if (rejectSameBunchPileup && !collision.selection_bit(o2::aod::evsel::kNoSameBunchPileup))
-        return;
+      return;
     registryData.fill(HIST("number_of_events_data_syst"), 5.5);
 
     // Require consistent FT0 vs PV z-vertex
     if (requireIsGoodZvtxFT0VsPV && !collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV))
-        return;
+      return;
     registryData.fill(HIST("number_of_events_data_syst"), 6.5);
 
     // Require TOF match for at least one vertex track
     if (requireIsVertexTOFmatched && !collision.selection_bit(o2::aod::evsel::kIsVertexTOFmatched))
-        return;
+      return;
     registryData.fill(HIST("number_of_events_data_syst"), 7.5);
 
     // Cut settings
@@ -1694,43 +1688,37 @@ struct AntinucleiInJets {
       0.079, 0.043, 0.067, 0.059, 0.032, 0.070, 0.048, 0.077, 0.062, 0.034,
       0.057, 0.055, 0.073, 0.038, 0.050, 0.075, 0.041, 0.061, 0.033, 0.069,
       0.035, 0.044, 0.076, 0.049, 0.037, 0.054, 0.072, 0.046, 0.058, 0.040,
-      0.068, 0.042, 0.056, 0.039, 0.047, 0.065, 0.051, 0.053, 0.063, 0.030
-    };
+      0.068, 0.042, 0.056, 0.039, 0.047, 0.065, 0.051, 0.053, 0.063, 0.030};
     static std::vector<double> maxDcazSyst = {
       0.064, 0.047, 0.032, 0.076, 0.039, 0.058, 0.043, 0.069, 0.050, 0.035,
       0.074, 0.061, 0.045, 0.033, 0.068, 0.055, 0.037, 0.071, 0.042, 0.053,
       0.077, 0.038, 0.065, 0.049, 0.036, 0.059, 0.044, 0.067, 0.041, 0.034,
       0.073, 0.052, 0.040, 0.063, 0.046, 0.031, 0.070, 0.054, 0.037, 0.062,
-      0.048, 0.035, 0.075, 0.051, 0.039, 0.066, 0.043, 0.060, 0.032, 0.056
-    };
+      0.048, 0.035, 0.075, 0.051, 0.039, 0.066, 0.043, 0.060, 0.032, 0.056};
     static std::vector<double> nSigmaItsMinSyst = {
       -2.91, -2.77, -3.03, -3.40, -2.69, -3.28, -2.96, -3.11, -3.36, -3.14,
       -2.99, -2.75, -3.17, -2.64, -2.72, -3.44, -2.87, -2.95, -3.00, -2.66,
       -2.93, -3.31, -3.05, -3.12, -3.21, -3.01, -2.89, -2.73, -3.26, -2.97,
       -2.81, -3.33, -2.68, -3.30, -2.78, -3.39, -2.84, -3.45, -2.92, -3.15,
-      -3.22, -2.58, -3.07, -2.86, -3.10, -2.76, -2.94, -3.25, -3.04, -2.82
-    };
+      -3.22, -2.58, -3.07, -2.86, -3.10, -2.76, -2.94, -3.25, -3.04, -2.82};
     static std::vector<double> nSigmaItsMaxSyst = {
       2.91, 2.77, 3.03, 3.40, 2.69, 3.28, 2.96, 3.11, 3.36, 3.14,
       2.99, 2.75, 3.17, 2.64, 2.72, 3.44, 2.87, 2.95, 3.00, 2.66,
       2.93, 3.31, 3.05, 3.12, 3.21, 3.01, 2.89, 2.73, 3.26, 2.97,
       2.81, 3.33, 2.68, 3.30, 2.78, 3.39, 2.84, 3.45, 2.92, 3.15,
-      3.22, 2.58, 3.07, 2.86, 3.10, 2.76, 2.94, 3.25, 3.04, 2.82
-    };
+      3.22, 2.58, 3.07, 2.86, 3.10, 2.76, 2.94, 3.25, 3.04, 2.82};
     static std::vector<double> minNsigmaTpcSyst = {
       -3.18, -2.86, -3.12, -2.91, -3.49, -2.57, -3.33, -2.98, -3.46, -2.70,
       -3.01, -2.65, -3.27, -3.40, -2.81, -3.10, -2.55, -3.22, -3.07, -2.77,
       -3.35, -2.68, -3.43, -2.88, -3.04, -2.53, -3.30, -2.79, -3.15, -2.66,
       -3.41, -2.75, -3.26, -2.61, -3.09, -2.54, -3.36, -2.95, -3.20, -2.58,
-      -3.44, -2.83, -3.11, -2.62, -3.28, -2.69, -3.23, -2.73, -3.39, -2.90
-    };
+      -3.44, -2.83, -3.11, -2.62, -3.28, -2.69, -3.23, -2.73, -3.39, -2.90};
     static std::vector<double> maxNsigmaTpcSyst = {
       3.18, 2.86, 3.12, 2.91, 3.49, 2.57, 3.33, 2.98, 3.46, 2.70,
       3.01, 2.65, 3.27, 3.40, 2.81, 3.10, 2.55, 3.22, 3.07, 2.77,
       3.35, 2.68, 3.43, 2.88, 3.04, 2.53, 3.30, 2.79, 3.15, 2.66,
       3.41, 2.75, 3.26, 2.61, 3.09, 2.54, 3.36, 2.95, 3.20, 2.58,
-      3.44, 2.83, 3.11, 2.62, 3.28, 2.69, 3.23, 2.73, 3.39, 2.90
-    };
+      3.44, 2.83, 3.11, 2.62, 3.28, 2.69, 3.23, 2.73, 3.39, 2.90};
 
     // Initialize ITS PID Response object
     o2::aod::ITSResponse itsResponse;
@@ -1743,12 +1731,12 @@ struct AntinucleiInJets {
         continue;
 
       // Loop over different cut settings
-      for (int isyst = 0 ; isyst < nSyst ; isyst++) {
+      for (int isyst = 0; isyst < nSyst; isyst++) {
 
         // Apply track selection
-        if (!passedTrackSelectionSyst(track,isyst))
+        if (!passedTrackSelectionSyst(track, isyst))
           continue;
-                        
+
         // Define variables
         double nsigmaTPCPr = track.tpcNSigmaPr();
         double nsigmaTOFPr = track.tofNSigmaPr();
@@ -1807,57 +1795,49 @@ struct AntinucleiInJets {
       0.079, 0.043, 0.067, 0.059, 0.032, 0.070, 0.048, 0.077, 0.062, 0.034,
       0.057, 0.055, 0.073, 0.038, 0.050, 0.075, 0.041, 0.061, 0.033, 0.069,
       0.035, 0.044, 0.076, 0.049, 0.037, 0.054, 0.072, 0.046, 0.058, 0.040,
-      0.068, 0.042, 0.056, 0.039, 0.047, 0.065, 0.051, 0.053, 0.063, 0.030
-    };
+      0.068, 0.042, 0.056, 0.039, 0.047, 0.065, 0.051, 0.053, 0.063, 0.030};
     static std::vector<double> maxDcazSyst = {
       0.064, 0.047, 0.032, 0.076, 0.039, 0.058, 0.043, 0.069, 0.050, 0.035,
       0.074, 0.061, 0.045, 0.033, 0.068, 0.055, 0.037, 0.071, 0.042, 0.053,
       0.077, 0.038, 0.065, 0.049, 0.036, 0.059, 0.044, 0.067, 0.041, 0.034,
       0.073, 0.052, 0.040, 0.063, 0.046, 0.031, 0.070, 0.054, 0.037, 0.062,
-      0.048, 0.035, 0.075, 0.051, 0.039, 0.066, 0.043, 0.060, 0.032, 0.056
-    };
+      0.048, 0.035, 0.075, 0.051, 0.039, 0.066, 0.043, 0.060, 0.032, 0.056};
     static std::vector<double> nSigmaItsMinSyst = {
       -2.91, -2.77, -3.03, -3.40, -2.69, -3.28, -2.96, -3.11, -3.36, -3.14,
       -2.99, -2.75, -3.17, -2.64, -2.72, -3.44, -2.87, -2.95, -3.00, -2.66,
       -2.93, -3.31, -3.05, -3.12, -3.21, -3.01, -2.89, -2.73, -3.26, -2.97,
       -2.81, -3.33, -2.68, -3.30, -2.78, -3.39, -2.84, -3.45, -2.92, -3.15,
-      -3.22, -2.58, -3.07, -2.86, -3.10, -2.76, -2.94, -3.25, -3.04, -2.82
-    };
+      -3.22, -2.58, -3.07, -2.86, -3.10, -2.76, -2.94, -3.25, -3.04, -2.82};
     static std::vector<double> nSigmaItsMaxSyst = {
       2.91, 2.77, 3.03, 3.40, 2.69, 3.28, 2.96, 3.11, 3.36, 3.14,
       2.99, 2.75, 3.17, 2.64, 2.72, 3.44, 2.87, 2.95, 3.00, 2.66,
       2.93, 3.31, 3.05, 3.12, 3.21, 3.01, 2.89, 2.73, 3.26, 2.97,
       2.81, 3.33, 2.68, 3.30, 2.78, 3.39, 2.84, 3.45, 2.92, 3.15,
-      3.22, 2.58, 3.07, 2.86, 3.10, 2.76, 2.94, 3.25, 3.04, 2.82
-    };
+      3.22, 2.58, 3.07, 2.86, 3.10, 2.76, 2.94, 3.25, 3.04, 2.82};
     static std::vector<double> minNsigmaTpcSyst = {
       -3.18, -2.86, -3.12, -2.91, -3.49, -2.57, -3.33, -2.98, -3.46, -2.70,
       -3.01, -2.65, -3.27, -3.40, -2.81, -3.10, -2.55, -3.22, -3.07, -2.77,
       -3.35, -2.68, -3.43, -2.88, -3.04, -2.53, -3.30, -2.79, -3.15, -2.66,
       -3.41, -2.75, -3.26, -2.61, -3.09, -2.54, -3.36, -2.95, -3.20, -2.58,
-      -3.44, -2.83, -3.11, -2.62, -3.28, -2.69, -3.23, -2.73, -3.39, -2.90
-    };
+      -3.44, -2.83, -3.11, -2.62, -3.28, -2.69, -3.23, -2.73, -3.39, -2.90};
     static std::vector<double> maxNsigmaTpcSyst = {
       3.18, 2.86, 3.12, 2.91, 3.49, 2.57, 3.33, 2.98, 3.46, 2.70,
       3.01, 2.65, 3.27, 3.40, 2.81, 3.10, 2.55, 3.22, 3.07, 2.77,
       3.35, 2.68, 3.43, 2.88, 3.04, 2.53, 3.30, 2.79, 3.15, 2.66,
       3.41, 2.75, 3.26, 2.61, 3.09, 2.54, 3.36, 2.95, 3.20, 2.58,
-      3.44, 2.83, 3.11, 2.62, 3.28, 2.69, 3.23, 2.73, 3.39, 2.90
-    };
+      3.44, 2.83, 3.11, 2.62, 3.28, 2.69, 3.23, 2.73, 3.39, 2.90};
     static std::vector<double> minNsigmaTofSyst = {
       -3.18, -2.86, -3.12, -2.91, -3.49, -2.57, -3.33, -2.98, -3.46, -2.70,
       -3.01, -2.65, -3.27, -3.40, -2.81, -3.10, -2.55, -3.22, -3.07, -2.77,
       -3.35, -2.68, -3.43, -2.88, -3.04, -2.53, -3.30, -2.79, -3.15, -2.66,
       -3.41, -2.75, -3.26, -2.61, -3.09, -2.54, -3.36, -2.95, -3.20, -2.58,
-      -3.44, -2.83, -3.11, -2.62, -3.28, -2.69, -3.23, -2.73, -3.39, -2.90
-    };
+      -3.44, -2.83, -3.11, -2.62, -3.28, -2.69, -3.23, -2.73, -3.39, -2.90};
     static std::vector<double> maxNsigmaTofSyst = {
       3.94, 3.62, 3.83, 3.15, 3.23, 3.49, 3.10, 3.78, 3.54, 3.36,
       3.91, 3.84, 3.72, 3.00, 3.63, 3.13, 3.68, 3.40, 3.97, 3.01,
       3.74, 3.25, 3.89, 3.08, 3.30, 3.48, 3.59, 3.16, 3.47, 3.31,
       3.92, 3.03, 3.43, 3.24, 3.11, 3.86, 3.60, 3.07, 3.21, 3.98,
-      3.14, 3.69, 3.56, 3.12, 3.28, 3.46, 3.34, 3.39, 3.05, 3.76
-    };
+      3.14, 3.69, 3.56, 3.12, 3.28, 3.46, 3.34, 3.39, 3.05, 3.76};
 
     // Loop over generated collisions
     for (const auto& collision : genCollisions) {
@@ -1939,10 +1919,10 @@ struct AntinucleiInJets {
         const auto mcparticle = track.mcParticle();
 
         // Loop over different cut settings
-        for (int isyst = 0 ; isyst < nSyst ; isyst++) {
+        for (int isyst = 0; isyst < nSyst; isyst++) {
 
           // Apply track selection
-          if (!passedTrackSelectionSyst(track,isyst))
+          if (!passedTrackSelectionSyst(track, isyst))
             continue;
 
           // Define variables