The expected 95% CL exclusion contour (-1$\sigma$) for a simplified model of dark-matter production involving an axial-vector operator, Dirac DM and couplings $g_{q}$ = 0.25, $g_{\chi}$ = 1 and $g_{l}$ = 0 as a function of the dark-matter mass $m_{\chi}$ and the mediator mass $m_{\mathrm{med}}$. The plane under the limit curve is excluded.

The expected 95% CL exclusion contour for a simplified model of dark-matter production involving an axial-vector operator, Dirac DM and couplings $g_{q}$ = 0.25, $g_{\chi}$ = 1 and $g_{l}$ = 0 as a function of the dark-matter mass $m_{\chi}$ and the mediator mass $m_{\mathrm{med}}$. The plane under the limit curve is excluded.

The observed 95% CL exclusion contour for a simplified model of dark-matter production involving an axial-vector operator, Dirac DM and couplings $g_{q}$ = 0.1, $g_{\chi}$ = 1 and $g_{l}$ = 0.1 as a function of the dark-matter mass $m_{\chi}$ and the mediator mass $m_{\mathrm{med}}$. The plane under the limit curve is excluded.

The expected 95% CL exclusion contour (+1$\sigma$) for a simplified model of dark-matter production involving an axial-vector operator, Dirac DM and couplings $g_{q}$ = 0.1, $g_{\chi}$ = 1 and $g_{l}$ = 0.1 as a function of the dark-matter mass $m_{\chi}$ and the mediator mass $m_{\mathrm{med}}$. The plane under the limit curve is excluded.

The expected 95% CL exclusion contour (-1$\sigma$) for a simplified model of dark-matter production involving an axial-vector operator, Dirac DM and couplings $g_{q}$ = 0.1, $g_{\chi}$ = 1 and $g_{l}$ = 0.1 as a function of the dark-matter mass $m_{\chi}$ and the mediator mass $m_{\mathrm{med}}$. The plane under the limit curve is excluded.

The expected 95% CL exclusion contour for a simplified model of dark-matter production involving an axial-vector operator, Dirac DM and couplings $g_{q}$ = 0.1, $g_{\chi}$ = 1 and $g_{l}$ = 0.1 as a function of the dark-matter mass $m_{\chi}$ and the mediator mass $m_{\mathrm{med}}$. The plane under the limit curve is excluded.

The observed 95% CL exclusion contour for a simplified model of dark-matter production involving a vector operator with couplings $g_{q}$ = 0.25, $g_{\chi}$ =1 and $g_{l}$ = 0 as a function of the dark-matter mass $m_{\chi}$ and the mediator mass $m_{\mathrm{med}}$. The plane under the limit curve is excluded.

The expected 95% CL exclusion contour ($+1\sigma$) for a simplified model of dark-matter production involving a vector operator with couplings $g_{q}$ = 0.25, $g_{\chi}$ =1 and $g_{l}$ = 0 as a function of the dark-matter mass $m_{\chi}$ and the mediator mass $m_{\mathrm{med}}$. The plane under the limit curve is excluded.

The expected 95% CL exclusion contour ($-1\sigma$) for a simplified model of dark-matter production involving a vector operator with couplings $g_{q}$ = 0.25, $g_{\chi}$ =1 and $g_{l}$ = 0 as a function of the dark-matter mass $m_{\chi}$ and the mediator mass $m_{\mathrm{med}}$. The plane under the limit curve is excluded.

The expected 95% CL exclusion contour for a simplified model of dark-matter production involving a vector operator with couplings $g_{q}$ = 0.25, $g_{\chi}$ =1 and $g_{l}$ = 0 as a function of the dark-matter mass $m_{\chi}$ and the mediator mass $m_{\mathrm{med}}$. The plane under the limit curve is excluded.

The observed 95% CL exclusion contour for a simplified model of dark-matter production involving for a vector operator with couplings $g_{q}$ = 0.1, $g_{\chi}$ = 1 and $g_{l}$ = 0.01 as a function of the dark-matter mass $m_{\chi}$ and the mediator mass $m_{\mathrm{med}}$. The plane under the limit curve is excluded.

The expected 95% CL exclusion contour ($+1\sigma$) for a simplified model of dark-matter production involving for a vector operator with couplings $g_{q}$ = 0.1, $g_{\chi}$ = 1 and $g_{l}$ = 0.01 as a function of the dark-matter mass $m_{\chi}$ and the mediator mass $m_{\mathrm{med}}$. The plane under the limit curve is excluded.

