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.

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

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

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

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

(Color online.) $J/\psi$ RAA versus pT for several centrality bins for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The statistical (systematic) uncertainties are shown with vertical bars (open boxes). The filled boxes about unity on the right show the size of the normalization uncertainty. PHENIX low-$p_{T}$ $J/\psi$ results [6] and STAR high-$p_{T}$ results in Cu+Cu collisions [15] are shown for comparison. The curves are from the predictions by Model I (Liu et al.) [11] and Model II (Zhao et al.) [12].

(Color online) $R_{AA}$ versus $N_{part}$ for high-$p_{T}$ $J/\psi$ in comparison with a) low-$p_{T}$ $J/\psi$ data from PHENIX [6], and high-$p_{T}$ $\pi^{\pm}$ from STAR [41, 42], and b) high-$p_{T}$ $J/\psi$ from CMS [8]. The statistical (systematic) uncertainties are shown in vertical bars (boxes). The shaded green band about unity shows the systematic uncertainties from $N_{bin}$ and the box about unity on the right shows the $R_{AA}$ normalization uncertainty from the statistical and global systematic uncertainties of the $p+p$ reference data. The percentage in the top horizontal axis of (a) and (b) refer to the collision centrality for high-$p_{T}$ $J/\psi$ data at the RHIC and LHC, respectively.

$\Eta/ \pi^0 ratios

$\pi^0$ nuclear modification factors

$\pi^0$ nuclear modification factors

$pi^0$ nuclear modification factors vs reaction plane, $\Delta \phi = 0^{\circ}-15^{\circ}$

$pi^0$ nuclear modification factors vs reaction plane, $\Delta \phi = 0^{\circ}-15^{\circ}$

$pi^0$ nuclear modification factors vs reaction plane, $\Delta \phi = 15^{\circ}-30^{\circ}$

$pi^0$ nuclear modification factors vs reaction plane, $\Delta \phi = 15^{\circ}-30^{\circ}$

$pi^0$ nuclear modification factors vs reaction plane, $\Delta \phi = 30^{\circ}-45^{\circ}$

$pi^0$ nuclear modification factors vs reaction plane, $\Delta \phi = 30^{\circ}-45^{\circ}$

$pi^0$ nuclear modification factors vs reaction plane, $\Delta \phi = 45^{\circ}-60^{\circ}$

$pi^0$ nuclear modification factors vs reaction plane, $\Delta \phi = 45^{\circ}-60^{\circ}$

$pi^0$ nuclear modification factors vs reaction plane, $\Delta \phi = 60^{\circ}-75^{\circ}$

$pi^0$ nuclear modification factors vs reaction plane, $\Delta \phi = 60^{\circ}-75^{\circ}$

$pi^0$ nuclear modification factors vs reaction plane, $\Delta \phi = 75^{\circ}-90^{\circ}$

$pi^0$ nuclear modification factors vs reaction plane, $\Delta \phi = 75^{\circ}-90^{\circ}$

Centrality dependence of $p_T$ ingetrated $R_{AA}(\Delta \phi)$

Centrality dependence of $p_T$ ingetrated $R_{AA}(\Delta \phi)$

Statistical correlation (in percentage) of the unfolded AFB between the mass bins in the muon channel for rapidity |y| of 0-1.

Statistical correlation (in percentage) of the unfolded AFB between the mass bins in the muon channel for rapidity |y| of 1-1.25.

Statistical correlation (in percentage) of the unfolded AFB between the mass bins in the muon channel for rapidity |y| of 1.25-1.5.

Statistical correlation (in percentage) of the unfolded AFB between the mass bins in the muon channel for rapidity |y| of 1.5-2.4.

Statistical correlation (in percentage) of the unfolded AFB between the mass bins in the electron channel for rapidity |y| of 0-1.

Statistical correlation (in percentage) of the unfolded AFB between the mass bins in the electron channel for rapidity |y| of 1-1.25.

Statistical correlation (in percentage) of the unfolded AFB between the mass bins in the electron channel for rapidity |y| of 1.25-1.5.

Statistical correlation (in percentage) of the unfolded AFB between the mass bins in the electron channel for rapidity |y| of 1.5-2.4.

Statistical correlation (in percentage) of the unfolded AFB between the mass bins in the electron+muon channel for rapidity |y| of 0-1.

Statistical correlation (in percentage) of the unfolded AFB between the mass bins in the electron+muon channel for rapidity |y| of 1-1.25.

