The d$^{2}\sigma$/dydp$_{\mathrm{T}}$ production cross section of D$^{0}$ for $8 < p_{\mathrm{T}} < 12$ GeV in ultraperipheral PbPb collisions.

Measured fiducial cross section of the VBF production mode.

Correlation matrix for fit configuration 1

Correlation matrix for fit configuration 2

Correlation matrix for fit configuration 3

The measured differences of mass differences between vector and ground states.

The measured ratios of mass differences between vector and ground states.

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.

Likelihood scan of cHuRe1122 versus cHuRe33. Other Wilson coefficients are fixed to zero.

Likelihood scan of cHdRe1122 versus cHdRe33. Other Wilson coefficients are fixed to zero.

Likelihood scan of cW versus cWtil. Other Wilson coefficients are fixed to zero.

Likelihood scan of cHqMRe1122 versus cHqMRe33. Other Wilson coefficients are profiled as well.

Likelihood scan of cHq3MRe1122 versus cHq3MRe33. Other Wilson coefficients are profiled as well.

Likelihood scan of cHuRe1122 versus cHuRe33. Other Wilson coefficients are profiled as well.

Likelihood scan of cHdRe1122 versus cHdRe33. Other Wilson coefficients are profiled as well.

Likelihood scan of cW versus cWtil. Other Wilson coefficients are profiled as well.

Signal efficiencies with Full Run 2 dimuon channel in 1b final state for different bbll (destructive interference) signal scenarios

Signal efficiencies with Full Run 2 dimuon channel in 2b final state for different bbll signal scenarios

Signal efficiencies with Full Run 2 dimuon channel in 2b final state for different bbll (destructive interference) signal scenarios

Signal efficiencies with Full Run 2 dielectron channel for different bbll signal scenarios

Signal efficiencies with Full Run 2 dielectron channel for different bbll (destructive interference) signal scenarios

Signal efficiencies with Full Run 2 dielectron channel in 1b final state for different bbll signal scenarios

Signal efficiencies with Full Run 2 dielectron channel in 1b final state for different bbll (destructive interference) signal scenarios

Signal efficiencies with Full Run 2 dielectron channel in 2b final state for different bbll signal scenarios

Signal efficiencies with Full Run 2 dielectron channel in 2b final state for different bbll (destructive interference) signal scenarios

Signal efficiencies with Full Run 2 dimuon channel for bsll signal

Signal efficiencies with Full Run 2 dielectron channel for bsll signal

Observed and expected background yields for different mass ranges in the dielectron channel in 0 b jet final state. The sum of all background contributions is shown as well as a breakdown into the four main categories. The quoted uncertainty includes both the statistical and the systematic components.

Observed and expected background yields for different mass ranges in the dielectron channel in 1 b jet final state. The sum of all background contributions is shown as well as a breakdown into the four main categories. The quoted uncertainty includes both the statistical and the systematic components.

Observed and expected background yields for different mass ranges in the dielectron channel in 2 b jet final state. The sum of all background contributions is shown as well as a breakdown into the four main categories. The quoted uncertainty includes both the statistical and the systematic components.

Observed and expected background yields for different mass ranges in the dimuon channel in 0 b jet final state. The sum of all background contributions is shown as well as a breakdown into the four main categories. The quoted uncertainty includes both the statistical and the systematic components.

Observed and expected background yields for different mass ranges in the dimuon channel in 1 b jet final state. The sum of all background contributions is shown as well as a breakdown into the four main categories. The quoted uncertainty includes both the statistical and the systematic components.

Observed and expected background yields for different mass ranges in the dimuon channel in 2 b jet final state. The sum of all background contributions is shown as well as a breakdown into the four main categories. The quoted uncertainty includes both the statistical and the systematic components.

Measured flavor ratio from DY+jets MC and data in 0b channel. The quoted uncertainty includes both the statistical and the systematic components. The flavor ratio from DY+jets MC deviates slightly from unity as the ratio is not corrected for mass-dependent acceptance times efficiency.

Measured flavor ratio from DY+jets MC and data in 1b+2b channel. The quoted uncertainty includes both the statistical and the systematic components. The flavor ratio from DY+jets MC deviates slightly from unity as the ratio is not corrected for mass-dependent acceptance times efficiency.

Lower limits on the energy scale ($\Lambda$) at 95% CL for bbll signal with different chirality and interference assumptions in the dielectron channel combining all b jet final states. The limits are obtained for $m_{ee}>300$ GeV.

Lower limits on the energy scale ($\Lambda$) at 95% CL for bbll signal with different chirality and interference assumptions in the dimuon channel combining all b jet final states. The limits are obtained for $m_{\mu\mu}>300$ GeV.

