$\lambda_{\phi}$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with NRQCD predictions in the Helicity and Collins-Soper frames.

$\lambda_{\theta \phi}$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with NRQCD predictions in the Helicity and Collins-Soper frames.

$\lambda_{\theta \phi}$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with NRQCD predictions in the Helicity and Collins-Soper frames.

$\bar{\lambda}$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with CSM and NRQCD predictions in the Helicity frame.

$\bar{\lambda}$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with CSM and NRQCD predictions in the Helicity frame.

Invariant $F$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with CSM and NRQCD predictions in the Helicity frame.

Invariant $F$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with CSM and NRQCD predictions in the Helicity frame.

$b\bar{b}$ invariant yield measured via B meson decay to like-sign dimuons as a function of pair $p_{T}$. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Au collisions as a function of centrality and rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in 3He+Au collisions as a function of centrality and rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Al MinBias collisions as a function of pT. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Au MinBias collisions as a function of pT. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in 3He+Au MInBias collisions as a function of pT. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Al collisions as a function of pT and centrality at backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Al collisions as a function of pT and centrality at forward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Au collisions as a function of pT and centrality at backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Au collisions as a function of pT and centrality at forward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in 3He+Au collisions as a function of pT and centrality at backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in 3He+Au collisions as a function of pT and centrality at forward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Al most central collisions as a function of pT is directly compared with p+Au collisions (Data listed in Table16) at forward and backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in d+Au and p+Au most central collisions as a function of pT is directly compared at forward and backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column. Please see PRC 87, 034904 for d+Au data points.

J/psi nuclear modification in 3He+Au most central collisions as a function of pT is directly compared to p+Au collisions (Data listed in Table 16) at forward and backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Al as a function of N_coll. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Au collisions as a function of N_coll. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in 3He+Au as a function of N_coll. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Al collisions as a function of nuclear thickness (T_A). The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in p+Au collisions as a function of nuclear thickness (T_A). The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi nuclear modification in 3He+Au collisions as a function of nuclear thickness (T_A). The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi mean pT squared in p+Al collisions as a function of N_coll. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi mean pT squared in p+Au collisions as a function of N_coll. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi mean pT squared in 3He+Au collisions as a function of N_coll. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields as a function of rapidity in small systems. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields in p+Al collisions as a function of centrality and rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields in p+Au collisions as a function of centrality and apidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields in 3He+Au collisions as a function of centrality and rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields in p+Al collisions as a function of pT and centrality at forward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields in p+Al collisions as a function of pT and centrality at backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields in p+Au collisions as a function of pT and centrality at forward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields in p+Au collisions as a function of pT and centrality at backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields in 3He+Au collisions as a function of pT and centrality at forward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

J/psi invariant yields in 3He+Au collisions as a function of pT and centrality at backward rapidity. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

Charged hadron $dN_{ch}/d\eta$ as a function of pseudorapidity in high-multiplicity 0%-5% central $p$+Al collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Charged hadron $dN_{ch}/d\eta$ as a function of pseudorapidity in various mutiplicity classes of $^3$He+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Charged hadron $dN_{ch}/d\eta$ as a function of pseudorapidity in various mutiplicity classes of $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Charged hadron $dN_{ch}/d\eta$ as a function of pseudorapidity in various mutiplicity classes of $p$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Charged hadron $dN_{ch}/d\eta$ as a function of pseudorapidity in various mutiplicity classes of $p$+Al collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Midrapidity charged hadron $dN_{ch}/d\eta$ per participating quark pair ($N_{qp}$/2) as a function of the number of participating quarks ($N_{qp}$).

Midrapidity charged hadron $dN_{ch}/d\eta$ per participating quark pair ($N_{qp}$/2) as a function of the number of participating quarks ($N_{qp}$).

Midrapidity charged hadron $dN_{ch}/d\eta$ per participating quark pair ($N_{qp}$/2) as a function of the number of participating quarks ($N_{qp}$).

Midrapidity charged hadron $dN_{ch}/d\eta$ per participating quark pair ($N_{qp}$/2) as a function of the number of participating quarks ($N_{qp}$).

