Unfolded distributions of $\Sigma p_{T}$ density in Z events, as a function of $p_{T}^{\mu\mu}$ in the towards ($\Delta\phi< 60^{\circ}$) region. Error bars represent the statistical and systematic uncertainties added in quadrature.

Unfolded distributions of $\Sigma p_{T}$ density in Z events, as a function of $p_{T}^{\mu\mu}$ in the transverse ($60^{\circ} <\Delta\phi< 120^{\circ}$) region. Error bars represent the statistical and systematic uncertainties added in quadrature.

Unfolded distributions of $\Sigma p_{T}$ density in Z events, as a function of $p_{T}^{\mu\mu}$ in the away ($\Delta\phi> 120^{\circ}$) region. Error bars represent the statistical and systematic uncertainties added in quadrature.

The measured $E_{T}^{miss}$ distribution in the Z/$\gamma ^{*}$($\rightarrow \mu \mu$)+jets control region, for the $E_{T}^{miss}$ > 250GeV inclusive selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit. The last bin of the distribution contains overflows.

The measured leading jet $p_{T}$ distribution in the Z/$\gamma ^{*}$($\rightarrow \mu \mu$)+jets control region, for the $E_{T}^{miss}$ > 250GeV inclusive selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit. The last bin of the distribution contains overflows.

The measured $E_{T}^{miss}$ distribution in the top control region, for the $E_{T}^{miss}$ > 250GeV inclusive selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit. The last bin of the distribution contains overflows.

The measured leading jet $p_{T}$ distribution in the top control region, for the $E_{T}^{miss}$ > 250GeV inclusive selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit. The last bin of the distribution contains overflows.

Measured distribution of the $E_{T}^{miss}$ for the $E_{T}^{miss}$ > 250GeV selection compared to the SM predictions. The latter are normalized with normalization factors as determined by the global fit that considers exclusive $E_{T}^{miss}$ regions. The last bin of the distribution contains overflows.

Measured distribution of the leading jet $p_{T}$ for the $E_{T}^{miss}$ > 250GeV selection compared to the SM predictions. The latter are normalized with normalization factors as determined by the global fit that considers exclusive $E_{T}^{miss}$ regions. The last bin of the distribution contains overflows.

Measured distribution of the leading jet $|\eta|$ for the $E_{T}^{miss}$ > 250GeV selection compared to the SM predictions. The latter are normalized with normalization factors as determined by the global fit that considers exclusive $E_{T}^{miss}$ regions. The last bin of the distribution contains overflows.

Measured distribution of the jet multiplicity for the $E_{T}^{miss}$ > 250GeV selection compared to the SM predictions. The latter are normalized with normalization factors as determined by the global fit that considers exclusive $E_{T}^{miss}$ regions. The last bin of the distribution contains overflows.

The expected $95\%$ CL exclusion limit for a simplified model of dark matter production involving an axial-vector operator, Dirac DM and couplings $g_{q} = 0.25$ and $g_{\chi} = 1$ as a function of the assumed mediator mass m$_{Z_{A}}$ and the dark matter mass m$_{\chi}$.

The measured $E_{T}^{miss}$ distribution in the W($\rightarrow \mu \nu$)+jets control region, for the $E_{T}^{miss}$ > 250GeV inclusive selection, compared to the background predictions. The latter include the global normalization factors extracted from the fit. The last bin of the distribution contains overflows.

The observed $95\%$ CL exclusion limit for a simplified model of dark matter production involving an axial-vector operator, Dirac DM and couplings $g_{q} = 0.25$ and $g_{\chi} = 1$ as a function of the assumed mediator mass m$_{Z_{A}}$ and the dark matter mass m$_{\chi}$.

The observed $90\%$ CL exclusion limit on the spin-dependent WIMP–proton scattering cross section in the context of the simplified model with axial-vector couplings, assuming minimal mediator width and the coupling values $g_{q} = 0.25$ and $g_{\chi} = 1$.

The expected $95\%$ CL exclusion limit for a simplified model of dark matter production involving a vector operator, Dirac DM and couplings $g_{q} = 0.25$ and $g_{\chi} = 1$ as a function of the assumed mediator mass m$_{Z_{V}}$ and the dark matter mass m$_{\chi}$.

The observed $95\%$ CL exclusion limit for a simplified model of dark matter production involving a vector operator, Dirac DM and couplings $g_{q} = 0.25$ and $g_{\chi} = 1$ as a function of the assumed mediator mass m$_{Z_{V}}$ and the dark matter mass m$_{\chi}$.

The expected and observed $95\%$ CL limits on the signal strength $\mu = \sigma^{95\% CL}/\sigma$ as a function of the mediator mass for a very light WIMP, in a model with spin-0 pseudoscalar mediator and $g_{q}=g_{\chi}=1.0$.

The expected and observed $95\%$ CL limits on the signal strength $\mu = \sigma^{95\% CL}/\sigma$ as a function of the WIMP mass for $m_{Z_{P}}=10$ GeV, in a model with spin-0 pseudoscalar mediator and $g_{q}=g_{\chi}=1.0$.

The expected exclusion contour at $95\%$ CL in the m$_{\eta}$–m$_{\chi}$ parameter plane for the coloured scalar mediator model, with minimal width and coupling set to $g=1$.

The observed exclusion contour at $95\%$ CL in the m$_{\eta}$–m$_{\chi}$ parameter plane for the coloured scalar mediator model, with minimal width and coupling set to $g=1$.

The expected excluded region at the $95\%$ CL in the ($\tilde{t}_{1}$,$\chi^{0}_{1}$) mass plane for the decay channel $\tilde{t}_{1} \rightarrow c + \chi^{0}_{1}$ (B = $100\%$).

The observed excluded region at the $95\%$ CL in the ($\tilde{t}_{1}$,$\chi^{0}_{1}$) mass plane for the decay channel $\tilde{t}_{1} \rightarrow c + \chi^{0}_{1}$ (B = $100\%$).

The expected excluded region at the $95\%$ CL in the ($\tilde{t}_{1}$,$\chi^{0}_{1}$) mass plane for the decay channel $\tilde{t}_{1} \rightarrow b + ff' + \chi^{0}_{1}$ (B = $100\%$).

The observed excluded region at the $95\%$ CL in the ($\tilde{t}_{1}$,$\chi^{0}_{1}$) mass plane for the decay channel $\tilde{t}_{1} \rightarrow b + ff' + \chi^{0}_{1}$ (B = $100\%$).

