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.

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{loose}}^{2\ell3j}$ signal region, which uses the invariant mass of the c-tagged jet and the geometrically closest b-tagged jet, $m_{\mathit{cb}}^{\mathrm{min}\Delta R}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{loose}}^{1\ell6ji}$ signal region, which uses the jet multiplicity, $N_{\mathrm{jets}}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{tight}}^{1\ell6ji}$ signal region, which uses the jet multiplicity, $N_{\mathrm{jets}}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{loose}}^{2\ell4ji}$ signal region, which uses the jet multiplicity, $N_{\mathrm{jets}}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{tight}}^{2\ell4ji}$ signal region, which uses the jet multiplicity, $N_{\mathrm{jets}}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Jet multiplicity and tagging criteria for the SRs and CRs in the single-lepton and dilepton channels. All single-lepton regions are defined in the 5-jet-exclusive and 6-jet-inclusive jet selections, all dilepton regions in the 3-jet-exclusive and 4-jet-inclusive jet selections. By construction, $\mathrm{SR}^{2\ell}_{\mathrm{tight}}$ only exists for the 4-jet-inclusive selection. Identical requirements on inclusive and exclusive working points, e.g., one c@22% and one c@11% in the $\mathrm{CR}_1^{1\ell}$, indicate a veto of additional tags at the inclusive working point.

Breakdown of the fiducial $t\bar{t}+{\geq}2c$ and $t\bar{t}+1c$ cross-section uncertainties. Systematic uncertainties are grouped into signal and background modeling, instrumental, and MC statistics categories. The table lists the fractional uncertainty in the measured cross-sections in percent and is estimated from the covariance matrix of the profile likelihood fit, following EPJC 84 (2024) 593. Individual groups of uncertainties can be added in quadrature to obtain the total uncertainty. JES and JER denote the jet energy scale and resolution uncertainties, respectively. The fractional uncertainties in the fitted $t\bar{t}+{\geq}1c$ and $t\bar{t}+\text{light}$ normalization factors are included in the data statistics category. For presentation purposes, all uncertainties are symmetrized.

Measured fiducial cross-section values in comparison with various NLO+PS predictions. The uncertainties in the predictions include independent and simultaneous variations of the $\mu_R$ and $\mu_F$ scales and uncertainties in the choice of the PDF set, but no uncertainties in the parton shower or hadronization. The uncertainties in the measured values include all statistical and systematic uncertainties. For presentation purposes, the quoted measurement and prediction uncertainties are symmetrized.

Measured and predicted values for the $t\bar{t}+{\geq}1b$, $t\bar{t}+{\geq}2c$, and $t\bar{t}+1c$ cross-sections and for the cross-section ratios of these processes to total $t\bar{t}+\text{jets}$ production. The quoted uncertainties in the measurements include statistical and systematic uncertainties. The predictions are taken from the nominal setup using inclusive $t\bar{t}$ Powheg+Pythia8 predictions for the $t\bar{t}+{\geq}2c$, $t\bar{t}+1c$, and $t\bar{t}+\text{light}$ processes, and $t\bar{t}+b\bar{b}$ Powheg+Pythia8 predictions for the $t\bar{t}+{\geq}1b$ process. The first use the 5FS, the latter the 4FS, where the $b$-quarks are treated as massive particles and are included in the matrix element calculation. The uncertainties in the predictions include independent and simultaneous variations of the $\mu_R$ and $\mu_F$ scales and uncertainties in the choice of the PDF set.

Nuisance parameters with a large impact on the fiducial measurement: ranking of the ten most impactful nuisance parameters for the $t\bar{t}+{\geq}2c$ cross-section measurement. The colored boxes indicate the impact of the nuisance parameters, $\Delta\mu$ (top axis), on the corresponding parameter of interest, $\mu$. They are calculated from the covariance matrix using the post-fit up (down) uncertainties of the parameter of interest and the nuisance parameter, $\sigma_{\mathrm{POI}}^{\mathrm{up}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{up}}$ ($\sigma_{\mathrm{POI}}^{\mathrm{down}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{down}}$), and their linear correlation coefficient, $\rho$, following EPJC 84 (2024) 593. The data points and error bars indicate the post-fit values and uncertainties of the nuisance parameters relative to their pre-fit values, $\theta_0$, in units of their pre-fit uncertainties, $\Delta \theta$ (bottom axis). A hollow circle indicates that the nuisance parameter has a Poisson constraint.