The expected 95% CL exclusion contour ($-1\sigma$) for a simplified model of dark-matter production involving for a vector operator with couplings $g_{q}$ = 0.1, $g_{\chi}$ = 1 and $g_{l}$ = 0.01 as a function of the dark-matter mass $m_{\chi}$ and the mediator mass $m_{\mathrm{med}}$. The plane under the limit curve is excluded.

The expected 95% CL exclusion contour for a simplified model of dark-matter production involving for a vector operator with couplings $g_{q}$ = 0.1, $g_{\chi}$ = 1 and $g_{l}$ = 0.01 as a function of the dark-matter mass $m_{\chi}$ and the mediator mass $m_{\mathrm{med}}$. The plane under the limit curve is excluded.

The 90% CL exclusion limit on the $\chi$-proton scattering cross section in a simplifed model of dark-matter production involving an axial-vector operator, Dirac DM and couplings $g_{q}$ = 0.25, $g_{\chi}$ = 1 and $g_{l}$ = 0 as a function of the dark-matter mass $m_{\chi}$.

The 90% CL exclusion limit on the $\chi$-nucleon scattering cross section in a simplifed model of dark-matter production involving a vector operator, Dirac DM and couplings $g_{q}$ = 0.25, $g_{\chi}$ = 1 and $g_{l}$ = 0 as a function of the dark-matter mass $m_{\chi}$.

The observed and expected 95% CL limits on M* for a dimension-7 operator EFT model with a contact interaction of type $\gamma\gamma\chi\chi$ as a function of dark-matter mass $m_{\chi}$.

The observed and expected limit at 95% CL on the production cross section of a $Z\gamma$ resonance as a function of its mass.

The observed 95% CL upper limits on signal strength $\mu$, the ratio of the experimental limit to the predicted signal cross section, 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 DM and the mediator particle masses. The bound on $\mu$ only applies for choices of couplings that yield the same kinematic distributions as the benchmark model.

Acceptance, in %, 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 for various $m_{\chi}$ and $m_{\textrm{med}}$ values. It is constant as a function of $m_{\chi}$ for a fixed $m_{\textrm{med}}$ in the region $m_{\textrm{med}}>2m_{\chi}$.

Fiducial reconstruction efficiencies in the three inclusive signal regions 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 for various $m_{\chi}$ and $m_{\textrm{med}}$ values. It is constant as a function of $m_{\chi}$ for a fixed $m_{\textrm{med}}$ in the region $m_{\textrm{med}}>2m_{\chi}$.

DM exclusion limit in the two-dimensional phase space of WIMP mass m_χ vs mediator mass m_med determined using the combined ee+μμ channel. Both the observed and expected limits are presented, and the 1σ uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_χ'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs, and the other grey line gives the perturbative limit (arXiv 1603.04156). The relic density line (arXiv 1603.04156) illustrates the combination of m_χ and m_med that would explain the observed DM relic density.

DM exclusion limit in the two-dimensional phase space of WIMP mass m_χ vs mediator mass m_med determined using the combined ee+μμ channel. Both the observed and expected limits are presented, and the 1σ uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_χ'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs, and the other grey line gives the perturbative limit (arXiv 1603.04156). The relic density line (arXiv 1603.04156) illustrates the combination of m_χ and m_med that would explain the observed DM relic density.

DM exclusion limit in the two-dimensional phase space of WIMP mass m_χ vs mediator mass m_med determined using the combined ee+μμ channel. Both the observed and expected limits are presented, and the 1σ uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_χ'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs, and the other grey line gives the perturbative limit (arXiv 1603.04156). The relic density line (arXiv 1603.04156) illustrates the combination of m_χ and m_med that would explain the observed DM relic density.

The 95% CL upper limits on the production cross-section of the ZH process with the prompt Z→ ee and Z→ μμ decays and the invisible Higgs boson decays as a function of m_H, obtained from the combined ee+μμ channel. The observed and expected limits are given, as well as the ±1σ and ±2σ error bands on the expected limits. The signal process is only modelled for the qq initial state, and the theory uncertainties on the signal prediction are not considered.