The expected exclusion plane at $95\%$ CL as a function of sbottom and neutralino masses for the decay channel $\tilde{b}_{1} \rightarrow b + \chi^{0}_{1}$ (B = $100\%$).

The observed exclusion plane at $95\%$ CL as a function of sbottom and neutralino masses for the decay channel $\tilde{b}_{1} \rightarrow b + \chi^{0}_{1}$ (B = $100\%$).

The expected exclusion region at $95\%$ CL as a function of squark mass and the squark-neutralino mass difference for $\tilde{q}_{1} → q + \chi^{0}_{1}$ (q =u,d,c,s).

The observed exclusion region at $95\%$ CL as a function of squark mass and the squark-neutralino mass difference for $\tilde{q}_{1} → q + \chi^{0}_{1}$ (q =u,d,c,s).

Expected and observed $95\%$ CL lower limits on the fundamental Planck scale in 4+n dimensions, M$_D$, as a function of the number of extra dimensions.

Expected and observed $95\%$ CL upper limit on the signal strength $\mu$ in the hypothesis of an axial-vector mediator, g$_{q}=0.25$, g$_{\chi}=1.0$ and minimal mediator width, as a function of the assumed mediator and DM masses.

Observed $90\%$ CL exclusion limit on the spin-dependent WIMP–neutron scattering cross section in the context of the simplified model with axial-vector couplings, assuming minimal mediator width and the coupling values $g_{q}=0.25$ and $g_{\chi}=1$.

Expected and observed $95\%$ CL upper limit on the signal strength $\mu$ in the hypothesis of a pseudoscalar mediator, $g_{q}=g_{\chi}=1.0$ and minimal mediator width, as a function of the assumed mediator and DM masses.

Statistical correlation between bins in data for the differential cross sections for the production of W⁺ bosons as a function of the inclusive jet multiplicity.

Statistical correlation between bins in data for the differential cross sections for the production of W^- bosons as a function of the inclusive jet multiplicity.

Differential cross section for the production of W bosons as a function of H_T for events with N_jets &ge; 1.

Statistical correlation between bins in data for the differential cross section for the production of W bosons as a function of H_T for events with N_jets &ge; 1.

Differential cross sections for the production of W⁺ bosons, W^- bosons and the W⁺/W^- cross section ratio as a function of the H_T for events with N_jets &ge; 1.

Statistical correlation between bins in data for the differential cross sections for the production of W⁺ bosons as a function of the H_T for events with N_jets &ge; 1.

Statistical correlation between bins in data for the differential cross sections for the production of W^- bosons as a function of the H_T for events with N_jets &ge; 1.

Differential cross section for the production of W bosons as a function of the W p_T for events with N_jets &ge; 1.

Statistical correlation between bins in data for the differential cross section for the production of W bosons as a function of the W p_T for events with N_jets &ge; 1.

Differential cross sections for the production of W⁺ bosons, W^- bosons and the W⁺/W^- cross section ratio as a function of the W p_T for events with N_jets &ge; 1.

Statistical correlation between bins in data for the differential cross sections for the production of W⁺ bosons as a function of the W p_T for events with N_jets &ge; 1.

Statistical correlation between bins in data for the differential cross sections for the production of W^- bosons as a function of the W p_T for events with N_jets &ge; 1.

Differential cross section for the production of W bosons as a function of the leading jet p_T for events with N_jets &ge; 1.

Statistical correlation between bins in data for the differential cross section for the production of W bosons as a function of the leading jet p_T for events with N_jets &ge; 1.

Differential cross sections for the production of W⁺ bosons, W^- bosons and the W⁺/W^- cross section ratio as a function of the leading jet p_T for events with N_jets &ge; 1.

Statistical correlation between bins in data for the differential cross sections for the production of W⁺ bosons as a function of the leading jet p_T for events with N_jets &ge; 1.

Statistical correlation between bins in data for the differential cross sections for the production of W^- bosons as a function of the leading jet p_T for events with N_jets &ge; 1.

Differential cross section for the production of W bosons as a function of the leading jet rapidity for events with N_jets &ge; 1.

Statistical correlation between bins in data for the differential cross section for the production of W bosons as a function of the leading jet rapidity for events with N_jets &ge; 1.

Differential cross sections for the production of W⁺ bosons, W^- bosons and the W⁺/W^- cross section ratio as a function of the leading jet rapidity for events with N_jets &ge; 1.

Statistical correlation between bins in data for the differential cross sections for the production of W⁺ bosons as a function of the leading jet rapidity for events with N_jets &ge; 1.

Statistical correlation between bins in data for the differential cross sections for the production of W^- bosons as a function of the leading jet rapidity for events with N_jets &ge; 1.

Differential cross section for the production of W bosons as a function of second leading jet p_T for events with N_jets &ge; 2.

Statistical correlation between bins in data for the differential cross section for the production of W bosons as a function of second leading jet p_T for events with N_jets &ge; 2.

Differential cross section for the production of W bosons as a function of second leading jet rapidity for events with N_jets &ge; 2.

Statistical correlation between bins in data for the differential cross section for the production of W bosons as a function of second leading jet rapidity for events with N_jets &ge; 2.

Differential cross section for the production of W bosons as a function of &Delta; R_jet1,jet2 for events with N_jets &ge; 2.

Statistical correlation between bins in data for the differential cross section for the production of W bosons as a function of &Delta; R_jet1,jet2 for events with N_jets &ge; 2.

Differential cross section for the production of W bosons as a function of dijet invariant mass for events with N_jets &ge; 2.

Statistical correlation between bins in data for the differential cross section for the production of W bosons as a function of dijet invariant mass for events with N_jets &ge; 2.

Cross section for the production of W bosons as a function of exclusive jet multiplicity.

Statistical correlation between bins in data for the cross section for the production of W bosons as a function of exclusive jet multiplicity.

Differential cross section for the production of W bosons as a function of the H_T for events with N_jets &ge; 2.

Statistical correlation between bins in data for the differential cross section for the production of W bosons as a function of the H_T for events with N_jets &ge; 2.

Differential cross sections for the production of W⁺ bosons, W^- bosons and the W⁺/W^- cross section ratio as a function of the H_T for events with N_jets &ge; 2.

Statistical correlation between bins in data for the differential cross sections for the production of W⁺ bosons as a function of the H_T for events with N_jets &ge; 2.

Statistical correlation between bins in data for the differential cross sections for the production of W^- bosons as a function of the H_T for events with N_jets &ge; 2.

Differential cross section for the production of W bosons as a function of the W p_T for events with N_jets &ge; 2.

Statistical correlation between bins in data for the differential cross section for the production of W bosons as a function of the W p_T for events with N_jets &ge; 2.

