Fully corrected assorted $\pi^{0}$-hadron conditional pair distributions for $d$+Au collisions centrality 0-88% and $p$+$p$ collisions. The trigger $\pi^0$s are within 5 < $p_{T,trig}$ < 10 GeV/$c$ and are correlated with hadrons with $p_{T,assoc}$ of 1.5-2 GeV/$c$, 2-3GeV/$c$, 3-4 GeV/$c$, and 4-5 GeV/$c$.

Near- and far-side widths as a function of $p_{T, assoc}$ for charged hadron azimuthal correlations from minimumbias $d$+Au collisions.

Near-side width, far-side width as a function of $p_{T,assoc}$ for a charged pion and neutral pion triggers from the $p_{T,trig}$ range of 5-10 GeV/$c$ in minimum bias $d$+Au collisions.

Near-side width, far-side width as a function of $p_{T,assoc}$ for a charged pion and neutral pion triggers from the $p_{T,trig}$ range of 5-10 GeV/$c$ in minimum bias $d$+Au collisions.

Near- and far-side conditional yields as a function of $p_{T, assoc}$ for charged hadron triggers (2.5−4 GeV/$c$) and associated charged hadrons from $d$+Au collisions.

Near- and far-side conditional yields as a function of $p_{T, assoc}$ for charged hadron triggers (4−6 GeV/$c$) and associated charged hadrons from $d$+Au collisions.

Near- and far-side conditional yields as a function of $p_{T, assoc}$ for charged pion triggers (5−10 GeV/$c$) and associated charged hadrons from minimum-bias $d$+Au collisions.

Near- and far-side conditional yields as a function of $p_{T, assoc}$ for neutral pion triggers (5−10 GeV/$c$) and associated charged hadrons from minimum-bias $d$+Au collisions.

Near- and far-side widths and conditional yields as a function of $N_{coll}$ for charged hadron triggers (2.5−4 GeV/$c$) and associated charged hadrons (1–2.5 GeV/$c$) from $d$+Au collisions.

Near- and far-side widths and conditional yields as a function of centrality for neutral pion triggers (5−10 GeV/$c$) and associated charged hadrons (2–3 GeV/$c$) from $d$+Au collisions.

The extracted $\sqrt{<j^2_T>}$ for minimum-bias $d$+Au collisions from all four dihadron correlations. The trigger $p_T$ ranges are 3 - 5 GeV/$c$, 5 - 10 GeV/$c$ and 5 - 10 GeV/$c$ for $h^{\pm}$ - $h^{\pm}$, $\pi^0$ - $h^{\pm}$ and $\pi^{\pm}$ - $h^{\pm}$ assorted-$p_T$ correlations.

The extracted $\sqrt{<j^2_T>}$ for minimum-bias $d$+Au collisions from all four dihadron correlations. The trigger $p_T$ ranges are 3 - 5 GeV/$c$, 5 - 10 GeV/$c$ and 5 - 10 GeV/$c$ for $h^{\pm}$ - $h^{\pm}$, $\pi^0$ - $h^{\pm}$ and $\pi^{\pm}$ - $h^{\pm}$ assorted-$p_T$ correlations.

The extracted $\sqrt{<j^2_T>}$ for minimum-bias $d$+Au collisions from all four dihadron correlations. The trigger $p_T$ ranges are 3 - 5 GeV/$c$, 5 - 10 GeV/$c$ and 5 - 10 GeV/$c$ for $h^{\pm}$ - $h^{\pm}$, $\pi^0$ - $h^{\pm}$ and $\pi^{\pm}$ - $h^{\pm}$ assorted-$p_T$ correlations.

The extracted $\sqrt{<j^2_T>}$ for minimum-bias $d$+Au collisions from all four dihadron correlations. The trigger $p_T$ ranges are 3 - 5 GeV/$c$, 5 - 10 GeV/$c$ and 5 - 10 GeV/$c$ for $h^{\pm}$ - $h^{\pm}$, $\pi^0$ - $h^{\pm}$ and $\pi^{\pm}$ - $h^{\pm}$ assorted-$p_T$ correlations.

$<sin^2(\phi_{jj})>$ for minimum bias $d$+Au collisions as function of associated particle $p_T$ for $\pi^{\pm}$ - $h^{\pm}$ and $\pi^0$ - $h^{\pm}$ correlations where the trigger particles have a $p_T$ between 5 - 10 GeV/$c$.