DM exclusion limit in the two-dimensional phase space of m_χ vs m_med set using the combined ee+μμ channel. The DM model assumes a vector mediator, a fermionic WIMP, and the coupling parameters g_q = 0.25 and g_χ = 1. Both the observed and expected limits are presented, and the 1σ uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_χ'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs. The relic density line [arXiv 1603.04156] illustrates the combination of m_χ and m_med that would explain the observed DM relic density.

DM exclusion limit in the two-dimensional phase space of m_χ vs m_med set using the combined ee+μμ channel. The DM model assumes a vector mediator, a fermionic WIMP, and the coupling parameters g_q = 0.25 and g_χ = 1. Both the observed and expected limits are presented, and the 1σ uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_χ'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs. The relic density line [arXiv 1603.04156] illustrates the combination of m_χ and m_med that would explain the observed DM relic density.

DM exclusion limit in the two-dimensional phase space of m_χ vs m_med set using the combined ee+μμ channel. The DM model assumes a vector mediator, a fermionic WIMP, and the coupling parameters g_q = 0.25 and g_χ = 1. Both the observed and expected limits are presented, and the 1σ uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_χ'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs. The relic density line [arXiv 1603.04156] illustrates the combination of m_χ and m_med that would explain the observed DM relic density.

DM exclusion limit in the two-dimensional phase space of m_χ vs m_med set using the combined ee+μμ channel. The DM model assumes a vector mediator, a fermionic WIMP, and the coupling parameters g_q = 0.25 and g_χ = 1. Both the observed and expected limits are presented, and the 1σ uncertainty band for the expected limits is also provided. Regions bounded by the limit curves are excluded at the 95% CL. The grey line labelled with "m_med = 2m_χ'' indicates the kinematic threshold where the mediator can decay on-shell into WIMPs. The relic density line [arXiv 1603.04156] illustrates the combination of m_χ and m_med that would explain the observed DM relic density.

Limits on the WIMP and proton scattering cross section as a function of WIMP mass, obtained for simplified dark matter models with an axial-vector mediator (left) and a vector mediator (right). The solid black line shows the observed limit at the 90% confidence level from this search. The left plot incorporates the latest results from these competitive direct-search experiments, LUX [arXiv 1705.03380], PICO-2L [Phys. Rev. Lett. 114 (2015) 231302], and PICO-60 [arXiv 1702.07666]. The right plot incorporates the latest results from the following experiments, CRESST-II [Eur. Phys. J. C76 (2016) 25], CDMSlite [Phys. Rev. Lett. 116 (2016) 071301], PandaX-II [Phys. Rev. Lett. 117 (2016) 121303], LUX [Phys. Rev. Lett. 118 (2017) 021303], and XENON1T [arXiv 1705.06655].

Limits on the WIMP and proton scattering cross section as a function of WIMP mass, obtained for simplified dark matter models with an axial-vector mediator (left) and a vector mediator (right). The solid black line shows the observed limit at the 90% confidence level from this search. The left plot incorporates the latest results from these competitive direct-search experiments, LUX [arXiv 1705.03380], PICO-2L [Phys. Rev. Lett. 114 (2015) 231302], and PICO-60 [arXiv 1702.07666]. The right plot incorporates the latest results from the following experiments, CRESST-II [Eur. Phys. J. C76 (2016) 25], CDMSlite [Phys. Rev. Lett. 116 (2016) 071301], PandaX-II [Phys. Rev. Lett. 117 (2016) 121303], LUX [Phys. Rev. Lett. 118 (2017) 021303], and XENON1T [arXiv 1705.06655].

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-W/Z signal region with the merged event topology after the profile-likelihood fit. This is shown for the 0b-tagged jet, low purity, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-W/Z signal region with the merged event topology after the profile-likelihood fit. This is shown for the 1b-tagged jet, high purity, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-W/Z signal region with the merged event topology after the profile-likelihood fit. This is shown for the 1b-tagged jet, low purity, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-W/Z signal region with the merged event topology after the profile-likelihood fit. This is shown for the 2b-tagged jet event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-W/Z signal region with the resolved event topology after the profile-likelihood fit. This is shown for the 0b-tagged jet event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-W/Z signal region with the resolved event topology after the profile-likelihood fit. This is shown for the 1b-tagged jet event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-W/Z signal region with the resolved event topology after the profile-likelihood fit. This is shown for the 2b-tagged jet event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the merged event topology for a Z' mass of 90 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, high purity, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the merged event topology for a Z' mass of 90 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, low purity, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the merged event topology for a Z' mass of 90 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, high purity, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the merged event topology for a Z' mass of 90 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, low purity, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the merged event topology for a Z' mass of 90 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 90 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 350 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 90 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 350 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 90 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 350 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

