@@ -50,7 +50,7 @@ using o2::constants::math::Rad2Deg;
5050struct NucleiAntineutronCex {
5151 // Slicing per colision
5252 Preslice<aod::McParticles> perMcByColl = aod::mcparticle::mcCollisionId;
53-
53+
5454 using TracksWCovMc = soa::Join<aod::TracksIU, aod::TracksExtra, aod::McTrackLabels, o2::aod::TracksCovIU>;
5555
5656 // === Cut values ===
@@ -106,11 +106,11 @@ struct NucleiAntineutronCex {
106106 histos.add (" pPz" , " p_{z};p_{z} (GeV/c);Entries" , kTH1F , {{100 , -10 ., 10 .}});
107107 histos.add (" pEta" , " Pseudorapidity;#eta;Entries" , kTH1F , {{100 , -10 ., 10 .}});
108108 histos.add (" pP_ITScuts" , " Momentum with ITS cuts;|p| (GeV/c);Entries" , kTH1F , {{100 , 0 ., 10 .}});
109-
109+
110110 // Process enum breakdown (secondary antiproton that anchors the SV)
111111 histos.add (" hProcEnumAP_CEX" , " procEnum of secondary #bar{p} (CEX);procEnum;Entries" , kTH1I , {{100 , -0.5 , 99.5 }});
112- histos.add (" hProcEnumAP_BG" , " procEnum of secondary #bar{p} (BG);procEnum;Entries" , kTH1I , {{100 , -0.5 , 99.5 }});
113-
112+ histos.add (" hProcEnumAP_BG" , " procEnum of secondary #bar{p} (BG);procEnum;Entries" , kTH1I , {{100 , -0.5 , 99.5 }});
113+
114114 // CEX pair from antineutron (MC)
115115 histos.add (" cexPairMcP" , " CEX pair total momentum;|p| (GeV/c);Entries" , kTH1F , {{100 , 0 ., 10 .}});
116116 histos.add (" cexPairMcPt" , " CEX pair p_{T};p_{T} (GeV/c);Entries" , kTH1F , {{100 , 0 ., 10 .}});
@@ -135,8 +135,8 @@ struct NucleiAntineutronCex {
135135 histos.add (" cexbg_pairmc_vtx" , " Background pair vertex;X (cm);Y (cm)" , kTH2F , {{200 , -60 ., 60 .}, {200 , -60 ., 60 .}});
136136 histos.add (" cexbg_pairmc_vtxz" , " Background secondary vertex Z;Z (cm);Entries" , kTH1F , {{200 , -60 ., 60 .}});
137137 histos.add (" cexbg_pairmc_pITScuts" , " Background momentum (ITS cuts);|p| (GeV/c);Entries" , kTH1F , {{100 , 0 ., 10 .}});
138-
139- // Pi0 events
138+
139+ // Pi0 events
140140 histos.add (" cexn_pairmc_p_pi0" , " Pair p / antineutron p for CEX + #pi^{0};p/p_{#bar{n}};Entries" , kTH1F , {{100 , 0 ., 2 .}});
141141
142142 // CEX pair from antineutron (TRK)
@@ -366,10 +366,10 @@ struct NucleiAntineutronCex {
366366 break ;
367367 }
368368 }
369-
369+
370370 if (pionPlus || pionMinus)
371371 continue ;
372-
372+
373373 // Check for neutral pion at the same secondary vertex
374374 bool pion0 = false ;
375375 for (const auto & particle4 : mcPartsThis) {
@@ -728,7 +728,7 @@ struct NucleiAntineutronCex {
728728
729729 const TVector3 pv2sv (secX - pvtxX, secY - pvtxY, secZ - pvtxZ);
730730 const double pairPointingAngleDeg = pv2sv.Angle (total_trk_pVec) * Rad2Deg;
731-
731+
732732 const double pvsvThetaDeg = pv2sv.Theta () * Rad2Deg;
733733
734734 double pvsvPhiDeg = pv2sv.Phi () * Rad2Deg;
@@ -797,14 +797,14 @@ struct NucleiAntineutronCex {
797797 histos.fill (HIST (" vtxfit_mc_d3D" ), d3d);
798798
799799 const bool isCex = (motherPdg == -kNeutron );
800-
800+
801801 // Nature of the process
802802 if (isCex) {
803803 histos.fill (HIST (" hProcEnumAP_CEX" ), static_cast <int >(procEnum));
804804 } else {
805805 histos.fill (HIST (" hProcEnumAP_BG" ), static_cast <int >(procEnum));
806806 }
807-
807+
808808 const float vtxfitDX = secX - antipVx;
809809 const float vtxfitDY = secY - antipVy;
810810 const float vtxfitDZ = secZ - antipVz;
@@ -891,8 +891,7 @@ struct NucleiAntineutronCex {
891891
892892 antipTrkItsNSigmaPr,
893893 antipTrkItsPidValid,
894- antipTrkTgl
895- );
894+ antipTrkTgl);
896895 }
897896 }
898897 // ==== end DCAFitter2 ====
0 commit comments