Forward-backward asymmetry, $N_{\ell}(+\eta_{\textrm{lab}})/ N_{\ell}(-\eta_{\textrm{lab}})$, for the negative leptons.

Forward-backward asymmetry, $N_{\ell}(+\eta_{\textrm{lab}})/ N_{\ell}(-\eta_{\textrm{lab}})$, for the positive leptons.

Forward-backward asymmetry, $N_{\ell}(+\eta_{\textrm{lab}})/ N_{\ell}(-\eta_{\textrm{lab}})$, for the charged-summed W bosons, as a function of the lepton pseudorapidity.

Interpretation of results in the T1tttt model. The colored regions show the upper limits (95\% CL) on the production cross section for $pp\rightarrow \tilde{g}\tilde{g},\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1$ in the $m_{\tilde{g}}$-$m_{\tilde{\chi}^0_1}$ plane. The curve shows the observed limit on the corresponding SUSY particle masses obtained by comparing the excluded cross section with $-1\sigma$ theoretical cross sections.

Interpretation of results in the T1tttt model. The colored regions show the upper limits (95\% CL) on the production cross section for $pp\rightarrow \tilde{g}\tilde{g},\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1$ in the $m_{\tilde{g}}$-$m_{\tilde{\chi}^0_1}$ plane. The curve shows the expected limit on the corresponding SUSY particle masses obtained by comparing the excluded cross section with theoretical cross sections.

Interpretation of results in the T1tttt model. The colored regions show the upper limits (95\% CL) on the production cross section for $pp\rightarrow \tilde{g}\tilde{g},\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1$ in the $m_{\tilde{g}}$-$m_{\tilde{\chi}^0_1}$ plane. The curve shows the expected limit on the corresponding SUSY particle masses corresponding to a $+1\sigma$ variation of the experimental uncertainty.

Interpretation of results in the T1tttt model. The colored regions show the upper limits (95\% CL) on the production cross section for $pp\rightarrow \tilde{g}\tilde{g},\tilde{g}\rightarrow t\bar{t}\tilde{\chi}^0_1$ in the $m_{\tilde{g}}$-$m_{\tilde{\chi}^0_1}$ plane. The curve shows the expected limit on the corresponding SUSY particle masses corresponding to a $-1\sigma$ variation of the experimental uncertainty.

Observed excluded region (95\% CL), shown in blue, in the $m_{\tilde{g}}$-$m_{\tilde{\chi}^0_1}$ plane for a model combining T5tttt, gluino pair production, followed by gluino decay to an on-shell top squark, together with a model for direct top squark pair production. The top squarks decay via the two-body process $\tilde{t}\rightarrow t\tilde{\chi}^0_1$. The neutralino and top squark masses are related by the constraint $m_{\tilde{t}} = m_{\tilde{\chi}^0_1} + 175$ GeV.

Expected excluded region (95\% CL), shown in blue, in the $m_{\tilde{g}}$-$m_{\tilde{\chi}^0_1}$ plane for a model combining T5tttt, gluino pair production, followed by gluino decay to an on-shell top squark, together with a model for direct top squark pair production. The top squarks decay via the two-body process $\tilde{t}\rightarrow t\tilde{\chi}^0_1$. The neutralino and top squark masses are related by the constraint $m_{\tilde{t}} = m_{\tilde{\chi}^0_1} + 175$ GeV.

Signal yields for the charge scenarios considered in the analysis, as well as their associateds uncertainties.

Signal yields for two charge scenarios considered in the analysis, as well as their associated uncertainties.

Signal yields for the charge scenarios considered in the analysis, as well as their associateds uncertainties.

Signal yields for two charge scenarios considered in the analysis, as well as their associated uncertainties.

Signal yields for the charge scenarios considered in the analysis, as well as their associateds uncertainties.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the early 2016 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 6%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text. Comparing with respect to the binomial fit starting at $N_{\text{hits}}^{\text{low dE/dx}} = 2$, and not $N_{\text{hits}}^{\text{low dE/dx}} = 1$, is needed to account for the fact that early 2016 data is more strongly affected low dE/dx by instrumental effects that widen the N hits distribution.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the early 2016 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 6%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text. Comparing with respect to the binomial fit starting at $N_{\text{hits}}^{\text{low dE/dx}} = 2$, and not $N_{\text{hits}}^{\text{low dE/dx}} = 1$, is needed to account for the fact that early 2016 data is more strongly affected low dE/dx by instrumental effects that widen the N hits distribution.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the late 2016 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 78%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the late 2016 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 78%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the 2017 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 65%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the 2017 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 65%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the 2018 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 9%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the 2018 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 9%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Exclusion region (hatched) at 95% CL in the FCP charge-mass plane for the considered signal. The expected exclusion is shown with the associated 1 and 2 standard deviations $\sigma$ bands. Signal points at charges 0.9, 0.8, 2/3, 0.5, and 1/3 e are connected by straight lines to guide the eye. This is a conservative interpolation. Previous exclusions from CMS [Phys. Rev. D 87 (2013) 092008), JHEP 07 (2013) 122] as well as OPAL [Phys. Lett. B 572 (2003) 8] are given for comparison.

