Fig. 6b: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Away and Toward regions after the subtraction of Number density $N_{\rm ch}$ and $\Sigma p_{\rm T}$ distributions in the transverse region for p-Pb collisions for $p_{\rm T} >$ 0.5 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. 7a: $<p_{\rm T}>$ as a function of $p_{\rm T}^{\rm trig}$ in Toward (left) and Away (right) regions after the subtraction of transverse region for pp collisions for $p_{\rm T} >$ 0.5 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. 7b: $<p_{\rm T}>$ as a function of $p_{\rm T}^{\rm trig}$ in Toward (left) and Away (right) regions after the subtraction of transverse region for p-Pb collisions for $p_{\rm T} >$ 0.5 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A1: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 0.15 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A2: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 1.0 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A3: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 0.15 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A4: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 1.0 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Single-differential production cross-sections for nonprompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainties are statistical, the second are correlated systematic uncertainties shared between intervals, and the last are uncorrelated systematic uncertainties.

Single-differential production cross-sections for prompt $J/\psi$ mesons as a function of $y$. The first uncertainties are statistical, the second are correlated systematic uncertainties shared between intervals, and the last are uncorrelated systematic uncertainties.

Single-differential production cross-sections for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainties are statistical, the second are correlated systematic uncertainties shared between intervals, and the last are uncorrelated systematic uncertainties.

Fraction of nonprompt $J/\psi$ mesons (in \%) in ($p_\text{T},y$) intervals. The first uncertainty is statistical and the second is systematic.

Nuclear modification factor $R_{p\text{Pb}}$ for prompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Nuclear modification factor $R_{p\text{Pb}}$ for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Relative changes of cross-sections (in \%), for a polarisation of $\lambda_{\theta}=-0.2$ rather than zero, in ($p_\text{T},y$) intervals.

Relative changes of cross-sections (in \%), for a polarisation of $\lambda_{\theta}=-1$ rather than zero, in ($p_\text{T},y$) intervals.

Relative changes of cross-sections (in \%), for a polarisation of $\lambda_{\theta}=+1$ rather than zero, in ($p_\text{T},y$) intervals.

The Z boson differential cross section as a function of |yZ|.

The Z boson differential cross section as a function of |yZ|.

The Z boson differential cross section as a function of |yZ|.

The Z boson differential cross section as a function of pTZ.

The Z boson differential cross section as a function of pTZ.

The Z boson differential cross section as a function of pTZ.

The TAA-normalized yields of Z bosons as a function of centrality.

The TAA-normalized yields of Z bosons as a function of centrality.

The TAA-normalized yields of Z bosons as a function of centrality.

Distributions of $\Delta\phi_{\mathrm{trk,Z}}$ in 30-50% centrality PbPb collisions at 5.02 TeV.

Distributions of $\Delta\phi_{\mathrm{trk,Z}}$ in 0-30% centrality PbPb collisions at 5.02 TeV.

PbPb - pp difference of $\Delta\phi_{\mathrm{trk,Z}}$ distributions for 70-90% centrality PbPb collisions at 5.02 TeV.

PbPb - pp difference of $\Delta\phi_{\mathrm{trk,Z}}$ distributions for 50-70% centrality PbPb collisions at 5.02 TeV.

PbPb - pp difference of $\Delta\phi_{\mathrm{trk,Z}}$ distributions for 30-50% centrality PbPb collisions at 5.02 TeV.

PbPb - pp difference of $\Delta\phi_{\mathrm{trk,Z}}$ distributions for 0-30% centrality PbPb collisions at 5.02 TeV.

Distributions of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ in pp collisions at 5.02 TeV.

Distributions of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ in 70-90% centrality PbPb collisions at 5.02 TeV.

Distributions of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ in 50-70% centrality PbPb collisions at 5.02 TeV.

Distributions of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ in 30-50% centrality PbPb collisions at 5.02 TeV.

Distributions of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ in 0-30% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ distributions for 70-90% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ distributions for 50-70% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ distributions for 30-50% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of $\xi^{\mathrm{trk,Z}}_{\mathrm{T}}$ distributions for 0-30% centrality PbPb collisions at 5.02 TeV.

Distributions of p$^{\mathrm{trk}}_{\mathrm{T}}$ in pp collisions at 5.02 TeV.

Distributions of p$^{\mathrm{trk}}_{\mathrm{T}}$ in pp collisions at 5.02 TeV.

