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.

Charged-hadron spectrum in the centrality interval 5-10% for p+Pb, divided by 〈TPPB〉. 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 〈TPPB〉. 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 〈TPPB〉. 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 〈TPPB〉. 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 〈TPPB〉. 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 〈TPPB〉. 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 〈TPPB〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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 〈TAA〉. 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.

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.

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

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

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.

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

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.

Double differential cross-section for $J/\psi$-from-$b$ mesons as a function of $p_\perp$ in bins of $y$. The first uncertainties are statistical, the second are the correlated systematic uncertainties shared between bins and the last are the uncorrelated systematic uncertainties.

The fraction of $J/\psi$-from-$b$ mesons (in %) in bins of the $J/\psi$ $p_\perp$ and $y$. The uncertainties are statistical only. The systematic uncertainties are negligible.

The fraction of $J/\psi$-from-$b$ mesons (in %) in bins of the $J/\psi$ $p_\perp$ and $y$. The uncertainties are statistical only. The systematic uncertainties are negligible.

Differential cross-sections as a function of $p_\perp$ integrated over $y$ for prompt $J/\psi$ mesons. The first uncertainties are statistical and the second (third) are uncorrelated (correlated) systematic uncertainties amongst bins.

Differential cross-sections as a function of $p_\perp$ integrated over $y$ for prompt $J/\psi$ mesons. The first uncertainties are statistical and the second (third) are uncorrelated (correlated) systematic uncertainties amongst bins.

Differential cross-sections as a function of $p_\perp$ integrated over $y$ for $J/\psi$-from-$b$ mesons. The first uncertainties are statistical and the second (third) are uncorrelated (correlated) systematic uncertainties amongst bins.

Differential cross-sections as a function of $p_\perp$ integrated over $y$ for $J/\psi$-from-$b$ mesons. The first uncertainties are statistical and the second (third) are uncorrelated (correlated) systematic uncertainties amongst bins.

Differential cross-sections as a function of $y$ integrated over $p_\perp$ for prompt $J/\psi$ mesons. The first uncertainties are statistical and the second (third) are the correlated (uncorrelated) systematic uncertainties.

Differential cross-sections as a function of $y$ integrated over $p_\perp$ for prompt $J/\psi$ mesons. The first uncertainties are statistical and the second (third) are the correlated (uncorrelated) systematic uncertainties.

Differential cross-sections as a function of $y$ integrated over $p_\perp$ for $J/\psi$-from-$b$ mesons. The first uncertainties are statistical and the second (third) are the correlated (uncorrelated) systematic uncertainties.

Differential cross-sections as a function of $y$ integrated over $p_\perp$ for $J/\psi$-from-$b$ mesons. The first uncertainties are statistical and the second (third) are the correlated (uncorrelated) systematic uncertainties.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ in bins of $y$ for prompt $J/\psi$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ in bins of $y$ for prompt $J/\psi$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ in bins of $y$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ in bins of $y$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $y$ integrated over $p_\perp$ for prompt $J/\psi$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $y$ integrated over $p_\perp$ for prompt $J/\psi$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $y$ integrated over $p_\perp$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $y$ integrated over $p_\perp$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ integrated over $y$ for prompt $J/\psi$ mesons.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ integrated over $y$ for prompt $J/\psi$ mesons.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ integrated over $y$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ integrated over $y$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

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.

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 in 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 in 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/7}$, 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/7}$, 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/7}$, 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/7}$, 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/7}$, 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/7}$, 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/7}$, 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/7}$, 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.

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

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.

Differential cross-section $\mathrm{d} \sigma ( pp \to ( \Upsilon \to \mu^+ \mu^- ) X ) / \mathrm{d}p_T$ [pb/(GeV/$c$)] for $2.0 < y < 4.5$.

Differential cross-section $\mathrm{d} \sigma ( pp \to ( \Upsilon \to \mu^+ \mu^- ) X ) / \mathrm{d}y$ [pb] for $p_T < 30$ GeV.

