Observed and expected 95 % CL exclusion limits on $\tan\beta$ as a function of $m_{H^{\pm}}$, shown in the context of the mh125 scenario, for $m_{H^{\pm}}>150$ GeV and $(1 \leq \tan\beta \leq 60)$. The surrounding shaded bands correspond to the 1$\sigma$ and 2$\sigma$ confidence intervals around the expected limit.

Observed and expected 95 % CL exclusions on $\tan\beta$ as a function of $m_{H^{\pm}}$ derived from exclusion limits on $\mathrm{\cal{B}}(t\to bH^+)\times \mathrm{\cal{B}}(H^+ \to \tau \nu)$, for Type I 2HDM low mass scenario. The surrounding shaded bands correspond to the 1$\sigma$ and 2$\sigma$ confidence intervals around the expected limit.

Observed and expected 95 % CL exclusions on $\tan\beta$ as a function of $m_{H^{\pm}}$ derived from exclusion limits on $\mathrm{\cal{B}}(t\to bH^+)\times \mathrm{\cal{B}}(H^+ \to \tau \nu)$, for Lepton Specific 2HDM low mass scenario. The surrounding shaded bands correspond to the 1$\sigma$ and 2$\sigma$ confidence intervals around the expected limit.

The observed and expected ($\pm2\sigma$) upper limits at 95% CL on the branching ratio for $H\rightarrow aa\rightarrow \gamma\gamma\tau\tau$ as a function of the resonance mass hypothesis $m_{a}$.

The polynomial parameterisation of the total signal selection efficiency as a function of $m_{a}$, defined as the fraction of the reconstructed and selected events out of the total number of events in the signal simulated samples.

Relative systematic uncertainties in the averaged cross-section for various differential distributions in the (a, b) inclusive and (c, d) collinear phase spaces. The upper solid line shows the total uncertainty in the measured cross-section in data, and includes correlations between the systematic components. The 'Others' category contains sub-dominant uncertainties arising from missing transverse momentum reconstruction and the jet-to-vertex fraction uncertainties.

Differential cross-section measurement in the inclusive phase space as function of (a) the minimum angular separation between the lepton and any jet with transverse momentum greater than 100 GeV (ΔR _min(ℓ,jet¹⁰⁰_i)), and (b) the ratio of W p_T to the closest-jet p_T (p_T^ℓν / p_T^{closest jet}). The measured cross-sections in data and their statistical uncertainty are shown by the solid dots and error bars, while the shaded band shows the systematic and statistical uncertainties added in quadrature. Errors bars on the theory prediction include theoretical uncertainties as discussed in the text. The bins around p_T^ℓν / p_T^{closest jet} equal to 1 are merged to be insensitive to the singularity that exists in the fixed-order MCFM calculation. Additionally, in (b), the central two bins from the MCFM prediction are merged due to a singularity in fixed-order calculations, which requires p_T resummation of the W + jets system. The two lower panels show the effect of NLO EW_virt as computed from Sherpa 2.2.11.

Differential cross-section measurement in the inclusive phase space as function of (a) the minimum angular separation between the lepton and any jet with transverse momentum greater than 100 GeV (ΔR _min(ℓ,jet¹⁰⁰_i)), and (b) the ratio of W p_T to the closest-jet p_T (p_T^ℓν / p_T^{closest jet}). The measured cross-sections in data and their statistical uncertainty are shown by the solid dots and error bars, while the shaded band shows the systematic and statistical uncertainties added in quadrature. Errors bars on the theory prediction include theoretical uncertainties as discussed in the text. The bins around p_T^ℓν / p_T^{closest jet} equal to 1 are merged to be insensitive to the singularity that exists in the fixed-order MCFM calculation. Additionally, in (b), the central two bins from the MCFM prediction are merged due to a singularity in fixed-order calculations, which requires p_T resummation of the W + jets system. The two lower panels show the effect of NLO EW_virt as computed from Sherpa 2.2.11.

Differential cross-section as a function of the invariant mass of the leading two jets (m_jj) in the inclusive, 2-jet phase space. The measured cross-sections in data and their statistical uncertainty are shown by the solid dots and error bars, while the shaded band shows the systematic and statistical uncertainties added in quadrature. Errors bars on the theory prediction include theoretical uncertainties as discussed in the text. The two lower panels show the effect of NLO EW_virt as computed from Sherpa 2.2.11.

Differential distributions at reconstruction level in the (a, c) electron or (b, d) muon channel for (a, b) inclusive and (c, d) collinear signal regions after the application of the background normalisation factors. The signal process is stacked above all background predictions. The bottom panel shows the ratio of the data to the total signal plus background prediction. The shaded band includes statistical and systematic uncertainties from signal and background processes added in quadrature.

Differential distributions at reconstruction level in the (a, c) electron or (b, d) muon channel for (a, b) inclusive and (c, d) collinear signal regions after the application of the background normalisation factors. The signal process is stacked above all background predictions. The bottom panel shows the ratio of the data to the total signal plus background prediction. The shaded band includes statistical and systematic uncertainties from signal and background processes added in quadrature.

Relative systematic uncertainties in the averaged cross-section for various differential distributions in the (a, b) inclusive and (c, d) collinear phase spaces. The upper solid line shows the total uncertainty in the measured cross-section in data, and includes correlations between the systematic components. The 'Others' category contains sub-dominant uncertainties arising from missing transverse momentum reconstruction and the jet-to-vertex fraction uncertainties.

Relative systematic uncertainties in the averaged cross-section for various differential distributions in the (a, b) inclusive and (c, d) collinear phase spaces. The upper solid line shows the total uncertainty in the measured cross-section in data, and includes correlations between the systematic components. The 'Others' category contains sub-dominant uncertainties arising from missing transverse momentum reconstruction and the jet-to-vertex fraction uncertainties.

