The correlation matrix corresponding to the statistical uncertainties on the differential cross-section ($d\sigma/dp_T$) fit results for $W^-$. To combine with $W^+$ results (Fig8a), use the rows and columns ordered as $W^+$ and then $W^-$. Assume no correlation in the statistical uncertainties between $W^+$ and $W^-$ (zero entries in the off-diagonal blocks).

The correlation matrix corresponding to all of the systematic uncertainties on the differential cross-section ($d\sigma/dp_T$) fit results for $W^+$ and $W^-$ combined together, with rows and columns ordered as $W^+$ and then $W^-$.

Data from Figure 2, panel d, upper panel, $v_{0}(p_{T})$ for $\eta_{gap}$ = 1, 0-5% Centrality

Data from Figure 2, panel d, lower panel, $v_{0}(p_{T})$ for $\eta_{gap}$ = 1, 50-60% Centrality

Data from Figure 3, panel a, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, 0-5% Centrality, $\eta_{gap}$=0

Data from Figure 3, panel a, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, 0-5% Centrality, $\eta_{gap}$=1

Data from Figure 3, panel a, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, 0-5% Centrality, $\eta_{gap}$=2

Data from Figure 3, panel a, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, 0-5% Centrality, $\eta_{gap}$=3

Data from Figure 3, panel b, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, 60-70% Centrality, $\eta_{gap}$=0

Data from Figure 3, panel b, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, 60-70% Centrality, $\eta_{gap}$=1

Data from Figure 3, panel b, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, 60-70% Centrality, $\eta_{gap}$=2

Data from Figure 3, panel b, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, 60-70% Centrality, $\eta_{gap}$=3

Data from Figure 4, panel a, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$=1, 0-5% Centrality

Data from Figure 4, panel a, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$=1, 5-10% Centrality

Data from Figure 4, panel a, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$=1, 10-20% Centrality

Data from Figure 4, panel a, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$=1, 20-30% Centrality

Data from Figure 4, panel a, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$=1, 30-40% Centrality

Data from Figure 4, panel a, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$=1, 40-50% Centrality

Data from Figure 4, panel a, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$=1, 50-60% Centrality

Data from Figure 4, panel a, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$=1, 60-70% Centrality

Data from Figure 4, panel b, $v_{0}(p_{T})/v_{0}$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$=1, 0-5% Centrality

Data from Figure 4, panel b, $v_{0}(p_{T})/v_{0}$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$=1, 5-10% Centrality

Data from Figure 4, panel b, $v_{0}(p_{T})/v_{0}$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$=1, 10-20% Centrality

Data from Figure 4, panel b, $v_{0}(p_{T})/v_{0}$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$=1, 20-30% Centrality

Data from Figure 4, panel b, $v_{0}(p_{T})/v_{0}$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$=1, 30-40% Centrality

Data from Figure 4, panel b, $v_{0}(p_{T})/v_{0}$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$=1, 40-50% Centrality

Data from Figure 4, panel b, $v_{0}(p_{T})/v_{0}$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$=1, 50-60% Centrality

Data from Figure 4, panel b, $v_{0}(p_{T})/v_{0}$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$=1, 60-70% Centrality

Data from Figure 5, panel a, $v_{0}(p_{T})/v_{0}$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$=1, 0-5% Centrality

Data from Figure 5, panel b, $v_{0}(p_{T})/v_{0}$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$=1, 0-5% Centrality

Data from Appendix, Figure 6, panel a, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$ = 1, 0-5% Centrality

Data from Appendix, Figure 6, panel b, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$ = 1, 10-20% Centrality

Data from Appendix, Figure 6, panel c, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$ = 1, 30-40% Centrality

Data from Appendix, Figure 6, panel d, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$ = 1, 60-70% Centrality

Data from Appendix, Figure 7, Zero Crossing Point of $v_{0}(p_{T})$ for $\eta_{gap}$ = 1

Data from Appendix, Figure 7, $\left\langle [p_{T}]\right\rangle$ for $p_{T}$ range of 0.5-10 GeV, $\eta_{gap}$ = 1

Data from Appendix, Figure 8, panel a, Closure of Sum Rule 1 for $\eta_{gap}$ = 1