Statistical correlation (in percentage) of the unfolded AFB between the mass bins in the electron+muon channel for rapidity |y| of 1.25-1.5.

Statistical correlation (in percentage) of the unfolded AFB between the mass bins in the electron+muon channel for rapidity |y| of 1.5-2.4.

Figure 2-a. Expected limit +1sigma.

Figure 2-a. Expected limit.

Figure 2-a. Expected limit -1sigma.

Figure 2-b. Observed limit +1sigma-th.

Figure 2-b. Observed limit.

Figure 2-b. Observed limit -1sigma-th.

Figure 2-b. Expected limit +1sigma.

Figure 2. Auxiliary material. Cross section for the Gbb and Gtt model.

Figure 3a. Auxiliary material. Theoretical uncertainty for the Gbb model.

Figure 3b. Auxiliary material. Theoretical uncertainty for the Gtt model.

Figure 4a. Auxiliary material. Number of MC events for the Gbb model.

Figure 4b. Auxiliary material. Number of MC events for the Gtt model.

Figure 5a. Auxiliary material. Acceptance for the Gbb model in the signal region SR4-L.

Figure 5b. Auxiliary material. Acceptance for the Gbb model in the signal region SR4-M.

Figure 5c. Auxiliary material. Acceptance for the Gbb model in the signal region SR4-T.

Figure 6a. Auxiliary material. Acceptance for the Gtt model in the signal region SR6-L.

Figure 6b. Auxiliary material. Acceptance for the Gtt model in the signal region SR6-T.

Figure 7a. Auxiliary material. Efficiency for the Gbb model in the signal region SR4-L.

Figure 7b. Auxiliary material. Efficiency for the Gbb model in the signal region SR4-M.

Figure 7c. Auxiliary material. Efficiency for the Gbb model in the signal region SR4-T.

Figure 8a. Auxiliary material. Efficiency for the Gtt model in the signal region SR6-L.

Figure 8b. Auxiliary material. Efficiency for the Gtt model in the signal region SR6-T.

Figure 9a. Auxiliary material. Experimental uncertainty for the Gbb model in the signal region SR4-L.

Figure 9b. Auxiliary material. Experimental uncertainty for the Gbb model in the signal region SR4-M.

Figure 9c. Auxiliary material. Experimental uncertainty for the Gbb model in the signal region SR4-T.

Figure 10a. Auxiliary material. Experimental uncertainty for the Gtt model in the signal region SR6-L.

Figure 10b. Auxiliary material. Experimental uncertainty for the Gtt model in the signal region SR6-T.

Figure 11a. Auxiliary material. Observed CLs for the Gbb model in the signal region SR4-L.

Figure 11b. Auxiliary material. Observed CLs for the Gbb model in the signal region SR4-M.

Figure 11c. Auxiliary material. Observed CLs for the Gbb model in the signal region SR4-T.

Figure 12a. Auxiliary material. Observed CLs for the Gtt model in the signal region SR6-L.

Figure 12b. Auxiliary material. Observed CLs for the Gtt model in the signal region SR6-T.

Figure 13a. Auxiliary material. Expected CLs for the Gbb model in the signal region SR4-L.

Figure 13b. Auxiliary material. Expected CLs for the Gbb model in the signal region SR4-M.

Figure 13c. Auxiliary material. Expected CLs for the Gbb model in the signal region SR4-T.

Figure 14a. Auxiliary material. Expected CLs for the Gtt model in the signal region SR6-L.

Figure 14b. Auxiliary material. Expected CLs for the Gtt model in the signal region SR6-T.

Figure 15a. Auxiliary material. Upper Limit for the Gbb model.

Figure 15b. Auxiliary material. Upper Limit for the Gtt model.

Figure 16a. Auxiliary material. Best Signal Region for the Gbb model.

Figure 16b. Auxiliary material. Best Signal Region for the Gtt model.

Normalized distributions of Tranverse Thrust Minor for 5 lower limit values of leading particle PT.

Normalized distributions of Tranverse Sphericity for 4 ranges of leading particle PT.

Normalized distributions of Tranverse Sphericity for 5 lower limit values of leading particle PT.

Mean Values of Thrust, Thrust Minor and Sphericity verses Multiplicity.

Mean Values of Thrust, Thrust Minor and Sphericity verses charged particle PT scalar sum.