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.

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{lead.~lep.}}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.~lep.}}$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production. Overflow events are included in the last bin. The largest observed value is 609 GeV.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.~lep.}}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.~lep.}}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $p_{\mathrm{T}}^{\mathrm{lead.~jet}}$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production. Overflow events are included in the last bin. The largest observed value is 1485 GeV.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{lead.~jet}}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{lead.~jet}}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $H_{\mathrm{T}}$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production. Overflow events are included in the last bin. The largest observed value is 2969 GeV.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $H_{\mathrm{T}}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $H_{\mathrm{T}}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $S_{\mathrm{T}}$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production. Overflow events are included in the last bin. The largest observed value is 3296 GeV.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $S_{\mathrm{T}}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $S_{\mathrm{T}}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $m_{\mathrm{T},e\mu}$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production. Overflow events are included in the last bin. The largest observed value is 4130 GeV.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $m_{\mathrm{T},e\mu}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $m_{\mathrm{T},e\mu}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $m_{e\mu}$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production. Overflow events are included in the last bin. The largest observed value is 3519 GeV.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $m_{e\mu}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $m_{e\mu}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $p_{\mathrm{T},e\mu}$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production. Overflow events are included in the last bin. The largest observed value is 1067 GeV.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T},e\mu}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T},e\mu}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $\Delta\phi(e,\mu)$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $\Delta\phi(e,\mu)$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $\Delta\phi(e,\mu)$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $y_{e\mu}$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $y_{e\mu}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $y_{e\mu}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $\cos\theta^*$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $\cos\theta^*$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $\cos\theta^*$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $n_{\mathrm{jet}}$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $n_{\mathrm{jet}}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $n_{\mathrm{jet}}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $m_{e\mu}$ for $p_{\mathrm{T}}^{\mathrm{lead.~jet}} > 200$ GeV. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production. Overflow events are included in the last bin. The largest observed value is 3519 GeV.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $m_{e\mu}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $m_{e\mu}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $\Delta\phi(e,\mu)$ for $p_{\mathrm{T}}^{\mathrm{lead.~jet}} > 200$ GeV. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $\Delta\phi(e,\mu)$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $\Delta\phi(e,\mu)$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $\Delta\phi(\mathrm{sub-lead.~lep.}, \mathrm{lead.~jet})$ for $p_{\mathrm{T}}^{\mathrm{lead.~lep.}} > 200$ GeV. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $\Delta\phi(\mathrm{sub-lead.~lep.}, \mathrm{lead.~jet})$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $\Delta\phi(\mathrm{sub-lead.~lep.}, \mathrm{lead.~jet})$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $\Delta R(\mathrm{sub-lead.~lep.}, \mathrm{lead.~jet})$ for $p_{\mathrm{T}}^{\mathrm{lead.~lep.}} > 200$ GeV. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $\Delta R(\mathrm{sub-lead.~lep.}, \mathrm{lead.~jet})$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $\Delta R(\mathrm{sub-lead.~lep.}, \mathrm{lead.~jet})$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.~lep.}} / p_{\mathrm{T}}^{\mathrm{lead.~lep.}}$ for $p_{\mathrm{T}}^{\mathrm{lead.~lep.}} > 200$ GeV. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.~lep.}} / p_{\mathrm{T}}^{\mathrm{lead.~lep.}}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.~lep.}} / p_{\mathrm{T}}^{\mathrm{lead.~lep.}}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.~lep.}} / p_{\mathrm{T}}^{\mathrm{lead.~jet}}$ for $p_{\mathrm{T}}^{\mathrm{lead.~lep.}} > 200$ GeV. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production. The largest observed value is 19.6.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.~lep.}} / p_{\mathrm{T}}^{\mathrm{lead.~jet}}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.~lep.}} / p_{\mathrm{T}}^{\mathrm{lead.~jet}}$

Distribution of the four-lepton invariant mass in the four-lepton final state for the ggF-MVA-low category.

Distribution of the four-lepton invariant mass in the four-lepton final state for the VBF-MVA-enriched category.

Distribution of the transverse mass in the llvv final state for the ggF-enriched eevv category.

Distribution of the transverse mass in the llvv final state for the ggF-enriched mumuvv category.

Distribution of the transverse mass in the llvv final state for the VBF-enriched eevv category.

Distribution of the transverse mass in the llvv final state for the VBF-enriched mumuvv category.

The upper limits at 95% CL on the cross section times branching ratio as a function of the heavy resonance mass for the ggF production mode

The upper limits at 95% CL on the cross section times branching ratio as a function of the heavy resonance mass for the VBF production mode

The upper limits at 95% CL on the cross section for the ggF production mode times branching ratio as a function of the heavy resonance mass for an additional heavy scalar assuming a width of 1% of its mass

The upper limits at 95% CL on the cross section for the ggF production mode times branching ratio as a function of the heavy resonance mass for an additional heavy scalar assuming a width of 5% of its mass

The upper limits at 95% CL on the cross section for the ggF production mode times branching ratio as a function of the heavy resonance mass for an additional heavy scalar assuming a width of 10% of its mass

The upper limits at 95% CL on the cross section for the ggF production mode times branching ratio as a function of the heavy resonance mass for an additional heavy scalar assuming a width of 15% of its mass

The upper limits at 95% CL on the cross section times branching ratio for a KK graviton produced with k/M_{PI} = 1

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.

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

Slope parameter T, dN/y, $N_\phi/N_h^-$. Systematic uncertainties on T are 12%, on dN/dy are 12%, and on $N_{\phi}/N_{h^-}$ are 22%.