@@ -375,21 +375,21 @@ struct AntinucleiInJets {
375375 registryCorr.add (" rho_netP_netD_ue" , " rho_netP_netD_ue" , HistType::kTH2F , {nAntideuteronsAxis, nAntiprotonsAxis});
376376 registryCorr.add (" rho_netP_netD_fullEvent" , " rho_netP_netD_fullEvent" , HistType::kTH2F , {nAntideuteronsAxis, nAntiprotonsAxis});
377377
378- // Efficiency histograms jet
378+ // Efficiency histograms jet
379379 registryCorr.add (" q1d_jet" , " q1d_jet" , HistType::kTH2F , {nAntideuteronsAxis, ptPerNucleonAxis});
380380 registryCorr.add (" q1p_jet" , " q1p_jet" , HistType::kTH2F , {nAntiprotonsAxis, ptPerNucleonAxis});
381381 registryCorr.add (" q1d_square_jet" , " q1d_square_jet" , HistType::kTH3F , {ptPerNucleonAxis, ptPerNucleonAxis, nBarD2Axis});
382382 registryCorr.add (" q1p_square_jet" , " q1p_square_jet" , HistType::kTH3F , {ptPerNucleonAxis, ptPerNucleonAxis, nBarP2Axis});
383383 registryCorr.add (" q1d_q1p_jet" , " q1d_q1p_jet" , HistType::kTH3F , {ptPerNucleonAxis, ptPerNucleonAxis, nBarDnBarPAxis});
384384
385- // Efficiency histograms UE
385+ // Efficiency histograms UE
386386 registryCorr.add (" q1d_ue" , " q1d_ue" , HistType::kTH2F , {nAntideuteronsAxis, ptPerNucleonAxis});
387387 registryCorr.add (" q1p_ue" , " q1p_ue" , HistType::kTH2F , {nAntiprotonsAxis, ptPerNucleonAxis});
388388 registryCorr.add (" q1d_square_ue" , " q1d_square_ue" , HistType::kTH3F , {ptPerNucleonAxis, ptPerNucleonAxis, nBarD2Axis});
389389 registryCorr.add (" q1p_square_ue" , " q1p_square_ue" , HistType::kTH3F , {ptPerNucleonAxis, ptPerNucleonAxis, nBarP2Axis});
390390 registryCorr.add (" q1d_q1p_ue" , " q1d_q1p_ue" , HistType::kTH3F , {ptPerNucleonAxis, ptPerNucleonAxis, nBarDnBarPAxis});
391391
392- // Efficiency histograms full event
392+ // Efficiency histograms full event
393393 registryCorr.add (" q1d_fullEvent" , " q1d_fullEvent" , HistType::kTH2F , {nAntideuteronsAxis, ptPerNucleonAxis});
394394 registryCorr.add (" q1p_fullEvent" , " q1p_fullEvent" , HistType::kTH2F , {nAntiprotonsAxis, ptPerNucleonAxis});
395395 registryCorr.add (" q1d_square_fullEvent" , " q1d_square_fullEvent" , HistType::kTH3F , {ptPerNucleonAxis, ptPerNucleonAxis, nBarD2Axis});
@@ -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 );
@@ -2182,14 +2182,12 @@ struct AntinucleiInJets {
21822182 if (track.pt () < ptOverAbins[0 ] || track.pt () >= ptOverAbins[nBins]) {
21832183 continue ;
21842184 }
2185- }
2186- else if (isDeuteron (track) && passedItsPidDeut) {
2185+ } else if (isDeuteron (track) && passedItsPidDeut) {
21872186 double ptPerNucleon = 0.5 * track.pt ();
21882187 if (ptPerNucleon < ptOverAbins[0 ] || ptPerNucleon >= ptOverAbins[nBins]) {
21892188 continue ;
21902189 }
2191- }
2192- else {
2190+ } else {
21932191 continue ;
21942192 }
21952193
@@ -2222,12 +2220,12 @@ struct AntinucleiInJets {
22222220 int netProtonFullEvent = nTotProtonFullEvent - nTotAntiprotonFullEvent;
22232221 int netDeuteronFullEvent = nTotDeuteronFullEvent - nTotAntideuteronFullEvent;
22242222 registryCorr.fill (HIST (" rho_fullEvent" ), nTotAntideuteronFullEvent, nTotAntiprotonFullEvent, multiplicity);
2225- registryCorr.fill (HIST (" rho_netP_netD_fullEvent" ), netDeuteronFullEvent, netProtonFullEvent);
2226-
2223+ registryCorr.fill (HIST (" rho_netP_netD_fullEvent" ), netDeuteronFullEvent, netProtonFullEvent);
2224+
22272225 // Fill efficiency histograms
22282226 for (int i = 0 ; i < nBins; i++) {
22292227 double ptAcenter_i = 0.5 * (ptOverAbins[i] + ptOverAbins[i + 1 ]);
2230-
2228+
22312229 registryCorr.fill (HIST (" q1d_fullEvent" ), nAntideuteronFullEvent[i], ptAcenter_i);
22322230 registryCorr.fill (HIST (" q1p_fullEvent" ), nAntiprotonFullEvent[i], ptAcenter_i);
22332231 for (int j = 0 ; j < nBins; j++) {
@@ -2239,8 +2237,7 @@ struct AntinucleiInJets {
22392237 }
22402238
22412239 // ==== end of processing of full event
2242-
2243-
2240+
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));
@@ -2338,14 +2335,12 @@ struct AntinucleiInJets {
23382335 if (track.pt () < ptOverAbins[0 ] || track.pt () >= ptOverAbins[nBins]) {
23392336 continue ;
23402337 }
2341- }
2342- else if (isDeuteron (track) && passedItsPidDeut) {
2338+ } else if (isDeuteron (track) && passedItsPidDeut) {
23432339 double ptPerNucleon = 0.5 * track.pt ();
23442340 if (ptPerNucleon < ptOverAbins[0 ] || ptPerNucleon >= ptOverAbins[nBins]) {
23452341 continue ;
23462342 }
2347- }
2348- else {
2343+ } else {
23492344 continue ;
23502345 }
23512346
@@ -2373,17 +2368,17 @@ struct AntinucleiInJets {
23732368 }
23742369 }
23752370 } // end of loop over constituents
2376-
2371+
23772372 // Fill correlation histograms
23782373 int netProtonJet = nTotProtonJet - nTotAntiprotonJet;
23792374 int netDeuteronJet = nTotDeuteronJet - nTotAntideuteronJet;
23802375 registryCorr.fill (HIST (" rho_jet" ), nTotAntideuteronJet, nTotAntiprotonJet, multiplicity);
2381- registryCorr.fill (HIST (" rho_netP_netD_jet" ), netDeuteronJet, netProtonJet);
2382-
2376+ registryCorr.fill (HIST (" rho_netP_netD_jet" ), netDeuteronJet, netProtonJet);
2377+
23832378 // Fill efficiency histograms
23842379 for (int i = 0 ; i < nBins; i++) {
23852380 double ptAcenter_i = 0.5 * (ptOverAbins[i] + ptOverAbins[i + 1 ]);
2386-
2381+
23872382 registryCorr.fill (HIST (" q1d_jet" ), nAntideuteronJet[i], ptAcenter_i);
23882383 registryCorr.fill (HIST (" q1p_jet" ), nAntiprotonJet[i], ptAcenter_i);
23892384 for (int j = 0 ; j < nBins; j++) {
@@ -2393,8 +2388,7 @@ struct AntinucleiInJets {
23932388 registryCorr.fill (HIST (" q1d_q1p_jet" ), ptAcenter_i, ptAcenter_j, nAntideuteronJet[i] * nAntiprotonJet[j]);
23942389 }
23952390 }
2396-
2397-
2391+
23982392 // Particle counters
23992393 std::vector<int > nAntiprotonUE (nBins, 0 );
24002394 std::vector<int > nAntideuteronUE (nBins, 0 );
@@ -2431,7 +2425,7 @@ struct AntinucleiInJets {
24312425 // Reject tracks that lie outside the maxConeRadius from both UE axes
24322426 if (deltaRUe1 > maxConeRadius && deltaRUe2 > maxConeRadius)
24332427 continue ;
2434-
2428+
24352429 // Particle identification using the ITS cluster size
24362430 bool passedItsPidProt (true ), passedItsPidDeut (true );
24372431 double nSigmaITSprot = static_cast <double >(itsResponse.nSigmaITS <o2::track::PID ::Proton>(track));
@@ -2449,14 +2443,12 @@ struct AntinucleiInJets {
24492443 if (track.pt () < ptOverAbins[0 ] || track.pt () >= ptOverAbins[nBins]) {
24502444 continue ;
24512445 }
2452- }
2453- else if (isDeuteron (track) && passedItsPidDeut) {
2446+ } else if (isDeuteron (track) && passedItsPidDeut) {
24542447 double ptPerNucleon = 0.5 * track.pt ();
24552448 if (ptPerNucleon < ptOverAbins[0 ] || ptPerNucleon >= ptOverAbins[nBins]) {
24562449 continue ;
24572450 }
2458- }
2459- else {
2451+ } else {
24602452 continue ;
24612453 }
24622454
@@ -2484,17 +2476,17 @@ struct AntinucleiInJets {
24842476 }
24852477 }
24862478 }
2487-
2479+
24882480 // Fill correlation histograms
24892481 int netProtonUE = nTotProtonUE - nTotAntiprotonUE;
24902482 int netDeuteronUE = nTotDeuteronUE - nTotAntideuteronUE;
24912483 registryCorr.fill (HIST (" rho_ue" ), nTotAntideuteronUE, nTotAntiprotonUE, multiplicity);
2492- registryCorr.fill (HIST (" rho_netP_netD_ue" ), netDeuteronUE, netProtonUE);
2493-
2484+ registryCorr.fill (HIST (" rho_netP_netD_ue" ), netDeuteronUE, netProtonUE);
2485+
24942486 // Fill efficiency histograms
24952487 for (int i = 0 ; i < nBins; i++) {
24962488 double ptAcenter_i = 0.5 * (ptOverAbins[i] + ptOverAbins[i + 1 ]);
2497-
2489+
24982490 registryCorr.fill (HIST (" q1d_ue" ), nAntideuteronUE[i], ptAcenter_i);
24992491 registryCorr.fill (HIST (" q1p_ue" ), nAntiprotonUE[i], ptAcenter_i);
25002492 for (int j = 0 ; j < nBins; j++) {
0 commit comments