The elliptic flow after nonflow subtraction, $v_{2}^{sub}/4$, for $f_0(980)$ as a function of $<KE_{T}>/4$ in pPb collision at 8.16 TeV.

The elliptic flow after nonflow subtraction, $v_{2}^{sub}/2$, for $f_0(980)$ as a function of $<p_{T}>/2$ in pPb collision at 8.16 TeV.

The elliptic flow after nonflow subtraction, $v_{2}^{sub}/4$, for $f_0(980)$ as a function of $<p_{T}>/4$ in pPb collision at 8.16 TeV.

$\Delta\phi$ projection of balance function in low $p_{T}$ in 0-10% centrality

$\Delta\phi$ projection of balance function in low $p_{T}$ in 30-40% centrality

1D $\Delta\eta$ projection of balance function in low $p_{T}$ in 70-80% centrality

$\Delta\eta$ projection of balance function in low $p_{T}$ in 0-40 $N_{trk}$

$\Delta\eta$ projection of balance function in low $p_{T}$ in 120-150 $N_{trk}$

$\Delta\eta$ projection of balance function in low $p_{T}$ in 270-300 $N_{trk}$

$\Delta\phi$ projection of balance function in low $p_{T}$ in 0-40 $N_{trk}$

$\Delta\phi$ projection of balance function in low $p_{T}$ in 120-150 $N_{trk}$

$\Delta\phi$ projection of balance function in low $p_{T}$ in 270-300 $N_{trk}$

The width of balance function in $<|\Delta\eta|>$ in low $p_{T}$ for PbPb collisions

Balance function width in low $p_{T}$ PbPb relative $<|\Delta\eta|>/<|\Delta\eta|>_{N_{ch}<65}$

The width of balance function in $<|\Delta\eta|>$ in low $p_{T}$ for pPb collisions

Balance function width in low $p_{T}$ in pPb relative $<|\Delta\eta|>/<|\Delta\eta|>_{N_{ch}<24}$

The width of the balance function in $<|\Delta\phi|>$ in low $p_{T}$ for PbPb collisions

Balance function width in low $p_{T}$ in PbPb for relative $<|\Delta\phi|>/<|\Delta\phi|>_{N_{ch}<65}$

The width of the balance function in $<|\Delta\phi|>$ in low $p_{T}$ for pPb collisions

Balance function width in low $p_{T}$ in pPb for relative $|\Delta\phi|/<|\Delta\phi|>_{N_{ch}<24}$

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in pPb for 0-40 $N_{trk}$

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in pPb for 120-150 $N_{trk}$

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in pPb for 270-300 $N_{trk}$

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in pPb for $N_{trk}$ 0-40

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in pPb for $N_{trk}$ 120-150

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in pPb for $N_{trk}$ 270-300

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in pPb for $N_{trk}$ 0-40

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in pPb for $N_{trk}$ 120-150

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in pPb for $N_{trk}$ 270-300

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for $N_{trk}$ 0-40

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for $N_{trk}$ 120-150

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for $N_{trk}$ 270-300

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for low $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for high $p_{T}$

Average width of the balance function in $\Delta \phi$ with $p_{T}$ for low $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for high $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for low $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for high $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for low $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for high $p_{T}$

Reconstruction level differential cross section spectla, obtained for the central acceptance, $|\eta| < 3$, for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology compared to the to the EPOS-LHC predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

The number of high purity tracksin the first $\eta$ bin after a gap of $4.5 \leq \Delta\eta^F < 5.0$ for events with the Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

The transeverse momentum of high purity tracksin the first $\eta$ bin after a gap of $4.5 \leq \Delta\eta^F < 5.0$ for events with the Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

The total energy of all Particle Flow candidatesin the first $\eta$ bin after a gap of $4.5 \leq \Delta\eta^F < 5.0$ for events with the Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

Unfolded diffraction enhanced differential cross section for events with Pomeron-Lead ($\mathrm{I\!P}\mathrm{Pb}$) topology, defined on the region $|\eta| < 5.2$, compared to hadron level predictions of the MC generators. The data are corrected for the contribution from events with undetectable energy in the HF calorimeter adjacent to the rapidity gap. The corrections are obtained using the EPOS-LHC MC sample. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffraction enhanced differential cross section for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology, defined on the region $|\eta| < 5.2$, compared to hadron level predictions of the MC generators. The data are corrected for the contribution from events with undetectable energy in the HF calorimeter adjacent to the rapidity gap. The corrections are obtained using the EPOS-LHC MC sample. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Lead ($\mathrm{I\!P}\mathrm{Pb}$) topology, compared to the to the hadron-level EPOS-LHC predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology, compared to the to the hadron-level EPOS-LHC predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Lead ($\mathrm{I\!P}\mathrm{Pb}$) topology, compared to the to the hadron-level QGSJET II predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology, compared to the to the hadron-level QGSJET II predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Reconstruction level diffraction enhanced differential cross section spectrum for Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology corrected for the contribution from events with undetectable energy in the HF calorimeter adjacent to the rapidity gap, compared with the spectrum obtained with events satisfying the ZDC veto requirement $E_{ZDC−} < 1 \mathrm{~TeV}$ which selects only the events without lead nuclear break up (without correction for HF undetectable energy). Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