Data from Appendix, Figure 8, panel b, Closure of Sum Rule 2 for $\eta_{gap}$ = 1

Data from Appendix, Figure 9, panel a, $v_{0}$ vs $N_{ch}$ for $\eta_{gap}$ = 1

Data from Appendix, Figure 9, panel b, $v_{0} / v^{5\%}_{0}$ vs $N_{ch}$ for $\eta_{gap}$ = 1

Data from Appendix, Figure 10, panel a, $v_{0}\sqrt{N_{ch}}$ vs $N_{ch}$ for $\eta_{gap}$ = 1

Data from Appendix, Figure 10, panel b, $v_{0}\sqrt{N_{ch}}$ vs Centrality for $\eta_{gap}$ = 1

Data from Appendix, Figure 11, $v_{0}(p_{T})\,\sqrt{\left\langle N_{\mathrm{ch}}\right\rangle}$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$ = 1, 0-5% Centrality

Data from Appendix, Figure 11, $v_{0}(p_{T})\,\sqrt{\left\langle N_{\mathrm{ch}}\right\rangle}$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$ = 1, 5-10% Centrality

Data from Appendix, Figure 11, $v_{0}(p_{T})\,\sqrt{\left\langle N_{\mathrm{ch}}\right\rangle}$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$ = 1, 10-20% Centrality

Data from Appendix, Figure 11, $v_{0}(p_{T})\,\sqrt{\left\langle N_{\mathrm{ch}}\right\rangle}$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$ = 1, 20-30% Centrality

Data from Appendix, Figure 11, $v_{0}(p_{T})\,\sqrt{\left\langle N_{\mathrm{ch}}\right\rangle}$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$ = 1, 30-40% Centrality

Data from Appendix, Figure 11, $v_{0}(p_{T})\,\sqrt{\left\langle N_{\mathrm{ch}}\right\rangle}$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$ = 1, 40-50% Centrality

Data from Appendix, Figure 11, $v_{0}(p_{T})\,\sqrt{\left\langle N_{\mathrm{ch}}\right\rangle}$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$ = 1, 50-60% Centrality

Data from Appendix, Figure 11, $v_{0}(p_{T})\,\sqrt{\left\langle N_{\mathrm{ch}}\right\rangle}$ for $p^{ref}_{T}$ = 0.5-2 GeV, $\eta_{gap}$ = 1, 60-70% Centrality

Data from Auxiliary, Figure 1, Top row, Left column, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-5 GeV, 0-5% Centrality, $\eta_{gap}$=0

Data from Auxiliary, Figure 1, Top row, Left column, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-5 GeV, 0-5% Centrality, $\eta_{gap}$=1

Data from Auxiliary, Figure 1, Top row, Left column, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-5 GeV, 0-5% Centrality, $\eta_{gap}$=2

Data from Auxiliary, Figure 1, Top row, Left column, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-5 GeV, 0-5% Centrality, $\eta_{gap}$=3

Data from Auxiliary, Figure 1, Top row, Right column, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-5 GeV, 60-70% Centrality, $\eta_{gap}$=0

Data from Auxiliary, Figure 1, Top row, Right column, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-5 GeV, 60-70% Centrality, $\eta_{gap}$=1

Data from Auxiliary, Figure 1, Top row, Right column, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-5 GeV, 60-70% Centrality, $\eta_{gap}$=2

Data from Auxiliary, Figure 1, Top row, Right column, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-5 GeV, 60-70% Centrality, $\eta_{gap}$=3

Data from Auxiliary, Figure 1, Bottom row, Left column, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-5 GeV, 0-5% Centrality, $\eta_{gap}$=0

Data from Auxiliary, Figure 1, Bottom row, Left column, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-5 GeV, 0-5% Centrality, $\eta_{gap}$=1

Data from Auxiliary, Figure 1, Bottom row, Left column, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-5 GeV, 0-5% Centrality, $\eta_{gap}$=2

Data from Auxiliary, Figure 1, Bottom row, Left column, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-5 GeV, 0-5% Centrality, $\eta_{gap}$=3

Data from Auxiliary, Figure 1, Bottom row, Right column, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-5 GeV, 60-70% Centrality, $\eta_{gap}$=0

