The 2D distribution of the $\Delta\Phi$ (cluster, $\mu_{trigger}$) vs the number of hits in clusters, in events from the out-of-time region.

Exclusion limits at 95% CL on the $B \rightarrow K \Phi$ as a function of the decay distance of the long-lived particle.

Exclusion limits at 95% CL on the $B \rightarrow K \Phi$ as a function of the decay distance of the long-lived particle.

Exclusion limits at 95% CL on the $B \rightarrow K \Phi$ as a function of the decay distance of the long-lived particle.

Exclusion limits at 95% CL on the $B \rightarrow K \Phi$ as a function of the decay distance of the long-lived particle.

Exclusion limits at 95% CL on the $B \rightarrow K \Phi$ as a function of the decay distance of the long-lived particle.

The estimated BR limits for a range of generator parameters: $0.3 < m_{LLP}$ < 3.0 GeV and 0 < $c\tau$ < 10000 mm.

The contour of the 2D distribution which includes the generation parameter values corresponding to the BR of 0.01.

The contour of the 2D distribution which includes the generation parameter values corresponding to the BR of 0.001.

Signal cut-flow efficiencies for the hadronic shower decay mode

Signal cut-flow efficiencies for the electromagnetic shower decay mode

Predictions of a leading order (LO) QCD simulation, normalized to an integrated luminosity of 138 fb$^{-1}$. The number of events are examined within bins of the four-jet mass and the ratio $\alpha$, which is the average dijet mass divided by the four-jet mass.

The 68% probability contour in the $m_{\mathrm{4j}}$ vs. $\overline{m}_{\mathrm{2j}}$ plane from a signal simulation of a diquark with a width of 1.5% and a mass of 8.4 TeV, decaying to a pair of vector-like quarks, each with a mass of 2.1 TeV.

The 68% probability contour in the $m_{\mathrm{4j}}$ vs. $\overline{m}_{\mathrm{2j}}$ plane from a signal simulation of a diquark with a width of 5% and a mass of 8.4 TeV, decaying to a pair of vector-like quarks, each with a mass of 2.1 TeV.

The 68% probability contour in the $m_{\mathrm{4j}}$ vs. $\overline{m}_{\mathrm{2j}}$ plane from a signal simulation of a diquark with a width of 10% and a mass of 8.4 TeV, decaying to a pair of vector-like quarks, each with a mass of 2.1 TeV.

The 68% probability contour in the $m_{\mathrm{4j}}$ vs. $\alpha$ plane from a signal simulation of a diquark with a width of 1.5% and a mass of 8.4 TeV, decaying to a pair of vector-like quarks, each with a mass of 2.1 TeV.

The 68% probability contour in the $m_{\mathrm{4j}}$ vs. $\alpha$ plane from a signal simulation of a diquark with a width of 5% and a mass of 8.4 TeV, decaying to a pair of vector-like quarks, each with a mass of 2.1 TeV.

The 68% probability contour in the $m_{\mathrm{4j}}$ vs. $\alpha$ plane from a signal simulation of a diquark with a width of 10% and a mass of 8.4 TeV, decaying to a pair of vector-like quarks, each with a mass of 2.1 TeV.

Signal differential distributions as a function of four-jet mass for $\alpha_{\mathrm{true}}$ = 0.25, diquark masses of 2, 5, 8.6 TeV and various widths, for all $\alpha$ bins inclusively. The integral of each distribution has been normalized to unity.

The product of acceptance, $A$, and efficiency, $\varepsilon$, of a resonant signal with $\alpha_{\mathrm{true}}$ = 0.25 vs. the diquark mass for various diquark widths, and for all $\alpha$ bins inclusively. The acceptance is defined as the fraction of generated events passing the kinematic selection criteria, while the efficiency is the fraction of signal events satisfying $m_{\mathrm{4j}} > 1.6$ TeV. We also show the signal acceptance alone in the curves where $\varepsilon = 1$.

The four-jet mass distribution in data for 0.22 < $\alpha$ < 0.24, fitted with three background-only functions (Dijet-3p, PowExp-3p and ModDijet-3p), each with three free parameters. Examples of predicted diquark resonances with $\alpha_{\mathrm{true}}$ = 0.25, $M_{\mathrm{S}}$ = 8.6 TeV, and $\Gamma/M_{\mathrm{S}}$ = 1.5%, 10% are also shown, with cross sections equal to the observed upper limits at 95% confidence level.

