$\Delta N_p$ multiplicity distributions in Au+Au collisions at $\sqrt{S_{NN}}=27$ GeV for 0-5 percent, 30-40 percent and 70-80 percent collision centralities at midrapidity.

$\Delta N_p$ multiplicity distributions in Au+Au collisions at $\sqrt{S_{NN}}=39$ GeV for 0-5 percent, 30-40 percent and 70-80 percent collision centralities at midrapidity.

$\Delta N_p$ multiplicity distributions in Au+Au collisions at $\sqrt{S_{NN}}=62.4$ GeV for 0-5 percent, 30-40 percent and 70-80 percent collision centralities at midrapidity.

$\Delta N_p$ multiplicity distributions in Au+Au collisions at $\sqrt{S_{NN}}=200$ GeV for 0-5 percent, 30-40 percent and 70-80 percent collision centralities at midrapidity.

Centrality dependence of the cumulants of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{S_{NN}}=7.7$ GeV.

Centrality dependence of the cumulants of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{S_{NN}}=11.5$ GeV.

Centrality dependence of the cumulants of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{S_{NN}}=19.6$ GeV.

Centrality dependence of the cumulants of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{S_{NN}}=27$ GeV.

Centrality dependence of the cumulants of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{S_{NN}}=39$ GeV.

Centrality dependence of the cumulants of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{S_{NN}}=62.4$ GeV.

Centrality dependence of the cumulants of $\Delta N_p$ distributions for Au+Au collisions at $\sqrt{S_{NN}}=200$ GeV.

Centrality dependence of $S\sigma$/Skellam and $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{S_{NN}}=7.7$ GeV.

Centrality dependence of $S\sigma$/Skellam and $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{S_{NN}}=11.5$ GeV.

Centrality dependence of $S\sigma$/Skellam and $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{S_{NN}}=19.6$ GeV.

Centrality dependence of $S\sigma$/Skellam and $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{S_{NN}}=27$ GeV.

Centrality dependence of $S\sigma$/Skellam and $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{S_{NN}}=39$ GeV.

Centrality dependence of $S\sigma$/Skellam and $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{S_{NN}}=62.4$ GeV.

Centrality dependence of $S\sigma$/Skellam and $\kappa\sigma^2$ for $\Delta N_p$ in Au+Au collisions at $\sqrt{S_{NN}}=200$ GeV.

Collision energy and centrality dependence of the net-proton $S\sigma$ and $\kappa\sigma^2$ from Au+Au and p+p collisions at RHIC.

Collision energy and centrality dependence of the net-proton $S\sigma$ and $\kappa\sigma^2$ from Au+Au and p+p collisions at RHIC.

Collision energy and centrality dependence of the net-proton $S\sigma$ and $\kappa\sigma^2$ from Au+Au and p+p collisions at RHIC.

Collision energy and centrality dependence of the net-proton $S\sigma$ and $\kappa\sigma^2$ from Au+Au and p+p collisions at RHIC.

Cumulants of net-proton distribution at $\sqrt{S_{NN}}=7.7$ GeV. (efficiency corrected).

Cumulants of net-proton distribution at $\sqrt{S_{NN}}=11.5$ GeV. (efficiency corrected).

Cumulants of net-proton distribution at $\sqrt{S_{NN}}=19.6$ GeV. (efficiency corrected).

Cumulants of net-proton distribution at $\sqrt{S_{NN}}=27$ GeV. (efficiency corrected).

Cumulants of net-proton distribution at $\sqrt{S_{NN}}=39$ GeV. (efficiency corrected).

Cumulants of net-proton distribution at $\sqrt{S_{NN}}=62.4$ GeV. (efficiency corrected).

Cumulants of net-proton distribution at $\sqrt{S_{NN}}=200$ GeV. (efficiency corrected).

Cumulants of proton distribution at $\sqrt{S_{NN}}=7.7$ GeV. (efficiency corrected).

Cumulants of proton distribution at $\sqrt{S_{NN}}=11.5$ GeV. (efficiency corrected).

Cumulants of proton distribution at $\sqrt{S_{NN}}=19.6$ GeV. (efficiency corrected).

Cumulants of proton distribution at $\sqrt{S_{NN}}=27$ GeV. (efficiency corrected).

Cumulants of proton distribution at $\sqrt{S_{NN}}=39$ GeV. (efficiency corrected).

Cumulants of proton distribution at $\sqrt{S_{NN}}=62.4$ GeV. (efficiency corrected).

Cumulants of proton distribution at $\sqrt{S_{NN}}=200$ GeV. (efficiency corrected).

Cumulants of anti-proton distribution at $\sqrt{S_{NN}}=7.7$ GeV. (efficiency corrected).

Cumulants of anti-proton distribution at $\sqrt{S_{NN}}=11.5$ GeV. (efficiency corrected).

Cumulants of anti-proton distribution at $\sqrt{S_{NN}}=19.6$ GeV. (efficiency corrected).