Nuisance parameters with a large impact on the fiducial measurement: ranking of the ten most impactful nuisance parameters for the $t\bar{t}+1c$ cross-section measurement. The colored boxes indicate the impact of the nuisance parameters, $\Delta\mu$ (top axis), on the corresponding parameter of interest, $\mu$. They are calculated from the covariance matrix using the post-fit up (down) uncertainties of the parameter of interest and the nuisance parameter, $\sigma_{\mathrm{POI}}^{\mathrm{up}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{up}}$ ($\sigma_{\mathrm{POI}}^{\mathrm{down}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{down}}$), and their linear correlation coefficient, $\rho$, following EPJC 84 (2024) 593. The data points and error bars indicate the post-fit values and uncertainties of the nuisance parameters relative to their pre-fit values, $\theta_0$, in units of their pre-fit uncertainties, $\Delta \theta$ (bottom axis). A hollow circle indicates that the nuisance parameter has a Poisson constraint.

Nuisance parameters with a large impact on the fiducial measurement: ranking of the ten most impactful nuisance parameters for the $t\bar{t}+{\geq}2c$ ratio to total $t\bar{t}+\text{jets}$ production. The colored boxes indicate the impact of the nuisance parameters, $\Delta\mu$ (top axis), on the corresponding parameter of interest, $\mu$. They are calculated from the covariance matrix using the post-fit up (down) uncertainties of the parameter of interest and the nuisance parameter, $\sigma_{\mathrm{POI}}^{\mathrm{up}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{up}}$ ($\sigma_{\mathrm{POI}}^{\mathrm{down}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{down}}$), and their linear correlation coefficient, $\rho$, following EPJC 84 (2024) 593. The data points and error bars indicate the post-fit values and uncertainties of the nuisance parameters relative to their pre-fit values, $\theta_0$, in units of their pre-fit uncertainties, $\Delta \theta$ (bottom axis). A hollow circle indicates that the nuisance parameter has a Poisson constraint.

Nuisance parameters with a large impact on the fiducial measurement: ranking of the ten most impactful nuisance parameters for the $t\bar{t}+1c$ ratio to total $t\bar{t}+\text{jets}$ production. The colored boxes indicate the impact of the nuisance parameters, $\Delta\mu$ (top axis), on the corresponding parameter of interest, $\mu$. They are calculated from the covariance matrix using the post-fit up (down) uncertainties of the parameter of interest and the nuisance parameter, $\sigma_{\mathrm{POI}}^{\mathrm{up}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{up}}$ ($\sigma_{\mathrm{POI}}^{\mathrm{down}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{down}}$), and their linear correlation coefficient, $\rho$, following EPJC 84 (2024) 593. The data points and error bars indicate the post-fit values and uncertainties of the nuisance parameters relative to their pre-fit values, $\theta_0$, in units of their pre-fit uncertainties, $\Delta \theta$ (bottom axis). A hollow circle indicates that the nuisance parameter has a Poisson constraint.

Efficiencies and rejection rates for the $b$- and $c$-jet tagging working points of the $b/c$-tagger, as well as the selection cuts on the discriminants, $\mathcal{D}_b'$ and $\mathcal{D}_c'$ that define the tagging bins. The primes indicate that a standard logistic function is applied to the discriminants. The b@70% and c@22% working points are inclusive of their respective tighter working points. The values were estimated from single-lepton and dilepton $t\bar{t}$ events simulated with Powheg + Pythia 8.

Efficiencies and rejection rates for the $b$- and $c$-jet tagging working points of the $b/c$-tagger, as well as the selection cuts on the discriminants, $\mathcal{D}_b'$ and $\mathcal{D}_c'$ that define the tagging bins. The primes indicate that a standard logistic function is applied to the discriminants. The b@70% and c@22% working points are inclusive of their respective tighter working points. The values were estimated from single-lepton and dilepton $t\bar{t}$ events simulated with Powheg + Pythia 8.

Pre-fit event yields in the twelve single-lepton and dilepton control regions (CRs) of the analysis. The quoted uncertainties include both statistical and systematic components, but all normalization factors are fixed to their pre-fit values with no uncertainties. All uncertainties are assumed to be uncorrelated. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes. The "Fakes" category includes events from misidentified leptons and non-prompt leptons, estimated through the data-driven matrix method in the single-lepton channel and through MC predictions in the dilepton channel.

Pre-fit event yields in the seven signal regions (SRs) of the analysis. The quoted uncertainties include both statistical and systematic components, but all normalization factors are fixed to their pre-fit values with no uncertainties. All uncertainties are assumed to be uncorrelated. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes. The "Fakes" category includes events from misidentified leptons and non-prompt leptons, estimated through the data-driven matrix method in the single-lepton channel and through MC predictions in the dilepton channel.