Ratio of double-differential cross-sections $( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(2S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y ) / ( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(1S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y )$.

Ratio of double-differential cross-sections $( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(3S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y ) / ( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(1S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y )$.

Ratio of double-differential cross-sections $( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(3S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y ) / ( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(2S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y )$.

The ratio of differential cross-sections $( \mathrm{d} \sigma ( pp \to ( \Upsilon(iS) \to \mu^+ \mu^- ) X ) / \mathrm{d}p_T ) / ( \mathrm{d} \sigma ( pp \to ( \Upsilon(jS) \to \mu^+ \mu^- ) X ) / \mathrm{d}p_T )$ for $2.0 < y < 4.5$.

The ratio of differential cross-sections $( \mathrm{d} \sigma ( pp \to ( \Upsilon(iS) \to \mu^+ \mu^- ) X ) / \mathrm{d}y ) / ( \mathrm{d} \sigma ( pp \to ( \Upsilon(jS) \to \mu^+ \mu^- ) X ) / \mathrm{d}y )$ for $p_T < 30$ GeV/$c$.

Double-differential cross-section $\mathrm{d}^2 \sigma ( pp \to ( \Upsilon \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d}y$ [pb/(GeV/$c$)].

Double-differential cross-section $\mathrm{d}^2 \sigma ( pp \to ( \Upsilon \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d}y$ [pb/(GeV/$c$)].

Double-differential cross-section $\mathrm{d}^2 \sigma ( pp \to ( \Upsilon \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d}y$ [pb/(GeV/$c$)].

Differential cross-section $\mathrm{d} \sigma ( pp \to ( \Upsilon \to \mu^+ \mu^- ) X ) / \mathrm{d}p_T$ [pb/(GeV/$c$)] for $2.0 < y < 4.5$.

Differential cross-section $\mathrm{d} \sigma ( pp \to ( \Upsilon \to \mu^+ \mu^- ) X ) / \mathrm{d}y$ [pb] for $p_T < 30$ GeV.

Ratio of double-differential cross-sections $( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(2S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y ) / ( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(1S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y )$.

Ratio of double-differential cross-sections $( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(3S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y ) / ( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(1S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y )$.

Ratio of double-differential cross-sections $( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(3S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y ) / ( \mathrm{d}^2 \sigma ( pp \to ( \Upsilon(2S) \to \mu^+ \mu^- ) X ) / \mathrm{d} p_T/\mathrm{d} y )$.

The ratio of differential cross-sections $( \mathrm{d} \sigma ( pp \to ( \Upsilon(iS) \to \mu^+ \mu^- ) X ) / \mathrm{d}p_T ) / ( \mathrm{d} \sigma ( pp \to ( \Upsilon(jS) \to \mu^+ \mu^- ) X ) / \mathrm{d}p_T )$ for $2.0 < y < 4.5$.

The ratio of differential cross-sections $( \mathrm{d} \sigma ( pp \to ( \Upsilon(iS) \to \mu^+ \mu^- ) X ) / \mathrm{d}y ) / ( \mathrm{d} \sigma ( pp \to ( \Upsilon(jS) \to \mu^+ \mu^- ) X ) / \mathrm{d}y )$ for $p_T < 30$ GeV/$c$.

The ratio of differential cross-section $\mathrm{d} \sigma ( pp \to ( \Upsilon \to \mu^+ \mu^- ) X ) / \mathrm{d}p_T$ [pb/(GeV/$c$)] for $\sqrt{s}=8$ and $\sqrt{s}=7$ TeV in $2.0 < y < 4.5$.

The ratio of differential cross-section $\mathrm{d} \sigma ( pp \to ( \Upsilon \to \mu^+ \mu^- ) X ) / \mathrm{d}y$ for $\sqrt{s}=8$ and $\sqrt{s}=7$ TeV and $p_T < 30$ GeV/$c$.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} p_T(D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} y(\Upsilon(1S))}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} y(\Upsilon(1S))}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} y(D^0)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} y(D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{\pi}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} \left| \Delta \phi \right| }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{\pi}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} \left| \Delta \phi \right| }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} \left| \Delta y \right| }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} \left| \Delta y \right| }$ for $2<y(\Upsilon(1S))<4.5$, for $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} p_T(\Upsilon(1S)D^0) }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} p_T(\Upsilon(1S)D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} y(\Upsilon(1S)D^0) }$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} y(\Upsilon(1S)D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} \mathcal{A}_T }$, where $\mathcal{A}_T = \frac{p_T(\Upsilon(1S))-p_T(D^0)}{p_T(\Upsilon(1S))-p_T(D^0)}$, for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} \mathcal{A}_T }$, where $\mathcal{A}_T = \frac{p_T(\Upsilon(1S))-p_T(D^0)}{p_T(\Upsilon(1S))-p_T(D^0)}$, for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^0)}{\mathrm{d} m(\Upsilon(1S)D^0) }$, for $2<y(\Upsilon(1S))<4.5$, $2<y(D^0)<4.5$, $p_T(D^0)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Normalized differential cross-section $\frac{1}{\sigma}\frac{ \mathrm{d}\sigma(\Upsilon(1S)D^+)}{\mathrm{d} m(\Upsilon(1S)D^+)}$ for $2<y(\Upsilon(1S))<4.5$, $2<y(D^+)<4.5$, $p_T(D^+)>1$ GeV/$c$. Only statistical uncertainties are quoted as systematic uncertainties are found to be negligible. The distribution is normalized to unity.

