The first measurements of the invariant differential cross sections of inclusive $\pi^0$ and $\eta$ meson production at mid-rapidity in proton-proton collisions at $\sqrt{s}=0.9$ TeV and $\sqrt{s}=7$ TeV are reported. The $\pi^0$ measurement covers the ranges $0.4<p_T<7$ GeV/$c$ and $0.3<p_T<25$ GeV/$c$ for these two energies, respectively. The production of $\eta$ mesons was measured at $\sqrt{s}=7$ TeV in the range $0.4<p_T<15$ GeV/$c$. Next-to-Leading Order perturbative QCD calculations, which are consistent with the $\pi^0$ spectrum at $\sqrt{s}=0.9$ TeV, overestimate those of $\pi^0$ and $\eta$ mesons at $\sqrt{s}=7$ TeV, but agree with the measured $\eta/\pi^0$ ratio at $\sqrt{s}=7$ TeV.
We present a measurement of the cross section and transverse single-spin asymmetry ($A_N$) for $\eta$ mesons at large pseudorapidity from $\sqrt{s}=200$~GeV $p^{\uparrow}+p$ collisions. The measured cross section for $0.5<p_T<5.0$~GeV/$c$ and $3.0<|\eta|<3.8$ is well described by a next-to-leading-order perturbative-quantum-chromodynamics calculation. The asymmetries $A_N$ have been measured as a function of Feynman-$x$ ($x_F$) from $0.2<|x_{F}|<0.7$, as well as transverse momentum ($p_T$) from $1.0<p_T<4.5$~GeV/$c$. The asymmetry averaged over positive $x_F$ is $\langle{A_{N}}\rangle=0.061{\pm}0.014$. The results are consistent with prior transverse single-spin measurements of forward $\eta$ and $\pi^{0}$ mesons at various energies in overlapping $x_F$ ranges. Comparison of different particle species can help to determine the origin of the large observed asymmetries in $p^{\uparrow}+p$ collisions.