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.

Differential cross sections of the Y(2S) meson as a function of pT for pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross sections of the Y(3S) meson as a function of pT for pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross sections of the Y(3S) meson as a function of pT for pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross sections of the Y(1S) meson as a function of pT for PbPb collisions. Global uncertainties arise from the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross sections of the Y(1S) meson as a function of pT for PbPb collisions. Global uncertainties arise from the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross sections of the Y(2S) meson as a function of pT for PbPb collisions. Global uncertainties arise from the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross sections of the Y(2S) meson as a function of pT for PbPb collisions. Global uncertainties arise from the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross sections of the Y(3S) meson as a function of pT for PbPb collisions.

Differential cross sections of the Y(3S) meson as a function of pT for PbPb collisions.

Differential cross sections of the Y(1S) meson as a function of rapidity for pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross sections of the Y(1S) meson as a function of rapidity for pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross sections of the Y(2S) meson as a function of rapidity for pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross sections of the Y(2S) meson as a function of rapidity for pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross sections of the Y(3S) meson as a function of rapidity for pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross sections of the Y(3S) meson as a function of rapidity for pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross sections of the Y(1S) meson as a function of pT for PbPb collisions. Global uncertainties arise from the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross sections of the Y(1S) meson as a function of pT for PbPb collisions. Global uncertainties arise from the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross sections of the Y(2S) meson as a function of pT for PbPb collisions. Global uncertainties arise from the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross sections of the Y(2S) meson as a function of pT for PbPb collisions. Global uncertainties arise from the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross sections of the Y(3S) meson as a function of pT for PbPb collisions.

Differential cross sections of the Y(3S) meson as a function of pT for PbPb collisions.

Nuclear modification factor of the Y(1S) meson as a function of pT. Global uncertainty combines the uncertainties from TAA, pp luminosity, and number of minimum bias events in PbPb collisions.

Nuclear modification factor of the Y(1S) meson as a function of pT. Global uncertainty combines the uncertainties from TAA, pp luminosity, and number of minimum bias events in PbPb collisions.

Nuclear modification factor of the Y(2S) meson as a function of pT. Global uncertainty combines the uncertainties from TAA, pp luminosity, and number of minimum bias events in PbPb collisions.

Nuclear modification factor of the Y(2S) meson as a function of pT. Global uncertainty combines the uncertainties from TAA, pp luminosity, and number of minimum bias events in PbPb collisions.

Nuclear modification factor of the Y(3S) meson as a function of pT.

Nuclear modification factor of the Y(3S) meson as a function of pT.

Nuclear modification factor of the Y(1S) meson as a function of rapidity. Global uncertainty combines the uncertainties from TAA, pp luminosity, and number of minimum bias events in PbPb collisions.

Nuclear modification factor of the Y(1S) meson as a function of rapidity. Global uncertainty combines the uncertainties from TAA, pp luminosity, and number of minimum bias events in PbPb collisions.

Nuclear modification factor of the Y(2S) meson as a function of rapidity. Global uncertainty combines the uncertainties from TAA, pp luminosity, and number of minimum bias events in PbPb collisions.

Nuclear modification factor of the Y(2S) meson as a function of rapidity. Global uncertainty combines the uncertainties from TAA, pp luminosity, and number of minimum bias events in PbPb collisions.

Nuclear modification factor of the Y(3S) meson as a function of rapidity.

Nuclear modification factor of the Y(3S) meson as a function of rapidity.

Nuclear modification factor of the Y(1S) meson as a function of collision centrality. Global uncertainty combines the pp luminosity and number of minimum bias events in PbPb collisions. The PP Y(1S) global uncertainty consists of the statistical and systematic uncertainties of pp yield for Y(1S) state.

Nuclear modification factor of the Y(1S) meson as a function of collision centrality. Global uncertainty combines the pp luminosity and number of minimum bias events in PbPb collisions. The PP Y(1S) global uncertainty consists of the statistical and systematic uncertainties of pp yield for Y(1S) state.

