The PHENIX experiment at the Relativistic Heavy Ion Collider has measured low mass vector meson, $\omega$, $\rho$, and $\phi$, production through the dimuon decay channel at forward rapidity ($1.2<|y|<2.2$) in $p$$+$$p$ collisions at $\sqrt{s}=200$ GeV. The differential cross sections for these mesons are measured as a function of both $p_T$ and rapidity. We also report the integrated differential cross sections over $1<p_T<7$ GeV/$c$ and $1.2<|y|<2.2$: $d\sigma/dy(\omega+\rho\rightarrow\mu\mu) = 80 \pm 6 \mbox{(stat)} \pm 12 \mbox{(syst)}$ nb and $d\sigma/dy(\phi\rightarrow\mu\mu) = 27 \pm 3 \mbox{(stat)} \pm 4 \mbox{(syst)}$ nb. These results are compared with midrapidity measurements and calculations.
Yields and phase space distributions of φ -mesons emitted from p+p (minimum bias trigger), p+Pb (at various centralities) and central Pb+Pb collisions are reported ( E beam =158 A GeV). The decay φ →K + K − was used for identification. The φ / π ratio is found to increase by a factor of 3.0±0.7 from inelastic p+p to central Pb+Pb. Significant enhancement in this ratio is also observed in subclasses of p+p events (characterized by high charged-particle multiplicity) as well as in the forward hemisphere of central p+Pb collisions. In Pb+Pb no shift or significant broadening of the φ -peak is seen.