Cumulants of anti-proton distribution at $\sqrt{S_{NN}}=27$ GeV. (efficiency corrected).

Cumulants of anti-proton distribution at $\sqrt{S_{NN}}=39$ GeV. (efficiency corrected).

Cumulants of anti-proton distribution at $\sqrt{S_{NN}}=62.4$ GeV. (efficiency corrected).

Cumulants of anti-proton distribution at $\sqrt{S_{NN}}=200$ GeV. (efficiency corrected).

Comparison of data to Monte Carlo for the distributions of two-photon invariant mass (left) and energy for the higher (center) and lower (right) energy photon.

Invariant mass distribution for the two-photon system with 7 < pT < 8 GeV/c.

Signal fractions calculated within the “peak region” of 0.1 < $M_{\gamma\gamma}$ < 0.2 $GeV/c^{2}$. Fractions for the full dataset as well as the subsets of longitudinal and transverse polarizations are shown as a function of pT (left), and the fractions for the transversely polarized data are shown as a function of $x_{F}$ integrated over 5 < pT < 12 GeV/c (right).

Signal fractions calculated within the “peak region” of 0.1 < $M_{\gamma\gamma}$ < 0.2 $GeV/c^{2}$. Fractions for the full dataset as well as the subsets of longitudinal and transverse polarizations are shown as a function of pT (left), and the fractions for the transversely polarized data are shown as a function of $x_{F}$ integrated over 5 < pT < 12 GeV/c (right).

Upper panel, the $\pi^{0}$ cross section (blue markers) is shown compared with an NLO pQCD cal- culation with three options for the scale parameter. The lower panel presents the ratio of the data to the pT-scale theory curve, as well as the ratio of the 2pT-scale and pT/2-scale theory curves to the pT-scale curve.

Upper panel, the $\pi^{0}$ cross section (blue markers) is shown compared with an NLO pQCD cal- culation with three options for the scale parameter. The lower panel presents the ratio of the data to the pT-scale theory curve, as well as the ratio of the 2pT-scale and pT/2-scale theory curves to the pT-scale curve.

The ALL results (blue markers) are presented with the DSSV prediction and the GRSV prediction using the best fit to polarized DIS ($\Delta$ g = std) and the maximum and minimum allowed values for gluon polarization.

The ALL results (blue markers) are presented with the DSSV prediction and the GRSV prediction using the best fit to polarized DIS ($\Delta$ g = std) and the maximum and minimum allowed values for gluon polarization.

The ALL results (blue markers) are presented with the DSSV prediction and the GRSV prediction using the best fit to polarized DIS ($\Delta$ g = std) and the maximum and minimum allowed values for gluon polarization.

The AN results are plotted versus $x_{F}$ integrated over 5 < pT < 12 GeV/c (left panel) and versus pT integrated over 0.06 < |$x_{F}$| < 0.27 (right panel).

The AN results are plotted versus $x_{F}$ integrated over 5 < pT < 12 GeV/c (left panel) and versus pT integrated over 0.06 < |$x_{F}$| < 0.27 (right panel).

The AN results are plotted versus $x_{F}$ integrated over 5 < pT < 12 GeV/c (left panel) and versus pT integrated over 0.06 < |x_{F}| < 0.27 (right panel).

The AN results are plotted versus $x_{F}$ integrated over 5 < pT < 12 GeV/c (left panel) and versus pT integrated over 0.06 < |$x_{F}$| < 0.27 (right panel).

The K0S cross section for 5 jet ET ranges.

The LAMBDA cross section for 5 jet ET ranges.

The correlation differences $\Delta\langle A^{2}\rangle=\delta\langle A^{2}_{UD}\rangle-\delta\langle A^{2}_{ LR}\rangle$ and $\Delta\langle A_{+}A_{-}\rangle=\delta\langle A_{+}A_{-}\rangle_{ UD}-\delta\langle A_{+}A_{-}\rangle_{ LR}$ as a function of $p_{T}$. The data are from 20-40% central (upper) and 0-20% central (lower) Au+Au collisions.

The charge asymmetry correlations $\delta\langle A^{2}\rangle$ (left panels) and $\delta\langle A_{+}A_{-}\rangle$ (center panels), and correlation differences $\Delta\langle A^{2}\rangle=\delta\langle A^{2}_{ UD}\rangle-\delta\langle A^{2}_{ LR}\rangle$ and $\Delta\langle A_{+}A_{-}\rangle=\delta\langle A_{+}A_{-}\rangle_{ UD}-\delta\langle A_{+}A_{-}\rangle_{ LR}$ (right panels), as a function of the azimuthal elliptic anisotropy of high-$p_{T}$ ($p_{T}>2$~GeV/$c$) particles (upper panels) and low-$p_{T}$ ($p_{T}<2$~GeV/$c$) particles (lower panels). Data are from 20-40% Au+Au collisions. The particle $p_{T}$ range of $0.15<p_{T}<2$~GeV/$c$ is used for both EP construction and asymmetry calculation.