Mean PT and RMS for J/PSI production from b-hadron decay (assuming unpolarised).

DSIG/DY for prompt (unpolarised) J/PSI and from b-hadron decay integrated over PT.

Double differential cross section in PT and YRAP for prompt J/PSI production assuming no polarisation. The first systematic error is the uncorrelated systematic error and the second is the correlated systematic error.

Double differential cross section in PT and YRAP for J/PSI from b-hadron decays assuming no polarisation. The first systematic error is the uncorrelated systematic error and the second is the correlated systematic error.

Double differential cross section in PT and YRAP for prompt J/PSI production assuming fully transversely polarised J/PSIs The first systematic error is the uncorrelated systematic error and the second is the correlated systematic error.

Double differential cross section in PT and YRAP for prompt J/PSI production assuming fully longitudinally polarised J/PSIs The first systematic error is the uncorrelated systematic error and the second is the correlated systematic error.

Fraction of J/PSIs from b-hadron decay in bins of PT and YRAP assuming no polaraisation. The first systematic error is uncorrelated between bins and the second is the uncertainty due to the unknown polarisation of the prompt J/PSI.

Covariance matrix of total experimental uncertainties (statistical and systematic uncertainties added in quadrature) of absolute double differential fiducial cross section. The bin index is PT_i + 10*y_j.

Inclusive cross-section for $Z$ boson production in bins of PHI*. The uncertainties are statistical, systematic, beam and luminosity.

Inclusive cross-section for $Z$ boson production in bins of muon pseudorapidity. The uncertainties are statistical, systematic, beam and luminosity.

$W^{\pm}$ to $Z$ boson cross-section ratio in bins of muon pseudorapidity. The uncertainties are statistical, systematic, beam and luminosity.

Ratio of $W^+$ to $W^-$ cross-section in bins of muon pseudorapidity. The uncertainties are statistical, systematic and beam.

Lepton charge asymmetry in bins of muon pseudorapidity. The uncertainties are statistical, systematic and beam.

Ratios of $W^{\pm}$ to $Z$ boson cross-section ratios at differnent centre-of-mass energy (8 TeV to 7 TeV) in bins of muon pseudorapidity. The uncertainties are statistical and systematic.

Inclusive cross-section for $Z$ boson production in bins of muon pseudorapidity at centre-of-mass energy of 7 TeV. The uncertainties are statistical, systematic, beam and luminosity.

$W^{\pm}$ to $Z$ boson cross-section ratio in bins of muon pseudorapidity at centre-of-mass energy of 7 TeV. The uncertainties are statistical, systematic, beam and luminosity.

