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.

Detailed list of the contribution of each source of uncertainty to the total uncertainty on the measured values of $\sigma(tq)$, $\sigma(\bar{t}q)$, $R_t$, and $\sigma(tq+\bar{t}q)$. The evaluation of the systematic uncertainties has a statistical uncertainty of $0.3\,\%$. Uncertainties contributing less than $1.0\,\%$ are marked with "$<1$" in the paper. To provide numerical values for this table in HEPdata, these uncertainties are approximated with $\pm 0.5\,\%$. This approximation is applied to all measurements for the following uncertainties$:$ JES statistical, JES physics modeling, JES mixed detector and modeling, JES close-by-jets, JES pileup, $b$-JES, jet vertex fraction, mistag efficiency and $W+\;$jets shape variation. For the measurement of $\sigma(tq)$ the approximation is applied in addition to the following uncertainties$:$ JES flavor response, $c$-tagging efficiency, $t\bar{t}$ generator + parton shower and $t\bar{t}$ ISR/FSR. For the measurement of $\sigma(\bar{t}q)$ the approximation is applied in addition to these uncertainties$:$ JES flavor response, $b/\bar{b}$ acceptance, and $t\bar{t}$ ISR/FSR. For the measurement of $R_t$ the approximation is applied in addition to these uncertainties$:$ JES detector, $b$-tagging efficiency, $c$-tagging efficiency, $b/\bar{b}$ acceptance and $tq$ scale variations. For the measurement of $\sigma(tq+\bar{t}q)$ the approximation is applied in addition to these uncertainties$:$ JES flavor response, $c$-tagging efficiency, $b/\bar{b}$ acceptance, $t\bar{t}$ generator + parton shower and $t\bar{t}$ ISR/FSR.

Differential t-channel top-quark production cross sections and normalized differential t-channel top-quark production cross sections as functions of ABS(YRAP(T)).

The cross sections for top-quark and top-antiquark production in the t-channel, together with the cross-section ratio.

Differential t-channel top-quark production cross sections and normalized differential t-channel top-quark production cross sections as functions of ABS(YRAP(TBAR)).

Parametrization factors for the $m_{t}$ dependence [see Eq. (4) in the paper] of $\sigma(tq)$, $\sigma(\bar tq)$, and $\sigma(t q+\bar t q)$.

The cross sections for top-quark and top-antiquark production in the t-channel, together with the cross-section ratio.

The value of $|V_{tb}|^{2}$ is extracted by dividing the measured value of $\sigma(t q+\bar t q)$ by the prediction of the approximate NNLO calculation. The experimental and theoretical uncertainties are added in quadrature.

Event yields for the 2-jet-$\ell^{+}$ and 2-jet-$\ell^{-}$ HPR channels. The expectation for the signal and background yields correspond to the result of the binned maximum-likelihood fit described in Sec. $\mathrm{VI~D}$ in the paper. The uncertainty of the expectations is the normalization uncertainty of each process after the fit, as described in Sec. $\mathrm{VII~F}$ in the paper.

Migration matrix relating the parton level $p_{\mathrm{T}}(t)$ shown on the $y$ axis and the reconstruction level $p_{\mathrm{T}}(\ell \nu b)$ shown on the $x$ axis for the top-quark $p_{\mathrm{T}}$ distribution.

Migration matrix relating the parton level $p_{\mathrm{T}}(\bar t)$ shown on the $y$ axis and the reconstruction level $p_{\mathrm{T}}(\ell \nu b)$ shown on the $x$ axis for the top-antiquark $p_{\mathrm{T}}$ distribution.

Migration matrix relating the parton level $|y(t)|$ shown on the $y$ axis and the reconstruction level $|y(\ell \nu b)|$ shown on the $x$ axis for the top-quark $|y|$ distribution.

Migration matrix relating the parton level $|y(\bar t)|$ shown on the $y$ axis and the reconstruction level $|y(\ell \nu b)|$ shown on the $x$ axis for the top-antiquark $|y|$ distribution.

Differential t-channel top-quark production cross section as a function of $p_{\mathrm{T}}(t)$ with the uncertainties for each bin given in percent.

Differential t-channel top-quark production cross section as a function of $p_{\mathrm{T}}(\bar t)$ with the uncertainties for each bin given in percent.

Differential t-channel top-quark production cross section as a function of $|y(t)|$ with the uncertainties for each bin given in percent.