Distributions of p$^{\mathrm{trk}}_{\mathrm{T}}$ in 70-90% centrality PbPb collisions at 5.02 TeV.

Distributions of p$^{\mathrm{trk}}_{\mathrm{T}}$ in 50-70% centrality PbPb collisions at 5.02 TeV.

Distributions of p$^{\mathrm{trk}}_{\mathrm{T}}$ in 30-50% centrality PbPb collisions at 5.02 TeV.

Distributions of p$^{\mathrm{trk}}_{\mathrm{T}}$ in 0-30% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of p$^{\mathrm{trk}}_{\mathrm{T}}$ distributions for 70-90% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of p$^{\mathrm{trk}}_{\mathrm{T}}$ distributions for 50-70% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of p$^{\mathrm{trk}}_{\mathrm{T}}$ distributions for 30-50% centrality PbPb collisions at 5.02 TeV.

PbPb / pp ratio of p$^{\mathrm{trk}}_{\mathrm{T}}$ distributions for 0-30% centrality PbPb collisions at 5.02 TeV.

Neutron multiplicity dependence of acoplanarity ($\alpha$) from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Neutron multiplicity dependence of acoplanarity ($\alpha$) from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Neutron multiplicity dependence of acoplanarity ($\alpha$) from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Neutron multiplicity dependence of average core acoplanarity ($<\alpha^{core}>$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Neutron multiplicity dependence of average invariant mass ($<m_{\mu\mu}>$) from process $\gamma\gamma$ to $\mu^+\mu^-$ in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 0n1n class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing larger neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 0nXn class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing larger neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 1nXn class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing larger neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 0n1n class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing smaller neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 0nXn class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing smaller neutron multiplicity.

Acoplanarity ($\alpha$) distributions from process $\gamma\gamma$ to $\mu^+\mu^-$ for 1nXn class in ultraperipheral PbPb at $\sqrt{s_{NN}}=5.02$ TeV. Dimuon rapidity is in the hemisphere containing smaller neutron multiplicity.

$p_{\rm T}$-differential yields of $\Lambda$(1520) (sum of particle and anti-particle states) in p--Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ $\mathrm{=}$ 5.02 TeV in multiplicity interval 20--40\%. The uncertainty 'sys,$p_{\rm T}$-correlated' indicates the systematic uncertainty after removing the contributions of $p_{\rm T}$-uncorrelated uncertainty.

$p_{\rm T}$-differential yields of $\Lambda$(1520) (sum of particle and anti-particle states) in p--Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ $\mathrm{=}$ 5.02 TeV in multiplicity interval 40--60\%. The uncertainty 'sys,$p_{\rm T}$-correlated' indicates the systematic uncertainty after removing the contributions of $p_{\rm T}$-uncorrelated uncertainty.

$p_{\rm T}$-differential yields of $\Lambda$(1520) (sum of particle and anti-particle states) in p--Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ $\mathrm{=}$ 5.02 TeV in multiplicity interval 60--100\%. The uncertainty 'sys,$p_{\rm T}$-correlated' indicates the systematic uncertainty after removing the contributions of $p_{\rm T}$-uncorrelated uncertainty.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to charged $\pi$ ($\pi^{+} + \pi^{-}$) at midrapidity as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in inelastic pp collisions at $\sqrt{s}$ $\mathrm{=}$ 7 TeV.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to charged $\pi$ ($\pi^{+} + \pi^{-}$) as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in p--Pb collisions at $\sqrt{s}$ $\mathrm{=}$ 5.02 TeV. The uncertainty 'sys,uncorrelated' indicates the multiplicity uncorrelated systematic uncertainty.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to K$^{-}$ at midrapidity as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in inelastic pp collisions at $\sqrt{s}$ $\mathrm{=}$ 7 TeV.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to K$^{-}$ as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in p--Pb collisions at $\sqrt{s}$ $\mathrm{=}$ 5.02 TeV. The uncertainty 'sys,uncorrelated' indicates the multiplicity uncorrelated systematic uncertainty.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to its ground state particle $\Lambda$ ($\Lambda + \bar{\Lambda}$) at midrapidity as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in inelastic pp collisions at $\sqrt{s}$ $\mathrm{=}$ 7 TeV.

Ratio of $p_{\rm T}$-integrated yields of $\Lambda$(1520) (sum of particle and anti-particle states) to its ground state particle $\Lambda$ ($\Lambda + \bar{\Lambda}$) as a function of $\langle {\rm d}N_{\rm ch}/{\rm d} \eta_{\rm lab} \rangle$ in p--Pb collisions at $\sqrt{s}$ $\mathrm{=}$ 5.02 TeV. The uncertainty 'sys,uncorrelated' indicates the multiplicity uncorrelated systematic uncertainty.

95% CL intervals on the double ratio of measured yields, $(N_{\varUpsilon(3S)} / N_{\varUpsilon(1S)})_{\mathrm{PbPb}} / (N_{\varUpsilon(3S)} / N_{\varUpsilon(1S)})_{pp}$, as a function of centrality, for upsilon $|y|<2.4$ and $p_T<30$GeV, and $p_{T}^{\mu}>4$GeV.

68% CL intervals on the double ratio of measured yields, $(N_{\varUpsilon(3S)} / N_{\varUpsilon(1S)})_{\mathrm{PbPb}} / (N_{\varUpsilon(3S)} / N_{\varUpsilon(1S)})_{pp}$, as a function of centrality, for upsilon $|y|<2.4$ and $p_T<30$GeV, and $p_{T}^{\mu}>4$GeV.

Differential cross section (multiplied by the dimuon branching fraction) of prompt J/$\psi$ mesons in pPb collisions at backward $y_{\mathrm{CM}}$.

Rapidity dependence of the cross section (multiplied by the dimuon branching fraction) for prompt J/$\psi$ mesons in the $p_{\mathrm{T}}$ intervals 6.5<$p_{\mathrm{T}}$<10 GeV and 10<$p_{\mathrm{T}}$<30 GeV in pp collisions.

Rapidity dependence of the cross section (multiplied by the dimuon branching fraction) for prompt J/$\psi$ mesons in the $p_{\mathrm{T}}$ intervals 6.5<$p_{\mathrm{T}}$<10 GeV and 10<$p_{\mathrm{T}}$<30 GeV in pPb collisions.

Transverse momentum dependence of $R_{\mathrm{pPb}}$ for prompt J/$\psi$ mesons in seven $y_{\mathrm{CM}}$ ranges.

Rapidity dependence of $R_{\mathrm{pPb}}$ for prompt J/$\psi$ mesons in $p_{\mathrm{T}}$ range 6.5<$p_{\mathrm{T}}$<10 GeV.

Rapidity dependence of $R_{\mathrm{pPb}}$ for prompt J/$\psi$ mesons in $p_{\mathrm{T}}$ range 10<$p_{\mathrm{T}}$<30 GeV.

Transverse momentum dependence of $R_{\mathrm{FB}}$ for prompt J/$\psi$ mesons in the 0<$y_{\mathrm{CM}}$<0.9 range.

Transverse momentum dependence of $R_{\mathrm{FB}}$ for prompt J/$\psi$ mesons in the 0.9<$y_{\mathrm{CM}}$<1.5 range.

Transverse momentum dependence of $R_{\mathrm{FB}}$ for prompt J/$\psi$ mesons in the 1.5<$y_{\mathrm{CM}}$<1.93 range.

Dependence of $R_{\mathrm{FB}}$ for prompt J/$\psi$ mesons on the hadronic activity in the event, given by the transverse energy deposited in the CMS detector at large pseudorapidities $E_{\mathrm{T}}^{\mathrm{HF}|\eta|>4}$.

Differential cross section (multiplied by the dimuon branching fraction) of nonprompt J/$\psi$ mesons in pp collisions at forward $y_{\mathrm{CM}}$.

Differential cross section (multiplied by the dimuon branching fraction) of nonprompt J/$\psi$ mesons in pp collisions at backward $y_{\mathrm{CM}}$.

Differential cross section (multiplied by the dimuon branching fraction) of nonprompt J/$\psi$ mesons in pPb collisions at forward $y_{\mathrm{CM}}$.

Differential cross section (multiplied by the dimuon branching fraction) of nonprompt J/$\psi$ mesons in pPb collisions at backward $y_{\mathrm{CM}}$.

Rapidity dependence of the cross section (multiplied by the dimuon branching fraction) for nonprompt J/$\psi$ mesons in the $p_{\mathrm{T}}$ intervals 6.5<$p_{\mathrm{T}}$<10 GeV and 10<$p_{\mathrm{T}}$<30 GeV in pp collisions.

Rapidity dependence of the cross section (multiplied by the dimuon branching fraction) for nonprompt J/$\psi$ mesons in the $p_{\mathrm{T}}$ intervals 6.5<$p_{\mathrm{T}}$<10 GeV and 10<$p_{\mathrm{T}}$<30 GeV in pPb collisions.

Transverse momentum dependence of $R_{\mathrm{pPb}}$ for nonprompt J/$\psi$ mesons in seven $y_{\mathrm{CM}}$ ranges.

Rapidity dependence of $R_{\mathrm{pPb}}$ for nonprompt J/$\psi$ mesons in $p_{\mathrm{T}}$ range 6.5<$p_{\mathrm{T}}$<10 GeV.

Rapidity dependence of $R_{\mathrm{pPb}}$ for nonprompt J/$\psi$ mesons in $p_{\mathrm{T}}$ range 10<$p_{\mathrm{T}}$<30 GeV.

Transverse momentum dependence of $R_{\mathrm{FB}}$ for nonprompt J/$\psi$ mesons in the 0<$y_{\mathrm{CM}}$<0.9 range.

Transverse momentum dependence of $R_{\mathrm{FB}}$ for nonprompt J/$\psi$ mesons in the 0.9<$y_{\mathrm{CM}}$<1.5 range.

Transverse momentum dependence of $R_{\mathrm{FB}}$ for nonprompt J/$\psi$ mesons in the 1.5<$y_{\mathrm{CM}}$<1.93 range.

Dependence of $R_{\mathrm{FB}}$ for nonprompt J/$\psi$ mesons on the hadronic activity in the event, given by the transverse energy deposited in the CMS detector at large pseudorapidities $E_{\mathrm{T}}^{\mathrm{HF}|\eta|>4}$.

Distribution of the transverse momentum ratio between the Z boson and the jet in PbPb collisions. The distribution is normalized by the number of Z events.

The mean value of the X_JZ distribution in pp collisions. The resolutions of the measured jet energy and azimuthal angle are smeared to match those of the PbPb sample.

The mean value of the X_JZ distribution in PbPb collisions.

The average number of jet partners per Z boson in pp collisions. The resolutions of the measured jet energy and azimuthal angle are smeared to match those of the PbPb sample.

The average number of jet partners per Z boson in PbPb collisions.

Distribution of the transverse momentum ratio between the Z boson and the jet in pp collisions. The distribution is normalized by the number of Z events. No resolution smearing is applied to the measured jet energy and azimuthal angle.

The $v_{2}$ result from SP method as a function of $p_{T}$ in 20-30\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}$ result from SP method as a function of $p_{T}$ in 30-40\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}$ result from SP method as a function of $p_{T}$ in 40-50\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}$ result from SP method as a function of $p_{T}$ in 50-60\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{3}$ result from SP method as a function of $p_{T}$ in 0-5\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{3}$ result from SP method as a function of $p_{T}$ in 5-10\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{3}$ result from SP method as a function of $p_{T}$ in 10-20\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{3}$ result from SP method as a function of $p_{T}$ in 20-30\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{3}$ result from SP method as a function of $p_{T}$ in 30-40\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{3}$ result from SP method as a function of $p_{T}$ in 40-50\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{3}$ result from SP method as a function of $p_{T}$ in 50-60\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}$ result from 4-, 6- and 8-particle cumulant methods as a function of $p_{T}$ in 5-10\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}$ result from 4-, 6- and 8-particle cumulant methods as a function of $p_{T}$ in 10-20\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}$ result from 4-, 6- and 8-particle cumulant methods as a function of $p_{T}$ in 20-30\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}$ result from 4-, 6- and 8-particle cumulant methods as a function of $p_{T}$ in 30-40\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}$ result from 4-, 6- and 8-particle cumulant methods as a function of $p_{T}$ in 40-50\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}$ result from 4-, 6- and 8-particle cumulant methods as a function of $p_{T}$ in 50-60\% centrality bin of PbPb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Shaded boxes represent systematic uncertainties.

The $v_{2}^{high}$ as a function of $v_{2}^{low}$ results from SP method in PbPb collisions at $sqrt{s_{NN}}$ = 5.02 TeV. Only statistical uncertainties are shown.

The $v_{2}^{high}$ as a function of $v_{2}^{low}$ results from 4-particle cumulant method in PbPb collisions at $sqrt{s_{NN}}$ = 5.02 TeV. Only statistical uncertainties are shown.