The combined differential $D^{*\pm}$-production cross section as a function of $Q^2$, with its uncorrelated and correlated uncertainties.

The combined differential $D^{*\pm}$-production cross section as a function of $y$, with its uncorrelated and correlated uncertainties.

The combined double-differential $D^{*\pm}$-production cross section as a function of $Q^2$ and $y$, with its uncorrelated and correlated uncertainties.

The effective exponential slope, $b_n$, obtained from the fit of double differential photoproduction cross sections ${\rm d^2}\sigma_{\gamma p}/{\rm d}x_L{\rm d}p_{T,n}^2$ to a single exponential function in bins of $x_L$. The first uncertainty represents the fit error from the statistical and uncorrelated systematic uncertainty and the second one is due to the correlated systematic uncertainty.

Energy dependence of the exclusive photoproduction of a $\rho^0$ meson associated with a leading neutron, $\gamma p \to \rho^0 n \pi^+$. The first uncertainty is statistical and the second is systematic. The global normalisation uncertainty of $4.4\%$ is not included. $\Phi_{\gamma}$ is the integral of the photon flux, Eq. (3) of paper, in a given $W_{\gamma p}$ bin.

Differential photoproduction cross section ${\rm d}\sigma_{\gamma p}/{\rm d}\eta$ for the exclusive process $\gamma p \to \rho^0 n \pi^+$ as a function of the $\rho^0$ pseudorapidity in the kinematic range $0.35 < x_L < 0.95$, $\theta_n < 0.75$ mrad and $20 < W_{\gamma p} < 100$ GeV. The statistical, uncorrelated and correlated systematic uncertainties, $\delta_{stat}$, $\delta_{sys}^{unc}$ and $\delta_{sys}^{cor}$ respectively, are given, which does not include the global normalisation error of $4.4\%$.

Differential photoproduction cross section ${\rm d}\sigma_{\gamma p}/{\rm d}t^{\prime}$ for the exclusive process $\gamma p \to \rho^0 n \pi^+$ as a function of the $\rho^0$ four-momentum transfer squared, $t^{\prime}$, in the kinematic range $0.35 < x_L < 0.95$, $\theta_n < 0.75$ mrad and $20 < W_{\gamma p} < 100$ GeV. The statistical, uncorrelated and correlated systematic uncertainties, $\delta_{stat}$, $\delta_{sys}^{unc}$ and $\delta_{sys}^{cor}$ respectively, are given, which does not include the global normalisation error of $4.4\%$.

Exponential slopes, $b_1$ and $b_2$, and the ratio $\sigma_1/\sigma_2$, obtained from the components of fit, Eq. (14) of paper, to the differential cross section ${\rm d}\sigma_{\gamma p}/{\rm d}t^{\prime}$ in bins of $x_L$. The errors represent the statistical and systematic uncertainties added in quadrature.

Exponential slopes, $b_1$ and $b_2$, and the ratio $\sigma_1/\sigma_2$, obtained from the components of fit, Eq. (14) of paper, to the differential cross section ${\rm d}\sigma_{\gamma p}/{\rm d}t^{\prime}$ in bins of $p_{T,n}^2$. The errors represent the statistical and systematic uncertainties added in quadrature.

Energy dependence of elastic $\rho^0$ photoproduction on the pion, $\gamma \pi^+ \to \rho^0 \pi^+$, extracted in the one-pion-exchange approximation using OPE1 sample. The first uncertainty represents the full experimental error and the second is the model error coming from the pion flux uncertainty (see text). $\Gamma_\pi(x_L)$ represents the value of the pion flux, Eqs. (5-6) of paper, integrated over the $p_{T,n}<0.2$ GeV range, at a given $x_L$.

Cross section of elastic $\rho^0$ photoproduction on the pion, $\gamma \pi^+ \to \rho^0 \pi^+$, extracted in the one-pion-exchange approximation using three different samples: full sample, OPE1 and OPE2. The first uncertainty represents the full experimental error and the second is the model error coming from the pion flux uncertainty (see text). $\Gamma_\pi$ represents the value of the pion flux, Eqs. (5-6) of paper, integrated over the corresponding $(x_L,p_{T,n})$ range.

The dihadron correlated yield normalized per radian per unit of pseudorapidity as function of $\Delta\eta$ in d+Au collisions on the away side (|$\Delta\phi$ - $\pi$| < $\pi$/3). Shown is the high FTPC-Au activity data. Trigger and associated particles have 1 < $p_T$ < 3 GeV/c and |$\eta$| < 1.