Differential t-channel top-quark production cross section as a function of $|y(\bar t)|$ with the uncertainties for each bin given in percent.

Normalized differential t-channel top-quark production cross section as a function of $p_{\mathrm{T}}(t)$ with the uncertainties for each bin given in percent.

Normalized differential t-channel top-quark production cross section as a function of $p_{\mathrm{T}}(\bar t)$ with the uncertainties for each bin given in percent.

Normalized differential t-channel top-quark production cross section as a function of $|y(t)|$ with the uncertainties for each bin given in percent.

Normalized differential t-channel top-quark production cross section as a function of $|y(\bar t)|$ with the uncertainties for each bin given in percent.

Comparison between the measured differential cross sections and the predictions from the NLO calculation using the MSTW2008 PDF set. For each variable and prediction a $\chi^{2}$ value is calculated with HERAfitter using the covariance matrix of each measured spectrum. The theory uncertainties of the predictions are treated as uncorrelated. The number of degrees of freedom (NDF) is equal to the number of bins in the measured spectrum.

Detailed list of the contribution of each source of uncertainty to the total relative uncertainty on the measured $\dfrac{\mathrm{d}\sigma}{\mathrm{d}p_{\mathrm{T}}(t)}$ distribution given in percent for each bin. The list includes only those uncertainties that contribute with more than $1\%$. The following uncertainties contribute to the total uncertainty with less than $1\%$ to each bin content$:$ JES detector, JES statistical, JES physics modeling, JES mixed detector and modeling, JES close-by jets, JES pileup, JES flavor composition, JES flavor response, jet-vertex fraction, $b/\bar{b}$ acceptance, $E_{\mathrm{T}}^{\mathrm{miss}}$ modeling, $W+$ jets shape variation, and $t \bar{t}$ generator. In cases when the uncertainty is report to be "$<1\%$" in the table of the paper the uncertainty is approximated by a value of $0.5\%$.

Detailed list of the contribution of each source of uncertainty to the total relative uncertainty on the measured $\dfrac{\mathrm{d}\sigma}{\mathrm{d}p_{\mathrm{T}}(\bar t)}$ distribution given in percent for each bin. The list includes only those uncertainties that contribute with more than $1\%$. The following uncertainties contribute to the total uncertainty with less than $1\%$ to each bin content$:$ JES detector, JES statistical, JES physics modeling, JES mixed detector and modeling, JES close-by jets, JES pileup, JES flavor composition, JES flavor response, $b-$JES, jet-vertex fraction, mistag efficiency, $b/\bar{b}$ acceptance, $E_{\mathrm{T}}^{\mathrm{miss}}$ modeling, $W+$ jets shape variation, and $t \bar{t}$ generator. In cases when the uncertainty is report to be "$<1\%$" in the table of the paper the uncertainty is approximated by a value of $0.5\%$.

Detailed list of the contribution of each source of uncertainty to the total relative uncertainty on the measured $\dfrac{\mathrm{d}\sigma}{\mathrm{d}|y(t)|}$ distribution given in percent for each bin. The list includes only those uncertainties that contribute with more than $1\%$. The following uncertainties contribute to the total uncertainty with less than $1\%$ to each bin content$:$ JES detector, JES statistical, JES physics modeling, JES mixed detector and modeling, JES close-by jets, JES pileup, JES flavor composition, JES flavor response, jet-vertex fraction, $b/\bar{b}$ acceptance, $E_{\mathrm{T}}^{\mathrm{miss}}$ modeling, $W+$ jets shape variation, $t \bar{t}$ generator, $t \bar{t}$ ISR/FSR, and unfolding. In cases when the uncertainty is report to be "$<1\%$" in the table of the paper the uncertainty is approximated by a value of $0.5\%$.

Detailed list of the contribution of each source of uncertainty to the total relative uncertainty on the measured $\dfrac{\mathrm{d}\sigma}{\mathrm{d}|y(\bar t)|}$ distribution given in percent for each bin. The list includes only those uncertainties that contribute with more than $1\%$. The following uncertainties contribute to the total uncertainty with less than $1\%$ to each bin content$:$ JES detector, JES statistical, JES physics modeling, JES mixed detector and modeling, JES close-by jets, JES pileup, JES flavor composition, JES flavor response, $b-$JES, jet-vertex fraction, $b/\bar{b}$ acceptance, mistag efficiency, $E_{\mathrm{T}}^{\mathrm{miss}}$ modeling, $W+$ jets shape variation, $t \bar{t}$ generator, $t \bar{t}$ ISR/FSR, and unfolding. In cases when the uncertainty is report to be "$<1\%$" in the table of the paper the uncertainty is approximated by a value of $0.5\%$.

