Summary of results for cross-section of non-prompt $\psi(2S)$ decaying to a muon pair for 13 TeV data in fb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $J/\psi$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $\psi(2S)$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Correction factors for an assumed angular dependence $\propto 1+\lambda_{\theta} \cos^2\theta $ in the helicity frame, for several values of the parameter $\lambda_{\theta}$. The correction factors were found to be the same (within 1% - 2%) for $J\psi$ and $\psi(2S)$ mesons, for prompt and non-prompt production mechanisms, and for the three rapidity intervals considered in this paper.

Left pannel. The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 2.3%) is not graphically represented. The blue line represents the average for $p_T > 18$ GeV. For comparison, the LHCb measurement [10.1103/PhysRevLett.124.122002] is also shown. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (2.3%)

Left pannel. The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 2.3%) is not graphically represented. The blue line represents the average for $p_T > 18$ GeV. For comparison, the LHCb measurement [10.1103/PhysRevLett.124.122002] is also shown. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (2.3%)

Left pannel. The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 2.3%) is not graphically represented. The blue line represents the average for $p_T > 18$ GeV. For comparison, the LHCb measurement [10.1103/PhysRevLett.124.122002] is also shown. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (2.3%)

Right pannel. The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 2.3%) is not graphically represented. The blue line represents the average for $p_T > 18$ GeV. For comparison, the LHCb measurement [10.1103/PhysRevLett.124.122002] is also shown. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (2.3%)

Right pannel. The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 2.3%) is not graphically represented. The blue line represents the average for $p_T > 18$ GeV. For comparison, the LHCb measurement [10.1103/PhysRevLett.124.122002] is also shown. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (2.3%)

Right pannel. The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 2.3%) is not graphically represented. The blue line represents the average for $p_T > 18$ GeV. For comparison, the LHCb measurement [10.1103/PhysRevLett.124.122002] is also shown. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (2.3%)

Left panel.The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 4.9\%) is not graphically represented. The blue line represents the average of all the points. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (4.9%)

Left panel.The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 5.7%) is not graphically represented. The blue line represents the average of all the points. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (5.7%)

Left panel.The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 5.7%) is not graphically represented. The blue line represents the average of all the points. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (5.7%)

Right panel.The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 4.9\%) is not graphically represented. The blue line represents the average of all the points. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (4.9%)

Right panel.The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 5.7%) is not graphically represented. The blue line represents the average of all the points. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (5.7%)

Right panel.The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 5.7%) is not graphically represented. The blue line represents the average of all the points. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (5.7%)

Charged-hadron spectrum in the centrality interval 5-10% for p+Pb, divided by &#9001;TPPB&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 10-20% for p+Pb, divided by &#9001;TPPB&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 20-30% for p+Pb, divided by &#9001;TPPB&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 30-40% for p+Pb, divided by &#9001;TPPB&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 40-60% for p+Pb, divided by &#9001;TPPB&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 60-90% for p+Pb, divided by &#9001;TPPB&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 0-90% for p+Pb, divided by &#9001;TPPB&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron cross-section in pp collisions. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 0-5% for Pb+Pb, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Charged-hadron spectrum in the centrality interval 5-10% for Pb+Pb, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Charged-hadron spectrum in the centrality interval 10-20% for Pb+Pb, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Charged-hadron spectrum in the centrality interval 20-30% for Pb+Pb, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Charged-hadron spectrum in the centrality interval 30-40% for Pb+Pb, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Charged-hadron spectrum in the centrality interval 40-50% for Pb+Pb, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Charged-hadron spectrum in the centrality interval 50-60% for Pb+Pb, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Charged-hadron spectrum in the centrality interval 60-80% for Pb+Pb, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Charged-hadron cross-section in pp collisions. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 0-5% for Xe+Xe, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 5-10% for Xe+Xe, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 10-20% for Xe+Xe, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 20-30% for Xe+Xe, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 30-40% for Xe+Xe, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 40-50% for Xe+Xe, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 50-60% for Xe+Xe, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Charged-hadron spectrum in the centrality interval 60-80% for Xe+Xe, divided by &#9001;TAA&#9002;. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-60% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Nuclear modification factor in centrality interval 5-10% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Nuclear modification factor in centrality interval 10-20% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Nuclear modification factor in centrality interval 20-30% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Nuclear modification factor in centrality interval 30-40% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Nuclear modification factor in centrality interval 40-50% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Nuclear modification factor in centrality interval 50-60% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Nuclear modification factor in centrality interval 60-80% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature. The systematic uncertainty on momentum bias is negligible at low pT; in such cases, it is omitted in the table below.

Nuclear modification factor in centrality interval 0-5% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-50% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 50-60% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-80% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-60% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-60% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-60% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-60% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-90% for p+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-50% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-50% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-50% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-50% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 50-60% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 50-60% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 50-60% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 50-60% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-80% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-80% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-80% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-80% for Pb+Pb. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 0-5% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 5-10% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 10-20% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 20-30% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 30-40% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-50% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-50% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 40-50% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 50-60% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 50-60% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 50-60% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-80% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-80% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Nuclear modification factor in centrality interval 60-80% for Xe+Xe. The systematic uncertainties are described in the section 7 of the paper. The total systematic uncertainties are determined by adding the contributions from all relevant sources in quadrature.

Final state radiation correction used in the $\phi_{\eta}^{*}$ cross-section measurement. The first uncertainty is statistical and the second is systematic.

Final state radiation correction used in the $y^{Z}-p_{T}^{Z}$ cross-section measurement. The first uncertainty is statistical and the second is systematic.

Final state radiation correction used in the $y^{Z}-\phi_{\eta}^{*}$ cross-section measurement. The first uncertainty is statistical and the second is systematic.

Correlation matrix of statistical uncertainty for one-dimensional $y^Z$ measurement.

Correlation matrix of statistical uncertainty for one-dimensional $p_{T}^{Z}$ measurement.

Correlation matrix of statistical uncertainty for one-dimensional $\phi_{\eta}^{*}$ measurement.

Correlation matrix of statistical uncertainty for two-dimensional $y^Z-p_{T}^{Z}$ measurement.

Correlation matrix of statistical uncertainty for two-dimensional $y^Z-\phi_{\eta}^{*}$ measurement.

Correlation matrix of efficiency uncertainty for one-dimensional $y^Z$ measurement.

Correlation matrix of efficiency uncertainty for one-dimensional $p_{T}^{Z}$ measurement.

Correlation matrix of efficiency uncertainty for one-dimensional $\phi_{\eta}^{*}$ measurement.

Correlation matrix of efficiency uncertainty for two-dimensional $y^Z-p_{T}^{Z}$ measurement.

Correlation matrix of efficiency uncertainty for two-dimensional $y^Z-\phi_{\eta}^{*}$ measurement.