The four-jet mass distribution in data for 0.24 < $\alpha$ < 0.26, fitted with three background-only functions (Dijet-3p, PowExp-3p and ModDijet-3p), each with three free parameters. Examples of predicted diquark resonances with $\alpha_{\mathrm{true}}$ = 0.25, $M_{\mathrm{S}}$ = 8.6 TeV, and $\Gamma/M_{\mathrm{S}}$ = 1.5%, 10% are also shown, with cross sections equal to the observed upper limits at 95% confidence level.

The four-jet mass distribution in data for 0.26 < $\alpha$ < 0.28, fitted with three background-only functions (Dijet-3p, PowExp-3p and ModDijet-3p), each with three free parameters. Examples of predicted diquark resonances with $\alpha_{\mathrm{true}}$ = 0.25, $M_{\mathrm{S}}$ = 8.6 TeV and $\alpha_{\mathrm{true}}$ = 0.29, $M_{\mathrm{S}}$ = 3.6 TeV, each shown for widths of $\Gamma/M_{\mathrm{S}}$ = 1.5% and 10% are also included, with cross sections equal to the observed upper limits at 95% confidence level.

The four-jet mass distribution in data for 0.28 < $\alpha$ < 0.30, fitted with three background-only functions (Dijet-3p, PowExp-3p and ModDijet-3p), each with three free parameters. Examples of predicted diquark resonances with $\alpha_{\mathrm{true}}$ = 0.25, $M_{\mathrm{S}}$ = 8.6 TeV and $\alpha_{\mathrm{true}}$ = 0.29, $M_{\mathrm{S}}$ = 3.6 TeV, each shown for widths of $\Gamma/M_{\mathrm{S}}$ = 1.5% and 10% are also included, with cross sections equal to the observed upper limits at 95% confidence level.

The four-jet mass distribution in data for 0.30 < $\alpha$ < 0.32, fitted with three background-only functions (Dijet-3p, PowExp-3p and ModDijet-3p), each with three free parameters. Examples of predicted diquark resonances with $\alpha_{\mathrm{true}}$ = 0.25, $M_{\mathrm{S}}$ = 8.6 TeV and $\alpha_{\mathrm{true}}$ = 0.29, $M_{\mathrm{S}}$ = 3.6 TeV, each shown for widths of $\Gamma/M_{\mathrm{S}}$ = 1.5% and 10% are also included, with cross sections equal to the observed upper limits at 95% confidence level.

The four-jet mass distribution in data for 0.32 < $\alpha$ < 0.34, fitted with three background-only functions (Dijet-3p, PowExp-3p and ModDijet-3p), each with three free parameters. Examples of predicted diquark resonances with $\alpha_{\mathrm{true}}$ = 0.25, $M_{\mathrm{S}}$ = 8.6 TeV, and $\Gamma/M_{\mathrm{S}}$ = 1.5%, 10% are also shown, with cross sections equal to the observed upper limits at 95% confidence level.

The inclusive four-jet mass distribution in data for $\alpha$ > 0.10, fitted with three background-only functions (Dijet-5p, PowExp-5p and ModDijet-5p), each with five free parameters. Examples of predicted diquark resonances with $\alpha_{\mathrm{true}}$ = 0.25, $M_{\mathrm{S}}$ = 8.6 TeV and $\alpha_{\mathrm{true}}$ = 0.29, $M_{\mathrm{S}}$ = 3.6 TeV, each shown for widths of $\Gamma/M_{\mathrm{S}}$ = 1.5% and 10% are also included, with cross sections equal to the observed upper limits at 95% confidence level.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.25$ and width of the initial resonance Y equal to 1.5%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 1.5%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.25$ and width of the initial resonance Y equal to 5%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 5%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.25$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.11$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.13$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.15$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.17$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.19$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.21$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.23$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.27$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.29$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.31$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.33$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.42$ and width of the initial resonance Y equal to 10%. The corresponding expected limits and their variations at the 1 and 2 standard deviation (s.d.) levels are also included. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate a mediator width equal to 10%.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.11$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.13$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.15$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.17$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.19$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.21$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.23$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.27$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.29$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.31$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.33$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for resonant production of paired dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}} / M_{\mathrm{Y}} = 0.42$ and widths of the initial resonance Y equal to 1.5, 5, and 10%. Limits are compared to predictions for scalar $\mathrm{S}_{\mathrm{uu}}$ and $\mathrm{S}_{\mathrm{dd}}$ diquarks with couplings to pairs of up and down quarks, $y_{\mathrm{uu}}$ and $y_{\mathrm{dd}}$, and to pairs of vector-like quarks, $y_{\chi}$ and $y_{\omega}$, set appropriately in order to generate the corresponding widths.