Detailed list of the contribution of each source of uncertainty to the total relative uncertainty on the measured $\dfrac{1}{\sigma}\dfrac{\mathrm{d}\sigma}{\mathrm{d}p_{\mathrm{T}}(t)}$ distribution given in percent for each bin. The list includes only those uncertainties that contribute with more than $1\%$. The JES $\eta$ intercalibration uncertainty has a sign switch from the first to the second bin. For the $tq$ generator $+$ parton shower uncertainty a sign switch is denoted with $\mp$. The following uncertainties contribute to the total uncertainty with less than $1\%$ to each bin content$:$ JES detector, JES statistical, JES physics modeling, JES mixed detector and modeling, JES close-by jets, JES pileup, JES flavor composition, JES flavor response, $b-$JES, jet-vertex fraction, $b/\bar{b}$ acceptance, $c-$tagging efficiency, $E_{\mathrm{T}}^{\mathrm{miss}}$ modeling, lepton uncertainties, $W+$ jets shape variation, and $t \bar{t}$ generator. In cases when the uncertainty is report to be "$<1\%$" in the table of the paper the uncertainty is approximated by a value of $0.5\%$.

Detailed list of the contribution of each source of uncertainty to the total relative uncertainty on the measured $\dfrac{1}{\sigma}\dfrac{\mathrm{d}\sigma}{\mathrm{d}p_{\mathrm{T}}(\bar t)}$ distribution given in percent for each bin. The list includes only those uncertainties that contribute with more than $1\%$. Sign switches within one uncertainty are denoted with $\mp$ and $\pm$. The following uncertainties contribute to the total uncertainty with less than $1\%$ to each bin content$:$ JES detector, JES statistical, JES physics modeling, JES mixed detector and modeling, JES close-by jets, JES pileup, JES flavor composition, JES flavor response, $b-$JES, jet-vertex fraction, $b/\bar{b}$ acceptance, mistag efficiency, $E_{\mathrm{T}}^{\mathrm{miss}}$ modeling, lepton uncertainties, $W+$ jets shape variation, and $t \bar{t}$ generator. In cases when the uncertainty is report to be "$<1\%$" in the table of the paper the uncertainty is approximated by a value of $0.5\%$.

Detailed list of the contribution of each source of uncertainty to the total relative uncertainty on the measured $\dfrac{1}{\sigma}\dfrac{\mathrm{d}\sigma}{\mathrm{d}|y(t)|}$ distribution given in percent for each bin. The list includes only those uncertainties that contribute with more than $1\%$. Sign switches within one uncertainty are denoted with $\mp$ and $\pm$. The following uncertainties contribute to the total uncertainty with less than $1\%$ to each bin content$:$ JES detector, JES statistical, JES physics modeling, JES mixed detector and modeling, JES close-by jets, JES pileup, JES flavor composition, JES flavor response, b-JES, jet-vertex fraction, $b/\bar{b}$ acceptance, $b-$tagging efficiency, $c-$tagging efficiency, mistag efficiency, $E_{\mathrm{T}}^{\mathrm{miss}}$ modeling, lepton uncertainties, $W+$jets shape variation, $t \bar{t}$ generator, $t \bar{t}$ISR/FSR, and unfolding. In cases when the uncertainty is report to be "$<1\%$" in the table of the paper the uncertainty is approximated by a value of $0.5\%$.

Detailed list of the contribution of each source of uncertainty to the total relative uncertainty on the measured $\dfrac{1}{\sigma}\dfrac{\mathrm{d}\sigma}{\mathrm{d}|y(\bar t)|}$ distribution given in percent for each bin. The list includes only those uncertainties that contribute with more than $1\%$. Sign switches within one uncertainty are denoted with $\mp$ and $\pm$. The following uncertainties contribute to the total uncertainty with less than $1\%$ to each bin content $:$ JES detector, JES statistical, JES physics modeling, JES mixed detector and modeling, JES close-by jets, JES pileup, JES flavor composition, JES flavor response, b-JES, jet energy resolution, jet-vertex fraction, $b/\bar{b}$ acceptance, $b-$tagging efficiency, $c-$ tagging efficiency, mistag efficiency, $E_{\mathrm{T}}^{\mathrm{miss}}$ modeling, lepton uncertainties, $W+$ jets shape variation, $t \bar{t}$ generator, $t \bar{t}$ ISR/FSR, and unfolding. In cases when the uncertainty is report to be "$<1\%$" in the table of the paper the uncertainty is approximated by a value of $0.5\%$.