The charge asymmetry correlations $\delta\langle A^{2}\rangle$ (left panels) and $\delta\langle A_{+}A_{-}\rangle$ (center panels), and correlation differences $\Delta\langle A^{2}\rangle=\delta\langle A^{2}_{ UD}\rangle-\delta\langle A^{2}_{ LR}\rangle$ and $\Delta\langle A_{+}A_{-}\rangle=\delta\langle A_{+}A_{-}\rangle_{ UD}-\delta\langle A_{+}A_{-}\rangle_{ LR}$ (right panels), as a function of the azimuthal elliptic anisotropy of high-$p_{T}$ ($p_{T}>2$~GeV/$c$) particles (upper panels) and low-$p_{T}$ ($p_{T}<2$~GeV/$c$) particles (lower panels). Data are from 20-40% Au+Au collisions. The particle $p_{T}$ range of $0.15<p_{T}<2$~GeV/$c$ is used for both EP construction and asymmetry calculation.

The wedge size dependences of charge multiplicity asymmetry correlations (upper panel), and their differences between out-of-plane and in-plane, $\Delta\langle A^{2}\rangle$ and $\Delta\langle A_{+}A_{-}\rangle$ (lower panel) for 20-40% central Au+Au collisions.

The wedge size dependences of charge multiplicity asymmetry correlations (upper panel), and their differences between out-of-plane and in-plane, $\Delta\langle A^{2}\rangle$ and $\Delta\langle A_{+}A_{-}\rangle$ (lower panel) for 20-40% central Au+Au collisions.

Charge multiplicity asymmetry correlations as a function of the wedge location, $\phi_{ w}$, in 20-40% central Au+Au collisions. The wedge size is $30^{\circ}$. The curves are the characteristic $\cos(2\phi_{ w}$) to guide the eye.

The wedge size dependence of the difference between the same-sign and opposite-sign, $\Delta\langle A^{2}\rangle-\Delta\langle A_{+}A_{-}\rangle$.

The values of $\Delta\langle A^{2}\rangle-\Delta\langle A_{+}A_{-}\rangle$, scaled by $N_{part}$, as a function of the measured average elliptic anisotropy $<v^{obs}_{2}}>$ in Au+Au collisions. The centrality bin number is labeled by each data point, 0 for 70-80% up to 8 for 0-5%.

$\Delta=\Delta\langle A^{2}\rangle-\Delta\langle A_{+}A_{-}\rangle$ as a function of $v^{obs}_{2}$, the event-by-event elliptical anisotropy of particle distributions relative to the second-harmonic event plane reconstructed from TPC tracks (left panel) and the first harmonic event plane reconstructed from the ZDC-SMD neutron signals (middle panel), in 20-40% central Au+Au collisions. Right panel: Average $\Delta$ for events with $|v^{obs}_{2}|<0.04$ relative to the TPC event plane as a function of centrality.

$\Delta=\Delta\langle A^{2}\rangle-\Delta\langle A_{+}A_{-}\rangle$ as a function of $v^{obs}_{2}$, the event-by-event elliptical anisotropy of particle distributions relative to the second-harmonic event plane reconstructed from TPC tracks (left panel) and the first harmonic event plane reconstructed from the ZDC-SMD neutron signals (middle panel), in 20-40% central Au+Au collisions. Right panel: Average $\Delta$ for events with $|v^{obs}_{2}|<0.04$ relative to the TPC event plane as a function of centrality.

$\Delta=\Delta\langle A^{2}\rangle-\Delta\langle A_{+}A_{-}\rangle$ as a function of $v^{obs}_{2}$, the event-by-event elliptical anisotropy of particle distributions relative to the second-harmonic event plane reconstructed from TPC tracks (left panel) and the first harmonic event plane reconstructed from the ZDC-SMD neutron signals (middle panel), in 20-40% central Au+Au collisions. Right panel: Average $\Delta$ for events with $|v^{obs}_{2}|<0.04$ relative to the TPC event plane as a function of centrality.

Left panel: $\Delta\langle A^{2}\rangle$ and $\Delta\langle A_{+}A_{-}\rangle$ as a function of $v^{obs}_{2}$ with a random EP in 20-40% Au+Au collisions. The lines depict the results with reconstructed EP lower right panel for comparison. Middle panel: The differences in $\Delta\langle A^{2}\rangle$ and $\Delta\langle A_{+}A_{-}\rangle$ between the reconstructed EP and random EP results, respectively. Right panel: $\Delta=\Delta\langle A^{2}\rangle-\Delta\langle A_{+}A_{-}\rangle$ as a function of $v^{obs}_{2}$ with a random EP (open triangles).