Differential cross sections for the production of W⁺ bosons, W^- bosons and the W⁺/W^- cross section ratio as a function of the W p_T for events with N_jets &ge; 2.

Statistical correlation between bins in data for the differential cross sections for the production of W⁺ bosons as a function of the W p_T for events with N_jets &ge; 2.

Statistical correlation between bins in data for the differential cross sections for the production of W^- bosons as a function of the W p_T for events with N_jets &ge; 2.

Differential cross section for the production of W bosons as a function of the leading jet p_T for events with N_jets &ge; 2.

Statistical correlation between bins in data for the differential cross section for the production of W bosons as a function of the leading jet p_T for events with N_jets &ge; 2.

Differential cross sections for the production of W⁺ bosons, W^- bosons and the W⁺/W^- cross section ratio as a function of the leading jet p_T for events with N_jets &ge; 2.

Statistical correlation between bins in data for the differential cross sections for the production of W⁺ bosons as a function of the leading jet p_T for events with N_jets &ge; 2.

Statistical correlation between bins in data for the differential cross sections for the production of W^- bosons as a function of the leading jet p_T for events with N_jets &ge; 2.

Differential cross section for the production of W bosons as a function of the electron &eta; for events with N_jets &ge; 0.

Statistical correlation between bins in data for the differential cross section for the production of W bosons as a function of the electron &eta; for events with N_jets &ge; 0.

Differential cross sections for the production of W⁺ bosons, W^- bosons and the W⁺/W^- cross section ratio as a function of the electron &eta; for events with N_jets &ge; 0.

Statistical correlation between bins in data for the differential cross sections for the production of W⁺ bosons as a function of the electron &eta; for events with N_jets &ge; 0.

Statistical correlation between bins in data for the differential cross sections for the production of W^- bosons as a function of the electron &eta; for events with N_jets &ge; 0.

Differential cross section for the production of W bosons as a function of the electron &eta; for events with N_jets &ge; 1.

Statistical correlation between bins in data for the differential cross section for the production of W bosons as a function of the electron &eta; for events with N_jets &ge; 1.

Differential cross sections for the production of W⁺ bosons, W^- bosons and the W⁺/W^- cross section ratio as a function of the electron &eta; for events with N_jets &ge; 1.

Statistical correlation between bins in data for the differential cross sections for the production of W⁺ bosons as a function of the electron &eta; for events with N_jets &ge; 1.

Statistical correlation between bins in data for the differential cross sections for the production of W^- bosons as a function of the electron &eta; for events with N_jets &ge; 1.

List of experimentally considered systematic uncertainties for the W+jets cross section measurement

Non-perturbative corrections for the cross section for the production of W bosons for different inclusive jet multiplicities.

Non-perturbative corrections for the differential cross sections for the production of W⁺ bosons and W^- bosons as a function of the inclusive jet multiplicity.

Non-perturbative corrections for the differential cross section for the production of W bosons as a function of H_T for events with N_jets &ge; 1.

Non-perturbative corrections for the differential cross sections for the production of W⁺ bosons and W^- bosons as a function of the H_T for events with N_jets &ge; 1.

Non-perturbative corrections for the differential cross section for the production of W bosons as a function of the W p_T for events with N_jets &ge; 1.

Non-perturbative corrections for the differential cross sections for the production of W⁺ bosons and W^- bosons as a function of the W p_T for events with N_jets &ge; 1.

Non-perturbative corrections for the differential cross section for the production of W bosons as a function of the leading jet p_T for events with N_jets &ge; 1.

Non-perturbative corrections for the differential cross sections for the production of W⁺ bosons and W^- bosons as a function of the leading jet p_T for events with N_jets &ge; 1.

Non-perturbative corrections for the differential cross section for the production of W bosons as a function of the leading jet rapidity for events with N_jets &ge; 1.

Non-perturbative corrections for the differential cross sections for the production of W⁺ bosons and W^- bosons as a function of the leading jet rapidity for events with N_jets &ge; 1.

Non-perturbative corrections for the differential cross section for the production of W bosons as a function of second leading jet p_T for events with N_jets &ge; 2.

Non-perturbative corrections for the differential cross section for the production of W bosons as a function of second leading jet rapidity for events with N_jets &ge; 2.

Non-perturbative corrections for the differential cross section for the production of W bosons as a function of &Delta; R_jet1,jet2 for events with N_jets &ge; 2.

Non-perturbative corrections for the differential cross section for the production of W bosons as a function of dijet invariant mass for events with N_jets &ge; 2.

Non-perturbative corrections for the cross section for the production of W bosons as a function of exclusive jet multiplicity.

Non-perturbative corrections for the differential cross section for the production of W bosons as a function of the H_T for events with N_jets &ge; 2.

Non-perturbative corrections for the differential cross sections for the production of W⁺ bosons and W^- bosons as a function of the H_T for events with N_jets &ge; 2.

Non-perturbative corrections for the differential cross section for the production of W bosons as a function of the W p_T for events with N_jets &ge; 2.

Non-perturbative corrections for the differential cross sections for the production of W⁺ bosons and W^- bosons as a function of the W p_T for events with N_jets &ge; 2.

Non-perturbative corrections for the differential cross section for the production of W bosons as a function of the leading jet p_T for events with N_jets &ge; 2.

Non-perturbative corrections for the differential cross sections for the production of W⁺ bosons and W^- bosons as a function of the leading jet p_T for events with N_jets &ge; 2.

NNLO/NLO k-factors determined with NNLO Njetti for the differential cross sections for the production of W⁺ bosons and W^- bosons as a function of the H_T for events with N_jets &ge; 1. These numbers were obtained with code described in Phys. Rev. Lett. 115 (2015) 062002 [arXiv:1504.02131].

NNLO/NLO k-factors determined with NNLO Njetti for the differential cross sections for the production of W⁺ bosons and W^- bosons as a function of the W p_T for events with N_jets &ge; 1. These numbers were obtained with code described in Phys. Rev. Lett. 115 (2015) 062002 [arXiv:1504.02131].

NNLO/NLO k-factors determined with NNLO Njetti for the differential cross sections for the production of W⁺ bosons and W^- bosons as a function of the leading jet p_T for events with N_jets &ge; 1. These numbers were obtained with code described in Phys. Rev. Lett. 115 (2015) 062002 [arXiv:1504.02131].

NNLO/NLO k-factors determined with NNLO Njetti for the differential cross sections for the production of W⁺ bosons and W^- bosons as a function of the leading jet rapidity for events with N_jets &ge; 1. These numbers were obtained with code described in Phys. Rev. Lett. 115 (2015) 062002 [arXiv:1504.02131].

