The JEC, JER, and total systematic uncertainties in unfolded cross sections as functions of transverse momentum, for 1.5<|y|<2. The total systematic uncertainty includes also the luminosity, jet identification and trigger efficiency uncertainties.

The unfolded measured particle-level inclusive jet cross section as functions of jet pT (markers), for |y|<0.5, compared to the NLO perturbative QCD prediction (histogram), using the CT14NLO PDF set, with muR = muF = HT, and corrected for the NP effects. The experimental and theoretical systematic uncertainties are shown.

The unfolded measured particle-level inclusive jet cross section as functions of jet pT (markers), for 0.5<|y|<1, compared to the NLO perturbative QCD prediction (histogram), using the CT14NLO PDF set, with muR = muF = HT, and corrected for the NP effects. The experimental and theoretical systematic uncertainties are shown.

The unfolded measured particle-level inclusive jet cross section as functions of jet pT (markers), for 1<|y|<1.5, compared to the NLO perturbative QCD prediction (histogram), using the CT14NLO PDF set, with muR = muF = HT, and corrected for the NP effects. The experimental and theoretical systematic uncertainties are shown.

The unfolded measured particle-level inclusive jet cross section as functions of jet pT (markers), for 1.5<|y|<2, compared to the NLO perturbative QCD prediction (histogram), using the CT14NLO PDF set, with muR = muF = HT, and corrected for the NP effects. The experimental and theoretical systematic uncertainties are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = pT, for |y|<0.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = pT, for 0.5<|y|<1. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = pT, for 1<|y|<1.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = pT, for 1.5<|y|<2. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = HT, for |y|<0.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = HT, for 0.5<|y|<1. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = HT, for 1<|y|<1.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = HT, for 1.5<|y|<2. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the CT14NNLO PDF set, with mu = HT, for |y|<0.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the CT14NNLO PDF set, with mu = HT, for 0.5<|y|<1. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the CT14NNLO PDF set, with mu = HT, for 1<|y|<1.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the CT14NNLO PDF set, with mu = HT, for 1.5<|y|<2. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the NNPDF31NNLO PDF set, with mu = HT, for |y|<0.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the NNPDF31NNLO PDF set, with mu = HT, for 0.5<|y|<1. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the NNPDF31NNLO PDF set, with mu = HT, for 1<|y|<1.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the NNPDF31NNLO PDF set, with mu = HT, for 1.5<|y|<2. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

The effect of aS(MZ) variation, for |y|<0.5. The NNLO theoretical cross section predictions using the NNPDF31NNLO PDF with m = HT, calculated for different choices of aS (0.108, 0.110, 0.112, 0.114, 0.116, 0.117, 0.118, 0.119, 0.120, 0.122, and 0.124), are divided by the benchmark NNLO prediction for aS = 0.118 and the same choice of PDF set, muR, and muF. Also shown is the experimental unfolded measurement divided by the same benchmark prediction. The width of the unity line corresponds to the statistical uncertainty from the MC integration for the determination of the NNLO prediction. The error bars on the unfolded data correspond to the total experimental statistical and systematic uncertainty added in quadrature.

The effect of aS(MZ) variation, for 0.5<|y|<1. The NNLO theoretical cross section predictions using the NNPDF31NNLO PDF with m = HT, calculated for different choices of aS (0.108, 0.110, 0.112, 0.114, 0.116, 0.117, 0.118, 0.119, 0.120, 0.122, and 0.124), are divided by the benchmark NNLO prediction for aS = 0.118 and the same choice of PDF set, muR, and muF. Also shown is the experimental unfolded measurement divided by the same benchmark prediction. The width of the unity line corresponds to the statistical uncertainty from the MC integration for the determination of the NNLO prediction. The error bars on the unfolded data correspond to the total experimental statistical and systematic uncertainty added in quadrature.

The effect of aS(MZ) variation, for 1<|y|<1.5. The NNLO theoretical cross section predictions using the NNPDF31NNLO PDF with m = HT, calculated for different choices of aS (0.108, 0.110, 0.112, 0.114, 0.116, 0.117, 0.118, 0.119, 0.120, 0.122, and 0.124), are divided by the benchmark NNLO prediction for aS = 0.118 and the same choice of PDF set, muR, and muF. Also shown is the experimental unfolded measurement divided by the same benchmark prediction. The width of the unity line corresponds to the statistical uncertainty from the MC integration for the determination of the NNLO prediction. The error bars on the unfolded data correspond to the total experimental statistical and systematic uncertainty added in quadrature.