The dihadron correlated yield normalized per radian per unit of pseudorapidity as function of $\Delta\eta$ in d+Au collisions on the away side (|$\Delta\phi$ - $\pi$| < $\pi$/3). Shown is the high-activity data after subtracting the unscaled low-activity data. Trigger and associated particles have 1 < $p_T$ < 3 GeV/c and |$\eta$| < 1.

The dihadron correlated yield normalized per radian per unit of pseudorapidity as function of $\Delta\eta$ in d+Au collisions on the away side (|$\Delta\phi$ - $\pi$| < $\pi$/3). Shown is the high-activity data after subtracting the unscaled low-activity data. Trigger and associated particles have 1 < $p_T$ < 3 GeV/c and |$\eta$| < 1.

The dihadron correlated yield normalized per radian per unit of pseudorapidity as function of $\Delta\eta$ in d+Au collisions on the near (|$\Delta\phi$| < $\pi$/3) side. Shown is the high-activity data after subtracting the unscaled low-activity data. Trigger and associated particles have 1 < $p_T$ < 3 GeV/c and |$\eta$| < 1.

The dihadron correlated yield normalized per radian per unit of pseudorapidity as function of $\Delta\eta$ in d+Au collisions on the near (|$\Delta\phi$| < $\pi$/3, solid circles) and away side (|$\Delta\phi$ - $\pi$| < $\pi$/3, open circles). Shown are the (a) low and (b) high FTPC-Au activity data, and the high-activity data after subtracting the (c) unscaled. Trigger and associated particles have 1 < $p_T$ < 3 GeV/c and |$\eta$| < 1.

(a) The near-side jetlike correlated yield obtained from Gaussian fit as in Fig. 1 as function of the uncorrected dN/d$\eta$ at midrapidity measured in the TPC. Two event selections are used: FTPC-Au multiplicity and ZDC-Au energy. The curve is the result from a HIJING calculation.

(a) The near-side jetlike correlated yield obtained from Gaussian fit as in Fig. 1 as function of the uncorrected dN/d$\eta$ at midrapidity measured in the TPC. Two event selections are used: FTPC-Au multiplicity and ZDC-Au energy. The curve is the result from a HIJING calculation.

$p_T^{(a)}$ dependence. The ratio of the correlated yields in high over low FTPC-Au multiplicity events as a function of $p_T^{(a)}$ $p_T^{(t)}$ where $p_T^{(a)}$ $p_T^{(t)}$ is fixed.

$p_T^{(a)}$ dependence. The ratio of the correlated yields in high over low FTPC-Au multiplicity events as a function of $p_T^{(a)}$ $p_T^{(t)}$ where $p_T^{(a)}$ $p_T^{(t)}$ is fixed.

(a) The dihadron correlated yield normalized per radian per unit of pseudorapidity as a function of $\Delta\phi$ in d+Au collisions at high (0-20%) FTPC-Au multiplicities. Trigger and associated particles are 1 < $p_T$ < 3 GeV/c within 0.5 < |$\Delta\eta$| < 0.7. ZYAM positions are indicated with arrows.

(a) The dihadron correlated yield normalized per radian per unit of pseudorapidity as a function of $\Delta\phi$ in d+Au collisions at low (40-100%) FTPC-Au multiplicities. Trigger and associated particles are 1 < $p_T$ < 3 GeV/c within 0.5 < |$\Delta\eta$| < 0.7. ZYAM positions are indicated with arrows.multiplicity versus $\Delta\phi$.

The raw associated yield at high FTPC-Au multiplicity minus the unscaled ZYAM-subtracted correlated yields at low FTPC-Au multiplicity versus $\Delta\phi$.

The raw associated yield at high FTPC-Au multiplicity minus the scaled ZYAM-subtracted correlated yields at low FTPC-Au multiplicity versus $\Delta\phi$.

Bin averaged hadron level differential cross sections as a function of the variable $z_{I\!\!P}$ for diffractive dijet DIS. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations and the radiative corrections ($1+\delta_{\text{rad}}$) are also listed. The overall normalisation uncertainty of $6\%$ is not included in the table.

Bin averaged hadron level differential cross sections as a function of the variable $x_{I\!\!P}$ for diffractive dijet DIS. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations and the radiative corrections ($1+\delta_{\text{rad}}$) are also listed. The overall normalisation uncertainty of $6\%$ is not included in the table.

Bin averaged hadron level differential cross sections as a function of the variable $y$ for diffractive dijet DIS. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations and the radiative corrections ($1+\delta_{\text{rad}}$) are also listed. The overall normalisation uncertainty of $6\%$ is not included in the table.