Differential cross-section as a function of (a) the leading p_T ^jet, (b) the inclusive jet multiplicity, (c) S_T, and (d) p^ℓν_T in the collinear phase-space. The measured cross-sections in data and their statistical uncertainty are shown by the solid dots and error bars, while the shaded band shows the systematic and statistical uncertainties added in quadrature. Errors bars on the theory prediction include theoretical uncertainties as discussed in the text. The two lower panels show the effect of NLO EW_virt as computed from Sherpa 2.2.11.

Differential cross-section as a function of (a) the leading p_T ^jet, (b) the inclusive jet multiplicity, (c) S_T, and (d) p^ℓν_T in the collinear phase-space. The measured cross-sections in data and their statistical uncertainty are shown by the solid dots and error bars, while the shaded band shows the systematic and statistical uncertainties added in quadrature. Errors bars on the theory prediction include theoretical uncertainties as discussed in the text. The two lower panels show the effect of NLO EW_virt as computed from Sherpa 2.2.11.

Differential cross-section as a function of (a) the leading p_T ^jet, (b) the inclusive jet multiplicity, (c) S_T, and (d) p^ℓν_T in the collinear phase-space. The measured cross-sections in data and their statistical uncertainty are shown by the solid dots and error bars, while the shaded band shows the systematic and statistical uncertainties added in quadrature. Errors bars on the theory prediction include theoretical uncertainties as discussed in the text. The two lower panels show the effect of NLO EW_virt as computed from Sherpa 2.2.11.

Differential cross-section as a function of (a) the leading p_T ^jet, (b) the inclusive jet multiplicity, (c) S_T, and (d) p^ℓν_T in the collinear phase-space. The measured cross-sections in data and their statistical uncertainty are shown by the solid dots and error bars, while the shaded band shows the systematic and statistical uncertainties added in quadrature. Errors bars on the theory prediction include theoretical uncertainties as discussed in the text. The two lower panels show the effect of NLO EW_virt as computed from Sherpa 2.2.11.

Measured fiducial cross-sections in each signal region. The measured cross-sections in data and their total uncertainty are shown by the solid dots and error band. Various theoretical predictions are overlaid and compared with the data in the lower panel. Uncertainties in predictions accurate to NLO precision include QCD scale and PDF uncertainties, while predictions that add NLO EW_virt (electroweak) corrections show uncertainties only due to different electroweak combination schemes. Uncertainties in predictions accurate at fixed-order NNLO precision are derived primarily from MCFM, do not include PDF variations, with an additional extrapolation from the Sherpa QCD with EW_virt prediction to assess the effects of higher-order electroweak corrections. On predictions with EW_virt corrections, the total uncertainty is shown as a shaded region, while the electroweak correction is overlaid with an error bar line.

Differential distributions at reconstructed-level comparing data with the background model in the (a) tt̄ muon, (b) Z + jets muon, and (c) multijet electron control selections. All backgrounds and signal samples are stacked to produce the figures. The data points are shown as solid dots with their corresponding statistical uncertainty, while the total background model uncertainty is shown by the shaded band. Errors bars on the theory prediction include theoretical uncertainties as discussed in the text.

Differential distributions at reconstructed-level comparing data with the background model in the multijet electron validation region for (a) the transverse momentum of the leading jet and (b) the transverse momentum of the lepton–neutrino system. All backgrounds and signal samples are stacked to produce the figures. The data points are shown as solid dots with their corresponding statistical uncertainty, while the total background model uncertainty is shown by the shaded band. Errors bars on the theory prediction include theoretical uncertainties as discussed in the text.

Differential distributions at reconstructed-level comparing data with the background model in the multijet electron validation region for (a) the transverse momentum of the leading jet and (b) the transverse momentum of the lepton–neutrino system. All backgrounds and signal samples are stacked to produce the figures. The data points are shown as solid dots with their corresponding statistical uncertainty, while the total background model uncertainty is shown by the shaded band. Errors bars on the theory prediction include theoretical uncertainties as discussed in the text.

Differential distributions at reconstructed-level comparing data with the background model in the (a) tt̄ muon, (b) Z + jets muon, and (c) multijet electron control selections. All backgrounds and signal samples are stacked to produce the figures. The data points are shown as solid dots with their corresponding statistical uncertainty, while the total background model uncertainty is shown by the shaded band. Errors bars on the theory prediction include theoretical uncertainties as discussed in the text.

Differential distributions at reconstructed-level comparing data with the background model in the (a) tt̄ muon, (b) Z + jets muon, and (c) multijet electron control selections. All backgrounds and signal samples are stacked to produce the figures. The data points are shown as solid dots with their corresponding statistical uncertainty, while the total background model uncertainty is shown by the shaded band. Errors bars on the theory prediction include theoretical uncertainties as discussed in the text.

Response matrices from the Sherpa 2.2.11 event generator for (a) ΔR _min(ℓ,jet¹⁰⁰_i) in the inclusive electron selection, and (b) transverse momentum of the leading jet in the collinear muon selection. The sum of the entries in each reconstructed bin (column) are normalised to unity.

Response matrices from the Sherpa 2.2.11 event generator for (a) ΔR _min(ℓ,jet¹⁰⁰_i) in the inclusive electron selection, and (b) transverse momentum of the leading jet in the collinear muon selection. The sum of the entries in each reconstructed bin (column) are normalised to unity.

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the resolved qqbb signal region.

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the resolved lvbb signal region.

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the merged lvbb signal region

Product of acceptance and efficiency for pp->tbH(->Wh) as function of the charged Higgs boson mass for the merged qqbb signal region

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 5jex3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 5jex4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 6jin3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the low-purity signal regions of the resolved qqbb 6jin4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 5jex3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 5jex4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 6jin3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the high-purity signal regions of the resolved qqbb 6jin4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 5jex3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 5jex4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 6jin3bex event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

