The correlation matrix corresponding to the statistical uncertainties on the differential cross-section ($d\sigma/dp_T$) fit results for $W^-$. To combine with $W^+$ results (Fig8a), use the rows and columns ordered as $W^+$ and then $W^-$. Assume no correlation in the statistical uncertainties between $W^+$ and $W^-$ (zero entries in the off-diagonal blocks).

The correlation matrix corresponding to all of the systematic uncertainties on the differential cross-section ($d\sigma/dp_T$) fit results for $W^+$ and $W^-$ combined together, with rows and columns ordered as $W^+$ and then $W^-$.

Observed 95% $\text{CL}_\text{S}$ upper limits on the production cross-section times acceptance times branching ratio to jets, $\sigma \cdot A \cdot \text{BR}$, of Gaussian-shaped signals of 5%, 10%, and 15% width relative to their peak mass, $m_G$. Also included are the corresponding expected upper limits predicted for the case the $m_{jj}$ distribution is observed to be identical to the background prediction in each bin and the $1\sigma$ and $2\sigma$ envelopes of outcomes expected for Poisson fluctuations around the background expectation. Limits are derived from the J100 signal region.

Observed 95% $\text{CL}_\text{S}$ upper limits on the the value of the coupling to quarks, $g_q$, of the $Z^\prime$ model as a function of the resonance mass $m_{Z^\prime}$. Also included are the corresponding expected upper limits predicted for the case the $m_{jj}$ distribution is observed to be identical to the background prediction in each bin and the $1\sigma$ and $2\sigma$ envelopes of outcomes expected for Poisson fluctuations around the background expectation. Limits are derived from the J50 signal region.

Observed 95% $\text{CL}_\text{S}$ upper limits on the the value of the coupling to quarks, $g_q$, of the $Z^\prime$ model as a function of the resonance mass $m_{Z^\prime}$. Also included are the corresponding expected upper limits predicted for the case the $m_{jj}$ distribution is observed to be identical to the background prediction in each bin and the $1\sigma$ and $2\sigma$ envelopes of outcomes expected for Poisson fluctuations around the background expectation. Limits are derived from the J100 signal region.

Acceptance of simulated $Z^\prime$ signal events given the analysis event selection as a function of $m_{Z^\prime}$. The solid black line corresponds to the acceptance without a lower $m_{jj}$ requirement applied, the solid blue line includes the $m_{jj} > 344$ GeV requirement of the J50 signal region, and the solid red line includes the $m_{jj} > 481$ GeV requirement of the J100 signal region.

Numbers of events in each of the two considered datasets after consecutively applying the analysis selections.

Observed upper limits at 95% CL on the product of the cross section and branching fraction squared, $\sigma (\tilde{\mathrm{t}} \, \bar{\tilde{\mathrm{t}}}) \, \mathcal{B}^{2} ( \tilde{\mathrm{t}} \rightarrow \mathrm{c} \tilde{\chi}_1^0 )$ (left panel), are shown using the color scale where the $\tilde{\mathrm{t}}$ mass is on the $x$-axis and the mass difference between the $\tilde{\mathrm{t}}$ and the LSP is on the $y$-axis. The expected lower mass limits (magenta line) together with their $\pm 1\sigma$ uncertainties (magenta dashed lines) and the observed lower mass limits (black line) are indicated for 100% branching fractions.

Chargino-neutralino production. Observed upper limits at 95% CL on the product of the cross section and the two branching fractions, $\sigma (\tilde{\chi}_1^\pm \tilde{\chi}_2^0) \, \mathcal{B} ( \tilde{\chi}_1^\pm \rightarrow \mathrm{W}^{\pm} \tilde{\chi}_1^0 ) \, \mathcal{B} ( \tilde{\chi}_2^0 \rightarrow \mathrm{Z} \tilde{\chi}_1^0 )$, are shown using the color scale where the $\tilde{\chi}_1^\pm/\tilde{\chi}_2^0$ mass is on the $x$-axis and the mass difference between the $\tilde{\chi}_1^\pm/\tilde{\chi}_2^0$ and the LSP is on the $y$-axis. For these results, based on the TChiWZ simplified model, the $\tilde{\chi}_1^\pm$ and $\tilde{\chi}_2^0$ masses are set equal. The expected lower mass limits (magenta line) together with their $\pm 1\sigma$ uncertainties (magenta dashed lines) and the observed lower mass limits (black line) are indicated for 100% branching fractions for Wino-like cross-sections (left panel) and for Higgsino-like cross-sections (right panel).

