diff --git a/PWGLF/TableProducer/Nuspex/nucleiAntineutronCex.cxx b/PWGLF/TableProducer/Nuspex/nucleiAntineutronCex.cxx index dfb60ad5070..a0e3c431029 100644 --- a/PWGLF/TableProducer/Nuspex/nucleiAntineutronCex.cxx +++ b/PWGLF/TableProducer/Nuspex/nucleiAntineutronCex.cxx @@ -50,7 +50,7 @@ using o2::constants::math::Rad2Deg; struct NucleiAntineutronCex { // Slicing per colision Preslice perMcByColl = aod::mcparticle::mcCollisionId; - + using TracksWCovMc = soa::Join; // === Cut values === @@ -106,11 +106,11 @@ struct NucleiAntineutronCex { histos.add("pPz", "p_{z};p_{z} (GeV/c);Entries", kTH1F, {{100, -10., 10.}}); histos.add("pEta", "Pseudorapidity;#eta;Entries", kTH1F, {{100, -10., 10.}}); histos.add("pP_ITScuts", "Momentum with ITS cuts;|p| (GeV/c);Entries", kTH1F, {{100, 0., 10.}}); - + // Process enum breakdown (secondary antiproton that anchors the SV) histos.add("hProcEnumAP_CEX", "procEnum of secondary #bar{p} (CEX);procEnum;Entries", kTH1I, {{100, -0.5, 99.5}}); - histos.add("hProcEnumAP_BG", "procEnum of secondary #bar{p} (BG);procEnum;Entries", kTH1I, {{100, -0.5, 99.5}}); - + histos.add("hProcEnumAP_BG", "procEnum of secondary #bar{p} (BG);procEnum;Entries", kTH1I, {{100, -0.5, 99.5}}); + // CEX pair from antineutron (MC) histos.add("cexPairMcP", "CEX pair total momentum;|p| (GeV/c);Entries", kTH1F, {{100, 0., 10.}}); histos.add("cexPairMcPt", "CEX pair p_{T};p_{T} (GeV/c);Entries", kTH1F, {{100, 0., 10.}}); @@ -135,8 +135,8 @@ struct NucleiAntineutronCex { histos.add("cexbg_pairmc_vtx", "Background pair vertex;X (cm);Y (cm)", kTH2F, {{200, -60., 60.}, {200, -60., 60.}}); histos.add("cexbg_pairmc_vtxz", "Background secondary vertex Z;Z (cm);Entries", kTH1F, {{200, -60., 60.}}); histos.add("cexbg_pairmc_pITScuts", "Background momentum (ITS cuts);|p| (GeV/c);Entries", kTH1F, {{100, 0., 10.}}); - - //Pi0 events + + // Pi0 events histos.add("cexn_pairmc_p_pi0", "Pair p / antineutron p for CEX + #pi^{0};p/p_{#bar{n}};Entries", kTH1F, {{100, 0., 2.}}); // CEX pair from antineutron (TRK) @@ -366,10 +366,10 @@ struct NucleiAntineutronCex { break; } } - + if (pionPlus || pionMinus) continue; - + // Check for neutral pion at the same secondary vertex bool pion0 = false; for (const auto& particle4 : mcPartsThis) { @@ -728,7 +728,7 @@ struct NucleiAntineutronCex { const TVector3 pv2sv(secX - pvtxX, secY - pvtxY, secZ - pvtxZ); const double pairPointingAngleDeg = pv2sv.Angle(total_trk_pVec) * Rad2Deg; - + const double pvsvThetaDeg = pv2sv.Theta() * Rad2Deg; double pvsvPhiDeg = pv2sv.Phi() * Rad2Deg; @@ -797,14 +797,14 @@ struct NucleiAntineutronCex { histos.fill(HIST("vtxfit_mc_d3D"), d3d); const bool isCex = (motherPdg == -kNeutron); - + // Nature of the process if (isCex) { histos.fill(HIST("hProcEnumAP_CEX"), static_cast(procEnum)); } else { histos.fill(HIST("hProcEnumAP_BG"), static_cast(procEnum)); } - + const float vtxfitDX = secX - antipVx; const float vtxfitDY = secY - antipVy; const float vtxfitDZ = secZ - antipVz; @@ -891,8 +891,7 @@ struct NucleiAntineutronCex { antipTrkItsNSigmaPr, antipTrkItsPidValid, - antipTrkTgl - ); + antipTrkTgl); } } // ==== end DCAFitter2 ====