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.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Post-fit SR yields of the high-mass channel. The upper panel shows the observed number of events, as well the post-fit background predictions in each region. The bottom panel shows the ratio of the observed data and the total background prediction. The shaded areas correspond to the total statistical and systematic uncertainties obtained after the fit and described in Section 6.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2016 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 2.3%.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2017 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 3.7%.

Pre-fit data and background (reweighted $2b$) predictions for each $4b$ SR $E_\text{T}^\text{miss}$ and $m_\text{eff}$ bin of the low-mass channel for the 2018 data-taking period. The bottom panel shows the significance of any differences between the observed $4b$ data and the background prediction. The $1\sigma$ and $2\sigma$ bands are shown in green and yellow, respectively. All systematics are included except the background normalization, which is 1.8%.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. Results from a previous ATLAS search using 24.3-36.1 fb$^{-1}$ [13] are shown by the solid (observed) and dashed (expected) blue lines. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass and high-mass channels. The low-mass channel is used for $m_{\tilde{H}}<250$ GeV while the high-mass channel is used for $m_{\tilde{H}}\ge250$ GeV. The plot shows the 95% CL observed (solid) and expected (dashed) upper limits on $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})$, assuming the theory cross section for higgsino pair production. The higgsinos are assumed to decay as $\tilde{H}\rightarrow h + \tilde{G}$ or $\tilde{H}\rightarrow Z + \tilde{G}$. The phase space above the lines is excluded.

Exclusion limits of the low-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the low-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the high-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Exclusion limits of the high-mass channel. The plot shows the observed (solid) and expected (dashed) 95% CL upper limits on the cross section of higgsino pair production, assuming a higgsino decay branching ratio of $\mathcal{B}(\tilde{H}\rightarrow h + \tilde{G})=100\%$. The theory cross section and its uncertainty are shown by the solid red line and red shading. The bottom panel shows the ratio of the limits to the theory cross section. The phase space above the lines is excluded.

Results of the background-only fit in the low-mass channel discovery region SR_LM_150. Both pre-fit and post-fit values are shown.

Results of the background-only fit in the low-mass channel discovery region SR_LM_300. Both pre-fit and post-fit values are shown.

The experimental efficiency of the low-mass channel for the exclusion and discovery signal regions as a function of higgsino mass. The experimental efficiency is defined as the number of events passing the detector-level event selections divided by the number of events passing the event selections for a perfect detector. The denominator is obtained by implementing particle-level event selections that emulate the detector-level selections. This treats the lack of availability of $b$-jet triggers as an inefficiency.

The particle-level acceptance for the low-mass exclusion and discovery signal regions, shown as a function of higgsino mass. The acceptance is defined as the fraction of signal events passing the particle-level event selection that emulates the detector-level selection. The acceptance calculation considers only those signal events where both higgsinos decay to Higgs bosons.

The experimental efficiency of the high-mass channel discovery regions as a function of higgsino mass. For each higgsino mass, the efficiency is shown for the SR-1 region corresponding to the mass. For masses above 1100 GeV, SR-1-1100 is used. The experimental efficiency is defined as the number of events passing the detector-level event selections divided by the number of events passing the event selections for a perfect detector. The denominator is obtained by implementing particle-level event selections that emulate the detector-level selections. The efficiency calculation considers only those signal events where both higgsinos decay to Higgs bosons.

The particle-level acceptance for the high-mass signal regions, shown as a function of higgsino mass. For each higgsino mass, the acceptance is shown for the SR-1 region corresponding to the mass. For masses above 1100 GeV, SR-1-1100 is used. The acceptance is defined as the fraction of signal events passing the particle-level event selection that emulates the detector-level selection. The acceptance calculation considers only those signal events where both higgsinos decay to Higgs bosons.

Cutflow for the low-mass channel for a representative 130 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 150 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 250 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 600 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 700 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 800 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 900 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 1000 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the low-mass channel for a representative 1100 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. The $b$-jet cut requires 4 or more $b$-jets with $p_\text{T}>40$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$, with the availability of $b$-jet triggers lowering the luminosity to 126 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 250 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 600 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 700 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 800 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 900 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1000 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1100 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1200 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1300 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1400 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

