Differential $\phi$ meson production cross section $\mathrm{d}^2\sigma/\mathrm{d}y\mathrm{d}\pt$ as a function of $\pt$ at $\sqrt{s}=5.02$ TeV (left) and $\sqrt{s}=13$ TeV (right) in several rapidity intervals. The boxes represent the systematic uncertainties, the error bars the statistical uncertainties Data are fitted with a Levy-Tsallis function.