Invariant mass distributions after combinatorial background subtraction for \kstar and \phim candidates in the multiplicity class 0--10$\%$ and transverse momentum range 2.2 $\leq$ $p_{\mathrm{T}}$ $<$ 3.0 GeV/$c$ in the rapidity interval \mbox{-0.3 $< y <$ 0} (panels (a) and (b)) and 0 $< y <$ 0.3 (panels (c) and (d)). The \kstar peak is described by a Breit-Wigner function whereas the \phim peak is fitted with a Voigtian function. The residual background is described by a polynomial function of order 2.