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 1.58%, 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 [17] event generator with the Monash 2013 tune [46]. The width of the line indicates the statistical uncertainty of the prediction. The blue dashed line shows a NLO calculation [14] 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.