The observed distribution of the dijet invariant mass in the $\ell$+jets decay channel in the 0 b tagged jet category. The vertical bars on the data points represent the statistical uncertainties. The integral in the $\rm{e}$($\mu$)+jets channel (including overflow events) is 1375 (2406).

The observed distribution of the dijet invariant mass in the $\ell$+jets decay channel in the 1 b tagged jet category. The vertical bars on the data points represent the statistical uncertainties. The integral in the $\rm{e}$($\mu$)+jets channel (including overflow events) is 61 (129).

The observed distribution of the dijet invariant mass in the $\ell$+jets decay channel in the $\geq$2 b tagged jet category. The vertical bars on the data points represent the statistical uncertainties. The integral in the $\rm{e}$($\mu$)+jets channel (including overflow events) is 19 (33).

The observed distribution of the the minimum angular distance $\Delta R$ between all pairs of jets $j$ and $j'$ in the $\ell$+jets decay channel in the 0 b tagged jet category. The vertical bars on the data points represent the statistical uncertainties. The integral in the $\rm{e}$($\mu$)+jets channel (including overflow events) is 1375 (2406).

The observed distribution of the the minimum angular distance $\Delta R$ between all pairs of jets $j$ and $j'$ in the $\ell$+jets decay channel in the 1 b tagged jet category. The vertical bars on the data points represent the statistical uncertainties. The integral in the $\rm{e}$($\mu$)+jets channel (including overflow events) is 61 (129).

The observed distribution of the the minimum angular distance $\Delta R$ between all pairs of jets $j$ and $j'$ in the $\ell$+jets decay channel in the $\geq$2 b tagged jet category. The vertical bars on the data points represent the statistical uncertainties. The integral in the $\rm{e}$($\mu$)+jets channel (including overflow events) is 19 (33).

The observed distribution of the jet multiplicity in the $\rm{e}^\pm \mu^\mp$ decay channel for events passing the dilepton criteria. The error bars indicate the statistical uncertainties. The last bin of the distributions contains the overflow events.

The observed distribution of the scalar $p_{\rm{T}}$ sum of all jets in the $\rm{e}^\pm \mu^\mp$ decay channel for events passing the dilepton criteria. The error bars indicate the statistical uncertainties. The last bin of the distributions contains the overflow events.

The observed distribution of the invariant mass of the lepton pair in the $\rm{e}^\pm \mu^\mp$ decay channel after requiring at least two jets. The error bars indicate the statistical uncertainties. The integral is 24. The last bin of the distributions contains the overflow events.

The observed distribution of the $p_{\rm{T}}$ of the lepton pair in the $\rm{e}^\pm \mu^\mp$ decay channel after requiring at least two jets. The error bars indicate the statistical uncertainties. The integral is 24. The last bin of the distributions contains the overflow events.

The observed distribution of the $p_{\rm{T}}^{\rm{miss}}$ in the $\mu^\pm \mu^\mp$ decay channel for events passing the dilepton criteria and Z boson veto. The error bars indicate the statistical uncertainties. The last bin of the distributions contains the overflow events.

The observed distribution of the invariant mass of the lepton pair in the $\mu^\pm \mu^\mp$ decay channel for events passing the dilepton criteria, Z boson veto, and $p_{\rm{T}}^{\rm{miss}}>35$ GeV requirement. The error bars indicate the statistical uncertainties. The last bin of the distributions contains the overflow events.

rapidity bin 1.5 < |Y| < 2.0 anti-kt R=0.4

rapidity bin 2.0 < |Y| < 2.5 anti-kt R=0.4

rapidity bin 2.5 < |Y| < 3.0 anti-kt R=0.4

rapidity bin 0 < y* < 0.5 anti-kt R=0.4

rapidity bin 0.5 < y* < 1.0 anti-kt R=0.4

rapidity bin 1.0 < y* < 1.5 anti-kt R=0.4

rapidity bin 1.5 < y* < 2.0 anti-kt R=0.4

rapidity bin 2.0 < y* < 2.5 anti-kt R=0.4

rapidity bin 2.5 < y* < 3.0 anti-kt R=0.4

Differential cross section dsigma(Z+c)/dpTZ times branching ratio and differential cross sections ratio dsigma(Z+c)/dpTZ / dsigma(Z+b)/dpTZ for three ranges of the transverse momentum of the Z boson and combining the dielectron and dimuon Z0 decay channels

Differential cross section dsigma(Z+c)/dpTjet times branching ratio and differential cross sections ratio dsigma(Z+c)/dpTZ / dsigma(Z+b)/dpTjet for three ranges of the transverse momentum of the jet and combining the dielectron and dimuon Z0 decay channels

Transverse momentum distribution of the selected muon-inside-a-jet for events with an identified muon among the jet constituents in the dielectron channel

Transverse momentum distribution of the selected muon-inside-a-jet for events with an identified muon among the jet constituents in the dimuon channel

The invariant mass distribution of three-prong secondary vertices for events selected in the D± mode, in the dielectron channel

The invariant mass distribution of three-prong secondary vertices for events selected in the D ± mode, in the dimuon channel

The invariant mass distribution of the three-track system composed of a two-prong secondary vertex and a primary particle for events selected in the D∗(2010)± mode, in the dielectron channel

The invariant mass distribution of the three-track system composed of a two-prong secondary vertex and a primary particle for events selected in the D ∗ (2010) ± mode, in the dimuon channel

Transverse momentum distribution of the c-tagged jet normalized to unity, in simulated W+c and Z+c samples and in W+c data events.

Number of reconstructed secondary vertices normalized to unity, in simulated W+c and Z+c samples and in W+c data events.

Distribution of the corrected secondary-vertex mass normalized to unity, in simulated W + c and Z + c samples, and in W + c data events.

Distribution of the JP discriminant for the D± mode normalized to unity, in simulated W + c and Z + c samples, and in W + c data events.

Distribution of the JP discriminant for the D∗(2010) ± mode normalized to unity, in simulated W + c and Z + c samples, and in W + c data events.

Distribution of the corrected secondary-vertex mass normalized to unity from simulated Z + b and eμ-tt data events

Corrected secondary-vertex mass distributions, after background subtraction in the dielectron channel for events selected in the semileptonic mode.

Corrected secondary-vertex mass distributions, after background subtraction in the dimuon channel for events selected in the semileptonic mode.

Background-subtracted distributions of the JP discriminant in the dielectron channel for Z + jets events with a D± → K∓ π± π± candidate.