The elliptic flow $v_2\{SP\}$ for charged hadron as a function of $p_T$ in PbPb collision at 5.02 TeV.

The elliptic flow $v_2\{4\}$ for $K_S^0$ as a function of $p_T$ in PbPb collision at 5.02 TeV.

The elliptic flow $v_2\{6\}$ for $K_S^0$ as a function of $p_T$ in PbPb collision at 5.02 TeV.

The elliptic flow $v_2\{8\}$ for $K_S^0$ as a function of $p_T$ in PbPb collision at 5.02 TeV.

The elliptic flow $v_2\{SP\}$ for $K_S^0$ as a function of $p_T$ in PbPb collision at 5.02 TeV.

The elliptic flow $v_2\{4\}$ for $\Lambda$ as a function of $p_T$ in PbPb collision at 5.02 TeV.

The elliptic flow $v_2\{6\}$ for $\Lambda$ as a function of $p_T$ in PbPb collision at 5.02 TeV.

The elliptic flow $v_2\{8\}$ for $\Lambda$ as a function of $p_T$ in PbPb collision at 5.02 TeV.

The elliptic flow $v_2\{SP\}$ for $\Lambda$ as a function of $p_T$ in PbPb collision at 5.02 TeV.

The elliptic flow $v_2\{4\}$ for charged hadron as a function of $p_T$ in pPb collision at 8.16 TeV.

The elliptic flow $v_2\{6\}$ for charged hadron as a function of $p_T$ in pPb collision at 8.16 TeV.

The elliptic flow $v_2\{8\}$ for charged hadron as a function of $p_T$ in pPb collision at 8.16 TeV.

The elliptic flow $v_2\{p-SP\}$ for charged hadron as a function of $p_T$ in pPb collision at 8.16 TeV.

The elliptic flow $v_2\{Pb-SP\}$ for charged hadron as a function of $p_T$ in pPb collision at 8.16 TeV.

The elliptic flow $v_2\{4\}$ for $K_S^0$ as a function of $p_T$ in pPb collision at 8.16 TeV.

The elliptic flow $v_2\{6\}$ for $K_S^0$ as a function of $p_T$ in pPb collision at 8.16 TeV.

The elliptic flow $v_2\{p-SP\}$ for $K_S^0$ as a function of $p_T$ in pPb collision at 8.16 TeV.

The elliptic flow $v_2\{Pb-SP\}$ for $K_S^0$ as a function of $p_T$ in pPb collision at 8.16 TeV.

The elliptic flow $v_2\{4\}$ for $\Lambda$ as a function of $p_T$ in pPb collision at 8.16 TeV.

The elliptic flow $v_2\{6\}$ for $\Lambda$ as a function of $p_T$ in pPb collision at 8.16 TeV.

The elliptic flow $v_2\{p-SP\}$ for $\Lambda$ as a function of $p_T$ in pPb collision at 8.16 TeV.

The elliptic flow $v_2\{Pb-SP\}$ for $\Lambda$ as a function of $p_T$ in pPb collision at 8.16 TeV.

The elliptic flow $v_2\{4,veto\}$ for charged hadron as a function of $p_T$ in pPb collision at 8.16 TeV.

The elliptic flow $v_2\{6,veto\}$ for charged hadron as a function of $p_T$ in pPb collision at 8.16 TeV.

The ratio $v_2\{4,veto\}/v_2\{4\}$ for charged hadron as a function of $p_T$ in pPb collision at 8.16 TeV.

The ratio $v_2\{6,veto\}/v_2\{6\}$ for charged hadron as a function of $p_T$ in pPb collision at 8.16 TeV.

The elliptic flow $v_2\{4,Sub\}$ for charged hadron as a function of $p_T$ in PbPb collision at 5.02 TeV.

The elliptic flow $v_2\{4,Sub\}$ for $K_S^0$ as a function of $p_T$ in PbPb collision at 5.02 TeV.