Chargino-neutralino production. Observed upper limits at 95% CL on the product of the cross section and the two branching fractions, $\sigma (\tilde{\chi}_1^\pm \tilde{\chi}_2^0) \, \mathcal{B} ( \tilde{\chi}_1^\pm \rightarrow \mathrm{W}^{\pm} \tilde{\chi}_1^0 ) \, \mathcal{B} ( \tilde{\chi}_2^0 \rightarrow \mathrm{Z} \tilde{\chi}_1^0 )$, are shown using the color scale where the $\tilde{\chi}_1^\pm/\tilde{\chi}_2^0$ mass is on the $x$-axis and the mass difference between the $\tilde{\chi}_1^\pm/\tilde{\chi}_2^0$ and the LSP is on the $y$-axis. For these results, based on the TChiWZ simplified model, the $\tilde{\chi}_1^\pm$ and $\tilde{\chi}_2^0$ masses are set equal. The expected lower mass limits (magenta line) together with their $\pm 1\sigma$ uncertainties (magenta dashed lines) and the observed lower mass limits (black line) are indicated for 100% branching fractions for Wino-like cross-sections (left panel) and for Higgsino-like cross-sections (right panel).

Chargino pair production. The left panel shows the observed upper limits at 95% CL on the product of the cross section and the branching fraction squared, $\sigma (\tilde{\chi}_1^+ \tilde{\chi}_1^- ) \, \mathcal{B}^{2} ( \tilde{\chi}_1^\pm \rightarrow \mathrm{W}^{\pm} \tilde{\chi}_1^0 ) $ are shown using the color scale where the $\tilde{\chi}_1^\pm$ mass is on the $x$-axis and the mass difference between the $\tilde{\chi}_1^\pm $ and the LSP is on the $y$-axis. The expected lower mass limits (magenta line) together with their $\pm 1\sigma$ uncertainties (magenta dashed lines) and the observed lower mass limits (black line) are indicated for 100% branching fractions for Wino-like cross-sections. The right panel shows the results for chargino pair production with decays as in the TChiSlepSnu model with democratic decay via an intermediate sneutrino or charged slepton ($\widetilde{\ell}^{\pm}_{\mathrm{L}}$) with mass halfway between the chargino and the lightest neutralino. These model predictions also assume Wino-like cross sections.

Chargino pair production. The left panel shows the observed upper limits at 95% CL on the product of the cross section and the branching fraction squared, $\sigma (\tilde{\chi}_1^+ \tilde{\chi}_1^- ) \, \mathcal{B}^{2} ( \tilde{\chi}_1^\pm \rightarrow \mathrm{W}^{\pm} \tilde{\chi}_1^0 ) $ are shown using the color scale where the $\tilde{\chi}_1^\pm$ mass is on the $x$-axis and the mass difference between the $\tilde{\chi}_1^\pm $ and the LSP is on the $y$-axis. The expected lower mass limits (magenta line) together with their $\pm 1\sigma$ uncertainties (magenta dashed lines) and the observed lower mass limits (black line) are indicated for 100% branching fractions for Wino-like cross-sections. The right panel shows the results for chargino pair production with decays as in the TChiSlepSnu model with democratic decay via an intermediate sneutrino or charged slepton ($\widetilde{\ell}^{\pm}_{\mathrm{L}}$) with mass halfway between the chargino and the lightest neutralino. These model predictions also assume Wino-like cross sections.

Slepton pair production. Observed 95% CL upper limits on the product of the cross section and branching fraction squared for direct slepton pair production followed by decay of both sleptons to the corresponding lepton and neutralino (color scale). Slepton $\tilde{l}_{\mathrm{L}/\mathrm{R}}$ indicates the scalar supersymmetric partner of left- and right-handed electrons and muons. The limit is shown as a function of the slepton mass and the mass difference between the slepton and the lightest neutralino. The regions to the left of the lines denote the regions excluded for a branching fraction of 100%. The median expected exclusion regions for 100% branching fraction are delimited by the dashed lines.

Slepton pair production. Observed 95% CL upper limits on the product of the cross section and branching fraction squared for direct slepton pair production followed by decay of both sleptons to the corresponding lepton and neutralino (color scale). Slepton $\tilde{l}_{\mathrm{L}/\mathrm{R}}$ indicates the scalar supersymmetric partner of left- and right-handed electrons and muons. The limit is shown as a function of the slepton mass and the mass difference between the slepton and the lightest neutralino. The regions to the left of the lines denote the regions excluded for a branching fraction of 100%. The median expected exclusion regions for 100% branching fraction are delimited by the dashed lines.

