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, $\sqrt{s}=8$ TeV and $\sqrt{s}=13$ TeV, measured in the $\mu^+\mu^-$ decay channel in the rapidity interval $2.5< y< 4$ (full symbols) and in the K$^+$K$^-$ channel at midrapidity (open symbols). The boxes represent the systematic uncertainties, the error bars the statistical uncertainties. The data points are fitted with a Levy-Tsallis function. The ratio between the fit function at forward rapidity and the data at midrapidity is plotted in the bottom panels.