Cutflow for the high-mass channel for a representative 1500 GeV signal. The preselection requires 4 or more jets with $p_\text{T}>25$ GeV and 2 or more $b$-jets with $p_\text{T}>25$ GeV. As the samples are generated with $\mathcal{B}(\tilde{H}\rightarrow h\tilde{G})$=50%, $\mathcal{B}(\tilde{H}\rightarrow Z\tilde{G})$=50% to allow for both decays to be studied, the $hh$ events selection is used to select the events where each of the higgsinos decays to a Higgs boson. Expected yields are normalized to a luminosity of 139 fb$^{-1}$. All selections are cumulative, with the exception of the SR cuts, which are each applied separately.

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_{\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.

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.

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.

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.

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.

Best fit values of $\kappa_{gZ}= \kappa_{g}\cdot\kappa_{Z} / \kappa_{H}$ and of the ratios of coupling modifiers, as defined in the most generic parameterisation described in the context of the $\kappa$ framework, from the combined analysis of the $\sqrt{s}$=7 and 8 TeV data. The results are shown for the combination of ATLAS and CMS and also separately for each experiment, together with their total uncertainties and their breakdown into the four components described in the text. The uncertainties in $\lambda_{tg}$ and $\lambda_{WZ}$, for which a negative solution is allowed, are calculated around the overall best fit value. The expected total uncertainties in the measurements are also shown. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Measured global signal strength $\mu$ and its total uncertainty, together with the breakdown of the uncertainty into its four components as defined in the text. The results are shown for the combination of ATLAS and CMS, and separately for each experiment. The expected uncertainty, with its breakdown, is also shown.

Measured signal strengths $\mu$ and their total uncertainties for different Higgs boson production processes. The results are shown for the combination of ATLAS and CMS, and separately for each experiment, for the combined $\sqrt{s}$=7 and 8 TeV data. The expected uncertainties in the measurements are also displayed. These results are obtained assuming that the Higgs boson branching fractions are the same as in the SM.

Measured signal strengths $\mu$ and their total uncertainties for different Higgs boson decay channels. The results are shown for the combination of ATLAS and CMS, and separately for each experiment, for the combined $\sqrt{s}$=7 and 8 TeV data. The expected uncertainties in the measurements are also displayed. These results are obtained assuming that the Higgs boson production process cross sections at $\sqrt{s}$ = 7 and 8 TeV are the same as in the SM.

Fit results allowing BSM loop couplings, assuming that $|\kappa_{V}| \le$ 1, where $\kappa_{V}$ denotes $\kappa_{Z}$ or $\kappa_{W}$, and that $\mathrm{B_{BSM}} \ge$ 0. The results for the combination of ATLAS and CMS are reported with their measured and expected uncertainties. Also shown are the results from each experiment. For the parameters with both signs allowed, the other 1$\sigma$ interval is shown on a second line. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Fit results allowing BSM loop couplings, assuming that there are no additional BSM contributions to the Higgs boson width, i.e. $\mathrm{B_{BSM}}=0$. The results for the combination of ATLAS and CMS are reported with their measured and expected uncertainties. Also shown are the results from each experiment. For the parameters with both signs allowed, the other 1$\sigma$ interval is shown on a second line. When a parameter is constrained and reaches a boundary, namely $\mathrm{B_{BSM}}$ = 0, the uncertainty is not indicated. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Fit results for the parameterisation assuming the absence of BSM particles in the loops ($\mathrm{B_{BSM}}$=0). The results with their measured and expected uncertainties are reported for the combination of ATLAS and CMS, together with the individual results from each experiment. For the parameters with both signs allowed, the other 1$\sigma$ CL interval is shown on a second line. When a parameter is constrained and reaches a boundary, namely $|\kappa_\mu|$ = 0, the uncertainty is not indicated. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Summary of fit results for a parameterisation probing the ratios of coupling modifiers for up-type versus down-type fermions. The results for the combination of ATLAS and CMS are reported together with their measured and expected uncertainties. Also shown are the results from each experiment. The parameter $\kappa_{uu}$ is positive definite since $\kappa_H$ is always assumed to be positive. For the parameter $\lambda_{du}$, for which both signs are allowed, the other 1$\sigma$ CL interval is shown on a second line. Negative values for the parameter $\lambda_{Vu}$ are excluded by more than 4$\sigma$.

Summary of fit results for a parameterisation probing the ratios of coupling modifiers for leptons versus quarks. The results for the combination of ATLAS and CMS are reported together with their measured and expected uncertainties. Also shown are the results from each experiment. The parameter $\kappa_{qq}$ is positive definite since $\kappa_H$ is always assumed to be positive. For the parameter $\lambda_{lq}$, for which there is no sensitivity to the sign, only the absolute values are shown. Negative values for the parameter $\lambda_{Vq}$ are excluded by more than 4$\sigma$.

Correlation matrix obtained from the fit combining the ATLAS and CMS data using the generic parameterisation with 23 parameters, $\sigma_i\cdot\mathrm{B}^f$ for each specific channel $i\to H\to f$. Only 20 parameters are shown because the other three, corresponding to the $H\to ZZ$ decay channel for the $WH$, $ZH$, and $ttH$ production processes, are not measured with a meaningful precision.

Correlation matrix obtained from the fit combining the ATLAS and CMS data using the generic parameterisation with nine parameters, $\sigma(gg\to H\to ZZ)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{ZZ}$.

Correlation matrix obtained from the fit combining the ATLAS and CMS data using the generic parameterisation with seven parameters, $\kappa_{gZ}=\kappa_{g}\cdot\kappa_{Z} / \kappa_{H}$ and the ratios of coupling modifiers.

Expected correlation matrix obtained from the fit combining ATLAS and CMS pre-fit Asimov data sets using the generic parameterisation with 23 parameters, $\sigma_i\cdot\mathrm{B}^f$ for each specific channel $i\to H\to f$. Only 20 parameters are shown because the other three, corresponding to the $H\to ZZ$ decay channel for the $WH$, $ZH$, and $ttH$ production processes, are not measured with a meaningful precision.

Expected correlation matrix obtained from the fit combining ATLAS and CMS pre-fit Asimov data sets using the generic parameterisation with nine parameters, $\sigma(gg\to H\to ZZ)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{ZZ}$.

Expected correlation matrix obtained from the fit combining ATLAS and CMS pre-fit Asimov data sets using the generic parameterisation with seven parameters, $\kappa_{gZ}=\kappa_{g}\cdot\kappa_{Z} / \kappa_{H}$ and the ratios of coupling modifiers.

ATLAS+CMS combined best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{g}$,$\kappa_{\gamma}$) plane.

ATLAS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{g}$,$\kappa_{\gamma}$) plane.

CMS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{g}$,$\kappa_{\gamma}$) plane.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{F}$,$\kappa_{V}$) plane.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\gamma\gamma$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow ZZ$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow WW$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\tau\tau$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow bb$ channel.

ATLAS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{F}$,$\kappa_{V}$) plane.

CMS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{F}$,$\kappa_{V}$) plane.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\gamma\gamma$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow ZZ$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow WW$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\tau\tau$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow bb$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow\gamma\gamma$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow ZZ$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow WW$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow\tau\tau$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow bb$ channel.