Data from Auxiliary, Figure 1, Bottom row, Right column, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-5 GeV, 60-70% Centrality, $\eta_{gap}$=1

Data from Auxiliary, Figure 1, Bottom row, Right column, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-5 GeV, 60-70% Centrality, $\eta_{gap}$=2

Data from Auxiliary, Figure 1, Bottom row, Right column, $v_{0}(p_{T})$ for $p^{ref}_{T}$ = 0.5-5 GeV, 60-70% Centrality, $\eta_{gap}$=3

Data from Auxiliary, Figure 2, $v_{0}(p_{T})$ For $\eta_{gap}$ = 0, and 0-5% Centrality

Data from Auxiliary, Figure 2, $v_{0}(p_{T})$ For $\eta_{gap}$ = 0, and 5-10% Centrality

Data from Auxiliary, Figure 2, $v_{0}(p_{T})$ For $\eta_{gap}$ = 0, and 10-20% Centrality

Data from Auxiliary, Figure 2, $v_{0}(p_{T})$ For $\eta_{gap}$ = 0, and 20-30% Centrality

Data from Auxiliary, Figure 2, $v_{0}(p_{T})$ For $\eta_{gap}$ = 0, and 30-40% Centrality

Data from Auxiliary, Figure 2, $v_{0}(p_{T})$ For $\eta_{gap}$ = 0, and 40-50% Centrality

Data from Auxiliary, Figure 2, $v_{0}(p_{T})$ For $\eta_{gap}$ = 0, and 50-60% Centrality

Data from Auxiliary, Figure 2, $v_{0}(p_{T})$ For $\eta_{gap}$ = 0, and 60-70% Centrality

Data from Auxiliary, Figure 3, $v_{0}(p_{T})$ For $\eta_{gap}$ = 1, and 0-5% Centrality

Data from Auxiliary, Figure 3, $v_{0}(p_{T})$ For $\eta_{gap}$ = 1, and 5-10% Centrality

Data from Auxiliary, Figure 3, $v_{0}(p_{T})$ For $\eta_{gap}$ = 1, and 10-20% Centrality

Data from Auxiliary, Figure 3, $v_{0}(p_{T})$ For $\eta_{gap}$ = 1, and 20-30% Centrality

Data from Auxiliary, Figure 3, $v_{0}(p_{T})$ For $\eta_{gap}$ = 1, and 30-40% Centrality

Data from Auxiliary, Figure 3, $v_{0}(p_{T})$ For $\eta_{gap}$ = 1, and 40-50% Centrality

Data from Auxiliary, Figure 3, $v_{0}(p_{T})$ For $\eta_{gap}$ = 1, and 50-60% Centrality

Data from Auxiliary, Figure 3, $v_{0}(p_{T})$ For $\eta_{gap}$ = 1, and 60-70% Centrality

Data from Auxiliary, Figure 4, $v_{0}(p_{T})$ For $\eta_{gap}$ = 2, and 0-5% Centrality

Data from Auxiliary, Figure 4, $v_{0}(p_{T})$ For $\eta_{gap}$ = 2, and 5-10% Centrality

Data from Auxiliary, Figure 4, $v_{0}(p_{T})$ For $\eta_{gap}$ = 2, and 10-20% Centrality

Data from Auxiliary, Figure 4, $v_{0}(p_{T})$ For $\eta_{gap}$ = 2, and 20-30% Centrality

Data from Auxiliary, Figure 4, $v_{0}(p_{T})$ For $\eta_{gap}$ = 2, and 30-40% Centrality

Data from Auxiliary, Figure 4, $v_{0}(p_{T})$ For $\eta_{gap}$ = 2, and 40-50% Centrality

Data from Auxiliary, Figure 4, $v_{0}(p_{T})$ For $\eta_{gap}$ = 2, and 50-60% Centrality

Data from Auxiliary, Figure 4, $v_{0}(p_{T})$ For $\eta_{gap}$ = 2, and 60-70% Centrality

Data from Auxiliary, Figure 5, $v_{0}(p_{T})$ For $\eta_{gap}$ = 3, and 0-5% Centrality

Data from Auxiliary, Figure 5, $v_{0}(p_{T})$ For $\eta_{gap}$ = 3, and 5-10% Centrality