Exclusion region (hatched) at 95% CL in the FCP charge-mass plane for the considered signal. The expected exclusion is shown with the associated 1 and 2 standard deviations $\sigma$ bands. Signal points at charges 0.9, 0.8, 2/3, 0.5, and 1/3 e are connected by straight lines to guide the eye. This is a conservative interpolation. Previous exclusions from CMS [Phys. Rev. D 87 (2013) 092008), JHEP 07 (2013) 122] as well as OPAL [Phys. Lett. B 572 (2003) 8] are given for comparison.

Particle-level distributions of ungroomed AK4 pTD2 in 120 < PT < 150 GeV in the central dijet region.

Particle-level distributions of ungroomed AK4 thrust in 120 < PT < 150 GeV in the Z+jet region.

Particle-level distributions of ungroomed AK4 thrust in 120 < PT < 150 GeV in the central dijet region.

Particle-level distributions of ungroomed AK4 width in 120 < PT < 150 GeV in the Z+jet region.

Particle-level distributions of ungroomed AK4 width in 120 < PT < 150 GeV in the central dijet region.

Particle-level distributions of ungroomed AK4 LHA in 120 < PT < 150 GeV in the Z+jet region.

Particle-level distributions of ungroomed AK4 LHA in 120 < PT < 150 GeV in the central dijet region.

Particle-level distributions of ungroomed AK4 LHA in 408 < PT < 1500 GeV in the Z+jet region.

Particle-level distributions of ungroomed AK4 LHA in 1000 < PT < 4000 GeV in the central dijet region.

Particle-level distributions of ungroomed AK8 LHA in 120 < PT < 150 GeV in the Z+jet region.

Particle-level distributions of ungroomed AK8 LHA in 120 < PT < 150 GeV in the central dijet region.

Particle-level distributions of ungroomed AK4 LHA (charged-only) in 120 < PT < 150 GeV in the Z+jet region.

Particle-level distributions of ungroomed AK4 LHA (charged-only) in 120 < PT < 150 GeV in the central dijet region.

Particle-level distributions of groomed AK4 LHA in 120 < PT < 150 GeV in the Z+jet region.

Particle-level distributions of groomed AK4 LHA in 120 < PT < 150 GeV in the central dijet region.

Mean of ungroomed LHA for AK4 jets as a function of PT in the Z+jet region.

Mean of ungroomed LHA for AK4 jets as a function of PT in the central dijet region.

95% CL intervals on the double ratio of measured yields, $(N_{\varUpsilon(3S)} / N_{\varUpsilon(1S)})_{\mathrm{PbPb}} / (N_{\varUpsilon(3S)} / N_{\varUpsilon(1S)})_{pp}$, as a function of centrality, for upsilon $|y|<2.4$ and $p_T<30$GeV, and $p_{T}^{\mu}>4$GeV.

68% CL intervals on the double ratio of measured yields, $(N_{\varUpsilon(3S)} / N_{\varUpsilon(1S)})_{\mathrm{PbPb}} / (N_{\varUpsilon(3S)} / N_{\varUpsilon(1S)})_{pp}$, as a function of centrality, for upsilon $|y|<2.4$ and $p_T<30$GeV, and $p_{T}^{\mu}>4$GeV.

RMS of raw subtracted $p_T$ for jets reconstructed in PbPb collisions (in different centrality bins) with the anti-k$_T$ algorythm with a distance parameter $R=0.3$, in $|\eta|<0.2$.

Inclusive jet cross section in pp collisions, for anti-k$_T$ jets with a distance parameter $R=0.2, 0.3, 0.4$, reconstructed in $|\eta|<0.2$. The '3.7%' luminosity uncertainty is a fully correlated uncertainty, and not included in the point-to-point uncertainty.

Inclusive jet spectra in different centrality bins of PbPb collisions, for anti-k$_T$ jets of distance parameter $R=0.2$, reconstructed in $|\eta|<0.2$. The PbPb spectra in each case are scaled by TAA. The uncertainty on the TAA is considered as global uncertainty (common to all points for a given centrality), and it is not included in the point-to-point uncertainties.

