Bayesian upper exclusion limits on the ratio $g_{W'}/g_{W}$ for an SSM-like W$\prime$ boson are shown. The unity coupling ratio (blue dotted curve) corresponds to the SSM common benchmark. The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian lower exclusion limits on the NUGIM G(221) mixing angle $\cot(\theta_{E})$ are shown as a function of the W$\prime$ boson mass. The theoretically excluded region is shaded in grey. The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian upper exclusion limits at 95% CL on the product of the production cross section and branching fraction of a QBH in an associated $ au$ lepton and neutrino final state. Minimum threshold masses $m_{th}$ of up to 6.6 TeV are excluded at 95% CL. The observed limit (solid line) is obtained from the intersection with the LO QBH cross section (blue dotted curve). The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian upper limits at 95% CL on the cross section of the process $pp\rightarrow\tau\nu$ mediated via LQ exchange in the t-channel. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The predicted LQ cross section at LO in the three coupling benchmark scenarios is depicted in different colors for $g_{U}=1$. The uncertainty bandy correspond to the sum in quadrature of PDF and scale variations. The first benchmark scenario considers only couplings to left-handed SM fermions (i.e. $\beta_{\text{R}}^{ij}=0$) and is referred to as "best fit LH". The second benchmark, referred to as "best fit LH+RH", considers $|\beta_{\text{R}}^{\text{b}\tau}|=1$ and all other $\beta_{\text{R}}^{ij}=0$. In the third "democratic" benchmark, equal couplings only to LH fermions are assumed, i.e. $\beta_{\text{L}}^{ij}=1$ and $\beta_{\text{R}}^{ij}=0$ for all $i$ and $j$.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the LH scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the LH+RH scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the democratic scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Bayesian upper exclusion limits at 95% CL on each of the Wilson coefficients described by the EFT model based on 2016-2018 data. The three different coupling types represent a left-handed vector coupling ($\epsilon^{cb}_{L}$), tensor-like coupling ($\epsilon^{cb}_{T}$), and scalar-tensor-like coupling ($\epsilon^{cb}_{S_{L}}$). The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Summary of exclusion limits (expected and observed) calculated at 95% CL for full Run-2 CMS data.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "year_binnumber", following the binning from Figure 4.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "year_binnumber", following the binning used in Figure 4.

Predicted signal yields for the 2017 data-taking period, corresponding to 41.3 fb$^{-1}$, after the application of each search requirement (cumulative) for various signal hypothesis. The requirements listed are presented as total efficiencies w.r.t. the previous selection step.

Centrality dependence of (a) $v_2${2} and (b) $v_2${4}. (a) The red points indicate no pseudorapidity gap whereas the magenta points indicate a pseudorapidity gap of |$\Delta\eta$| > 2.0. (b) The black points indicate $v_2${4} with no pseudorapidity gap, the blue points indicate a two-subevent method with |$\Delta\eta$| > 2.0 but where some short-range pairs are allowed, and the red points indicate a two-subevent method with |$\Delta\eta$| > 2.0 where no short-range pairs are allowed.

Centrality dependence of (a) $v_2${2} and (b) $v_2${4}. (a) The red points indicate no pseudorapidity gap whereas the magenta points indicate a pseudorapidity gap of |$\Delta\eta$| > 2.0. (b) The black points indicate $v_2${4} with no pseudorapidity gap, the blue points indicate a two-subevent method with |$\Delta\eta$| > 2.0 but where some short-range pairs are allowed, and the red points indicate a two-subevent method with |$\Delta\eta$| > 2.0 where no short-range pairs are allowed.

Multi-particle $v_2$ as a function of centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The magenta open diamonds indicate the $v_2${SP }, the blue open squares indicate $v_2${4}, the black open circles indicate $v_2${6}, and the green filled diamonds indicate $v_2${8}.

Multi-particle $v_2$ as a function of centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The magenta open diamonds indicate the $v_2${SP }, the blue open squares indicate $v_2${4}, the black open circles indicate $v_2${6}, and the green filled diamonds indicate $v_2${8}.

Multi-particle $v_2$ as a function of centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The magenta open diamonds indicate the $v_2${SP }, the blue open squares indicate $v_2${4}, the black open circles indicate $v_2${6}, and the green filled diamonds indicate $v_2${8}.

Multi-particle $v_2$ as a function of centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The magenta open diamonds indicate the $v_2${SP }, the blue open squares indicate $v_2${4}, the black open circles indicate $v_2${6}, and the green filled diamonds indicate $v_2${8}.

Cumulant method estimate of $\sigma_{v_{2}}$ / $\langle v_{2} \rangle$ as a function of centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The data are shown as black open squares. The same calculation as done in data is done in ampt, shown as a solid green line. Calculations of $\sigma_{\epsilon_{2}}$ /$ \langle\epsilon_{2}\rangle$ performed in the Monte Carlo Glauber model are shown as blue lines. The solid blue line is the Monte Carlo Glauber calculation done using the same estimate as the data, the dashed blue line is the direct calculation of the moments of the MC Glauber $\epsilon_{2}$ distribution.

Centrality dependence of $v_3${2, |$\Delta\eta$| > 2} in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV, shown as magenta diamonds. The systematic uncertainty is indicated as a shaded magenta band. Also shown as black squares are $\sqrt{\langle v^{2}_{3}\rangle}$ as determined from the folding analysis, which is shown in the next section.

Centrality dependence of $v_3${2, |$\Delta\eta$| > 2} in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV, shown as magenta diamonds. The systematic uncertainty is indicated as a shaded magenta band. Also shown as black squares are $\sqrt{\langle v^{2}_{3}\rangle}$ as determined from the folding analysis, which is shown in the next section.

Centrality dependence of $c_3${4} for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. (a) Calculations using both arms: $c_3${4} (black circles), $c_3${4}$_{ab|ab}$ (blue diamonds), $c_3${4}$_{aa|bb}$ (red squares), and comparison to STAR [36] (black stars). (b) Comparison of $c_3${4} determined using both arms (open symbols) and a single arm (closed symbols). Note that the open black circles are the same in (a) and (b).

Centrality dependence of $c_3${4} for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. (a) Calculations using both arms: $c_3${4} (black circles), $c_3${4}$_{ab|ab}$ (blue diamonds), $c_3${4}$_{aa|bb}$ (red squares), and comparison to STAR [36] (black stars). (b) Comparison of $c_3${4} determined using both arms (open symbols) and a single arm (closed symbols). Note that the open black circles are the same in (a) and (b).

Centrality dependence of $c_3${4} for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. (a) Calculations using both arms: $c_3${4} (black circles), $c_3${4}$_{ab|ab}$ (blue diamonds), $c_3${4}$_{aa|bb}$ (red squares), and comparison to STAR [36] (black stars). (b) Comparison of $c_3${4} determined using both arms (open symbols) and a single arm (closed symbols). Note that the open black circles are the same in (a) and (b).

Folding results for (a) $v_2$, (b) $\sigma_{v_{2}}$ , and (c) $\sigma_{v_{2}}$ / $\langle v_{2}\rangle$. The black lines above and below the points indicate the systematic uncertainties. The red (green) boxes indicate the statistical uncertainties at the 68.27% (95.45%) confidence level. In the case of $\langle v_{2}\rangle$, the statistical uncertainties at the 68.27% confidence level are too small to be seen, and the uncertainties at the 95.45% confidence level are visible but noticeably smaller than the marker size. Shown as blue squares are the same quantities as determined using the cumulant based calculation—these points are slightly offset in the x-coordinate to improve visibility.

Folding results for (a) $v_2$, (b) $\sigma_{v_{2}}$ , and (c) $\sigma_{v_{2}}$ / $\langle v_{2}\rangle$. The black lines above and below the points indicate the systematic uncertainties. The red (green) boxes indicate the statistical uncertainties at the 68.27% (95.45%) confidence level. In the case of $\langle v_{2}\rangle$, the statistical uncertainties at the 68.27% confidence level are too small to be seen, and the uncertainties at the 95.45% confidence level are visible but noticeably smaller than the marker size. Shown as blue squares are the same quantities as determined using the cumulant based calculation—these points are slightly offset in the x-coordinate to improve visibility.

Folding results for (a) $v_2$, (b) $\sigma_{v_{2}}$ , and (c) $\sigma_{v_{2}}$ / $\langle v_{2}\rangle$. The black lines above and below the points indicate the systematic uncertainties. The red (green) boxes indicate the statistical uncertainties at the 68.27% (95.45%) confidence level. In the case of $\langle v_{2}\rangle$, the statistical uncertainties at the 68.27% confidence level are too small to be seen, and the uncertainties at the 95.45% confidence level are visible but noticeably smaller than the marker size. Shown as blue squares are the same quantities as determined using the cumulant based calculation—these points are slightly offset in the x-coordinate to improve visibility.

Folding results for (a) $\langle v_{3}\rangle$, (b) $\sigma_{v_{3}}$ , and (c) $\sigma_{v_{3}}$ /$\langle v_{3}\rangle$. The black lines above and below the points indicate the systematic uncertainties. The red (green) boxes indicate the statistical uncertainties at the 68.27% (95.45%) confidence level. The $\sigma_{v_{3}}$ /$\langle v_{3}\rangle$ values are all ≈ 0.52, the apparent limiting value of this quantity for the Bessel-Gaussian distribution.

Folding results for (a) $\langle v_{3}\rangle$, (b) $\sigma_{v_{3}}$ , and (c) $\sigma_{v_{3}}$ /$\langle v_{3}\rangle$. The black lines above and below the points indicate the systematic uncertainties. The red (green) boxes indicate the statistical uncertainties at the 68.27% (95.45%) confidence level. The $\sigma_{v_{3}}$ /$\langle v_{3}\rangle$ values are all ≈ 0.52, the apparent limiting value of this quantity for the Bessel-Gaussian distribution.

Folding results for (a) $\langle v_{3}\rangle$, (b) $\sigma_{v_{3}}$ , and (c) $\sigma_{v_{3}}$ /$\langle v_{3}\rangle$. The black lines above and below the points indicate the systematic uncertainties. The red (green) boxes indicate the statistical uncertainties at the 68.27% (95.45%) confidence level. The $\sigma_{v_{3}}$ /$\langle v_{3}\rangle$ values are all ≈ 0.52, the apparent limiting value of this quantity for the Bessel-Gaussian distribution.

The second-order Fourier coefficients $v_2 (p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for the centrality classes marked in each panel.

The second-order Fourier coefficients $v_2 (p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for the centrality classes marked in each panel.

The second-order Fourier coefficients $v_2 (p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for the centrality classes marked in each panel.

The second-order Fourier coefficients $v_2 (p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for the centrality classes marked in each panel.

The second-order Fourier coefficients $v_2 (p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for the centrality classes marked in each panel.

The third-order Fourier coefficients $v_3(p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for 0%-30% centrality.

The third-order Fourier coefficients $v_3(p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for 0%-30% centrality.

The third-order Fourier coefficients $v_3(p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for 0%-30% centrality.

The first-order Fourier coefficients $v_1(p_T)$ for charge-combined $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ identified hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for 10%-50% centrality.

The first-order Fourier coefficients $v_1(p_T)$ for charge-combined $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ identified hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for 10%-50% centrality.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Cu + Cu collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Cu + Cu collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled third-order Fourier coefficients $v_3(p_T)/\epsilon_3$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled third-order Fourier coefficients $v_3(p_T)/\epsilon_3$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled third-order Fourier coefficients $v_3(p_T)/\epsilon_3$ for charged hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled third-order Fourier coefficients $\frac{v_3}{($_3 N_{part}^{1/3})}$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled third-order Fourier coefficients $\frac{v_3}{(e_3 N_{part}^{1/3})}$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled third-order Fourier coefficients $\frac{v_3}{(e_3 N_{part}^{1/3})}$ for charged hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

The second-order Fourier coefficients $v_2$($p_T$ ) for charged hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV in comparison with hydrodynamics calculations for the centrality classes marked in each panel.

The third-order Fourier coefficients $v_3(p_T)$ for charged hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV in comparison with hydrodynamics calculations for the centrality classes marked in each panel.

The ratio of $\psi(2S)$ to $J/\psi(1S)$ cross sections, as measured in the dimuon decay channel in pp collisions at 200 GeV over the rapidity interval $1.2<|y|<2.2$, as a function of transverse momentum $p_{T}$.

The double ratios $[\sigma_{\psi(2S)}/\sigma_{\psi(1S)}]_{pA}/[\sigma_{\psi(2S)}/\sigma_{\psi(1S)}]_{pp}$. The global systematic uncertainty is $\pm 14\%$.

The charged-particle multiplicity $dN_{ch}/d\eta$ determined from AMPT simulations and the nuclear overlap $S_T$ from Glauber simulations for the data shown in Fig. 5. Here we assume $dN/d\eta$ = $\frac{3}{2}dN_{ch}/d\eta$.

The charged-particle multiplicity $dN_{ch}/d\eta$ determined from AMPT simulations and the nuclear overlap $S_T$ from Glauber simulations for the data shown in Fig. 5. Here we assume $dN/d\eta$ = $\frac{3}{2}dN_{ch}/d\eta$.

The charged-particle multiplicity $dN_{ch}/d\eta$ determined from AMPT simulations and the nuclear overlap $S_T$ from Glauber simulations for the data shown in Fig. 5. Here we assume $dN/d\eta$ = $\frac{3}{2}dN_{ch}/d\eta$.

Data table for $\sigma_{\psi(2S)}/\sigma_{\psi(1S)}$ ratios, measured in the dimuon decay channel. Here "fwd" refers to the forward rapidity interval $1.2<y<2.2$, "bwd" refers to the backward rapidity interval $-2.2<y<-1.2$, and "comb" refers to the rapidity interval $1.2<|y|<2.2$.

Efficiency corrected cumulants of net-charge distributions as a function of $\langle N_{part} \rangle$ from Au+Au collisions at different collision energies.

$\langle N_{part} \rangle$ dependence of efficiency corrected $\mu / \sigma^2$ of net-charge distributions for Au+Au collisions at different collision energies.

$\langle N_{part} \rangle$ dependence of efficiency corrected $S \sigma$ of net-charge distributions for Au+Au collisions at different collision energies.

$\langle N_{part} \rangle$ dependence of efficiency corrected $\kappa \sigma^2$ of net-charge distributions for Au+Au collisions at different collision energies.

$\langle N_{part} \rangle$ dependence of efficiency corrected $S \sigma^3 / \mu$ of net-charge distributions for Au+Au collisions at different collision energies.

The energy dependence of efficiency corrected $\mu / \sigma^2$, $S \sigma$, $\kappa \sigma^2$, and $S \sigma^3 / \mu$ of netcharge distributions for central (0%–5%) Au+Au collisions.

The energy dependence of the chemical freeze-out parameter $\mu_B$.

Away-side jet widths from a Gaussian fit by $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta in Au+Au collisions.

Away-side $I_{AA}$ for the entire away-side $|\Delta \phi - \pi| < \pi /2$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Away-side $I_{AA}$ for the entire away-side $|\Delta \phi - \pi| < \pi /2$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Away-side $I_{AA}$ for a narrow “head” $|\Delta \phi - \pi| < \pi /6$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Away-side $I_{AA}$ for a narrow “head” $|\Delta \phi - \pi| < \pi /6$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Supplemental near-side jet widths from a Gaussian fit by $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta in $p+p$ collisions.

Supplemental near-side jet widths from a Gaussian fit by $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta in Au+Au collisions.

Supplemental near-side jet widths from a Gaussian fit by $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta in Au+Au collisions.

Supplemental away-side $I_{AA}$ for $|\Delta \phi - \pi| < \pi /3$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Supplemental away-side $I_{AA}$ for $|\Delta \phi - \pi| < \pi /3$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

$p+p$ yield vs $p_T^a$ supplemental tables.

Au+Au yield vs $p_T^a$ supplemental tables.

Au+Au yield vs $p_T^a$ supplemental tables.

\Delta v_1 between proton and anti-proton vs rapidity in Au+Au 200 GeV 50-80% centrality.

\Delta v_1 between proton and anti-proton vs rapidity in Zr+Zr and Ru+Ru (combined) 200 GeV 50-80% centrality.

\Delta v_1 between proton and anti-proton vs rapidity in Au+Au 27 GeV 50-80% centrality.

\Delta dv1/dy as a function of centrality in Au+Au 200 GeV.

\Delta dv1/dy as a function of centrality in Ru+Ru and Zr+Zr 200 GeV

\Delta dv1/dy as a function of centrality in Au+Au 27 GeV.

Isoscalar WIMP-nucleon elastic coupling limit for Operator 9

Isoscalar WIMP-nucleon elastic coupling limit for Operator 12

Isoscalar WIMP-nucleon elastic coupling limit for Operator 7

Isoscalar WIMP-nucleon elastic coupling limit for Operator 13

Isoscalar WIMP-nucleon elastic coupling limit for Operator 3

Isoscalar WIMP-nucleon elastic coupling limit for Operator 5

Isoscalar WIMP-nucleon elastic coupling limit for Operator 11

Isoscalar WIMP-nucleon elastic coupling limit for Operator 14

Isoscalar WIMP-nucleon elastic coupling limit for Operator 4

Isoscalar WIMP-nucleon elastic coupling limit for Operator 10

Isoscalar WIMP-nucleon elastic coupling limit for Operator 1

Isoscalar WIMP-nucleon elastic coupling limit for Operator 6

Isoscalar WIMP-nucleon elastic coupling limit for Operator 15

Isovector WIMP-nucleon elastic coupling limit for Operator 3

Isovector WIMP-nucleon elastic coupling limit for Operator 10

Isovector WIMP-nucleon elastic coupling limit for Operator 13

Isovector WIMP-nucleon elastic coupling limit for Operator 12

Isovector WIMP-nucleon elastic coupling limit for Operator 8

Isovector WIMP-nucleon elastic coupling limit for Operator 1

Isovector WIMP-nucleon elastic coupling limit for Operator 4

Isovector WIMP-nucleon elastic coupling limit for Operator 5

Isovector WIMP-nucleon elastic coupling limit for Operator 7

Isovector WIMP-nucleon elastic coupling limit for Operator 11

Isovector WIMP-nucleon elastic coupling limit for Operator 15

Isovector WIMP-nucleon elastic coupling limit for Operator 9

Isovector WIMP-nucleon elastic coupling limit for Operator 14

Isovector WIMP-nucleon elastic coupling limit for Operator 6

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 12

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 3

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 9

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 7

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 10

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 15

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 8

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 4

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 5

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 6

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 13

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 11

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 14

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 1

Isovector WIMP-nucleon inelastic coupling limit for Operator 13

Isovector WIMP-nucleon inelastic coupling limit for Operator 1

Isovector WIMP-nucleon inelastic coupling limit for Operator 4

Isovector WIMP-nucleon inelastic coupling limit for Operator 7

Isovector WIMP-nucleon inelastic coupling limit for Operator 9

Isovector WIMP-nucleon inelastic coupling limit for Operator 3

Isovector WIMP-nucleon inelastic coupling limit for Operator 11

Isovector WIMP-nucleon inelastic coupling limit for Operator 5

Isovector WIMP-nucleon inelastic coupling limit for Operator 10

Isovector WIMP-nucleon inelastic coupling limit for Operator 8

Isovector WIMP-nucleon inelastic coupling limit for Operator 14

Isovector WIMP-nucleon inelastic coupling limit for Operator 15

Isovector WIMP-nucleon inelastic coupling limit for Operator 6

Isovector WIMP-nucleon inelastic coupling limit for Operator 12

Forward-rapidity $J/\psi$ —> $\mu^+\mu^-$ $d$+Au Nuclear Dependence at $\sqrt{s}$ = 200 GeV. (sys. A, B systematics are relative, i.e. they multiply the $R_{dAu}$ value. The sysA uncertainty includes statistical uncertainties as well as point-to-point uncorrelated systematics, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

Mid rapidity $d$+Au —> $e^+e^-$ $J/\psi$ $R_{dAu}$ at $\sqrt{s}$=200 GeV. (Sys. A, B systematics are absolute, i.e. they add/subtract directly from $R_{dAu}$. The sys. A uncertainty includes statistical uncertainties as well as point-to-point uncorrelated systematics, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

$d$+Au vs centrality vs $y$. (All uncertainties are absolute. The sys. A uncertainty includes both the statistical uncertainty and the point-to-point uncorrelated systematic, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

Forward-rapidity $J/\psi$ —> $\mu^+\mu^-$ $d$+Au Nuclear Dependence at $\sqrt{s}$ = 200 GeV. (sys. A, B systematics are relative, i.e. they multiply the $R_{dAu}$ value. The sysA uncertainty includes statistical uncertainties as well as point-to-point uncorrelated systematics, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

Mid rapidity $d$+Au —> $e^+e^-$ $J/\psi$ $R_{dAu}$ at $\sqrt{s}$=200 GeV. (Sys. A, B systematics are absolute, i.e. they add/subtract directly from $R_{dAu}$. The sys. A uncertainty includes statistical uncertainties as well as point-to-point uncorrelated systematics, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

Forward-rapidity $J/\psi$ —> $\mu^+\mu^-$ $d$+Au Nuclear Dependence at $\sqrt{s}$ = 200 GeV. (sys. A, B systematics are relative, i.e. they multiply the $R_{CP}$ value. The sys. A uncertainty includes statistical uncertainties as well as point-to-point uncorrelated systematics, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

Mid rapidity $d$+Au —> $e^+e^-$ $J/\psi$ $R_{CP}$ at $\sqrt{s}$=200 GeV. (Sys. A, B systematics are absolute, i.e. they add/subtract directly from $R_{CP}$. The sys. A uncertainty includes statistical uncertainties as well as point-to-point uncorrelated systematics, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

Comparison between data and background prediction before the fit (Pre-Fit) for the mass of the FCNC top-quark candidate in SR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The four FCNC LH signals are also shown separately, normalized to five times the cross-section corresponding to the most stringent observed branching ratio limits. The first (last) bin in all distributions includes the underflow (overflow). The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction before the fit (Pre-Fit) for the mass of the SM top-quark candidate in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The four FCNC LH signals are also shown separately, normalized to five times the cross-section corresponding to the most stringent observed branching ratio limits. The first (last) bin in all distributions includes the underflow (overflow). The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction before the fit (Pre-Fit) for the transverse momentum of the Z boson candidate in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The four FCNC LH signals are also shown separately, normalized to five times the cross-section corresponding to the most stringent observed branching ratio limits. The first (last) bin in all distributions includes the underflow (overflow). The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{1}$ discriminant in the mass sideband CR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{2}^{u}$ discriminant in the mass sideband CR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{1}$ discriminant in SR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{2}^{u}$ discriminant in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{1}$ discriminant in the mass sideband CR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{2}^{c}$ discriminant in the mass sideband CR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{1}$ discriminant in SR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{2}^{c}$ discriminant in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the leading lepton $p_{T}$ in the $t\bar{t}$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The last bin includes the overflow. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the third lepton $p_{T}$ in the $t\bar{t}$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The last bin includes the overflow. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the leading lepton $p_{T}$ in the $t\bar{t}Z$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The last bin includes the overflow. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{1}$ discriminant in the $t\bar{t}Z$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

The observed 90% C.L upper limit on neutrino effective millicharge (\delta_{Q}[e_{0}]) in SR1. The data contains observed upper limit, median sensitivity and 1\sigma and 2\sigma sensitivity range.

The observerd 90% C.L upper limit on solar axion g_{ae} coupling constant in valid across the mass range for SR1. The data contains observed upper limit, median sensitivity and 1\sigma and 2\sigma sensitivity range as a function of Solar Axion Mass .

The observed 90% C.L. upper limit on axion like particle (ALP) coupling constant (g_{ae}) in SR1. ALP mass spanning from 1 keV to 17 keV is considered. The data contains observed upper limit, median sensitivity and 1\sigma and 2\sigma sensitivity range as a function of ALP Mass.

The observed 90% C.L. upper limit on hidden photon (HP) coupling constant squared (\kappa^2) in SR1. HP mass spanning from 1 keV to 17 keV is considered. The data contains observed upper limit, median sensitivity and 1\sigma and 2\sigma sensitivity range as a function of HP Mass.

The observed 90% C.L. upper limit for the spin-independent (cm^2) WIMP cross-section vs. WIMP mass (GeV/c^2), using the Migdal signal pathway. The result is obtained for WIMPs undergoing the Migdal effect from 0.5 to 9 GeV/c2. The data contains observed upper limit, median sensitivity and 1\sigma and 2\sigma sensitivity range as a function of WIMP Mass.

The observed 90% C.L. upper limit for the spin-dependent WIMP cross-section (cm^2) vs. WIMP mass (GeV/c^2) for coupling on neutrons. The result is obtained for WIMPs undergoing the Migdal effect from 0.5 to 9 GeV/c2. The data contains observed upper limit, median sensitivity and 1\sigma and 2\sigma sensitivity range as a function of WIMP Mass.

The observed 90% C.L. upper limit for the spin-dependent WIMP cross-section (cm^2) vs. WIMP mass (GeV/c^2) for coupling on protons. The result is obtained for WIMPs undergoing the Migdal effect from 0.5 to 9 GeV/c2. The data contains observed upper limit, median sensitivity and 1\sigma and 2\sigma sensitivity range as a function of WIMP Mass.

Data selection efficiency as a function of nuclear recoil energy

Data points used in analysis in log_10(S2)-S1 space

High $p_T$ ($p_T$ > 6 GeV/$c$) integrated $v_2$ vs $N_{part}$ for (a) $\pi^0$ (b) inclusive photon, and (c) direct photon. Results are shown with both reaction plane detectors$:$ (solid black circles) BBC and (solid red squares) RXN. Each point represents a 10% wide centrality bin from 60%-0%.

Inclusive unlike-sign $\mu \mu$ pair yield from $p$+$p$ collisions at $\sqrt{s}$ = 200 GeV as a function of mass in different $p_T$ slices for the south and north muon arms and the ratio of data to expected sources.

Inclusive unlike-sign $\mu \mu$ pair yield from $p$+$p$ collisions at $\sqrt{s}$ = 200 GeV as a function of mass in different $p_T$ slices for the south and north muon arms and the ratio of data to expected sources.

Inclusive like-sign $\mu \mu$ pair yield from $p$+$p$ collisions at $\sqrt{s}$ = 200 GeV as a function of mass in different $p_T$ slices for the south and north muon arms and the ratio of data to expected sources.

Inclusive like-sign $\mu \mu$ pair yield from $p$+$p$ collisions at $\sqrt{s}$ = 200 GeV as a function of mass in different $p_T$ slices for the south and north muon arms and the ratio of data to expected sources.

Inclusive like-sign $\mu \mu$ pair yield from $p$+$p$ collisions at $\sqrt{s}$ = 200 GeV as a function of mass in different $p_T$ slices for the south and north muon arms and the ratio of data to expected sources.

The $\mu \mu$ pair data in the charm mass region as a function of $\Delta \phi$.

The $\mu \mu$ pair data in the charm mass region as a function of $\Delta \phi.$

The $\mu \mu$ pair data in the charm mass region as a function of pair $p_T$.

The $\mu \mu$ pair data in the charm mass region as a function of pair $p_T$

The like-sign $\mu \mu$ pair data in the bottom mass region as a function of $\Delta \phi$.

The like-sign $\mu \mu$ pair data in the bottom mass region as a function of $\Delta \phi$.

The like-sign $\mu \mu$ pair data in the bottom mass region as a function of pair $p_T$.

The like-sign $\mu \mu$ pair data in the bottom mass region as a function of pair $p_T$.

The unlike-sign $\mu \mu$ pair data used to determine the Drell-Yan contribution as a function of mass.

The unlike-sign $\mu \mu$ pair data used to determine the Drell-Yan contribution as a function of mass.

The unlike-sign $\mu \mu$ pair data used to determine the Drell-Yan contribution as a function of pair $p_T$.

The unlike-sign $\mu \mu$ pair data used to determine the Drell-Yan contribution as a function of pair $p_T$.

Summary of arm-averaged relative systematic uncertainties for the total bottom cross section $\sigma _{b\bar{b}}$, the differential Drell-Yan cross section $d^{2}\sigma _{DY\implies \mu\mu}/dmdy$, and the $b\bar{b}$ ($c\bar{c}$) differential yields $dN_{b\bar{b}(c\bar{c})\implies \mu\mu}/d\Delta\phi$. The systematic uncertainty type is indicated in the second column and is applicable only to the differential measurements. The uncertainties for the differential measurements vary with azimuthal opening angle, pair $p_{T}$ , or mass. Asymmetric uncertainties are quoted in bracketed values. For the $c\bar{c}$ measurement, the regions $\Delta\phi$ < $\pi$/2, $p_{T}$ < 0.5 GeV/$c$ and $p_{T}$ > 2.0 GeV/$c$ are excluded because the yield approaches zero and relative systematic uncertainties diverge. With these regions excluded, the difference between the systematic uncertainties of all measurements for the south and north muon arms differs by no greater than 2% for all systematic uncertainties sources.

The differential yield $dN/d\phi$ of unlike-sign muon pairs from charm with mass 1.5-2.5 GeV/$c^2$ in the ideal muon arm acceptance, as a function of the pair azimuthal opening angle. Statistical (stat.) and types B and C systematic errors (sys. B & sys. C) are all $10^{-9}\times$[rad$^{-1}$].

The differential yield $dN/d\phi$ of like-sign muon pairs from bottom with mass 3.5-10.0 GeV/$c^2$ in the ideal muon arm acceptance, as a function of the pair azimuthal opening angle. Statistical (stat.) and types B and C systematic errors (sys. B & sys. C) are all $10^{-9}\times$[rad$^{-1}$].

The differential yield $dN/dp_{T}$ of unlike-sign muon pairs from bottom with mass 1.5-2.5 GeV/$c^2$ in the ideal muon arm acceptance, as a function of the pair transverse momentum. Statistical (stat.) and types B and C systematic errors (sys. B & sys. C) are all $10^{-9}\times$[$c$/GeV].

The differential yield $dN/dp_{T}$ of like-sign muon pairs from bottom with mass 3.5-10.0 GeV/$c^2$ in the ideal muon arm acceptance, as a function of the pair transverse momentum. Statistical (stat.) and types B and C systematic errors (sys. B & sys. C) are all $10^{-9}\times$[$c$/GeV].

$\sigma_{b\bar{b}}$ from fit using different models. Only statistical uncertainties are shown.

The differential Drell-Yan cross section $\frac{d^{2}\sigma}{dmdy}$ as a function of the muon pair mass, where the muon pair rapidity |$y\mu\mu$| is between 1.2 and 2.2. Statistical (stat.) and types B and C systematic errors (sys. B & sys. C) are all $10^{-9}\times$[pb $c^2$/GeV]. In the data table below, the value of mass correspond to the center of the stated range.

The differential Drell-Yan cross section $\frac{1}{2\pi p_{T}} \frac{d^{2}\sigma}{dydp_{T}}$ as a function of the muon pair transverse momentum, where the muon pair mass $m_{\mu\mu}$ is between 4.8 and 8.2 GeV/$c^2$ and the muon pair rapidity |$y\mu\mu$| is between 1.2 and 2.2. Statistical (stat.) and types B and C systematic errors (sys. B & sys. C) are all $10^{-9}\times$[pb $c^2$/GeV]. In the data table below, the value of p_T correspond to the center of the stated range.

The relative nuclear-modification factor for U+U and Au+Au as a function of collision centrality.

The ratio of the invariant yield for U+U to that for Au+Au, as a function of collision centrality.

Probability of finding at least one TAR jet, where the p_T-leading TAR jet passes the m_Wcand and D₂^&beta;=1 requirements, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=2100 GeV, with the preselections applied.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=500 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=1000 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=1700 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=2100 GeV, with the preselections applied that do not pass the requirements of the merged category.

Observed exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected+1&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected-1&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected+2&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected-2&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=0.5 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=0.5 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=1 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=1 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=1.7 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=1.7 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=2.1 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=2.1 TeV, shown in dashed blue.

Data overlaid on SM background yields stacked in each SR and CR category after the fit to data ('Post-fit'). The yields in the SR are broken down into their contributions to the individual bins. The maximum-likelihood estimators are set to the conditional values of the CR-only fit, and propagated to SR and CRs.

Dominant sources of uncertainty for three dark Higgs scenarios after the fit to data. The uncertainties are quantified in terms of their contribution to the fitted signal uncertainty that is expressed relative to the theory prediction. Three representative dark Higgs signal scenarios with g_q=0.25, g_&chi;=1.0, sin&theta;=0.01 and m_&chi;=200 GeV are considered, which are indicated using the (m_Z', m_s) format in units of GeV in the table columns.

Cumulative efficiencies in the merged category for three representative dark Higgs signal scenarios with g_q=0.25, g_&chi;=1.0, sin&theta;=0.01, m_Z' = 1 TeV, and m_&chi;=200 GeV considering s&rarr;W(&ell;&nu;)W(qq) decays only.

Cumulative efficiencies in the resolved category for three representative dark Higgs signal scenarios with g_q=0.25, g_&chi;=1.0, sin&theta;=0.01, m_Z' = 1 TeV, and m_&chi;=200 GeV considering s&rarr;W(&ell;&nu;)W(qq) decays only.

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={300, 350, 400, 500, 750} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={1000, 1700} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={2100, 2500, 2900, 3300} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Figure 4a

Figure 4b

Figure 4c

Figure 5a

Figure 5b

Figure 5c

Figure 7

Figure 8a

Figure 8b

Figure 8c

Figure 9a

Figure 9b

Figure 9c

Figure 10a

Figure 10b

Figure 10c

$P_{\rm H}$ measurements are shown as a function of hyperon $p_{\rm T}$ at $\sqrt{s_{\rm NN}}$=19.6 and 27 GeV. Statistical uncertainties are represented as lines while systematic uncertainties are represented as boxes. There is no observed dependence of $P_{\rm H}$ on $p_{\rm T}$ at $\sqrt{s_{\rm NN}}$=19.6 or 27 GeV, consistent with previous observations.

$P_{\rm H}$ measurements are shown as a function of hyperon $y$ at $\sqrt{s_{\rm NN}}$=19.6 and 27 GeV. Statistical uncertainties are represented as lines while systematic uncertainties are represented as boxes. The data set at $\sqrt{s_{\rm NN}}$=19.6 GeV takes advantage of STAR upgrades to reach larger $|y|$.

HFE (electrons from semileptonic decays of heavy-flavor hadrons) invariant yields in different centrality intervals of Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV. The vertical bars and the boxes represent statistical and systematic uncertainties, respectively. The horizontal bars indicate the bin width.

HFE (electrons from semileptonic decays of heavy-flavor hadrons) $R_{\rm AA}$ (red circles) as a function of $p_{\rm T}$ in different centrality intervals of Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV, compared with STAR (yellow stars) and PHENIX (green squares) published results, and Duke ((modified Langevin transport model, blue line) and PHSD (parton-hadron-string dynamics model, orange line) model calculations. Vertical bars and boxes around data points represent combined statistical and systematic uncertainties from both Au+Au and $p$+$p$ measurements, respectively. Boxes at unity show the global uncertainties, which for this analysis include the 8% global uncertainty on $p$+$p$ reference and the $N_{\rm coll}$ uncertainties. The left box is for PHENIX and the right one for STAR.

HFE (electrons from semileptonic decays of heavy-flavor hadrons) $R_{\rm AA}$ (red circles) as a function of $N_{\rm part}$ in Au+Au collisions at $\sqrt{s_{\rm NN}}=200$ GeV, compared with PHENIX measurements (green squares), and Duke (blue line) and PHSD (orange line) model calculations. Vertical bars and boxes around data points represent statistical and systematic uncertainties from Au+Au measurements, respectively. The gray band represents the $N_{\rm coll}$ uncertainties. The boxes at unity show the global uncertainties including the total uncertainties of the $p$+$p$ reference.

Lévy scale parameter $R$ versus average $m_T$ of the pair. The graphical representation of statistical and systematic uncertainties is the same as in Fig. 4.

Lévy index parameter $\alpha$ versus average $m_T$ of the pair. Statistical and systematic uncertainties are indicated similarly to Fig. 4. The horizontal line, %\aplha$ = 1.207, represents the 0%-30% centrality average value of %\aplha$.

Inverse square of the Lévy scale parameter 1/$R^2$ versus average $m_T$ of the pair. Statistical and systematic uncertainties shown as bars and boxes, respectively.

New scale parameter $\hat{R}$ versus average $m_T$ of the pair, with a linear fit. Statistical and systematic uncertainties shown as bars and boxes, respectively.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Dependence of $\overline{p}$-$\overline{d}$ correlation on pseudorapidity acceptance of $\overline{p}$ and $\overline{d}$ selection in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. Results are for 20.0--30.0$\%$ collision centrality.

Dependence of $\overline{p}$-$\overline{d}$ correlation on pseudorapidity acceptance of $\overline{p}$ and $\overline{d}$ selection in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. Results are for 50.0--60.0$\%$ collision centrality.

HBT parameters of charge-combined kaon pairs, shown as value $\pm$ statistical uncertainty [absolute value] $\pm$ systematic uncertainty [%] for the centrality bins shown in Fig. 3.

$R_o/R_s$ of positive pion pairs shown as value $\pm$ statistical uncertainty [absolute value] $\pm$ systematic uncertainty [%] for the centrality bins shown in Fig. 4.

$R_o/R_s$ of negative pion pairs shown as value $\pm$ statistical uncertainty [absolute value] $\pm$ systematic uncertainty [%] for the centrality bins shown in Fig. 4.

$R_o/R_s$ of charge-combined kaon pairs shown as value $\pm$ statistical uncertainty [absolute value] $\pm$ systematic uncertainty [%] for the centrality bins shown in Fig. 4.

Azimuthal angle dependence of HBT radii of charged kaons, shown as value $\pm$ statistical uncertainty $\pm$ systematic uncertainty for the centrality bins shown in Fig. 8.

Azimuthal angle dependence of HBT radii of charged pions, shown as value $\pm$ statistical uncertainty $\pm$ systematic uncertainty for the 0%–20% and 20%–60% centrality bins.

Oscillation amplitudes relative to the event plane for charged pions, shown as value $\pm$ statistical uncertainty $\pm$ systematic uncertainty for the 0%–20% and 20%–60% centrality bins.

Oscillation amplitudes relative to the event plane for charged kaons, shown as value $\pm$ statistical uncertainty $\pm$ systematic uncertainty for the 0%–20% and 20%–60% centrality bins.

Summary of extracted parameters in the blast-wave model fit for two fitting conditions. The $T_f$ and $\rho_0$ parameters were obtained by fits to $p_T$ spectra, and the other parameters were obtained by a simultaneous fit to $v_2$ and the HBT radii.

$J/\psi$ $R_{CP}$ for 0-20%, 20-40%, and 40-60% (central) relative to 60-86% (peripheral) Au+Au collisions at 62.4 GeV. Type A errors are point-to-point uncorrelated, Type B errors are point-to-point correlated, and Type C are global systematic from the uncertainty in the peripheral $<N_{coll}>$ value.

$J/\psi$ $R_{CP}$ for 0-40% (central) relative to 40-86% (peripheral) Au+Au collisions at 39 GeV. Type A errors are point-to-point uncorrelated, Type B errors are point-to-point correlated, and Type C are global systematic from the uncertainty in the peripheral $<N_{coll}>$ value.

$J/\psi$ $R_{AA}$ at $\sqrt{s_{NN}}$ = 39 and 62.4 GeV. Type A errors are point-to-point uncorrelated, Type B errors are point-to-point correlated, and Type C are global systematic from the uncertainty in the peripheral $<N_{coll}>$ value.

$J/\psi$ $R_{AA}$ at $\sqrt{s_{NN}}$ = 39 and 62.4 GeV.

Comparison of neutral pion $A_N$ as function of $x_F$ from $\sqrt{s}$ = 62.4 GeV.

Fractional composition of electromagnetic clusters in the MPC at $\sqrt{s}$ = 200 GeV for $x_F$ > 0.4.

The $A_N$ of electromagnetic clusters at $\sqrt{s}$ = 200 GeV as a function of $x_F$ for $3.1 < |\eta| < 3.5$.

The $A_N$ of electromagnetic clusters at $\sqrt{s}$ = 200 GeV as a function of $x_F$ for $3.5 < |\eta| < 3.8$.

The $A_N$ of electromagnetic clusters at $\sqrt{s}$ = 200 GeV at large |$x_F$| > 0.4 in forward/backward rapidities as function of $p_T$.

Transverse single spin asymmetries as function of $p_T$ and $x_F$ at forward/backward rapidities.

Transverse single spin asymmetries $A_N$ of $\pi^0$ mesons at $\sqrt{s}$ = 200 GeV at midrapidity as function of $p_T$.

Transverse single spin asymmetries $A_N$ of $\pi^0$ mesons at $\sqrt{s}$ = 200 GeV at midrapidity as function of $p_T$.

Transverse single spin asymmetries $A_N$ of $\eta$ mesons at $\sqrt{s}$ = 200 GeV at midrapidity as function of $p_T$.

Transverse single spin asymmetries $A_N$ of $\eta$ mesons at $\sqrt{s}$ = 200 GeV at midrapidity as function of $p_T$.

Invariant transverse momentum spectra of $\omega$ production in $p$+$p$ and $d$+Au collisions at $\sqrt{s}$=200 GeV.

Invariant transverse momentum spectra of $\omega$ production in $p$+$p$ and $d$+Au collisions at $\sqrt{s}$=200 GeV.

Invariant transverse momentum spectra of $\omega$ production in Cu+Cu (left) and Au+Au (right) collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel for three centrality bins and minimum bias.

Invariant transverse momentum spectra of $\omega$ production in Cu+Cu (left) and Au+Au (right) collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel for three centrality bins and minimum bias.

Invariant transverse momentum spectra of $\omega$ production in Cu+Cu (left) and Au+Au (right) collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel for three centrality bins and minimum bias.

Invariant transverse momentum spectra of $\omega$ production in Cu+Cu (left) and Au+Au (right) collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel for three centrality bins and minimum bias.

Invariant transverse momentum spectra of $\omega$ production in Cu+Cu (left) and Au+Au (right) collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel for three centrality bins and minimum bias.

Invariant transverse momentum spectra of $\omega$ production in Cu+Cu (left) and Au+Au (right) collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel for three centrality bins and minimum bias.

Measured ${\omega}/{\pi}$ ratio as a function of $p_T$ in $p$ + $p$ collisions at $\sqrt{s}$=200 GeV.

Measured ${\omega}/{\pi}$ ratio as a function of $p_T$ in $p$ + $p$ collisions at $\sqrt{s}$=200 GeV.

Measured ${\omega}/{\pi}$ ratio as a function of $p_T$ in $p$ + $p$ collisions at $\sqrt{s}$=200 GeV.

${\omega}/{\pi}$ ratio versus transverse momentum in $d$+Au (0–88%) for $\omega \rightarrow {e}^{+}{e}^{-}$.

${\omega}/{\pi}$ ratio versus transverse momentum in $d$+Au (0–88%) for $\omega \rightarrow {\pi}^{+}{\pi}^{-}\pi^0$.

${\omega}/{\pi}$ ratio versus transverse momentum in $d$+Au (0–88%) for $\omega \rightarrow {\pi}^{0}\gamma$.

${\omega}/{\pi}$ ratio versus transverse momentum in Cu+Cu (0–94%) for $\omega \rightarrow {\pi}^{0}\gamma$.

${\omega}/{\pi}$ ratio versus transverse momentum in Au+Au (0–92%) for $\omega \rightarrow {\pi}^{0}\gamma$.

${\omega}/{\pi}$ ratio versus transverse momentum in Au+Au (0–92%) for $\omega \rightarrow {\pi}^{0}\gamma$.

Nuclear modification factor, RdAu, measured for the ${R}_{dAu}$ in 0–88, 0–20, 20-40, and 40-60% centrality bins in $d$+Au collisions at $\sqrt{s}$=200 GeV.

Nuclear modification factor, RdAu, measured for the ${R}_{dAu}$ in 60–88% centrality bins in $d$+Au collisions at $\sqrt{s}$=200 GeV.

Nuclear modification factor, RdAu, measured for the ${R}_{dAu}$ in 0–88, 20-40, and 40–60% centrality bins in $d$+Au collisions at $\sqrt{s}$=200 GeV.

Nuclear modification factor, RdAu, measured for the ${R}_{dAu}$ in 0–20 and 60–88% centrality bins in $d$+Au collisions at $\sqrt{s}$=200 GeV.

Nuclear modification factor, RdAu, measured for the ${R}_{dAu}$ in 0–88, 0–20, 20-40, 40-60, and 60–88% centrality bins in $d$+Au collisions at $\sqrt{s}$=200 GeV.

$R_{AA}$ of the $\omega$ in Cu+Cu collisions from the $\omega \rightarrow \pi^0\gamma$ decay channel.

$R_{AA}$ of the $\omega$ in Au+Au collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel.

$R_{AA}$ of the $\omega$ in Au+Au collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel.

$R_{AA}$ of the $\omega$ in Au+Au collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel.

$R_{AA}$ of the $\omega$ in Au+Au collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel.

$R_{AA}$ of the $\omega$ in and Au+Au collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel.

($R_{AA}$) of the $\omega \rightarrow \pi^+\pi^-\pi^0$ collision for the $\omega$ meson integrated over the range $p_T$ > 7 $\mathrm{GeV}$/$c$ as a function of the number participating nucleons ($N_{part}$).

($R_{AA}$) of the $\omega \rightarrow \pi^0 \gamma$ collision for the $\omega$ meson integrated over the range $p_T$ > 7 $\mathrm{GeV}$/$c$ as a function of the number participating nucleons ($N_{part}$).

$R_{AA}$ of the $\omega \rightarrow \pi^0 \gamma$ collision for the $\omega$ meson integrated over the range $p_T$ > 7 $\mathrm{GeV}$/$c$ as a function of the number participating nucleons ($N_{part}$).

Parameters of the Tsallis fit with parameter $n$ constrained to a fixed linear dependence on mass (for mesons). Cross sections are in $\mu$b for $J/\psi$ and $\psi^{\prime}$ and in mb for all other particles.

Constant and linear fits to the Tsallis parameter $T$ of mesons with fixed parameter $n$. The last column (Prob.) gives the probability estimated by the $\chi^2$/$n.d.f.$ of the fit.

Parameters of the Tsallis fit with Eq. 8 in the paper with parameters $n$ and $T$ constrained to have a fixed linear dependence on mass (for mesons). Cross sections are in $\mu$b for $J/\psi$ and $\psi^{\prime}$ and in mb for all other particles.

Cross sections in mb and $<p_T>$ in GeV/$c$ of different particles in $p$+$p$ collisions at $\sqrt{s}$ = 200 GeV. PHENIX columns show the values obtained by fits to the experiment spectra with the Tsallis functional form as described in the text. One should state explicitly that these values do not supersede values given in the ”Published” column by the experiments, in their publications listed in the last column, or elsewhere. An additional 9.7% systematic uncertainty should be added to all $d\sigma$/$dy$ values listed in the column “PHENIX” to account for the trigger uncertainties. Values in the column “Published” are also given without these systematic uncertainties. The column “S.M.” is the prediction of the statistical model discussed in the text. The characteristic widths of the particle spectra are $\langle p_T \rangle$ for all species except for $J/\psi$ and $\psi^{\prime}$ for which the values given in the Table are $\langle p_T^2 \rangle$. For $\psi^{\prime}$ the integration is done in the $p_T$ regionbelow 5 GeV/$c$. All errors are the combined statistical and systematic uncertainties.

Cross sections in mb and $<p_T>$ in GeV/$c$ of different particles in $p$+$p$ collisions at $\sqrt{s}$ = 200 GeV. PHENIX columns show the values obtained by fits to the experiment spectra with the Tsallis functional form as described in the text. One should state explicitly that these values do not supersede values given in the ”Published” column by the experiments, in their publications listed in the last column, or elsewhere. An additional 9.7% systematic uncertainty should be added to all $d\sigma$/$dy$ values listed in the column “PHENIX” to account for the trigger uncertainties. Values in the column “Published” are also given without these systematic uncertainties. The column “S.M.” is the prediction of the statistical model discussed in the text. The characteristic widths of the particle spectra are $\langle p_T \rangle$ for all species except for $J/\psi$ and $\psi^{\prime}$ for which the values given in the Table are $\langle p_T^2 \rangle$. For $\psi^{\prime}$ the integration is done in the $p_T$ regionbelow 5 GeV/$c$. All errors are the combined statistical and systematic uncertainties.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

The invariant differential cross section ${d^2 \sigma}$/${2\pi p_T dy p_T}$ of $\omega$ meson production measured in the indicated decay channel. Type A uncertainties are $p_T$ independent. Type C & B systematic uncertainties are residual with Type B having an unknown $p_T$ dependence. For each $y$-coordinate in each $p_T$ bin was varied by the same amount according to the uncertainty Type C.

The experimental sensitivity and observed limit on the number of dark photon candidates as a function of the assumed dark photon mass.

Limits on the $U-\gamma$ mixing parameters from PHENIX at 90%, 85% and 80% confidence levels.

$\pi^0$-h$^{\pm}$ $\sqrt{\langle p^2_{out} \rangle}$ values from fits to the $\Delta \phi$ correlations.

$\pi^0$ per trigger yields as a function of $p_{out}$.

$\pi^0$ per trigger yields as a function of $p_{out}$.

$\pi^0$ per trigger yields as a function of $p_{out}$.

Direct photon per trigger yields as a function of $p_{out}$.

$\pi^0$ per-trigger yields as a function of $\Delta \phi$.

Direct photon per-trigger yields as a function of $\Delta \phi$.

Gaussian widths from fits to the $p_{out}$ distributions.

Gaussian widths from fits to the $p_{out}$ distributions.

Frame-independent angular decay coefficient $\tilde{\lambda}$ as a function of transverse momentum for the four reference frames and three $p_T$ bins. The numbers in the CS frame for the $p_T$ = 2-3 GeV/$c$ bin are 90% confidence level upper limits.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The per-trigger yields are shown as a function of $\Delta\phi$ in several $p_T^{trig}$ $\otimes$ $p_T^{assoc}$ bins.

The $p_{out}$ distributions are shown for dihadron and direct photon-hadron correlations, binned in $p^{trig}_T$ in $x_E$.

The $p_{out}$ distributions are shown for dihadron and direct photon-hadron correlations, binned in $p^{trig}_T$ in $x_E$.

The $p_{out}$ distributions are shown for dihadron and direct photon-hadron correlations, binned in $p^{trig}_T$ in $x_E$.

The $p_{out}$ distributions are shown for dihadron and direct photon-hadron correlations, binned in $p^{trig}_T$ in $x_E$.

The $p_{out}$ distributions are shown for dihadron and direct photon-hadron correlations, binned in $p^{trig}_T$ in $x_E$.

The $p_{out}$ distributions are shown for dihadron and direct photon-hadron correlations, binned in $p^{trig}_T$ in $x_E$.

The $p_{out}$ distributions are shown for dihadron and direct photon-hadron correlations, binned in $p^{trig}_T$ in $x_E$.

The $p_{out}$ distributions are shown for dihadron and direct photon-hadron correlations, binned in $p^{trig}_T$ in $x_E$.

The $p_{out}$ distributions are shown for dihadron and direct photon-hadron correlations, binned in $p^{trig}_T$ in $x_E$.

The nonperturbative away-side jet widths as a function of $p^{trig}_T$ at $\sqrt{s}$ = 200 GeV for both direct photon-hadron and dihadron correlations.

The nonperturbative away-side jet widths as a function of $p^{trig}_T$ at $\sqrt{s}$ = 200 GeV for both direct photon-hadron and dihadron correlations.

The $p_{out}$ distributions are shown in several bins of $x_E$, integrated over a range of $p_T^{trig}$.

The $p_{out}$ distributions are shown in several bins of $x_E$, integrated over a range of $p_T^{trig}$.

The $p_{out}$ distributions are shown in several bins of $x_E$, integrated over a range of $p_T^{trig}$.

The $p_{out}$ distributions are shown in several bins of $x_E$, integrated over a range of $p_T^{trig}$.

The $p_{out}$ distributions are shown in several bins of $x_E$, integrated over a range of $p_T^{trig}$.

The $p_{out}$ distributions are shown in several bins of $x_E$, integrated over a range of $p_T^{trig}$.

The $p_{out}$ distributions are shown in several bins of $x_E$, integrated over a range of $p_T^{trig}$.

The $p_{out}$ distributions are shown in several bins of $x_E$, integrated over a range of $p_T^{trig}$.

The Gaussian widths of $p_{out}$ as a function of $x_E$ are shown for both $\pi^0$ and direct photon triggered correlations.

The nonperturbative Gaussian widths are shown as a function of $x_E$ for $\sqrt{s}$ = 510 GeV. The points for $\sqrt{s}$ = 200 GeV are given in figure 8 and its associated table.

The nonperturbative Gaussian widths are shown as a function of $x_E$ for $\sqrt{s}$ = 510 GeV. The points for $\sqrt{s}$ = 200 GeV are given in figure 8 and its associated table.

$J_{dA}$ versus $x^{frag}_{Au}$ for $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for different centrality classes.

$J_{dA}$ versus $x^{frag}_{Au}$ for $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for different centrality classes.

$J_{dA}$ versus $x^{frag}_{Au}$ for $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for different centrality classes.

$p_T$ versus $R_{dA}$ for $\pi^0$ and range of 3.0 < $\eta$ < 3.8.

PHENIX $\pi^0$ invariant cross-section for $p$+$p$ collisions for range 3.0 < $\eta$ < 3.8.

Direct photon cross section. (a) The invariant cross sections of the direct photon in $p+p$ [3, 4] and $d$+Au collisions. The $p+p$ fit result with the empirical parameterization described in the text is shown as well as NLO pQCD calculations, and the scaled $p+p$ fit is compared with the $d$+Au data. The closed and open symbols show the results from the virtual photon and $\pi_{0}$-tagging methods, respectively. The asterisk symbols show the result from the statistical subtraction method for $d$+Au data, overlapping with the virtual photon result in 3 < $p_{T}$ < 5 GeV/c. The values in the table are equal to this mean value. The bars and bands represent the point-to-point (ptp.) and $p_{T}$-correlated (cor.) uncertainties, respectively. (b) The $p+p$ data over the fit. The uncertainties of the fit due to both point-to-point (ptp.) and pT -correlated uncertainties of the data are summed quadratically, and the sum is shown as dotted lines. The NLO pQCD calculations divided by the fit are also shown.

Direct photon cross section. (a) The invariant cross sections of the direct photon in $p+p$ [3, 4] and $d$+Au collisions. The $p+p$ fit result with the empirical parameterization described in the text is shown as well as NLO pQCD calculations, and the scaled $p+p$ fit is compared with the $d$+Au data. The closed and open symbols show the results from the virtual photon and $\pi_{0}$-tagging methods, respectively. The asterisk symbols show the result from the statistical subtraction method for $d$+Au data, overlapping with the virtual photon result in 3 < $p_{T}$ < 5 GeV/c. The values in the table are equal to this mean value. The bars and bands represent the point-to-point (ptp.) and $p_{T}$-correlated (cor.) uncertainties, respectively. (b) The $p+p$ data over the fit. The uncertainties of the fit due to both point-to-point (ptp.) and pT -correlated uncertainties of the data are summed quadratically, and the sum is shown as dotted lines. The NLO pQCD calculations divided by the fit are also shown.

(a) The invariant cross sections of the direct photon in $p+p$ [3, 4] and $d$+Au collisions. The $p+p$ fit result with the empirical parameterization described in the text is shown as well as NLO pQCD calculations, and the scaled $p+p$ fit is compared with the $d$+Au data. The closed and open symbols show the results from the virtual photon and $\pi_{0}$-tagging methods, respectively. The asterisk symbols show the result from the statistical subtraction method for $d$+Au data, overlapping with the virtual photon result in 3 < $p_{T}$ < 5 GeV/c. The values in the table are equal to this mean value. The bars and bands represent the point-to-point (ptp.) and $p_{T}$-correlated (cor.) uncertainties, respectively. (b) The $p+p$ data over the fit. The uncertainties of the fit due to both point-to-point (ptp.) and pT -correlated uncertainties of the data are summed quadratically, and the sum is shown as dotted lines. The NLO pQCD calculations divided by the fit are also shown.

$R_{dA}$ ($d$+Au data/scaled $p+p$ fit). Nuclear modification factor for $d$+Au, $R_{dA}$, as a function of $p_{T}$ . The closed and open symbols show the results from the virtual- and real-photon measurements, respectively. The values in the table are equal to this mean value. The bars and bands represent the point-to-point (ptp.) and $p_{T}$-correlated (cor.) uncertainties, respectively. The box on the right shows the uncertainty of $T_{dA}$ for $d$+Au. The curves indicate the theoretical calculations [24] with different combinations of the CNM effects such as the Cronin enhancement, isospin effect, nuclear shadowing and initial state energy loss.

$R_{dA}$ ($d$+Au data/scaled $p+p$ fit). Nuclear modification factor for $d$+Au, $R_{dA}$, as a function of $p_{T}$ . The closed and open symbols show the results from the virtual- and real-photon measurements, respectively. The values in the table are equal to this mean value. The bars and bands represent the point-to-point (ptp.) and $p_{T}$-correlated (cor.) uncertainties, respectively. The box on the right shows the uncertainty of $T_{dA}$ for $d$+Au. The curves indicate the theoretical calculations [24] with different combinations of the CNM effects such as the Cronin enhancement, isospin effect, nuclear shadowing and initial state energy loss.

$R_{AA}$ (Au+Au MB data/scaled $p+p$ fit). Nuclear modification factors for Au+Au (MB) and $d$+Au as a function of $p_{T}$ . The triangle symbols show results from the (closed) virtual [1] and (open) real photon [5] measurements, respectively. The values in the table are equal to this mean value. The bars and bands represent the point-to-point (ptp.) and $p_{T}$-correlated (cor.) uncertainties, respectively. The (+) symbols for $R_{dA}$ for $p_{T}$ < 5 GeV/$c$ illustrates the difference in magnitude for $R_{AA}$ between Au+Au and $d$+Au collisions.

$R_{AA}$ (Au+Au MB data/scaled $p+p$ fit). Nuclear modification factors for Au+Au (MB) and $d$+Au as a function of $p_{T}$ . The triangle symbols show results from the (closed) virtual [1] and (open) real photon [5] measurements, respectively. The values in the table are equal to this mean value. The bars and bands represent the point-to-point (ptp.) and $p_{T}$-correlated (cor.) uncertainties, respectively. The (+) symbols for $R_{dA}$ for $p_{T}$ < 5 GeV/$c$ illustrates the difference in magnitude for $R_{AA}$ between Au+Au and $d$+Au collisions.

$p_T$ dependence of $v_2$ for charged hadrons for several centrality selections as indicated.

$p_T$ dependence of $v_4$ for charged hadrons for several centrality selections as indicated.

$v_2$ vs. $N_{part}$ and $v_4$ vs. $N_{part}$ for charged hadrons for several $p_T$ selections as indicated.

The ratio ${v_4}/{(v_2)^2}$ vs. $N_{part}$ for the same $p_T$ selections.

Invariant cross section of heavy-flavor electrons in $p$$+$$p$ collisions at $\sqrt{s_{_{NN}}}$ = 62.4 GeV. Candidate (inclusive), photonic, and heavy-flavor electron $v_{2}$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV for 20%-40% centrality. Fit parameters to p+p data in Fig 13.

Integrated invariant yield per binary collision vs $N_{coll}$ for heavy-flavor electrons in Au$+$Au collisions at $\sqrt{s_{_{NN}}}$ = 62.4 GeV for the indicated $p_{T}^{e}$ ranges. The data point at $N_{coll}$ =1 is for $p$$+$$p$ collisions.

Integrated invariant yield per binary collision vs $N_{coll}$ for heavy-flavor electrons in Au$+$Au collisions at $\sqrt{s_{_{NN}}}$ = 62.4 GeV for the indicated $p_{T}^{e}$ ranges. The data point at $N_{coll}$ =1 is for $p$$+$$p$ collisions.

The $R_{AA}$ for electrons from heavy-flavor decays in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV for the indicated centralities. The error bars (boxes) represent the statistical (systematic) uncertainties. The global uncertainty due to the uncertainty in Ncoll for each centrality is given by the box on the right side of each plot.

The $R_{AA}$ for electrons from heavy-flavor decays in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV for the indicated centralities. The error bars (boxes) represent the statistical (systematic) uncertainties. The global uncertainty due to the uncertainty in Ncoll for each centrality is given by the box on the right side of each plot.

The $R_{AA}$ for electrons from heavy-flavor decays in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV for the indicated centralities. The error bars (boxes) represent the statistical (systematic) uncertainties. The global uncertainty due to the uncertainty in Ncoll for each centrality is given by the box on the right side of each plot.

The $R_{AA}$ for electrons from heavy-flavor decays in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV for the indicated centralities. The error bars (boxes) represent the statistical (systematic) uncertainties. The global uncertainty due to the uncertainty in Ncoll for each centrality is given by the box on the right side of each plot.

The $R_{AA}$ values for electrons from heavy-flavor decay in Au+Au collisions at 62.4 GeV with the $R_{AA}$ results from $d$+Au, Cu+Cu, and Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV (data from [9, 21]). The error bars (boxes) represent the statistical (systematic) uncertainties. The $p_T$ ranges are as indicated: (a) 1 < $p_T$ < 3 GeV/c and (b) 3 < $p_T$ < 5 GeV/$c$.

The $R_{AA}$ values for electrons from heavy-flavor decay in Au+Au collisions at 62.4 GeV with the $R_{AA}$ results from $d$+Au, Cu+Cu, and Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV (data from [9, 21]). The error bars (boxes) represent the statistical (systematic) uncertainties. The $p_T$ ranges are as indicated: (a) 1 < $p_T$ < 3 GeV/c and (b) 3 < $p_T$ < 5 GeV/$c$.

Heavy-flavor electron RCP between centrality 0\%–20\% and 40\%–60\% in Au$+$Au collisions at $\sqrt{s_{_{NN}}}$ = 62.4 GeV. The error bars (boxes) show the statistical (systematic) uncertainties.

Candidate (inclusive), photonic, and heavy-flavor electron $v_{2}$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV for 20%-40% centrality.

Measured $\pi^0$ $R_{AA}$ as a function of $p_T$ for the 0-10% most central collisions at $\sqrt{s_{NN}}$ = 22.4 GeV in comparison to a jet quenching calculation. The error (tot.) includes the quadratic sum of the statistical uncertainties and the point-to-point uncorrelated and correlated systematic uncertainties.

Measured $\pi^0$ $R_{AA}$ as a function of $p_T$ for the 0-10% most central collisions at $\sqrt{s_{NN}}$ = 62.4 GeV in comparison to a jet quenching calculation. The error (tot.) includes the quadratic sum of the statistical uncertainties and the point-to-point uncorrelated and correlated systematic uncertainties.

Measured $\pi^0$ $R_{AA}$ as a function of $p_T$ for the 0-10% most central collisions at $\sqrt{s_{NN}}$ = 200.0 GeV in comparison to a jet quenching calculation. The error (tot.) includes the quadratic sum of the statistical uncertainties and the point-to-point uncorrelated and correlated systematic uncertainties.

The average $R_{AA}$ in the interval 2.5 < $p_T$ < 3.5 GeV/$c$ as a function of centrality for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 22.4 GeV. The error (sys.) includes the normalization and $<N_{coll}>$ uncertainties for a typical $N_{coll}$ uncertainty of 12%.

The average $R_{AA}$ in the interval 2.5 < $p_T$ < 3.5 GeV/$c$ as a function of centrality for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV. The error (sys.) includes the normalization and $<N_{coll}>$ uncertainties for a typical $N_{coll}$ uncertainty of 12%.

The average $R_{AA}$ in the interval 2.5 < $p_T$ < 3.5 GeV/$c$ as a function of centrality for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV. The error (sys.) includes the normalization and $<N_{coll}>$ uncertainties for a typical $N_{coll}$ uncertainty of 12%.

Direct photon $p_T$ spectra in different centrality bins.

Direct photon $p_T$ spectra after subtraction of the $N_{coll}$ scaled $p$+$p$ contribution in different centrality bins.

Integrated thermal photon yields as a function of $N_{part}$ for different lower $p_T$ integration limits.

Cross section of midrapidity charged-hadron production from $p$ + $p$ collisions at $\sqrt{s}$ = 62.4 GeV as a function of $p_T$ for positive and negative hadrons.

Double-helicity asymmetries and the statistical uncertainties as a function of $p_T$ for positive and negative inclusive charged hadrons from $p$ + $p$ collisions at $\sqrt{s}$ = 62.4 GeV.

Charge-combined, invariant cross section of muons from open heavy-flavor decays as a function of $p_T$ in $p$+$p$ collisions at $\sqrt{s}$ = 200 GeV at forward rapidity.