Distributions of the mWh observable in the signal regions of the resolved lvbb 6jin4bin event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The distributions are presented after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contribution assuming $m_{H^{\pm}}$ = 700 GeV, normalised to the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) of 0.064 pb, is shown as a dashed histogram.

The mWh distributions in the Low-NN-score SRs of the merged qqbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Low-NN-score SRs of the merged qqbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged qqbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged qqbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Low-NN-score SRs of the merged lvbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Low-NN-score SRs of the merged lvbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Medium-NN-score SRs of the merged lvbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the Medium-NN-score SRs of the merged lvbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged lvbb 0b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

The mWh distributions in the High-NN-score SRs of the merged lvbb 1b event categories. The term ‘Others’ summarises events from tHjb, tWh, tttt, and SM Vh production. The background prediction is shown after a background-only maximum-likelihood fit to data. The individual background uncertainty does not take into account the possible correlations between the nuisance parameter. The expected signal contributions assuming $m_{H^{\pm}}$ = 900 GeV and $m_{H^{\pm}}$ = 2000 GeV, normalised the expected limit of the cross-section times branching ratio ($\sigma_{sig}$ × B) values of 0.036 pb and 2.7 fb respectively, are shown as dashed histograms.

Cutflow table for a few selected signal samples after applying the selection requirements of the resolved analysis. At the initial cut stage, all events from the tbH ± → W ± h(→ b b) production process (including those from the fully hadronic decay channel) are considered.

Cutflow table for a few selected signal samples after applying the selection requirements of the merged analysis. At the initial cut stage, all events from the tbH ± → W ± h(→ bb) production process (including those from the fully hadronic decay channel) are considered.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM1.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM1.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM1.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM1.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM1.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM1.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM2.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM2.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM2.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM2.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM2.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM2.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM3.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM3.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM3.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM3.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM3.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-HM3.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp1.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp1.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp1.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp1.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp1.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp1.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp2.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp2.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp2.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp2.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp2.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp2.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp3.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp3.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp3.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp3.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp3.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp3.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp-1c.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to charm quarks and neutralinos when only considering SR-Comp-1c.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to top or charm quarks and 1 GeV neutralinos for different branching fraction assumptions.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to top or charm quarks and 1 GeV neutralinos for different branching fraction assumptions.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to top or charm quarks and 1 GeV neutralinos for different branching fraction assumptions.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to top or charm quarks and 1 GeV neutralinos for different branching fraction assumptions.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to top or charm quarks and 1 GeV neutralinos for different branching fraction assumptions.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to top or charm quarks and 1 GeV neutralinos for different branching fraction assumptions.

Expected exclusion limit at 95% CL for pair production of top squarks decaying to top or charm quarks and 200 GeV neutralinos for different branching fraction assumptions.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to top or charm quarks and 200 GeV neutralinos for different branching fraction assumptions.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to top or charm quarks and 200 GeV neutralinos for different branching fraction assumptions.

Observed exclusion limit at 95% CL for pair production of top squarks decaying to top or charm quarks and 200 GeV neutralinos for different branching fraction assumptions.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of top squarks decaying to top or charm quarks and 200 GeV neutralinos for different branching fraction assumptions.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of top squarks decaying to top or charm quarks and 200 GeV neutralinos for different branching fraction assumptions.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of up-type scalar LQs coupled to first generation leptons and second generation quarks.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of up-type scalar LQs coupled to first generation leptons and second generation quarks.

Expected exclusion limit at 95% CL for pair production of up-type scalar LQs coupled to first generation leptons and second generation quarks.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of up-type scalar LQs coupled to first generation leptons and second generation quarks.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of up-type scalar LQs coupled to first generation leptons and second generation quarks.

Observed exclusion limit at 95% CL for pair production of up-type scalar LQs coupled to first generation leptons and second generation quarks.

Expected exclusion limit $(-1\sigma)$ at 95% CL for pair production of up-type scalar LQs coupled to second generation leptons and quarks.

Expected exclusion limit $(+1\sigma)$ at 95% CL for pair production of up-type scalar LQs coupled to second generation leptons and quarks.

Expected exclusion limit at 95% CL for pair production of up-type scalar LQs coupled to second generation leptons and quarks.

Observed exclusion limit $(+1\sigma)$ at 95% CL for pair production of up-type scalar LQs coupled to second generation leptons and quarks.

Observed exclusion limit $(-1\sigma)$ at 95% CL for pair production of up-type scalar LQs coupled to second generation leptons and quarks.

Observed exclusion limit at 95% CL for pair production of up-type scalar LQs coupled to second generation leptons and quarks.

Observed upper limit on the cross-section at 95% CL for pair production of top squarks decaying into charm quarks and neutralinos for the fit with the best expected CL$_\mathrm{S}$.

Observed upper limit on the cross-section for pair production of up-type scalar LQs coupled to first generation leptons and second generation quarks.

Observed upper limit on the cross-section for pair production of up-type scalar LQs coupled to second generation leptons and quarks.

Expected upper limit on the cross-section $(-1\sigma)$ for $U(1)$ vector LQ models in the minimal coupling scenario with $\beta_{23}=1$, corresponding to a branching fraction $\mathcal{B}(\mathrm{vLQ}^\mathrm{u,min}_{23}\rightarrow c\nu) = 0.5$.

Expected upper limit on the cross-section $(+1\sigma)$ for $U(1)$ vector LQ models in the minimal coupling scenario with $\beta_{23}=1$, corresponding to a branching fraction $\mathcal{B}(\mathrm{vLQ}^\mathrm{u,min}_{23}\rightarrow c\nu) = 0.5$.