The effect of aS(MZ) variation, for 1.5<|y|<2. The NNLO theoretical cross section predictions using the NNPDF31NNLO PDF with m = HT, calculated for different choices of aS (0.108, 0.110, 0.112, 0.114, 0.116, 0.117, 0.118, 0.119, 0.120, 0.122, and 0.124), are divided by the benchmark NNLO prediction for aS = 0.118 and the same choice of PDF set, muR, and muF. Also shown is the experimental unfolded measurement divided by the same benchmark prediction. The width of the unity line corresponds to the statistical uncertainty from the MC integration for the determination of the NNLO prediction. The error bars on the unfolded data correspond to the total experimental statistical and systematic uncertainty added in quadrature.

Unfolded correlation matrix for |y|<0.5.

Unfolded correlation matrix for 0.5<|y|<1.

Unfolded correlation matrix for 1<|y|<1.5.

Unfolded correlation matrix for 1.5<|y|<2.

$\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}$

Near-side ratios of V^{0}-triggered yields to charged-hadron triggered (h$-$h) ones, ${{Y}^{K^{0}_{S}-\mathrm{h}}_{\Delta\varphi}}/{{Y}^{\mathrm{h-h}}_{\Delta\varphi}}$ and ${{Y}^{(\Lambda+\overline{\Lambda})-\mathrm{h}}_{\Delta\varphi}}/{{Y}^{\mathrm{h-h}}_{\Delta\varphi}}$ in Pb-Pb collisions.

Away-side ratios of V^{0}-triggered yields to charged-hadron triggered (h$-$h) ones, ${{Y}^{K^{0}_{S}-\mathrm{h}}_{\Delta\varphi}}/{{Y}^{\mathrm{h-h}}_{\Delta\varphi}}$ and ${{Y}^{(\Lambda+\overline{\Lambda})-\mathrm{h}}_{\Delta\varphi}}/{{Y}^{\mathrm{h-h}}_{\Delta\varphi}}$ in Pb-Pb collisions.

Near-side ratios of V^{0}-triggered yields to charged-hadron triggered (h--h) ones, ${{Y}^{K^{0}_{S}-\mathrm{h}}_{\Delta\varphi}}/{{Y}^{\mathrm{h-h}}_{\Delta\varphi}}$ and ${{Y}^{(\Lambda+\overline{\Lambda})-\mathrm{h}}_{\Delta\varphi}}/{{Y}^{\mathrm{h-h}}_{\Delta\varphi}}$ in pp collisions.

Away-side ratios of V^{0}-triggered yields to charged-hadron triggered (h--h) ones, ${{Y}^{K^{0}_{S}-\mathrm{h}}_{\Delta\varphi}}/{{Y}^{\mathrm{h-h}}_{\Delta\varphi}}$ and ${{Y}^{(\Lambda+\overline{\Lambda})-\mathrm{h}}_{\Delta\varphi}}/{{Y}^{\mathrm{h-h}}_{\Delta\varphi}}$ in Pb-Pb collisions.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with 400 < $p_{T1}$ < 800 GeV and for an azimuthal separation between the two leading jets of $0 < \Delta\Phi_{1,2} < 150^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with 400 < $p_{T1}$ < 800 GeV and for an azimuthal separation between the two leading jets of $150 < \Delta\Phi_{1,2} < 170^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with 400 < $p_{T1}$ < 800 GeV and for an azimuthal separation between the two leading jets of $170 < \Delta\Phi_{1,2} < 180^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with $p_{T1}$ > 800 and for an azimuthal separation between the two leading jets of $0 < \Delta\Phi_{1,2} < 150^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with $p_{T1}$ > 800 and for an azimuthal separation between the two leading jets of $150 < \Delta\Phi_{1,2} < 170^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with $p_{T1}$ > 800 and for an azimuthal separation between the two leading jets of $170 < \Delta\Phi_{1,2} < 180^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Measured transverse momentum of the leading $p_{T}$ jet ($p_{T1}$). The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Measured transverse momentum of the subleading $p_{T}$ jet ($p_{T2}$). The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Measured transverse momentum of the third leading $p_{T}$ jet ($p_{T3}$). The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Measured transverse momentum of the fourth leading $p_{T}$ jet ($p_{T4}$). The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity distribution in TUnfold binning ($N_{jets},\Delta\phi_{1,2},p_{T1}$) as indicated in the XML file provided as additional resource. The uncertainties follow the notation of Table 1.