Data from Auxiliary, Figure 5, $v_{0}(p_{T})$ For $\eta_{gap}$ = 3, and 10-20% Centrality

Data from Auxiliary, Figure 5, $v_{0}(p_{T})$ For $\eta_{gap}$ = 3, and 20-30% Centrality

Data from Auxiliary, Figure 5, $v_{0}(p_{T})$ For $\eta_{gap}$ = 3, and 30-40% Centrality

Data from Auxiliary Figure 5, $v_{0}(p_{T})$ For $\eta_{gap}$ = 3, and 40-50% Centrality

Data from Auxiliary Figure 5, $v_{0}(p_{T})$ For $\eta_{gap}$ = 3, and 50-60% Centrality

Data from Auxiliary Figure 5, $v_{0}(p_{T})$ For $\eta_{gap}$ = 3, and 60-70% Centrality

The Charged-hadron yields as a function of pT in different y selections in p+Pb collisions.

The K^{0}_{S} yields as a function of pT in different y selections in Pb+Pb photo-nuclear collisions.

The #Lambda yields as a function of pT in different y selections in Pb+Pb photo-nuclear collisions.

The #Xi^{-} yields as a function of pT in different y selections in Pb+Pb photo-nuclear collisions.

The K^{0}_{S} yields as a function of pT in different y selections in p+Pb collisions.

The #Lambda yields as a function of pT in different y selections in p+Pb collisions.

The #Xi^{-} yields as a function of pT in different y selections in p+Pb collisions.

The Charged-hadron and identified-hadron yields as a function of y in Pb+Pb photo-nuclear collisions.

The Charged-hadron and identified-hadron yields as a function of y in Pb+Pb photo-nuclear collisions.

The Charged-hadron and identified-hadron yields as a function of y in p+Pb collisions.

The Charged-hadron and identified-hadron yields as a function of y in p+Pb collisions.

The meanpT of charged hadrons and identified hadrons as a function of Nchrec in Pb+Pb photo-nuclear collisions.

The meanpT of charged hadrons and identified hadrons as a function of Nchrec in Pb+Pb photo-nuclear collisions.

The meanpT of charged hadrons and identified hadrons as a function of Nchrec in Pb+Pb photo-nuclear collisions.

The meanpT of charged hadrons and identified hadrons as a function of Nchrec in p+Pb collisions.

The meanpT of charged hadrons and identified hadrons as a function of Nchrec in p+Pb collisions.

The baryon to meson ratio as a function of pT in Pb+Pb photo-nuclear collisions.

The baryon to meson ratio as a function of pT in Pb+Pb photo-nuclear collisions.

The baryon to meson ratio as a function of pT in p+Pb collisions.

The baryon to meson ratio as a function of pT in p+Pb collisions.

The ratios of identified-strange-hadron yield to charged-hadron yields as a function of Nchrec in Pb+Pb photo-nuclear collisions.

The ratios of identified-strange-hadron yield to charged-hadron yields as a function of Nchrec in Pb+Pb photo-nuclear collisions.

The ratios of identified-strange-hadron yield to charged-hadron yields as a function of Nchrec in Pb+Pb photo-nuclear collisions.

The ratios of identified-strange-hadron yield to charged-hadron yields as a function of Nchrec in p+Pb collisions.

The ratios of identified-strange-hadron yield to charged-hadron yields as a function of Nchrec in p+Pb collisions.

Breakdown of relative systematic uncertainties in the measured ${t\bar{t}}$ cross-section. The quoted uncertainties are obtained by repeating the fit with a group of nuisance parameters fixed to their fitted values and subtracting in quadrature the resulting total uncertainty from the uncertainty of the complete fit. However, the total uncertainty is not the quadratic sum of the grouped impacts, as this approach neglects the correlation among the different groups.

The figure shows the measurement of the top quark pair production cross-section, $\sigma_{t\bar{t}}$, in lead-lead (Pb+Pb) collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV, based on 1.9 nb$^{-1}$ of data.