Background-subtracted distributions of the JP discriminant in the dimuon channel for Z + jets events with a D± → K∓ π± π± candidate.

Background-subtracted distributions of the JP discriminant in the dielectron channel for Z +jets events with a D∗(2010)± → D0 π± → K∓ π± π± candidate.

Background-subtracted distributions of the JP discriminant in the dimuon channel for Z +jets events with a D∗(2010)± → D0 π± → K∓ π± π± candidate.

ontributions to the systematic uncertainty in the measured Z+c cross section and in the (Z+c)/(Z+b) cross section ratio.

Differential Z + c cross section as a function of the transverse momentum of the Z boson

(Z +c)/(Z +b) cross section ratio as a function of the transverse momentum of teh Z boson

Differential Z + c cross section as a function of the transverse momentum of the jet

(Z +c)/(Z +b) cross section ratio as a function of the transverse momentum of the jet.

Figure 3 (upper). Distribution of $M_{T2}(ll)$ in a control region enriched in $t\bar{t}$ events defined by requiring at least 2 jets, at least one of them b-tagged, and missing transverse momentum below 80 GeV. For the plot, simulated yields are normalized to data using the yields in the $M_{T2}(ll)$ < 100 GeV region.

Figure 4 (upper). Expected and observed yields in the five $t\bar{t}Z$ control regions, which are defined by different requirements on the number of reconstructed jets and b jets, before the fit.

Figure 4 (lower). Expected and observed yields in the five $t\bar{t}Z$ control regions, which are defined by different requirements on the number of reconstructed jets and b jets, after the fit.

Figure 5 (left). Distribution of $M_{T2}(ll)$ of SF events falling within in the Z boson mass window with at least two jets and no b jet, $p_{T}^{miss} > 80$ GeV and $M_{T2}(ll) > 100$ GeV.

Figure 5 (center). Distribution of $M_{T2}(blbl)$ of SF events falling within in the Z boson mass window with at least two jets and no b jet, $p_{T}^{miss} > 80$ GeV and $M_{T2}(ll) > 100$ GeV.

Figure 5 (right). Distribution of $p_{T}^{miss}$ of SF events falling within in the Z boson mass window with at least two jets and no b jet, $p_{T}^{miss} > 80$ GeV and $M_{T2}(ll) > 100$ GeV.

Figure 6. Event yields in the 13 Drell-Yan and multiboson control regions for events with SF leptons falling within the Z boson mass window and no b jets, after renormalizing with the scale factors obtained from the fit procedure described in the text.

Figure 7 (left). Distributions of $M_{T2}(ll)$ for observed events in the dimuon channel compared to the predicted SM background.

Figure 7 (center). Distributions of $M_{T2}(ll)$ for observed events in the dielectron channel compared to the predicted SM background.

Figure 7 (right). Distributions of $M_{T2}(ll)$ for observed events in the electron-muon channel compared to the predicted SM background.

Figure 8 (left). Distributions of $M_{T2}(ll)$ for all lepton flavors for the defined selection.

Figure 8 (center). Distributions of $M_{T2}(blbl)$ for all lepton flavors for the defined selection, after requiring $M_{T2}(ll) > 100$ GeV.

Figure 8 (right). Distributions of $p_{T}^{miss}$ for all lepton flavors for the defined selection, after requiring $M_{T2}(ll) > 100$ GeV.

Figure 9 (upper). Predicted backgrounds and observed yields in the dielectron and dimuon search regions. The yields are post fit and correspond to the numbers quoted in Tab. 4 of the publication.

Figure 9 (lower). Predicted backgrounds and observed yields in the electron-muon search regions. The yields are post fit and correspond to the numbers quoted in Tab. 4 of the publication.

Figure 10. Predicted backgrounds and observed yields in all search regions combined. The yields are post fit and correspond to the numbers quoted in Tab. 4 of the publication.

Figure 11 (upper). Observed limits for the T2tt model in the top squark - LSP mass plane. The numbers indicate the 95% CL upper limit on the cross section times the square of the branching fraction at each point of the plane.

Figure 11 (upper). Observed exclusion region at 95% CL assuming 100% branching fraction.

Figure 11 (upper). Expected exclusion region at 95% CL assuming 100% branching fraction.

Figure 11 (lower). Observed limits for the T2bW model in the top squark - LSP mass plane. The numbers indicate the 95% CL upper limit on the cross section times the square of the branching fraction at each point of the plane.

Figure 11 (lower). Observed exclusion region at 95% CL assuming 100% branching fraction.

Figure 11 (lower). Expected exclusion region at 95% CL assuming 100% branching fraction.

Figure 12 (upper left). Observed limits for the T8bbllnunu model in the top squark - LSP mass plane. The numbers indicate the 95% CL upper limit on the cross section times the square of the branching fraction at each point of the plane. $m_{\tilde{\chi}^{\pm}_{1}} = 0.5 (m_{\tilde{t}_{1}} + m_{\tilde{\chi}^{0}_{1}})$, $m_{\tilde{l}} = 0.05 (m_{\tilde{\chi}^{\pm}_{1}} - m_{\tilde{\chi}^{0}_{1}}) + m_{\tilde{\chi}^{0}_{1}}$

Figure 12 (upper left). Observed exclusion region at 95% CL assuming 100% branching fraction. $m_{\tilde{\chi}^{\pm}_{1}} = 0.5 (m_{\tilde{t}_{1}} + m_{\tilde{\chi}^{0}_{1}})$, $m_{\tilde{l}} = 0.05 (m_{\tilde{\chi}^{\pm}_{1}} - m_{\tilde{\chi}^{0}_{1}}) + m_{\tilde{\chi}^{0}_{1}}$

Figure 12 (upper left). Expected exclusion region at 95% CL assuming 100% branching fraction. $m_{\tilde{\chi}^{\pm}_{1}} = 0.5 (m_{\tilde{t}_{1}} + m_{\tilde{\chi}^{0}_{1}})$, $m_{\tilde{l}} = 0.05 (m_{\tilde{\chi}^{\pm}_{1}} - m_{\tilde{\chi}^{0}_{1}}) + m_{\tilde{\chi}^{0}_{1}}$

Figure 12 (upper left). Observed limits for the T8bbllnunu model in the top squark - LSP mass plane. The numbers indicate the 95% CL upper limit on the cross section times the square of the branching fraction at each point of the plane. $m_{\tilde{\chi}^{\pm}_{1}} = 0.5 (m_{\tilde{t}_{1}} + m_{\tilde{\chi}^{0}_{1}})$, $m_{\tilde{l}} = 0.5 (m_{\tilde{\chi}^{\pm}_{1}} - m_{\tilde{\chi}^{0}_{1}}) + m_{\tilde{\chi}^{0}_{1}}$

