# Figure 7

 Invariant cross sections for neutral meson production shown together with NLO pQCD predictions using PDFs MSTW08 (CTEQ6M5) with FFs DSS14 (AESSS) for $\pi^{0}$ ($\eta$) as well as PYTHIA8.210 calculations, for which two different tunes were available. The data points were fitted using a TCM fit, Eq.~\ref{eq:Bylinkin}, and a Levy-Tsallis fit, Eq.~\ref{eq:Tsallis}.