Invariant cross section of $p+p\rightarrow\omega+X$ production at midrapidity in pp collisions at $\sqrt{s}= 13$ TeV compared to theoretical predictions. The cross section is parametrized using a Levy-Tsallis function (dashed grey line), where the fit parameters are given in Table 2. The systematic uncertainty of the measurement is denoted by boxes, excluding the normalization uncertainty of 2.2%, which is shown separately as a grey box in the bottom panel of the figure. To account for the finite width of each $p_{\rm T}$-interval, the data points are shifted according to the underlying spectrum, as described by the fit. The red line represents the theoretical prediction obtained using the PYTHIA8.2  event generator with the Monash 2013 tune . The width of the line indicates the statistical uncertainty of the prediction. The blue dashed line shows a NLO calculation  incorporating $\omega$ fragmentation based on a broken SU(3) model, where all scales are chosen to be $\mu = p_{\rm T}$ and the blue shaded band denotes the scale uncertainty The bottom panel shows the ratio of the data, as well as the theoretical predictions to the Levy-Tsallis parametrization.