The PHENIX collaboration presents a systematic study of $\pi^0$ production from $p$ $+$ $p$, $p$ $+$Al, $p$ $+$Au, $d$ $+$Au, and $^{3}$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. Measurements were performed with different centrality selections as well as the total inelastic, 0%--100%, selection for all collision systems. For 0%--100% collisions, the nuclear modification factors, $R_{xA}$, are consistent with unity for $p_T$ above 8 GeV/$c$, but exhibit an enhancement in peripheral collisions and a suppression in central collisions. The enhancement and suppression characteristics are similar for all systems for the same centrality class. It is shown that for high-$p_T$-$\pi^0$ production, the nucleons in the $d$ and $^3$He interact mostly independently with the Au nucleus and that the counter intuitive centrality dependence is likely due to a physical correlation between multiplicity and the presence of a hard scattering process. These observations disfavor models where parton energy loss has a significant contribution to nuclear modifications in small systems. Nuclear modifications at lower $p_T$ resemble the Cronin effect -- an increase followed by a peak in central or inelastic collisions and a plateau in peripheral collisions. The peak height has a characteristic ordering by system size as $p$ $+$Au $>$ $d$ $+$Au $>$ $^{3}$He$+$Au $>$ $p$ $+$Al. For collisions with Au ions, current calculations based on initial state cold nuclear matter effects result in the opposite order, suggesting the presence of other contributions to nuclear modifications, in particular at lower $p_T$.
Differential cross section of $\pi^0$ in p+p collisions at $\sqrt{s}$ = 200 GeV
Invariant yield of $\pi^0$ from (a) p+Al, (b) p+Au, (c) d+Au, and (d) $^{3}$HeAu in different centrality selections at $\sqrt{s}$ = 200 GeV
Nuclear modification factors from inelastic (a) p+Al, (b) p+Au, (c) d+Au, and (d) $^{3}$HeAu collisions at $\sqrt{s}$ = 200 GeV. The right boxes are the $N_{coll}$ uncertainties from the Glauber model, while the left box represents the overall normalization uncertainty from p+p collisions
The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the differential cross section of $\phi$(1020) meson production at forward rapidity in $p$$+$$p$ collisions at $\sqrt{s}=$510 GeV via the dimuon decay channel. The integrated cross section in the rapidity and $p_T$ ranges $1.2<|y|<2.2$ and $2<p_T<7$ GeV/$c$ is $\sigma_\phi=2.79 \pm 0.20\,{\rm (stat)} \pm 0.17\,{\rm (syst)} \pm 0.34\, {\rm (norm)} \times 10^{-2}$~mb. The energy dependence of $\sigma_\phi$ ($1.2<|y|<2.2$; $2<p_T<5$ GeV/$c$) is studied using the PHENIX measurements at $\sqrt{s}=$200 and 510 GeV and the Large-Hadron-Collider measurements at $\sqrt{s}=$2.76 and 7 TeV. The experimental results are compared to various event generator predictions (pythia6, pythia8, phojet, ampt, epos3, and epos-lhc).
The $\phi$-meson-production cross section d$\sigma_{\phi}$/dy in $p$ + $p$ collisions at $\sqrt{s}$ = 510 GeV integrated in the transverse-momentum range 2 < $p_T$ < 7 GeV/$c$.
The $\phi$-meson-production cross section d$\sigma_{\phi}$/dy in $p$ + $p$ collisions at $\sqrt{s}$ = 510 GeV integrated in the transverse-momentum range 2 < $p_T$ < 7 GeV/$c$.
The $\phi$-meson-differential-production cross section d${}^{2}$$\sigma_{\phi}/dp_T dy$ for 1.2 < |y| < 2.2 in $p$ + $p$ collisions at $\sqrt{s}$ = 510 GeV.
The fraction of $J/\psi$ mesons which come from B-meson decay, $\textrm{F}_{B{\rightarrow}J/\psi}$, is measured for J/$\psi$ rapidity \mbox{$1.2<|y|<2.2$} and $p_T>0$ in $p$+$p$ and Cu+Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV with the PHENIX detector. The extracted fraction is $\textrm{F}_{B{\rightarrow}J/\psi}$ = 0.025 $\pm$ 0.006(stat) $\pm$ 0.010(syst) for $p$+$p$ collisions. For Cu+Au collisions, $\textrm{F}_{B{\rightarrow}J/\psi}$ is 0.094 $\pm$ 0.028(stat) $\pm$ 0.037(syst) in the Au-going direction ($-2.2<y<-1.2$) and 0.089 $\pm$ 0.026(stat) $\pm$ 0.040(syst) in the Cu-going direction ($1.2<y<2.2$). The nuclear modification factor, $R_{\rm CuAu}$, of B mesons in Cu+Au collisions is consistent with binary scaling of measured yields in $p$+$p$ at both forward and backward rapidity.
Differential cross section for $p$+$p$ $\rightarrow$ $b\bar{b}$ at $\sqrt{s}$ = 200 GeV.
Fraction F$_{B \rightarrow J/\psi}$ of $B$-meson decays in the inclusive $J/\psi$ sample in $p$+$p$ and Cu+Au collisions at $sqrt{s_{NN}}$ = 200 GeV versus rapidity along with a theoretical estimation based on fixed-order plus next-to-leading logs (FONLL) for the $B \rightarrow J/\psi$ cross section and Color-Evaporation-Model (CEM) for the prompt $J/\psi$.
Fraction F$_{B \rightarrow J/\psi}$ of $B$-meson decays in the inclusive $J/\psi$ sample in $p$+$p$ and Cu+Au collisions at $sqrt{s_{NN}}$ = 200 GeV versus rapidity along with a theoretical estimation based on fixed-order plus next-to-leading logs (FONLL) for the $B \rightarrow J/\psi$ cross section and Color-Evaporation-Model (CEM) for the prompt $J/\psi$.