Fig. A1: Number density $N_{ch}$ (left) and $\\Sigma p_{T}$ (right) distributions as a function of pleadingT and the comparisons to MC predictions in Toward (top), Transverse (middle), and Away (bottom) regions for $p_{T}^{track} >$ 0.5 GeV/$c$. The shaded areas in the upper panels represent the systematic uncertainties and vertical error bars indicate statistical uncertainties. In the lower panels, the shaded areas are the sum in quadrature of statistical and systematic uncertainties from the upper panels. No uncertainties are given for the MC calculations.

Fig. A2: Number density $N_{ch}$ (left) and $\\Sigma p_{T}$ (right) distributions as a function of pleadingT and the comparisons to MC predictions in Toward (top), Transverse (middle), and Away (bottom) regions for $p_{T}^{track} >$ 1.0 GeV/$c$. The shaded areas in the upper panels represent the systematic uncertainties and vertical error bars indicate statistical uncertainties. In the lower panels, the shaded areas are the sum in quadrature of statistical and systematic uncertainties from the upper panels. No uncertainties are given for the MC calculations.

Charged jet differential cross sections without UE subtraction in pp collisions at $\sqrt{s}$ = 5.02 TeV with the leading track bias. All jets must contain at least one track with $p_{T}$ > 5 GeV/$c$. Statistical uncertainties are displayed as vertical error bars. The total systematic uncertainties are shown as shaded bands around the data points. Data are scaled to enhance visibility

The second-order Fourier coefficients, $V_{3\Delta}(2, |\Delta\eta| > 2)$, as a function of $N_{offline}^{trk}$ for charged particles, after correcting for back-to-back jet correlations, estimated from the 10 $\leq$ $N_{offline}^{trk}$ < 20 range.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $N_{offline}^{trk}$ for charged particles.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $N_{offline}^{trk}$ for charged particles.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $N_{offline}^{trk}$ for charged particles.

The triangular flow after correcting for back-to-back jet correlations, $v_{3}^{sub}(2, |\Delta\eta| > 2)$, as a function of $N_{offline}^{trk}$ for charged particles.

The triangular flow after correcting for back-to-back jet correlations, $v_{3}^{sub}(2, |\Delta\eta| > 2)$, as a function of $N_{offline}^{trk}$ for charged particles.

The triangular flow after correcting for back-to-back jet correlations, $v_{3}^{sub}(2, |\Delta\eta| > 2)$, as a function of $N_{offline}^{trk}$ for charged particles.

The elliptic flow, $v_{2}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for charged particles.

The elliptic flow, $v_{2}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for charged particles.

The elliptic flow, $v_{2}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for charged particles.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for charged particles.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for charged particles.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for charged particles.

The elliptic flow, $v_{2}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for charged particles.

The elliptic flow, $v_{2}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for $K^{0}_{S}$.

The elliptic flow, $v_{2}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for $\Lambda/\bar{\Lambda}$.

The elliptic flow, $v_{2}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for charged particles.

The elliptic flow, $v_{2}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for $K^{0}_{S}$.

The elliptic flow, $v_{2}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for $\Lambda/\bar{\Lambda}$.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for charged particles.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for $K^{0}_{S}$.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $p_{T}$ for $\Lambda/\bar{\Lambda}$.

The elliptic flow per constituent quark after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)/n_{q}$, as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ for $K^{0}_{S}$.

The elliptic flow per constituent quark after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ for $\Lambda/\bar{\Lambda}$.

The four-particle cumulant, $c_{2}(4)$, as a function of $N_{offline}^{trk}$ for charged particles.

The four-particle cumulant, $c_{2}(4)$, as a function of $N_{offline}^{trk}$ for charged particles.

The four-particle cumulant, $c_{2}(4)$, as a function of $N_{offline}^{trk}$ for charged particles.

The six-particle cumulant, $c_{2}(6)$, as a function of $N_{offline}^{trk}$ for charged particles.

The elliptic flow after correcting for back-to-back jet correlations, $v_{2}^{sub}(2, |\Delta\eta| > 2)$, as a function of $N_{offline}^{trk}$ for charged particles.