Expected exclusion contours at 95% C.L. of the dark matter mediator particles (m_chi, m_Zp) for the combined mono-W and mono-Z search in the frame of the simplified model with the Dirac DM particles and couplings g_SM = 0.25 and g_DM = 1.

Observed exclusion contours at 95% C.L. of the dark matter mediator particles (m_chi, m_Zp) for the combined mono-W and mono-Z search in the frame of the simplified model with the Dirac DM particles and couplings g_SM = 0.25 and g_DM = 1.

Exclusion contours at 95% C.L. of the cross section as a function of Z' mass (mZp), for the mono-Z' dark fermion model, with the light dark sector scenario, with parameters g_SM = 0.1 and g_DM = 1.

Exclusion contours at 95% C.L. of the cross section as a function of Z' mass (mZp), for the mono-Z' dark fermion model, with the heavy dark sector scenario, with parameters g_SM = 0.1 and g_DM = 1.

Exclusion contours at 95% C.L. of the cross section as a function of Z' mass (mZp), for the mono-Z' dark Higgs model, with the light dark sector scenario, with parameters g_SM = 0.1 and g_DM = 1.

Exclusion contours at 95% C.L. of the cross section as a function of Z' mass (mZp), for the mono-Z' dark Higgs model, with the heavy dark sector scenario, with parameters g_SM = 0.1 and g_DM = 1.

Exclusion contours at 95% C.L. of the product of the couplings g_SM g_DM as a function of Z' mass (mZp), for the mono-Z' dark fermion model, with the light dark sector scenario, with parameter g_DM = 1.

Exclusion contours at 95% C.L. of the product of the couplings g_SM g_DM as a function of Z' mass (mZp), for the mono-Z' dark fermion model, with the heavy dark sector scenario, with parameter g_DM = 1.

Exclusion contours at 95% C.L. of the product of the couplings g_SM g_DM as a function of Z' mass (mZp), for the mono-Z' dark Higgs model, with the light dark sector scenario, with parameter g_DM = 1.

Exclusion contours at 95% C.L. of the product of the couplings g_SM g_DM as a function of Z' mass (mZp), for the mono-Z' dark Higgs model, with the heavy dark sector scenario, with parameter g_DM = 1.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the merged event topology for a Z' mass of 80 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, low purity, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the merged event topology for a Z' mass of 80 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, low purity, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the merged event topology for a Z' mass of 80 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, high purity, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the merged event topology for a Z' mass of 80 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, high purity, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the merged event topology for a Z' mass of 80 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 80 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 80 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 80 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the merged event topology for a Z' mass of 100 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, low purity, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the merged event topology for a Z' mass of 100 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, low purity, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the merged event topology for a Z' mass of 100 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, high purity, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the merged event topology for a Z' mass of 100 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, high purity, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the merged event topology for a Z' mass of 100 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 100 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 100 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 100 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 110 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 110 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 110 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 120 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 120 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 120 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 130 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 130 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 130 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 140 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 140 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 140 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 150 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 150 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 150 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 160 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 160 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 160 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 170 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 170 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 170 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 180 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 180 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 180 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 190 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 190 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 190 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 200 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 200 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 200 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 250 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 250 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 250 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 300 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 300 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 300 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 400 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 400 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 400 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 450 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 450 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 450 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 500 GeV after the profile-likelihood fit. This is shown for the 0b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 500 GeV after the profile-likelihood fit. This is shown for the 1b-tagged jet, event category.

The observed and expected MET distributions with 36.1fb-1 of data with sqrt(s) = 13 TeV in the mono-Z' signal region with the resolved event topology for a Z' mass of 500 GeV after the profile-likelihood fit. This is shown for the 2b-tagged jet, event category.

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.

The ratios (μ) of the 95% C.L. upper limits on the combined s→ W^±W^∓ and s→ ZZ cross section to simplified model expectations for the m_Z'=1 TeV scenario, for various m_s hypotheses. The observed limits (solid line) are consistent with the expectation under the SM-only hypothesis (dashed line) within uncertainties (filled band), except for a small excess for m_s=160 GeV, discussed in the text.