Statistical correlation matrices for the differential cross section as a function of $p_{\mathrm{T}}(t)$.

Statistical correlation matrices for the differential cross section as a function of $p_{\mathrm{T}}(\bar t)$.

Statistical correlation matrices for the differential cross section as a function of $|y(t)|$.

Statistical correlation matrices for the differential cross section as a function of $|y(\bar t)|$.

Statistical correlation matrices for the normalized differential cross section as a function of $p_{\mathrm{T}}(t)$.

Statistical correlation matrices for the normalized differential cross section as a function of $p_{\mathrm{T}}(\bar t)$.

Statistical correlation matrices for the normalized differential cross section as a function of $|y(t)|$.

Statistical correlation matrices for the normalized differential cross section as a function of $|y(\bar t)|$.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Differential cross section for prompt D+ production.

Differential cross section for prompt D*+ production.

Differential cross section for prompt D*+ production.

Integrated cross sections for prompt charm particle production extraplated to PT=0. The second (sys) error is the uncertainty in the branching ratio The third (sys) error is the uncertainty in the extrapolation procedure.

Integrated cross sections for prompt charm particle production extraplated to PT=0. The second (sys) error is the uncertainty in the branching ratio The third (sys) error is the uncertainty in the extrapolation procedure.

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

Measured differential cross section with associated uncertainties as a function of absolute rapidity of diphoton system. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of multiplicity of jets with transverse momentum pT(jet) > 30 GeV. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of multiplicity of jets with transverse momentum pT(jet) > 30 GeV. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of multiplicity of jets with transverse momentum pT(jet) > 30 GeV. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of multiplicity of jets with transverse momentum pT(jet) > 30 GeV. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of the leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of the leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of absolute rapidity of the leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of absolute rapidity of the leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of scalar transverse momentum sum of all jets. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of scalar transverse momentum sum of all jets. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of second leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of second leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of rapidity separation of leading two jets. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of rapidity separation of leading two jets. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of azimuthal angle between diphoton and dijet systems. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of azimuthal angle between diphoton and dijet systems. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of cosine of the decay angle in the Collins-Soper frame. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of cosine of the decay angle in the Collins-Soper frame. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of cosine of the decay angle in the Collins-Soper frame. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of cosine of the decay angle in the Collins-Soper frame. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of azimuthal angle between the two leading jets. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of azimuthal angle between the two leading jets. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of diphoton thrust pT. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of diphoton thrust pT. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of rapidity separation of the two photons. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of rapidity separation of the two photons. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of tau of highest-tau jet (see paper for description). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of tau of highest-tau jet (see paper for description). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of scalar sum of tau for all jets (see paper for description). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of scalar sum of tau for all jets (see paper for description). Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of absolute rapidity of second leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of absolute rapidity of second leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of third leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of third leading jet. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of dijet invariant mass. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of dijet invariant mass. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of the combined diphoton and dijet system. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of transverse momentum of the combined diphoton and dijet system. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of cosine of the decay angle in the Collins-Soper frame in bins of diphoton transverse momentum. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of cosine of the decay angle in the Collins-Soper frame in bins of diphoton transverse momentum. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of diphoton transverse momentum in jet multiplicity bins. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of diphoton transverse momentum in jet multiplicity bins. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of leading jet transverse momentum for exclusive one-jet events. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Measured differential cross section with associated uncertainties as a function of leading jet transverse momentum for exclusive one-jet events. Each systematic uncertainty sources is fully uncorrelated with the other sources and fully correlated across bins, except for the background modelling systematics for which an uncorrelated treatment across bins is more appropriate.

Diphoton kinematic acceptances in percent for gluon fusion for each fiducial region/variable bin studied in this paper, defined as the probability to fulfil the diphoton kinematic criteria: pT(2gam) < 0.35 (0.25) for the leading (subleading) photon and |eta(gam)|<2.37. The factors are evaluated using the POWHEG event generator with MPI modelling and hadronisation turned off. Consistent results for the diphoton variables are obtained by HRes 2.2. Uncertainties are taken from PDF variations. QCD scale varaitions have a negligible impact on these factors.