Expected upper limit on the cross-section for $U(1)$ vector LQ models in the minimal coupling scenario with $\beta_{23}=1$, corresponding to a branching fraction $\mathcal{B}(\mathrm{vLQ}^\mathrm{u,min}_{23}\rightarrow c\nu) = 0.5$.

Observed upper limit on the cross-section for $U(1)$ vector LQ models in the minimal coupling scenario with $\beta_{23}=1$, corresponding to a branching fraction $\mathcal{B}(\mathrm{vLQ}^\mathrm{u,min}_{23}\rightarrow c\nu) = 0.5$.

Expected upper limit on the cross-section $(-1\sigma)$ for $U(1)$ vector LQ models in the Yang-Mills coupling scenario with $\beta_{23}=1$, corresponding to a branching fraction $\mathcal{B}(\mathrm{vLQ}^\mathrm{u,YM}_{23}\rightarrow c\nu) = 0.5$.

Expected upper limit on the cross-section $(+1\sigma)$ for $U(1)$ vector LQ models in the Yang-Mills coupling scenario with $\beta_{23}=1$, corresponding to a branching fraction $\mathcal{B}(\mathrm{vLQ}^\mathrm{u,YM}_{23}\rightarrow c\nu) = 0.5$.

Expected upper limit on the cross-section for $U(1)$ vector LQ models in the Yang-Mills coupling scenario with $\beta_{23}=1$, corresponding to a branching fraction $\mathcal{B}(\mathrm{vLQ}^\mathrm{u,YM}_{23}\rightarrow c\nu) = 0.5$.

Observed upper limit on the cross-section for $U(1)$ vector LQ models in the Yang-Mills coupling scenario with $\beta_{23}=1$, corresponding to a branching fraction $\mathcal{B}(\mathrm{vLQ}^\mathrm{u,YM}_{23}\rightarrow c\nu) = 0.5$.

Distribution of $R_{\mathrm{ISR}}$ in SR-Comp-1c showing the data and the expected backgrounds, after simultaneously fitting all the control regions. The pre-fit yields for a representative $m(\tilde{t}_1,\tilde{\chi}_1^0)=(450,430)$ GeV signal is included. Processes with top quarks and multiple vector bosons are included in "Other". The right-most bin does not include the overflow entry but the $x$-axis range is chosen to include all observed data.

Distribution of $E_\mathrm{T}^\mathrm{miss}\mathrm{ Sig.}$ in SR-HM1 showing the data and the expected backgrounds, after simultaneously fitting all the control regions. The pre-fit yields for a representative $m(\tilde{t}_1,\tilde{\chi}_1^0)=(1000,1)$ GeV signal is included. Processes with top quarks and multiple vector bosons are included in "Other". The right-most bin does not include the overflow entry but the $x$-axis range is chosen to include all observed data.

Distribution of $m_{\mathrm{T}}^{\mathrm{min}}(c)$ in SR-HM2 showing the data and the expected backgrounds, after simultaneously fitting all the control regions. The pre-fit yields for a representative $m(\tilde{t}_1,\tilde{\chi}_1^0)=(750,450)$ GeV signal is included. Processes with top quarks and multiple vector bosons are included in "Other". The right-most bin does not include the overflow entry but the $x$-axis range is chosen to include all observed data.

Distribution of $R_{\mathrm{ISR}}$ in SR-Comp1 showing the data and the expected backgrounds, after simultaneously fitting all the control regions. The pre-fit yields for a representative $m(\tilde{t}_1,\tilde{\chi}_1^0)=(600,550)$ GeV signal is included. Processes with top quarks and multiple vector bosons are included in "Other". The right-most bin does not include the overflow entry but the $x$-axis range is chosen to include all observed data.

Distribution of $R_{\mathrm{ISR}}$ in SR-Comp2 showing the data and the expected backgrounds, after simultaneously fitting all the control regions. The pre-fit yields for a representative $m(\tilde{t}_1,\tilde{\chi}_1^0)=(550,470)$ GeV signal is included. Processes with top quarks and multiple vector bosons are included in "Other". The right-most bin does not include the overflow entry but the $x$-axis range is chosen to include all observed data.

Distribution of $R_{\mathrm{ISR}}$ in SR-Comp3 showing the data and the expected backgrounds, after simultaneously fitting all the control regions. The pre-fit yields for a representative $m(\tilde{t}_1,\tilde{\chi}_1^0)=(550,375)$ GeV signal is included. Processes with top quarks and multiple vector bosons are included in "Other". The right-most bin does not include the overflow entry but the $x$-axis range is chosen to include all observed data.

Acceptance across the SUSY signal grid for SR-Comp-1c.

Acceptance across the SUSY signal grid for SR-HM1.

Acceptance across the SUSY signal grid for SR-HM2.

Acceptance across the SUSY signal grid for SR-HM3.

Acceptance across the SUSY signal grid for SR-HM-Disc.

Acceptance across the SUSY signal grid for SR-Comp1.

Acceptance across the SUSY signal grid for SR-Comp2.

Acceptance across the SUSY signal grid for SR-Comp3.

Acceptance across the $\mathrm{LQ}^\mathrm{u}_{21}$ signal grid for SR-HM1.

Acceptance across the $\mathrm{LQ}^\mathrm{u}_{21}$ signal grid for SR-HM2.

Acceptance across the $\mathrm{LQ}^\mathrm{u}_{21}$ signal grid for SR-HM3.

Acceptance across the $\mathrm{LQ}^\mathrm{u}_{21}$ signal grid for SR-HM-Disc.

Acceptance across the $\mathrm{LQ}^\mathrm{u}_{22}$ signal grid for SR-HM1.

Acceptance across the $\mathrm{LQ}^\mathrm{u}_{22}$ signal grid for SR-HM2.