Data from Figure 1, panel c, $k_{2}$ vs $N_{ch}$ for Pb+Pb collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 1, panel c, $k_{2}$ vs $N_{ch}$ for Xe+Xe collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 1, panel d, $k_{3}$ vs $N_{ch}$ for Pb+Pb collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 1, panel d, $k_{3}$ vs $N_{ch}$ for Xe+Xe collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 2, panel a, $N_{ch}k_{2}$ vs $N_{ch}$ for Pb+Pb collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 2, panel a, $N_{ch}k_{2}$ vs $N_{ch}$ for Xe+Xe collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 2, panel b, $N^{2}_{ch}k_{3}$ vs $N_{ch}$ for Pb+Pb collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 2, panel b, $N^{2}_{ch}k_{3}$ vs $N_{ch}$ for Xe+Xe collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 2, panel c, $\Gamma$ vs $N_{ch}$ for Pb+Pb collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 2, panel c, $\Gamma$ vs $N_{ch}$ for Xe+Xe collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 3, panel a, $\left\langle[p_{T}]\right\rangle / \left\langle[p_{T}]\right\rangle^{5\%}$ vs $N_{ch}/N^{5\%}_{ch}$ for Pb+Pb collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 3, panel a, $\left\langle[p_{T}]\right\rangle / \left\langle[p_{T}]\right\rangle^{5\%}$ vs $N_{ch}/N^{5\%}_{ch}$ for Xe+Xe collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 3, panel b, $k_{2}/k^{5\%}_{2} vs N_{ch}/N^{5\%}_{ch}$ for Pb+Pb collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 3, panel b, $k_{2}/k^{5\%}_{2} vs N_{ch}/N^{5\%}_{ch}$ for Xe+Xe collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 3, panel c, $k_{3}/k^{5\%}_{3} vs N_{ch}/N^{5\%}_{ch}$ for Pb+Pb collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 3, panel c, $k_{3}/k^{5\%}_{3} vs N_{ch}/N^{5\%}_{ch}$ for Xe+Xe collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 3, panel d, $\Gamma/\Gamma^{5\%}$ vs $N_{ch}/N^{5\%}_{ch}$ for Pb+Pb collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 3, panel d, $\Gamma/\Gamma^{5\%}$ vs $N_{ch}/N^{5\%}_{ch}$ for Xe+Xe collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 4, $\Delta p_{T}/\left\langle[p_{T}]\right\rangle_{0-1\%}$ vs $\Delta N_{ch}/\left\langle N_{ch}\right\rangle_{0-1\%}$ for Pb+Pb collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 4, $\Delta p_{T}/\left\langle[p_{T}]\right\rangle_{0-1\%}$ vs $\Delta N_{ch}/\left\langle N_{ch}\right\rangle_{0-1\%}$ for Xe+Xe collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Figure 4, $\Delta p_{T}/\left\langle[p_{T}]\right\rangle_{0-1\%}$ vs $\Delta N_{ch}/\left\langle N_{ch}\right\rangle_{0-1\%}$ for Pb+Pb collisions, 0.5 $ <p_{T}< $ 2 GeV/c, $|\eta|< $ 2.5

Data from Figure 4, $\Delta p_{T}/\left\langle[p_{T}]\right\rangle_{0-1\%}$ vs $\Delta N_{ch}/\left\langle N_{ch}\right\rangle_{0-1\%}$ for Xe+Xe collisions, 0.5 $ <p_{T}< $ 2 GeV/c, $|\eta|< $ 2.5

Data from Second Supplementary Figure, Fig 6, panel a, Correlation between FCal $\Sigma E_{T}$ and $N_{ch}$ in Pb+Pb collisions

Data from Second Supplementary Figure, Fig 6, panel b, Correlation between $[p_{T}]$ - $\left\langle[p_{T}]\right\rangle_{0-1\%}$ vs $\Delta N_{ch}$/$\left\langle N_{ch}\right\rangle_{0-1\%}$ in Pb+Pb collisions

Data from Third Supplementary Figure, Fig 7, Normalized Event distribution of $N_{ch}$ /$N^{5\%}_{ch}$ for Pb+Pb collisions

Data from Third Supplementary Figure, Fig 7, Normalized Event distribution of $N_{ch}$ / $N^{5\%}_{ch}$ for Xe+Xe collisions

