diff --git a/PWGLF/TableProducer/Nuspex/nucleiAntineutronCex.cxx b/PWGLF/TableProducer/Nuspex/nucleiAntineutronCex.cxx index 32a896db20c..4ba3316a0b7 100644 --- a/PWGLF/TableProducer/Nuspex/nucleiAntineutronCex.cxx +++ b/PWGLF/TableProducer/Nuspex/nucleiAntineutronCex.cxx @@ -49,8 +49,8 @@ struct NucleiAntineutronCex { using TracksWCovMc = soa::Join; // === Cut values === - static constexpr double kIts2MinR = 2.2; // ITS2 min radius [cm] - static constexpr double kIts2MaxR = 39.0; // ITS2 max radius [cm] + static constexpr double kIts2MinR = 4.5; // ITS2 min radius (exluding IB) [cm] + static constexpr double kIts2MaxR = 48.0; // ITS2 max radius [cm] static constexpr double kIts2MaxVz = 39.0; // ITS2 max |vz| [cm] static constexpr double kAccMaxEta = 1.2; // acceptance |eta| static constexpr double kAccMaxVz = 5.3; // acceptance |vz| [cm] @@ -609,7 +609,7 @@ struct NucleiAntineutronCex { histos.fill(HIST("pItsPidValid"), pTrkItsPidValid); histos.fill(HIST("pTgl"), pTrkTgl); } - if (motherPdg == -kNeutron) { + if (motherPdg != -kNeutron) { histos.fill(HIST("pItsNsigmaPr_bg"), pTrkItsNSigmaPr); histos.fill(HIST("pItsPidValid_bg"), pTrkItsPidValid); histos.fill(HIST("pTgl_bg"), pTrkTgl); @@ -715,7 +715,7 @@ struct NucleiAntineutronCex { const float dOpenAngle = trkangleDeg - mcangleDeg; // Closest ITS layer: Radius need to be checked - static const std::array rLayers = {2.2, 2.8, 3.6, 19.6, 24.0, 29.0, 35.0}; + static const std::array rLayers = {2.3, 3.1, 3.9, 24.0, 30.0, 42.0, 48.0}; int16_t svNearestLayerId = -1; float svDeltaRToLayer = 1e9f; for (int i = 0; i < static_cast(rLayers.size()); ++i) {