Acceptance across the $\mathrm{LQ}^\mathrm{u}_{22}$ signal grid for SR-HM3.

Acceptance across the $\mathrm{LQ}^\mathrm{u}_{22}$ signal grid for SR-HM-Disc.

Acceptance for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM1.

Acceptance for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM2.

Acceptance for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM3.

Acceptance for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM-Disc.

Acceptance for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM1.

Acceptance for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM2.

Acceptance for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM3.

Acceptance for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM-Disc.

Efficiency for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM1. Due to the limited statistical power of the signal samples, efficiencies can fluctuate and can drop to zero.

Efficiency for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM2. Due to the limited statistical power of the signal samples, efficiencies can fluctuate and can drop to zero.

Efficiency for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM3. Due to the limited statistical power of the signal samples, efficiencies can fluctuate and can drop to zero.

Efficiency for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM-Disc. Due to the limited statistical power of the signal samples, efficiencies can fluctuate and can drop to zero.

Efficiency for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM1. Due to the limited statistical power of the signal samples, efficiencies can fluctuate and can drop to zero.

Efficiency for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM2. Due to the limited statistical power of the signal samples, efficiencies can fluctuate and can drop to zero.

Efficiency for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM3. Due to the limited statistical power of the signal samples, efficiencies can fluctuate and can drop to zero.

Efficiency for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM-Disc. Due to the limited statistical power of the signal samples, efficiencies can fluctuate and can drop to zero.

Acceptance times efficiency across the SUSY signal grid for SR-Comp-1c.

Acceptance times efficiency across the SUSY signal grid for SR-HM1.

Acceptance times efficiency across the SUSY signal grid for SR-HM2.

Acceptance times efficiency across the SUSY signal grid for SR-HM3.

Acceptance times efficiency across the SUSY signal grid for SR-HM-Disc.

Acceptance times efficiency across the SUSY signal grid for SR-Comp1.

Acceptance times efficiency across the SUSY signal grid for SR-Comp2.

Acceptance times efficiency across the SUSY signal grid for SR-Comp3.

Acceptance times efficiency across the $\mathrm{LQ}^\mathrm{u}_{21}$ signal grid for SR-HM1.

Acceptance times efficiency across the $\mathrm{LQ}^\mathrm{u}_{21}$ signal grid for SR-HM2.

Acceptance times efficiency across the $\mathrm{LQ}^\mathrm{u}_{21}$ signal grid for SR-HM3.

Acceptance times efficiency across the $\mathrm{LQ}^\mathrm{u}_{21}$ signal grid for SR-HM-Disc.

Acceptance times efficiency across the $\mathrm{LQ}^\mathrm{u}_{22}$ signal grid for SR-HM1.

Acceptance times efficiency across the $\mathrm{LQ}^\mathrm{u}_{22}$ signal grid for SR-HM2.

Acceptance times efficiency across the $\mathrm{LQ}^\mathrm{u}_{22}$ signal grid for SR-HM3.

Acceptance times efficiency across the $\mathrm{LQ}^\mathrm{u}_{22}$ signal grid for SR-HM-Disc.

Acceptance times efficiency for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM1.

Acceptance times efficiency for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM2.

Acceptance times efficiency for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM3.

Acceptance times efficiency for the $\mathrm{vLQ}^\mathrm{u,min}_{23}$ signal model in SR-HM-Disc.

Acceptance times efficiency for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM1.

Acceptance times efficiency for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM2.

Acceptance times efficiency for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM3.

Acceptance times efficiency for the $\mathrm{vLQ}^\mathrm{u,YM}_{23}$ signal model in SR-HM-Disc.

Cutflow table for SR-Comp-1c using the $m(\tilde{t}_{1},\tilde{\chi}_{1}^{0}) = (450,430)$ GeV signal point as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-Comp1 using the $m(\tilde{t}_{1},\tilde{\chi}_{1}^{0}) = (600,550)$ GeV signal point as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-Comp2 using the $m(\tilde{t}_{1},\tilde{\chi}_{1}^{0}) = (550,470)$ GeV signal point as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-Comp3 using the $m(\tilde{t}_{1},\tilde{\chi}_{1}^{0}) = (550,375)$ GeV signal point as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-HM1 using the $m(\tilde{t}_{1},\tilde{\chi}_{1}^{0}) = (1000,1)$ GeV signal point as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>40$ GeV, however, the $E_\mathrm{T}^\mathrm{miss}$ is calculated with jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-HM2 using the $m(\tilde{t}_{1},\tilde{\chi}_{1}^{0}) = (750,450)$ GeV signal point as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>40$ GeV, however, the $E_\mathrm{T}^\mathrm{miss}$ is calculated with jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-HM3 using the $m(\tilde{t}_{1},\tilde{\chi}_{1}^{0}) = (750,450)$ GeV signal point as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>40$ GeV, however, the $E_\mathrm{T}^\mathrm{miss}$ is calculated with jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-HM-Disc using the $m(\tilde{t}_{1},\tilde{\chi}_{1}^{0}) = (750,450)$ GeV signal point as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>40$ GeV, however, the $E_\mathrm{T}^\mathrm{miss}$ is calculated with jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-HM1 using the $m(\mathrm{LQ}_{21}^\mathrm{u}) = 600$ GeV signal point, with a 50% branching fraction to $c\nu$, as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>40$ GeV, however, the $E_\mathrm{T}^\mathrm{miss}$ is calculated with jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-HM2 using the $m(\mathrm{LQ}_{21}^\mathrm{u}) = 600$ GeV signal point, with a 50% branching fraction to $c\nu$, as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>40$ GeV, however, the $E_\mathrm{T}^\mathrm{miss}$ is calculated with jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-HM3 using the $m(\mathrm{LQ}_{21}^\mathrm{u}) = 600$ GeV signal point, with a 50% branching fraction to $c\nu$, as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>40$ GeV, however, the $E_\mathrm{T}^\mathrm{miss}$ is calculated with jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Cutflow table for SR-HM-Disc using the $m(\mathrm{LQ}_{21}^\mathrm{u}) = 600$ GeV signal point, with a 50% branching fraction to $c\nu$, as the benchmark. All jet-related variables use jets with $p_\mathrm{T}>40$ GeV, however, the $E_\mathrm{T}^\mathrm{miss}$ is calculated with jets with $p_\mathrm{T}>20$ GeV. The "Preselection" requirement includes: $E_\mathrm{T}^\mathrm{miss}>250$ GeV; events containing electrons or muons are vetoed; at least two jets must be present of which at least one must be identified as a $c$-tagged jet; events with jets that have been identified as $b$-tagged are discarded; the minimum azimuthal angular ($\phi$) distance between the (up to) four leading jets and the $\vec{p}_\mathrm{T}^\mathrm{miss}$ is required to be greater than 0.4.