Figure 12 (upper left). Observed exclusion region at 95% CL assuming 100% branching fraction. $m_{\tilde{\chi}^{\pm}_{1}} = 0.5 (m_{\tilde{t}_{1}} + m_{\tilde{\chi}^{0}_{1}})$, $m_{\tilde{l}} = 0.5 (m_{\tilde{\chi}^{\pm}_{1}} - m_{\tilde{\chi}^{0}_{1}}) + m_{\tilde{\chi}^{0}_{1}}$

Figure 12 (upper left). Expected exclusion region at 95% CL assuming 100% branching fraction. $m_{\tilde{\chi}^{\pm}_{1}} = 0.5 (m_{\tilde{t}_{1}} + m_{\tilde{\chi}^{0}_{1}})$, $m_{\tilde{l}} = 0.5 (m_{\tilde{\chi}^{\pm}_{1}} - m_{\tilde{\chi}^{0}_{1}}) + m_{\tilde{\chi}^{0}_{1}}$

Figure 12 (lower). Observed limits for the T8bbllnunu model in the top squark - LSP mass plane. The numbers indicate the 95% CL upper limit on the cross section times the square of the branching fraction at each point of the plane. $m_{\tilde{\chi}^{\pm}_{1}} = 0.5 (m_{\tilde{t}_{1}} + m_{\tilde{\chi}^{0}_{1}})$, $m_{\tilde{l}} = 0.95 (m_{\tilde{\chi}^{\pm}_{1}} - m_{\tilde{\chi}^{0}_{1}}) + m_{\tilde{\chi}^{0}_{1}}$

Figure 12 (lower). Observed exclusion region at 95% CL assuming 100% branching fraction. $m_{\tilde{\chi}^{\pm}_{1}} = 0.5 (m_{\tilde{t}_{1}} + m_{\tilde{\chi}^{0}_{1}})$, $m_{\tilde{l}} = 0.95 (m_{\tilde{\chi}^{\pm}_{1}} - m_{\tilde{\chi}^{0}_{1}}) + m_{\tilde{\chi}^{0}_{1}}$

Figure 12 (lower). Expected exclusion region at 95% CL assuming 100% branching fraction. $m_{\tilde{\chi}^{\pm}_{1}} = 0.5 (m_{\tilde{t}_{1}} + m_{\tilde{\chi}^{0}_{1}})$, $m_{\tilde{l}} = 0.95 (m_{\tilde{\chi}^{\pm}_{1}} - m_{\tilde{\chi}^{0}_{1}}) + m_{\tilde{\chi}^{0}_{1}}$

Figure 13. Observed limits for the T2tt model in the top squark - LSP mass plane combining the dilepton final state with single-lepton and the all-hadronic final states. The numbers indicate the 95% CL upper limit on the cross section times the square of the branching fraction at each point of the plane.

Figure 13. Observed exclusion region at 95% CL assuming 100% branching fraction.

Figure 13. Expected exclusion region at 95% CL assuming 100% branching fraction.

Figure 14. Observed limits for the T2bW model in the top squark - LSP mass plane combining the dilepton final state with single-lepton and the all-hadronic final states. The numbers indicate the 95% CL upper limit on the cross section times the square of the branching fraction at each point of the plane.

Figure 14. Observed exclusion region at 95% CL assuming 100% branching fraction.

Figure 14. Expected exclusion region at 95% CL assuming 100% branching fraction.

Additional Figure 3. Observed limits for the T2tt model in the top squark - LSP mass plane using aggregate signal regions. The numbers indicate the 95% CL upper limit on the cross section times the square of the branching fraction at each point of the plane.

Additional Figure 3. Observed exclusion region at 95% CL assuming 100% branching fraction.

Additional Figure 3. Expected exclusion region at 95% CL assuming 100% branching fraction.

Table 5. Ratios of $\mu=\frac{\sigma$}{\sigma_{theory}}$ of the 95% CL expected and observed limits to simplified model expectations for different DM particle masses and mediator masses.

Table 6. Expected and observed event yields, summed over all lepton flavors, in aggregate signal regions.

Table 7 (upper). Covariance matrix for the background prediction in aggregate signal regions.

Table 7 (lower). Correlation matrix for the background prediction in aggregate signal regions.

Additional Figure 1. Correlation matrix of the 13 signal regions, split into same-flavor (SF) and different-flavor (DF) regions.

Additional Figure 2. Covariance matrix of the 13 signal regions, split into same-flavor (SF) and different-flavor (DF) regions.

Limit on the signal strength of the DM signal in a simplified model with a scalar mediator.

Limit on the signal strength of the DM signal in a simplified model with a pseudoscalar mediator.

Expected and observed 95% CL upper limits on the product of the production cross section and the branching fraction, $\sigma$(qq -> ZH) B(H -> inv.), as a function of the SM-like Higgs boson mass. The limits consider only quark-induced Higgs boson production.

The 95% CL upper limits on the Wilson coefficient $\lambda(1\mathrm{TeV} / \Lambda_{U})^{d_{U}-1}$ of the unparticle-quark coupling operator.

The 95% CL lower limits on the reduced Planck mass $M_{D}$ as a function of the number of extra dimensions $d$.

The covariance matrix of the bin contents of the background fit. Because the matrix is symmetric by definiton, only the upper half is provided here for presentational purposes. The full matrix can be obtained by mirroring the given half along the diagonal.

Fit results in SRt1, SRt2 and SRt3 for an integrated luminosity of $36.1 fb^{-1}$. The background normalisation parameters are obtained from the background-only fit in the CRs and are applied to the SRs. Small backgrounds are indicated as Others. Benchmark signal models yields are given for each SR. The uncertainties on the yields include all systematic uncertainties.

95% CL upper limits on the visible cross-section ($\langle\epsilon\mathcal{A}\sigma\rangle^{\rm obs}_{95}$) and on the number of BSM events ($S^{\rm obs}_{95}$ ). The third column ($S^{\rm exp}_{95}$) shows the 95% CL upper limit on the number of signal events, given the expected number (and $\pm 1\sigma$ excursions on the expectation) of background events. The last column indicates the discovery $p$-value ($p(s = 0)$) and Z (the number of equivalent Gaussian standard deviations).