Left panel: $\Delta\langle A^{2}\rangle$ and $\Delta\langle A_{+}A_{-}\rangle$ as a function of $v^{obs}_{2}$ with a random EP in 20-40% Au+Au collisions. The lines depict the results with reconstructed EP lower right panel for comparison. Middle panel: The differences in $\Delta\langle A^{2}\rangle$ and $\Delta\langle A_{+}A_{-}\rangle$ between the reconstructed EP and random EP results, respectively. Right panel: $\Delta=\Delta\langle A^{2}\rangle-\Delta\langle A_{+}A_{-}\rangle$ as a function of $v^{obs}_{2}$ with a random EP (open triangles).

Left panel: $\Delta\langle A^{2}\rangle$ and $\Delta\langle A_{+}A_{-}\rangle$ as a function of $v^{obs}_{2}$ with a random EP in 20-40% Au+Au collisions. The lines depict the results with reconstructed EP lower right panel for comparison. Middle panel: The differences in $\Delta\langle A^{2}\rangle$ and $\Delta\langle A_{+}A_{-}\rangle$ between the reconstructed EP and random EP results, respectively. Right panel: $\Delta=\Delta\langle A^{2}\rangle-\Delta\langle A_{+}A_{-}\rangle$ as a function of $v^{obs}_{2}$ with a random EP (open triangles).

The correlation differences, $\langle A^{2}_{ LR}\rangle-\langle A^{2}_{ UD}\rangle$ and $\langle A_{+}A_{-}\rangle_{ LR}-\langle A_{+}A_{-}\rangle_{ UD}$, scaled by the number of participants $(\pi/4)^2 N_{part}$. Also shown are $\langle \cos(\phi_{\alpha}+\phi_{\beta}-2\phi_{c})\rangle/v_{2,c} \approx \langle \cos(\phi_{\alpha}+\phi_{\beta}-2\psi_{ RP}) \rangle$ of same- and opposite-sign particle pairs ($\alpha$ and $\beta$) calculated from particles used for the charge asymmetry correlations with particle $c$ from those used for EP construction. The upper panel shows the default results where particles are divided according to $\eta$; the lower panel shows results with particles divided randomly into two halves and compared to published correlators.

The correlation differences, $\langle A^{2}_{ LR}\rangle-\langle A^{2}_{ UD}\rangle$ and $\langle A_{+}A_{-}\rangle_{ LR}-\langle A_{+}A_{-}\rangle_{ UD}$, scaled by the number of participants $(\pi/4)^2 N_{part}$. Also shown are $\langle \cos(\phi_{\alpha}+\phi_{\beta}-2\phi_{c})\rangle/v_{2,c} \approx \langle \cos(\phi_{\alpha}+\phi_{\beta}-2\psi_{ RP}) \rangle$ of same- and opposite-sign particle pairs ($\alpha$ and $\beta$) calculated from particles used for the charge asymmetry correlations with particle $c$ from those used for EP construction. The upper panel shows the default results where particles are divided according to $\eta$; the lower panel shows results with particles divided randomly into two halves and compared to published correlators.

Upper panel: Positively charged particle multiplicity distributions versus the azimuthal angle, $\phi$, in 30-40% central Au+Au collisions. The data are shown separately for $\eta>0$ and $\eta<0$. The magnetic polarities of the STAR magnet have been summed and the $p_{T}$ is integrated over the range $0.15<p_{T}<2$~GeV/$c$. The inverse of these distributions, properly normalized, were used to correct for the track efficiency versus $\phi$. Lower panel: Reconstructed event-plane azimuthal angle distributions in 30-40% central Au+Au collisions. Particles within $0.15 < p_{T} < 2$~GeV/$c$ from $\eta<0$ and $\eta>0$ were used separately to reconstruct the EP.

Upper panel: Positively charged particle multiplicity distributions versus the azimuthal angle, $\phi$, in 30-40% central Au+Au collisions. The data are shown separately for $\eta>0$ and $\eta<0$. The magnetic polarities of the STAR magnet have been summed and the $p_{T}$ is integrated over the range $0.15<p_{T}<2$~GeV/$c$. The inverse of these distributions, properly normalized, were used to correct for the track efficiency versus $\phi$. Lower panel: Reconstructed event-plane azimuthal angle distributions in 30-40% central Au+Au collisions. Particles within $0.15 < p_{T} < 2$~GeV/$c$ from $\eta<0$ and $\eta>0$ were used separately to reconstruct the EP.

The relative charge asymmetry correlations, $\langle A_{+}A_{-}\rangle_{ UD}/\langle A_{+}A_{-}\rangle_{ LR}$, as a function of the number of participants, $N_{part}$, for four combinations of $\eta$ ranges used for EP reconstruction and asymmetry calculation.

