The \pt-differential yield of \fzero in \pp collisions at \five is compared with predictions from the HERWIG 7.2 event generator and with a coalescence calculation based on the AMPT model . The statistical and systematic uncertainties on data (full black markers) are shown as bars and boxes, respectively. The middle and bottom panels show the model to data ratios. The grey boxes at unity represent the sum in quadrature of the statistical and systematic uncertainty on the data. The ratios of $\rm{u\bar{u}s\bar{s}}$ tetraquark and $\rm{K^+K^-}$ molecule configurations from AMPT model predictions to data are multiplied by 100 to improve visibility. In all three panels, the uncertainties associated with the models are statistical ones.