Post-fit event yields in the twelve single-lepton and dilepton control regions (CRs) of the analysis. The quoted uncertainties include both statistical and systematic components, and post-fit correlations between uncertainties are taken into account. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes. The "Fakes" category includes events from misidentified leptons and non-prompt leptons, estimated through the data-driven matrix method in the single-lepton channel and through MC predictions in the dilepton channel.

Post-fit event yields in the seven signal regions (SRs) of the analysis. The quoted uncertainties include both statistical and systematic components, and post-fit correlations between uncertainties are taken into account. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes. The "Fakes" category includes events from misidentified leptons and non-prompt leptons, estimated through the data-driven matrix method in the single-lepton channel and through MC predictions in the dilepton channel.

Pre-fit background composition of the SR$3\ell$ Prod. The table shows the event yields as opposed to just the percentages of the relevant background processes.

Post-fit plot of $H_\text{T}(\text{jets})$ in the SR$2\ell$ Dec from a signal-blinded background-only fit.

Post-fit plot of $m(t_\text{SM}, b\text{-jet}_0)$ in the SR$2\ell$ Prod from a signal-blinded background-only fit.

Post-fit plot of $m(\ell_\text{OS},\ell_\text{SS,1})$ in the SR$3\ell$ Dec from a signal-blinded background-only fit.

Post-fit plot of $m(\ell_\text{OS},\ell_\text{SS,1})$ in the SR$3\ell$ Prod from a signal-blinded background-only fit.

Post-fit plot of $D_\text{NN}(tHc)$ in the SR$2\ell$ Dec from the full fit to data.

Post-fit plot of $D_\text{NN}(tHc)$ in the SR$2\ell$ Prod from the full fit to data.

Post-fit plot of $D_\text{NN}(tHc)$ in the SR$3\ell$ Dec from the full fit to data.

Post-fit plot of $D_\text{NN}(tHc)$ in the SR$3\ell$ Prod from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ HF$e$ from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ HF$\mu$ from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ $t\bar{t}V$ from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_2)$ in the CR$3\ell$ HF$e$ from the full fit to data.

Post-fit plot of $p_\text{T}(\ell_2)$ in the CR$3\ell$ HF$\mu$ from the full fit to data.

Post-fit plot of $p_\text{T}(b\text{-jet}_0)$ in the CR$3\ell$ $t\bar{t}W$ from the full fit to data.

Post-fit plot of $p_\text{T}(b\text{-jet}_0)$ in the CR$3\ell$ $t\bar{t}Z$ from the full fit to data.

Observed and expected upper exclusion limits on the branching ratio $\mathcal{B}(t\to Hu)$ for different analyses and their statistical combination.

Observed and expected upper exclusion limits on the branching ratio $\mathcal{B}(t\to Hc)$ for different analyses and their statistical combination.

Post-fit normalisation factors of free-floating background processes and the signal normalisation.

Post-fit predicted and observed yields in all $2\ell$SS signal and control regions. Pre-fit signal contributions for a signal cross section equivalent to $\mathcal{B}(t\to Hq)=0.1\,\%$ are given as well.

Post-fit predicted and observed yields in all $3\ell$ signal and control regions. Pre-fit signal contributions for a signal cross section equivalent to $\mathcal{B}(t\to Hq)=0.1\,\%$ are given as well.

Expected upper limits on $\mathcal{B}(t\to Hq)$ for the nominal fit and alternative fit configurations. One contains the full phase space but only considers statistical uncertainties. Two other configurations consider the full set of systematic uncertainties, but only encompass one final state.

Expected and observed upper limits on $\mathcal{B}(t\to Hq)$ and $|C_{u\phi}^{qt,tq}|$ for the full fit containing all systematic uncertainties.

Pre-fit plot of $H_\text{T}(\text{jets})$ in the SR$2\ell$ Dec from a signal-blinded background-only fit.

Pre-fit plot of $m(t_\text{SM}, b\text{-jet}_0)$ in the SR$2\ell$ Prod from a signal-blinded background-only fit.

Pre-fit plot of $m(\ell_\text{OS},\ell_\text{SS,1})$ in the SR$3\ell$ Dec from a signal-blinded background-only fit.

Pre-fit plot of $m(\ell_\text{OS},\ell_\text{SS,1})$ in the SR$3\ell$ Prod from a signal-blinded background-only fit.

Post-fit plot of $D_\text{NN}(tHu)$ in the SR$2\ell$ Dec from the full fit to data.