$<sin^2(\phi_{jj})>$ for minimum bias $d$+Au collisions as function of associated particle $p_T$ for $\pi^{\pm}$ - $h^{\pm}$ and $\pi^0$ - $h^{\pm}$ correlations where the trigger particles have a $p_T$ between 5 - 10 GeV/$c$.

$<sin^2(\phi_{jj})>$ for minimum bias $d$+Au collisions as function of trigger particle $p_T$ for $\pi^{\pm}$ - $h^{\pm}$ and $\pi^0$ - $h^{\pm}$ correlations.

$<sin^2(\phi_{jj})>$ for minimum bias $d$+Au collisions as function of trigger particle $p_T$ for $\pi^{\pm}$ - $h^{\pm}$ and $\pi^0$ - $h^{\pm}$ correlations.

Near-side and far-side $p_{out}$ distributions for minimum bias $d$+Au collisions obtained from $\pi^{\pm}$ - $h^{\pm}$ correlation. The trigger is 5 < $p_{T,trig}$ < 10 GeV/$c$, the associated particle is 0.5 < $p_{T,assoc}$ < 5.0 GeV/$c$.

Conditional yield as a function of $x_E$ for near-side and far-side for $\pi^{\pm}$ - $h^{\pm}$ and $\pi^0$ - $h^{\pm}$ from minimum-bias $d$+Au collisions.

Conditional yield as a function of $x_E$ for near-side and far-side for $\pi^{\pm}$ - $h^{\pm}$ and $\pi^0$ - $h^{\pm}$ from minimum-bias $d$+Au collisions.

Conditional yield as a function of $x_E$ for near-side and far-side correlations for $\pi^{\pm}$ - $h^{\pm}$ correlations from minimum-bias $d$+Au collisions. The trigger pions are 5 < $p_{T,trig}$ < 6 GeV/$c$.

Conditional yield as a function of $x_E$ for near-side and far-side correlation for $\pi^{\pm}$ - $h^{\pm}$ correlation for several different trigger $p_T$s for minimum-bias $d$+Au collisions.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ for $\pi^{\pm}$ - $h^{\pm}$ correlation for minimum-bias $d$+Au collisions.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ for $\pi^{\pm}$ - $h^{\pm}$ correlation for minimum-bias $d$+Au collisions.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ for $\pi^{\pm}$ - $h^{\pm}$ correlation for minimum-bias $d$+Au collisions.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ for $\pi^{\pm}$ - $h^{\pm}$ correlation for minimum-bias $d$+Au collisions.

The comparison of the $\sqrt{<j^2_T>}$ values between $d$+Au and $p$+$p$ for the $\pi^{\pm}$ - $h^{\pm}$ correlations and $\pi^0$ - $h^{\pm}$ correlations. The trigger pion range is 5 - 10 GeV/$c$.

The comparison of the $\sqrt{<j^2_T>}$ values between $d$+Au and $p$+$p$ for the $\pi^{\pm}$ - $h^{\pm}$ correlations and $\pi^0$ - $h^{\pm}$ correlations. The trigger pion range is 5 - 10 GeV/$c$.

Quadrature difference between minimum-bias $d$+Au and $p$+$p$ $<sin^2(\phi_{jj})>$ values.

Quadrature difference between minimum-bias $d$+Au and $p$+$p$ $<sin^2(\phi_{jj})>$ values.

The comparison of the $p_{out}$ distribution at the near-side and far-side between central $d$+Au collisions and $p$+$p$ collisions. Results are obtained for $\pi^{\pm}$ - $h^{\pm}$ correlations with the associated hadron range 0.5 < $p_{T,assoc}$ < 5 GeV/$c$ and trigger pion range of 5 < $p_{T,trig}$ < 10 GeV/$c$.

The comparison of the $x_E$ distribution from $\pi^{\pm}$ - $h^{\pm}$ correlation at the near-side and far-side between minimum bias $d$+Au collisions and $p$+$p$ collisions. The trigger $\pi^{\pm}$ are from 5 - 10 GeV/$c$.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ from $p$+$p$ collisions, triggers are $\pi^{\pm}$ from 5 - 10 GeV/$c$.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ from $p$+$p$ collisions, triggers are $\pi^{\pm}$ from 5 - 10 GeV/$c$.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ from $p$+$p$ collisions, triggers are $\pi^{\pm}$ from 5 - 10 GeV/$c$.

