\pt-differential inclusive \psis production cross section at midrapidity in \pp collisions at \thirteen is compared with ALICE results at forward rapidity  (left) and with the model calculations from ICEM  and NRQCD  (right). Both calculations include a non-prompt contribution from beauty-hadron decays via FONLL . Data points are plotted according to the Lafferty–Wyatt prescription . The lower panels show the ratios of the data points or model calculations to the power-law fit to the midrapidity data.