Update antinucleiInJets.cxx

Author: FrancescaCasillo

PWGLF/Tasks/Nuspex/antinucleiInJets.cxx

@@ -3704,41 +3704,41 @@ struct AntinucleiInJets {
       registryMC.fill(HIST("genEventsCoalescenceCorr"), 2.5);
 
       // Build deuterons
-      for (const auto& proton : protonCandidates) {
-        if (proton.used)
-          continue;
-
-        for (const auto& neutron : neutronCandidates) {
-          if (neutron.used)
-            continue;
-
-          // Physics consistency check
-          if (proton.pdgCode * neutron.pdgCode < 0)
-            continue;
-
-          if (passDeuteronCoalescence(proton, neutron, coalescenceMomentum, mRand)) {
-
-            neutron.used = true;
-            proton.used = true;
-
-            int sign = (proton.pdgCode > 0) ? +1 : -1;
-            int deuteronPdg = sign * o2::constants::physics::Pdg::kDeuteron;
-
-            double pxDeut = proton.px + neutron.px;
-            double pyDeut = proton.py + neutron.py;
-            double pzDeut = proton.pz + neutron.pz;
-            double energyDeut = std::sqrt(pxDeut * pxDeut + pyDeut * pyDeut + pzDeut * pzDeut + massDeut * massDeut);
-            LorentzVector pd(pxDeut, pyDeut, pzDeut, energyDeut);
-            if (pd.Eta() >= minEta && pd.Eta() <= maxEta && (0.5 * pd.Pt()) >= MinPtPerNucleon) {
-              // Store Deuteron
-              finalDeuterons.push_back({pxDeut, pyDeut, pzDeut, deuteronPdg, proton.mcIndex, false});
+        for (size_t iP = 0; iP < protonCandidates.size(); ++iP) {
+            if (protonCandidates[iP].used)
+                continue;
+            
+            for (size_t iN = 0; iN < neutronCandidates.size(); ++iN) {
+                if (neutronCandidates[iN].used)
+                    continue;
+                
+                // Physics consistency check
+                if (protonCandidates[iP].pdgCode * neutronCandidates[iN].pdgCode < 0)
+                    continue;
+                
+                if (passDeuteronCoalescence(protonCandidates[iP], neutronCandidates[iN], coalescenceMomentum, mRand)) {
+                    
+                    neutronCandidates[iN].used = true;
+                    protonCandidates[iP].used = true;
+                    
+                    int sign = (protonCandidates[iP].pdgCode > 0) ? +1 : -1;
+                    int deuteronPdg = sign * o2::constants::physics::Pdg::kDeuteron;
+                    
+                    double pxDeut = protonCandidates[iP].px + neutronCandidates[iN].px;
+                    double pyDeut = protonCandidates[iP].py + neutronCandidates[iN].py;
+                    double pzDeut = protonCandidates[iP].pz + neutronCandidates[iN].pz;
+                    double energyDeut = std::sqrt(pxDeut * pxDeut + pyDeut * pyDeut + pzDeut * pzDeut + massDeut * massDeut);
+                    LorentzVector pd(pxDeut, pyDeut, pzDeut, energyDeut);
+                    if (pd.Eta() >= minEta && pd.Eta() <= maxEta && (0.5 * pd.Pt()) >= MinPtPerNucleon) {
+                        // Store Deuteron
+                        finalDeuterons.push_back({pxDeut, pyDeut, pzDeut, deuteronPdg, protonCandidates[iP].mcIndex, false});
+                    }
+                    
+                    break;
+                }
             }
-
-            break;
-          }
         }
-      }
-
+        
       // Add unused protons to final vectors
       for (const auto& proton : protonCandidates) {
         if (!proton.used) {