The fitted fractional change in $dN/d_{x_E}$ per unit $p_{T,trig}$ of the far-side conditional yield for different ranges of $x_E$.

Centrality dependent near and far $CY(p_T)$ for $\pi^{\pm}$ - $h^{\pm}$ correlations and $\pi^0$ - $h^{\pm}$ correlations.

Centrality dependent near and far $CY(p_T)$ for $\pi^{\pm}$ - $h^{\pm}$ correlations and $\pi^0$ - $h^{\pm}$ correlations.

Centrality dependence to the ratio of near and far $CY(p_T)$ for $d$+Au to $p$+$p$.

Centrality dependence to the ratio of near and far $CY(p_T)$ for $d$+Au to $p$+$p$.

Near- and far-side widths as a function of $p_{T, trig}$ for charged pion triggers and associated charged hadrons (2–4.5 Gev/$c$) from minimum-bias $d$+Au collisions.

Near- and far-side widths as a function of $p_{T, trig}$ for neutral pion triggers and associated charged hadrons (2–4.5 Gev/$c$) from minimum-bias $d$+Au collisions.

<p>Charged hadron spectra for centrality selected d+Au collisions.</p>

<p>Charged hadron spectra for centrality selected d+Au collisions.</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. All centralities.</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. Centrality A (Central).</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. Centrality B.</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. Centrality C.</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. Centrality D (Peripheral).</p>

<p>$R_{dAu}$ as a function of $p_T$.</p>

<p>$R_{dAu}$ as a function of $p_T$.</p>

<p>$R_{dAu}$ as a function of $p_T$.</p>

<p>$R_{dAu}$ as a function of $p_T$.</p>

<p>$R_{NAu}$ as a function of $p_T$. Centrality A (Peripheral).</p>

<p>$R_{NAu}$ as a function of $p_T$. Centrality B.</p>

<p>$R_{dAu}$ as a function of $p_T$. Centrality C.</p>

<p>$R_{dAu}$ as a function of $p_T$. Centrality D (Central).</p>

<p>$R_{dAu}$ values averaged in three momentum ranges as functions of $\nu−1$. Horizontal bars show the uncertainty in the value of $\nu$ for each centrality class.</p>

<p>$R_{NAu}$ values averaged in three momentum ranges as functions of $\nu−1$. Horizontal bars show the uncertainty in the value of $\nu$ for each centrality class.</p>

<p>Charged hadron spectra for centrality selected d+Au collisions.</p>

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.).

Data selection efficiency as a function of nuclear recoil energy

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

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;).

Nuclear modification factor $R_{AA}$ in 10% most central $Au$+$Au$ collisions.

Nuclear modification factor $R_{dA}$ for ($h^+$+$h^-$)/2 in minimum bias $d$+$Au$.

$R_{AA}$ in the 10% most central $Au$+$Au$ collisions.

$R_{dA}$ for ($h^+$+$h^-$)/2 measurement in $d$+$Au$.

$R_{dA}$ for $\pi^0$ measurement in $d$+$Au$.

Ratio of $R_{dA}$ in minimum bias $d$+$Au$ and $R_{pA}$ from neutron tagged $d$+$Au$ collisions for ($h^+$+$h^-$)/2.

Ratio of $R_{dA}$ in minimum bias $d$+$Au$ and $R_{pA}$ from neutron tagged $d$+$Au$ collisions for $\pi^0$

Additional information containing number of events which were used to reconstruct the numvers matching to Figure 1 and 2.

Additional information containing number of events which were used to reconstruct the numvers matching to Figure 1 and 2.

Additional information containing number of events which were used to reconstruct the numvers matching to Figure 1 and 2.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Au+Au collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 62.4 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 22.5 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 22.5 GeV Cu+Cu collisions.

The uncorrected multiplicity distributions of charged hadrons with 0.2 < $p_T$ < 2.0 GeV/$c$ for 200 GeV $p$+$p$ collisions.