Lower limits on the energy scale ($\Lambda$) at 95% CL for bbll signal with different chirality and interference assumptions in the combined (dimuon+dielectron) channel combining all b jet final states. The limits are obtained for $m_{\ell\ell}>300$ GeV.

Observed and Expected upper limits on the product of the production cross section for the bsll signal in the dielectron and dimuon channels for 0 b-tagged jets.

Observed and Expected upper limits on the product of the production cross section for the bsll signal in the dielectron and dimuon channels for 1 b-tagged jets.

Cutflows and signal efficiencies for the Stealth SYY SUSY model in the $1\ell$ channel corresponding to two values of $m_{\tilde{t}}$.

Cutflows and signal efficiencies for the RPV SUSY model in the $2\ell$ channel corresponding to two values of $m_{\tilde{t}}$.

Cutflows and signal efficiencies for the Stealth SYY SUSY model in the $2\ell$ channel corresponding to two values of $m_{\tilde{t}}$.

Signal efficiencies for the RPV SUSY model in the $0\ell$ channel corresponding to the four ABCD regions

Signal efficiencies for the SYY SUSY model in the $0\ell$ channel corresponding to the four ABCD regions

Signal efficiencies for the RPV SUSY model in the $1\ell$ channel corresponding to the four ABCD regions

Signal efficiencies for the SYY SUSY model in the $1\ell$ channel corresponding to the four ABCD regions

Signal efficiencies for the RPV SUSY model in the $2\ell$ channel corresponding to the four ABCD regions

Signal efficiencies for the SYY SUSY model in the $2\ell$ channel corresponding to the four ABCD regions

Expected and observed 95% CL upper limit on the top squark pair production cross section as a function of the top squark mass for the RPV SUSY model in the $0\ell$ channel

Expected and observed 95% CL upper limit on the top squark pair production cross section as a function of the top squark mass for the RPV SUSY model in the $1\ell$ channel

Expected and observed 95% CL upper limit on the top squark pair production cross section as a function of the top squark mass for the RPV SUSY model in the $2\ell$ channel

Expected and observed 95% CL upper limit on the top squark pair production cross section as a function of the top squark mass for the RPV SUSY model in the Combo channel

The $N_{jets}$ distributions from the background-only fits to data in the $0\ell$. The signal distribution overlaid corresponds to the RPV model with $m_{\tilde{t}} = 400$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $1\ell$. The signal distribution overlaid corresponds to the RPV model with $m_{\tilde{t}} = 400$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $2\ell$. The signal distribution overlaid corresponds to the RPV model with $m_{\tilde{t}} = 400$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $0\ell$. The signal distribution overlaid corresponds to the RPV model with $m_{\tilde{t}} = 800$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $1\ell$. The signal distribution overlaid corresponds to the RPV model with $m_{\tilde{t}} = 800$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $2\ell$. The signal distribution overlaid corresponds to the RPV model with $m_{\tilde{t}} = 800$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

Expected and observed 95% CL upper limit on the top squark pair production cross section as a function of the top squark mass for the SYY SUSY model in the $0\ell$ channel

Expected and observed 95% CL upper limit on the top squark pair production cross section as a function of the top squark mass for the SYY SUSY model in the $1\ell$ channel

Expected and observed 95% CL upper limit on the top squark pair production cross section as a function of the top squark mass for the SYY SUSY model in the $2\ell$ channel

Expected and observed 95% CL upper limit on the top squark pair production cross section as a function of the top squark mass for the SYY SUSY model in the Combo channel

The $N_{jets}$ distributions from the background-only fits to data in the $0\ell$. The signal distribution overlaid corresponds to the Stealth SYY model with $m_{\tilde{t}} = 400$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $1\ell$. The signal distribution overlaid corresponds to the Stealth SYY model with $m_{\tilde{t}} = 400$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $2\ell$. The signal distribution overlaid corresponds to the Stealth SYY model with $m_{\tilde{t}} = 400$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $0\ell$. The signal distribution overlaid corresponds to the Stealth SYY model with $m_{\tilde{t}} = 800$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $1\ell$. The signal distribution overlaid corresponds to the Stealth SYY model with $m_{\tilde{t}} = 800$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $2\ell$. The signal distribution overlaid corresponds to the Stealth SYY model with $m_{\tilde{t}} = 800$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

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

Fiducial differential cross-sections as a function of $p_{\mathrm{T}}^{\mathrm{sub-lead.\,lep.}}$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

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

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

Fiducial differential cross-sections as a function of $p_{\mathrm{T},e\mu}$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

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

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

Fiducial differential cross-sections as a function of $|y_{e\mu}|$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

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

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

Fiducial differential cross-sections as a function of $m_{e\mu}$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

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

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

Fiducial differential cross-sections as a function of $|\Delta\phi_{e\mu}|$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $|\Delta\phi_{e\mu}|$.

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $|\Delta\phi_{e\mu}|$.