Inclusive jet spectra in different centrality bins of PbPb collisions, for anti-k$_T$ jets of distance parameter $R=0.3$, reconstructed in $|\eta|<0.2$. The PbPb spectra in each case are scaled by TAA. The uncertainty on the TAA is considered as global uncertainty (common to all points for a given centrality), and it is not included in the point-to-point uncertainties.

Inclusive jet spectra in different centrality bins of PbPb collisions, for anti-k$_T$ jets of distance parameter $R=0.4$, reconstructed in $|\eta|<0.2$. The PbPb spectra in each case are scaled by TAA. The uncertainty on the TAA is considered as global uncertainty (common to all points for a given centrality), and it is not included in the point-to-point uncertainties.

Inclusive jet RAA as a function of jet $p_T$, for anti-k$_T$ jets with distance parameter $R=0.2, 0.3, 0.4$, with $|\eta|<2$, for 0-5\% PbPb centrality. The '5.5%' combined uncertainty on the pp luminosity and PbPb TAA is considered as global uncertainty (common to all points), and it is not included in the point-to-point uncertainties.

Inclusive jet RAA as a function of jet $p_T$, for anti-k$_T$ jets with distance parameter $R=0.2, 0.3, 0.4$, with $|\eta|<2$, for 5-10\% PbPb centrality. The '5.9%' combined uncertainty on the pp luminosity and PbPb TAA is considered as global uncertainty (common to all points), and it is not included in the point-to-point uncertainties.

Inclusive jet RAA as a function of jet $p_T$, for anti-k$_T$ jets with distance parameter $R=0.2, 0.3, 0.4$, with $|\eta|<2$, for 10-30\% PbPb centrality. The '6.9%' combined uncertainty on the pp luminosity and PbPb TAA is considered as global uncertainty (common to all points), and it is not included in the point-to-point uncertainties.

Inclusive jet RAA as a function of jet $p_T$, for anti-k$_T$ jets with distance parameter $R=0.2, 0.3, 0.4$, with $|\eta|<2$, for 30-50\% PbPb centrality. The '10%' combined uncertainty on the pp luminosity and PbPb TAA is considered as global uncertainty (common to all points), and it is not included in the point-to-point uncertainties.

Inclusive jet RAA as a function of jet $p_T$, for anti-k$_T$ jets with distance parameter $R=0.2, 0.3, 0.4$, with $|\eta|<2$, for 50-70\% PbPb centrality. The '15%' combined uncertainty on the pp luminosity and PbPb TAA is considered as global uncertainty (common to all points), and it is not included in the point-to-point uncertainties.

Inclusive jet RAA as a function of jet $p_T$, for anti-k$_T$ jets with distance parameter $R=0.2, 0.3, 0.4$, with $|\eta|<2$, for 70-90\% PbPb centrality. The '18%' combined uncertainty on the pp luminosity and PbPb TAA is considered as global uncertainty (common to all points), and it is not included in the point-to-point uncertainties.

Inclusive jet RAA as a function of average Npart, for anti-k$_T$ jets with distance parameter $R=0.2, 0.3, 0.4$, reconstructed in the range $80< p_T < 90$GeV/c and $|\eta|<2$.

Inclusive jet RAA as a function of average Npart, for anti-k$_T$ jets with distance parameter $R=0.2, 0.3, 0.4$, reconstructed in the range $130< p_T < 150$GeV/c and $|\eta|<2$.

Inclusive jet RAA as a function of jet $p_T$, for anti-k$_T$ jets with distance parameter $R=0.2, 0.3$, reconstructed in the range $|\eta|<2$, for 0-10\% PbPb centrality. The '5.7%' combined uncertainty on the pp luminosity and PbPb TAA is considered as global uncertainty (common to all points), and it is not included in the point-to-point uncertainties.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of rapidity.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of $p_{\rm T}$.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of centrality at high $p_{\rm T}$, 6.5$<p_{\rm T}<$30 GeV/c, in three different $|y|$ regions.

Prompt J/$\psi$ $R_{\rm AA}$ as a function of centrality, at forward rapidity, 1.6$<|y|<$2.4, for two different $p_{\rm T}$ regions.

Nonprompt J/$\psi$ $v_{2}$ as a function of $p_{\rm T}$.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of rapidity.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of $p_{\rm T}$.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

Nonprompt J/$\psi$ $R_{\rm AA}$ as a function of centrality. The average ${N}_{\rm part}$ values correspond to events flatly distributed across centrality.

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(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(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(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 Υ(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 Υ(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 Υ(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.

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.

Observed and expected likelihood scans $f_{\Lambda1}^{Z\gamma}\cos\phi_{\Lambda1}^{Z\gamma}$. See Section 2 of the paper for more details.