Correlation coefficients of differential cross-section measurements as a function of lepton $\eta$. The beam energy and luminosity uncertainties, which are fully correlated between cross-section measurements, are excluded. This table lists bin-to-bin correlations for $W^{+}$.

Correlation coefficients of differential cross-section measurements as a function of lepton $\eta$. The beam energy and luminosity uncertainties, which are fully correlated between cross-section measurements, are excluded. This table lists bin-to-bin correlations for $W^{-}$.

Correlation coefficients of differential cross-section measurements as a function of lepton $\eta$. The beam energy and luminosity uncertainties, which are fully correlated between cross-section measurements, are excluded. This table lists differential cross-section correlations between $W^{-}$ and $W^{+}$ bins.

Correlation coefficients of differential cross-section measurements as a function of $y_{Z}$. The beam energy and luminosity uncertainties, which are fully correlated between cross-section measurements, are excluded.

Correlation coefficients of differential cross-section measurements as a function of transverse momentum. The beam energy and luminosity uncertainties, which are fully correlated between cross-section measurements, are excluded.

Correlation coefficients of differential cross-section measurements as a function of $\phi^{*}$. The beam energy and luminosity uncertainties, which are fully correlated between cross-section measurements, are excluded.

Correlation coefficients between differential cross-section measurements as a function $y_{Z}$ and $W^{+}$ boson muon $\eta$. The beam energy and luminosity uncertainties, which are fully correlated between cross-section measurements, are excluded.

Correlation coefficients between differential cross-section measurements as a function $y_{Z}$ and $W^{-}$ boson muon $\eta$. The beam energy and luminosity uncertainties, which are fully correlated between cross-section measurements, are excluded.

Correlation coefficients of differential cross-section measurements as a function of $\eta^{\mu^{+}}$ of the Z boson. The beam energy and luminosity uncertainties, which are fully correlated between cross-section measurements, are excluded.

Correlation coefficients of differential cross-section measurements as a function of $\eta^{\mu^{-}}$ of the Z boson. The beam energy and luminosity uncertainties, which are fully correlated between cross-section measurements, are excluded.

Correlation coefficients of differential cross-section measurements as a function of $\eta^{\mu^{+}}$. The beam energy and luminosity uncertainties, which are fully correlated between cross-section measurements, are excluded.

Correlation coefficients of differential cross-section measurements as a function of $\eta^{\mu^{-}}$. The beam energy and luminosity uncertainties, which are fully correlated between cross-section measurements, are excluded.

No description provided.

No description provided.

No description provided.

The nuclear modification factors $R_{p+A}^{\rm{FONLL}}(p_{\rm{T}})$ of $B^{+}$ measured in $p$ + Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. The statistical and systematic uncertainties on the $p$ + Pb data are shown, followed by upper and lower systematic uncertainties from the FONLL predictions.

The nuclear modification factors $R_{p+A}^{\rm{FONLL}}(p_{\rm{T}})$ of $B^{0}$ measured in $p$ + Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. The statistical and systematic uncertainties on the $p$ + Pb data are shown, followed by upper and lower systematic uncertainties from the FONLL predictions.

The nuclear modification factors $R_{p+A}^{\rm{FONLL}}(p_{\rm{T}})$ of $B_{s}^{0}$ measured in $p$ + Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. The statistical and systematic uncertainties on the $p$ + Pb data are shown, followed by upper and lower systematic uncertainties from the FONLL predictions.

The measured $y_{\rm{CM}}$ (rapidity in the center-of-mass frame) differential production cross section of $B^{+}$ in $p$ + Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV, together with the cross section calculated by the FONLL model.

The nuclear modification factors $R_{p+A}^{\rm{FONLL}}(y_{\rm{CM}})$ of $B^{+}$ measured in $p$ + Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV. The statistical and systematic uncertainties on the $p$ + Pb data are shown, followed by upper and lower systematic uncertainties from the FONLL predictions.