Data from Fourth Supplementary Figure, Fig 8, panel a, $\gamma$ vs $N_{ch}$ for Pb+Pb collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Fourth Supplementary Figure, Fig 8, panel a, $\gamma$ vs $N_{ch}$ for Xe+Xe collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Fourth Supplementary Figure, Fig 8, panel b, $\gamma/\gamma^{5\%}$ vs $N_{ch}/N^{5\%}_{ch}$ for Pb+Pb collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Fourth Supplementary Figure, Fig 8, panel b, $\gamma/\gamma^{5\%}$ vs $N_{ch}/N^{5\%}_{ch}$ for Xe+Xe collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Auxiliary Fig 1 panel a, Correlation between FCal $\Sigma E_{T}$ and $N_{ch}$ in Xe+Xe collisions

Data from Auxiliary Fig 1 panel b, Correlation between $[p_{T}]$ - $\left\langle[p_{T}]\right\rangle_{0-1\%}$ vs $\Delta N_{ch}$/$\left\langle N_{ch}\right\rangle_{0-1\%}$ in Xe+Xe collisions

Data from Auxiliary Fig 2 panel b, $N_{ch}$ /$N^{5\%}_{ch}$ vs Centrality [\%] for Pb+Pb collisions

Data from Auxiliary Fig 2 panel b, $N_{ch}$ /$N^{5\%}_{ch}$ vs Centrality [\%] for Xe+Xe collisions

Data from Auxiliary Fig 3 panel b, $[p_{T}]$ vs $N_{ch}$ for Xe+Xe collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Auxiliary Figure 4, $\Delta p_{T}/\left\langle[p_{T}]\right\rangle_{0-0.5\%}$ vs $\Delta N_{ch}/\left\langle N_{ch}\right\rangle_{0-0.5\%}$ for Pb+Pb collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Auxiliary Figure 4, $\Delta p_{T}/\left\langle[p_{T}]\right\rangle_{0-0.5\%}$ vs $\Delta N_{ch}/\left\langle N_{ch}\right\rangle_{0-0.5\%}$ for Xe+Xe collisions, 0.5 $ <p_{T}< $ 5 GeV/c, $|\eta|< $ 2.5

Data from Auxiliary Figure 4, $\Delta p_{T}/\left\langle[p_{T}]\right\rangle_{0-0.5\%}$ vs $\Delta N_{ch}/\left\langle N_{ch}\right\rangle_{0-0.5\%}$ for Pb+Pb collisions, 0.5 $ <p_{T}< $ 2 GeV/c, $|\eta|< $ 2.5

Data from Auxiliary Figure 4, $\Delta p_{T}/\left\langle[p_{T}]\right\rangle_{0-0.5\%}$ vs $\Delta N_{ch}/\left\langle N_{ch}\right\rangle_{0-0.5\%}$ for Xe+Xe collisions, 0.5 $ <p_{T}< $ 2 GeV/c, $|\eta|< $ 2.5

The ratios between the hydrodynamic probes and their Taylor expansions. The ratios are plotted as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

The ratios between the hydrodynamic probes and their Taylor expansions. The ratios are plotted as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

The magnitudes of the measured skewness $\gamma_1^\textrm{exp}$ and "corrected" skewness $\gamma_{1,\textrm{corr}}^\textrm{exp}$ as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

The magnitudes of the measured kurtosis $\gamma_2^\textrm{exp}$ and "corrected" kurtosis $\gamma_{2,\textrm{corr}}^\textrm{exp}$ as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

The magnitudes of the measured superskewness $\gamma_3^\textrm{exp}$ and "corrected" superskewness $\gamma_{3,\textrm{corr}}^\textrm{exp}$ as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

The $v_{3}$ values as a function of centrality for prompt and nonprompt J/$\psi$ mesons.The kinematic range is 6.5 $< p_{\text{T}} <$ 50 GeV/c and $|y| < 2.4$.

The $v_{2}$ values of prompt J/$\psi$ mesons as a function of $p_{\text{T}}$ in the 10–60% centrality range. The results for 3 $< p_{\text{T}} <$ 6.5 and 6.5 $< p_{\text{T}} <$ 50 GeV/c are studied in the rapidity range of 1.6 $< |y| <$ 2.4 and $|y| <$ 2.4, respectively.

