The pseudoscalar production cross section via gluon-gluon fusion was obtained using MadGraph5_aMC@NLO 2.6.2 [1]}. We have used an interface with SPheno where masses, decay widths and branching fractions are calculated and linked to MadGraph5 for every single generated point. The relevance of requiring an agreement with QFV observables is emphasized in this plot, where grey points represent scenarios where at least one of such observables was found to deviate from its SM value by more than 2\(\sigma \) . The blue points denote scenarios for which a full 1\(\sigma \)  agreement was found for all QFV observables while the darker shades of blue also survive all remaining theoretical and experimental constraints. As was to be expected, the total signal strength for this channel diminishes as the pseudoscalar mass increases – see in particular the sharp drop-off around \(m_{A_1} \sim 375\)  GeV, that is twice the top mass. Indeed, for larger pseudoscalar
masses, a new decay channel \(A_1 \rightarrow t\bar{t}\)  becomes kinematically allowed and tends to reduce the branching ratio for \(A_1 \rightarrow \tau \bar{\tau }\) . Notice, however, that there is a number of points for lower \(m_{A_1}\)  which can almost be probed by the current CMS bounds.