Observed local $p$-value for a four-jet resonance, Y, decaying to a pair of dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}}/M_{\mathrm{Y}}$ = 0.25, and various widths of Y superimposed.

Observed local $p$-value for a four-jet resonance, Y, decaying to a pair of dijet resonances, X, with $\alpha_{\mathrm{true}} = M_{\mathrm{X}}/M_{\mathrm{Y}}$ = 0.29, and various widths of Y superimposed.

The table presents the cumulative cutflow for a signal with $M_{\mathrm{S}} = 2.0$ TeV, $M_{\chi} = 0.5$ TeV and different width hypotheses ($\Gamma/M_{\mathrm{S}} =$ 1.5, 5, and 10%). Each row corresponds to the number of signal events that survive all cuts up to and including the one listed. The percentage in parentheses shows the efficiency of the current cut alone, defined as the ratio of the number of events that survive all cuts up to and including this one to the number of events that survived all previous cuts.

The table presents the cumulative cutflow for a signal with $M_{\mathrm{S}} = 5.0$ TeV, $M_{\chi} = 1.25$ TeV and different width hypotheses ($\Gamma/M_{\mathrm{S}} =$ 1.5, 5, and 10%). Each row corresponds to the number of signal events that survive all cuts up to and including the one listed. The percentage in parentheses shows the efficiency of the current cut alone, defined as the ratio of the number of events that survive all cuts up to and including this one to the number of events that survived all previous cuts.

The table presents the cumulative cutflow for a signal with $M_{\mathrm{S}} = 8.6$ TeV, $M_{\chi} = 2.15$ TeV and different width hypotheses ($\Gamma/M_{\mathrm{S}} =$ 1.5, 5, and 10%). Each row corresponds to the number of signal events that survive all cuts up to and including the one listed. The percentage in parentheses shows the efficiency of the current cut alone, defined as the ratio of the number of events that survive all cuts up to and including this one to the number of events that survived all previous cuts.

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|y_{e u}|$, versus $S_{e u}$ mass for scalar LQs coupled to electrons.

Observed and Expected UL exclusions on the $BR(H\to S)$ of generic signals with $m_{A'} = 1.0\;GeV$ and $BR(A' \rightarrow \pi\pi) = 1.0$.

The observed data in the dielectron channel and the fitted signal-plus-background templates, shown for the $S_{e d}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Observed and Expected UL exclusions on the number of 'signal' yields in sideband A from the model agnostic fit, assuming no 'SUEP-like' signal contamination in the CRs.

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|y_{e d}|$, versus $S_{e d}$ mass for scalar LQs coupled to electrons.

Postfit values for SR in the B-Only CR+SR and CR-Only fits with the prefit signal overlayed for reference

The observed data in the dimuon channel and the fitted signal-plus-background templates, shown for the $S_{\mu u}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Postfit values for PR in the B-Only CR+SR and CR-Only fits with the prefit signal overlayed for reference

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|y_{\mu u}|$, versus $S_{\mu u}$ mass for scalar LQs coupled to muons.

Postfit values for CRDY in the B-Only CR+SR and CR-Only fits with the prefit signal overlayed for reference

The observed data in the dimuon channel and the fitted signal-plus-background templates, shown for the $S_{\mu d}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Postfit values for CRTT in the B-Only CR+SR and CR-Only fits with the prefit signal overlayed for reference

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|y_{\mu d}|$, versus $S_{\mu d}$ mass for scalar LQs coupled to muons.