Correlation matrix at particle level for the measured jet multiplicity in TUnfold binning ($N_{jets},\Delta\phi_{1,2},p_{T1}$) as indicated in the XML file provided as additional resource.

Jet $p_{T}$ distributions in TUnfold binning ($p_{T1},p_{T2},p_{T3},p_{T4}$) as indicated in the XML file provided as additional resource. The uncertainties follow the notation of Table 1.

Correlation matrix at particle level for the measured jet $p_{T}$ distributions in TUnfold binning ($p_{T1},p_{T2},p_{T3},p_{T4}$) as indicated in the XML file provided as additional resource.

Inclusive muon $v_{2}^{\mu}$ as a function of $p_{\mathrm{T}}$ is measured by two-particle cumulant method at backward rapidities in high-multiplicity (0$-$20%) p$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$= 8.16 TeV. The event activity is estimated with the V0M estimator.

Inclusive muon $v_{2}^{\mu}$ as a function of $p_{\mathrm{T}}$ is measured by two-particle cumulant method at forward rapidities in high-multiplicity (0$-$20%) p$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$= 8.16 TeV. The event activity is estimated with the CL1 estimator.

Inclusive muon $v_{2}^{\mu}$ as a function of $p_{\mathrm{T}}$ is measured by two-particle cumulant method at forward rapidities in high-multiplicity (0$-$20%) p$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$= 8.16 TeV. The event activity is estimated with the ZDC estimator.

Inclusive muon $v_{2}^{\mu}$ as a function of $p_{\mathrm{T}}$ is measured by two-particle cumulant method at backward rapidities in high-multiplicity (0$-$20%) p$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$= 8.16 TeV. The event activity is estimated with the CL1 estimator.

Inclusive muon $v_{2}^{\mu}$ as a function of $p_{\mathrm{T}}$ is measured by two-particle cumulant method at backward rapidities in high-multiplicity (0$-$20%) p$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$= 8.16 TeV. The event activity is estimated with the ZDC estimator.

Inclusive muon $v_{2}^{\mu}$ as a function of $p_{\mathrm{T}}$ is measured by two-particle correlation method at forward rapidities in high-multiplicity (0$-$20%) p$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$= 8.16 TeV. The event activity is estimated with the CL1 estimator.

Inclusive muon $v_{2}^{\mu}$ as a function of $p_{\mathrm{T}}$ is measured by two-particle correlation method at forward rapidities in high-multiplicity (0$-$20%) p$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$= 8.16 TeV. The event activity is estimated with the ZDC estimator.

Inclusive muon $v_{2}^{\mu}$ as a function of $p_{\mathrm{T}}$ is measured by two-particle corrrelation method at backward rapidities in high-multiplicity (0$-$20%) p$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$= 8.16 TeV. The event activity is estimated with the CL1 estimator.

Inclusive muon $v_{2}^{\mu}$ as a function of $p_{\mathrm{T}}$ is measured by two-particle correlation method at backward rapidities in high-multiplicity (0$-$20%) p$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$= 8.16 TeV. The event activity is estimated with the ZDC estimator.

Negative log-likelihood difference in $\text{c}_{\varphi\text{Q}}^{3}$ where the other Wilson coefficients are fixed to 0.

Negative log-likelihood difference in $\text{c}_{\varphi\text{t}}$ where the other Wilson coefficients are fixed to 0.

Negative log-likelihood difference in $\text{c}_{\varphi\text{tb}}$ where the other Wilson coefficients are fixed to 0.

Negative log-likelihood difference in $\text{c}_{\text{tW}}$ where the other Wilson coefficients are fixed to 0.

Negative log-likelihood difference in $\text{c}_{\text{bW}}$ where the other Wilson coefficients are fixed to 0.

Negative log-likelihood difference in $\text{c}_{\text{tZ}}$ where the other Wilson coefficients are fixed to 0.