The asymmetry correlations, $\langle A^{2}_{+,{ UD}}\rangle$ (left panel), $\langle A^{2}_{-,{ UD}}\rangle$ (middle panel), and $\langle A_{+}A_{-}\rangle_{ UD}$ (right panel), multiplied by the number of participants $N_{part}$, before and after the corrections for the $\phi$-dependent acceptance $\times$ efficiency. The colored curves are the corresponding statistical fluctuations and detector effects. The results are separated for $\eta>0$ and $\eta<0$. Only $\eta>0$ is shown for $\langle A^{2}_{+,{ UD}}\rangle$ and $\eta<0$ for $\langle A^{2}_{-,{ UD}}\rangle$ for clarity.

The asymmetry correlations, $\langle A^{2}_{+,{ UD}}\rangle$ (left panel), $\langle A^{2}_{-,{ UD}}\rangle$ (middle panel), and $\langle A_{+}A_{-}\rangle_{ UD}$ (right panel), multiplied by the number of participants $N_{part}$, before and after the corrections for the $\phi$-dependent acceptance $\times$ efficiency. The colored curves are the corresponding statistical fluctuations and detector effects. The results are separated for $\eta>0$ and $\eta<0$. Only $\eta>0$ is shown for $\langle A^{2}_{+,{ UD}}\rangle$ and $\eta<0$ for $\langle A^{2}_{-,{ UD}}\rangle$ for clarity.

The asymmetry correlations, $\langle A^{2}_{+,{ UD}}\rangle$ (left panel), $\langle A^{2}_{-,{ UD}}\rangle$ (middle panel), and $\langle A_{+}A_{-}\rangle_{ UD}$ (right panel), multiplied by the number of participants $N_{part}$, before and after the corrections for the $\phi$-dependent acceptance $\times$ efficiency. The colored curves are the corresponding statistical fluctuations and detector effects. The results are separated for $\eta>0$ and $\eta<0$. Only $\eta>0$ is shown for $\langle A^{2}_{+,{ UD}}\rangle$ and $\eta<0$ for $\langle A^{2}_{-,{ UD}}\rangle$ for clarity.

Left panel: Statistical fluctuation and detector effects in the charge asymmetry variance (multiplied by the number of participants $N_{part}$) from the $\eta<0$ region. The dotted curve shows the $1/N$ approximation, the dashed curve shows the statistical fluctuations $\langle A^{2}_{-,{ LR,stat}}\rangle$ obtained by the (50-50) method (see the text). The solid curve connecting the solid points shows the net effect of the statistical fluctuations and detector non-uniformities, $\langle A^{2}_{-,{ LR,stat+det}}\rangle$, obtained by the adding-$\pi$ method, and the crosses shows the same but obtained from the ``scramble'' method. See the text for details. Middle panel: The statistical fluctuation effects relative to the $1/N$ approximation (thin curves) and the effects due to the imperfect detector (thick curves). The dashed curves are for the region $\eta>0$ and the solid curves are for $\eta<0$. Right panel: The ratios of the statistical fluctuation and detector effects $\langle A^{2}_{+,{ LR,stat+det}}\rangle/\langle A^{2}_{-,{ LR,stat+det}}\rangle$ and $\langle A^{2}_{+,{ UD,stat+det}}\rangle/\langle A^{2}_{+,{ LR,stat+det}}\rangle$, separately for the $\eta>0$ and $\eta<0$ regions. The $\phi$-acceptance correction was applied for the results in this figure.

Left panel: Statistical fluctuation and detector effects in the charge asymmetry variance (multiplied by the number of participants $N_{part}$) from the $\eta<0$ region. The dotted curve shows the $1/N$ approximation, the dashed curve shows the statistical fluctuations $\langle A^{2}_{-,{ LR,stat}}\rangle$ obtained by the (50-50) method (see the text). The solid curve connecting the solid points shows the net effect of the statistical fluctuations and detector non-uniformities, $\langle A^{2}_{-,{ LR,stat+det}}\rangle$, obtained by the adding-$\pi$ method, and the crosses shows the same but obtained from the ``scramble'' method. See the text for details. Middle panel: The statistical fluctuation effects relative to the $1/N$ approximation (thin curves) and the effects due to the imperfect detector (thick curves). The dashed curves are for the region $\eta>0$ and the solid curves are for $\eta<0$. Right panel: The ratios of the statistical fluctuation and detector effects $\langle A^{2}_{+,{ LR,stat+det}}\rangle/\langle A^{2}_{-,{ LR,stat+det}}\rangle$ and $\langle A^{2}_{+,{ UD,stat+det}}\rangle/\langle A^{2}_{+,{ LR,stat+det}}\rangle$, separately for the $\eta>0$ and $\eta<0$ regions. The $\phi$-acceptance correction was applied for the results in this figure.

