The EW-Zy jj differential cross-section in the Signal Region as a function of the leading lepton pT. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The EW-Zy jj differential cross-section in the Signal Region as a function of the leading photon pT. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The EW-Zy jj differential cross-section in the Signal Region as a function of the leading jet pT. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The EW-Zy jj differential cross-section in the Signal Region as a function of the Zy system pT. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The EW-Zy jj differential cross-section in the Signal Region as a function of the dijet invariant mass. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The EW-Zy jj differential cross-section in the Signal Region as a function of the dijet rapidity difference. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The EW-Zy jj differential cross-section in the Signal Region as a function of the Zy and dijet azimuthal difference. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The total Zy jj differential cross-section in the Extended SR as a function of the leading lepton pT. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The total Zy jj differential cross-section in the Extended SR as a function of the leading photon pT. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The total Zy jj differential cross-section in the Extended SR as a function of the leading jet pT. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The total Zy jj differential cross-section in the Extended SR as a function of the Z boson pT. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The total Zy jj differential cross-section in the Extended SR as a function of the Zy system pT. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The total Zy jj differential cross-section in the Extended SR as a function of the dijet invariant mass. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The total Zy jj differential cross-section in the Extended SR as a function of the dijet rapidity difference. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The total Zy jj differential cross-section in the Extended SR as a function of the Zy and dijet azimuthal difference. The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

The total Zy jj differential cross-section in the Extended SR as a function of the centrality of the system ζ(Zy). The lower panels show the ratios of the MC predictions to the data. The band around unfolded data represents the total uncertainty (including statistical uncertainty) and takes into account the correlations as obtained from the fit. The hatched area represents the uncertainty in the prediction. Events beyond the upper limit of the histogram are included in the last bin.

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.

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{lead.~lep.}}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.~lep.}}$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production. Overflow events are included in the last bin. The largest observed value is 609 GeV.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.~lep.}}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.~lep.}}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $p_{\mathrm{T}}^{\mathrm{lead.~jet}}$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production. Overflow events are included in the last bin. The largest observed value is 1485 GeV.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{lead.~jet}}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{lead.~jet}}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $H_{\mathrm{T}}$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production. Overflow events are included in the last bin. The largest observed value is 2969 GeV.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $H_{\mathrm{T}}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $H_{\mathrm{T}}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $S_{\mathrm{T}}$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production. Overflow events are included in the last bin. The largest observed value is 3296 GeV.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $S_{\mathrm{T}}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $S_{\mathrm{T}}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $m_{\mathrm{T},e\mu}$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production. Overflow events are included in the last bin. The largest observed value is 4130 GeV.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $m_{\mathrm{T},e\mu}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $m_{\mathrm{T},e\mu}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $m_{e\mu}$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production. Overflow events are included in the last bin. The largest observed value is 3519 GeV.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $m_{e\mu}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $m_{e\mu}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $p_{\mathrm{T},e\mu}$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production. Overflow events are included in the last bin. The largest observed value is 1067 GeV.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T},e\mu}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T},e\mu}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $\Delta\phi(e,\mu)$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $\Delta\phi(e,\mu)$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $\Delta\phi(e,\mu)$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $y_{e\mu}$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $y_{e\mu}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $y_{e\mu}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $\cos\theta^*$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $\cos\theta^*$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $\cos\theta^*$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $n_{\mathrm{jet}}$. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $n_{\mathrm{jet}}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $n_{\mathrm{jet}}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $m_{e\mu}$ for $p_{\mathrm{T}}^{\mathrm{lead.~jet}} > 200$ GeV. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production. Overflow events are included in the last bin. The largest observed value is 3519 GeV.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $m_{e\mu}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $m_{e\mu}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $\Delta\phi(e,\mu)$ for $p_{\mathrm{T}}^{\mathrm{lead.~jet}} > 200$ GeV. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $\Delta\phi(e,\mu)$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $\Delta\phi(e,\mu)$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $\Delta\phi(\mathrm{sub-lead.~lep.}, \mathrm{lead.~jet})$ for $p_{\mathrm{T}}^{\mathrm{lead.~lep.}} > 200$ GeV. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $\Delta\phi(\mathrm{sub-lead.~lep.}, \mathrm{lead.~jet})$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $\Delta\phi(\mathrm{sub-lead.~lep.}, \mathrm{lead.~jet})$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $\Delta R(\mathrm{sub-lead.~lep.}, \mathrm{lead.~jet})$ for $p_{\mathrm{T}}^{\mathrm{lead.~lep.}} > 200$ GeV. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $\Delta R(\mathrm{sub-lead.~lep.}, \mathrm{lead.~jet})$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $\Delta R(\mathrm{sub-lead.~lep.}, \mathrm{lead.~jet})$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.~lep.}} / p_{\mathrm{T}}^{\mathrm{lead.~lep.}}$ for $p_{\mathrm{T}}^{\mathrm{lead.~lep.}} > 200$ GeV. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.~lep.}} / p_{\mathrm{T}}^{\mathrm{lead.~lep.}}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.~lep.}} / p_{\mathrm{T}}^{\mathrm{lead.~lep.}}$