Post-fit plot of $D_\text{NN}(tHu)$ in the SR$2\ell$ Prod from the full fit to data.

Post-fit plot of $D_\text{NN}(tHu)$ in the SR$3\ell$ Dec from the full fit to data.

Post-fit plot of $D_\text{NN}(tHu)$ in the SR$3\ell$ Prod from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHc)$ in the SR$2\ell$ Dec from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHc)$ in the SR$2\ell$ Prod from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHc)$ in the SR$3\ell$ Dec from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHc)$ in the SR$3\ell$ Prod from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHu)$ in the SR$2\ell$ Dec from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHu)$ in the SR$2\ell$ Prod from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHu)$ in the SR$3\ell$ Dec from the full fit to data.

Pre-fit plot of $D_\text{NN}(tHu)$ in the SR$3\ell$ Prod from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ HF$e$ from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ HF$\mu$ from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_1)$ in the CR$2\ell$ $t\bar{t}V$ from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_2)$ in the CR$3\ell$ HF$e$ from the full fit to data.

Pre-fit plot of $p_\text{T}(\ell_2)$ in the CR$3\ell$ HF$\mu$ from the full fit to data.

Pre-fit plot of $p_\text{T}(b\text{-jet}_0)$ in the CR$3\ell$ $t\bar{t}W$ from the full fit to data.

Pre-fit plot of $p_\text{T}(b\text{-jet}_0)$ in the CR$3\ell$ $t\bar{t}Z$ from the full fit to data.

Ranking of fit nuisance parameters according to their impact on the post-fit $tHu$ signal normalisation when fixed to $\pm1\sigma$

Ranking of fit nuisance parameters according to their impact on the post-fit $tHc$ signal normalisation when fixed to $\pm1\sigma$

Expected upper exclusion limits on the branching ratio $\mathcal{B}(t\to Hu)$ for each individual final state and the full analysis.

Expected upper exclusion limits on the branching ratio $\mathcal{B}(t\to Hc)$ for each individual final state and the full analysis.

$D^+ D^0$~mass distributions for selected candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^+ D^+$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^+ D^0 \pi^+$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Background-subtracted distributions for the multiplicity of tracks reconstructed in the vertex detector for $T_{cc}^+\to D^0 D^0 \pi^+$ signal, low-mass $D^0\bar{D}^0$ and $D^0D^0$ pairs. The binning scheme is chosen to have an approximately uniform distribution for $D^0\bar{D}^0$ pairs. The distributions for the $D^0\bar{D}^0$ and $D^0D^0$ pairs are normalised to the same yields as the $T_{cc}^+\to D^0 D^0 \pi^+$ signal. For better visualisation, the points are slightly displaced from the bin centres. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Background-subtracted transverse momentum spectra for $T_{cc}^+\to D^0 D^0 \pi^+$ signal, low-mass $D^0\bar{D}^0$ and $D^0D^0$ pairs. The binning scheme is chosen to have an approximately uniform distribution for $D^0\bar{D}^0$ pairs. The distributions for the $D^0\bar{D}^0$ and $D^0D^0$ pairs are normalised to the same yields as the $T_{cc}^+\to D^0 D^0 \pi^+$ signal. For better visualisation, the points are slightly displaced from the bin centres. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^0 \bar{D}^0$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^0 D^0$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $\bar{D}^0D^0\pi^+$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^0D^0\pi^+$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^0D^-$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Mass distributions for selected $D^0D^+$ candidates with the $D^0$ background subtracted. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

Distribution of $D^0 D^0 \pi^+$ mass where the contribution of the non-$D^0$ background has been statistically subtracted by assigning the a weight to every candidate.

Mass distribution for $D^0 \pi^+$ pairs from selected $D^0 D^0 \pi^+$ candidates with a mass below the $D^{*+}D^0$ mass threshold with non-$D^0$ background subtracted by assigning the a weight to every candidate.

$D^0 D^0$~mass distributions for selected candidates with the $D^0$ background subtracted by assigning the a weight to every candidate.

$D^+ D^0$~mass distributions for selected candidates with the $D^0$ background subtracted by assigning the a weight to every candidate.

Mass distributions for selected $D^+ D^+$ candidates with the $D^0$ background subtracted by assigning the a weight to every candidate.

Mass distributions for selected $D^+ D^0 \pi^+$ candidates with the $D^0$ background subtracted by assigning the a weight to every candidate.

Spline interpolation of the weight function (depending of $\cos \theta_{P_c}$) applied to $\Lambda_{b}^{0}$ candidates. The weight function is determined as the inverse of the density of $\Lambda_{b}^{0}$ candidates in the narrow pentaquark peak region while the interpolation is done between the bin centers.

