J/$\psi$ polarization parameters, measured for Pb--Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV, in the helicity and Collins-Soper reference frames in the rapidity interval $2.5<y<4$.The first uncertainty is statistical and the second systematic.

J/$\psi$ polarization parameters, measured for Pb--Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV, in the helicity and Collins-Soper reference frames in the rapidity interval $2.5<y<4$.The first uncertainty is statistical and the second systematic.

J/$\psi$ polarization parameters, measured for Pb--Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV, in the helicity and Collins-Soper reference frames in the rapidity interval $2.5<y<4$.The first uncertainty is statistical and the second systematic.

$\Upsilon(1{\rm S})$ polarization parameters in the helicity and Collins-Soper reference frames measured in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV in the rapidity interval $2.5<y<4$ and for transverse momentum $p_{\rm T}<15$ GeV/$c$. The first uncertainty is statistical and the second systematic.

$\Upsilon(1{\rm S})$ polarization parameters in the helicity and Collins-Soper reference frames measured in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV in the rapidity interval $2.5<y<4$ and for transverse momentum $p_{\rm T}<15$ GeV/$c$. The first uncertainty is statistical and the second systematic.

$\Upsilon(1{\rm S})$ polarization parameters in the helicity and Collins-Soper reference frames measured in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV in the rapidity interval $2.5<y<4$ and for transverse momentum $p_{\rm T}<15$ GeV/$c$. The first uncertainty is statistical and the second systematic.

$\Upsilon(1{\rm S})$ polarization parameters in the helicity and Collins-Soper reference frames measured in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV in the rapidity interval $2.5<y<4$ and for transverse momentum $p_{\rm T}<15$ GeV/$c$. The first uncertainty is statistical and the second systematic.

$\Upsilon(1{\rm S})$ polarization parameters in the helicity and Collins-Soper reference frames measured in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV in the rapidity interval $2.5<y<4$ and for transverse momentum $p_{\rm T}<15$ GeV/$c$. The first uncertainty is statistical and the second systematic.

$\Upsilon(1{\rm S})$ polarization parameters in the helicity and Collins-Soper reference frames measured in Pb--Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV in the rapidity interval $2.5<y<4$ and for transverse momentum $p_{\rm T}<15$ GeV/$c$. The first uncertainty is statistical and the second systematic.

Rapidity distribution of $\bar{p}$. Combined statitiscal uncertainty and systematic uncertainty from PID contramination. Systematic uncertainties from the track reconstruction efficiency are $\pm$25%.

Average transverse momentum of $p+\bar{p}$. Systematic uncertainties are $\pm$6%.

Average transverse momentum vs. centrality. Systematic error bands in figure are $\pm$6% for <pT>.

$dN/dy$ vs. centrality. Systematic error bands in figure are $\pm$25% for dN/dy.

$dN/dy$ vs. centrality. Systematic error bands in figure are $\pm$25% for dN/dy.

Average transverse momentum vs. the number of negative hadrons. Systematic error bands in figure are $\pm$6% for <pT>.

$dN/dy$ vs. the number of negative hadrons. Systematic error bands in figure are $\pm$25% for dN/dy.

$dN/dy$ vs. the number of negative hadrons. Systematic error bands in figure are $\pm$25% for dN/dy.

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.

Symmetrized pair-density net ratios $\widehat{r}[X(p_{t1}),X(p_{t2})]-1$ for all nonidentified charged primary particles for peripheral Au-Au collision events at $\sqrt{s_{NN}}=130$ GeV.

The values of $\Delta \langle p^{2}_{\rm T}\rangle$ as a function of centrality at backward ($-4.46 < y_{\rm cms} < -2.96$) and forward ($2.03 < y_{\rm cms} < 3.53$) centre-of-mass rapidity. The first uncertainty is statistical, the second and third ones are the systematic uncertainties. The third is related to the uncertainty on the $\langle p^{2}_{\rm T} \rangle$ in pp collisions and is correlated over centrality.

The nuclear modification factor ($Q_{\rm pPb}$) of inclusive J/$\psi$ as function of $\langle N_{\rm coll} \rangle$ at backward ($-4.46 < y_{\rm cms} < -2.96$) and forward ($2.03 < y_{\rm cms} < 3.53$) centre-of-mass rapidity. The first uncertainty is statistical, the second and third ones are the systematic uncertainties. The third uncertainty is fully correlated over centrality.