Measured fiducial cross section for $pp\rightarrow W^+W^-$+jets production for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.~lep.}} / p_{\mathrm{T}}^{\mathrm{lead.~jet}}$ for $p_{\mathrm{T}}^{\mathrm{lead.~lep.}} > 200$ GeV. The second column contains the results obtained with a fiducial particle phase space that includes a veto on $b$-jets. This alternative result is obtained from the nominal result by the application of bin-wise correction that is calculated as the ratio of the predicted differential cross-section in the nominal analysis phase space and the predicted cross-section for a phase space that includes a veto on events with $b$-jets with $p_{\mathrm{T}} > 20$ GeV. Also shown are the Standard Model predictions for $q\bar{q} \rightarrow WW$, obtained from Sherpa 2.2.2, MadGraph 2.3.3 + Pythia 8.212 using FxFx merging, and Powheg MiNLO + Pythia 8.244. These predictions are supplemented by the Sherpa 2.2.2 + OpenLoops simulation of $gg\rightarrow WW$. Finally, the prediction from MATRIX is given, which includes nNLO QCD and NLO EW corrections to $WW$+jet production. The largest observed value is 19.6.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.~lep.}} / p_{\mathrm{T}}^{\mathrm{lead.~jet}}$

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.~lep.}} / p_{\mathrm{T}}^{\mathrm{lead.~jet}}$

Cutflow for two representative signal samples used in this analysis. The excited tau mass $m_{\tau^\ast} = 2.75$ TeV and the contact interaction scale $\Lambda=10$ TeV. The LQ mass $m_\mathrm{LQ} = 1.3$ TeV. The event yields include all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Acceptance x efficiency of the $\tau^\ast$ signal SR selection

Acceptance x efficiency of the LQ signal SR selection

Expected 95% CL exclusion contours in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for boosted.

Observed 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for boosted.

Expected 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for boosted.

Observed 95% CL exclusion contours for the assumpition of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for boosted.

Expected 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for boosted.

Observed 95% CL exclusion contours for the assumption of Majorana type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for resolved.

Expected 95% CL exclusion contours for the assumption of Majorana type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for resolved.

Observed 95% CL exclusion contours for the assumption of Majorana type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for resolved.

Expected 95% CL exclusion contours for the assumption of Majorana type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for resolved.

Observed 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for resolved.

Expected 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for resolved.

Observed 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for resolved.

Expected 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for resolved.

The $m_{eJ}$ distribution in the boosted one electron channel.

The $m_{eeJ}$ distribution in the boosted two electron channel.

The $m_{\mu\mu J}$ distribution in the boosted two muon channel.

The $m_{jj}$ distribution in the resolved electron channel.

The $m_{jj}$ distribution in the resolved muon channel.

The $m_{eejj}$ distribution in the resolved electron channel.

The $m_{\mu\mu jj}$ distribution in the resolved muon channel.

The $h_{T}$ distribution in the resolved electron channel.

The $h_{T}$ distribution in the resolved muon channel.

Acceptance $\times$ efficiency for bSR1e.

Acceptance $\times$ efficiency for bSR2e.

Acceptance $\times$ efficiency for bSR2mu.

Acceptance $\times$ efficiency for rSROS2e.

Acceptance $\times$ efficiency for rSROS2mu.

Acceptance $\times$ efficiency for rSRSS2e.

Acceptance $\times$ efficiency for rSRSS2mu.

Expected and observed 95% CL upper limits for $\sigma(pp\to W_{R})\times B(W_{R}\to eeqq)$ as a function of $m(W_{R})$ for the special case of $m(W_{R}) = 2\times m(N_{R})$. The theoretical cross section corresponds to either the Dirac or Majorana neutrino interpretation by ignoring the charge of the leptons.

Expected and observed 95% CL upper limits for $\sigma(pp\to W_{R})\times B(W_{R} \to \mu\mu qq)$ as a function of $m(W_{R})$ for the special case of $m(W_{R}) = 2\times m(N_{R})$. The theoretical cross section corresponds to either the Dirac or Majorana neutrino interpretation by ignoring the charge of the leptons.

Expected and observed 95% CL upper limits for $\sigma(pp\to W_{R})\times B(W_{R} \to eeqq)$ as a function of $m(W_{R})$ for the special case of $m(N_{R})=200$ GeV. The theoretical cross section corresponds to either the Dirac or Majorana neutrino interpretation by ignoring the charge of the leptons.