Comparison of the data with the post-fit SM prediction of the $E_{\mathrm T}^{\mathrm{miss}}$ distribution in SRb1. The last bins include overflows, where applicable. All signal region requirements except the one on the distribution shown are applied. The signal region requirement on the distribution shown is indicated by an arrow. The bottom panel shows the ratio of the data over the prediction. The band includes all systematic uncertainties.

Comparison of the data with the post-fit SM prediction of the $\cos{\theta}^*_{bb}$ distribution in SRb2. The last bins include overflows, where applicable. All signal region requirements except the one on the distribution shown are applied. The signal region requirement on the distribution shown is indicated by an arrow. The bottom panel shows the ratio of the data over the prediction. The band includes all systematic uncertainties.

Comparison of the data with the post-fit SM prediction of the $m_{\mathrm T}^{\mathrm{b,min}}$ distribution in SRt1. The last bins include overflows, where applicable. All signal region requirements except the one on the distribution shown are applied. The signal region requirement on the distribution shown is indicated by an arrow. The bottom panel shows the ratio of the data over the prediction. The band includes all systematic uncertainties.

Comparison of the data with the post-fit SM prediction of the $E_{\mathrm T}^{\mathrm{miss,sig}}$ distribution in SRt2. The last bins include overflows, where applicable. All signal region requirements except the one on the distribution shown are applied. The signal region requirement on the distribution shown is indicated by an arrow. The bottom panel shows the ratio of the data over the prediction. The band includes all systematic uncertainties.

Comparison of the data with the post-fit SM prediction of the $\xi^{+}_{\ell\ell}$ distribution in SRt3. The last bins include overflows, where applicable. All signal region requirements except the one on the distribution shown are applied. The signal region requirement on the distribution shown is indicated by an arrow. The bottom panel shows the ratio of the data over the prediction. The band includes all systematic uncertainties.

Expected exclusion limits for colour-neutral scalar model in the SRb2 as a function of the mediator mass for a DM mass of 1 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$.

Observed exclusion limits for colour-neutral scalar model in SRb2 as a function of the mediator mass for a DM mass of 1 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$.

Expected exclusion limits for colour-neutral scalar model in the SRt1/SRt2 as a function of the mediator mass for a DM mass of 1 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$. To derive the results for the fully hadronic $t\bar{t}$ final state the region among SRt1 and SRt2 providing the best expected sensitivity is used.

Observed exclusion limits for colour-neutral scalar model in SRt1/SRt2 as a function of the mediator mass for a DM mass of 1 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$. To derive the results for the fully hadronic $t\bar{t}$ final state the region among SRt1 and SRt2 providing the best expected sensitivity is used.

Expected exclusion limits for colour-neutral scalar model in SRt3 as a function of the mediator mass for a DM mass of 1 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$.

Observed exclusion limits for colour-neutral scalar model in SRt3 as a function of the mediator mass for a DM mass of 1 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$.

Expected exclusion limits for colour-neutral pseudo-scalar model in SRb2 as a function of the mediator mass for a DM mass of 1 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$.

Observed exclusion limits for colour-neutral pseudo-scalar model in SRb2 as a function of the mediator mass for a DM mass of 1 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$.

Expected exclusion limits for colour-neutral pseudo-scalar model in SRt1/SRt2 as a function of the mediator mass for a DM mass of 1 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$. To derive the results for the fully hadronic $t\bar{t}$ final state the region among SRt1 and SRt2 providing the best expected sensitivity is used.

Observed exclusion limits for colour-neutral pseudo-scalar model in SRt1/SRt2 as a function of the mediator mass for a DM mass of 1 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$. To derive the results for the fully hadronic $t\bar{t}$ final state the region among SRt1 and SRt2 providing the best expected sensitivity is used.

Expected limits for colour-neutral pseudo-scalar model in SRt3 as a function of the mediator mass for a DM mass of 1 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$.

Observed exclusion limits for colour-neutral pseudo-scalar model in SRt3 as a function of the mediator mass for a DM mass of 1 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$.

Expected exclusion limits for colour-neutral scalar model in SRt1/SRt2 as a function of the DM mass for a mediator mass of 10 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$. To derive the results for the fully hadronic $t\bar{t}$ final state the region among SRt1 and SRt2 providing the best expected sensitivity is used.

Observed exclusion limits for colour-neutral scalar model in SRt1/SRt2 as a function of the DM mass for a mediator mass of 10 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$. To derive the results for the fully hadronic $t\bar{t}$ final state the region among SRt1 and SRt2 providing the best expected sensitivity is used.

Expected exclusion limits for colour-neutral scalar model in SRt3 as a function of the DM mass for a mediator mass of 10 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$.

Observed exclusion limits for colour-neutral scalar model in SRt3 as a function of the DM mass for a mediator mass of 10 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$.

Expected exclusion limits for colour-neutral pseudo-scalar model in SRt1/SRt2 as a function of the DM mass for a mediator mass of 10 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$. To derive the results for the fully hadronic $t\bar{t}$ final state the region among SRt1 and SRt2 providing the best expected sensitivity is used.

Observed exclusion limits for colour-neutral pseudo-scalar model in SRt1/SRt2 as a function of the DM mass for a mediator mass of 10 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$. To derive the results for the fully hadronic $t\bar{t}$ final state the region among SRt1 and SRt2 providing the best expected sensitivity is used.

Expected exclusion limits for colour-neutral pseudo-scalar model in SRt3 as a function of the DM mass for a mediator mass of 10 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$.

Observed xclusion limits for colour-neutral pseudo-scalar model in SRt3 as a function of the DM mass for a mediator mass of 10 GeV. The limits are calculated at 95% C.L. and are expressed in terms of the ratio of the excluded cross-section over the nominal cross-section for a coupling assumption of $g = g_\chi = g_\nu = 1$.

Comparison of the 90% CL limits on the spin-independent DM–nucleon cross-section versus mediator mass between these results and the direct-detection experiments, in the context of the colour-neutral simplified model with scalar mediator. The black line indicates the exclusion contour derived from the observed limits of SRt3. Values inside the contour are excluded.

Expected exclusion limits for colour-charged scalar mediators ($b$-FDM) as a function of the mediator and DM masses for $36.1fb^{-1}$ of data. The limits are calculated at 95% C.L for a coupling assumption $\lambda_b$ yielding the measured relic density.

Exclusion observed limits for colour-charged scalar mediators ($b-$FDM) as a function of the mediator and DM masses for $36.1fb^{-1}$ of data. The limits are calculated at 95% C.L for a coupling assumption $\lambda_b$ yielding the measured relic density.