Bin averaged hadron level differential cross sections as a function of the variable $Q^2$ for diffractive dijet DIS. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations and the radiative corrections ($1+\delta_{\text{rad}}$) are also listed. The overall normalisation uncertainty of $6\%$ is not included in the table.

Bin averaged hadron level differential cross sections for diffractive dijet DIS as a function of the variable $E_T^{*\text{jet1}}$. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations and the the radiative corrections ($1+\delta_{\text{rad}}$) are also listed. The overall normalisation uncertainty of $6\%$ is not included in the table.

Bin averaged hadron level differential cross sections for diffractive dijet DIS as a function of the variable $M_X$. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations and the the radiative corrections ($1+\delta_{\text{rad}}$) are also listed. The overall normalisation uncertainty of $6\%$ is not included in the table.

Bin averaged hadron level differential cross sections for diffractive dijet DIS as a function of the variable $\vert\Delta\eta^{\text{jets}}\vert$. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations and the the radiative corrections ($1+\delta_{\text{rad}}$) are also listed. The overall normalisation uncertainty of $6\%$ is not included in the table.

Bin averaged hadron level differential cross sections for diffractive dijet DIS as a function of the variable $\langle\eta^{\text{jets}}\rangle$. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations and the the radiative corrections ($1+\delta_{\text{rad}}$) are also listed. The overall normalisation uncertainty of $6\%$ is not included in the table.

Bin averaged hadron level differential diffractive dijet $ep$ cross sections as a function of the variable $z_{I\!\!P}$ for the dijet photoproduction kinematic range. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations are given. The overall normalisation uncertainty of $6\%$ is not included in the table.

Bin averaged hadron level differential diffractive dijet $ep$ cross sections as a function of the variable $x_{I\!\!P}$ for the dijet photoproduction kinematic range. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations are given. The overall normalisation uncertainty of $6\%$ is not included in the table.

Bin averaged hadron level differential diffractive dijet $ep$ cross sections as a function of the variable $y$ for the dijet photoproduction kinematic range. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations are given. The overall normalisation uncertainty of $6\%$ is not included in the table.

Bin averaged hadron level differential diffractive dijet $ep$ cross sections as a function of the variable $x_\gamma$ for the dijet photoproduction kinematic range. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations are given. The overall normalisation uncertainty of $6\%$ is not included in the table.

Bin averaged hadron level differential diffractive dijet $ep$ cross sections in the photoproduction kinematic range as a function of the variable $E_T^{*\text{jet1}}$. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations are given. The overall normalisation uncertainty of $6\%$ is not included in the table.

Bin averaged hadron level differential diffractive dijet $ep$ cross sections in the photoproduction kinematic range as a function of the variable $M_X$. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations are given. The overall normalisation uncertainty of $6\%$ is not included in the table.

Bin averaged hadron level differential diffractive dijet $ep$ cross sections in the photoproduction kinematic range as a function of the variable $\vert\Delta\eta^{\text{jets}}\vert$. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations are given. The overall normalisation uncertainty of $6\%$ is not included in the table.

Bin averaged hadron level differential diffractive dijet $ep$ cross sections in the photoproduction kinematic range as a function of the variable $\langle\eta^{\text{jets}}\rangle$. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations are given. The overall normalisation uncertainty of $6\%$ is not included in the table.

Ratios of differential diffractive dijet $ep$ cross sections, measured in photoproduction, to measurements in DIS as a function of the variable $\vert\Delta\eta^{\text{jets}}\vert$. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations are given.

Ratios of differential diffractive dijet $ep$ cross sections, measured in photoproduction, to measurements in DIS as a function of the variable $y$. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations are given.

Ratios of differential diffractive dijet $ep$ cross sections, measured in photoproduction, to measurements in DIS as a function of the variable $z_{I\!\!P}$. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations are given.

Ratios of differential diffractive dijet $ep$ cross sections, measured in photoproduction, to measurements in DIS as a function of the variable $E_T^{*\text{jet1}}$. The hadronisation correction factors ($1+\delta_{\text{hadr}}$) applied to the NLO calculations are given.

The measured differential cross sections $\rm{d}\sigma/\rm{d}|\eta|$ of $D^{*\pm}$ and $D^\pm$ production with $20<p_T<100\,$GeV. The first and second errors are the statistical and systematic uncertainties, respectively. The systematic uncertainty fractions corresponding to the tracking ($\delta_2$) uncertainties (Table 2 of the paper) are strongly correlated. The fully correlated uncertainties linked with the luminosity measurement ($3.5\%$) and branching fractions ($1.5\%$ and $2.1\%$ for $D^{*\pm}$ and $D^\pm$, respectively) are not shown.