The $v_{2}$ values as functions of $p_{\text{T}}$ for prompt J/$\psi$ and prompt $\psi$(2S) mesons in the 10–60% centrality range. The results for $p_{\text{T}}$ < 6.5 and $p_{\text{T}}$ > 6.5 GeV/c are studied in the rapidity range of 1.6 $< |y| <$ 2.4 and $|y| <$ 2.4, respectively.

The $v_{2}$ values as functions of centrality for prompt J/$\psi$ and prompt $\psi$(2S) mesons. The kinematic range is 6.5 $< p_{\text{T}} <$ 50 GeV/c and $|y| <$ 2.4.

The $v_{3}$ values as functions of $p_{\text{T}}$ for prompt J/$\psi$ and prompt $\psi$(2S) mesons in the 10–60% centrality range. The results for $p_{\text{T}}$ < 6.5 and $p_{\text{T}}$ > 6.5 GeV/c are studied in the rapidity range of 1.6 $< |y| <$ 2.4 and $|y| <$ 2.4, respectively.

The $v_{3}$ values as functions of centrality for prompt J/$\psi$ and prompt $\psi$(2S) mesons. The kinematic range is 6.5 $< p_{\text{T}} <$ 50 GeV/c and $|y| <$ 2.4.

Measured $R_\text{AA}$ for the $\Upsilon(3S)$ state as functions of PbPb collision centrality, integrated over the full kinematic range $p_\text{T}$ < 30 GeV/c and |y| < 2.4. The global uncertainty "PP MB" represents the pp luminosity and PbPb $N_\text{MB}$ combined uncertainties, whereas the global uncertainty "PP 3S" corresponds to the uncertainty on the $\Upsilon(3S)$ pp yields.

Measured $R_\text{AA}$ for the$\Upsilon(2S)$ state in the 0–90% centrality interval, integrated over the full kinematic range $p_\text{T}$ < 30 GeV/c and |y| < 2.4. The global uncertainty "PP MB" represents the pp luminosity and PbPb $N_\text{MB}$ combined uncertainties, whereas the global uncertainty "PP 2S" corresponds to the uncertainty on the $\Upsilon(2S)$ pp yields.

Measured $R_\text{AA}$ for the$\Upsilon(2S)$ state in the 0–90% centrality interval, integrated over the full kinematic range $p_\text{T}$ < 30 GeV/c and |y| < 2.4. The global uncertainty "PP MB" represents the pp luminosity and PbPb $N_\text{MB}$ combined uncertainties, whereas the global uncertainty "PP 2S" corresponds to the uncertainty on the $\Upsilon(2S)$ pp yields.

Measured $R_\text{AA}$ for the $\Upsilon(3S)$ state in the 0–90% centrality interval, integrated over the full kinematic range $p_\text{T}$ < 30 GeV/c and |y| < 2.4. The global uncertainty "PP MB" represents the pp luminosity and PbPb $N_\text{MB}$ combined uncertainties, whereas the global uncertainty "PP 3S" corresponds to the uncertainty on the$\Upsilon(3S)$ pp yields.

Measured $R_\text{AA}$ for the $\Upsilon(3S)$ state in the 0–90% centrality interval, integrated over the full kinematic range $p_\text{T}$ < 30 GeV/c and |y| < 2.4. The global uncertainty "PP MB" represents the pp luminosity and PbPb $N_\text{MB}$ combined uncertainties, whereas the global uncertainty "PP 3S" corresponds to the uncertainty on the$\Upsilon(3S)$ pp yields.

Measured $R_\text{AA}$ for the $\Upsilon(2S)$ state as functions of transverse momentum $p_\text{T}$, integrated over the rapidity range of |y| < 2.4 and the 0–90% centrality interval. The global uncertainty combines the uncertainties of $T_\text{AA}$ , pp luminosity, and PbPb $N_\text{MB}$.

Measured $R_\text{AA}$ for the $\Upsilon(2S)$ state as functions of transverse momentum $p_\text{T}$, integrated over the rapidity range of |y| < 2.4 and the 0–90% centrality interval. The global uncertainty combines the uncertainties of $T_\text{AA}$ , pp luminosity, and PbPb $N_\text{MB}$.