The ratios (μ) of the 95% C.L. upper limits on the combined s→ W^±W^∓ and s→ ZZ cross section to simplified model expectations for the m_Z'=1.7 TeV scenario, for various m_s hypotheses. The observed limits (solid line) are consistent with the expectation under the SM-only hypothesis (dashed line) within uncertainties (filled band), except for a small excess for m_s=160 GeV, discussed in the text.

Observed upper limits at 95% C.L. on σ(pp → s χχ) × B(s→ VV) for m_Z'=0.5 TeV signal points. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the σ(pp → s χχ) × B(s→ VV) process for m_Z'=0.5 TeV, shown in dashed blue.

Observed upper limits at 95% C.L. on σ(pp → s χχ) × B(s→ VV) for m_Z'=1 TeV signal points. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the σ(pp → s χχ) × B(s→ VV) process for m_Z'=1 TeV, shown in dashed blue.

Observed upper limits at 95% C.L. on σ(pp → s χχ) × B(s→ VV) for m_Z'=1.7 TeV signal points. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the σ(pp → s χχ) × B(s→ VV) process for m_Z'=1.7 TeV, shown in dashed blue.

SM background post-fit yields stacked in each SR and CR category and E_T^miss bin and data overlaid with the maximum likelihood estimators set to the conditional values of the combined signal and control region fit. The hatched uncertainty band shown includes simulation statistics uncertainties, experimental systematic uncertainties, and V+jets theory modelling systematic uncertainties. Pre-fit uncertainties cover differences between the data and pre-fit background prediction.

Cumulative efficiencies for the merged category for signal samples with m_s=160 GeV (a), m_s=235 GeV (b) and m_s=310 GeV (c), each with m_Z'=1 TeV. The dark Higgs candidate selection includes stringent jet substructure requirements and typically at most one candidate is present in signal events. Here, Δ φ_{jets_1,2,3 E_T^miss} is the smallest azimuthal angle between the E_T^miss and any of the three highest-p_T (leading) small-R jets.

Cumulative efficiencies for the intermediate category for signal samples with m_s=160 GeV (a), m_s=235 GeV (b) and m_s=310 GeV (c), each with m_Z'=1 TeV. The TAR+Comb algorithm reconstructs the dark Higgs candidate from a TAR jet with m^TAR>60 GeV that is supplemented by up to two additional small-R jets within ΔR_cone=2.5 of the TAR jet. Here, Δ φ_{jets_1,2,3 E_T^miss} is the smallest azimuthal angle between the E_T^miss and any of the three highest-p_T (leading) small-R jets. For details see text.

The product of acceptance and efficiency (A × ϵ), defined as the number of signal events satisfying the full set of selection criteria in the merged or intermediate signal regions, divided by the total number of generated signal events, for the s(W^±W^∓) dark Higgs signal points with dark Higgs boson mass m_s and Z' boson mass m_Z'.

The product of acceptance and efficiency (A × ϵ), defined as the number of signal events satisfying the full set of selection criteria in the merged or intermediate signal regions, divided by the total number of generated signal events, for the s(ZZ) dark Higgs signal points with dark Higgs boson mass m_s and Z' boson mass m_Z'.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 0L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 1L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Only signals simulating the tW+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_a$ vs. $m_{H^{\pm}}$ and assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 150 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The observed exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

The expected exclusion contour at 95% CL as a function of the $m_{H^{\pm}}$ vs. tan$\beta$ and assuming $m_a$ = 250 $\mathrm{GeV}$, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Masses that are within the contours are excluded. Signals simulating the tW+DM + tt+DM final states are considered in this contour. These exclusion contours are derived using the 2L channel only.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.7$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Only signals simulating the tW+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_a$ vs. $ m_{H^{\pm}}$ signal grid assuming tan$\beta$ = 1, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 150 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

Model dependent upper limit on the cross section for the $m_{H^{\pm}}$ vs. tan$\beta$ signal grid assuming $m_a$ = 250 GeV, $m_{\mathrm{DM}} = 10 \mathrm{GeV}$, $g_{\chi} = 1$ and sin$\theta = 0.35$. Signals simulating the tW+DM + tt+DM final states are considered. Upper limits with large $\mu_{\mathrm{sig}}$ for the observed limit are capped at 500.