The jet-like yield in $\mid\Delta\eta\mid<$0.78 as a function of $p_T^{trigger}$ for $K_S^0$-h and $\Lambda$-h correlations for 1.5 GeV/$c$ < $p_T^{associated}$ < $p_T^{trigger}$ in minimum bias $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Statistical errors listed in table. 5.4% systematic error on all points.

The jet-like yield in $\mid\Delta\eta\mid<$0.78 as a function of $p_T^{trigger}$ for $K_S^0$-h and $\Lambda$-h correlations for 1.5 GeV/$c$ < $p_T^{associated}$ < $p_T^{trigger}$ in 0-60% Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV. Statistical errors listed in table. 5.4% systematic error on all points.

The jet-like yield in $\mid\Delta\eta\mid<$0.78 as a function of $p_T^{trigger}$ for $K_S^0$-h and $\Lambda$-h correlations for 1.5 GeV/$c$ < $p_T^{associated}$ < $p_T^{trigger}$ in 40-80% Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Statistical errors listed in table. 5.4% systematic error on all points.

The jet-like yield in $\mid\Delta\eta\mid<$0.78 as a function of $p_T^{trigger}$ for $K_S^0$-h and $\Lambda$-h correlations for 1.5 GeV/$c$ < $p_T^{associated}$ < $p_T^{trigger}$ in 0-12% Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. Statistical errors listed in table. 5.4% systematic error on all points.

Centrality dependence of the jet-like yield of $K_S^0$-h and $\Lambda$-h correlations for 3 < $p_T^{trigger}$ < 6 GeV/$c$ and 1.5 GeV/$c$ < $p_T^{associated}$ < $p_T^{trigger}$ in Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV. 5.4% systematic error on all points.

Centrality dependence of the jet-like yield of $K_S^0$-h and $\Lambda$-h correlations for 3 < $p_T^{trigger}$ < 6 GeV/$c$ and 1.5 GeV/$c$ < $p_T^{associated}$ < $p_T^{trigger}$ in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. 5.4% systematic error on all points.

Centrality dependence of the jet-like yield of $K_S^0$-h and $\Lambda$-h correlations for 3 < $p_T^{trigger}$ < 6 GeV/$c$ and 1.5 GeV/$c$ < $p_T^{associated}$ < $p_T^{trigger}$ in $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. 5.4% systematic error on all points.

Centrality dependence of the jet-like yield of $K_S^0$-h and $\Lambda$-h correlations for 3 < $p_T^{trigger}$ < 6 GeV/$c$ and 1.5 GeV/$c$ < $p_T^{associated}$ < $p_T^{trigger}$ in minimum bias $d$+Au and 0-60% central Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV. 5.4% systematic error on all points.

Correlation matrix of systematic uncertainties for $\sigma^{\pm}$. The statistical and integrated luminosity uncertainties are not included. The first and the last eleven indices correspond to the systematic uncertainites of $W^+$ ($\sigma^+_{1-11}$) and $W^-$ ($\sigma^-_{1-11}$) cross sections, respectively, in 11 muon pseudorapidity bins. The values are expressed as percentages.

Correlation matrix of statistical uncertainties for $\sigma^{\pm}$. The first and the last eleven indices correspond to the statistical uncertainites of $W^+$ ($\sigma^+_{1-11}$) and $W^-$ ($\sigma^-_{1-11}$) cross sections, respectively, in 11 muon pseudorapidity bins. The values are expressed as percentages.

Correlation matrix of systematic uncertainties for $\mathcal{A}$. The values are expressed as percentages.

$R_{pPb}$ as a function of $p_{T}$ integrated over rapidity range $-1.8<y*<1.3$ for the 0--90% centrality interval for the three geometric models: (a) Glauber, (b) Glauber-Gribov with $\omega_{\sigma}=0.11$ and (c) Glauber-Gribov with $\omega_{\sigma}=0.2$.

$R_{pPb}$ values as a function of $p_{T}$ in the 0--90% centrality interval averaged over $|y*|<0.5$.The $\langle T_{Pb} \rangle$ value for the ATLAS centrality correction is calculated with the Glauber model. The total systematic uncertainties, which include the uncertainty in $\langle T_{Pb} \rangle$, are indicated by lines of the same colour. Strict quantitative agreement is not expected as each measurement uses different rapidity intervals for the centrality determination and apply different event selection criteria to reject diffractive collisions.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $-1.8<y*<1.3$ for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $-1.8<y*<1.3$ for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $-1.8<y*<1.3$ for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{CP}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $|\eta|<2.3$ for seven centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{CP}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $|\eta|<2.3$ for seven centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{CP}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $|\eta|<2.3$ for seven centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