Diphoton kinematic acceptances in percent for gluon fusion for each fiducial region/variable bin studied in this paper, defined as the probability to fulfil the diphoton kinematic criteria: pT(2gam) < 0.35 (0.25) for the leading (subleading) photon and |eta(gam)|<2.37. The factors are evaluated using the POWHEG event generator with MPI modelling and hadronisation turned off. Consistent results for the diphoton variables are obtained by HRes 2.2. Uncertainties are taken from PDF variations. QCD scale varaitions have a negligible impact on these factors.

Isolation efficiencies in percent for gluon fusion H -> 2gam for each fiducial region/variable bin measured in this analysis. The isolation efficiency is defined as the probability for both photons to fulfil the isolation criteria (ETiso < 14 GeV as described in the text) for events that pass the diphoton kinematic criteria. Uncertainties are assigned in the same way as for the non-perturbative correction factors: by varying the fragmentation and underlying event modelling. These factors can be multiplied by the kinematic acceptance factors (see table~ ef{tab:fid_acceptance}) to extrapolate an inclusive gluon fusion Higgs prediction to the fiducial volume used in this analysis.

Isolation efficiencies in percent for gluon fusion H -> 2gam for each fiducial region/variable bin measured in this analysis. The isolation efficiency is defined as the probability for both photons to fulfil the isolation criteria (ETiso < 14 GeV as described in the text) for events that pass the diphoton kinematic criteria. Uncertainties are assigned in the same way as for the non-perturbative correction factors: by varying the fragmentation and underlying event modelling. These factors can be multiplied by the kinematic acceptance factors (see table~ ef{tab:fid_acceptance}) to extrapolate an inclusive gluon fusion Higgs prediction to the fiducial volume used in this analysis.

Non-perturbative correction factors in percent accounting for the impact of hadronisation and the underlying event activity for all measured variables and fiducial regions. Uncertainties are evaluated by deriving these factors using different generators and tunes as described in the text.

Non-perturbative correction factors in percent accounting for the impact of hadronisation and the underlying event activity for all measured variables and fiducial regions. Uncertainties are evaluated by deriving these factors using different generators and tunes as described in the text.

Fiducial cross sections in fb from WH, ZH and ttH combined, in each variable bin and fiducial region.

Fiducial cross sections in fb from WH, ZH and ttH combined, in each variable bin and fiducial region.

Statistical bin-to-bin correlations between the 24 bins used in the analysis described in Phys.Lett. B753 (2016) 69-85 (arXiv:1508.02507). See Figure 3 of this paper for the bins used.

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.

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

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

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

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

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

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

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

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

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

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

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

Measured cross section in bins of $p_{\rm{T}}^{\rm{j}1}$.

Measured fractions in bins of $p_{\rm{T}}^{\rm{H}}$.

Measured fractions in bins of $|y^{\rm{H}}|$.

Measured fractions in bins of $N_{\rm{jets}}$.

Measured fractions in bins of $p_{\rm{T}}^{\rm{j1}}$.

Correlation matrix for the total uncertainty of the differential cross-section measurement in bins of $p_{\rm{T}}^{\rm{H}}$.

Correlation matrix for the total uncertainty of the differential cross-section measurement in bins of $|y^{\rm{H}}|$.

Correlation matrix for the total uncertainty of the differential cross-section measurement in bins of $N_{\rm{jets}}$.

Correlation matrix for the total uncertainty of the differential cross-section measurement in bins of $p_{\rm{T}}^{\rm{j1}}$.

Correlation matrix for the total uncertainty of the differential shape measurement in bins of $p_{\rm{T}}^{\rm{H}}$.

Correlation matrix for the total uncertainty of the differential shape measurement in bins of $|y^{\rm{H}}|$.

Correlation matrix for the total uncertainty of the differential shape measurement in bins of $N_{\rm{jets}}$.

Correlation matrix for the total uncertainty of the differential shape measurement in bins of $p_{\rm{T}}^{\rm{j1}}$.

The total $pp\to H$ cross section in the $H\to\gamma\gamma$ and $H\to ZZ^*\to 4\ell$ channels and their combination.

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.

No description provided.

No description provided.