Left panel: Statistical fluctuation and detector effects in the charge asymmetry variance (multiplied by the number of participants $N_{part}$) from the $\eta<0$ region. The dotted curve shows the $1/N$ approximation, the dashed curve shows the statistical fluctuations $\langle A^{2}_{-,{ LR,stat}}\rangle$ obtained by the (50-50) method (see the text). The solid curve connecting the solid points shows the net effect of the statistical fluctuations and detector non-uniformities, $\langle A^{2}_{-,{ LR,stat+det}}\rangle$, obtained by the adding-$\pi$ method, and the crosses shows the same but obtained from the ``scramble'' method. See the text for details. Middle panel: The statistical fluctuation effects relative to the $1/N$ approximation (thin curves) and the effects due to the imperfect detector (thick curves). The dashed curves are for the region $\eta>0$ and the solid curves are for $\eta<0$. Right panel: The ratios of the statistical fluctuation and detector effects $\langle A^{2}_{+,{ LR,stat+det}}\rangle/\langle A^{2}_{-,{ LR,stat+det}}\rangle$ and $\langle A^{2}_{+,{ UD,stat+det}}\rangle/\langle A^{2}_{+,{ LR,stat+det}}\rangle$, separately for the $\eta>0$ and $\eta<0$ regions. The $\phi$-acceptance correction was applied for the results in this figure.

The asymmetry correlations, $\langle A^{2}_{+,{ UD}}\rangle$ (left panel), $\langle A^{2}_{-,{ UD}}\rangle$ (middle panel), and $\langle A_{+}A_{-}\rangle_{ UD}$ (right panel), multiplied by the number of participants $N_{part}$, separately for $\eta>0$ and $\eta<0$. The star and triangle depict the results from $d$+Au collisions. The net effects of the statistical fluctuations and detector non-uniformities are shown as the curves.

The asymmetry correlations, $\langle A^{2}_{+,{ UD}}\rangle$ (left panel), $\langle A^{2}_{-,{ UD}}\rangle$ (middle panel), and $\langle A_{+}A_{-}\rangle_{ UD}$ (right panel), multiplied by the number of participants $N_{part}$, separately for $\eta>0$ and $\eta<0$. The star and triangle depict the results from $d$+Au collisions. The net effects of the statistical fluctuations and detector non-uniformities are shown as the curves.

The asymmetry correlations, $\langle A^{2}_{+,{ UD}}\rangle$ (left panel), $\langle A^{2}_{-,{ UD}}\rangle$ (middle panel), and $\langle A_{+}A_{-}\rangle_{ UD}$ (right panel), multiplied by the number of participants $N_{part}$, separately for $\eta>0$ and $\eta<0$. The star and triangle depict the results from $d$+Au collisions. The net effects of the statistical fluctuations and detector non-uniformities are shown as the curves.

The statistical and detector effect subtracted asymmetry correlations, $\delta\langle A^{2}_{+,0^{\circ}\pm\Delta\phi_{w}}\rangle$ (left panel), $\delta\langle A^{2}_{-,0^{\circ}\pm\Delta\phi_{w}}\rangle$ (middle panel), and $\delta\langle A_{+}A_{-}\rangle_{ LR}$ (right panel), multiplied by the number of participants $N_{part}$, before and after the corrections for the $\phi$-dependent acceptance $\times$ efficiency.

The statistical and detector effect subtracted asymmetry correlations, $\delta\langle A^{2}_{+,0^{\circ}\pm\Delta\phi_{w}}\rangle$ (left panel), $\delta\langle A^{2}_{-,0^{\circ}\pm\Delta\phi_{w}}\rangle$ (middle panel), and $\delta\langle A_{+}A_{-}\rangle_{ LR}$ (right panel), multiplied by the number of participants $N_{part}$, before and after the corrections for the $\phi$-dependent acceptance $\times$ efficiency.

The statistical and detector effect subtracted asymmetry correlations, $\delta\langle A^{2}_{+,0^{\circ}\pm\Delta\phi_{w}}\rangle$ (left panel), $\delta\langle A^{2}_{-,0^{\circ}\pm\Delta\phi_{w}}\rangle$ (middle panel), and $\delta\langle A_{+}A_{-}\rangle_{ LR}$ (right panel), multiplied by the number of participants $N_{part}$, before and after the corrections for the $\phi$-dependent acceptance $\times$ efficiency.

Dynamical charge asymmetry variance after removing the statistical and detector effects, $\delta\langle A^{2}_{\pm,{ LR}}\rangle=\langle A^{2}_{\pm,{ LR}}\rangle-\langle A^{2}_{\pm,{ LR,stat+det}}\rangle$, scaled by the number of participants $N_{part}$. The results are consistent between positive and negative $\eta$ regions and between positive and negative charges.

The second-harmonic event-plane resolution of sub-events ($\eta<0$ and $\eta>0$) as a function of the number of participants.

Charge multiplicity asymmetry correlations, $\delta\langle A^{2}\rangle$ (left panel) and $\delta\langle A_{+}A_{-}\rangle$ (middle panel), and their differences between UD\ and LR\ (right panel) as a function of the event-plane resolution $\epsilon_{ EP}(f)=\sqrt{\mean{\cos2(\psi_{{ EP},\eta>0}(f)-\psi_{{ EP},\eta<0}(f))}}$ in 20-30% central Au+Au collisions. The solid lines are free linear fits to the data, while the dashed lines are linear fits with the intercept fixed to zero at $\epsilon_{ EP}(f)=0$.