The $p_{\rm T}$-differential inclusive J/$\psi$ $Q_{\rm pPb}$ for six centrality classes at backward ($-4.46 < y_{\rm cms} < -2.96$) centre-of-mass rapidity. The first uncertainty is statistical, the second one and the third ones are the uncorrelated and correlated systematic uncertainties, respectively. The third uncertainty is fully correlated over $p_{\rm T}$. The fourth uncertainty is the correlated systematic which is common for all centrality classes. In Figure 7 correlated systematic uncertainties are added in quadrature.

The $p_{\rm T}$-differential inclusive J/$\psi$ $Q_{\rm pPb}$ for six centrality classes at forward ($2.03 < y_{\rm cms} < 3.53$) centre-of-mass rapidity. The first uncertainty is statistical, the second one and the third ones are the uncorrelated and correlated systematic uncertainties, respectively. The third uncertainty is fully correlated over $p_{\rm T}$. The fourth uncertainty is the correlated systematic which is common for all centrality classes. In Figure 8 correlated systematic uncertainties are added in quadrature.

The ratio of inclusive J/$\psi$ $Q_{\rm pPb}$ between peripheral to central collisions ($Q_{\rm PC}$) as function of $p_{\rm T}$ for six centrality classes at backward ($-4.46 < y_{\rm cms} < -2.96$) centre-of-mass rapidity. The first uncertainty is statistical, the second one and the third ones are the systematic uncertainties. The third uncertainty is fully correlated over $p_{\rm T}$.

The ratio of inclusive J/$\psi$ $Q_{\rm pPb}$ between peripheral to central collisions ($Q_{\rm PC}$) as a function of $p_{\rm T}$ for six centrality classes at forward ($2.03 < y_{\rm cms} < 3.53$) centre-of-mass rapidity. The first uncertainty is statistical, the second and the third ones are the systematic uncertainties. The third uncertainty is fully correlated over $p_{\rm T}$.

B.R.$_{\psi(\rm 2S)\rightarrow\mu^{+}\mu^{-}}{\sigma_{\psi(\rm 2S)}}$/B.R.$_{J/\psi\rightarrow\mu^{+}\mu^{-}}{\sigma_{\rm J/\psi}}$ as a function of $\langle N_{\rm coll} \rangle$ at backward ($-4.46 < y_{\rm cms} < -2.96$) centre-of-mass rapidity. The first uncertainty is statistical, the second and the third ones are the systematic uncertainties. The third uncertainty is fully correlated over centrality.

B.R.$_{\psi(\rm 2S)\rightarrow\mu^{+}\mu^{-}}{\sigma_{\psi(\rm 2S)}}$/B.R.$_{J/\psi\rightarrow\mu^{+}\mu^{-}}{\sigma_{\rm J/\psi}}$ as a function of $\langle N_{\rm coll} \rangle$ at forward ($2.03 < y_{\rm cms} < 3.53$) centre-of-mass rapidity. The first uncertainty is statistical, the second and the third ones are the systematic uncertainties. The third uncertainty is fully correlated over centrality.

Double ratio $[\sigma_{\psi(\rm 2S)}/\sigma_{\rm J/\psi}]_{\rm pPb}/[\sigma_{\psi(\rm 2S)}/\sigma_{\rm J/\psi}]_{\rm pp}$ as a function of $\langle N_{\rm coll} \rangle$ at backward ($-4.46 < y_{\rm cms} < -2.96$) centre-of-mass rapidity. The first uncertainty is statistical, the second and the third ones are the systematic uncertainties. The third uncertainty is fully correlated over centrality.

Double ratio $[\sigma_{\psi(\rm 2S)}/\sigma_{\rm J/\psi}]_{\rm pPb}/[\sigma_{\psi(\rm 2S)}/\sigma_{\rm J/\psi}]_{\rm pp}$ as a function of $\langle N_{\rm coll} \rangle$ at forward ($2.03 < y_{\rm cms} < 3.53$) centre-of-mass rapidity. The first uncertainty is statistical, the second and the third ones are the systematic uncertainties. The third uncertainty is fully correlated over centrality.

Inclusive $\psi$(2S) $Q_{\rm pPb}$ as a function of $\langle N_{\rm coll} \rangle$ at backward ($-4.46 < y_{\rm cms} < -2.96$) centre-of-mass rapidity. The first uncertainty is statistical, the second is the uncorrelated systematic, the third is the correlated systematic for $\psi$(2S) or J/$\psi$ only while the fourth one is the correlated systematic shared among $\psi$(2S) and J/$\psi$.