Negative log-likelihood difference in $\text{c}_{\text{t}\varphi}$ where the other Wilson coefficients are simultaneously fit.

Negative log-likelihood difference in $\text{c}_{\varphi\text{Q}}^{-}$ where the other Wilson coefficients are simultaneously fit.

Negative log-likelihood difference in $\text{c}_{\varphi\text{Q}}^{3}$ where the other Wilson coefficients are simultaneously fit.

Negative log-likelihood difference in $\text{c}_{\varphi\text{t}}$ where the other Wilson coefficients are simultaneously fit.

Negative log-likelihood difference in $\text{c}_{\varphi\text{tb}}$ where the other Wilson coefficients are simultaneously fit.

Negative log-likelihood difference in $\text{c}_{\text{tW}}$ where the other Wilson coefficients are simultaneously fit.

Negative log-likelihood difference in $\text{c}_{\text{bW}}$ where the other Wilson coefficients are simultaneously fit.

Negative log-likelihood difference in $\text{c}_{\text{tZ}}$ where the other Wilson coefficients are simultaneously fit.

Negative log-likelihood difference in $\text{c}_{\varphi\text{t}}, \text{c}_{\varphi\text{Q}}^{-}$ where the other Wilson coefficients are fixed to 0.

Negative log-likelihood difference in $\text{c}_{\varphi\text{Q}}^{3}, \text{c}_{\varphi\text{Q}}^{-}$ where the other Wilson coefficients are fixed to 0.

Negative log-likelihood difference in $\text{c}_{\text{tW}}, \text{c}_{\text{tZ}}$ where the other Wilson coefficients are fixed to 0.

$f_{0}(980)$ to pions ($(\pi^{+}+\pi^{-})/2$), $p_{\rm T}$-integrated yield ratio at midrapidity as a function of $\langle {\rm d}N_{\rm ch}/{\rm d}\eta \rangle$. The uncertainty 'stat,syst' indicates the total uncertainty (stat+syst) on the measurement. The branching ratio correction amounts to BR = (46 $\pm$ 6)% [ Phys. Rev. Lett. 111 no. 6, (2013) 062001] assuming dominance of $\pi\pi$ and KK channel has been applied to the $p_{\rm T}$-differential yields of $f_{0}(980)$. Here, the branching ratio relative uncertainty (13%) for $f_{0}(980)$ is not included.

Ratio of measured PSI(2S) over J/PSI(1S) cross section in charged-particle multiplicity intervals and integrated in multiplicity.

Ratio of measured PSI(2S) over J/PSI(1S) cross section in charged-particle multiplicity intervals and integrated in multiplicity.

Ratio of measured PSI(2S) over J/PSI(1S) cross section in charged-particle multiplicity intervals and integrated in multiplicity.

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($-0.9 < y < -0.6$).

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($-1.2 < y < -0.9$).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($0.0 < y < 0.3$).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($-0.3 < y < 0.0$).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($-0.6 < y < -0.3$).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($-0.9 < y < -0.6$).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV ($-1.2 < y < -0.9$).

d$N$/d$y$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (V0A multiplicity classes).

d$N$/d$y$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (V0A multiplicity classes).

$\langle p_\mathrm{T}\rangle$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (V0A multiplicity classes).

$\langle p_\mathrm{T}\rangle$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (V0A multiplicity classes).

(d$N$/$\mathrm{d}y$)/(d$N$/$\mathrm{d}y)_{y=0}$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (0-100 %).

(d$N$/$\mathrm{d}y$)/(d$N$/$\mathrm{d}y)_{y=0}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (0-100 %).

$\langle p_\mathrm{T}\rangle_\mathrm{y}$/$\langle p_\mathrm{T}\rangle_\mathrm{y=0}$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (0-100 %).

$\langle p_\mathrm{T}\rangle_\mathrm{y}$/$\langle p_\mathrm{T}\rangle_\mathrm{y=0}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (0-100%).

$y_{asym}$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (V0A multiplicity classes).

$y_{asym}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (V0A multiplicity classes).

$Q_{CP}$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (0-10 %/40-100 %).

$Q_{CP}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (0-10 %/40-100%).

$Q_{CP}$ of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (10-40 %/40-100 %).

$Q_{CP}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV (10-40 %/40-100%).