Expected and observed 95% CL upper limits for $\sigma(pp\to W_{R})\times B(W_{R} \to \mu \mu qq)$ as a function of $m(W_{R})$ for the special case of $m(N_{R})=200$ GeV. The theoretical cross section corresponds to either the Dirac or Majorana neutrino interpretation by ignoring the charge of the leptons.

Expected and observed 95% CL upper limits for $\sigma(pp \to W_{R}) \times B(W_{R} \to eeqq)$ as a function of $m(W_{R})$ for the Majorana neutrino interpretations for the special case of $m(W_{R}) = 2\times m(N_{R})$.

Expected and observed 95% CL upper limits for $\sigma(pp \to W_{R}) \times B(W_{R} \to \mu\mu qq)$ as a function of $m(W_{R})$ for the Majorana neutrino interpretations for the special case of $m(W_{R}) = 2 \times m(N_{R})$.

Expected and observed 95% CL upper limits for $\sigma(pp \to W_{R}) \times B(W_{R} \to eeqq)$ as a function of $m(W_{R})$ for the Dirac neutrino interpretations for the special case of $m(W_{R}) = 2 \times m(N_{R})$.

Expected and observed 95% CL upper limits for $\sigma(pp \to W_{R}) \times B(W_{R} \to \mu \mu qq)$ as a function of $m(W_{R})$ for the Dirac neutrino interpretations for the special case of $m(W_{R}) = 2 \times m(N_{R})$.

Background-only fit to the observed numbers of events (black points) as a function of the diphoton invariant mass $m_{\gamma\gamma}$.

Expected and observed 95% CL upper limits on the signal fiducial cross section $\sigma_{\textrm{fid}}$, assuming 100% branching ratio for ALP decay into two photons, as a function of the hypothetical ALP mass $m_{\textrm{X}}$. The 1$\sigma$ and 2$\sigma$ confidence intervals are shown by the coloured bands.

Expected and observed 95% CL upper limits on the the ALP coupling constant, assuming 100% branching ratio for ALP decay into two photons, as a function of the hypothetical ALP mass $m_{\textrm{X}}$. The 1$\sigma$ and 2$\sigma$ confidence intervals are shown by the coloured bands. Contours of the ALP natural width $\Gamma$ are illustrated by the smooth blue solid lines.

Fiducial region acceptance of signal events for coupling constant f$^{-1}=0.05$ TeV$^{-1}$. EL, SD, and DD stand for the exclusive, single-dissociative, and double-dissociative processes. The acceptances are parameterised as $p_0 + p^{i}_{1} e^{p_{2}^{i} m_{\textrm{X}}} + p_{3}^{i} e^{p_{4}^{i} m_{\textrm{X}}} $ , where $i$ specifies the process type, exclusive, single-dissociative, or double-dissociative.

Correction factors $A_{\textrm{X}}/A_{0}$ and $C_{\textrm{X}}$ for the exclusive signal events with ALP coupling constant f$^{-1}$=0.05 TeV$^{-1}$, where acceptance $A_{\textrm{X}}$ and efficiency $C_{\textrm{X}}$ are defined in arXiv:2102.13405. $A_{\textrm{X}}$/$A_{0}$ and $C_{\textrm{X}}$ are parameterised as $p_0 + p^{i}_{1} e^{p_{2}^{i} m_{\textrm{X}}} + p_{3}^{i} e^{p_{4}^{i} m_{\textrm{X}}}$, where $i$ specifies the process type.

Simulated ALP signal cutflow for $m_{\rm{X}}$=300 GeV and $m_{\rm{X}}$=1200 GeV. The yields are normalized to a luminosity of 14.6 fb$^{-1}$. Selections are applied for each side after the selection on acoplanarity. The union of the sets of selected events is taken finally.

Di-muon mass distribution for the $\mu\mu\mu$ signal region for data and expected background. The expected signal distribution for $m_a = 35$ GeV is shown assuming $\sigma(t\bar{t}a)\times \text{Br}(a\rightarrow\mu\mu) = 4$ fb. Rare background processes include $ZZ+$jets, $WWZ$, $WZZ$, $ZZZ$, and $t\bar{t}t\bar{t}$.

Fit of the Double Crystal Ball width as a function of the $a$-boson mass in the $e\mu\mu$ channel.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $t\bar{t}a$.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $H^{+}\rightarrow aW$, considering a top pair cross section of $\sigma(pp\rightarrow t\bar{t})= 833$ pb for the charged Higgs boson mass, $m_{H^{+}} = 120$ GeV.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $H^{+}\rightarrow aW$, considering a top pair cross section of $\sigma(pp\rightarrow t\bar{t})= 833$ pb for the charged Higgs boson mass, $m_{H^{+}} = 140$ GeV.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $H^{+}\rightarrow aW$, considering a top pair cross section of $\sigma(pp\rightarrow t\bar{t})= 833$ pb for the charged Higgs boson mass, $m_{H^{+}} = 160$ GeV.