Measured total $Z$-boson cross-section for different datasets. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured single differential cross-sections in interval regions of $y^{Z}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured single differential cross-sections in interval regions of $p_{T}^{Z}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured single differential cross-sections in interval regions of $\phi_{\eta}^{*}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured double differential cross-sections in interval regions of $y^{Z}$ and $p_{T}^{Z}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Measured double differential cross-sections in interval regions of $y^{Z}$ and $\phi_{\eta}^{*}$. The first uncertainty is statistical, the second systematic, and the third is due to the luminosity.

Systematic uncertainties in the single differential cross-sections in interval regions of $y^{Z}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

Systematic uncertainties in the single differential cross-sections in interval regions of $p_{T}^{Z}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

Systematic uncertainties in the single differential cross-sections in interval regions of $\phi_{\eta}^{*}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

Systematic uncertainties in the double differential cross-sections in interval regions of $y^{Z}$ and $p_{T}^{Z}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

Systematic uncertainties in the double differential cross-sections in interval regions of $y^{Z}$ and $\phi_{\eta}^{*}$, presented in percentage. The contributions from efficiency (Eff), background (BKG), final state radiation (FSR), closure test (Closure), and alignment and calibration (Alignment) are shown.

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.

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].

Measured rapidity differential cross section of coherent Psi(2S) photoproduction in Pb–Pb UPCs events.

Nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of the average number of participants $\langle N_{\mathrm{part}} \rangle$ or as a function of the collision centrality.

Nuclear modification factor of $\Upsilon(2\mathrm{S})$ as a function of the average number of participants $\langle N_{\mathrm{part}} \rangle$ or as a function of the collision centrality.

Ratio of $\Upsilon(2\mathrm{S})$ and $\Upsilon(1\mathrm{S})$ yields (cf. equation 2 in the article for the definition) as a function of the average number of participants $\langle N_{\mathrm{part}} \rangle$ or as a function of the collision centrality. The global uncertainty is the quadratic sum of the branching ratio uncertainties.

Relative nuclear modification factor or double yield ratio between $\Upsilon(2\mathrm{S})$ and $\Upsilon(1\mathrm{S})$ as a function of the average number of participants $\langle N_{\mathrm{part}} \rangle$ or as a function of the collision centrality. The global uncertainty corresponds to the systematic uncertainty on the cross-section ratio in proton–proton collisions.

Nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of transverse momentum for the 0–90% centrality interval.

Nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of rapidity for the 0–90% centrality interval.

Nuclear modification factor of $\Upsilon(2\mathrm{S})$ as a function of rapidity for the 0–90% centrality interval.