The distributions of $m_{\mathrm{b1},\mathrm{W-tagged}}$ in the 0L inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{T}}^{\mathrm{b,E_{\mathrm{T}^{\mathrm{miss}}}}}$ in the 0L inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $N_{\mathrm{W-tagged}}$ in the 0L inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{b1},\mathrm{\cancel{b1}}}$ in the hadronic top inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{b1},\mathrm{\cancel{b1}}}$ in the leptonic top inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

The distributions of $m_{\mathrm{b1},\mathrm{\cancel{b1}}}$ in the leptonic top inclusive signal region. For each bin yields for the data and total SM prediction are provided. The SM prediction is provided with the total uncertainty, including the MC statistical uncertainty, detector-related systematic uncertainties and theoretical uncertainties. The rightmost bin includes overflow events.

Cutflow for the reference point $(\it{m}_{\mathrm{H^{\pm}}}, \it{m}_{a}, tan\beta, sin\theta )=$ (500,100,1,0.7) , (800,150,20,0.7), (600,250,30,0.7), (1000,400,1,0.7) in 0L regions. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(\it{m}_{\mathrm{H^{\pm}}}, \it{m}_{a}, tan\beta, sin\theta )=$ (500,100,1,0.7) , (800,150,20,0.7), (600,250,30,0.7), (1000,400,1,0.7) in 1L leptonic top regions. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Cutflow for the reference point $(\it{m}_{\mathrm{H^{\pm}}}, \it{m}_{a}, tan\beta, sin\theta )=$ (500,100,1,0.7) , (800,150,20,0.7), (600,250,30,0.7), (1000,400,1,0.7) in 1L hadronic top regions. Results are shown including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.7. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 0L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_a$--m$_{H^{\pm}}$ assuming tan$\beta$ = 1, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 150 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Signal acceptance in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the acceptance given in the table is multiplied by factor of $10^{3}$

Signal efficiency in the 1L region for 2HDM+a model DM signals on the plane defined by m$_{H^{\pm}}$--tan$\beta$ assuming m$_a$ = 250 GeV, m$_{\chi}$= 10 GeV and sin$\theta$ = 0.35. Please mind that the efficiency given in the table is multiplied by factor of $10^{2}$

Probability of finding at least one TAR jet, where the p_T-leading TAR jet passes the m_Wcand and D₂^β=1 requirements, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=2100 GeV, with the preselections applied.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=500 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=1000 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=1700 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=2100 GeV, with the preselections applied that do not pass the requirements of the merged category.

Observed exclusion contour at 95% C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_χ=1, m_χ=200 GeV, and sinθ=0.01.

Expected exclusion contour at 95% C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_χ=1, m_χ=200 GeV, and sinθ=0.01.

Expected+1σ exclusion contour at 95% C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_χ=1, m_χ=200 GeV, and sinθ=0.01.

Expected-1σ exclusion contour at 95% C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_χ=1, m_χ=200 GeV, and sinθ=0.01.

Expected+2σ exclusion contour at 95% C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_χ=1, m_χ=200 GeV, and sinθ=0.01.

Expected-2σ exclusion contour at 95% C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_χ=1, m_χ=200 GeV, and sinθ=0.01.

Observed upper limits at 95% C.L. on σ(pp → s χχ) × B(s → W^± W^∓) for m_Z'=0.5 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the σ(pp → s χχ) × B(s → W^± W^∓) process for m_Z'=0.5 TeV, shown in dashed blue.

Observed upper limits at 95% C.L. on σ(pp → s χχ) × B(s → W^± W^∓) for m_Z'=1 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the σ(pp → s χχ) × B(s → W^± W^∓) process for m_Z'=1 TeV, shown in dashed blue.

Observed upper limits at 95% C.L. on σ(pp → s χχ) × B(s → W^± W^∓) for m_Z'=1.7 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the σ(pp → s χχ) × B(s → W^± W^∓) process for m_Z'=1.7 TeV, shown in dashed blue.

Observed upper limits at 95% C.L. on σ(pp → s χχ) × B(s → W^± W^∓) for m_Z'=2.1 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the σ(pp → s χχ) × B(s → W^± W^∓) process for m_Z'=2.1 TeV, shown in dashed blue.

Data overlaid on SM background yields stacked in each SR and CR category after the fit to data ('Post-fit'). The yields in the SR are broken down into their contributions to the individual bins. The maximum-likelihood estimators are set to the conditional values of the CR-only fit, and propagated to SR and CRs.