Midrapidity charged hadron $dN_{ch}/d\eta$ per participating quark pair ($N_{qp}$/2) as a function of the number of participating quarks ($N_{qp}$).

Elliptic flow $v_2$ as a function of pseudorapidity in high-multiplicity 0%-5% central $p$+Al, $p$+Au, $d$+Au, and $^3$He+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The centrality dependence of e+e− invariant mass spectra within the STAR acceptance from Au+Au collisions and U+U collisions for pair pT < 0.15 GeV/c. The vertical bars on data points depict the statistical uncertainties, while the systematic uncertainties are shown as gray boxes. The hadronic cocktail yields from U+U collisions are ∼5%–12% higher than those from Au+Au collisions in given centrality bins; thus only cocktails for Au+Au collisions are shown here as solid lines, with shaded bands representing the systematic uncertainties for clarity.

The centrality dependence of e+e− invariant mass spectra within the STAR acceptance from Au+Au collisions and U+U collisions for pair pT < 0.15 GeV/c. The vertical bars on data points depict the statistical uncertainties, while the systematic uncertainties are shown as gray boxes. The hadronic cocktail yields from U+U collisions are ∼5%–12% higher than those from Au+Au collisions in given centrality bins; thus only cocktails for Au+Au collisions are shown here as solid lines, with shaded bands representing the systematic uncertainties for clarity.

The centrality dependence of e+e− invariant mass spectra within the STAR acceptance from Au+Au collisions and U+U collisions for pair pT < 0.15 GeV/c. The vertical bars on data points depict the statistical uncertainties, while the systematic uncertainties are shown as gray boxes. The hadronic cocktail yields from U+U collisions are ∼5%–12% higher than those from Au+Au collisions in given centrality bins; thus only cocktails for Au+Au collisions are shown here as solid lines, with shaded bands representing the systematic uncertainties for clarity.

The corresponding ratios of data over cocktail for Figure 1.

The corresponding ratios of data over cocktail for Figure 1.

The corresponding ratios of data over cocktail for Figure 1.

The corresponding ratios of data over cocktail for Figure 1.

The corresponding ratios of data over cocktail for Figure 1.

The corresponding ratios of data over cocktail for Figure 1.

The e+e− pair pT distributions within the STAR acceptance for different mass regions in 60%–80% Au+Au and U+U collisions compared to cocktails. The systematic uncertainties of the data are shown as gray boxes. The gray bands depict the systematic uncertainties of the cocktails.

The e+e− pair pT distributions within the STAR acceptance for different mass regions in 60%–80% Au+Au and U+U collisions compared to cocktails. The systematic uncertainties of the data are shown as gray boxes. The gray bands depict the systematic uncertainties of the cocktails.

The e+e− pair pT distributions within the STAR acceptance for different mass regions in 60%–80% Au+Au and U+U collisions compared to cocktails. The systematic uncertainties of the data are shown as gray boxes. The gray bands depict the systematic uncertainties of the cocktails.

The e+e− pair pT distributions within the STAR acceptance for different mass regions in 60%–80% Au+Au and U+U collisions compared to cocktails. The systematic uncertainties of the data are shown as gray boxes. The gray bands depict the systematic uncertainties of the cocktails.

The e+e− pair pT distributions within the STAR acceptance for different mass regions in 60%–80% Au+Au and U+U collisions compared to cocktails. The systematic uncertainties of the data are shown as gray boxes. The gray bands depict the systematic uncertainties of the cocktails.

The e+e− pair pT distributions within the STAR acceptance for different mass regions in 60%–80% Au+Au and U+U collisions compared to cocktails. The systematic uncertainties of the data are shown as gray boxes. The gray bands depict the systematic uncertainties of the cocktails.

The low-pT (pT < 0.15 GeV/c) e+e− excess mass spectra (data − cocktail) within the STAR acceptance in 60%–80% for Au + Au and U + U collisions, compared with a broadened rho meson model calculation [8]. The contributions of rho and phi from the photonuclear process are shown, as are the contributions of photon-photon process from two models [33,34]. The model calculations are for Au + Au collisions in the corresponding centrality bins.