Interpolated production cross sections and cross-section ratio for $\Upsilon(1\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ and $\Upsilon(2\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ in proton–proton collisions, integrated over the rapidity interval 2.5–4.0 and $p_{\mathrm{T}} <$ 15 GeV/$c$. Results used for the determination of the integrated nuclear modification factors and of the nuclear modification factors as a function of centrality (Figures 4 and 5 (right)).

Interpolated $p_{\mathrm{T}}$-differential cross section for $\Upsilon(1\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ in proton–proton collisions. Results used for the determination of the nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of its transverse momentum (Figure 6).

Interpolated rapidity-differential cross section for $\Upsilon(1\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ in proton–proton collisions. Results used for the determination of the nuclear modification factor of $\Upsilon(1\mathrm{S})$ as a function of rapidity (Figure 7).

Interpolated rapidity-differential cross section for $\Upsilon(2\mathrm{S}) \rightarrow \mu^{+}\mu^{-}$ in proton–proton collisions. Results used for the determination of the nuclear modification factor of $\Upsilon(2\mathrm{S})$ as a function of rapidity (Figure 7).

$^3\overline{\mathrm{H}} /\,^3\mathrm{H}$ ratio in $\mathrm{INEL}>0$ p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Transverse momentum spectra obtained from the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for $\mathrm{INEL}>0$ events in p$-$Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Transverse momentum spectra obtained from the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for events in 0$-$10% multiplicity class in p$-$Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Transverse momentum spectra obtained from the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for events in 10$-$20% multiplicity class in p$-$Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Transverse momentum spectra obtained from the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for events in 20$-$40% multiplicity class in p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Transverse momentum spectra obtained from the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for events in 40$-$100% multiplicity class in p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Mean transverse momentum obtained from the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ as a function of the mean charged-particle multiplicity density in p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$ TeV

$(^3\mathrm{He} +\,^3\overline{\mathrm{He}}) / (\mathrm{p} + \overline{\mathrm{p}})$ ratio in p$-$Pb as a function of the mean charged-particle multiplicity density

$(^3\mathrm{H} +\,^3\overline{\mathrm{H}}) / (\mathrm{p} + \overline{\mathrm{p}})$ ratio in p$-$Pb as a function of the mean charged-particle multiplicity density

Coalescence parameter $B_3$ calculated using the average of $\mathrm{INEL}>0$ $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ yields

Coalescence parameter $B_3$ calculated using the average of $\mathrm{INEL}>0$ $^3\mathrm{H}$ and $^3\overline{\mathrm{H}}$ yields

Ratio of (anti-)$^3$H and (anti-)$^3$He yields

Coalescence parameter $B_3$ calculated using the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for $\mathrm{INEL}>0$ events in p$-$Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Coalescence parameter $B_3$ calculated using the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for events in 0$-$10% multiplicity class in p$-$Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Coalescence parameter $B_3$ calculated using the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for events in 10$-$20% multiplicity class in p$-$Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Coalescence parameter $B_3$ calculated using the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for events in 20$-$40% multiplicity class in p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Coalescence parameter $B_3$ calculated using the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for events in 40$-$100% multiplicity class in p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Integrated production yield obtained from the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ corrected for the spin degeneracy factor $(2J +1)$ for minimum-bias p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Upper limit on the integrated production yield of $^4\overline{\mathrm{He}}$ at 90$\%$ confidence level corrected for the spin degeneracy factor $(2J +1)$ for minimum-bias p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Integrated production yield obtained from the average of $^3\mathrm{He}$ and $^3\overline{\mathrm{He}}$ for the different multiplicity classes and $\mathrm{INEL} > 0$ p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

Integrated production yield obtained from the average of $^3\mathrm{H}$ and $^3\overline{\mathrm{H}}$ for $\mathrm{INEL} > 0$ p--Pb collisions at $\sqrt{s_\mathrm{NN}} = 5.02$

$\Upsilon(1S)$ production cross-section in $p$Pb and Pb$p$, as a function of $y^*$. The uncertainty is the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ production cross-section in $p$Pb and Pb$p$, as a function of $p_{T}$. The uncertainty is the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ production cross-section in $p$Pb and Pb$p$, as a function of $y^*$. The uncertainty is the sum in quadrature of the statistical and systematic components.

Scaled $pp$ differential cross-section in $p_{T}$ at $\sqrt{s_{NN}}=8.16$ TeV. The first uncertainty is statistical, the second is systematic, which includes the systematic uncertainty from the $pp$ measurement and that estimated by changing the interpolation function.

Scaled $pp$ differential cross-section in $y$ at $\sqrt{s_{NN}}=8.16$ TeV. The first uncertainty is statistical, the second is systematic, which includes the systematic uncertainty from the $pp$ measurement and that estimated by changing the interpolation function.

$\Upsilon(1S)$ nuclear modification factor, $R_{p{\rm Pb}}^{\Upsilon(1S)}$, in $p$Pb and Pb$p$ as a function of $p_{T}$ integrated over $y^*$ in the range $1.5 < y^* < 4.0$ for $p$Pb and $-5.0 < y^* < -2.5$ for Pb$p$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(1S)$ nuclear modification factor, $R_{p{\rm Pb}}^{\Upsilon(1S)}$, in $p$Pb and Pb$p$ as a function of $y^*$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV$/c$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ nuclear modification factor, $R_{p{\rm Pb}}^{\Upsilon(2S)}$, in $p$Pb and Pb$p$ as a function of $p_{T}$ integrated over $y^*$ in the range $1.5 < y^* < 4.0$ for $p$Pb and $-5.0 < y^* < -2.5$ for Pb$p$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ nuclear modification factor, $R_{p{\rm Pb}}^{\Upsilon(2S)}$, in $p$Pb and Pb$p$ as a function of $y^*$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV$/c$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(1S)$ forward-to-backward ratio, $R_{{\rm FB}}^{\Upsilon(1S)}$, as a function of $p_{T}$ integrated over $\vert y^* \vert$ in the range $2.5 < \vert y^* \vert < 4.0$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(1S)$ forward-to-backward ratio, $R_{{\rm FB}}^{\Upsilon(1S)}$, as a function of $\vert y^* \vert$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV/c. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ to $\Upsilon(1S)$ ratio, in $p$Pb and Pb$p$ as a function of $p_{T}$ integrated over $y^*$ in the range $1.5 < y^* < 4.0$ for $p$Pb and $-5.0 < y^* < -2.5$ for Pb$p$. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(2S)$ to $\Upsilon(1S)$ ratio, in $p$Pb and Pb$p$ as a function of $y^*$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV/c. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

$\Upsilon(1S)$ to non-prompt $J/\psi$, in $p$Pb and Pb$p$ as a function of $y^*$ integrated over $p_{T}$ in the range $0 < p_{T} < 25$ GeV/c. The quoted uncertainties are the sum in quadrature of the statistical and systematic components.

The covariance matrix for the kinematic bins matrix of the $D_s^+$ production asymmetry measured using the combined $\sqrt{s} =7$ and 8 TeV data sets (corresponding to Table 1). The statistical and systematic uncertainties are combined. The kinematic bins are indexed as follows: the three rapidity bins, y, in the range from 2.0 to 4.5 are indexed with $i_y$ ranging from 0 to 2 with the following bin edges: 2.0, 3.0, 3.5, 4.5. The four $p_{\rm T}$ bins in the range from 2.5 to 25.0 GeV/c are indexed with $i_{p_T}$ from 0 to 3 with the following bin edges: 2.5, 4.7, 6.5, 8.5, 25.0 GeV/c. The index of a kinematic bins is then given by: $i_{\rm bin} = 4 \times i_y + i_{p_T}$.

The covariance matrix for the kinematic bins matrix of the $D_s^+$ production asymmetry measured using the $\sqrt{s} =7$ TeV data set (corresponding to Table 2). The statistical and systematic uncertainties are combined. The kinematic bins are indexed as follows: the three rapidity bins, y, in the range from 2.0 to 4.5 are indexed with $i_y$ ranging from 0 to 2 with the following bin edges: 2.0, 3.0, 3.5, 4.5. The four $p_{\rm T}$ bins in the range from 2.5 to 25.0 GeV/c are indexed with $i_{p_T}$ from 0 to 3 with the following bin edges: 2.5, 4.7, 6.5, 8.5, 25.0 GeV/c. The index of a kinematic bins is then given by: $i_{\rm bin} = 4 \times i_y + i_{p_T}$.

The covariance matrix for the kinematic bins matrix of the $D_s^+$ production asymmetry measured using the $\sqrt{s} =8$ TeV data set (corresponding to Table 3). The statistical and systematic uncertainties are combined. The kinematic bins are indexed as follows: the three rapidity bins, y, in the range from 2.0 to 4.5 are indexed with $i_y$ ranging from 0 to 2 with the following bin edges: 2.0, 3.0, 3.5, 4.5. The four $p_{\rm T}$ bins in the range from 2.5 to 25.0 GeV/c are indexed with $i_{p_T}$ from 0 to 3 with the following bin edges: 2.5, 4.7, 6.5, 8.5, 25.0 GeV/c. The index of a kinematic bins is then given by: $i_{\rm bin} = 4 \times i_y + i_{p_T}$.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

Differential cross-sections multiplied by dimuon branching fractions for $\Upsilon$(1S), $\Upsilon$(2S) and $\Upsilon$(3S) states versus $p_T$ integrated over $y$ between 2.0 and 4.5.

Differential cross-sections multiplied by dimuon branching fractions for $\Upsilon$(1S), $\Upsilon$(2S) and $\Upsilon$(3S) states versus $y$ integrated over $p_T$ from 0 to 15 GeV/c.

Ratios of double-differential cross-sections times dimuon branching fractions for $\Upsilon$(2S) to $\Upsilon$(1S).

Ratios of double-differential cross-sections times dimuon branching fractions for $\Upsilon$(3S) to $\Upsilon$(1S).

Ratios of differential cross-sections times dimuon branching fractions versus $p_T$ over $y$ for $\Upsilon$(2S) to $\Upsilon$(1S) and $\Upsilon$(3S) to $\Upsilon$(1S).

Ratios of differential cross-sections times dimuon branching fractions versus $y$ over $p_T$ for $\Upsilon$(2S) to $\Upsilon$(1S) and $\Upsilon$(3S) to $\Upsilon$(1S).

The cross-section ratios between 13 TeV and 8 TeV in different bins of $p_T$ and $y$ for $\Upsilon$(1S).

The cross-section ratios between 13 TeV and 8 TeV in different bins of $p_T$ and $y$ for $\Upsilon$(2S).

The cross-section ratios between 13 TeV and 8 TeV in different bins of $p_T$ and $y$ for $\Upsilon$(3S).

Ratios of differential cross-sections between 13 TeV and 8 TeV measurements versus $p_T$ over $y$ for $\Upsilon$(1S), $\Upsilon$(2S) and $\Upsilon$(3S).

Ratios of differential cross-sections between 13 TeV and 8 TeV measurements versus $y$ over $p_T$ for $\Upsilon$(1S), $\Upsilon$(2S) and $\Upsilon$(3S).

Double-differential cross section times the dimuon branching fraction of the Y(2S) meson for different ranges of pT in bins of |y| and for the full |y| < 1.2 range, for the unpolarized decay hypothesis. The global uncertainty in the integrated luminosity of 2.3% is not included in the systematic uncertainties.

Double-differential cross section times the dimuon branching fraction of the Y(3S) meson for different ranges of pT in bins of |y| and for the full |y| < 1.2 range, for the unpolarized decay hypothesis. The global uncertainty in the integrated luminosity of 2.3% is not included in the systematic uncertainties.

Multiplicative scaling factors to obtain the J/psi differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, +0.10, -1) from the unpolarized cross section measurements given in Table 1.

Multiplicative scaling factors to obtain the psi(2S) differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, +0.03, -1) from the unpolarized cross section measurements given in Table 2.

Multiplicative scaling factors to obtain the Y(1S) differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, k, -1) from the unpolarized cross section measurements given in Table 3. The parameter k corresponds to a linear interpolation of the CMS measured value as a function of pT for pT < 50 GeV. For pT > 50 GeV, where no measurements exist, k is taken as the average of all the measured values for pT < 50 GeV.

Multiplicative scaling factors to obtain the Y(2S) differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, k, -1) from the unpolarized cross section measurements given in Table 3. The parameter k corresponds to a linear interpolation of the CMS measured value as a function of pT for pT < 50 GeV. For pT > 50 GeV, where no measurements exist, k is taken as the average of all the measured values for pT < 50 GeV.

