ATLAS and CMS have recently reported enhancements in the top-antitop production rate near threshold, a region where non-perturbative QCD dynamics associated with toponium formation become relevant. We investigate how this behaviour is modified in the presence of a neutral pseudoscalar that couples to gluons and top quarks, using an effective description that consistently incorporates perturbative Standard Model and new physics contributions, their interference and non-perturbative threshold effects. We show that the combined effect of those ingredients markedly shapes the viable region of the pseudoscalar parameter space, particularly for narrow resonances with masses close to twice the top mass. While Standard Model threshold effects could explain a sizeable part of the measured enhancements, the current data remain compatible with additional contributions from pseudoscalar interactions.

We present version 1.11 of MadAnalysis 5, which extends the software package in several major ways to improve the handling of efficiency tables, the computation of observables in different reference frames and the calculation of statistical limits and/or significance. We detail how these improvements, whose development was motivated by the desire to implement two Run 2 LHC analyses targeting signatures with soft leptons and missing energy and exhibiting mild excesses (ATLAS-SUSY-2018-16 and ATLAS-SUSY-2019-09), have been implemented by both direct extensions of the code and integrations with third-party software. We then document the implementation and validation of these analyses, demonstrating their utility along with the improved statistics capabilities of MadAnalysis 5 through an investigation of the Next-to-Minimal Supersymmetric Standard Model in the context of a larger set of overlapping excesses in channels with soft leptons/jets and missing transverse energy.