The mean from the NBD fit as a function of $N_{part}$ for 200 GeV Au+Au collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 62.4 GeV Au+Au collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 200 GeV Cu+Cu collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 62.4 GeV Cu+Cu collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The mean from the NBD fit as a function of $N_{part}$ for 22.5 GeV Cu+Cu collisions over the range 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of centrality for 200 GeV Au+Au collisions in the range 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 200 GeV Au+Au collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 62.4 GeV Au+Au collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 200 GeV Cu+Cu collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 62.4 GeV Cu+Cu collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Fluctuations expressed as the scaled variance as a function of $N_{part}$ for 22.5 GeV Cu+Cu collisions for 0.2 < $p_T$ < 2.0 GeV/$c$.

Scaled variance for 200 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Au+Au collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 200 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

Scaled variance for 62.4 GeV Cu+Cu collisions plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Au+Au collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 200 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The inverse of the parameter $k_{NBD}$ from the Negative Binomial Distribution fits for 62.4 GeV Cu+Cu collisions. The fluctuations are plotted as a function of $N_{part}$ for several $p_T$ ranges.

The scaled variance as a funciton of $N_{part}$ for 200 GeV Au+Au collisions in the range 0.2 < $p_T$ < 2.0 GeV/$c$.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

The correlation of the clan model parameters $\bar{n_c}$ and $\bar{N_c}$ for all of the collision species measured as a function of centrality.

Fit results based on $k(\delta_{\eta})$=$1/{{2\alpha\xi}/{\delta_{\eta}}}$ ($\xi << \delta_{\eta}$).

$F_{p_T}$ (in percent) of non-random fluctuations as a function of the $p_T$ range over which $M_{p_T}$ is calculated, 0.2 GeV/$c$ < $p_T$ < $p^{max}_T$, for the 20-25% centrality class ($N_{part}$ = 181.6).

Invariant cross section of $\omega$ production in $p$+$p$ and $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV measured in $\omega \rightarrow \pi^0\pi^+\pi^-$ and $\omega \rightarrow \pi^0\gamma$ decay channels.

Invariant cross section of $\omega$ production in $p$+$p$ and $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV measured in $\omega \rightarrow \pi^0\pi^+\pi^-$ and $\omega \rightarrow \pi^0\gamma$ decay channels.

Invariant cross section of $\omega$ production in $p$+$p$ and $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV measured in $\omega \rightarrow \pi^0\pi^+\pi^-$ and $\omega \rightarrow \pi^0\gamma$ decay channels.

Measured $\omega/{\pi^0}$ vs $p_T$ in $p$+$p$ collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Measured $\omega/{\pi^0}$ vs $p_T$ in $p$+$p$ collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Measured $\omega/{\pi^0}$ vs $p_T$ in $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Measured $\omega/{\pi^0}$ vs $p_T$ in $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Measured $R_{dA}$ vs $p_T$ for neutral mesons in $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for minimum bias.

Measured $R_{dA}$ vs $p_T$ for neutral mesons in $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for minimum bias.

Measured $R_{dA}$ vs $p_T$ for neutral mesons in $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for 0-20% centrality.

Measured $R_{dA}$ vs $p_T$ for neutral mesons in $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV for 0-20% centrality.

Reconstructed $\omega$ mass measured in the $\pi^0\pi^+\pi^-$ decay channel vs $p_T$ in $p$+$p$ collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Reconstructed $\omega$ mass measured in the $\pi^0\pi^+\pi^-$ decay channel vs $p_T$ in $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Invariant $\pi^0$ yields at midrapidity as a function of $p_T$ for minimum bias and nine centralities in $Au\ +\ Au$ at $\sqrt{s_{NN}} = 200\ GeV$ [0%–10% (80%–92%) is most central (peripheral)]. The labels "uncorr." and "corr." include systematic errors that are uncorrelated and correlated point-to-point, respectively.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nuclear modification factor $R_{AA}(p_T)$ for $\pi^0$ in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. The error bars include all point-to-point experimental ($p+p$, Au+Au) errors. The label "sys(uncorr)" refers to the uncorrelated point-to-point uncertainty (which includes statistical errors), while "sys(norm)" is the normalization uncertainty.

Nucelar modification factor $R_{CP}$ for ($p+\bar{p}$)/2.