Measured $R_\text{AA}$ for the $\Upsilon(3S)$ state as functions of transverse momentum $p_\text{T}$, integrated over the rapidity range of |y| < 2.4 and the 0–90% centrality interval. The global uncertainty combines the uncertainties of $T_\text{AA}$ , pp luminosity, and PbPb $N_\text{MB}$.

Measured $R_\text{AA}$ for the $\Upsilon(3S)$ state as functions of transverse momentum $p_\text{T}$, integrated over the rapidity range of |y| < 2.4 and the 0–90% centrality interval. The global uncertainty combines the uncertainties of $T_\text{AA}$ , pp luminosity, and PbPb $N_\text{MB}$.

The double ratios of $\Upsilon(3S)$ /$\Upsilon(2S)$ as functions of PbPb collision centrality, integrated over the full kinematic range $p_\text{T}$ < 30 GeV/c and |y| < 2.4. The global uncertainty represents the combined systematic and statistical uncertainties from pp data.

The double ratios of $\Upsilon(3S)$ /$\Upsilon(2S)$ as functions of PbPb collision centrality, integrated over the full kinematic range $p_\text{T}$ < 30 GeV/c and |y| < 2.4. The global uncertainty represents the combined systematic and statistical uncertainties from pp data.

The double ratios of $\Upsilon(3S)$ /$\Upsilon(2S)$ in the 0–90% centrality interval, integrated over the full kinematic range $p_\text{T}$ < 30 GeV/c and |y| < 2.4. The global uncertainty represents the combined systematic and statistical uncertainties from pp data.

The double ratios of $\Upsilon(3S)$ /$\Upsilon(2S)$ in the 0–90% centrality interval, integrated over the full kinematic range $p_\text{T}$ < 30 GeV/c and |y| < 2.4. The global uncertainty represents the combined systematic and statistical uncertainties from pp data.

The double ratios of $\Upsilon(3S)$ /$\Upsilon(2S)$ as functions of transverse momentum $p_\text{T}$, integrated over the rapidity range of |y| < 2.4 and the 0–90% centrality interval.

The double ratios of $\Upsilon(3S)$ /$\Upsilon(2S)$ as functions of transverse momentum $p_\text{T}$, integrated over the rapidity range of |y| < 2.4 and the 0–90% centrality interval.

Yields for $\Upsilon(2S)$ mesons in PbPb collisions in centrality $0--90\%$ and $|y| < 2.4$, corrected for acceptance and efficiency, and normalized by the nuclear thickness function $\langle T_{AA} \rangle$ and the number of minimum bias events $N_{MB}$. The numbers are in units of pb$/GeV/c$.

Yields for $\Upsilon(2S)$ mesons in PbPb collisions in $p_\text{T} < 30GeV/c$ and $|y| < 2.4$, corrected for acceptance and efficiency, and normalized by the nuclear thickness function $\langle T_{AA} \rangle$ and the number of minimum bias events $N_{MB}$. The numbers are in units of pb$/GeV/c$.

Yields for $\Upsilon(3S)$ mesons in PbPb collisions in centrality $0--90\%$ and $|y| < 2.4$, corrected for acceptance and efficiency, and normalized by the nuclear thickness function $\langle T_{AA} \rangle$ and the number of minimum bias events $N_{MB}$. The numbers are in units of pb$/GeV/c$.

Yields for $\Upsilon(3S)$ mesons in PbPb collisions in $p_\text{T} < 30GeV/c$ and $|y| < 2.4$, corrected for acceptance and efficiency, and normalized by the nuclear thickness function $\langle T_{AA} \rangle$ and the number of minimum bias events $N_{MB}$. The numbers are in units of pb$/GeV/c$.

Yields for $\Upsilon(1S)$ mesons in PbPb collisions in $p_\text{T} < 30GeV/c$ and $|y| < 2.4$, corrected for acceptance and efficiency, and normalized by the nuclear thickness function $\langle T_{AA} \rangle$ and the number of minimum bias events $N_{MB}$. The numbers are in units of pb$/GeV/c$. The values and global integrated luminosity uncertainty are quoted from data taken from CMS in 2018 LHC PbPb run.