Inclusive $\psi$(2S) $Q_{\rm pPb}$ as a function of $\langle N_{\rm coll} \rangle$ at forward ($2.03 < y_{\rm cms} < 3.53$) centre-of-mass rapidity. The first uncertainty is statistical, the second is the uncorrelated systematic, the third is the correlated systematic for $\psi$(2S) or J/$\psi$ only while the fourth one is the correlated systematic shared among $\psi$(2S) and J/$\psi$.

Inclusive J/$\psi$ cross sections extrapolated from existing measurements at backward ($-4.46 < y < -2.96$) and forward ($2.03 < y < 3.53$) rapidity in pp collisions. The first and the second uncertainties are the uncorrelated and correlated systematic uncertainties, respectively. The correlated uncertainty is correlated over rapidity. This table refers to Fig. 1 of ALICE-PUBLIC-2020-007.

Inclusive $\psi$(2S) cross sections extrapolated from existing measurements at backward ($-4.46 < y < -2.96$) and forward ($2.03 < y < 3.53$) rapidity in pp collisions . The first and the second uncertainties are the uncorrelated and correlated systematic uncertainties, respectively. The correlated uncertainty is fully correlated over rapidity. This table refers Fig. 1 of ALICE-PUBLIC-2020-007.

The rapidity differential inclusive J/$\psi$ cross sections extrapolated from existing measurements for backward ($-4.46 < y < -2.96$) rapidity interval in pp collisions. The first and the second uncertainties are the uncorrelated and correlated systematic uncertainties, respectively. The correlated uncertainty is correlated over rapidity. This table refers Fig. 2 of ALICE-PUBLIC-2020-007.

The rapidity differential inclusive J/$\psi$ cross sections extrapolated from existing measurements for forward ($2.03 < y < 3.53$) rapidity interval in pp collisions. The first and the second uncertainties are the uncorrelated and correlated systematic uncertainties, respectively. The correlated uncertainty is correlated over rapidity. This table refers Fig. 2 of ALICE-PUBLIC-2020-007.

The $p_{\rm T}$-differential inclusive J/$\psi$ cross sections extrapolated from existing measurements at backward ($-4.46 < y < -2.96$) rapidity in pp collisions. The first and the second uncertainties are the uncorrelated and correlated systematic uncertainties, respectively. The correlated uncertainty is fully correlated over $p_{\rm T}$. This table refers Fig. 3 of ALICE-PUBLIC-2020-007.

The $p_{\rm T}$-differential inclusive J/$\psi$ cross sections extrapolated from existing measurements at forward ($2.03 < y < 3.53$) rapidity in pp collisions. The first and the second uncertainties are the uncorrelated and correlated systematic uncertainties, respectively. The correlated uncertainty is fully correlated over $p_{\rm T}$. This table refers Fig. 3 of ALICE-PUBLIC-2020-007.

The values of $\langle p_{\rm T} \rangle$ and $\langle p^{2}_{\rm T} \rangle$ of inclusive J/$\psi$ at backward ($-4.46 < y < -2.96$) and forward ($2.03 < y < 3.53$) rapidity calculated from extrapolated cross sections in pp collisions. The uncertainty is the systematic uncertainty. The systematic uncertainty is obtained as the quadratic sum of the uncorrelated and correlated systematic uncertainties. This table refers Fig. 4 of ALICE-PUBLIC-2020-007.

Transverse momentum spectra of PI− in rapidity slices produced in the 5% most central Be+Be collisions at 75A GeV/c.

Transverse momentum spectra of PI− in rapidity slices produced in the 5% most central Be+Be collisions at 150A GeV/c.

Transverse mass spectra 1/mT d2n/(dydmT) of negatively charged pions produced in central Be+Be collisions at the SPS energies.

Rapidity spectra of PI− produced in the 5% most central Be+Be collisions at the SPS energies.

Collision energy dependence of the ratio of the width sigma_y of the rapidity distribution to the beam rapidity y_beam for central Be+Be.

Energy dependence of the mean transverse mass of PI- measured at mid-rapidity in central Be+Be collisions.

The "kink" plot showing the ratio of pion multiplicity <PI> to number of wounded nucleons <W> versus the Fermi energy variable F=s^1/4_NN from central Be+Be collisions.

The $p_{T}$ dependance of $v_{4}^{Inclusive}$, $v_{4}^{Linear}$, $v_{4}^{Mode-coupled}$, $\chi_{4,22}$ and $\rho_{4,22}$

The $p_{T}$ dependance of $v_{5}^{Inclusive}$, $v_{5}^{Linear}$, $v_{5}^{Mode-coupled}$, $\chi_{5,23}$ and $\rho_{5,23}$