The low-pT (pT < 0.15 GeV/c) e+e− excess mass spectra (data − cocktail) within the STAR acceptance in 60%–80% for Au + Au and U + U collisions, compared with a broadened rho meson model calculation [8]. The contributions of rho and phi from the photonuclear process are shown, as are the contributions of photon-photon process from two models [33,34]. The model calculations are for Au + Au collisions in the corresponding centrality bins.

The low-pT (pT < 0.15 GeV/c) e+e− excess mass spectra (data − cocktail) within the STAR acceptance in 40%–60% for Au + Au and U + U collisions, compared with a broadened rho meson model calculation [8]. The contributions of rho and phi from the photonuclear process are shown, as are the contributions of photon-photon process from two models [33,34]. The model calculations are for Au + Au collisions in the corresponding centrality bins.

The low-pT (pT < 0.15 GeV/c) e+e− excess mass spectra (data − cocktail) within the STAR acceptance in 40%–60% for Au + Au and U + U collisions, compared with a broadened rho meson model calculation [8]. The contributions of rho and phi from the photonuclear process are shown, as are the contributions of photon-photon process from two models [33,34]. The model calculations are for Au + Au collisions in the corresponding centrality bins.

The centrality dependence of integrated excess yields in the mass regions of 0.4–0.76, 0.76–1.2, and 1.2–2.6 GeV=c2 in Au + Au and U + U collisions. The centrality dependence of hadronic cocktail yields in the mass region of 0.76–1.2 GeV/c2 in both collisions is also shown for comparison. The systematic uncertainties are shown as gray boxes.

The centrality dependence of integrated excess yields in the mass regions of 0.4–0.76, 0.76–1.2, and 1.2–2.6 GeV=c2 in Au + Au and U + U collisions. The centrality dependence of hadronic cocktail yields in the mass region of 0.76–1.2 GeV/c2 in both collisions is also shown for comparison. The systematic uncertainties are shown as gray boxes.

The centrality dependence of integrated excess yields in the mass regions of 0.4–0.76, 0.76–1.2, and 1.2–2.6 GeV=c2 in Au + Au and U + U collisions. The centrality dependence of hadronic cocktail yields in the mass region of 0.76–1.2 GeV/c2 in both collisions is also shown for comparison. The systematic uncertainties are shown as gray boxes.

The centrality dependence of integrated excess yields in the mass regions of 0.4–0.76, 0.76–1.2, and 1.2–2.6 GeV=c2 in Au + Au and U + U collisions. The centrality dependence of hadronic cocktail yields in the mass region of 0.76–1.2 GeV/c2 in both collisions is also shown for comparison. The systematic uncertainties are shown as gray boxes.

The centrality dependence of integrated excess yields in the mass regions of 0.4–0.76, 0.76–1.2, and 1.2–2.6 GeV=c2 in Au + Au and U + U collisions. The centrality dependence of hadronic cocktail yields in the mass region of 0.76–1.2 GeV/c2 in both collisions is also shown for comparison. The systematic uncertainties are shown as gray boxes.

The centrality dependence of integrated excess yields in the mass regions of 0.4–0.76, 0.76–1.2, and 1.2–2.6 GeV/c2 in Au + Au and U + U collisions. The centrality dependence of hadronic cocktail yields in the mass region of 0.76–1.2 GeV/c2 in both collisions is also shown for comparison. The systematic uncertainties are shown as gray boxes.

The p2T distributions of excess yields within the STAR acceptance in the mass regions of 0.4–0.76 GeV/c2 in 60%–80% Au + Au and U + U collisions. The systematic uncertainties are shown as gray boxes. The solid and dotted lines are exponential fits to the data in Au + Au and U + U collisions, respectively. The dot-dashed and dot-dot- dashed lines represent the p2T distributions for the photon-photon process from two models [33,34] within the STAR acceptance in 60%–80% Au + Au collisions. The dashed lines illustrate the corresponding p2T distributions for e+e− pairs from the model [33] traversing 1 fm in a constant magnetic field of 10^14 T perpendicular to the beam line.