Dominant sources of uncertainty for three dark Higgs scenarios after the fit to data. The uncertainties are quantified in terms of their contribution to the fitted signal uncertainty that is expressed relative to the theory prediction. Three representative dark Higgs signal scenarios with g_q=0.25, g_χ=1.0, sinθ=0.01 and m_χ=200 GeV are considered, which are indicated using the (m_Z', m_s) format in units of GeV in the table columns.

Cumulative efficiencies in the merged category for three representative dark Higgs signal scenarios with g_q=0.25, g_χ=1.0, sinθ=0.01, m_Z' = 1 TeV, and m_χ=200 GeV considering s→W(ℓν)W(qq) decays only.

Cumulative efficiencies in the resolved category for three representative dark Higgs signal scenarios with g_q=0.25, g_χ=1.0, sinθ=0.01, m_Z' = 1 TeV, and m_χ=200 GeV considering s→W(ℓν)W(qq) decays only.

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={300, 350, 400, 500, 750} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={1000, 1700} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={2100, 2500, 2900, 3300} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The j+b-jet invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The 2b-jet invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The j+e invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The b-jet+e invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The j+$\gamma$ invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The j+$\mu$ invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The b-jet+$\mu$ invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The b-jet+$\gamma$ invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

The 95% CL upper limits on the cross section times acceptance ($A$), efficiency ($\epsilon$), and branching ratio ($B$) for Gaussian-shaped signals with different signal widths. The limits are calculated for events with at least one lepton with pT > 60 GeV in the 10 pb AR. Two width hypotheses are shown; $\sigma X / mX = 0$ and $\sigma X / mX = 0.15$. In both cases, the detector resolution for jets is included in the simulation of signal samples. The $\pm1 \sigma$ and $\pm2 \sigma$ bands around the expected limit are shown for $\sigma X / mX = 0$ signals. Mass points are spaced 5% apart, relative to the preceding point, starting at 0.3 TeV.

Distributions of the anomaly score from the AE for data and five benchmark BSM models. Their legends, from top to bottom, are; (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$, with the $H^{+}$ mass of 2 TeV; (2) a 6 TeV Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a 4 TeV $Z'$ boson decaying to a composite lepton $E$ with a mass of 0.5 TeV; (4) the sequential standard model $W' \rightarrow WZ' \rightarrow \ell\nu q\bar{q}$, with the $W'$ mass of 6.2 TeV and the $Z'$ mass of 6 TeV; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$ with a mass of 6 TeV. The distributions for the BSM models are smoothed to remove fluctuations due to low MC event counts. The vertical lines indicate the start of the three anomaly regions (ARs). The labels of the three ARs indicate the visible cross section for hypothetical processes yielding the same number of events as observed in the 140 $fb^{-1}$ dataset. The AE is applied to preselected events without any requirements on invariant mass distributions.

Invariant mass distributions of $j+Y$ for $M_{jY}$ > 0.3 TeV after preselection along with the fit from Eq.(1). The fit is represented by red lines, and the associated uncertainties are indicated by yellow bands. The lower panels show the bin-by-bin significances of deviations from the fit, calculated as $(d_{\textit{i}} - f_{\textit{i}}) / \delta_{\textit{i}}$, where $d_{\textit{i}}$ is the data yield, $f_{\textit{i}}$ is the fit value, and $\delta_{i}$ is the uncertainty of data in the $\textit{i}$-th bin

Distributions of the anomaly score for data and several anomaly scenarios. The example BSM model (shown with the dashed blue lines) is the sequential standard model $W' \rightarrow WZ' \rightarrow \ell\nu q\bar{q}$. The mass of $W'$ is set to 2.2 TeV and the mass of the $Z'$ is set to 2 TeV. This model leads to the final state of one lepton, two jets, and small missing transverse energy that is similar to the SM backgrounds. The other histograms represent events from the same model with artificial modifications to represent anomalous events; (1) "Anomaly 1" is the case where all jets beyond the second jet ($N_{j}$ > 2) are replaced with photons; (2) "Anomaly 2" is the case where all jets beyond the second jet are replaced with $b$-jets; (3) "Anomaly 3" is the case of low-multiplicity events where all jets beyond the second jet are removed. The histograms are normalized to unit. The left peak near log(Loss) = -9 visible in Anomaly 1 and 2 is from the events without a third jet.