Top panels: Transverse momentum spectrum of \kstar (left) and \phim (right) as a function of \pt for the NSD events, measured in the rapidity interval $-$0.5 $ y$ $$ 0 for p--Pb collisions at \snn = 8.16 TeV. The statistical and systematic uncertainties are shown as bars and boxes, respectively. The NSD spectrum is compared with the predictions from EPOS-LHC~, DPMJET~ and HIJING~. Bottom panels: The ratios of \pt spectra from model to data. The shaded bands around unity describe the statistical and systematic uncertainties of the data point.