Nucelar modification factor $R_{CP}$ for $\pi^0$.

Charged hadron to $\pi^0$ ratio in central (0-10%) Au+Au collisions.

Charged hadron to $\pi^0$ ratio in peripheral (60-92%) Au+Au collisions. The peripheral data points are offset by +130 MeV/$c$ for clarity.

Systematic error on normalization a) common to both centralities b) that can move the $h$/$\pi^0$ ratios for each centrality independent of the other.

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$ spectra of charged hadrons for minimum bias collisions along with spectra for 9 centrality classes derived from the pseudo-rapidity region $|\eta|$ < 0.18. Stat. stands for statistical error, syst. stands for the systematic errors and occ. stands for occupancy error.

$p_T$ spectra of charged hadrons for minimum bias collisions along with spectra for 9 centrality classes derived from the pseudo-rapidity region $|\eta|$ < 0.18. Stat. stands for statistical error, syst. stands for the systematic errors and occ. stands for occupancy error.

Centrality dependence of $<p^{trunc}_T>$. Data points for truncated mean $p_T$ calculation for three $p_T$ ranges. The integral region is from 8 GeV/$c$ > $p_T$ > $p_{T_{min}}$ for the most peripheral bin. Since the data stop around 6.5 GeV/$c$, truncated mean $p_T$ is lower.

Ratio of charged hadron yields per nucleon-nucleon collision between central (0-10%) and peripheral (60-92%) Au + Au collisions. Stat. stands for statistical error, syst. stands for the systematic errors and occ. stands for occupancy error.

$R_{AA}$ for $(h^++h^-)/2$ and $\pi^0$ as function of $p_{T}$ for minimum bias and 9 centrality classes. Stat. stands for statistical error from Au+Au. The asymmetric uncertainty from the fit of the $\pi^0$ mass peak is denoted by fit. $\sigma_{ENBC}$ stands for error on number of binary collisions. Error from $pp$ parametrization that is not $p_T$ dependent. This error is not included in the calculation of the Upper bound of the systematic error and the Lower bound of the systematic error.

$R_{AA}$ for $(h^++h^-)/2$ and $\pi^0$ as function of $p_{T}$ for minimum bias and 9 centrality classes. Stat. stands for statistical error from Au+Au. The asymmetric uncertainty from the fit of the $\pi^0$ mass peak is denoted by fit. $\sigma_{ENBC}$ stands for error on number of binary collisions. Error from $pp$ parametrization that is not $p_T$ dependent. This error is not included in the calculation of the Upper bound of the systematic error and the Lower bound of the systematic error.

$R_{AA}$ for $(h^++h^-)/2$ and $\pi^0$ as function of $p_{T}$ for minimum bias and 9 centrality classes. Stat. stands for statistical error from Au+Au. The asymmetric uncertainty from the fit of the $\pi^0$ mass peak is denoted by fit. $\sigma_{ENBC}$ stands for error on number of binary collisions. Error from $pp$ parametrization that is not $p_T$ dependent. This error is not included in the calculation of the Upper bound of the systematic error and the Lower bound of the systematic error.

Charged hadron to $\pi^0$ ratios for minimum bias events and 9 centrality classes.

Charged hadron to $\pi^0$ ratios for minimum bias events and 9 centrality classes.

Charged hadron to $\pi^0$ ratios for minimum bias events and 9 centrality classes.

Charged hadron to $\pi^0$ ratios for minimum bias events and 9 centrality classes.

Charged hadron to $\pi^0$ ratios for minimum bias events and 9 centrality classes.

Charged hadron to $\pi^0$ ratios for minimum bias events and 9 centrality classes.

Au + Au yield integrated for $p_T$ > 4.5 (GeV/$c$) over the $N$ + $N$ yield, normalized using either $N_{coll}$ or $N_{part}$, plot as a function of $<N_{part}>$. Statistical and systematic error from Au+Au and $N$+$N$ reference includes: a) $NN$ reference error with 10.4% normalization error subtracted, this error is 11.2%, b) systematic error on Au+Au is 8.5% (correction), and c) 5% occupancy correction (it is centrality dependent). Error on $N_{coll}$+10.4% error on $N$+$N$ reference.