HT inclusive distribution in SRZ.

Jet multiplicity distribution in SRZ.

B-tagged jet multiplicity distribution in SRZ.

Dilepton invariant mass distribution in SRMedium.

Dilepton invariant mass distribution in SRMedium.

Dilepton invariant mass in SRHigh.

Expected 95% exclusion contour for the gluino simplified model with fixed N1 mass and decays to quarks and Z bosons in SRZ.

Observed 95% exclusion contour for the gluino simplified model with fixed N1 mass and decays to quarks and Z bosons in SRZ.

Expected 95% exclusion contour for the squark simplified model with fixed N1 mass and decays to quarks and Z bosons in SRZ.

Observed 95% exclusion contour for the squark simplified model with fixed N1 mass and decays to quarks and Z bosons in SRZ.

Expected 95% exclusion contour for the gluino simplified model with varying N1 mass and decays to quarks and Z bosons in SRZ.

Observed 95% exclusion contour for the gluino simplified model with varying N1 mass and decays to quarks and Z bosons in SRZ.

Expected 95% exclusion contour for the gluino simplified model with sleptons in the edge SRs.

Observed 95% exclusion contour for the gluino simplified model with sleptons in the edge SRs.

Expected 95% exclusion contour for the gluino simplified model with decays to on- or off-shell Z bosons in the edge SRs.

Observed 95% exclusion contour for the gluino simplified model with decays to on- or off-shell Z bosons in the edge SRs.

Expected 95% exclusion contour for the gluino simplified model with decays to on- or off-shell Z bosons in SRZ.

Observed 95% exclusion contour for the gluino simplified model with decays to on- or off-shell Z bosons in SRZ.

Cutflow table for three benchmark signal points from the $\tilde{g}$--$\chi_{2}^{0}$ on-shell grid in SRZ for the $ee$ and $\mu\mu$ channels separately. The number of generated events is N RAW, while N WEIGHTED corresponds to the number of events expected in 14.7/fb, after all MC simulation corrections have been applied.

Cutflow table for three benchmark signal points for the slepton model in SR-low, SR-medium and SR-high for the $ee$ and $\mu\mu$ channels separately. The number of generated events is N RAW, while N WEIGHTED corresponds to the number of events expected in 14.7/fb, after all MC simulation corrections have been applied. The event generation was filtered to require at least two leptons (electrons or muons) with $p_{T}>5$ GeV and $|\eta|<2.8$, with 20000 events required to pass the filter. The corresponding number of unfiltered events is provided in the ``Pre-filter events'' row of the table.

Signal acceptance for the gluino simplified model with decays to quarks and Z bosons in SRZ.

Signal efficiency for the gluino simplified model with decays to quarks and Z bosons in SRZ.

Signal acceptance for the gluino simplified model with decays via sleptons in SRLow.

Signal efficiency for the gluino simplified model with decays via sleptons in SRLow.

Signal acceptance for the gluino simplified model with decays via sleptons in SRMedium.

Signal efficiency for the gluino simplified model with decays via sleptons in SRMedium.

Signal acceptance for the gluino simplified model with decays via sleptons in SRHigh.

Signal efficiency for the gluino simplified model with decays via sleptons in SRHigh.

Upper limits on the signal cross-section at 95% CL for the gluino simplified model with fixed N1 mass and decays to quarks and Z bosons in SRZ.

Upper limits on the signal cross-section at 95% CL for the squark simplified model with fixed N1 mass and decays to quarks and Z bosons in SRZ.

Upper limits on the signal cross-section at 95% CL for the gluino simplified model with varying N1 mass and decays to quarks and Z bosons in SRZ.

Upper limits on the signal cross-section at 95% CL for the gluino simplified model with decays via sleptons.

Upper limits on the signal cross-section at 95% CL for the gluino simplified model with decays to on-shell or off-shell Z bosons.

Observed and predicted W+W- transverse mass distribution in the E-E channel. Also tabulated are the predictions for a RS graviton of mass 1000 GeV and a bulk RS graviton with mass 600 GeV.

Observed and predicted W+W- transverse mass distribution in the E-MU channel. Also tabulated are the predictions for a RS graviton of mass 1000 GeV and a bulk RS graviton with mass 600 GeV.

Correlations in the PT distribution.

Normalised fiducial cross section in bins of the Mass of the WZ.

Correlations in the Mass distribution.

The meff_incl distribution in Signal Region C.

The meff_incl distribution in Signal Region D.

The meff_incl distribution in Signal Region E.