Per-selection efficiencies on Generic - $m_{A'} = 1.0\;GeV$ $BR(A' \rightarrow \pi\pi) = 1.0$ Hadronic - $m_{A'} = 0.7\;GeV$ and $BR(A' \rightarrow ee) = BR(A' \rightarrow \mu\mu) = 0.15$ $BR(A' \rightarrow \pi\pi) = 0.7$ Leptonic - $m_{A'} = 0.5\;GeV$ and $BR(A' \rightarrow ee) = BR(A' \rightarrow \mu\mu) = 0.2$ and $BR(A' \rightarrow \pi\pi) = 0.6$ samples. Selections: 2 OSSF leptons, one SUEP cluster, $60\leq Z_{m}\leq 120$ GeV, $p_{T}^Z\geq 25$ GeV, veto events with b-tagged jets, $p_{T}^{SUEP} > 60\;GeV$, and $dR^{Ak4}_{Ak15} < 1.5$. Efficiencies are the ratio of the histogram size at the selection step to the previous step (except twoLepton vs. initial). The final column is the total cumulative efficiency (dR vs. initial). Note that these selection numbers were all calculated using samples generated with 600,000 events for each signal point considered.

The observed data in the dielectron channel and the fitted signal-plus-background templates, shown for the $V_{e u}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Interpolated observed and expected UL exclusions on $BR(H\to S) = 1$ for generic signals with $m_{A'} = 1.0\;GeV$ and $BR(A' \rightarrow \pi\pi) = 1.0$. The discreate heatmap points can be found in the UL Generic (Observed & Expected) figure in this sandbox

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|g_{e u}|$, versus $V_{e u}$ mass for vector LQs coupled to electrons.

The observed data in the dielectron channel and the fitted signal-plus-background templates, shown for the $V_{e d}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|g_{e d}|$, versus $V_{e d}$ mass for vector LQs coupled to electrons.

The observed data in the dimuon channel and the fitted signal-plus-background templates, shown for the $V_{\mu u}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|g_{\mu u}|$, versus $V_{\mu u}$ mass for vector LQs coupled to muons.