Multiplicative scaling factors to obtain the Y(3S) differential cross sections for different polarization scenarios (Lambda_{Theta} = +1, k, -1) from the unpolarized cross section measurements given in Table 3. The parameter k corresponds to a linear interpolation of the CMS measured value as a function of pT for pT < 50 GeV. For pT > 50 GeV, where no measurements exist, k is taken as the average of all the measured values for pT < 50 GeV.

Ratio of differential cross section times branching fractions of the prompt psi(2S) to J/psi mesons assuming unpolarized productions as a function of pT for |y| < 1.2.

Ratios of differential cross section times branching fractions of the prompt Y(2S) to Y(1S) and Y(3S) to Y(1S) mesons assuming unpolarized productions as a function of pT for |y| < 1.2.

The frame-invariant polarization parameter $\tilde{\lambda}$ measured in the helicity frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\theta}$ measured in the helicity frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\theta\phi}$ measured in the helicity frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\phi}$ measured in the helicity frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

The frame-invariant polarization parameter $\tilde{\lambda}$ measured in the helicity frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_\theta$ measured in the Collins-Soper frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\theta\phi}$ measured in the Collins-Soper frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_\phi$ measured in the Collins-Soper frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

The frame-invariant polarization parameter $\tilde{\lambda}$ measured in the Collins-Soper frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_\theta$ measured in the Collins-Soper frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\theta\phi}$ measured in the Collins-Soper frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_\phi$ measured in the Collins-Soper frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

The frame-invariant polarization parameter $\tilde{\lambda}$ measured in the Collins-Soper frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_\theta$ measured in the Gottfried-Jackson frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\theta\phi}$ measured in the Gottfried-Jackson frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_\phi$ measured in the Gottfried-Jackson frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

The frame-invariant polarization parameter $\tilde{\lambda}$ measured in the Gottfried-Jackson frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_\theta$ measured in the Gottfried-Jackson frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\theta\phi}$ measured in the Gottfried-Jackson frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_\phi$ measured in the Gottfried-Jackson frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

The frame-invariant polarization parameter $\tilde{\lambda}$ measured in the Gottfried-Jackson frame for the $\Upsilon(1S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the polarization parameter $\lambda_\theta$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(1S)$ produced at $\sqrt{s}=7$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the polarization parameter $\lambda_{\theta\phi}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(1S)$ produced at $\sqrt{s}=7$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the polarization parameter $\lambda_{\phi}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(1S)$ produced at $\sqrt{s}=7$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the frame-invariant polarization parameter $\tilde{\lambda}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(1S)$ produced at $\sqrt{s}=7$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the polarization parameter $\lambda_{\theta}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(1S)$ produced at $\sqrt{s}=8$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the polarization parameter $\lambda_{\theta\phi}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(1S)$ produced at $\sqrt{s}=8$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the polarization parameter $\lambda_{\phi}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(1S)$ produced at $\sqrt{s}=8$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the frame-invariant polarization parameter $\tilde{\lambda}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(1S)$ produced at $\sqrt{s}=8$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\theta}$ measured in the helicity frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\theta\phi}$ measured in the helicity frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\phi}$ measured in the helicity frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The frame-invariant polarization parameter $\tilde{\lambda}$ measured in the helicity frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\theta}$ measured in the helicity frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\theta\phi}$ measured in the helicity frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\phi}$ measured in the helicity frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The frame-invariant polarization parameter $\tilde{\lambda}$ measured in the helicity frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\theta}$ measured in the Collins-Soper frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\theta\phi}$ measured in the Collins-Soper frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\phi}$ measured in the Collins-Soper frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The frame-invariant polarization parameter $\tilde{\lambda}$ measured in the Collins-Soper frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\theta}$ measured in the Collins-Soper frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\theta\phi}$ measured in the Collins-Soper frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\phi}$ measured in the Collins-Soper frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The frame-invariant polarization parameter $\tilde{\lambda}$ measured in the Collins-Soper frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\theta}$ measured in the Gottfried-Jackson frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\theta\phi}$ measured in the Gottfried-Jackson frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\phi}$ measured in the Gottfried-Jackson frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The frame-invariant polarization parameter $\tilde{\lambda}$ measured in the Gottfried-Jackson frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\theta}$ measured in the Gottfried-Jackson frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\theta\phi}$ measured in the Gottfried-Jackson frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\phi}$ measured in the Gottfried-Jackson frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The frame-invariant polarization parameter $\tilde{\lambda}$ measured in the Gottfried-Jackson frame for the $\Upsilon(2S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

Values of the polarization parameter $\lambda_\theta$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(2S)$ produced at $\sqrt{s}=7$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the polarization parameter $\lambda_{\theta\phi}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(2S)$ produced at $\sqrt{s}=7$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the polarization parameter $\lambda_{\phi}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(2S)$ produced at $\sqrt{s}=7$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the frame-invariant polarization parameter $\tilde{\lambda}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(2S)$ produced at $\sqrt{s}=7$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the polarization parameter $\lambda_{\theta}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(2S)$ produced at $\sqrt{s}=8$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the polarization parameter $\lambda_{\theta\phi}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(2S)$ produced at $\sqrt{s}=8$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the polarization parameter $\lambda_{\phi}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(2S)$ produced at $\sqrt{s}=8$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the frame-invariant polarization parameter $\tilde{\lambda}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(2S)$ produced at $\sqrt{s}=8$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\theta}$ measured in the helicity frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\theta\phi}$ measured in the helicity frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\phi}$ measured in the helicity frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The frame-invariant polarization parameter $\tilde{\lambda}$ measured in the helicity frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\theta}$ measured in the helicity frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\theta\phi}$ measured in the helicity frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\phi}$ measured in the helicity frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The frame-invariant polarization parameter $\tilde{\lambda}$ measured in the helicity frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\theta}$ measured in the Collins-Soper frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\theta\phi}$ measured in the Collins-Soper frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\phi}$ measured in the Collins-Soper frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The frame-invariant polarization parameter $\tilde{\lambda}$ measured in the Collins-Soper frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\theta}$ measured in the Collins-Soper frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\theta\phi}$ measured in the Collins-Soper frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\phi}$ measured in the Collins-Soper frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The frame-invariant polarization parameter $\tilde{\lambda}$ measured in the Collins-Soper frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\theta}$ measured in the Gottfried-Jackson frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\theta\phi}$ measured in the Gottfried-Jackson frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\phi}$ measured in the Gottfried-Jackson frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The frame-invariant polarization parameter $\tilde{\lambda}$ measured in the Gottfried-Jackson frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=7\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\theta}$ measured in the Gottfried-Jackson frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The polarization parameter $\lambda_{\theta\phi}$ measured in the Gottfried-Jackson frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