Au + Au yield integrated for $p_T$ > 4.5 (GeV/$c$) over the $N$ + $N$ yield, normalized using either $N_{coll}$ or $N_{part}$, plot as a function of $<N_{part}>$. Statistical and systematic error from Au+Au and $N$+$N$ reference includes: a) $NN$ reference error with 10.4% normalization error subtracted, this error is 11.2%, b) systematic error on Au+Au is 8.5% (correction), and c) 5% occupancy correction (it is centrality dependent). Error on $N_{coll}$+10.4% error on $N$+$N$ reference.

Invariant yields from the $p$+$p$ reference data. The normalization uncertainty is labeled norm, the uncertainty from the fit is labeled fit, and the uncertainty $\pi^0$ is from the $\pi^0$ ratio.

$\gamma p \rightarrow J/\psi p$ differential cross sections 10.36–11.44 GeV beam energy range, average $t$ and beam energy in bins of $t$. The first cross section uncertainties are statistical, and the second are systematic. The overall average beam energy is 10.82 GeV. There is an additional fully correlated systematic uncertainty of 19.5% on the total cross section, not included here.

Parameter $\alpha$ from the fit to the data.

Results of the measurements by PHENIX at $\sqrt{s_{NN}}$ = 200 GeV.

Results of the measurements by PHENIX at $\sqrt{s_{NN}}$ = 130 GeV.

Results of the measurements by PHENIX at $\sqrt{s_{NN}}$ = 19.6 GeV.

Ratios of measured quantities at 200 GeV/130 GeV and 200 GeV/19.6 GeV. The number of $N_P$ is the average between two energies.

Ratios of measured quantities at 200 GeV/130 GeV and 200 GeV/19.6 GeV. The number of $N_P$ is the average between two energies.

200, 130, and 19.6 GeV are RHIC average values, and 17.2, 8.7, and 4.8 GeV are SPS average values of $dN_{ch}$/$d\eta$/($0.5N_p$) at different $\sqrt{s_{NN}}$. An additional 5% error (recalc.) is added for collision energies 4.8 - 17.6 GeV for the uncertainty related to recalculation to the Center of Mass system.

200, 130, and 19.6 GeV are RHIC average values, and 17.2, 8.7, and 4.8 GeV are SPS average values of $dN_{ch}$/$d\eta$/($0.5N_p$) at different $\sqrt{s_{NN}}$. An additional 5% error (recalc.) is added for collision energies 4.8 - 17.6 GeV for the uncertainty related to recalculation to the Center of Mass system.

200, 130, and 19.6 GeV are RHIC average values, and 17.2, 8.7, and 4.8 GeV are SPS average values of $dN_{ch}$/$d\eta$/($0.5N_p$) at different $\sqrt{s_{NN}}$. An additional 5% error (recalc.) is added for collision energies 4.8 - 17.6 GeV for the uncertainty related to recalculation to the Center of Mass system.

200, 130, and 19.6 GeV are RHIC average values, and 17.2, 8.7, and 4.8 GeV are SPS average values of $dN_{ch}$/$d\eta$/($0.5N_p$) at different $\sqrt{s_{NN}}$. An additional 5% error (recalc.) is added for collision energies 4.8 - 17.6 GeV for the uncertainty related to recalculation to the Center of Mass system.

$dN/dy$ and $T$ for different centrality bins.

Blast wave model parameters as a function of centrality from fitting $\pi^±$, $K^±$, $p$, and $\bar{p}$ spectra.

$m_T$ spectra of $\phi$ mesons in different centrality bins.

Centrality dependence of the $\phi$ mass centroid and $\phi$ intrinsic width.

Minimum-bias $m_T$ spectra of $\phi$ for different subsystem combinations.

Minimum-bias $m_T$ spectra of $\phi$ for different subsystem combinations.

Minimum-bias $m_T$ spectra of $\phi$ for different subsystem combinations.

Minimum-bias $m_T$ spectra of $\phi$ for different subsystem combinations.

$\phi$ $m_T$ spectra in different centrality bins.

Centrality dependence of $\phi$ yield at midrapidity.

Centrality dependence of $\phi$ yield at midrapidity.

Centrality dependence of the inverse slope, $T$.

Centrality dependence of particle ratios for $K^+$/$\pi^+$ and $K^-$/$\pi^-$.

Centrality dependence of particle ratios for $\phi$/$\pi$ and $\phi$/$K$.