Comparison of the data with the post-fit Monte Carlo prediction of sub leading jet $p_{T}$ in SRb1. The last bin includes overflows. The bottom panel shows the ratio of the data over the Monte Carlo prediction. The band includes all systematic uncertainties.

Comparison of the data with the post-fit Monte Carlo prediction of $H_{\mathrm T}^{ratio}$ distribution in SRb2. The last bin includes overflows. The bottom panel shows the ratio of the data over the Monte Carlo prediction. The band includes all systematic uncertainties.

Comparison of the data with the post-fit Monte Carlo prediction of $\Delta R_{bb}$ in SRt1. The last bin includes overflows. The bottom panel shows the ratio of the data over the Monte Carlo prediction. The band includes all systematic uncertainties.

Comparison of the data with the post-fit Monte Carlo prediction of $M_{\mathrm T}^{b,min}$ in SRt2. The last bin includes overflows. The bottom panel shows the ratio of the data over the Monte Carlo prediction. The band includes all systematic uncertainties.

Comparison of the data with the post-fit Monte Carlo prediction of $\Delta \phi_{boost}$ in SRt3. The last bin includes overflows. The bottom panel shows the ratio of the data over the Monte Carlo prediction. The band includes all systematic uncertainties.

Acceptance of the SRb1 selection of the $b$-FDM model signal samples

Efficiency of the SRb1 selection of the $b$-FDM model signal samples

Acceptance of the SRb2 selection of the colour-neutral simplified model samples with scalar mediator decaying into dark matter pairs produced in association with bottom-quark pairs.

Efficiency of the SRb2 selection to the colour-neutral simplified model samples with scalar mediator decaying into dark matter pairs produced in association with bottom-quark pairs

Efficiency of the SRb2 selection to the colour-neutral simplified model samples with pseudo-scalar mediator decaying into dark matter pairs produced in association with bottom-quark pairs

Acceptance of SRb2 selection to the colour-neutral simplified model samples with pseudo-scalar mediator decaying into dark matter pairs produced in association with bottom-quark pairs

Acceptance of the SRt2 selections to the color-neutral simplified model samples with scalar mediator decaying into dark matter pairs produced in association with top quarks.

Acceptance of the SRt1 selections to the color-neutral simplified model samples with scalar mediator decaying into dark matter pairs produced in association with top quarks.

Efficiency of the SRt2 selections to the color-neutral simplified model samples with scalar mediator decaying into dark matter pairs produced in association with top quarks.

Efficiency of the SRt1 selections to the color-neutral simplified model samples with scalar mediator decaying into dark matter pairs produced in association with top quarks.

Acceptance of the SRt2 selections to the color-neutral simplified model samples with pseudo-scalar mediator decaying into dark matter pairs produced in association with top quarks.

Acceptance of the SRt1 selections to the color-neutral simplified model samples with pseudo-scalar mediator decaying into dark matter pairs produced in association with top quarks.

Efficiency of the SRt2 selections to the color-neutral simplified model samples with pseudo-scalar mediator decaying into dark matter pairs produced in association with top quarks.

Efficiency of the SRt1 selections to the color-neutral simplified model samples with pseudo-scalar mediator decaying into dark matter pairs produced in association with top quarks.

Acceptance of the SRt3 selections to the color-neutral simplified model samples with scalar mediator decaying into dark matter pairs produced in association with top quarks.

Efficiency of the SRt3 selections to the color-neutral simplified model samples with scalar mediator decaying into dark matter pairs produced in association with top quarks.

Acceptance of the SRt3 selections to the color-neutral simplified model samples with pseudo-scalar mediator decaying into dark matter pairs produced in association with top quarks.

Efficiency of the SRt3 selections to the color-neutral simplified model samples with pseudo-scalar mediator decaying into dark matter pairs produced in association with top quarks.

Number of signal events selected at different stages of the SRb1 selections for the $b$-FDM benchmark model $m(\phi_b,\chi)=(1000,35)$GeV

Number of signal events selected at different stages of the SRb2 selections for the $b\bar{b} +\phi$ benchmark model $m(\phi,\chi)=(20,1)$GeV, $g=1$.

Number of signal events selected at different stages of the SRt1 selections for the $t \bar{t} +\phi$ benchmark model $m(\phi,\chi)=(20,1)$GeV, $g=1$.

Number of signal events selected at different stages of the SRt2 selections for the $t \bar{t} +\phi$ benchmark model $m(\phi,\chi)=(20,1)$GeV, $g=1$.

Number of signal events selected at different stages of the SRt1 selections for the $t \bar{t} +a$ benchmark model $m(a,\chi)=(20,1)$GeV, $g=1$.

Number of signal events selected at different stages of the SRt2 selections for the $t \bar{t} +a$ benchmark model $m(\phi,\chi)=(20,1)$GeV, $g=1$.

Number of signal events selected at different stages of the SRt3 selections for the $t \bar{t} +\phi$ benchmark model $m(\phi,\chi)=(20,1)$GeV, $g=1$.

Figure 9a (upper left). The 95% CL upper limit on the production cross section of the T1tttt simplified models as a function of the gluino and LSP masses.

Figure 9a (upper left). Observed exclusion region at 95% CL assuming 100% branching fraction.

Figure 9a (upper left). Expected exclusion region at 95% CL assuming 100% branching fraction.

Figure 9b (upper right). The 95% CL upper limit on the production cross section of the T1ttbb simplified models as a function of the gluino and LSP masses.

Figure 9b (upper right). Observed exclusion region at 95% CL.

Figure 9b (upper right). Expected exclusion region at 95%.

Figure 9c (bottom left). The 95% CL upper limit on the production cross section of the T5tttt simplified models as a function of the gluino and LSP masses. imits are not given for the T5tttt model for mχ~01< 50 GeV for the reason stated in the text.

Figure 9c (bottom left). Observed exclusion region at 95% CL assuming 100% branching fraction.

Figure 9c (bottom left). Expected exclusion region at 95% CL assuming 100% branching fraction.

Figure 9d (bottom right). The 95% CL upper limit on the production cross section of the T5ttcc simplified models as a function of the gluino and LSP masses.

Figure 9d (bottom right). Observed exclusion region at 95% CL.

Figure 9d (bottom right). Expected exclusion region at 95% CL.

The observed number of events and the total background prediction for search regions with Nt = 1 and Nb = 1.

The observed number of events and the total background prediction for search regions with Nt = 1 and Nb ≥ 2.

The observed number of events and the total background prediction for search regions with Nt ≥ 2.

Selected T2tt signal yields for the aggregate search bins.

Selected gluino mediated signal yields for the aggregate search bins.