Distribution of the reconstructed $\tau\tau$ invariant mass ($m_{\tau\tau}$) for all events in the VBF_1 signal region for $p_{\text{T}}^{H}>200$ GeV. The observed Higgs boson signal corresponds to $(\sigma\times B)/(\sigma\times B)_{\text{SM}}\,=\,0.99$. Entries with values above the $x$-axis range are shown in the last bin of each distribution. The prediction for each sample is determined from the likelihood fit performed to measure the total $pp\rightarrow H\rightarrow\tau\tau$ cross-section.

Distribution of the reconstructed $\tau\tau$ invariant mass ($m_{\tau\tau}$) for all events in the VBF_0 signal region for $m_{jj}<1$ TeV. The observed Higgs boson signal corresponds to $(\sigma\times B)/(\sigma\times B)_{\text{SM}}\,=\,0.99$. Entries with values above the $x$-axis range are shown in the last bin of each distribution. The prediction for each sample is determined from the likelihood fit performed to measure the total $pp\rightarrow H\rightarrow \tau\tau$ cross-section.

Distribution of the reconstructed $\tau\tau$ invariant mass ($m_{\tau\tau}$) for all events in the VBF_0 signal region for $m_{jj}>1$ TeV. The observed Higgs boson signal corresponds to $(\sigma\times B)/(\sigma\times B)_{\text{SM}}\,=\,0.99$. Entries with values above the $x$-axis range are shown in the last bin of each distribution. The prediction for each sample is determined from the likelihood fit performed to measure the total $pp\rightarrow H \rightarrow \tau\tau$ cross-section.

Distribution of the reconstructed $\tau\tau$ invariant mass ($m_{\tau\tau}$) for all events in the VBF_1 signal region for $m_{jj}<1$ TeV. The observed Higgs boson signal corresponds to $(\sigma\times B)/(\sigma\times B)_{\text{SM}}\,=\,0.99$. Entries with values above the $x$-axis range are shown in the last bin of each distribution. The prediction for each sample is determined from the likelihood fit performed to measure the total $pp\rightarrow \rightarrow \tau\tau$ cross-section.

Distribution of the reconstructed $\tau\tau$ invariant mass ($m_{\tau\tau}$) for all events in the VBF_1 signal region for $m_{jj}>1$ TeV. The observed Higgs boson signal corresponds to $(\sigma\times B)/(\sigma\times B)_{\text{SM}}\,=\,0.99$. Entries with values above the $x$-axis range are shown in the last bin of each distribution. The prediction for each sample is determined from the likelihood fit performed to measure the total $pp\rightarrow H\rightarrow \tau\tau$ cross-section.

Comparison of the observed and expected event yields in the ttH categories. The window category contains events with $m_{\tau\tau} \in [100, 150]$ GeV, while events outside of this mass region are included in the sideband region. The observed Higgs boson signal corresponds to $(\sigma\times B)/(\sigma\times B)_{\text{SM}}\,=\,0.99$. The prediction for each sample is determined from the likelihood fit performed to measure the total $pp\rightarrow H\rightarrow\tau\tau$ cross-section.

The measured values of $\sigma_{H}\times B(H\rightarrow\tau\tau)$ relative to the SM expectations in the total (labelled as 'Combined') and per-production-mode fit. The total $\pm1\sigma$ uncertainty in each measurement is indicated by the error bar. The total systematic uncertainty in each measurement is also indicated.

The measured values for $\sigma_{H}\times B(H\rightarrow\tau\tau)$ relative to the SM expectations in the simplified template cross-section measurement. The total $\pm1\sigma$ uncertainty in each measurement is indicated by the error bar. The total systematic uncertainty in each measurement is also indicated.

Upper limits at 95% CL on simplified template cross-sections in the individual ttH $p_{\text{T}}^{H}$ bins relative to their SM expectation in the $H\rightarrow\tau\tau$ decay channel, derived using the CL$_\text{s}$ method. The observed upper limits, together with the expected limits both in the background-only hypothesis and in the SM hypothesis (dotted red lines) are included. In the case of the expected limits in the background-only hypothesis, one- and two-standard-deviation uncertainty bands are also shown.

Measured fiducial differential cross-sections as a function of $p_\text{T}(j_0)$. The measurements are compared with particle-level SM predictions from the POWHEG+PYTHIA8, POWHEG+HERWIG7 and MadGraph5_aMC@NLO+Pythia8 generators for the combined VBF, ggF, $VH$, and $t\bar{t}H$ Higgs boson production modes. The shaded box around each data point shows the statistical uncertainty, while the total uncertainty is indicated by the error bar. The contribution from the ggF, $VH$, and $t\bar{t}H$ production modes as predicted by POWHEG+PYTHIA8 is also shown.