Slepton pair production. Observed 95% CL upper limits on the product of the cross section and branching fraction squared for direct slepton pair production followed by decay of both sleptons to the corresponding lepton and neutralino (color scale). Slepton $\tilde{l}_{\mathrm{L}/\mathrm{R}}$ indicates the scalar supersymmetric partner of left- and right-handed electrons and muons. The limit is shown as a function of the slepton mass and the mass difference between the slepton and the lightest neutralino. The regions to the left of the lines denote the regions excluded for a branching fraction of 100%. The median expected exclusion regions for 100% branching fraction are delimited by the dashed lines.

Slepton pair production. Observed 95% CL upper limits on the cross section times branching fraction squared for direct selectron pair production (left) and smuon pair production (right) followed by decay of both sleptons to the corresponding lepton and neutralino (color scale). The limits are shown as a function of the slepton mass and the mass difference between the slepton and the lightest neutralino for the three different simplified possibilities of only RR, only LL, and both RR and LL where it is assumed that the R and L masses are identical. The regions to the left of the lines denote the regions excluded for a branching fraction of 100%. Median expected limits for 100% branching fraction are delimited by the dashed lines.

Slepton pair production. Observed 95% CL upper limits on the cross section times branching fraction squared for direct selectron pair production (left) and smuon pair production (right) followed by decay of both sleptons to the corresponding lepton and neutralino (color scale). The limits are shown as a function of the slepton mass and the mass difference between the slepton and the lightest neutralino for the three different simplified possibilities of only RR, only LL, and both RR and LL where it is assumed that the R and L masses are identical. The regions to the left of the lines denote the regions excluded for a branching fraction of 100%. Median expected limits for 100% branching fraction are delimited by the dashed lines.

Slepton pair production. Observed 95% CL upper limits on the cross section times branching fraction squared for direct selectron pair production (left) and smuon pair production (right) followed by decay of both sleptons to the corresponding lepton and neutralino (color scale). The limits are shown as a function of the slepton mass and the mass difference between the slepton and the lightest neutralino for the three different simplified possibilities of only RR, only LL, and both RR and LL where it is assumed that the R and L masses are identical. The regions to the left of the lines denote the regions excluded for a branching fraction of 100%. Median expected limits for 100% branching fraction are delimited by the dashed lines.

Slepton pair production. Observed 95% CL upper limits on the cross section times branching fraction squared for direct selectron pair production (left) and smuon pair production (right) followed by decay of both sleptons to the corresponding lepton and neutralino (color scale). The limits are shown as a function of the slepton mass and the mass difference between the slepton and the lightest neutralino for the three different simplified possibilities of only RR, only LL, and both RR and LL where it is assumed that the R and L masses are identical. The regions to the left of the lines denote the regions excluded for a branching fraction of 100%. Median expected limits for 100% branching fraction are delimited by the dashed lines.

Slepton pair production. Observed 95% CL upper limits on the cross section times branching fraction squared for direct selectron pair production (left) and smuon pair production (right) followed by decay of both sleptons to the corresponding lepton and neutralino (color scale). The limits are shown as a function of the slepton mass and the mass difference between the slepton and the lightest neutralino for the three different simplified possibilities of only RR, only LL, and both RR and LL where it is assumed that the R and L masses are identical. The regions to the left of the lines denote the regions excluded for a branching fraction of 100%. Median expected limits for 100% branching fraction are delimited by the dashed lines.

Slepton pair production. Observed 95% CL upper limits on the cross section times branching fraction squared for direct selectron pair production (left) and smuon pair production (right) followed by decay of both sleptons to the corresponding lepton and neutralino (color scale). The limits are shown as a function of the slepton mass and the mass difference between the slepton and the lightest neutralino for the three different simplified possibilities of only RR, only LL, and both RR and LL where it is assumed that the R and L masses are identical. The regions to the left of the lines denote the regions excluded for a branching fraction of 100%. Median expected limits for 100% branching fraction are delimited by the dashed lines.

The observed and expected upper limits at 95% CLs on $\sigma (\mathrm{pp} \to \mathrm{H}+\text{X}) {\mathcal{B}} (\mathrm{H} \to \mathrm{a}_1\,\mathrm{a}_1)$, relative to $\sigma_\text{SM}$, as a function of $m_{\mathrm{a}_1}$ for Type III 2HD+S model, $\tan\beta = 2$.

The observed and expected upper limits at 95% CLs on $\sigma (\mathrm{pp} \to \mathrm{H}+\text{X}) {\mathcal{B}} (\mathrm{H} \to \mathrm{a}_1\,\mathrm{a}_1)$, relative to $\sigma_\text{SM}$, as a function of $m_{\mathrm{a}_1}$ for Type IV 2HD+S model, $\tan\beta = 0.5$.