The p2T distributions of excess yields within the STAR acceptance in the mass regions of 0.4–0.76 GeV/c2 in 60%–80% Au + Au and U + U collisions. The systematic uncertainties are shown as gray boxes. The solid and dotted lines are exponential fits to the data in Au + Au and U + U collisions, respectively. The dot-dashed and dot-dot- dashed lines represent the p2T distributions for the photon-photon process from two models [33,34] within the STAR acceptance in 60%–80% Au + Au collisions. The dashed lines illustrate the corresponding p2T distributions for e+e− pairs from the model [33] traversing 1 fm in a constant magnetic field of 10^14 T perpendicular to the beam line.

The p2T distributions of excess yields within the STAR acceptance in the mass regions of 0.76–1.2 GeV/c2 in 60%–80% Au + Au and U + U collisions. The systematic uncertainties are shown as gray boxes. The solid and dotted lines are exponential fits to the data in Au + Au and U + U collisions, respectively. The dot-dashed and dot-dot- dashed lines represent the p2T distributions for the photon-photon process from two models [33,34] within the STAR acceptance in 60%–80% Au + Au collisions. The dashed lines illustrate the corresponding p2T distributions for e+e− pairs from the model [33] traversing 1 fm in a constant magnetic field of 10^14 T perpendicular to the beam line.

The p2T distributions of excess yields within the STAR acceptance in the mass regions of 0.76–1.2 GeV/c2 in 60%–80% Au + Au and U + U collisions. The systematic uncertainties are shown as gray boxes. The solid and dotted lines are exponential fits to the data in Au + Au and U + U collisions, respectively. The dot-dashed and dot-dot- dashed lines represent the p2T distributions for the photon-photon process from two models [33,34] within the STAR acceptance in 60%–80% Au + Au collisions. The dashed lines illustrate the corresponding p2T distributions for e+e− pairs from the model [33] traversing 1 fm in a constant magnetic field of 10^14 T perpendicular to the beam line.

The p2T distributions of excess yields within the STAR acceptance in the mass regions of 1.2–2.6 GeV/c2 in 60%–80% Au + Au and U + U collisions. The systematic uncertainties are shown as gray boxes. The solid and dotted lines are exponential fits to the data in Au + Au and U + U collisions, respectively. The dot-dashed and dot-dot- dashed lines represent the p2T distributions for the photon-photon process from two models [33,34] within the STAR acceptance in 60%–80% Au + Au collisions. The dashed lines illustrate the corresponding p2T distributions for e+e− pairs from the model [33] traversing 1 fm in a constant magnetic field of 10^14 T perpendicular to the beam line.

The p2T distributions of excess yields within the STAR acceptance in the mass regions of 1.2–2.6 GeV/c2 in 60%–80% Au + Au and U + U collisions. The systematic uncertainties are shown as gray boxes. The solid and dotted lines are exponential fits to the data in Au + Au and U + U collisions, respectively. The dot-dashed and dot-dot- dashed lines represent the p2T distributions for the photon-photon process from two models [33,34] within the STAR acceptance in 60%–80% Au + Au collisions. The dashed lines illustrate the corresponding p2T distributions for e+e− pairs from the model [33] traversing 1 fm in a constant magnetic field of 10^14 T perpendicular to the beam line.

The corresponding sqrt(p2T) of excess yields. The vertical bars on data points are the combined statistical and systematic uncertainties.

The corresponding sqrt(p2T) of excess yields. The vertical bars on data points are the combined statistical and systematic uncertainties.

Data/simulation comparisons of the relative jet yields as functions of jet transverse momentum for the endcap

Data/simulation comparisons of the jet yields vs jet neutral energy fraction (NEF), shown separately for jets in different pseudorapidity ranges. The points represent the data (solid circle for barrel, open square and open cross for endcap jets at two different pseudorapidity ranges), and the histogram is the simulation.

Data/simulation comparisons of dijet yields as a function of invariant mass (left) for all the accepted events.

Data/simulation comparisons of dijet yields as a function of the pseudorapidity gap between the jets.

Data/simulation comparisons of dijet yields as a function of azimuthal opening angle between the jets.

Data/simulation comparisons of the underlying event $\delta M$ vs underlying event corrected dijet invariant mass.

The dijet $p_{T}$ imbalance distribution.

The jet invariant mass distribution.

