Measurement of neutral mesons in<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mi>p</mml:mi><mml:mo>+</mml:mo><mml:mi>p</mml:mi></mml:math>collisions at<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:msqrt><mml:mi>s</mml:mi></mml:msqrt><mml:mo>=</mml:mo><mml:mn>200</mml:mn><mml:mtext> </mml:mtext><mml:mtext> </mml:mtext><mml:mi>GeV</mml:mi></mml:math>and scaling properties of hadron production

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the invariant differential cross section for production of ${K}_{S}^{0}$, $\ensuremath{\omega}$, ${\ensuremath{\eta}}^{\ensuremath{'}}$, and $\ensuremath{\phi}$ mesons in $p+p$ collisions at $\sqrt{s}=200\text{ }\text{ }\mathrm{GeV}$. Measurements of $\ensuremath{\omega}$ and $\ensuremath{\phi}$ production in different decay channels give consistent results. New results for the $\ensuremath{\omega}$ are in agreement with previously published data and extend the measured ${p}_{T}$ coverage. The spectral shapes of all hadron transverse momentum distributions measured by PHENIX are well described by a Tsallis distribution functional form with only two parameters, $n$ and $T$, determining the high-${p}_{T}$ and characterizing the low-${p}_{T}$ regions of the spectra, respectively. The values of these parameters are very similar for all analyzed meson spectra, but with a lower parameter $T$ extracted for protons. The integrated invariant cross sections calculated from the fitted distributions are found to be consistent with existing measurements and with statistical model predictions.