The observed and expected 95% CL upper limits on $\sigma (\mathrm{pp} \to \mathrm{H}+\text{X}) {\mathcal{B}} (\mathrm{H} \to \mathrm{a}_1\,\mathrm{a}_1)$, relative to $\sigma_\text{SM}$, as a function of $\tan\beta$ for the Type II 2HD+S model for $m_{\mathrm{a}_1} = 5$ GeV.

The observed and expected 95% CL upper limits on $\sigma (\mathrm{pp} \to \mathrm{H}+\text{X}) {\mathcal{B}} (\mathrm{H} \to \mathrm{a}_1\,\mathrm{a}_1)$, relative to $\sigma_\text{SM}$, as a function of $\tan\beta$ for the Type II 2HD+S model for $m_{\mathrm{a}_1} = 8$ GeV.

The observed and expected 95% CL upper limits on $\sigma (\mathrm{pp} \to \mathrm{H}+\text{X}) {\mathcal{B}} (\mathrm{H} \to \mathrm{a}_1\,\mathrm{a}_1)$, relative to $\sigma_\text{SM}$, as a function of $\tan\beta$ for the Type II 2HD+S model for $m_{\mathrm{a}_1} = 12$ GeV.

The observed and expected 95% CL upper limits on $\sigma (\mathrm{pp} \to \mathrm{H}+\text{X}) {\mathcal{B}} (\mathrm{H} \to \mathrm{a}_1\,\mathrm{a}_1)$, relative to $\sigma_\text{SM}$, as a function of $\tan\beta$ for the Type II 2HD+S model for $m_{\mathrm{a}_1} = 15$ GeV.

The observed and expected 95% CL upper limits on $\sigma (\mathrm{pp} \to \mathrm{H}+\text{X}) {\mathcal{B}} (\mathrm{H} \to \mathrm{a}_1\,\mathrm{a}_1)$, relative to $\sigma_\text{SM}$, as a function of $\tan\beta$ for the Type III 2HD+S model for $m_{\mathrm{a}_1} = 5$ GeV.

The observed and expected 95% CL upper limits on $\sigma (\mathrm{pp} \to \mathrm{H}+\text{X}) {\mathcal{B}} (\mathrm{H} \to \mathrm{a}_1\,\mathrm{a}_1)$, relative to $\sigma_\text{SM}$, as a function of $\tan\beta$ for the Type III 2HD+S model for $m_{\mathrm{a}_1} = 8$ GeV.

The observed and expected 95% CL upper limits on $\sigma (\mathrm{pp} \to \mathrm{H}+\text{X}) {\mathcal{B}} (\mathrm{H} \to \mathrm{a}_1\,\mathrm{a}_1)$, relative to $\sigma_\text{SM}$, as a function of $\tan\beta$ for the Type III 2HD+S model for $m_{\mathrm{a}_1} = 12$ GeV.

The observed and expected 95% CL upper limits on $\sigma (\mathrm{pp} \to \mathrm{H}+\text{X}) {\mathcal{B}} (\mathrm{H} \to \mathrm{a}_1\,\mathrm{a}_1)$, relative to $\sigma_\text{SM}$, as a function of $\tan\beta$ for the Type III 2HD+S model for $m_{\mathrm{a}_1} = 15$ GeV.

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.

The 95% CL limits on the cross-section $\sigma(pp \rightarrow Z' \rightarrow \chi \chi$) times branching ratio B in fb with all statistical and systematic uncertainties, for the $R_{\text{inv}}=$0.6 signal points.

Data yield in the PFN signal region.

Yield of data in the ANTELOPE signal region in the binning used for BumpHunter fitting.

Observed profile likelihood as a function of $\sigma\times\mathcal{B}_{H\rightarrow WW^{*}}$ normalised by the SM expectation for the single-channel measurements

Two-dimensional likelihood scan of the measured values of $\sigma_{ZH}\times\mathcal{B}_{H\rightarrow WW^{*}}$ vs. $\sigma_{WH}\times\mathcal{B}_{H\rightarrow WW^{*}}$.

The observed values of the background normalisation factors for the combined 2-POI fit. The uncertainties correspond to the total of all statistical and systematic sources.

Measured cross sections times the $H\rightarrow WW^{*}$ branching ratio for the $p_T^V$ scheme.

Measured cross sections times the $H\rightarrow WW^{*}$ branching ratio for the STXS scheme. Note: the uncertainties of "-0" are due to the confidence interval reaching the minimum allowed value of the POI.

Correlation matrix of the POIs and normalisation factors for the 2-POI inclusive analysis

Correlation matrix of the parameters of interest, normalisation factors and nuissance parameters for the $p_T^V$ scheme

Correlation matrix of the parameters of interest, normalisation factors and nuissance parameters for the STXS scheme

Event data from Z-dominated SR