Nuclear modification factor of the Y(2S) meson as a function of collision centrality. Global uncertainty combines the pp luminosity and number of minimum bias events in PbPb collisions. The PP Y(2S) global uncertainty consists of the statistical and systematic uncertainties of pp yield for Y(2S) state.

Nuclear modification factor of the Y(2S) meson as a function of collision centrality. Global uncertainty combines the pp luminosity and number of minimum bias events in PbPb collisions. The PP Y(2S) global uncertainty consists of the statistical and systematic uncertainties of pp yield for Y(2S) state.

Confidence intervals on the nuclear modification factors of the Y(3S) mesons as a function of collision centrality.

Confidence intervals on the nuclear modification factors of the Y(3S) mesons as a function of collision centrality.

Confidence intervals on the nuclear modification factors of the Y(2S) mesons as a function of collision centrality.

Confidence intervals on the nuclear modification factors of the Y(2S) mesons as a function of collision centrality.

Nuclear modification factors of the Y(1S) and Y(2S) mesons in the integrated pT, rapidity and centrality.

Nuclear modification factors of the Y(1S) and Y(2S) mesons in the integrated pT, rapidity and centrality.

Differential cross section of J/psi mesons from b hadrons as a function of dimuon pT in pp and PbPb collisions. The PbPb cross sections are normalised by TAA for direct comparison. Global uncertainties arise from the integrated luminosity uncertainty in pp collisions, and the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross section of prompt J/psi mesons as a function of dimuon rapidty in pp and PbPb collisions. The PbPb cross sections are normalised by TAA for direct comparison. Global uncertainties arise from the integrated luminosity uncertainty in pp collisions, and the number of minimum bias events and TAA uncertainties for PbPb collisions.

Differential cross section of J/psi mesons from b hadrons as a function of dimuon rapidty in pp and PbPb collisions. The PbPb cross sections are normalised by TAA for direct comparison. Global uncertainties arise from the integrated luminosity uncertainty in pp collisions, and the number of minimum bias events and TAA uncertainties for PbPb collisions.

Nuclear modification factor of prompt J/psi mesons as a function of dimuon rapidity. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of collision centrality. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of dimuon pT. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of dimuon pT in the mid-rapidity interval. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of dimuon pT in the most forward rapidity interval. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of collision centrality in the mid- and most forward rapidity intervals. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of pT in three centrality bins. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of collision centrality in the most forward rapidity range at moderate and high pT. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of collision centrality in the mid-rapidity range. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt psi(2S) mesons as a function of collision centrality in the mid-rapidity range. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

95% CL intervals on the nuclear modification factor of prompt psi(2S) mesons as a function of collision centrality in the mid-rapidity range. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt J/psi and psi(2S) mesons as a function of dimuon pT in the mid-rapidity range. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of collision centrality in the forward rapidity range. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt psi(2S) mesons as a function of collision centrality in the forward rapidity range. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

95% CL intervals on the nuclear modification factor of prompt psi(2S) mesons as a function of collision centrality in the forward rapidity range. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of prompt J/psi mesons as a function of dimuon pT in the forward rapidity range. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of prompt psi(2S) mesons as a function of dimuon pT in the forward rapidity range. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

95% CL intervals on the nuclear modification factor of prompt psi(2S) mesons as a function of dimuon pT in the forward rapidity range. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of dimuon rapidity. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of collision centrality. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of dimuon pT. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of dimuon pT in the mid-rapidity interval. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of dimuon pT in the most forward rapidity interval. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of collision centrality in the mid- and most forward rapidity intervals. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of pT in three centrality bins. Global uncertainties arise from the integrated luminosity uncertainty of the pp data, and the number of minimum bias events and TAA uncertainties of the PbPb data.

