@@ -2141,11 +2141,11 @@ struct AntinucleiInJets {
21412141
21422142 // Multiplicity percentile
21432143 const float multiplicity = collision.centFT0M ();
2144-
2144+
21452145 // pt/A bins
21462146 std::vector<double > ptOverAbins = {0.3 , 0.4 , 0.5 , 0.6 , 0.7 , 0.8 , 0.9 , 1.0 };
21472147 const int nBins = ptOverAbins.size () - 1 ;
2148-
2148+
21492149 // Particle counters
21502150 std::vector<int > nAntiprotonFullEvent (nBins, 0 );
21512151 std::vector<int > nAntideuteronFullEvent (nBins, 0 );
@@ -2223,11 +2223,11 @@ struct AntinucleiInJets {
22232223 int netDeuteronFullEvent = nTotDeuteronFullEvent - nTotAntideuteronFullEvent;
22242224 registryCorr.fill (HIST (" rho_fullEvent" ), nTotAntideuteronFullEvent, nTotAntiprotonFullEvent, multiplicity);
22252225 registryCorr.fill (HIST (" rho_netP_netD_fullEvent" ), netDeuteronFullEvent, netProtonFullEvent);
2226-
2226+
22272227 // Fill efficiency histograms
22282228 for (int i = 0 ; i < nBins; i++) {
22292229 double ptAcenter_i = 0.5 * (ptOverAbins[i] + ptOverAbins[i + 1 ]);
2230-
2230+
22312231 registryCorr.fill (HIST (" q1d_fullEvent" ), nAntideuteronFullEvent[i], ptAcenter_i);
22322232 registryCorr.fill (HIST (" q1p_fullEvent" ), nAntiprotonFullEvent[i], ptAcenter_i);
22332233 for (int j = 0 ; j < nBins; j++) {
@@ -2238,9 +2238,6 @@ struct AntinucleiInJets {
22382238 }
22392239 }
22402240
2241- // ==== end of processing of full event
2242-
2243-
22442241 // Loop over reconstructed tracks (refactoring: this part can be incorporated above)
22452242 int id (-1 );
22462243 std::vector<fastjet::PseudoJet> fjParticles;
@@ -2320,7 +2317,7 @@ struct AntinucleiInJets {
23202317 // Apply DCA selections
23212318 if (std::fabs (track.dcaXY ()) > maxDcaxy || std::fabs (track.dcaZ ()) > maxDcaz)
23222319 continue ;
2323-
2320+
23242321 // Particle identification using the ITS cluster size
23252322 bool passedItsPidProt (true ), passedItsPidDeut (true );
23262323 double nSigmaITSprot = static_cast <double >(itsResponse.nSigmaITS <o2::track::PID::Proton>(track));
@@ -2373,13 +2370,13 @@ struct AntinucleiInJets {
23732370 }
23742371 }
23752372 } // end of loop over constituents
2376-
2373+
23772374 // Fill correlation histograms
23782375 int netProtonJet = nTotProtonJet - nTotAntiprotonJet;
23792376 int netDeuteronJet = nTotDeuteronJet - nTotAntideuteronJet;
23802377 registryCorr.fill (HIST (" rho_jet" ), nTotAntideuteronJet, nTotAntiprotonJet, multiplicity);
23812378 registryCorr.fill (HIST (" rho_netP_netD_jet" ), netDeuteronJet, netProtonJet);
2382-
2379+
23832380 // Fill efficiency histograms
23842381 for (int i = 0 ; i < nBins; i++) {
23852382 double ptAcenter_i = 0.5 * (ptOverAbins[i] + ptOverAbins[i + 1 ]);
@@ -2393,8 +2390,7 @@ struct AntinucleiInJets {
23932390 registryCorr.fill (HIST (" q1d_q1p_jet" ), ptAcenter_i, ptAcenter_j, nAntideuteronJet[i] * nAntiprotonJet[j]);
23942391 }
23952392 }
2396-
2397-
2393+
23982394 // Particle counters
23992395 std::vector<int > nAntiprotonUE (nBins, 0 );
24002396 std::vector<int > nAntideuteronUE (nBins, 0 );
@@ -2431,7 +2427,7 @@ struct AntinucleiInJets {
24312427 // Reject tracks that lie outside the maxConeRadius from both UE axes
24322428 if (deltaRUe1 > maxConeRadius && deltaRUe2 > maxConeRadius)
24332429 continue ;
2434-
2430+
24352431 // Particle identification using the ITS cluster size
24362432 bool passedItsPidProt (true ), passedItsPidDeut (true );
24372433 double nSigmaITSprot = static_cast <double >(itsResponse.nSigmaITS <o2::track::PID::Proton>(track));
@@ -2484,17 +2480,17 @@ struct AntinucleiInJets {
24842480 }
24852481 }
24862482 }
2487-
2483+
24882484 // Fill correlation histograms
24892485 int netProtonUE = nTotProtonUE - nTotAntiprotonUE;
24902486 int netDeuteronUE = nTotDeuteronUE - nTotAntideuteronUE;
24912487 registryCorr.fill (HIST (" rho_ue" ), nTotAntideuteronUE, nTotAntiprotonUE, multiplicity);
24922488 registryCorr.fill (HIST (" rho_netP_netD_ue" ), netDeuteronUE, netProtonUE);
2493-
2489+
24942490 // Fill efficiency histograms
24952491 for (int i = 0 ; i < nBins; i++) {
24962492 double ptAcenter_i = 0.5 * (ptOverAbins[i] + ptOverAbins[i + 1 ]);
2497-
2493+
24982494 registryCorr.fill (HIST (" q1d_ue" ), nAntideuteronUE[i], ptAcenter_i);
24992495 registryCorr.fill (HIST (" q1p_ue" ), nAntiprotonUE[i], ptAcenter_i);
25002496 for (int j = 0 ; j < nBins; j++) {
0 commit comments