Fig. 2. (Color online.) Event-by-event photon multiplicity distributions (solid circles) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV} .$ The distributions for top $0-5 \%$ central $\mathrm{Au}+$ Au collisions and top $0-10 \%$ central $\mathrm{Cu}+\mathrm{Cu}$ collisions are also shown (open circles). The photon multiplicity distributions for central collisions are observed to be Gaussian (solid line). Only statistical errors are shown. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 2. (Color online.) Event-by-event photon multiplicity distributions (solid circles) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV} .$ The distributions for top $0-5 \%$ central $\mathrm{Au}+$ Au collisions and top $0-10 \%$ central $\mathrm{Cu}+\mathrm{Cu}$ collisions are also shown (open circles). The photon multiplicity distributions for central collisions are observed to be Gaussian (solid line). Only statistical errors are shown. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 2. (Color online.) Event-by-event photon multiplicity distributions (solid circles) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV} .$ The distributions for top $0-5 \%$ central $\mathrm{Au}+$ Au collisions and top $0-10 \%$ central $\mathrm{Cu}+\mathrm{Cu}$ collisions are also shown (open circles). The photon multiplicity distributions for central collisions are observed to be Gaussian (solid line). Only statistical errors are shown. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several centrality classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several centrality classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several centrality classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 3. (Color online.) Photon pseudorapidity distributions for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s}_{\mathrm{NN}}}=62.4$ and $200\mathrm{GeV}$. The results for several classes are shown. The solid curves are results of HIJING simulations for central $(0-5 \%$ for $\mathrm{Au}+\mathrm{Au}$ and $0-10 \%$ for $\mathrm{Cu}+\mathrm{Cu})$ and $30-40 \%$ mid-central collisions. The errors shown are systematic, statistical errors are negligible in comparison. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 4. (Color online.) Top panel: The number of photons divided by $\left\langle N_{\text{part}}\right\rangle / 2$ as a function of average number of participating nucleons for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s} _{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$ for $-3.7<\eta<-2.3 .$ Errors shown are systematic only and include those for $\left\langle N_{\text {part }}\right\rangle .$ Results from HIJING are shown as lines (solid for $\mathrm{Au}+\mathrm{Au}$ and dashed for $\mathrm{Cu}+\mathrm{Cu}$ ). Bottom panel: Same as above, for both photons and charged particles from PHOBOS [10]. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 4. (Color online.) Top panel: The number of photons divided by $\left\langle N_{\text{part}}\right\rangle / 2$ as a function of average number of participating nucleons for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ at $\sqrt{\mathrm{s} _{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$ for $-3.7<\eta<-2.3 .$ Errors shown are systematic only and include those for $\left\langle N_{\text {part }}\right\rangle .$ Results from HIJING are shown as lines (solid for $\mathrm{Au}+\mathrm{Au}$ and dashed for $\mathrm{Cu}+\mathrm{Cu}$ ). Bottom panel: Same as above, for both photons and charged particles from PHOBOS [10]. NOTE: For points with invisible error bars, the point size was considered as an absolute upper limit for the uncertainty.

Fig. 5. (Color online.) Photon pseudorapidity distributions normalized by the average number of participating nucleon pairs for different collision centralities are plotted as a function of pseudorapidity shifted by the beam rapidity (-5.36 for $200 \mathrm{GeV}$ and -4.19 for $62.4 \mathrm{GeV}$ ) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ collisions at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$. Errors are systematic only, statistical errors are negligible in comparison. For clarity of presentation, results for only four centralities are shown. The $\mathrm{Cu}+$ Cu data are shifted by 0.1 unit in $\eta-y_{\text {beam }} .$ The solid line is a second order polynomial fit to the data (see text for details).

Fig. 5. (Color online.) Photon pseudorapidity distributions normalized by the average number of participating nucleon pairs for different collision centralities are plotted as a function of pseudorapidity shifted by the beam rapidity (-5.36 for $200 \mathrm{GeV}$ and -4.19 for $62.4 \mathrm{GeV}$ ) for $\mathrm{Au}+\mathrm{Au}$ and $\mathrm{Cu}+\mathrm{Cu}$ collisions at $\sqrt{s_{\mathrm{NN}}}=62.4$ and $200 \mathrm{GeV}$. Errors are systematic only, statistical errors are negligible in comparison. For clarity of presentation, results for only four centralities are shown. The $\mathrm{Cu}+$ Cu data are shifted by 0.1 unit in $\eta-y_{\text {beam }} .$ The solid line is a second order polynomial fit to the data (see text for details).