The elliptic flow $v_2\{4,Sub\}$ for $\Lambda$ as a function of $p_T$ in PbPb collision at 5.02 TeV.

The elliptic flow $v_2\{4,Sub\}$ for $K_S^0$ as a function of $p_T$ in pPb collision at 8.16 TeV.

The elliptic flow $v_2\{4,Sub\}$ for $\Lambda$ as a function of $p_T$ in pPb collision at 8.16 TeV.

The elliptic flow $v_2\{4,Sub\}$ for charged hadron as a function of $p_T$ in pPb collision at 8.16 TeV.

The ratio $v_2\{4,Sub\}/v_2\{4\}$ for charged hadron as a function of $p_T$ in PbPb collision at 5.02 TeV.

The ratio $v_2\{4,Sub\}/v_2\{4\}$ for $K_S^0$ as a function of $p_T$ in PbPb collision at 5.02 TeV.

The ratio $v_2\{4,Sub\}/v_2\{4\}$ for $\Lambda$ as a function of $p_T$ in PbPb collision at 5.02 TeV.

The ratio $v_2\{4,Sub\}/v_2\{4\}$ for charged hadron as a function of $p_T$ in pPb collision at 8.16 TeV.

The ratio $v_2\{4,Sub\}/v_2\{4\}$ for $K_S^0$ as a function of $p_T$ in pPb collision at 8.16 TeV.

The ratio $v_2\{4,Sub\}/v_2\{4\}$ for $\Lambda$ as a function of $p_T$ in pPb collision at 8.16 TeV.

The flow fluctuation $\sigma/\langle v_{2} \rangle$ for charged hadron as a function of $p_T$ in PbPb collision at 5.02 TeV.

The flow fluctuation $\sigma/\langle v_2 \rangle$ for $K_S^0$ as a function of $p_T$ in PbPb collision at 5.02 TeV.

The flow fluctuation $\sigma/\langle v_2 \rangle$ for $\Lambda$ as a function of $p_T$ in PbPb collision at 5.02 TeV.

The flow fluctuation $\sigma/\langle v_2 \rangle$ for $K_S^0$ as a function of $p_T$ in pPb collision at 8.16 TeV.

The flow fluctuation $\sigma/\langle v_2 \rangle$ for $\Lambda$ as a function of $p_T$ in pPb collision at 8.16 TeV.

The flow fluctuation $\sigma/\langle v_2 \rangle$ using $v_2\{4\}$ for charged hadron as a function of $p_T$ in pPb collision at 8.16 TeV.

The flow fluctuation $\sigma/\langle v_2 \rangle$ using $v_2\{6\}$ for charged hadron as a function of $p_T$ in pPb collision at 8.16 TeV.

The flow fluctuation $\sigma/\langle v_2 \rangle$ using $v_2\{4,veto\}$ for charged hadron as a function of $p_T$ in pPb collision at 8.16 TeV.

The $v_2\{6\}/v_2\{4\}$ ratio for charged hadron as a function of $p_T$ in PbPb collision at 5.02 TeV.

The $v_2\{6\}/v_2\{4\}$ ratio for $K_S^0$ as a function of $p_T$ in PbPb collision at 5.02 TeV.

The $v_2\{6\}/v_2\{4\}$ ratio for $\Lambda$ as a function of $p_T$ in PbPb collision at 5.02 TeV.

The $v_2\{8\}/v_2\{4\}$ ratio for charged hadron as a function of $p_T$ in PbPb collision at 5.02 TeV.

The $v_2\{8\}/v_2\{4\}$ ratio for $K_S^0$ as a function of $p_T$ in PbPb collision at 5.02 TeV.

The $v_2\{8\}/v_2\{4\}$ ratio for $\Lambda$ as a function of $p_T$ in PbPb collision at 5.02 TeV.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $\phi^*$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential fiducial cross section (without the acceptance correction) for the DY process measured in the muon channel, as a function of $\phi^*$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Integrated fiducial cross section for the DY process measured in the muon channel. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of dimuon invariant mass. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of rapidity in the centre-of-mass frame for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of rapidity in the centre-of-mass frame for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $p_{\textrm{T}}$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $\phi^*$ for $15<m_{\mu\mu}<60$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Differential cross section for the DY process measured in the muon channel, as a function of $\phi^*$ for $60<m_{\mu\mu}<120$ GeV. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Integrated cross section for the DY process measured in the muon channel. The quoted error is the quadratic sum of the statistical and systematic uncertainties. The global normalisation uncertainty of 3.5% is listed separately.