Fiducial differential cross-sections as a function of $\cos(\theta^{*})$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $\cos(\theta^{*})$.

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $\cos(\theta^{*})$.

Fiducial differential cross-sections as a function of $E_{\mathrm{T}}^{\mathrm{miss}}$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $E_{\mathrm{T}}^{\mathrm{miss}}$.

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $E_{\mathrm{T}}^{\mathrm{miss}}$.

Fiducial differential cross-sections as a function of $H_{\mathrm{T}}^{\mathrm{lep}+\mathrm{MET}}$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

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

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

Fiducial differential cross-sections as a function of $m_{\mathrm{T},e\mu}$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

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

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

Fiducial differential cross-sections as a function of the number of jets. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable Number of jets ($p_{\mathrm{T}} > $ 30 GeV).

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable Number of jets ($p_{\mathrm{T}} > $ 30 GeV).

Fiducial differential cross-sections as a function of $S_{\mathrm{T}}$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

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

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

Fiducial differential cross-sections for $p_{\mathrm{T},e\mu}$ in the region with 85 $<m_{e\mu}<$ 150, for the nNNLO MATRIX prediction and various four flavour PDF sets. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The measurement is compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1 using four different four-flavour NNLO PDF sets: NNPDF3.0, NNPDF3.1, CT18, and MSHT20.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T},e\mu}$ (85 $<m_{e\mu}<$ 150).

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T},e\mu}$ (85 $<m_{e\mu}<$ 150).

Fiducial differential cross-sections for $|y_{e\mu}|$ in the region with 85 $<m_{e\mu}<$ 150, for the nNNLO MATRIX prediction and various four flavour PDF sets. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The measurement is compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1 using four different four-flavour NNLO PDF sets: NNPDF3.0, NNPDF3.1, CT18, and MSHT20.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $|y_{e\mu}|$ (85 $<m_{e\mu}<$ 150).

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $|y_{e\mu}|$ (85 $<m_{e\mu}<$ 150).

Fiducial differential cross-sections for $p_{\mathrm{T},e\mu}$ in the region with 150 $<m_{e\mu}<$ 250, for the nNNLO MATRIX prediction and various four flavour PDF sets. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The measurement is compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1 using four different four-flavour NNLO PDF sets: NNPDF3.0, NNPDF3.1, CT18, and MSHT20.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T},e\mu}$ (150 $<m_{e\mu}<$ 250).

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T},e\mu}$ (150 $<m_{e\mu}<$ 250).

Fiducial differential cross-sections for $|y_{e\mu}|$ in the region with 150 $<m_{e\mu}<$ 250, for the nNNLO MATRIX prediction and various four flavour PDF sets. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The measurement is compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1 using four different four-flavour NNLO PDF sets: NNPDF3.0, NNPDF3.1, CT18, and MSHT20.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $|y_{e\mu}|$ (150 $<m_{e\mu}<$ 250).

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $|y_{e\mu}|$ (150 $<m_{e\mu}<$ 250).

Fiducial differential cross-sections for $p_{\mathrm{T},e\mu}$ in the region with $m_{e\mu}>$ 250, for the nNNLO MATRIX prediction and various four flavour PDF sets. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The measurement is compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1 using four different four-flavour NNLO PDF sets: NNPDF3.0, NNPDF3.1, CT18, and MSHT20.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T},e\mu}$ ($m_{e\mu}>$ 250).

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T},e\mu}$ ($m_{e\mu}>$ 250).

Fiducial differential cross-sections for $|y_{e\mu}|$ in the region with $m_{e\mu}>$ 250, for the nNNLO MATRIX prediction and various four flavour PDF sets. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The measurement is compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1 using four different four-flavour NNLO PDF sets: NNPDF3.0, NNPDF3.1, CT18, and MSHT20.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $|y_{e\mu}|$ ($m_{e\mu}>$ 250).

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $|y_{e\mu}|$ ($m_{e\mu}>$ 250).

Measured value of the charge asymmetry, $A_{C}$.

Charge asymmetry $A_C$ measured as a function of the absolute difference of the absolute lepton rapidities $||\eta_{l+}|-|\eta_{l-}||$. The measured values are shown as points, whose error bars represent the statistical uncertainty, while the solid bands indicate the size of the total uncertainty. In the bottom panel, the central value of the charge asymmetry divided by the total uncertainty of the $A_C$ measurement is shown. The measurement is compared with the NNLO+PS predictions from Powheg MiNNLO+Pythia8 with vertical lines denoting PDF uncertainties.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $A_{C}$ vs. $||\eta_{l+}|-|\eta_{l-}||$.

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $A_{C}$ vs. $||\eta_{l+}|-|\eta_{l-}||$.