Centrality dependence of $\phi$ yields at different collision energies.

Centrality dependence of $\phi$ yields at different collision energies.

$\phi$/$\pi$ and $\phi$/$K^-$ ratios as a function of collision energies.

Inverse slope parameteres $T$ in $m_T$ spectra.

$N_{coll}$ scaled central to peripheral ratio $R_{CP}$ for $\phi$. The total systematic error is 19%, but 7% is uncorrelated between the $\phi$ and $p$ $R_{CP}$.

Collision centrality, energy, and system size dependence of shape parameters.

Collision centrality, energy, and system size dependence of shape parameters.

Collision centrality, energy, and system size dependence of shape parameters.

The ratio of the hadronic decay photon $v_2$ over inclusive photon $v_2$ ($v_2^{b.g}/v_2^{inclusive \gamma}$) compared with the direct photon excess ratio $R = (N_{direct \gamma} + N_{b.g})/N_{b.g}$.

The ratio of the hadronic decay photon $v_2$ over inclusive photon $v_2$ ($v_2^{b.g}/v_2^{inclusive \gamma}$) compared with the direct photon excess ratio $R = (N_{direct \gamma} + N_{b.g})/N_{b.g}$.

Measured widths of the near- and away-angle $\pi^0$ - $h^{\pm}$ correlation peaks for various trigger momenta.

Measured widths of the near- and away-angle $\pi^0$ - $h^{\pm}$ correlation peaks for various trigger momenta.

Measured widths of the near- and away-angle $\pi^0$ - $h^{\pm}$ correlation peaks for various trigger momenta.

The $\sigma_N$ and $\sigma_A$ values in GeV/$c$.

The $\sqrt{<p^2_{out}>}$ values in GeV/$c$.

The $\sqrt{<p^2_{out}>}$ values in GeV/$c$.

The $\sqrt{<p^2_{out}>}$ values in GeV/$c$.

The near and away side conditional yield per number of triggers for 1.4 < $p_{Ta}$ < 5.0 GeV/$c$.

The $\sqrt{<j^2_T>}$ values in GeV/$c$.

The $\sqrt{<j^2_T>}$ values in GeV/$c$.

The $\sqrt{<j^2_T>}$ values in GeV/$c$.

The $\sqrt{<j^2_T>}$ values in GeV/$c$.

Measured $x_E$ distributions associated with various transverse momenta of the trigger $\pi^0$.

Measured $x_E$ distributions associated with various transverse momenta of the trigger $\pi^0$.

Measured $x_E$ distributions associated with various transverse momenta of the trigger $\pi^0$.

Measured $x_E$ distributions associated with various transverse momenta of the trigger $\pi^0$.

Measured $x_E$ distributions associated with various transverse momenta of the trigger $\pi^0$.

Extracted values of $D(x)$ parameters according from the fit to the LEP data and power $n$ of the unmeasured final state parton spectra $\Sigma_q(\bar{p_T})$ extracted from the fit to the single inclusive $\pi^0$ invariant cross section for corresponding fragmentation and fixed values of $\sqrt{<k^2_T>}$ = 2.5 GeV/$c$.

The $\bar{x^-1_h}$ $<z_t>$ $\sqrt{<k^2_T>}$ $\sqrt{<k^2_T>}$ values as a function of $p_{Ta}$ for two different trigger $\pi^0$ transverse momentum bins.

The $\bar{x^-1_h}$ $<z_t>$ $\sqrt{<k^2_T>}$ $\sqrt{<k^2_T>}$ values as a function of $p_{Ta}$ for two different trigger $\pi^0$ transverse momentum bins.

Values of $\bar{x^-1_h}$ $<z_t>$ $\sqrt{<k^2_T>}$ and $\sqrt{<k^2_T>}$ for various trigger particle $p_{Tt}$ and associated momenta in 1.4 < $p_{Ta}$ < 5.0 GeV/$c$ region.

The $<z_t>$ and $\bar{x_h}$ values with $p_{Tt}$

The $<z_t>$ and $\bar{x_h}$ values with $p_{Ta}$ for two trigger $\pi^0$ momenta bins.

The $<z_t>$ and $\bar{x_h}$ values with $p_{Ta}$ for two trigger $\pi^0$ momenta bins.