Forward-backward ratios for $15<m_{\mu\mu}<60$ GeV. Quoted uncertainties are the quadratic sum of the statistical and systematic uncertainties.

Forward-backward ratios for $60<m_{\mu\mu}<120$ GeV. Quoted uncertainties are the quadratic sum of the statistical and systematic uncertainties.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of dimuon invariant mass.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of rapidity in the centre-of-mass frame for $15<m_{\mu\mu}<60$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of rapidity in the centre-of-mass frame for $60<m_{\mu\mu}<120$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $p_{\textrm{T}}$ for $15<m_{\mu\mu}<60$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $p_{\textrm{T}}$ for $60<m_{\mu\mu}<120$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $\phi^*$ for $15<m_{\mu\mu}<60$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $\phi^*$ for $60<m_{\mu\mu}<120$ GeV.

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of rapidity in the centre-of-mass frame, putting together low and high masses. Each bin is labeled lm_i or hm_i, where lm stands for $15<m_{\mu\mu}<60$ GeV, hm stands for $60<m_{\mu\mu}<120$ GeV, and i is the corresponding bin number in rapidity in the centre-of-mass frame. The global normalisation uncertainty of 3.5% is not included. Data from Figure [3, 4, 4].

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $p_{\textrm{T}}$, putting together low and high masses. Each bin is labeled lm_i or hm_i, where lm stands for $15<m_{\mu\mu}<60$ GeV, hm stands for $60<m_{\mu\mu}<120$ GeV, and i is the corresponding bin number in $p_{\textrm{T}}$. The global normalisation uncertainty of 3.5% is not included. Data from Figure [3, 4, 4].

Correlation matrix for the systematic uncertainties, excluding integrated luminosity, as a function of $\phi^*$, putting together low and high masses. Each bin is labeled lm_i or hm_i, where lm stands for $15<m_{\mu\mu}<60$ GeV, hm stands for $60<m_{\mu\mu}<120$ GeV, and i is the corresponding bin number in $\phi^*$. The global normalisation uncertainty of 3.5% is not included. Data from Figure [3, 4, 4].

Results of elliptic flow, corrected for short range correlations, for prompt neutral D mesons, as a function of multiplicity for $|y_{lab}|< 1$, with 6$ < p_{T} < $8 GeV in pp collisions at 13 TeV.

Results of elliptic flow, corrected for short range correlations, for prompt neutral D mesons, as a function of multiplicity for $|y_{lab}|< 1$, with 2$ < p_{T} < $4 GeV in pPb collisions at 8.16 TeV.

Results of elliptic flow, corrected for short range correlations, for prompt neutral D mesons, as a function of multiplicity for $|y_{lab}|< 1$, with 4$ < p_{T} < $6 GeV in pPb collisions at 8.16 TeV.

Results of elliptic flow, corrected for short range correlations, for prompt neutral D mesons, as a function of multiplicity for $|y_{lab}|< 1$, with 6$ < p_{T} < $8 GeV in pPb collisions at 8.16 TeV.

Results of elliptic flow, corrected for short range correlations, for nonprompt neutral D mesons, as a function of transverse momenta for $|y_{lab}|< 1$, with $185 \leq N^{offline}_{trk} < 250$ in pPb collisions at 8.16 TeV.

$v_3\{4\}$ as a function of $N_{trk}^{offline}$ in PbPb collisions at $\sqrt{s_{NN}} = 5.02$ TeV.

$v_2\{4\}$ as a function of $N_{trk}^{offline}$ in pPb collisions at $\sqrt{s_{NN}} = 8.16$ TeV.

$v_2\{6\}$ as a function of $N_{trk}^{offline}$ in pPb collisions at $\sqrt{s_{NN}} = 8.16$ TeV.

$v_2\{6\}$ as a function of $N_{trk}^{offline}$ in pPb collisions at $\sqrt{s_{NN}} = 8.16$ TeV.

$v_3\{4\}$ as a function of $N_{trk}^{offline}$ in pPb collisions at $\sqrt{s_{NN}} = 8.16$ TeV.

$v_2\{4\}/v_2\{2\}$ as a function of $N_{trk}^{offline}$ in PbPb collisions at $\sqrt{s_{NN}} = 5.02$ TeV.

$v_3\{4\}/v_3\{2\}$ as a function of $N_{trk}^{offline}$ in PbPb collisions at $\sqrt{s_{NN}} = 5.02$ TeV.

$v_2\{4\}/v_2\{2\}$ as a function of $N_{trk}^{offline}$ in pPb collisions at $\sqrt{s_{NN}} = 8.16$ TeV.