Nuclear modification factor of nonprompt J/psi mesons as a function of collision centrality in the most forward rapidity range at moderate and high pT. Global uncertainties arise from the full pp statistical and systematic uncertainties (including integrated luminosity uncertainty), and the number of minimum bias events uncertainty of the PbPb data.

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

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

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

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

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

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

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

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

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

Values of $g_1^{NS}$ for the COMPASS data in ($x$, $Q^2$) bins.

Parameters of two solutions of the QCD fit. The assumed functional form at $Q^2_0 = 1 ($GeV$/c)^2$ is $q(x) = \eta x^\alpha (1-x)^\beta (1 + \gamma x)/N$ where N is the normalisation. The parametrisation is used for the singlet and gluon helicity distribution and the two non-singlet distributions $\Delta q_3$ and $\Delta q_8$.

Normalisations of the different data sets used in the QCD fit.

Differential cross section for Y(2S) states as a function of their transverse momentum and per unit of rapidity in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(3S) states as a function of their transverse momentum and per unit of rapidity in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(3S) states as a function of their transverse momentum and per unit of rapidity in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(1S) states as a function of their transverse momentum and per unit of rapidity in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(1S) states as a function of their transverse momentum and per unit of rapidity in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(2S) states as a function of their transverse momentum and per unit of rapidity in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(2S) states as a function of their transverse momentum and per unit of rapidity in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(1S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(1S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(2S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(2S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(3S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(3S) states as a function of their rapidity and integrated over transverse momentum in pp collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 3.7% are not displayed.

Differential cross section for Y(1S) states as a function of their rapidity and integrated over transverse momentum in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(1S) states as a function of their rapidity and integrated over transverse momentum in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(2S) states as a function of their rapidity and integrated over transverse momentum in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Differential cross section for Y(2S) states as a function of their rapidity and integrated over transverse momentum in PbPb collisions. Statistical (systematic) uncertainties are displayed as error bars (boxes). Global relative uncertainties of 6.5% are not displayed.

Nuclear modification factor for Υ(1S) states in PbPb collisions as a function of pT. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(1S) states in PbPb collisions as a function of pT. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(2S) states in PbPb collisions as a function of pT. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(2S) states in PbPb collisions as a function of pT. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(1S) states in PbPb collisions as a function of |y|. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(1S) states in PbPb collisions as a function of |y|. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(2S) states in PbPb collisions as a function of |y|. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factor for Υ(2S) states in PbPb collisions as a function of |y|. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainty (7.5%) is displayed as a grey box at unity.

Nuclear modification factors for Υ(1S) meson production in PbPb collisions, as a function of centrality, displayed as the average number of participating nucleons. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainties from the PbPb data (3.2%) or from the pp reference (6.3%) are displayed at unity as empty and filled red black boxes, respectively.

Nuclear modification factors for Υ(1S) meson production in PbPb collisions, as a function of centrality, displayed as the average number of participating nucleons. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainties from the PbPb data (3.2%) or from the pp reference (6.3%) are displayed at unity as empty and filled red black boxes, respectively.

Nuclear modification factors for Υ(2S) meson production in PbPb collisions, as a function of centrality, displayed as the average number of participating nucleons. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainties from the PbPb data (3.2%) or from the pp reference (6.9%) are displayed at unity as empty and filled black black boxes, respectively.

Nuclear modification factors for Υ(2S) meson production in PbPb collisions, as a function of centrality, displayed as the average number of participating nucleons. Statistical (systematic) uncertainties are displayed as error bars (boxes), while the global (fully correlated) uncertainties from the PbPb data (3.2%) or from the pp reference (6.9%) are displayed at unity as empty and filled black black boxes, respectively.

Nuclear modification factors for Υ(1S), Υ(2S) and Υ(3S) meson production in PbPb collisions, integrated over centrality. For the Υ(3S), the upper limit derived on the nuclear modification factor is given.

Nuclear modification factors for Υ(1S), Υ(2S) and Υ(3S) meson production in PbPb collisions, integrated over centrality. For the Υ(3S), the upper limit derived on the nuclear modification factor is given.

95% CL intervals on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of pT, for the forward rapidity analysis bin.

Double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the midrapidity analysis bin.

95% CL intervals on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the midrapidity analysis bin.

Double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the forward rapidity analysis bin.

95% CL intervals on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, as a function of centrality, for the forward rapidity analysis bin.

Double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, integrated over centrality, for the midrapidity and forward rapidity analysis bins.

95% CL interval on the double ratio of measured yields, $(N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{\mathrm{PbPb}} / (N_{\psi\mathrm{(2S)}} / N_{J/\psi})_{pp}$, integrated over centrality, for the midrapidity and forward rapidity analysis bins.

Measured integrated cross sections for the $Z\gamma\gamma$ process for neutrino final state at $\sqrt{s} = 8$ TeV in the extended fiducial regions defined in the paper, table 5. The parton-to-particle correction factors are also shown, which are defined as the ratio of the cross sections at parton-level to the cross sections at particle-level.

The differential cross sections as a function of $E_{\mathrm{T}}^{\gamma}$ for the $pp \rightarrow \ell^{+}\ell^{-}\gamma$ process in the inclusive $N_{\mathrm{jets}} \geq 0$ extended fiducial region. The parton-to-particle correction factors are also shown, which are defined as the ratio of the cross sections at parton-level to the cross sections at particle-level.

The differential cross sections as a function of $E_{\mathrm{T}}^{\gamma}$ for the $pp \rightarrow \ell^{+}\ell^{-}\gamma$ process in the exclusive $N_{\mathrm{jets}} = 0$ extended fiducial region. The parton-to-particle correction factors are also shown, which are defined as the ratio of the cross sections at parton-level to the cross sections at particle-level.

The differential cross sections as a function of $E_{\mathrm{T}}^{\gamma}$ for the $pp \rightarrow \nu\bar{\nu}\gamma$ process in the inclusive $N_{\mathrm{jets}} \geq 0$ extended fiducial region. The parton-to-particle correction factors are also shown, which are defined as the ratio of the cross sections at parton-level to the cross sections at particle-level.

The differential cross sections as a function of $E_{\mathrm{T}}^{\gamma}$ for the $pp \rightarrow \nu\bar{\nu}\gamma$ process in the exclusive $N_{\mathrm{jets}} = 0$ extended fiducial region. The parton-to-particle correction factors are also shown, which are defined as the ratio of the cross sections at parton-level to the cross sections at particle-level.

The differential cross sections as a function of $m_{\ell^{+}\ell^{-}\gamma}$ for the $pp \rightarrow \ell^{+}\ell^{-}\gamma$ process in the inclusive $N_{\mathrm{jets}} \geq 0$ extended fiducial region. The parton-to-particle correction factors are also shown, which are defined as the ratio of the cross sections at parton-level to the cross sections at particle-level.

The differential cross sections as a function of $m_{\ell^{+}\ell^{-}\gamma}$ for the $pp \rightarrow \ell^{+}\ell^{-}\gamma$ process in the exclusive $N_{\mathrm{jets}} = 0$ extended fiducial region. The parton-to-particle correction factors are also shown, which are defined as the ratio of the cross sections at parton-level to the cross sections at particle-level.

The cross sections as a function of $N_{\mathrm{jets}}$ for the $pp \rightarrow \ell^{+}\ell^{-}\gamma$ process in the extended fiducial region. The parton-to-particle correction factors are also shown, which are defined as the ratio of the cross sections at parton-level to the cross sections at particle-level.

Measured fiducial cross section in all $\ell^+\ell^-\ell'^+\ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding luminosity uncertainty, the second is the luminosity uncertainty.

Cross section extrapolated to total phase space and all Z boson decays. The first systematic uncertainty is the combined systematic uncertainty excluding luminosity uncertainty, the second is the luminosity uncertainty.