Sample correlations in $\Delta\phi$ ($|\Delta\eta|<$ 1.78) for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-95% $d$+Au at $\sqrt{s_{NN}}$ = 200 GeV, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV, 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV, and 0-12% central Au+Au at $\sqrt{s_{NN}}$ = 200 GeV. The data are averaged between positive and negative $\Delta\phi$ and reflected in the plot.

Background subtracted sample correlations for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ on the near-side for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-95% $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV, 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV and 0-12% central Au+Au at $\sqrt{s_{NN}}$ = 200 GeV. The dependence of the jet-like correlation is shown as a function of $\Delta\eta$ ($|\Delta\phi|<$ 0.78). The data are averaged between positive and negative $\Delta\eta$.

Background subtracted sample correlations for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ on the near-side for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-95% $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV, 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV and 0-12% central Au+Au at $\sqrt{s_{NN}}$ = 200 GeV. The dependence of the jet-like correlation is shown as a function of $\Delta\phi$ ($|\Delta\eta|<$ 1.78). The data are averaged between positive and negative $\Delta\phi$.

Dependence of jet-like yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV. 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of jet-like yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV. 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of jet-like yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for $d$+Au at $\sqrt{s_{NN}}$ = 200 GeV. 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of jet-like yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV. 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of jet-like yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV. 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of jet-like yield on $p_T^{\mathrm{trigger}}$ for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, 0-95% $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV, 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV and 0-12% central Au+Au at $\sqrt{s_{NN}}$ = 200 GeV. The 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of jet-like yield on $p_T^{\mathrm{associated}}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ for 0-60% Cu+Cu and 0-80% Au+Au collisions at $\sqrt{s_{NN}}$ = 62.4 GeV and 0-95% $d$+Au, 0-60% Cu+Cu, 0-12% Au+Au and 40-80% Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The 5% systematic error due to the uncertainty on the associated particle's efficiency is not shown.

Dependence of the widths in $\Delta\phi$ on $p_T^{\mathrm{trigger}}$ for 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-95% $d$+Au, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV and $\sqrt{s_{NN}}$ = 200 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, and 0-12% and 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\phi$ on $p_T^{\mathrm{associated}}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ for 0-95% $d$+Au, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV and $\sqrt{s_{NN}}$ = 200 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, and 0-12% and 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\phi$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of the widths in $\Delta\phi$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of the widths in $\Delta\phi$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-95% $d$+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\phi$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\phi$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\eta$ on $p_T^{\mathrm{trigger}}$ for 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-95% $d$+Au, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV and $\sqrt{s_{NN}}$ = 200 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, and 0-12% and 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\eta$ on $p_T^{\mathrm{associated}}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ for 0-95% $d$+Au, 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV and $\sqrt{s_{NN}}$ = 200 GeV, 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV, and 0-12% and 40-80% Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\eta$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of the widths in $\Delta\eta$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-80% Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of the widths in $\Delta\eta$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-95% $d$+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\eta$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for 0-60% Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of the widths in $\Delta\eta$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of ridge yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of ridge yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of ridge yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of ridge yield on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Ratio of the ridge and jet-like yields on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV.

Ratio of the ridge and jet-like yields on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV.

Ratio of the ridge and jet-like yields on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV.

Ratio of the ridge and jet-like yields on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of $V_{3\Delta}/V_{2\Delta}$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of $V_{3\Delta}/V_{2\Delta}$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV.

Dependence of $V_{3\Delta}/V_{2\Delta}$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for $d$+Au at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of $V_{3\Delta}/V_{2\Delta}$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Cu+Cu at $\sqrt{s_{NN}}$ = 200 GeV.

Dependence of $V_{3\Delta}/V_{2\Delta}$ on $N_{part}$ for 3 $<$ $p_T^{trigger}$ $<$ 6 GeV/$c$ and 1.5 GeV/$c$ $<$ $p_T^{associated}$ $<$ $p_T^{trigger}$ for Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

$\left.\frac{\mathrm{d}N}{\mathrm{d}\eta}\right\vert_{200\mathrm{Ge\!V}}\,/\,\left.\frac{\mathrm{d}N}{\mathrm{d}\eta}\right\vert_{130\mathrm{Ge\!V}}$ versus $\eta$ for $x^{\pm}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$ for , per centrality

$\frac{\mathrm{d}N}{\mathrm{d}y}\frac{2}{N_{\mathrm{part}}}$ versus $N_{\mathrm{part}}$ for $x^{\pm}$ in $\mathrm{Au}-\mathrm{Au}$ at $\sqrt{s_{\mathrm{NN}}}=200\,\mathrm{Ge\!V}$