Invariant differential cross section of \piz (left) and \et (right) versus transverse momentum for \pp collisions at \sT. The data are parametrized with a modified \ac{TCM} model (see \Eq{eq:TwoComponentModel}) and compared to predictions from PYTHIA8 Monash, PYTHIA8 Ropes, EPOS LHC and predictions from NLO pQCD calculations using recent PDFs and FFs. Ratio plots of the data and model calculations to the modified \ac{TCM} fit of the data are shown in the lower panels. Statistical error bars are represented by vertical bars, and systematic uncertainties are shown as boxes.