Measured fiducial differential cross-sections as a function of $p_\text{T}^{H}$. The measurements are compared with particle-level SM predictions from the POWHEG+PYTHIA8, POWHEG+HERWIG7 and MadGraph5_aMC@NLO+Pythia8 generators for the combined VBF, ggF, $VH$, and $t\bar{t} H$ Higgs boson production modes. The shaded box around each data point shows the statistical uncertainty, while the total uncertainty is indicated by the error bar. The contribution from the ggF, $VH$, and $t\bar{t} H$ production modes as predicted by POWHEG+PYTHIA8 is also shown.

Measured fiducial differential cross-sections as a function of $\Delta\phi_{jj}^{\text{signed}}$. The measurements are compared with particle-level SM predictions from the POWHEG+PYTHIA8, POWHEG+HERWIG7 and MadGraph5_aMC@NLO+Pythia8 generators for the combined VBF, ggF, $VH$, and $t\bar{t}H$ Higgs boson production modes. The shaded box around each data point shows the statistical uncertainty, while the total uncertainty is indicated by the error bar. The contribution from the ggF, $VH$, and $t\bar{t} H$ production modes as predicted by POWHEG+PYTHIA8 is also shown.

Measured fiducial differential cross-sections as a function of $\Delta\phi_{jj}^{\text{signed}}$ vs $p_\text{T}^\text{H}$. The measurements are compared with particle-level SM predictions from the POWHEG+PYTHIA8, POWHEG+HERWIG7 and MadGraph5_aMC@NLO+Pythia8 generators for the combined VBF, ggF, $VH$, and $t\bar{t}H$ Higgs boson production modes. The shaded box around each data point shows the statistical uncertainty, while the total uncertainty is indicated by the error bar. The contribution from the ggF, $VH$, and $t\bar{t}H$ production modes as predicted by POWHEG+PYTHIA8 is also shown.

The unfolded data compared with the nominal SM predictions from POWHEG+PYTHIA8, as well as the predictions for several values of Wilson coefficients in the SMEFT, for $\Delta\phi_{jj}^{\text{signed}}$.

The unfolded data compared with the nominal SM predictions from POWHEG+PYTHIA8, as well as the predictions for several values of Wilson coefficients in the SMEFT, for $\Delta\phi_{jj}^{\text{signed}}$.

The unfolded data compared with the nominal SM predictions from POWHEG+PYTHIA8, as well as the predictions for several values of Wilson coefficients in the SMEFT, for $\Delta\phi_{jj}^{\text{signed}}$ vs $p_\text{T}^{H}$.

The unfolded data compared with the nominal SM predictions from POWHEG+PYTHIA8, as well as the predictions for several values of Wilson coefficients in the SMEFT, for $p_\text{T}^{H}$.

The 95% expected and observed confidence intervals for each of the six considered Wilson coefficients. Each coefficient is treated individually and both the linear and linear + quadratic models are considered when evaluating the sensitivity to the terms that enter at $\mathcal{O}(\Lambda^{-4})$. For the CP-even operators the $\Delta\phi_{jj}^{\text{signed}}$ distribution is used to extract the confidence interval. The limits are computed at a new-physics scale $\Lambda=1$ TeV.

The 95% expected and observed confidence intervals for each of the six considered Wilson coefficients. Each coefficient is treated individually and both the linear and linear + quadratic models are considered when evaluating the sensitivity to the terms that enter at $\mathcal{O}(\Lambda^{-4})$. For the CP-odd operators the $\Delta\phi_{jj}^{\text{signed}}$ vs $p_\text{T}^{H}$ distribution is used to extract the confidence interval. The limits are computed at a new-physics scale $\Lambda=1$ TeV.

Two-dimensional observed and expected limits at 68% and 95% confidence level obtained from $c_{HW}$ and $c_{HB}$ using only the SM-dimension-6 interference in the SMEFT basis. The limits are computed at a new-physics scale $\Lambda$ = 1 TeV. In this case, the $\Delta\phi_{jj}^{signed}$ distribution is used. IMPORTANT: HEPData requires the same number of rows for each column. However, each curve of on the plot has a different number of points. For that reason, the curves are padded with data of "0.0, 0.0". The "0.0, 0.0" points should be omitted in analysis.

Two-dimensional observed and expected limits at 68% and 95% confidence level obtained from $c_{H\tilde{W}}$ and $c_{H\tilde{B}}$ using only the SM-dimension-6 interference in the SMEFT basis. The limits are computed at a new-physics scale $\Lambda$ = 1 TeV. In this purely CP-odd case, the $\Delta\phi_{jj}^{signed}$ vs $p_{T}^{H}$ distribution is used. IMPORTANT: HEPData requires the same number of rows for each column. However, each curve of on the plot has a different number of points. For that reason, the curves are padded with data of "0.0, 0.0". The "0.0, 0.0" points should be omitted in analysis.

Two-dimensional observed and expected limits at 68% and 95% confidence level obtained from $c_{HW}$ and $c_{H\tilde{W}}$ using only the SM-dimension-6 interference in the SMEFT basis. The limits are computed at a new-physics scale $\Lambda$ = 1 TeV. In this case, the $\Delta\phi_{jj}^{signed}$ distribution is used. IMPORTANT: HEPData requires the same number of rows for each column. However, each curve of on the plot has a different number of points. For that reason, the curves are padded with data of "0.0, 0.0". The "0.0, 0.0" points should be omitted in analysis.

The expected values of $\sigma_{H}\times B(H\rightarrow\tau\tau)$ relative to the SM expectations in the total (labelled as 'Combined') and per-production-mode fit. The total $\pm1\sigma$ uncertainty in each measurement is indicated by the error bar. The total systematic uncertainty in each measurement is also indicated.