The polarization parameter $\lambda_{\phi}$ measured in the Gottfried-Jackson frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second is systematic.

The frame-invariant polarization parameter $\tilde{\lambda}$ measured in the Gottfried-Jackson frame for the $\Upsilon(3S)$ state in different bins of $p_{T}^{\Upsilon}$ and three rapidity ranges using data collected at $\sqrt{s}=8\,\mathrm{TeV}$. The first uncertainty is statistical and the second represents the systematic uncertainty.

Values of the polarization parameter $\lambda_\theta$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(3S)$ produced at $\sqrt{s}=7$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the polarization parameter $\lambda_{\theta\phi}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(3S)$ produced at $\sqrt{s}=7$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the polarization parameter $\lambda_{\phi}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(3S)$ produced at $\sqrt{s}=7$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the frame-invariant polarization parameter $\tilde{\lambda}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(3S)$ produced at $\sqrt{s}=7$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the polarization parameter $\lambda_{\theta}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(3S)$ produced at $\sqrt{s}=8$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the polarization parameter $\lambda_{\theta\phi}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(3S)$ produced at $\sqrt{s}=8$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the polarization parameter $\lambda_{\phi}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(3S)$ produced at $\sqrt{s}=8$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

Values of the frame-invariant polarization parameter $\tilde{\lambda}$ measured in the helicity(SH), Collins-Soper(CS) and Gottfried-Jackson(TH) frames for the $\Upsilon(3S)$ produced at $\sqrt{s}=8$ TeV in the rapidity range $2.2 < y^\Upsilon < 4.5$ and different bins of $p_{T}^{\Upsilon}$. The first quoted uncertainty is statistical and the second is systematic.

The cross section measurement as a function of the inclusive jet multiplicity, for jet multiplicities of up to 7.

The differential cross section measurement as a function of the transverse momentum of the first leading jet.

The differential cross section measurement as a function of the transverse momentum of the first leading jet.

The differential cross section measurement as a function of the transverse momentum of the second leading jet.

The differential cross section measurement as a function of the transverse momentum of the second leading jet.

The differential cross section measurement as a function of the transverse momentum of the third leading jet.

The differential cross section measurement as a function of the transverse momentum of the third leading jet.

The differential cross section measurement as a function of the transverse momentum of the fourth leading jet.

The differential cross section measurement as a function of the transverse momentum of the fourth leading jet.

The differential cross section measurement as a function of the scalar sum of the transverse momenta of jets (HT) for events with one or more jets.

The differential cross section measurement as a function of the scalar sum of the transverse momenta of jets (HT) for events with one or more jets.

The differential cross section measurement as a function of the scalar sum of the transverse momenta of jets (HT) for events with two or more jets.

The differential cross section measurement as a function of the scalar sum of the transverse momenta of jets (HT) for events with two or more jets.

The differential cross section measurement as a function of the scalar sum of the transverse momenta of jets (HT) for events with three or more jets.

The differential cross section measurement as a function of the scalar sum of the transverse momenta of jets (HT) for events with three or more jets.

The differential cross section measurement as a function of the scalar sum of the transverse momenta of jets (HT) for events with four or more jets.

The differential cross section measurement as a function of the scalar sum of the transverse momenta of jets (HT) for events with four or more jets.

The differential cross section measurement as a function of the transverse momentum of the dijet system of the two highest pT jets for events with two or more jets.

The differential cross section measurement as a function of the transverse momentum of the dijet system of the two highest pT jets for events with two or more jets.

The differential cross section measurement as a function of the transverse momentum of the dijet system of the two highest pT jets for events with three or more jets.

The differential cross section measurement as a function of the transverse momentum of the dijet system of the two highest pT jets for events with three or more jets.

The differential cross section measurement as a function of the transverse momentum of the dijet system of the two highest pT jets for events with four or more jets.

The differential cross section measurement as a function of the transverse momentum of the dijet system of the two highest pT jets for events with four or more jets.

The differential cross section measurement as a function of the dijet invariant mass of the two highest pT jets for events with two or more jets.

The differential cross section measurement as a function of the dijet invariant mass of the two highest pT jets for events with two or more jets.

The differential cross section measurement as a function of the dijet invariant mass of the two highest pT jets for events with three or more jets.

The differential cross section measurement as a function of the dijet invariant mass of the two highest pT jets for events with three or more jets.

The differential cross section measurement as a function of the dijet invariant mass of the two highest pT jets for events with four or more jets.

The differential cross section measurement as a function of the dijet invariant mass of the two highest pT jets for events with four or more jets.

The differential cross section measurement as a function of the absolute value of the pseudorapidity of the first leading jet.

The differential cross section measurement as a function of the absolute value of the pseudorapidity of the first leading jet.

The differential cross section measurement as a function of the absolute value of the pseudorapidity of the second leading jet.

The differential cross section measurement as a function of the absolute value of the pseudorapidity of the second leading jet.