Charge multiplicity asymmetry correlations, $\delta\langle A^{2}\rangle$ (left panel) and $\delta\langle A_{+}A_{-}\rangle$ (middle panel), and their differences between UD\ and LR\ (right panel) as a function of the event-plane resolution $\epsilon_{ EP}(f)=\sqrt{\mean{\cos2(\psi_{{ EP},\eta>0}(f)-\psi_{{ EP},\eta<0}(f))}}$ in 20-30% central Au+Au collisions. The solid lines are free linear fits to the data, while the dashed lines are linear fits with the intercept fixed to zero at $\epsilon_{ EP}(f)=0$.

TCharge multiplicity asymmetry correlations, $\delta\langle A^{2}\rangle$ (left panel) and $\delta\langle A_{+}A_{-}\rangle$ (middle panel), and their differences between UD\ and LR\ (right panel) as a function of the event-plane resolution $\epsilon_{ EP}(f)=\sqrt{\mean{\cos2(\psi_{{ EP},\eta>0}(f)-\psi_{{ EP},\eta<0}(f))}}$ in 20-30% central Au+Au collisions. The solid lines are free linear fits to the data, while the dashed lines are linear fits with the intercept fixed to zero at $\epsilon_{ EP}(f)=0$.

The TPC second-harmonic (upper left) and ZDC-SMD first harmonic (upper right) event-plane resolution as a function of $v^{obs}_{2}$ in 20-40% central Au+Au collisions. The TPC second-harmonic (lower left) and ZDC-SMD first harmonic (lower right) event-plane resolution for events with $|v^{obs}_{2}|<0.04$ as a function of centrality.

The TPC second-harmonic (upper left) and ZDC-SMD first harmonic (upper right) event-plane resolution as a function of $v^{obs}_{2}$ in 20-40% central Au+Au collisions. The TPC second-harmonic (lower left) and ZDC-SMD first harmonic (lower right) event-plane resolution for events with $|v^{obs}_{2}|<0.04$ as a function of centrality.

The TPC second-harmonic (upper left) and ZDC-SMD first harmonic (upper right) event-plane resolution as a function of $v^{obs}_{2}$ in 20-40% central Au+Au collisions. The TPC second-harmonic (lower left) and ZDC-SMD first harmonic (lower right) event-plane resolution for events with $|v^{obs}_{2}|<0.04$ as a function of centrality.

The TPC second-harmonic (upper left) and ZDC-SMD first harmonic (upper right) event-plane resolution as a function of $v^{obs}_{2}$ in 20-40% central Au+Au collisions. The TPC second-harmonic (lower left) and ZDC-SMD first harmonic (lower right) event-plane resolution for events with $|v^{obs}_{2}|<0.04$ as a function of centrality.

Jet-hadron correlations after background subtraction. Shown with Gaussian fits to jet peaks and systematic uncertanty bands Au+Au(4-6 GeV).

Gaussian widths of the awayside jet peaks in p+p(10-15 GeV).

Gaussian widths of the awayside jet peaks in Au+Au(10-15 GeV).

Gaussian widths of the awayside jet peaks in p+p(20-40 GeV).

Gaussian widths of the awayside jet peaks in Au+Au(20-40 GeV).

Awayside momentum difference(DAA) shown for 10-15 GeV.

Awayside momentum difference(DAA) shown for 10-15 GeV.

(Color online) Modulated sign correlations (msc) compared to the three-point correlator versus centrality for Au+Au collisions at $\sqrt{s_{NN}}$= 200 GeV. Shown with triangles is the msc, Eq. 5. The systematic uncertainties will be shown in detail in Fig. 7. Diamonds and circles show the three-point correlator, Eq. 1, and the grey bars reflect the conditions of $\Delta p_{T}$ > 0.15 GeV/c and $\Delta\eta$ > 0.15 applied to the three-point correlator, to be discussed in the text. For comparison, the model calculation of THERMINATOR [21] is also shown.

(Color online) The msc split into 2 composite parts versus centrality for Au+Au collisions at $\sqrt{s_{NN}}$= 200 GeV. Shaded areas represent the systematic uncertainty due to the event plane determination. For comparison with the $\Delta$msc term, we also put $−\nu_{2}/N$ and the model calculation of MEVSIM [24], to be described in the text.

(Color online) The msc split into 2 composite parts versus centrality for Au+Au collisions at $\sqrt{s_{NN}}$= 200 GeV. Shaded areas represent the systematic uncertainty due to the event plane determination. For comparison with the $\Delta$msc term, we also put $−\nu_{2}/N$ and the model calculation of MEVSIM [24], to be described in the text.