The distributions of the parton $x_1$ and $x_2$, which has been weighted by the partonic $\hat{a}_{LL}$, from PYTHIA detector level simulations at $\sqrt{s}$ = 200 GeV for different jet pseudorapidity ranges.

The distributions of the parton $x_1$ and $x_2$, which has been weighted by the partonic $\hat{a}_{LL}$, from PYTHIA detector level simulations at $\sqrt{s}$ = 200 GeV for different jet pseudorapidity ranges.

The distributions of the parton $x_1$ and $x_2$, which has been weighted by the partonic $\hat{a}_{LL}$, from PYTHIA detector level simulations at $\sqrt{s}$ = 200 GeV for different jet pseudorapidity ranges.

$A_{LL}$ as a function of parton-level invariant mass for dijets with the East barrel-endcap.

$A_{LL}$ as a function of parton-level invariant mass for dijets with the West barrel-endcap.

$A_{LL}$ as a function of parton-level invariant mass for dijets with the endcap-endcap.

$\pi^0$ spectra from figure 3a from 0-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ spectra from figure 3a from 20-40% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ spectra from figure 3a from 40-60% semi-peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ spectra supplemental to figure 3a from 60-90% peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

Ratio of $\pi^0$ spectra from the PbSc to the combined spectra from figure 3b from 0-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\pi^0$ spectra from the PbGl to the combined spectra from figure 3b from 0-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\pi^0$ spectra from the PbSc to the combined spectra from figure 3c from 20-40% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\pi^0$ spectra from the PbGl to the combined spectra from figure 3c from 20-40% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\pi^0$ spectra from the PbSc to the combined spectra from figure 3d from 40-60% semi-peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\pi^0$ spectra from the PbGl to the combined spectra from figure 3d from 40-60% semi-peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

$\eta$ spectra from figure 3e from minimum bias Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 3e from 0-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 3e from 20-40% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 3e from 40-60% semi-peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra supplemental to figure 3e from 60-90% peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

Ratio of $\eta$ spectra from the PbSc to the combined spectra from figure 3f from 0-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\eta$ spectra from the PbGl to the combined spectra from figure 3f from 0-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\eta$ spectra from the PbSc to the combined spectra from figure 3g from 20-40% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\eta$ spectra from the PbGl to the combined spectra from figure 3g from 20-40% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\eta$ spectra from the PbSc to the combined spectra from figure 3h from 40-60% semi-peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

Ratio of $\eta$ spectra from the PbGl to the combined spectra from figure 3h from 40-60% semi-peripheral Au+Cu collisions. Type A uncertainties are uncorrelated point-to-point. Systematic noncorr. uncertainties are uncorrelated between PbSc and PbGl subsystems.

$\eta/\pi^0$ ratio from figure 4 from minimum bias Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 4 from 0-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 4 from 20-40% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 4 from 40-60% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio supplemental to figure 4 from 60-90% peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\pi^0$ $R_{AB}$ from figure 5a from 0-10% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AB}$ from figure 5b from 10-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AB}$ from figure 5c from 0-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ $R_{AB}$ from figure 5a from 0-20% central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AB}$ from figure 5d from 20-40% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ $R_{AB}$ from figure 5b from 20-40% semi-central Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AB}$ from figure 5e from 40-60% semi-peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ $R_{AB}$ from figure 5c from 40-60% semi-peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AB}$ from figure 5f from minimum bias Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ $R_{AB}$ from figure 5f from minimum bias Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AB}$ supplemental to figure 5 from 60-90% peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ $R_{AB}$ supplemental to figure 5 from 60-90% peripheral Cu+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AB}$ supplemental to figure 6a from 20-30% central Au+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AB}$ supplemental to figure 6b from 0-10% central Cu+Cu collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AB}$ supplemental to figure 6b from 40-50% semi-central Au+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ for $p_T>5$ GeV/$c$ in Cu+Au collisions as a function of centrality from figure 7a. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ for $p_T>5$ GeV/$c$ in Au+Au collisions as a function of centrality from figure 7a. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ for $p_T>10$ GeV/$c$ in Cu+Au collisions as a function of centrality from figure 7b. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ for $p_T>10$ GeV/$c$ in Au+Au collisions as a function of centrality from figure 7b. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.