The expected values for $\sigma_{H}\times B(H\rightarrow\tau\tau)$ relative to the SM expectations in the simplified template cross-section measurement. The total $\pm1\sigma$ uncertainty in each measurement is indicated by the error bar. The total systematic uncertainty in each measurement is also indicated.

The measured values for $\sigma_{H}\times B(H\rightarrow\tau\tau)$ in the simplified template cross-section measurement. The total $\pm1\sigma$ uncertainty in the measurement is indicated by the black error bars, with the individual contribution from the statistical uncertainty in blue.

The measured correlations between differential fiducial cross-sections, $\sigma$, and the normalisation factors, $\alpha$, for the $Z(\to\tau\tau) + \text{jets}$ and top backgrounds for $p_\text{T}(j_0)$. The bins denoted as $\sigma_{1,2,3,4}$ correspond to the measurements in the four bins listed elsewhere for each observable.

The measured correlations between differential fiducial cross-sections, $\sigma$, and the normalisation factors, $\alpha$, for the $Z(\to\tau\tau) + \text{jets}$ and top backgrounds for $p_\text{T}^{H}}$. The bins denoted as $\sigma_{1,2,3,4}$ correspond to the measurements in the four bins listed elsewhere for each observable.

The measured correlations between differential fiducial cross-sections, $\sigma$, and the normalisation factors, $\alpha$, for the $Z(\to\tau\tau) + \text{jets}$ and top backgrounds for $\Delta\phi_{jj}^{\text{signed}}$. The bins denoted as $\sigma_{1,2,3,4}$ correspond to the measurements in the four bins listed elsewhere for each observable.

The measured correlations between differential fiducial cross-sections, $\sigma$, and the normalisation factors, $\alpha$, for the $Z(\to\tau\tau) + \text{jets}$ and top backgrounds for $\Delta\phi_{jj}^{\text{signed}}$ vs $p_\text{T}^{H}$. The bins denoted as $\sigma_{1,2,3,4}$ correspond to the measurements in the four bins listed elsewhere for each observable.

Fig. 6b: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Away and Toward regions after the subtraction of Number density $N_{\rm ch}$ and $\Sigma p_{\rm T}$ distributions in the transverse region for p-Pb collisions for $p_{\rm T} >$ 0.5 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. 7a: $<p_{\rm T}>$ as a function of $p_{\rm T}^{\rm trig}$ in Toward (left) and Away (right) regions after the subtraction of transverse region for pp collisions for $p_{\rm T} >$ 0.5 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. 7b: $<p_{\rm T}>$ as a function of $p_{\rm T}^{\rm trig}$ in Toward (left) and Away (right) regions after the subtraction of transverse region for p-Pb collisions for $p_{\rm T} >$ 0.5 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A1: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 0.15 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A2: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 1.0 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A3: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 0.15 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Fig. A4: Number density $N_{\rm ch}$ (left) and $\Sigma p_{\rm T}$ (right) distributions as a function of $p_{\rm T}^{\rm trig}$ in Transverse, Away, and Toward regions for $p_{\rm T} >$ 1.0 GeV/$c$. The shaded areas and the error bars around the data points represent the systematic and statistical uncertainties, respectively.

Differential cross section times dimuon branching fraction of Y(1S) as a function of $y^{Y}_{CM}$ in pPb collisions. The global uncertainty arises from the integrated luminosity uncertainty in pPb collisions.

Differential cross section times dimuon branching fraction of Y(2S) as a function of $y^{Y}_{CM}$ in pPb collisions. The global uncertainty arises from the integrated luminosity uncertainty in pPb collisions.

Differential cross section times dimuon branching fraction of Y(3S) as a function of $y^{Y}_{CM}$ in pPb collisions. The global uncertainty arises from the integrated luminosity uncertainty in pPb collisions.

Differential cross section times dimuon branching fraction of Y(1S) as a function of pT in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(2S) as a function of pT in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(3S) as a function of pT in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(1S) as a function of $|y^{Y}_{CM}|$ in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(2S) as a function of $|y^{Y}_{CM}|$ in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Differential cross section times dimuon branching fraction of Y(3S) as a function of $|y^{Y}_{CM}|$ in pp collisions. The global uncertainty arises from the integrated luminosity uncertainty in pp collisions.

Nuclear modification factor of Y(1S) as a function of pT. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) as a function of pT. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) as a function of pT. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(1S) as a function of $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) as a function of $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) as a function of $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(1S) at forward and backward $y^{Y}_{CM}$ for pT < 6 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) at forward and backward $y^{Y}_{CM}$ for pT < 6 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) at forward and backward $y^{Y}_{CM}$ for pT < 6 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(1S) at forward and backward $y^{Y}_{CM}$ for 6 < pT < 30 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(2S) at forward and backward $y^{Y}_{CM}$ for 6 < pT < 30 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

Nuclear modification factor of Y(3S) at forward and backward $y^{Y}_{CM}$ for 6 < pT < 30 GeV/c. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

RFB of Y(1S) versus $N^{|\eta_{lab}|<2.4}_{tracks}$.

RFB of Y(2S) versus $N^{|\eta_{lab}|<2.4}_{tracks}$.

RFB of Y(3S) versus $N^{|\eta_{lab}|<2.4}_{tracks}$.

RFB of Y(1S) versus $E^{|\eta_{lab}|>4}_{T}$.

RFB of Y(2S) versus $E^{|\eta_{lab}|>4}_{T}$.

RFB of Y(3S) versus $E^{|\eta_{lab}|>4}_{T}$.

Nuclear modification factor of Y(1S), Y(2S), and Y(3S) integrated over pT and $y^{Y}_{CM}$. The global uncertainty arises from the integrated luminosity uncertainties in pPb and pp collisions.

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.