The differential cross section measurement as a function of the absolute value of the pseudorapidity of the third leading jet.

The differential cross section measurement as a function of the absolute value of the pseudorapidity of the third leading jet.

The differential cross section measurement as a function of the absolute value of the pseudorapidity of the fourth leading jet.

The differential cross section measurement as a function of the absolute value of the pseudorapidity of the fourth leading jet.

The differential cross section measurement as a function of the rapidity difference between the two highest pT jets for events with two or more jets.

The differential cross section measurement as a function of the rapidity difference between the two highest pT jets for events with two or more jets.

The differential cross section measurement as a function of the rapidity difference between the two highest pT jets for events with three or more jets.

The differential cross section measurement as a function of the rapidity difference between the two highest pT jets for events with three or more jets.

The differential cross section measurement as a function of the rapidity difference between the two highest pT jets for events with four or more jets.

The differential cross section measurement as a function of the rapidity difference between the two highest pT jets for events with four or more jets.

The differential cross section measurement as a function of the rapidity difference between the first and the third highest pT jets for events with three or more jets.

The differential cross section measurement as a function of the rapidity difference between the first and the third highest pT jets for events with three or more jets.

The differential cross section measurement as a function of the rapidity difference between the second and the third highest pT jets for events with three or more jets.

The differential cross section measurement as a function of the rapidity difference between the second and the third highest pT jets for events with three or more jets.

The differential cross section measurement as a function of the rapidity difference between the most forward and the most backward jets for events with two or more jets.

The differential cross section measurement as a function of the rapidity difference between the most forward and the most backward jets for events with two or more jets.

The differential cross section measurement as a function of the rapidity difference between the most forward and the most backward jets for events with three or more jets.

The differential cross section measurement as a function of the rapidity difference between the most forward and the most backward jets for events with three or more jets.

The differential cross section measurement as a function of the rapidity difference between the most forward and the most backward jets for events with four or more jets.

The differential cross section measurement as a function of the rapidity difference between the most forward and the most backward jets for events with four or more jets.

The differential cross section measurement as a function of the azimuthal angle separation between the two highest pT jets for events with two or more jets.

The differential cross section measurement as a function of the azimuthal angle separation between the two highest pT jets for events with two or more jets.

The differential cross section measurement as a function of the azimuthal angle separation between the most forward and the most backward jets for events with two or more jets.

The differential cross section measurement as a function of the azimuthal angle separation between the most forward and the most backward jets for events with two or more jets.

The differential cross section measurement as a function of DeltaR separation between the two highest pT jets for events with two or more jets.

The differential cross section measurement as a function of DeltaR separation between the two highest pT jets for events with two or more jets.

The differential cross section measurement as a function of the azimuthal angle separation between the first leading jet and the muon.

The differential cross section measurement as a function of the azimuthal angle separation between the first leading jet and the muon.

The differential cross section measurement as a function of the azimuthal angle separation between the second leading jet and the muon.

The differential cross section measurement as a function of the azimuthal angle separation between the second leading jet and the muon.

The differential cross section measurement as a function of the azimuthal angle separation between the third leading jet and the muon.

The differential cross section measurement as a function of the azimuthal angle separation between the third leading jet and the muon.

The differential cross section measurement as a function of the azimuthal angle separation between the fourth leading jet and the muon.

The differential cross section measurement as a function of the azimuthal angle separation between the fourth leading jet and the muon.

Average number of jets as a function of HT for events with one or more jets.

Average number of jets as a function of HT for events with one or more jets.

Average number of jets as a function of HT for events with two or more jets.

Average number of jets as a function of HT for events with two or more jets.

Average number of jets as a function of the rapidity difference between the two highest pT jets for events with two or more jets.

Average number of jets as a function of the rapidity difference between the two highest pT jets for events with two or more jets.

Average number of jets as a function of the rapidity difference between the most forward and the most backward jets for events with two or more jets.

Average number of jets as a function of the rapidity difference between the most forward and the most backward jets for events with two or more jets.