(Color online) Three-point correlations split up into out-of-plane $(\langle sin(\Delta\phi_{\alpha}) sin(\Delta\phi_{\beta})\rangle)$ and in-plane $(\langle cos(\Delta\phi_{\alpha}) cos(\Delta\phi_{\beta})\rangle)$ composite parts for 40 − 60% Au+Au collisions at $\sqrt{s_{NN}}$ GeV. (a)shows the correlations versus $\langle\eta\rangle$ = $(\eta_{\alpha}+\eta_{\beta})/2$. (b) shows the correlations versus $|\Delta\eta|$ = $|\eta_{\alpha}-\eta_{\beta}|$. Statistical errors are smaller than the symbol size. Systematic errors are given by the shaded bands and apply only to the difference of in-plane and out-of-plane parts.

(Color online) Three-point correlations split up into out-of-plane $(\langle sin(\Delta\phi_{\alpha}) sin(\Delta\phi_{\beta})\rangle)$ and in-plane $(\langle cos(\Delta\phi_{\alpha}) cos(\Delta\phi_{\beta})\rangle)$ composite parts for 40 − 60% Au+Au collisions at $\sqrt{s_{NN}}$ GeV. (a)shows the correlations versus $\langle\eta\rangle$ = $(\eta_{\alpha}+\eta_{\beta})/2$. (b) shows the correlations versus $|\Delta\eta|$ = $|\eta_{\alpha}-\eta_{\beta}|$. Statistical errors are smaller than the symbol size. Systematic errors are given by the shaded bands and apply only to the difference of in-plane and out-of-plane parts.

(Color online) Three-point correlations split up into out-of-plane and in-plane composite parts for 40 − 60% Au+Au collisions at $\sqrt{s_{NN}}$= 200 GeV. (a) shows the correlatios versus $\langle p_{T}\rangle$ = $(p_{T,\alpha}+p_{T,\beta})/2$. (b) shows the correlations versus $|\Delta p_{T}|$ = $|p_{T,\alpha}-p_{T,\beta}|$. Statistical errors are smaller than the symbol size. Systematic errors are given by the shaded bands and apply only to the difference of in-plane and out-of-plane parts.

(Color online) Three-point correlations split up into out-of-plane and in-plane composite parts for 40 − 60% Au+Au collisions at $\sqrt{s_{NN}}$= 200 GeV. (a) shows the correlatios versus $\langle p_{T}\rangle$ = $(p_{T,\alpha}+p_{T,\beta})/2$. (b) shows the correlations versus $|\Delta p_{T}|$ = $|p_{T,\alpha}-p_{T,\beta}|$. Statistical errors are smaller than the symbol size. Systematic errors are given by the shaded bands and apply only to the difference of in-plane and out-of-plane parts.

Transverse mass dependence of Gaussian radii Rout, Rside and Rlong and of experimental values of the Gaussian fit parameter lambda for midrapidity kaon pairs from the 30% most central Au+Au collisions at sqrt(sNN)=200 GeV. The PHENIX datapoints of the plot in the paper are from Phys. Rev. Lett. 103, 142301 (2009).

Amplitude (A) as a function of centrality for Like Sign (LS) and Charge Independent (CI) particles for Au+Au collision at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of pseudorapidity for 0-5% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of pseudorapidity for 5-10% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of pseudorapidity for 10-20% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of pseudorapidity for 20-30% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of pseudorapidity for 30-40% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of pseudorapidity for 40-60% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of transverse momentum for the wideGaussian method and for fhe event plane in the TPC for 0-5% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of transverse momentum for the wideGaussian method and for fhe event plane in the TPC for 30-40% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of transverse momentum for 0-10% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of transverse momentum for 10-20% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of transverse momentum for 20-30% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of transverse momentum for 30-40% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of transverse momentum for 40-50% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of transverse momentum for 50-60% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The third harmonic coefficient as a function of centrality from different methods of analysis for Au+Au collisions at sqrt(snn)=200 GEV.

The square of the third harmonic coefficient as a function of pseudorapidity separation (Deta) for Au+Au collisions at sqrt(snn)=200 GEV.

The fourth power of the third harmonic coefficient from four-particle cumulants is plotted as a function of centrality for Au+Au collisions at sqrt(snn)=200 GEV.

The ratio of the fourth power of the third harmonic coefficient from four-particle cumulants and the fourth power of the third harmonic flow is plotted as a function of centrality for Au+Au collisions at sqrt(snn)=200 GEV.

The v_2{TPC} as a function of transverse momentum for 0-5% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The v_2{TPC} as a function of transverse momentum for 20-30% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The v_2{TPC} as a function of transverse momentum for 30-40% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The v_3{TPC} as a function of transverse momentum for 0-5% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The v_3{TPC} as a function of transverse momentum for 20-30% centrality Au+Au collisions at sqrt(snn)=200 GEV.

The v_3{TPC} as a function of transverse momentum for 30-40% centrality Au+Au collisions at sqrt(snn)=200 GEV.