Comparison of the expected background with data in a region dominated by non-prompt background. The region is defined exactly like the $e\mu\mu$ signal region but requiring the two selected muons to have the same electrical charge.

Scan of the observed p-value as a function of $m_a$ for the background-only hypothesis.

Scan of the observed p-value as a function of $m_a$ for the background-only hypothesis for $e\mu\mu$ channel.

Scan of the observed p-value as a function of $m_a$ for the background-only hypothesis for $\mu\mu\mu$ channel.

Acceptance times efficiency for the $t\bar{t}a$, $a\rightarrow \mu\mu$ signal as a function of the $m_a$ hypothesis.

Acceptance times efficiency for the $t \rightarrow bH^+$, $H^+ \rightarrow W^+ a$, $a\rightarrow \mu\mu$ signal as a function of the $m_a$ and $m_{H^+}$ hypotheses, for the $e\mu\mu$ signal region. The upper left region of the plot corresponds to signals for which on-shell decays of $H^+ \rightarrow Wa$ are not possible. The CP conjugate process is also included.

Acceptance times efficiency for the $t \rightarrow bH^+ $, $H^+ \rightarrow W^+ a$, $a\rightarrow \mu\mu$ signal as a function of the $m_a$ and $m_{H^+}$ hypotheses, for the $\mu\mu\mu$ signal region. The upper left region of the plot corresponds to signals for which on-shell decays of $H^+ \rightarrow Wa$ are not possible. The CP conjugate process is also included.

Expected lower limits on $\tan\beta$ for the $t\rightarrow H^{+}b, H^{+}\rightarrow aW, a\rightarrow\mu\mu$ signal in the context of the 2HDM as a function of $m_{H^{+}}$ and $m_a$.

Observed lower limits on $\tan\beta$ for the $t\rightarrow H^{+}b, H^{+}\rightarrow aW, a\rightarrow\mu\mu$ signal in the context of the 2HDM as a function of $m_{H^{+}}$ and $m_a$.

Cutflow for two signal mass points for $t\bar{t}a$, $m_{a} = 20$ GeV and $m_{a} = 60$ GeV, and one signal mass point for $H^+ \rightarrow Wa$ for the electron channel $e\mu\mu$. For the signal yields, a cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Cutflow for the dominant backgrounds estimated from simulation ($t\bar{t}Z$, $WZ$, $t\bar{t}H$) and data, for the electron channel $e\mu\mu$. The yields are presented before the profile likelihood fit.

Cutflow for the signal mass point $t\bar{t}a$, $m_{a} = 20$ GeV, as well as the dominant backgrounds estimated from simulation ($t\bar{t}Z$, $WZ$, $t\bar{t}H$) and data, for the corresponding signal mass hypothesis, for the muon channel $\mu\mu\mu$. For the signal yields, a cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Cutflow for the signal mass point $t\bar{t}a$, $m_{a} = 60$ GeV, as well as the dominant backgrounds estimated from simulation ($t\bar{t}Z$, $WZ$, $t\bar{t}H$) and data, for the corresponding signal mass hypothesis, for the muon channel $\mu\mu\mu$. For the signal yields, a cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Cutflow for $H^+ \rightarrow W a$ with $m_{H^{+}} = 120$ GeV, $m_{a} = 30$ GeV for the muon channel $\mu\mu\mu$. A cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Total yields in the $e\mu\mu$ and $\mu\mu\mu$ signal regions for the expected background and data after the profile likelihood fit for a $m_a = 30$ GeV signal mass hypothesis.

The $\Delta p_{\mathrm{T}}$ of photon-tagged jets as a function of jet $p_{\mathrm{T}}$.

The double $R_{\mathrm{AA}}$ ratio of $R_{\mathrm{AA}}^{\gamma\text{-jet}}/R_{\mathrm{AA}}^{\text{inclusive jet}}$ as a function of jet $p_{\mathrm{T}}$.

The differential cross-section of photon-tagged jets as a function of jet $p_{\mathrm{T}}$ in pp collisions.

The yields of photon-tagged jets as a function of jet $p_{\mathrm{T}}$ in Pb+Pb collisions for different centrality intervals.

The nuclear modification factor of photon-tagged jets as a function of jet $p_{\mathrm{T}}$ for different centrality intervals.

The $\Delta p_{\mathrm{T}}$ of photon-tagged jets as a function of jet $p_{\mathrm{T}}$.

The $R_{\mathrm{CP}}$ of photon-tagged jets as a function of jet $p_{\mathrm{T}}$.