Differential production cross-sections for prompt $D^{+} + D^{-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

Differential production cross-sections for prompt $D_{s}^{+} + D_{s}^{-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

Differential production cross-sections for prompt $D_{s}^{+} + D_{s}^{-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

Differential production cross-sections for prompt $D^{*+} + D^{*-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

Differential production cross-sections for prompt $D^{*+} + D^{*-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/5}$, for prompt $D^{0} + \bar{D}^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/5}$, for prompt $D^{0} + \bar{D}^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/5}$, for prompt $D^{+} + D^{-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/5}$, for prompt $D^{+} + D^{-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/5}$, for prompt $D_{s}^{+} + D_{s}^{-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/5}$, for prompt $D_{s}^{+} + D_{s}^{-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/5}$, for prompt $D^{*+} + D^{*-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/5}$, for prompt $D^{*+} + D^{*-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D^{+}$ and $D^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D^{+}$ and $D^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D_{s}^{+}$ and $D^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D_{s}^{+}$ and $D^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D^{*+}$ and $D^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D^{*+}$ and $D^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D_{s}^{+}$ and $D^{+}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D_{s}^{+}$ and $D^{+}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D^{*+}$ and $D^{+}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D^{*+}$ and $D^{+}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D_{s}^{+}$ and $D^{*+}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D_{s}^{+}$ and $D^{*+}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

Ratios of B+ differential production cross sections at 13 TeV and at 7 TeV as a function of |yB| for 10 < ptB < 100 GeV or 17 < ptB < 100 GeV. The calculations from FONLL and PYTHIA are provided as well.

B+ differential production cross sections DSIG/DPT for |yB|< 1.45 or |yB|< 2.1, at 13 TeV. The calculations from FONLL and PYTHIA are provided.

B+ differential production cross sections DSIG/DETARAP for 10 < ptB < 100 GeV or 17 < ptB < 100 GeV, at 13 TeV. The calculations from FONLL and PYTHIA are provided.

B+ production cross section for the phase-space region of 10 <= pTB < 17 GeV and |yB| < 1.45, and 17 <= pTB < 100 GeV and |yB| < 2.1. The calculations from FONLL and PYTHIA are provided.

Normalized distribution of the variable $\Delta^{p_{\mathrm{T}}}_{13}$, defined in Eq (16) of the paper, in data after unfolding to particle level.

Normalized distribution of the variable $\Delta^{p_{\mathrm{T}}}_{23}$, defined in Eq (16) of the paper, in data after unfolding to particle level.

Normalized distribution of the variable $\Delta^{p_{\mathrm{T}}}_{14}$, defined in Eq (16) of the paper, in data after unfolding to particle level.

Normalized distribution of the variable $\Delta^{p_{\mathrm{T}}}_{24}$, defined in Eq (16) of the paper, in data after unfolding to particle level.

Normalized distribution of the variable $\Delta\phi_{12}$, defined in Eq (16) of the paper, in data after unfolding to particle level.

Normalized distribution of the variable $\Delta\phi_{13}$, defined in Eq (16) of the paper, in data after unfolding to particle level.

Normalized distribution of the variable $\Delta\phi_{23}$, defined in Eq (16) of the paper, in data after unfolding to particle level.

Normalized distribution of the variable $\Delta\phi_{14}$, defined in Eq (16) of the paper, in data after unfolding to particle level.

Normalized distribution of the variable $\Delta\phi_{24}$, defined in Eq (16) of the paper, in data after unfolding to particle level.

Normalized distribution of the variable $\Delta y_{12}$, defined in Eq (16) of the paper, in data after unfolding to particle level.

Normalized distribution of the variable $\Delta y_{34}$, defined in Eq (16) of the paper, in data after unfolding to particle level.

Normalized distribution of the variable $\Delta y_{13}$, defined in Eq (16) of the paper, in data after unfolding to particle level.

Normalized distribution of the variable $\Delta y_{23}$, defined in Eq (16) of the paper, in data after unfolding to particle level.

Normalized distribution of the variable $\Delta y_{14}$, defined in Eq (16) of the paper, in data after unfolding to particle level.

Normalized distribution of the variable $\Delta y_{24}$, defined in Eq (16) of the paper, in data after unfolding to particle level.

Normalized distribution of the variable $\phi_{1+2} - \phi_{3+4}$, defined in Eq (16) of the paper, in data after unfolding to particle level.

Normalized distribution of the variable $\phi_{1+3} - \phi_{2+4}$, defined in Eq (16) of the paper, in data after unfolding to particle level.

Normalized distribution of the variable $\phi_{1+4} - \phi_{2+3}$, defined in Eq (16) of the paper, in data after unfolding to particle level.

The FSR correction applied as a function of $\phi ^ * _ \eta$ for electrons.

The FSR correction applied as a function of the boson transverse momentum.

The measured differential cross-sections as a function of the boson rapidity for muons. The first uncertainty is due to the size of the dataset, the second is due to experimental systematic uncertainties, and the third is due to the luminosity.

The measured differential cross-sections as a function of the boson rapidity for electrons. The first uncertainty is due to the size of the dataset, the second is due to experimental systematic uncertainties, and the third is due to the luminosity.

The measured differential cross-sections as a function of $\phi ^ * _ \eta$ for muons. The first uncertainty is due to the size of the dataset, the second is due to experimental systematic uncertainties, and the third is due to the luminosity.

The measured differential cross-sections as a function of $\phi ^ * _ \eta$ for electrons. The first uncertainty is due to the size of the dataset, the second is due to experimental systematic uncertainties, and the third is due to the luminosity.

The measured differential cross-sections as a function of pT . The first uncertainty is due to the size of the dataset, the second is due to experimental systematic uncertainties, and the third is due to the luminosity.

The correlation matrix for the differential cross-section measurement as a function of Z boson rapidity, for the dimuon final state, excluding the luminosity uncertainty, which is fully correlated between bins.

The correlation matrix for the differential cross-section measurements as a function of the Z boson rapidity, for the dielectron final state, excluding the luminosity uncertainty, which is fully correlated between bins.

The correlation matrix for the differential cross-section measurement as a function of $\phi ^ * _ \eta$ , for the dimuon final state, excluding the luminosity uncertainty, which is fully correlated between bins.

The correlation matrix for the differential cross-section measurements as a function of $\phi ^ * _ \eta$ , for the dielectron final state, excluding the luminosity uncertainty, which is fully correlated between bins.

The correlation matrix for the differential cross-section measurements as a function of boson pT , for the dimuon final state, excluding the luminosity uncertainty, which is fully correlated between bins.

pT-differential cross section of inclusive Dzero mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV in the rapidity interval -0.96<y<0.04. Branching ratio of D0->Kpi : 0.0388.

pT-differential cross section of prompt Dzero mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV in the rapidity interval -0.96<y<0.04. Branching ratio of D0->Kpi : 0.0388. Data points for pt<2 GeV/c from analysis "without vertexing". Data points for pt>2 GeV/c from the analysis "with vertexing" taken from PRL 113 (2014) 232301 (http://hepdata.cedar.ac.uk/view/ins1296081).

First column: production cross sections per unit of rapidity for prompt D0 mesons, inclusive D0 mesons (no feed-down subtraction) and charm quarks in -0.96 < ycm < 0.04 in p-Pb collisions at 5.02 TeV. For D0 mesons, the second (sys) error is from the luminosity uncertainty, the third (sys) error is from the branching-ratio uncertainties. For charm quarks, the second (sys) error is from the luminosity uncertainty, the third (sys) error is from the Fragmentation Function uncertainties, the fourth (sys) error is from the rapidity shapes of D0 mesons and single charm quarks. Second column: value of <pT> of prompt D0 mesons in p-Pb collisions at 5.02 TeV in the rapidity interval -0.96 < y < 0.04. The first uncertainty is statistical, the second is the systematic uncertainty.

Ratio of D+ to D0 meson yield in p-Pb collisions at sqrt{sNN}=5.02 TeV in -0.96<y< 0.04 as a function of pT. Branching ratio of D+->Kpipi : 0.0913. Branching ratio of D0->Kpi : 0.0388.

Ratio of D*+ to D0 meson yield in p-Pb collisions at sqrt{sNN}=5.02 TeV in -0.96<y< 0.04 as a function of pT. Branching ratio of D*+->D0pi->Kpipi : 0.0388*0.677. Branching ratio of D0->Kpi : 0.0388.

Ratio of Ds to D0 meson yield in p-Pb collisions at sqrt{sNN}=5.02 TeV in -0.96<y< 0.04 as a function of pT. Branching ratio of Ds->phipi->KKpi : 0.0224. Branching ratio of D0->Kpi : 0.0388.

Ratio of Ds to D+ meson yield in p-Pb collisions at sqrt{sNN}=5.02 TeV in -0.96<y< 0.04 as a function of pT. Branching ratio of Ds->phipi->KKpi : 0.0224. Branching ratio of D+->Kpipi : 0.0913.

cross section of prompt Dzero mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 2<pt<5 GeV/c. Branching ratio of D0->Kpi : 0.0388. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dzero mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 5<pt<8 GeV/c. Branching ratio of D0->Kpi : 0.0388. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dzero mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 8<pt<16 GeV/c. Branching ratio of D0->Kpi : 0.0388. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dplus mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 2<pt<5 GeV/c. Branching ratio of D+->Kpipi : 0.0913. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dplus mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 5<pt<8 GeV/c. Branching ratio of D+->Kpipi : 0.0913. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dplus mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 8<pt<16 GeV/c. Branching ratio of D+->Kpipi : 0.0913. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dstar mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 2<pt<5 GeV/c. Branching ratio of D*+->D0pi->Kpipi : 0.0388*0.677. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dstar mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 5<pt<8 GeV/c. Branching ratio of D*+->D0pi->Kpipi : 0.0388*0.677. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

cross section of prompt Dstar mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV vs rapidity in the pt interval 8<pt<16 GeV/c. Branching ratio of D*+->D0pi->Kpipi : 0.0388*0.677. The second (sys) error is the systematic uncertainty from the B feed-down contribution. The first (sys) error is the systematic uncertainty from the other sources.

Nuclear modification factor RpPb of D0 mesons in p-Pb collisions at sqrt{sNN}=5.02 TeV in -0.96<y<0.04 as a function of pT. Data point for 0<pt<1 GeV/c from analysis "without vertexing". Data points for pt>1 GeV/c from the analysis "with vertexing" taken from PRL 113 (2014) 232301 (http://hepdata.cedar.ac.uk/view/ins1296081).

Pt-integrated nuclear modification factor RpA of prompt D0 mesons in -0.96 < ycm < 0.04. The first error is the statistical uncertainty, the second is the systematic uncertainty.

$R_{pPb}$ as a function of $p_{T}$ integrated over rapidity range $-1.8<y*<1.3$ for the 0--90% centrality interval for the three geometric models: (a) Glauber, (b) Glauber-Gribov with $\omega_{\sigma}=0.11$ and (c) Glauber-Gribov with $\omega_{\sigma}=0.2$.

$R_{pPb}$ values as a function of $p_{T}$ in the 0--90% centrality interval averaged over $|y*|<0.5$.The $\langle T_{Pb} \rangle$ value for the ATLAS centrality correction is calculated with the Glauber model. The total systematic uncertainties, which include the uncertainty in $\langle T_{Pb} \rangle$, are indicated by lines of the same colour. Strict quantitative agreement is not expected as each measurement uses different rapidity intervals for the centrality determination and apply different event selection criteria to reject diffractive collisions.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $-1.8<y*<1.3$ for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $-1.8<y*<1.3$ for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $-1.8<y*<1.3$ for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{CP}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $|\eta|<2.3$ for seven centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{CP}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $|\eta|<2.3$ for seven centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{CP}$ as a function of $p_{T}$ extracted from the invariant yields integrated over $|\eta|<2.3$ for seven centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $p_{T}$ extracted from the invariant yields for six rapidity ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

$R_{pPb}$ as a function of $y*$ extracted from the invariant yields for eight $p_{T}$ ranges, for eight centrality intervals, and for different geometrical models used to calculate $\langle T_{Pb} \rangle$.

Inclusive Jet Cross Section for |rapidity| 1.5 TO 2.0 as a function of the jet transverse momentum. Jets are clustered with the anti-kt algorithm ( R = 0.7). The (sys) error is the total systematic error, including the luminosity uncertainty of 2.7%.

Inclusive Jet Cross Section for |rapidity| 2.0 TO 2.5 as a function of the jet transverse momentum. Jets are clustered with the anti-kt algorithm ( R = 0.7). The (sys) error is the total systematic error, including the luminosity uncertainty of 2.7%.

Inclusive Jet Cross Section for |rapidity| 2.5 TO 3.0 as a function of the jet transverse momentum. Jets are clustered with the anti-kt algorithm ( R = 0.7). The (sys) error is the total systematic error, including the luminosity uncertainty of 2.7%.

Inclusive Jet Cross Section for |rapidity| 3.2 TO 4.7 as a function of the jet transverse momentum. Jets are clustered with the anti-kt algorithm ( R = 0.7). The (sys) error is the total systematic error, including the luminosity uncertainty of 2.7%.

Inclusive Jet Cross Section for |rapidity| < 0.5 as a function of the jet transverse momentum. Jets are clustered with the anti-kt algorithm ( R = 0.4). The (sys) error is the total systematic error, including the luminosity uncertainty of 2.7%.

Inclusive Jet Cross Section for |rapidity| 0.5 TO 1.0 as a function of the jet transverse momentum. Jets are clustered with the anti-kt algorithm ( R = 0.4). The (sys) error is the total systematic error, including the luminosity uncertainty of 2.7%.

Inclusive Jet Cross Section for |rapidity| 1.0 TO 1.5 as a function of the jet transverse momentum. Jets are clustered with the anti-kt algorithm ( R = 0.4). The (sys) error is the total systematic error, including the luminosity uncertainty of 2.7%.

Inclusive Jet Cross Section for |rapidity| 1.5 TO 2.0 as a function of the jet transverse momentum. Jets are clustered with the anti-kt algorithm ( R = 0.4). The (sys) error is the total systematic error, including the luminosity uncertainty of 2.7%.

Inclusive Jet Cross Section for |rapidity| 2.0 TO 2.5 as a function of the jet transverse momentum. Jets are clustered with the anti-kt algorithm ( R = 0.4). The (sys) error is the total systematic error, including the luminosity uncertainty of 2.7%.

Inclusive Jet Cross Section for |rapidity| 2.5 TO 3.0 as a function of the jet transverse momentum. Jets are clustered with the anti-kt algorithm ( R = 0.4). The (sys) error is the total systematic error, including the luminosity uncertainty of 2.7%.

Inclusive Jet Cross Section for |rapidity| 3.2 TO 4.7 as a function of the jet transverse momentum. Jets are clustered with the anti-kt algorithm ( R = 0.4). The (sys) error is the total systematic error, including the luminosity uncertainty of 2.7%.

Centrality dependence of the Psi(2S) Q_pPb in the backward rapidity range (-4.46<y_cms<-2.96). The first uncertainty is statistical, the second is a systematic one. The third systematic uncertainty is fully correlated over centrality, while the fourth one is fully correlated over centrality and correlated between the JPsi and the Psi(2S) charmonium states.

Centrality dependence of the Psi(2S) Q_pPb in the forward rapidity range (2.03<y_cms<3.53). The first uncertainty is statistical, the second is the systematic one. The third systematic uncertainty is fully correlated over centrality, while the fourth one is fully correlated over centrality and correlated between the JPsi and the Psi(2S) charmonium states.

Psi(2S)/J/Psi)_pA/(Psi(2S)/J/Psi)_pp double ratio in the backward rapidity range (-4.46<y_cms<-2.96) plotted as a function of the crossing time. The first uncertainty is statistical, the second is the systematic one.

Psi(2S)/J/Psi)_pA/(Psi(2S)/J/Psi)_pp double ratio in the forward rapidity range (2.03<y_cms<3.53) plotted as a function of the crossing time. The first uncertainty is statistical, the second is the systematic one.