The elliptic flow, $v_{2}(4)$, as a function of $N_{offline}^{trk}$ for charged particles.

The elliptic flow, $v_{2}(6)$, as a function of $N_{offline}^{trk}$ for charged particles.

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and $N_{offline}^{trk}<$ 35 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and $N_{offline}^{trk}<$ 35 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and $N_{offline}^{trk}<$ 35 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and $N_{offline}^{trk}<$ 35 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and $N_{offline}^{trk}<$ 35 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 0.1 $<p_{T}<$ 1.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 80 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 0.1 $<p_{T}<$ 1.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 90 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 80 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 90 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 80 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 90 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 80 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 90 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 0.1 $<p_{T}<$ 1.0 $GeV/c$ and 80 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 0.1 $<p_{T}<$ 1.0 $GeV/c$ and 90 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and 80 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and 90 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and 80 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and 90 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and 80 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and 90 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 0.1 $<p_{T}<$ 1.0 $GeV/c$ and $N_{offline}^{trk}>$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 0.1 $<p_{T}<$ 1.0 $GeV/c$ and $N_{offline}^{trk}>$ 110 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and $N_{offline}^{trk}>$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and $N_{offline}^{trk}>$ 110 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and $N_{offline}^{trk}>$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and $N_{offline}^{trk}>$ 110 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and $N_{offline}^{trk}>$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and $N_{offline}^{trk}>$ 110 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Associated yield for the near side of the correlation function averaged over 2 $<|\Delta\eta|<$ 4 and integrated over the region $|\Delta\Phi| < \Delta\Phi_{ZYAM}$ for pp data at $\sqrt{s} =$ 13 $TeV$. The associated yield as a function of $p_{T}$ for events with $N^{offline}_{trk} \geq$ 105. The $p_{T}$ value for each $p_{T}$ bin is the average $p_{T}$ value.

Associated yield for the near side of the correlation function averaged over 2 $<|\Delta\eta|<$ 4 for pp data at $\sqrt{s} =$ 7 $TeV$. The associated yield as a function of $p_{T}$ for events with $N^{offline}_{trk} \geq$ 110. The $p_{T}$ value for each $p_{T}$ bin is the average $p_{T}$ value.

Associated yield for the near side of the correlation function averaged over 2 $<|\Delta\eta|<$ 4 and integrated over the region $|\Delta\Phi| < \Delta\Phi_{ZYAM}$ for pp data at $\sqrt{s} =$ 13 $TeV$. The associated yield as a function of $N_{trk}^{offline}$ for events with 1.0 $< p_{T} <$ 2.0 GeV/c. The $N_{trk}^{offline}$ value for each $N_{trk}^{offline}$ bin is the average $N_{trk}^{offline}$ value.

Associated yield for the near side of the correlation function averaged over 2 $<|\Delta\eta|<$ 4 for pp data at $\sqrt{s} =$ 7 $TeV$. The associated yield as a function of $N_{trk}^{offline}$ for events with 1.0 $< p_{T} <$ 2.0 GeV/c. The $N_{trk}^{offline}$ value for each $N_{trk}^{offline}$ bin is the average $N_{trk}^{offline}$ value.

Associated yield for the near side of the correlation function averaged over 2 $<|\Delta\eta|<$ 4 pPb data at $\sqrt{s} =$ 5.02 $TeV$. The associated yield as a function of $N_{trk}^{offline}$ for events with 1.0 $< p_{T} <$ 2.0 GeV/c. The $N_{trk}^{offline}$ value for each $N_{trk}^{offline}$ bin is the average $N_{trk}^{offline}$ value.

Associated yield for the near side of the correlation function averaged over 2 $<|\Delta\eta|<$ 4 PbPb data at $\sqrt{s} =$ 2.76 $TeV$. The associated yield as a function of $N_{trk}^{offline}$ for events with 1.0 $< p_{T} <$ 2.0 GeV/c. The $N_{trk}^{offline}$ value for each $N_{trk}^{offline}$ bin is the average $N_{trk}^{offline}$ value.