$v_3\{4\}/v_3\{2\}$ as a function of $N_{trk}^{offline}$ in pPb collisions at $\sqrt{s_{NN}} = 8.16$ TeV.

$v_2\{6\}/v_2\{4\}$ as a function of $v_2\{4\}/v_2\{2\}$ in pPb collisions at $\sqrt{s_{NN}} = 8.16$ TeV.

$v_2\{8\}/v_2\{6\}$ as a function of $v_2\{4\}/v_2\{2\}$ in pPb collisions at $\sqrt{s_{NN}} = 8.16$ TeV.

The elliptic flow, $v_{2}$, for $\Omega^{-}$ as a function of $p_{T}$ in pPb collision at 8.16 TeV.

The elliptic flow, $v_{2}$, for $D^{0}$ mesons as a function of $p_{T}$ in pPb collision at 8.16 TeV.

The elliptic flow after correcting back-to-back jet contribution, $v_{2}^{sub}$, for $K^{0}_{S}$ as a function of $p_{T}$ in pPb collision at 8.16 TeV.

The elliptic flow after correcting back-to-back jet contribution, $v_{2}^{sub}$, for $\Lambda$ as a function of $p_{T}$ in pPb collision at 8.16 TeV.

The elliptic flow after correcting back-to-back jet contribution, $v_{2}^{sub}$, for $\Xi^{-}$ as a function of $p_{T}$ in pPb collision at 8.16 TeV.

The elliptic flow after correcting back-to-back jet contribution, $v_{2}^{sub}$, for $\Omega^{-}$ as a function of $p_{T}$ in pPb collision at 8.16 TeV.

The elliptic flow after correcting back-to-back jet contribution, $v_{2}^{sub}$, for $D^{0}$ mesons as a function of $p_{T}$ in pPb collision at 8.16 TeV.

The elliptic flow per constituent quark after correcting back-to-back jet contribution, $v_{2}^{sub}/n_{q}$, for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ in pPb collision at 8.16 TeV.

The elliptic flow per constituent quark after correcting back-to-back jet contribution, $v_{2}^{sub}/n_{q}$, for $\Lambda$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ in pPb collision at 8.16 TeV.

The elliptic flow per constituent quark after correcting back-to-back jet contribution, $v_{2}^{sub}/n_{q}$, for $\Xi^{-}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ in pPb collision at 8.16 TeV.

The elliptic flow per constituent quark after correcting back-to-back jet contribution, $v_{2}^{sub}/n_{q}$, for $\Omega^{-}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ in pPb collision at 8.16 TeV.

The elliptic flow per constituent quark after correcting back-to-back jet contribution, $v_{2}^{sub}/n_{q}$, for $D^{0}$ meson as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ in pPb collision at 8.16 TeV.

The elliptic flow after correcting back-to-back jet contribution, $v_{2}^{sub}$, for $K^{0}_{S}$ as a function of $p_{T}$ in PbPb collision at 5.02 TeV.

The elliptic flow after correcting back-to-back jet contribution, $v_{2}^{sub}$, for $\Lambda$ as a function of $p_{T}$ in PbPb collision at 5.02 TeV.

The elliptic flow after correcting back-to-back jet contribution, $v_{2}^{sub}$, for $\Xi^{-}$ as a function of $p_{T}$ in PbPb collision at 5.02 TeV.

The elliptic flow after correcting back-to-back jet contribution, $v_{2}^{sub}$, for $\Omega^{-}$ as a function of $p_{T}$ in PbPb collision at 5.02 TeV.

The elliptic flow per constituent quark after correcting back-to-back jet contribution, $v_{2}^{sub}/n_{q}$, for $K^{0}_{S}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ in PbPb collision at 5.02 TeV.

The elliptic flow per constituent quark after correcting back-to-back jet contribution, $v_{2}^{sub}/n_{q}$, for $\Lambda$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ in PbPb collision at 5.02 TeV.

The elliptic flow per constituent quark after correcting back-to-back jet contribution, $v_{2}^{sub}/n_{q}$, for $\Xi^{-}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ in PbPb collision at 5.02 TeV.

The elliptic flow per constituent quark after correcting back-to-back jet contribution, $v_{2}^{sub}/n_{q}$, for $\Omega^{-}$ as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ in PbPb collision at 5.02 TeV.

The elliptic flow per constituent quark after correcting back-to-back jet contribution, $v_{2}^{sub}/n_{q}$, for $D^{0}$ meson as a function of transverse kinetic energy per constituent quark $KE_{T}/n_{q}$ in PbPb collision at 5.02 TeV.