Charge asymmetry $A_C$ measured as a function of the invariant mass of the dilepton system $m_{e\mu}$. The measured values are shown as points, whose error bars represent the statistical uncertainty, while the solid bands indicate the size of the total uncertainty. In the bottom panel, the central value of the charge asymmetry divided by the total uncertainty of the $A_C$ measurement is shown. The measurement is compared with the NNLO+PS predictions from Powheg MiNNLO+Pythia8 with vertical lines denoting PDF uncertainties.

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

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

Measured value of the asymmetry in the CP-odd observable $O_{z}$, $A_{CP}$.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the 1 b jet, 0 forward jet, $t\leq 0$ signal region in the single lepton channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the 1 b jet, 0 forward jet, $t>0$ signal region in the single lepton channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the 1 b jet, $\geq 1$ forward jet, $t\leq 0$ signal region in the single lepton channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the 1 b jet, $\geq 1$ forward jet, $t>0$ signal region in the single lepton channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the $\geq 2$ b jet, $t\leq 0$ signal region in the single lepton channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the $\geq 2$ b jet, $t>0$ signal region in the single lepton channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit NN output score distribution of the $\geq 2$ b jet, different flavor signal region in the dileptonic channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit NN output score distribution of the $\geq 2$ b jet, same flavor signal region in the dileptonic channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit NN output score distribution of the 1 b jet, different flavor signal region in the dileptonic channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit NN output score distribution of the 1 b jet, same flavor signal region in the dileptonic channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The model-independent 95\% CL limits on production cross section for new physics processes for the scalar interaction. The expected limit is shown by the black dashed line with the 68 and 95\% CL uncertainty bands in green and yellow, respectively, while the observed limit is shown by the solid black line. Theoretical LO cross section values for the DM model and their associated uncertainties are also presented (grey line).

The model-independent 95\% CL limits on production cross section for new physics processes for the pseudoscalar interaction. The expected limit is shown by the black dashed line with the 68 and 95\% CL uncertainty bands in green and yellow, respectively, while the observed limit is shown by the solid black line. Theoretical LO cross section values for the DM model and their associated uncertainties are also presented (grey line).

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the QCD control region in the all hadronic channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the tt(1l) control region in the all hadronic channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the W+jets control region in the all hadronic channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the Z(2l) control region in the all hadronic channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the W+jets control region in the single lepton channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the tt(2l) control region in the single lepton channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit $(N_{\text{jet}},N_{\text{bjet}})$ distribution of the ttZ control region in the dileptonic channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit NN distribution of the Drell-Yan control region in the dileptonic channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

Predicted signal yields for the 2016, 2017, and 2018 data-taking periods, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{\phi}) = (1,50)$ GeV scalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{\phi}) = (1,500)$ GeV scalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{a}) = (1,50)$ GeV pseudoscalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{a}) = (1,500)$ GeV pseudoscalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{\phi}) = (1,50)$ GeV scalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{\phi}) = (1,500)$ GeV scalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{a}) = (1,50)$ GeV pseudoscalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{a}) = (1,500)$ GeV pseudoscalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{\phi}) = (1,50)$ GeV scalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{\phi}) = (1,500)$ GeV scalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{a}) = (1,50)$ GeV pseudoscalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{a}) = (1,500)$ GeV pseudoscalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{\phi}) = (1,50)$ GeV scalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{\phi}) = (1,500)$ GeV scalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{a}) = (1,50)$ GeV pseudoscalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{a}) = (1,500)$ GeV pseudoscalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{\phi}) = (1,50)$ GeV scalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{\phi}) = (1,500)$ GeV scalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{a}) = (1,50)$ GeV pseudoscalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{a}) = (1,500)$ GeV pseudoscalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{\phi}) = (1,50)$ GeV scalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{\phi}) = (1,500)$ GeV scalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{a}) = (1,50)$ GeV pseudoscalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{a}) = (1,500)$ GeV pseudoscalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Cross sections calculated to leading order (LO) accuracy for the tt+DM scalar mediator signal process considering all possible decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the tt+DM pseudoscalar mediator signal process considering all possible decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the tt+DM scalar mediator signal process considering dileptonic decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the tt+DM pseudoscalar mediator signal process considering dileptonic decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the t+DM t-channel scalar mediator signal process considering all possible decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the t+DM tW-channel scalar mediator signal process considering all possible decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the t+DM t-channel pseudoscalar mediator signal process considering all possible decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the t+DM tW-channel pseudoscalar mediator signal process considering all possible decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the t+DM tW-channel scalar mediator signal process considering dileptonic decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the t+DM tW-channel pseudoscalar mediator signal process considering dileptonic decays from the visible particles in the final state.