The observed data in the dimuon channel and the fitted signal-plus-background templates, shown for the $V_{\mu d}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|g_{\mu d}|$, versus $V_{\mu d}$ mass for vector LQs coupled to muons.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the SR (t-pass) and SR (t-fail) after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. A representative mono-top signal (vector coupling scenario) with a mediator mass of 1 TeV, a DM candidate mass of 150 GeV, and a cross section of 1 pb is overlaid as an orange line. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the SR (t-pass) and SR (t-fail) after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. A representative mono-top signal (vector coupling scenario) with a mediator mass of 1 TeV, a DM candidate mass of 150 GeV, and a cross section of 1 pb is overlaid as an orange line. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the W(e$\nu$) (t-pass) and W(e$\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the W(e$\nu$) (t-pass) and W(e$\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the W($\mu\nu$) (t-pass) and W($\mu\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the W($\mu\nu$) (t-pass) and W($\mu\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the Z(ee) (t-pass) and Z(ee) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the Z(ee) (t-pass) and Z(ee) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the Z($\mu\mu$) (t-pass) and Z($\mu\mu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the Z($\mu\mu$) (t-pass) and Z($\mu\mu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the t$\bar{t}$ (e$\nu$) (t-pass) and t$\bar{t}$ (e$\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the t$\bar{t}$ (e$\nu$) (t-pass) and t$\bar{t}$ (e$\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the tt ($\mu\nu$) (t-pass) and tt ($\mu\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the tt ($\mu\nu$) (t-pass) and tt ($\mu\nu$) (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the $\gamma$ (t-pass) and $\gamma$ (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Postfit distributions of the magnitude of the hadronic recoil $R_{\mathrm{T}}$ in the $\gamma$ (t-pass) and $\gamma$ (t-fail) CRs after a fit of the background model to the data. The last bin of each distribution also contains events with $R_{\mathrm{T}}$ > 1000 GeV. The background processes are stacked together. The gray band represents the statistical and postfit systematic uncertainties in the predicted background yields after the fit.

Upper limits at 95% CL on $\sigma\mathcal{B}$ of mono-top production presented in the two- dimensional plane spanned by the mediator and DM candidate masses for a mediator mass between 200 and 2250 GeV and a DM candidate mass between 1 and 1125 GeV only consider- ing on-shell decays of the mediator to the DM candidates. The mediator has vector couplings to quarks and DM candidates in the upper plot and axial-vector couplings in the lower plot. The median expected exclusion range is indicated by a black solid line, demonstrating the search sensitivity of the analysis. The 68% probability interval of the expected exclusion is shown in black dashed lines. Contours of theory predictions for constant values of $\sigma\mathcal{B}$ are shown in gray dashed lines. The observed exclusion contour of mediator and DM candidate masses is rep- resented by the red solid line. The exclusion contour obtained from measurements of the DM relic density $\Omega_{\mathrm{nbm}}h^2$ by the Planck Collaboration is shown in the gray solid line.

Upper limits at 95% CL on $\sigma\mathcal{B}$ of mono-top production presented in the two- dimensional plane spanned by the mediator and DM candidate masses for a mediator mass between 200 and 2250 GeV and a DM candidate mass between 1 and 1125 GeV only consider- ing on-shell decays of the mediator to the DM candidates. The mediator has vector couplings to quarks and DM candidates in the upper plot and axial-vector couplings in the lower plot. The median expected exclusion range is indicated by a black solid line, demonstrating the search sensitivity of the analysis. The 68% probability interval of the expected exclusion is shown in black dashed lines. Contours of theory predictions for constant values of $\sigma\mathcal{B}$ are shown in gray dashed lines. The observed exclusion contour of mediator and DM candidate masses is rep- resented by the red solid line. The exclusion contour obtained from measurements of the DM relic density $\Omega_{\mathrm{nbm}}h^2$ by the Planck Collaboration is shown in the gray solid line.

Distributions of the number of hits in the cluster (Nhits) for the DT category in the signal region (SR). The last histogram bin contains all overflow events.

Distributions of the number of hits in the cluster (Nhits) for the DT category in the out-of-time (OOT) region. The last histogram bin contains all overflow events.

Distributions of the number of hits in the cluster (Nhits) for the CSC category in the signal region (SR). The last histogram bin contains all overflow events.

Distributions of the number of hits in the cluster (Nhits) for the CSC category in the out-of-time (OOT) region. The last histogram bin contains all overflow events.

Distributions of the number of hits in the cluster (Nhits) for the DT category in the out-of-time (OOT) region. The last histogram bin contains all overflow events.

The 95% CL observed and expected upper limits on the VLL production cross section as a function of the VLL mass for the pseudoscalar mean proper decay length c$\tau$ = 0.025 m. The pseudoscalar mass is 2 GeV.

Distributions of the number of hits in the cluster (Nhits) for the CSC category in the out-of-time (OOT) region. The last histogram bin contains all overflow events.

The 95% CL observed and expected upper limits on the VLL production cross section as a function of the pseudoscalar mean proper decay length c$\tau$ for VLL Mass of 700 GeV. The pseudoscalar mass is 2 GeV.

The 95% CL observed and expected upper limits on the VLL production cross section as a function of the VLL mass for the pseudoscalar mean proper decay length c$\tau$ = 0.025 m. The pseudoscalar mass is 2 GeV.

The 95% CL observed upper limits on the VLL production cross section as a function of the VLL mass and the pseudoscalar mean proper decay length c$\tau$. The pseudoscalar mass is 2 GeV.

The 95% CL observed and expected upper limits on the VLL production cross section as a function of the pseudoscalar mean proper decay length c$\tau$ for VLL Mass of 700 GeV. The pseudoscalar mass is 2 GeV.

The 95% CL observed exclusion contour on the VLL production cross section as a function of the VLL mass and the pseudoscalar mean proper decay length c$\tau$. The pseudoscalar mass is 2 GeV.

The 95% CL observed upper limits on the VLL production cross section as a function of the VLL mass and the pseudoscalar mean proper decay length c$\tau$. The pseudoscalar mass is 2 GeV.

Selection efficiencies in the CSC category signal region (SR) for different signal hypothesis. The pseudoscalar mass is 2 GeV.

The 95% CL observed exclusion contour on the VLL production cross section as a function of the VLL mass and the pseudoscalar mean proper decay length c$\tau$. The pseudoscalar mass is 2 GeV.

Selection efficiencies in the DT category signal region (SR) for different signal hypothesis. The pseudoscalar mass is 2 GeV.

Selection efficiencies in the CSC category signal region (SR) for different signal hypothesis. The pseudoscalar mass is 2 GeV.

Selection efficiencies in the DT category signal region (SR) for different signal hypothesis. The pseudoscalar mass is 2 GeV.

The cluster reconstruction efficiency, including both DT and CSC clusters, as a function of the simulated r and |z| decay positions of the pseudoscalar into photons in events with MET > 200 GeV, for a VLL mass of 300 GeV and a pseudoscalar mass of 2 GeV, and a range of ctau values uniformly distributed between 0.01 and 0.1 m.

The cluster reconstruction efficiency, including both DT and CSC clusters, as a function of the simulated r and |z| decay positions of the pseudoscalar into photons in events with MET > 200 GeV, for a VLL mass of 500 GeV and a pseudoscalar mass of 2 GeV, and a range of ctau values uniformly distributed between 0.01 and 0.1 m.

Expected and observed upper limits for the b-quark-associated Higgs boson production cross section times branching fraction of the decay into a b quark pair at 95% CL as functions of $m_\phi$ for the 2018 FH category. The vertical dashed lines indicate the boundaries of usage of the different fit ranges, as reflected in the rightmost column of Table 2.

Expected and observed upper limits for the b-quark-associated Higgs boson production cross section times branching fraction of the decay into a b quark pair at 95% CL as functions of $m_\phi$, corresponding to the Run 2 combination. The vertical dashed line separates the mass range where only the 2017 SL category contributes on its left, from the region where also the 2017 FH and 2018 FH categories contribute on its right. The expected limits from the 2017 SL, 2017 FH, and 2018 FH datasets as well as from the previously published result based on the 2016 dataset are also shown as colored lines.

Interpretation in the $M_h^{125}$ scenario of the MSSM: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of the mass of the CP-odd boson, $m_{A}$. The higgsino mass parameter has been set to $\mu = +1$ TeV. The hashed area indicates the parameter region in which the mass of the lightest MSSM Higgs boson does not coincide with 125 GeV within a margin of 3 GeV.

Interpretation in the $M_h^{125}$ scenario of the MSSM: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of the mass of the CP-odd boson, $m_{A}$. The higgsino mass parameter has been set to $\mu = -1$ TeV. The hashed area indicates the parameter region in which the mass of the lightest MSSM Higgs boson does not coincide with 125 GeV within a margin of 3 GeV.

Interpretation in the $M_h^{125}$ scenario of the MSSM: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of the mass of the CP-odd boson, $m_{A}$. The higgsino mass parameter has been set to $\mu = -2$ TeV. The hashed area indicates the parameter region in which the mass of the lightest MSSM Higgs boson does not coincide with 125 GeV within a margin of 3 GeV.

Interpretation in the $M_h^{125}$ scenario of the MSSM: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of the mass of the CP-odd boson, $m_{A}$. The higgsino mass parameter has been set to $\mu = -3$ TeV. The hashed area indicates the parameter region in which the mass of the lightest MSSM Higgs boson does not coincide with 125 GeV within a margin of 3 GeV.

Interpretation in the $m_h^{mod+}$ scenario of the MSSM: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of the mass of the CP-odd boson, $m_{A}$. The hashed area indicates the parameter region in which the mass of the lightest MSSM Higgs boson does not coincide with 125 GeV within a margin of 3 GeV.

Interpretation in the hMSSM scenario of the MSSM: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of the mass of the CP-odd boson, $m_{A}$.

Interpretation in 2HDM scenarios: observed and expected upper limits at 95% CL on the parameter tan(beta) as a function of $m_{A}$ for cos(beta-alpha) = 0.1, for the 2HDM $\it Type-II$ scenario.

Interpretation in 2HDM scenarios: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of cos(beta-alpha) for masses of $m_{A} = m_{H} = 300$ GeV, for the 2HDM $\it Type-II$ scenario.

Interpretation in 2HDM scenarios: observed and expected upper limits at 95% CL on the parameter tan(beta) as a function of $m_{A}$ for cos(beta-alpha) = 0.1, for the 2HDM $\it Flipped$ scenario.

Interpretation in 2HDM scenarios: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of cos(beta-alpha) for masses of $m_{A} = m_{H} = 300$ GeV, for the 2HDM $\it Flipped$ scenario.

Interpretation in the 2HDM flipped scenario: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of cos(beta-alpha) for mass of $m_{A} = m_{H} = 140$ GeV.

Interpretation in the 2HDM flipped scenario: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of cos(beta-alpha) for mass of $m_{A} = m_{H} = 600$ GeV.

Interpretation in the 2HDM flipped scenario: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of cos(beta-alpha) for mass of $m_{A} = m_{H} = 900$ GeV.

Interpretation in the 2HDM flipped scenario: observed and expected upper limits at 95% CL on the parameter tan(beta) as functions of cos(beta-alpha) for mass of $m_{A} = m_{H} = 1200$ GeV.

Simulated signal shapes for three representative values of the Higgs boson mass $m_\phi$ in the 2017 SL channel.

Simulated signal shapes for three representative values of the Higgs boson mass $m_\phi$ in the 2017 FH channel.

Simulated signal shapes for three representative values of the Higgs boson mass $m_\phi$ in the 2018 FH channel.

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the SL category

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the SL category

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the SL category

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2017 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2018 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2018 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2018 analysis in the FH category

Background-only fits to the M12 distribution in each fit range of the 2018 analysis in the FH category

The contour for the excluded mass--lifetime region for stau pair production obtained with the di-track search. All masses and lifetimes shown that are below the curve and above 200 GeV are excluded by the observed data (while the expected exclusion is between the upper curve down to 210 GeV for lifetimes above 3000 ns). The sensitivity extends indefinitely to longer lifetimes.

The contour for the excluded mass--lifetime region for stau pair production obtained with the di-track search. All masses and lifetimes shown that are below the curve and above 200 GeV are excluded by the observed data (while the expected exclusion is between the upper curve down to 210 GeV for lifetimes above 3000 ns). The sensitivity extends indefinitely to longer lifetimes.

The contour for the excluded mass--lifetime region for stau pair production obtained with the di-track search. All masses and lifetimes shown that are below the curve and above 200 GeV are excluded by the observed data (while the expected exclusion is between the upper curve down to 210 GeV for lifetimes above 3000 ns). The sensitivity extends indefinitely to longer lifetimes.

The distribution of data and predicted background in the di-track Exclusion-SR. The observed data events are indicated as magenta circles if they are inside the mass-compatibility angle (shown as grey lines) and as orange diamonds if they are outside, while the blue area is the mass distribution of the expected background. The last bins include overflow events

The distribution of data and predicted background in the di-track Exclusion-SR. The observed data events are indicated as magenta circles if they are inside the mass-compatibility angle (shown as grey lines) and as orange diamonds if they are outside, while the blue area is the mass distribution of the expected background. The last bins include overflow events

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with 3 ns lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with 3 ns lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with 3 ns lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with 3 ns lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with 3 ns lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with 3 ns lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with 10 ns lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with 10 ns lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with 10 ns lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with 10 ns lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with 10 ns lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with 10 ns lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with 30 ns lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with 30 ns lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with 30 ns lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with 30 ns lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with 30 ns lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with 30 ns lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with a detector-stable lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with a detector-stable lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with a detector-stable lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with a detector-stable lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with a detector-stable lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus mass for staus with a detector-stable lifetime.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 200 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 200 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 200 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 200 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 200 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 200 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 300 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 300 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 300 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 300 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 300 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 300 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 400 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 400 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 400 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 400 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 400 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 400 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 500 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 500 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 500 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 500 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 500 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 500 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 600 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 600 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 600 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 600 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 600 GeV staus.

The 95\% CL exclusion limits from the di-track search on cross-section versus lifetime for 600 GeV staus.

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Example description

Lower limits on tab$eta$ as a function of m(A) in the mhEFT(125) scenario. Values below the black solid line are excluded at $95%$ CL

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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