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.

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.

Expected (1 sigma band) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Expected (2 sigma band) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Expected (2 sigma band) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Observed exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Observed exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Observed (+1 sigma) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Observed (+1 sigma) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Observed (-1 sigma) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Observed (-1 sigma) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Exclusion limits on the production cross section as a function of m($\tilde{t}$) are shown for several values of $\tau(\tilde{t})$ along with the nominal signal production cross section and its theoretical uncertainty.

Exclusion limits on the production cross section as a function of m($\tilde{t}$) are shown for several values of $\tau(\tilde{t})$ along with the nominal signal production cross section and its theoretical uncertainty.

Parameterized event selection efficiencies for the $E_{T}^{miss}$ Trigger SR. The event-level efficiencies for each SR are extracted for all events passing the acceptance of the corresponding SR.

Parameterized event selection efficiencies for the $E_{T}^{miss}$ Trigger SR. The event-level efficiencies for each SR are extracted for all events passing the acceptance of the corresponding SR.

Parameterized event selection efficiencies for the Muon Trigger SR. The event-level efficiencies for each SR are extracted for all events passing the acceptance of the corresponding SR.

Parameterized event selection efficiencies for the Muon Trigger SR. The event-level efficiencies for each SR are extracted for all events passing the acceptance of the corresponding SR.

Parameterized muon-level reconstruction efficiencies as a function of the muon $p_{T}$ and $d_{0}$. The muon-level efficiencies are extracted using muons passing the muon acceptance criteria.

Parameterized muon-level reconstruction efficiencies as a function of the muon $p_{T}$ and $d_{0}$. The muon-level efficiencies are extracted using muons passing the muon acceptance criteria.

Parameterized vertex-level reconstruction efficiencies as a function of the radial position of the truth vertex. The efficiency is calculated independent of the muons originating from this truth vertex.

Parameterized vertex-level reconstruction efficiencies as a function of the radial position of the truth vertex. The efficiency is calculated independent of the muons originating from this truth vertex.

Parameterized vertex-level reconstruction efficiencies as a function of the radial position of the truth vertex. The efficiency is calculated only for truth vertices which have a muon originating from them which is matched to a reconstructed muon.

Parameterized vertex-level reconstruction efficiencies as a function of the radial position of the truth vertex. The efficiency is calculated only for truth vertices which have a muon originating from them which is matched to a reconstructed muon.

The $p_{T}$ distribution of all muons originating from LLP decays in the samples used to calculate and validate the efficiencies.

The $p_{T}$ distribution of all muons originating from LLP decays in the samples used to calculate and validate the efficiencies.

The invariant mass and multiplicity of selected decay products of all truth vertices used in the calculation and validation of the reconstructed efficiencies.

The invariant mass and multiplicity of selected decay products of all truth vertices used in the calculation and validation of the reconstructed efficiencies.

Upper limits at the 95% CL on the cross section for SSM $W^\prime$ production and decay to the muon+neutrino channel as a function of the $W^\prime$ pole mass.

Combined (electron and muon channels) upper limits at the 95% CL on the cross section for SSM $W^\prime$ production and decay to a single lepton generation as a function of the $W^\prime$ pole mass.

Observed upper limits at the 95% CL on the cross section for generic $W^\prime$ production and decay to the electron+neutrino channel as a function of the $W^\prime$ pole mass.

Observed upper limits at the 95% CL on the cross section for generic $W^\prime$ production and decay to the muon+neutrino channel as a function of the $W^\prime$ pole mass.

Combined (electron and muon channels) observed upper limits at the 95% CL on the cross section for generic $W^\prime$ production and decay to a single lepton generation as a function of the $W^\prime$ pole mass.

Observed upper limits at the 95% CL on the visible cross section in the electron+neutrino channel as a function of the transverse mass threshold.

Observed upper limits at the 95% CL on the visible cross section in the muon+neutrino channel as a function of the transverse mass threshold.

Product of acceptance and efficiency for the electron and muon selections as a function of the $W^\prime$ pole mass.

Expected upper limits at the 95% CL on the cross section for generic $W^\prime$ production and decay to the electron+neutrino channel as a function of the $W^\prime$ pole mass.

Expected upper limits at the 95% CL on the cross section for generic $W^\prime$ production and decay to the muon+neutrino channel as a function of the $W^\prime$ pole mass.

Combined (electron and muon channels) expected upper limits at the 95% CL on the cross section for generic $W^\prime$ production and decay to a single lepton generation as a function of the $W^\prime$ pole mass.

Distribution of the resonance mass in the c-tag Signal Region of the $ ce$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the b-tag Signal Region of the $ ce$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the untagged Signal Region of the $ c\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the c-tag Signal Region of the $ c\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the b-tag Signal Region of the $ c\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 0-tag Signal Region of the $ be$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 1-tag Signal Region of the $ be$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 2-tag Signal Region of the $ be$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 0-tag Signal Region of the $ b\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 1-tag Signal Region of the $ b\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 2-tag Signal Region of the $ b\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into electrons, shown as a function of $m_{LQ}$ for the $qe$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into muons, shown as a function of $m_{LQ}$ for the $q\mu$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into electrons, shown as a function of $m_{LQ}$ for the $ce$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into muons, shown as a function of $m_{LQ}$ for the $c\mu$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into electrons, shown as a function of $m_{LQ}$ for the $be$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into muons, shown as a function of $m_{LQ}$ for the $b\mu$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $qe$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $q\mu$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $ce$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $c\mu$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $be$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $b\mu$ channel.

The signal selection efficiency x acceptance summed over all signal regions, for all masses and LQ decay channels considered.

The observed and expected limits for all masses and LQ decay channels considered.

Expected exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Plus 1$\sigma$ uncertainty, varying the signal cross section within its uncertainty, on the observed exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Plus 1$\sigma$ uncertainty, varying the signal cross section within its uncertainty, on the observed exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Minus 1$\sigma$ uncertainty, varying the signal cross section within its uncertainty, on the observed exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Minus 1$\sigma$ uncertainty, varying the signal cross section within its uncertainty, on the observed exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Plus 1$\sigma$ uncertainty, due to uncertainties in the background prediction and experimental uncertainties affecting the signal, on the expected exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Plus 1$\sigma$ uncertainty, due to uncertainties in the background prediction and experimental uncertainties affecting the signal, on the expected exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Minus 1$\sigma$ uncertainty, due to uncertainties in the background prediction and experimental uncertainties affecting the signal, on the expected exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Minus 1$\sigma$ uncertainty, due to uncertainties in the background prediction and experimental uncertainties affecting the signal, on the expected exclusion contour on C1N2 production assuming on-shell $W/Z$ decays as a function of the C1/N2 and N1 masses, and derived from the combined fit of low-mass and ISR regions.

Upper limits on observed wino-bino simplified model signal cross section $\sigma_\text{obs}^\text{95}$.

Upper limits on observed wino-bino simplified model signal cross section $\sigma_\text{obs}^\text{95}$.

Upper limits on expected wino-bino simplified model signal cross section $\sigma_\text{exp}^\text{95}$.

Upper limits on expected wino-bino simplified model signal cross section $\sigma_\text{exp}^\text{95}$.

Signal acceptance in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal acceptance in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal efficiency in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal efficiency in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal acceptance in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal acceptance in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal efficiency in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal efficiency in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} \geq 100$ GeV.

Signal acceptance in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

Signal acceptance in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

Signal efficiency in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

Signal efficiency in SR-low, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

Signal acceptance in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

Signal acceptance in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

Signal efficiency in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

Signal efficiency in SR-ISR, for signals with $m(\widetilde{\chi}^{\pm}_{1}/\widetilde{\chi}^{0}_{2}) - m\widetilde{\chi}^{0}_{1} < 100$ GeV.

Upper limits on observed (expected) wino-bino simplified model signal cross section $\sigma_\text{obs(exp)}^\text{95}$.

Upper limits on observed (expected) wino-bino simplified model signal cross section $\sigma_\text{obs(exp)}^\text{95}$.

The observed upper limits on the model cross-section in units of pb for simplified models with ${\tilde{\tau}}_L {\tilde{\tau}}_L$ only production. Three points at $M({\tilde{\chi}}^{0}_{1})=200GeV$ were removed from the plot but kept in the table because they overlapped with the plot's legend and are far from the exclusion contour.

The observed 95\% CL exclusion contours for the combined fit of SR-lowMass and SR-highMass for simplified models with combined ${\tilde{\tau}}^{+}_{R,L} {\tilde{\tau}}^{-}_{R,L}$ production.

The observed 95\% CL exclusion contours for the combined fit of SR-lowMass and SR-highMass for simplified models with combined ${\tilde{\tau}}^{+}_{R,L} {\tilde{\tau}}^{-}_{R,L}$ production.

The expected 95% CL exclusion contours for the combined fit of SR-lowMass and SR-highMass for simplified models with combined ${\tilde{\tau}}^{+}_{R,L} {\tilde{\tau}}^{-}_{R,L}$ production.

The expected 95% CL exclusion contours for the combined fit of SR-lowMass and SR-highMass for simplified models with combined ${\tilde{\tau}}^{+}_{R,L} {\tilde{\tau}}^{-}_{R,L}$ production.

The observed 95% CL exclusion contours for the combined fit of SR-lowMass and SR-highMass for simplified models with ${\tilde{\tau}}_L {\tilde{\tau}}_L$ only production.

The observed 95% CL exclusion contours for the combined fit of SR-lowMass and SR-highMass for simplified models with ${\tilde{\tau}}_L {\tilde{\tau}}_L$ only production.

The expected 95% CL exclusion contours for the combined fit of SR-lowMass and SR-highMass for simplified models with ${\tilde{\tau}}_L {\tilde{\tau}}_L$ only production.

The expected 95% CL exclusion contours for the combined fit of SR-lowMass and SR-highMass for simplified models with ${\tilde{\tau}}_L {\tilde{\tau}}_L$ only production.

Observed 95% CL exclusion limits for simplified models with direct stau pair production in SR-lowMass.

Observed 95% CL exclusion limits for simplified models with direct stau pair production in SR-lowMass.

Expected 95% CL exclusion limits for simplified models with direct stau pair production in SR-lowMass.

Expected 95% CL exclusion limits for simplified models with direct stau pair production in SR-lowMass.

Observed 95% CL exclusion limits for simplified models with direct stau pair production in SR-highMass.

Observed 95% CL exclusion limits for simplified models with direct stau pair production in SR-highMass.

Expected 95% CL exclusion limits for simplified models with direct stau pair production in SR-highMass.

Expected 95% CL exclusion limits for simplified models with direct stau pair production in SR-highMass.

Signal acceptance in SR highMass for combined stau final states

Signal acceptance in SR highMass for combined stau final states

Signal acceptance in SR lowMass for combined stau final states

Signal acceptance in SR lowMass for combined stau final states

Signal efficiency in SR highMass for combined stau final states

Signal efficiency in SR highMass for combined stau final states

Signal efficiency in SR lowMass for combined stau final states

Signal efficiency in SR lowMass for combined stau final states

Signal acceptance*efficiency in SR highMass for combined stau final states

Signal acceptance*efficiency in SR highMass for combined stau final states

Signal acceptance*efficiency in SR lowMass for combined stau final states

Signal acceptance*efficiency in SR lowMass for combined stau final states

Cutflow for two reference points (${\tilde{\tau}}^{+}_{R,L} {\tilde{\tau}}^{-}_{R,L}$ production) in SR. The column labelled $N_{weighted}$ shows the results including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$, while $N_{raw}$ in brackets shows the results for the generated number of events. The quoted uncertainties are statistical only. The "Generator filter" includes the requirements that two $\tau$ in the event have ${p}_{T} > 15$ GeV and $|\eta| <$ 2.6. The "Baseline Cut" includes the requirement of two baseline $\tau$ with a minimum value at 0.01 of the boosted decision tree discriminant (JetBDTSigTransMin $>$ 0.01) and ${p}_{T, \tau_{1}} > 50$ GeV and ${p}_{T, \tau_{2}} > 40$ GeV. At the step "Trigger & offline cuts", the following requirements are applied: the event is recorded using the asymmetric di-$\tau$ trigger (di-$\tau$ $E_{T}^{miss}$ trigger) in SR-lowMass (SR-highMass), and the lepton $p_{T}$ and $E_{T}^{miss}$ are required at plateau.

Cutflow for two reference points (${\tilde{\tau}}^{+}_{R,L} {\tilde{\tau}}^{-}_{R,L}$ production) in SR. The column labelled $N_{weighted}$ shows the results including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$, while $N_{raw}$ in brackets shows the results for the generated number of events. The quoted uncertainties are statistical only. The "Generator filter" includes the requirements that two $\tau$ in the event have ${p}_{T} > 15$ GeV and $|\eta| <$ 2.6. The "Baseline Cut" includes the requirement of two baseline $\tau$ with a minimum value at 0.01 of the boosted decision tree discriminant (JetBDTSigTransMin $>$ 0.01) and ${p}_{T, \tau_{1}} > 50$ GeV and ${p}_{T, \tau_{2}} > 40$ GeV. At the step "Trigger & offline cuts", the following requirements are applied: the event is recorded using the asymmetric di-$\tau$ trigger (di-$\tau$ $E_{T}^{miss}$ trigger) in SR-lowMass (SR-highMass), and the lepton $p_{T}$ and $E_{T}^{miss}$ are required at plateau.

Observed and expected numbers of events in the control and signal regions where all control and signal region bins are included as constraints in the likelihood. The expected event yields of SM processes are given after the background-only fit. The entries marked as "--" are negligible. The uncertainties correspond to the sum in quadrature of statistical and systematic uncertainties. The correlation of systematic uncertainties among control regions and among background processes is fully taken into account.

Observed and expected numbers of events in the control and signal regions where all control and signal region bins are included as constraints in the likelihood. The expected event yields of SM processes are given after the background-only fit. The entries marked as "--" are negligible. The uncertainties correspond to the sum in quadrature of statistical and systematic uncertainties. The correlation of systematic uncertainties among control regions and among background processes is fully taken into account.

The post-fit $m_{T2}$ distribution for SR-lowMass. The stacked histograms show the expected SM backgrounds. The multi-jet contribution is estimated from data using the ABCD method. The contributions of multi-jet and $W$+jets events are scaled with the corresponding normalization factors derived from the background-only fit. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distributions from the SUSY reference points are also shown as dashed lines. The last bin includes the overflow events.

The post-fit $m_{T2}$ distribution for SR-lowMass. The stacked histograms show the expected SM backgrounds. The multi-jet contribution is estimated from data using the ABCD method. The contributions of multi-jet and $W$+jets events are scaled with the corresponding normalization factors derived from the background-only fit. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distributions from the SUSY reference points are also shown as dashed lines. The last bin includes the overflow events.

The post-fit $m_{T2}$ distribution for SR-highMass. The stacked histograms show the expected SM backgrounds. The multi-jet contribution is estimated from data using the ABCD method. The contributions of multi-jet and $W$+jets events are scaled with the corresponding normalization factors derived from the background-only fit. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distributions from the SUSY reference points are also shown as dashed lines. The last bin includes the overflow events.

The post-fit $m_{T2}$ distribution for SR-highMass. The stacked histograms show the expected SM backgrounds. The multi-jet contribution is estimated from data using the ABCD method. The contributions of multi-jet and $W$+jets events are scaled with the corresponding normalization factors derived from the background-only fit. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distributions from the SUSY reference points are also shown as dashed lines. The last bin includes the overflow events.

The $m_{T2}$ post-fit distributions in the multi-jet background validation region VR-F (lowMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The $m_{T2}$ post-fit distributions in the multi-jet background validation region VR-F (lowMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The $m_{T2}$ post-fit distributions in the multi-jet background validation VR-F (highMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The $m_{T2}$ post-fit distributions in the multi-jet background validation VR-F (highMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The $E_{T}^{miss}$ post-fit distributions in the multi-jet background validation region VR-F (lowMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The $E_{T}^{miss}$ post-fit distributions in the multi-jet background validation region VR-F (lowMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The $E_{T}^{miss}$ post-fit distributions in the multi-jet background validation region VR-F (highMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The $E_{T}^{miss}$ post-fit distributions in the multi-jet background validation region VR-F (highMass). The stacked histograms show the contribution of each relevant SM process. The multi-jet shape is taken from VR-E in the ABCD method and the normalization is determined by the transfer factor $T$ and rescaled by a correction factor determined by the fit. The hatched bands represent the combined statistical and systematic uncertainties in the sum of the SM backgrounds shown. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The pre-fit $m_{T2}$ distribution in the $WCR$. The SM backgrounds other than multi-jet production are estimated from MC simulation. The multi-jet contribution is estimated from data using the OS--SS method. The hatched bands represent the combined statistical and systematic uncertainties of the total SM background. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The pre-fit $m_{T2}$ distribution in the $WCR$. The SM backgrounds other than multi-jet production are estimated from MC simulation. The multi-jet contribution is estimated from data using the OS--SS method. The hatched bands represent the combined statistical and systematic uncertainties of the total SM background. For illustration, the distributions of the SUSY reference points are also shown as dashed lines.

The post-fit yields in the $WVR$, $TVRs$, $ZVRs$ and $VVVRs$. The SM backgrounds other than multi-jet production are estimated from MC simulation. The multi-jet contribution is negligible and is estimated from data using the ABCD method, using CRs obtained with the same technique used for the SRs. The hatched bands represent the combined statistical and systematic uncertainties of the total SM background. For illustration, the distributions of the SUSY reference points are also shown as dashed lines. The lower panels show the ratio of data to the SM background estimate.

The post-fit yields in the $WVR$, $TVRs$, $ZVRs$ and $VVVRs$. The SM backgrounds other than multi-jet production are estimated from MC simulation. The multi-jet contribution is negligible and is estimated from data using the ABCD method, using CRs obtained with the same technique used for the SRs. The hatched bands represent the combined statistical and systematic uncertainties of the total SM background. For illustration, the distributions of the SUSY reference points are also shown as dashed lines. The lower panels show the ratio of data to the SM background estimate.

Observed and predicted mTtot distribution in the b-veto category of the 1l1tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table.The last bin includes overflows. The combined prediction for A and H bosons with masses of 400, 1000 and 1500 GeV and $\tan\beta$ = 6, 12 and 25 respectively in the mh125 scenario are also provided. The combined prediction for A and H bosons with masses of 1000 and 1500 GeV is scaled by 100 in the paper figure, but not in the HepData table.

Observed and predicted mTtot distribution in the b-tag category of the 1l1tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 400, 1000 and 1500 GeV and $\tan\beta$ = 6, 12 and 25 respectively in the mh125 scenario are also provided. The combined prediction for A and H bosons with masses of 1000 and 1500 GeV is scaled by 100 in the paper figure, but not in the HepData table.

Observed and predicted mTtot distribution in the b-tag category of the 1l1tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 400, 1000 and 1500 GeV and $\tan\beta$ = 6, 12 and 25 respectively in the mh125 scenario are also provided. The combined prediction for A and H bosons with masses of 1000 and 1500 GeV is scaled by 100 in the paper figure, but not in the HepData table.

Observed and predicted mTtot distribution in the b-tag category of the 1l1tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 400, 1000 and 1500 GeV and $\tan\beta$ = 6, 12 and 25 respectively in the mh125 scenario are also provided. The combined prediction for A and H bosons with masses of 1000 and 1500 GeV is scaled by 100 in the paper figure, but not in the HepData table.

Observed and predicted mTtot distribution in the b-tag category of the 1l1tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 400, 1000 and 1500 GeV and $\tan\beta$ = 6, 12 and 25 respectively in the mh125 scenario are also provided. The combined prediction for A and H bosons with masses of 1000 and 1500 GeV is scaled by 100 in the paper figure, but not in the HepData table.

Observed and predicted mTtot distribution in the b-veto category of the 2tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 400, 1000 and 1500 GeV and $\tan\beta$ = 6, 12 and 25 respectively in the mh125 scenario are also provided. The combined prediction for A and H bosons with masses of 1000 and 1500 GeV is scaled by 100 in the paper figure, but not in the HepData table.

Observed and predicted mTtot distribution in the b-veto category of the 2tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 400, 1000 and 1500 GeV and $\tan\beta$ = 6, 12 and 25 respectively in the mh125 scenario are also provided. The combined prediction for A and H bosons with masses of 1000 and 1500 GeV is scaled by 100 in the paper figure, but not in the HepData table.

Observed and predicted mTtot distribution in the b-veto category of the 2tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 400, 1000 and 1500 GeV and $\tan\beta$ = 6, 12 and 25 respectively in the mh125 scenario are also provided. The combined prediction for A and H bosons with masses of 1000 and 1500 GeV is scaled by 100 in the paper figure, but not in the HepData table.

Observed and predicted mTtot distribution in the b-veto category of the 2tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 400, 1000 and 1500 GeV and $\tan\beta$ = 6, 12 and 25 respectively in the mh125 scenario are also provided. The combined prediction for A and H bosons with masses of 1000 and 1500 GeV is scaled by 100 in the paper figure, but not in the HepData table.

Observed and predicted mTtot distribution in the b-tag category of the 2tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 400, 1000 and 1500 GeV and $\tan\beta$ = 6, 12 and 25 respectively in the mh125 scenario are also provided. The combined prediction for A and H bosons with masses of 1000 and 1500 GeV is scaled by 100 in the paper figure, but not in the HepData table.

Observed and predicted mTtot distribution in the b-tag category of the 2tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 400, 1000 and 1500 GeV and $\tan\beta$ = 6, 12 and 25 respectively in the mh125 scenario are also provided. The combined prediction for A and H bosons with masses of 1000 and 1500 GeV is scaled by 100 in the paper figure, but not in the HepData table.

Observed and predicted mTtot distribution in the b-tag category of the 2tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 400, 1000 and 1500 GeV and $\tan\beta$ = 6, 12 and 25 respectively in the mh125 scenario are also provided. The combined prediction for A and H bosons with masses of 1000 and 1500 GeV is scaled by 100 in the paper figure, but not in the HepData table.

Observed and predicted mTtot distribution in the b-tag category of the 2tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 400, 1000 and 1500 GeV and $\tan\beta$ = 6, 12 and 25 respectively in the mh125 scenario are also provided. The combined prediction for A and H bosons with masses of 1000 and 1500 GeV is scaled by 100 in the paper figure, but not in the HepData table.

Observed and expected 95% CL upper limits on the gluon-gluon fusion Higgs boson production cross section times ditau branching fraction as a function of the Higgs boson mass.

Observed and expected 95% CL upper limits on the gluon-gluon fusion Higgs boson production cross section times ditau branching fraction as a function of the Higgs boson mass.

Observed and expected 95% CL upper limits on the gluon-gluon fusion Higgs boson production cross section times ditau branching fraction as a function of the Higgs boson mass.

Observed and expected 95% CL upper limits on the gluon-gluon fusion Higgs boson production cross section times ditau branching fraction as a function of the Higgs boson mass.

Observed and expected 95% CL upper limits on the b-associated Higgs boson production cross section times ditau branching fraction as a function of the boson mass.

Observed and expected 95% CL upper limits on the b-associated Higgs boson production cross section times ditau branching fraction as a function of the boson mass.

Observed and expected 95% CL upper limits on the b-associated Higgs boson production cross section times ditau branching fraction as a function of the boson mass.

Observed and expected 95% CL upper limits on the b-associated Higgs boson production cross section times ditau branching fraction as a function of the boson mass.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. No theoretical uncertainty is considered when computing these limits.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

Acceptance times efficiency for a scalar boson produced by gluon-gluon fusion as a function of the scalar boson mass.

Acceptance times efficiency for a scalar boson produced by gluon-gluon fusion as a function of the scalar boson mass.

Acceptance times efficiency for a scalar boson produced by gluon-gluon fusion as a function of the scalar boson mass.

Acceptance times efficiency for a scalar boson produced by gluon-gluon fusion as a function of the scalar boson mass.

Acceptance times efficiency for a scalar boson produced by b-associated production as a function of the scalar boson mass.

Acceptance times efficiency for a scalar boson produced by b-associated production as a function of the scalar boson mass.

Acceptance times efficiency for a scalar boson produced by b-associated production as a function of the scalar boson mass.

Acceptance times efficiency for a scalar boson produced by b-associated production as a function of the scalar boson mass.

Observed 95% CL upper limits on the scalar boson production cross section times ditau branching fraction as a function of the scalar boson mass and the fraction of the b-associated production. The limits are calculated from a statistical combination of the 1l1tau_h and 2tau_h channels.

Observed 95% CL upper limits on the scalar boson production cross section times ditau branching fraction as a function of the scalar boson mass and the fraction of the b-associated production. The limits are calculated from a statistical combination of the 1l1tau_h and 2tau_h channels.

Observed 95% CL upper limits on the scalar boson production cross section times ditau branching fraction as a function of the scalar boson mass and the fraction of the b-associated production. The limits are calculated from a statistical combination of the 1l1tau_h and 2tau_h channels.

Observed 95% CL upper limits on the scalar boson production cross section times ditau branching fraction as a function of the scalar boson mass and the fraction of the b-associated production. The limits are calculated from a statistical combination of the 1l1tau_h and 2tau_h channels.

Expected 95% CL upper limits on the scalar boson production cross section times ditau branching fraction as a function of the scalar boson mass and the fraction of the b-associated production. The limits are calculated from a statistical combination of the 1l1tau_h and 2tau_h channels.

Expected 95% CL upper limits on the scalar boson production cross section times ditau branching fraction as a function of the scalar boson mass and the fraction of the b-associated production. The limits are calculated from a statistical combination of the 1l1tau_h and 2tau_h channels.

Expected 95% CL upper limits on the scalar boson production cross section times ditau branching fraction as a function of the scalar boson mass and the fraction of the b-associated production. The limits are calculated from a statistical combination of the 1l1tau_h and 2tau_h channels.

Expected 95% CL upper limits on the scalar boson production cross section times ditau branching fraction as a function of the scalar boson mass and the fraction of the b-associated production. The limits are calculated from a statistical combination of the 1l1tau_h and 2tau_h channels.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 200 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 200 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 200 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 200 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 250 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 250 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 250 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 250 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 300 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 300 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 300 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 300 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 350 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 350 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 350 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 350 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 400 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 400 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 400 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 400 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 500 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 500 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 500 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 500 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 600 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 600 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 600 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 600 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 700 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 700 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 700 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 700 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 800 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 800 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 800 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 800 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1000 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1000 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1000 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1000 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1200 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1200 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1200 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1200 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1500 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1500 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1500 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1500 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2000 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2000 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2000 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2000 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2500 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2500 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2500 GeV signal mass point is shown in the HEPData table.

Observed two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2500 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 200 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 200 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 200 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 250 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 250 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 250 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 300 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 300 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 300 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 350 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 350 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 350 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 400 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 400 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 400 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 500 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 500 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 500 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 600 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 600 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 600 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 700 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 700 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 700 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 800 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 800 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 800 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1000 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1000 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1000 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1200 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1200 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1200 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1500 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1500 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1500 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2000 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2000 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2000 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2500 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2500 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2500 GeV signal mass point is shown in the HEPData table.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\chi})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\chi})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\chi})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\chi})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\chi})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\chi})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\chi})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\chi})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\chi})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\chi})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\chi})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\chi})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\tau})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\tau})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\tau})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\tau})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\tau})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\tau})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\tau})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\tau})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\tau})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\tau})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(\widetilde{\tau})$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(\widetilde{\tau})$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60. The theoretical uncertainty of signal cross section is considered.

The observed 95% CL upper limits with one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The theoretical uncertainty of signal cross section is considered.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}(alignment)$ scenario. The lowest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 1.0. The highest value of $\tan\beta$ considered by the $M_{h}^{125}(alignment)$ scenario is 20.0. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The theoretical uncertainty of signal cross section is considered.

Three-body selection. Distributions of $M_{\Delta}^R$ in (a,b) $SR_{W}^{3-body}$ and (c,d) $SR_{T}^{3-body}$ for (left) same-flavour and (right) different-flavour events satisfying the selection criteria of the given SR, except the one for the presented variable, after the background fit. The contributions from all SM backgrounds are shown as a histogram stack. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. The hatched bands represent the total statistical and systematic uncertainty. The rightmost bin of each plot includes overflow events. Reference top squark pair production signal models are overlayed for comparison. Red arrows in the upper panels indicate the signal region selection criteria. The bottom panels show the ratio of the observed data to the total SM background prediction, with hatched bands representing the total uncertainty in the background prediction; red arrows show data outside the vertical-axis range.

Three-body selection. Distributions of $M_{\Delta}^R$ in (a,b) $SR_{W}^{3-body}$ and (c,d) $SR_{T}^{3-body}$ for (left) same-flavour and (right) different-flavour events satisfying the selection criteria of the given SR, except the one for the presented variable, after the background fit. The contributions from all SM backgrounds are shown as a histogram stack. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. The hatched bands represent the total statistical and systematic uncertainty. The rightmost bin of each plot includes overflow events. Reference top squark pair production signal models are overlayed for comparison. Red arrows in the upper panels indicate the signal region selection criteria. The bottom panels show the ratio of the observed data to the total SM background prediction, with hatched bands representing the total uncertainty in the background prediction; red arrows show data outside the vertical-axis range.

Three-body selection. Distributions of $M_{\Delta}^R$ in (a,b) $SR_{W}^{3-body}$ and (c,d) $SR_{T}^{3-body}$ for (left) same-flavour and (right) different-flavour events satisfying the selection criteria of the given SR, except the one for the presented variable, after the background fit. The contributions from all SM backgrounds are shown as a histogram stack. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. The hatched bands represent the total statistical and systematic uncertainty. The rightmost bin of each plot includes overflow events. Reference top squark pair production signal models are overlayed for comparison. Red arrows in the upper panels indicate the signal region selection criteria. The bottom panels show the ratio of the observed data to the total SM background prediction, with hatched bands representing the total uncertainty in the background prediction; red arrows show data outside the vertical-axis range.

Four-body selection. (a) distributions of $E_{T}^{miss}$ in $SR^{4-body}_{Small\,\Delta m}$ and (b) distribution of $R_{2\ell 4j}$ in $SR^{4-body}_{Large\,\Delta m}$ for events satisfying the selection criteria of the given SR, except the one for the presented variable, after the background fit. The contributions from all SM backgrounds are shown as a histogram stack. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. The hatched bands represent the total statistical and systematic uncertainty. The rightmost bin of each plot includes overflow events. Reference top squark pair production signal models are overlayed for comparison. Red arrows in the upper panel indicate the signal region selection criteria. The bottom panels show the ratio of the observed data to the total SM background prediction, with hatched bands representing the total uncertainty in the background prediction; red arrows show data outside the vertical-axis range.

Four-body selection. (a) distributions of $E_{T}^{miss}$ in $SR^{4-body}_{Small\,\Delta m}$ and (b) distribution of $R_{2\ell 4j}$ in $SR^{4-body}_{Large\,\Delta m}$ for events satisfying the selection criteria of the given SR, except the one for the presented variable, after the background fit. The contributions from all SM backgrounds are shown as a histogram stack. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. The hatched bands represent the total statistical and systematic uncertainty. The rightmost bin of each plot includes overflow events. Reference top squark pair production signal models are overlayed for comparison. Red arrows in the upper panel indicate the signal region selection criteria. The bottom panels show the ratio of the observed data to the total SM background prediction, with hatched bands representing the total uncertainty in the background prediction; red arrows show data outside the vertical-axis range.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the Observed limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100\% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100\% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100\% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100\% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100\% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}_1^0$ with 100\% branching ratio, in the (a) $m(\tilde{t}_1)$--$m(\tilde{\chi}_1^0)$ and (b) $m(\tilde{t}_1)$--$\Delta m(\tilde{t}_1,\tilde{\chi}_1^0)$ planes. The dashed lines and the shaded bands are the expected limits and their $\pm1\sigma$ uncertainties. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the mediator mass for a DM particle mass of $m(\chi)=1$ GeV. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross-section to the nominal cross-section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines shows the observed (expected) exclusion limits.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the mediator mass for a DM particle mass of $m(\chi)=1$ GeV. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross-section to the nominal cross-section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines shows the observed (expected) exclusion limits.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the mediator mass for a DM particle mass of $m(\chi)=1$ GeV. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross-section to the nominal cross-section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines shows the observed (expected) exclusion limits.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the mediator mass for a DM particle mass of $m(\chi)=1$ GeV. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross-section to the nominal cross-section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines shows the observed (expected) exclusion limits.

Two-body selection. Background fit results for $\mathrm{CR}^{\mathrm{2-body}}_{t\bar{t}}$, $\mathrm{CR}^{\mathrm{2-body}}_{t\bar{t}Z}$, $\mathrm{VR}^{\mathrm{2-body}}_{t\bar{t}, DF}$, $\mathrm{VR}^{\mathrm{2-body}}_{t\bar{t}, SF}$ and $\mathrm{VR}^{\mathrm{2-body}}_{t\bar{t} Z}$. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. Combined statistical and systematic uncertainties are given. Entries marked `--' indicate a negligible background contribution (less than 0.001 events). The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Three-body selection. Background fit results for $\mathrm{CR}^{\mathrm{3-body}}_{t\bar{t}}$, $\mathrm{CR}^{\mathrm{3-body}}_{VV}$, $\mathrm{CR}^{\mathrm{2-body}}_{t\bar{t}Z}$, $\mathrm{VR}^{\mathrm{3-body}}_{VV}$, $\mathrm{VR(1)}^{\mathrm{3-body}}_{t\bar{t}}$ and $\mathrm{VR(2)}^{\mathrm{3-body}}_{t\bar{t}}$. ''Others'' includes contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$ processes. Combined statistical and systematic uncertainties are given. Entries marked `--' indicate a negligible background contribution (less than 0.001 events). The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Four-body selection. Background fit results for $\mathrm{CR}^{\mathrm{4-body}}_{t\bar{t}}$,$\mathrm{CR}^{\mathrm{4-body}}_{VV}$, $\mathrm{VR}^{\mathrm{4-body}}_{t\bar{t}}$, $VR^{4-body}_{VV}$ and $\mathrm{VR}^{\mathrm{4-body}}_{VV,lll}$. The ''Others'' category contains the contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$. Combined statistical and systematic uncertainties are given. Entries marked `--' indicate a negligible background contribution (less than 0.001 events). The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Two-body selection. Background fit results for the different-flavour leptons binned SRs. The ''Others'' category contains the contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$. Combined statistical and systematic uncertainties are given. Entries marked `--' indicate a negligible background contribution (less than 0.001 events). The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Two-body selection. Background fit results for the same-flavour leptons binned SRs. The ''Others'' category contains the contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$. Combined statistical and systematic uncertainties are given. The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Three-body selection. Observed event yields and background fit results for the three-body selection SRs. The ''Others'' category contains contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$. Combined statistical and systematic uncertainties are given. Entries marked `--' indicate a negligible background contribution (less than 0.001 events). The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Four-body selection. Observed event yields and background fit results for SR$^{\mathrm{4-body}}_{\mathrm{Small}\,\Delta m}$ and SR$^{\mathrm{4-body}}_{\mathrm{Large}\,\Delta m}$. The ''Others'' category contains the contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$. Combined statistical and systematic uncertainties are given. The individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Exclusion limits contours (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}^0_1$ with 100% branching ratio in $\tilde{t}_1--\tilde{\chi}^0_1$ masses planes. The dashed lines and the shaded bands are the expected limit and its $\pm 1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The exclusion limits contours for the two-body, three-body and four-body selections are respectively shown in blue, green and red.

Exclusion limits contours (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}^0_1$ with 100% branching ratio in $\tilde{t}_1--\tilde{\chi}^0_1$ masses planes. The dashed lines and the shaded bands are the expected limit and its $\pm 1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The exclusion limits contours for the two-body, three-body and four-body selections are respectively shown in blue, green and red.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow t \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b W \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b W \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b W \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm 1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b W \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty. The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty.The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limit contour (95% CL) for a simplified model assuming $\tilde{t}_1$ pair production, decaying via $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}_1^0$ with 100% branching ratio, in $\tilde{t}_1$--$\tilde{\chi}_1^0$ masses plane. The dashed lines and the shaded bands are the expected limit and its $\pm1\sigma$ uncertainty.The thick solid lines are the observed limits for the central value of the signal cross-section. The expected and observed limits do not include the effect of the theoretical uncertainties in the signal cross-section. The dotted lines show the effect on the observed limit when varying the signal cross-section by $\pm1\sigma$ of the theoretical uncertainty. The observed (a) and expected (b) CLs values are respectively shown.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the DM particle mass for a mediator mass of 10 GeV. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross-section to the nominal cross-section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines shows the observed (expected) exclusion limits.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the DM particle mass for a mediator mass of 10 GeV. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross-section to the nominal cross-section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines shows the observed (expected) exclusion limits.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the DM particle mass for a mediator mass of 10 GeV. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross-section to the nominal cross-section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines shows the observed (expected) exclusion limits.

Exclusion limits for (a) $t\bar{t} + \phi $ scalar and (b) $t\bar{t} + a $ pseudoscalar models as a function of the DM particle mass for a mediator mass of 10 GeV. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross-section to the nominal cross-section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines shows the observed (expected) exclusion limits.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection efficiency (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Three-body selection efficiency (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection efficiency (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection efficiency (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection efficiency (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Four-body selection Efficiency (a) SR$^{4-body}_{Small \Delta m}$ , (b) $SR^{4-body}_{Large \Delta m}$ for a simplified model assuming $\tilde{t}_1$ pair production.

Four-body selection Efficiency (a) SR$^{4-body}_{Small \Delta m}$ , (b) $SR^{4-body}_{Large \Delta\ m}$ for a simplified model assuming $\tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} +\phi$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ t \tilde{t} +\phi$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ t \tilde{t} +\phi$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ t \tilde{t} +\phi$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ t \tilde{t} +\phi$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $ t \tilde{t} +\phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + \phi$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-DF$^{2-body}_{[110,120)}$, (b) SR-DF1$^{2-body}_{[120,140)}$, (c) SR-DF2$^{2-body}_{[140,160)}$, (d) SR-DF3$^{2-body}_{[160,180)}$, (e) SR-DF4$^{2-body}_{[180,220)}$, (f) SR-DF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) SR-SF$^{2-body}_{[110,120)}$, (b) SR-SF1$^{2-body}_{[120,140)}$, (c) SR-SF2$^{2-body}_{[140,160)}$, (d) SR-SF3$^{2-body}_{[160,180)}$, (e) SR-SF4$^{2-body}_{[180,220)}$, (f) SR-SF5$^{2-body}_{[220,\infty)}$ for a simplified model assuming $t \tilde{t} + a$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Two-body selection acceptance (a) $SR^{2-body}_{[110,\infty)}$ , (b) $SR^{2-body}_{[120,\infty)}$ , (c) $SR^{2-body}_{[140,\infty)}$ , (d) $SR^{2-body}_{[160,\infty)}$ , (e) $SR^{2-body}_{[180,\infty)}$ , (f) $SR^{2-body}_{[200,\infty)}$ , (g) $SR^{2-body}_{[220,\infty)}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection acceptance (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection acceptance (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection acceptance (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Three-body selection acceptance (a) SR-DF$^{3-body}_{t}$, (b) SR-SF$^{3-body}_{t}$, (c) SR-DF$^{3-body}_{W}$, (d) SR-SF$^{3-body}_{W}$ for a simplified model assuming $ \tilde{t}_1$ pair production.

Four-body selection acceptance (a) SR$^{4-body}_{Small \Delta m}$ , (b) $SR^{4-body}_{Large \Delta m}$ for a simplified model assuming $\tilde{t}_1$ pair production.

Four-body selection acceptance (a) SR$^{4-body}_{Small \Delta m}$ , (b) $SR^{4-body}_{Large \Delta m}$ for a simplified model assuming $\tilde{t}_1$ pair production.

Two-body selection The numbers indicate the observed upper limits on the signal strenght for (a) a simplified model assuming $\tilde{t}_1$ pair production, (b) for $t\tilde{t} + a $ pseudoscalar models, (c) for $t\tilde{t} + \phi $ scalar models. In Figure (a), the red line corresponds to the observed limit.

Two-body selection The numbers indicate the observed upper limits on the signal strenght for (a) a simplified model assuming $\tilde{t}_1$ pair production, (b) for $t\tilde{t} + a $ pseudoscalar models, (c) for $t\tilde{t} + \phi $ scalar models. In Figure (a), the red line corresponds to the observed limit.

Two-body selection The numbers indicate the observed upper limits on the signal strenght for (a) a simplified model assuming $\tilde{t}_1$ pair production, (b) for $t\tilde{t} + a $ pseudoscalar models, (c) for $t\tilde{t} + \phi $ scalar models. In Figure (a), the red line corresponds to the observed limit.

Three-body selection The numbers indicate the upper limits on the signal strenght for a simplified model assuming $\tilde{t}_1$ pair production. For comparison, the red line corresponds to the observed limit.

Four-body selection The numbers indicate the upper limits on the signal strenght for a simplified model assuming $\tilde{t}_1$ pair production. For comparison, the red line corresponds to the observed limit.

Two-body selection The numbers indicate the upper limits on the signal cross-section for (a) a simplified model assuming $\tilde{t}_1$ pair production, (b) for $t\tilde{t} + a $ pseudoscalar models, (c) for $t\tilde{t} + \phi $ scalar models. In Figure (a), the red line corresponds to the observed limit.

Two-body selection The numbers indicate the upper limits on the signal cross-section for (a) a simplified model assuming $\tilde{t}_1$ pair production, (b) for $t\tilde{t} + a $ pseudoscalar models, (c) for $t\tilde{t} + \phi $ scalar models. In Figure (a), the red line corresponds to the observed limit.

Two-body selection The numbers indicate the upper limits on the signal cross-section for (a) a simplified model assuming $\tilde{t}_1$ pair production, (b) for $t\tilde{t} + a $ pseudoscalar models, (c) for $t\tilde{t} + \phi $ scalar models. In Figure (a), the red line corresponds to the observed limit.

Three-body selection The numbers indicate the upper limits on the signal cross-section for a simplified model assuming $\tilde{t}_1$ pair production. For comparison, the red line corresponds to the observed limit.

Four-body selection The numbers indicate the upper limits on the signal cross-section for a simplified model assuming $\tilde{t}_1$ pair production. For comparison, the red line corresponds to the observed limit.

Two-body selection. Background fit results for the $inclusive$ SRs. The Others category contains the contributions from $VVV$, $t\bar{t} t$, $t\bar{t}t\bar{t}$, $t\bar{t} W$, $t\bar{t} WW$, $t\bar{t} WZ$, $t\bar{t} H$, and $tZ$. Combined statistical and systematic uncertainties are given. Note that the individual uncertainties can be correlated, and do not necessarily add up quadratically to the total background uncertainty.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow t^{(*)}\tilde{\chi}^0_1$ with $m(\tilde{t}_1)=600~ GeV$ and $m(\tilde{\chi}^0_1)=400~ GeV$ in the SRs for the two-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the scalar signal model $t\bar{t} + \phi $ with $m(\phi)=150~ GeV$ and $m(\chi)=1~ GeV$ in the SRs for the two-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the pseudoscalar signal model $t\bar{t} + a $ with $m(a)=150~ GeV$ and $m(\chi)=1~ GeV$ in the SRs for the two-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow bW\tilde{\chi}^0_1$ with $m(\tilde{t}_1)=550~ GeV$ and $m(\tilde{\chi}^0_1)=385~ GeV$ in the SRs for the three-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow bW\tilde{\chi}^0_1$ with $m(\tilde{t}_1)=550~ GeV$ and $m(\tilde{\chi}^0_1)=400~ GeV$ in the SRs for the three-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow bW\tilde{\chi}^0_1$ with $m(\tilde{t}_1)=550~ GeV$ and $m(\tilde{\chi}^0_1)=430~ GeV$ in the SRs for the three-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow bW\tilde{\chi}^0_1$ with $m(\tilde{t}_1)=550~ GeV$ and $m(\tilde{\chi}^0_1)=460~ GeV$ in the SRs for the three-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}^0_1$ with $m(\tilde{t}_1)=400~ GeV$ and $m(\tilde{\chi}^0_1)=380~ GeV$ in the SRs for the four-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}^0_1$ with $m(\tilde{t}_1)=460~ GeV$ and $m(\tilde{\chi}^0_1)=415~ GeV$ in the SRs for the four-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Cut flow for the simplified signal model $\tilde{t}_1 \rightarrow b l \nu \tilde{\chi}^0_1$ with $m(\tilde{t}_1)=400~ GeV$ and $m(\tilde{\chi}^0_1)=320~ GeV$ in the SRs for the four-body selection. The number of events is normalized to the cross-section and to an integrated luminosity of $139~fb^{-1}$.

Summary of the observed 95% CL upper limits on $\sigma_{ZH} BR(H\rightarrow aa \rightarrow (b\bar{b})(b\bar{b}))$ for the resolved analysis.

Efficiency and acceptance for simulated $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ samples in two signal regions (SR) of the analysis, one with two $a\to b\bar{b}$ candidates in the High Purity Category (HPC), and the other with one $a\to b\bar{b}$ candidate in the High Purity Category (HPC) and one in the Low Purity Category (LPC).

Efficiency and acceptance for simulated $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ samples in two signal regions (SR) of the analysis, one with two $a\to b\bar{b}$ candidates in the High Purity Category (HPC), and the other with one $a\to b\bar{b}$ candidate in the High Purity Category (HPC) and one in the Low Purity Category (LPC).

Event yields for a simulated $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ sample with $m_a = 17.5\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$. Cut 0 corresponds to the initial number of events. Cut 1 requires the single lepton trigger. Cut 2 requires 2 identified leptons. Cut 3 requires the Z-boson mass window. Cut 4 requires 2 reconstructed $a\to b\bar{b}$ candidates. Cut 5a requires 2 identified $a\to b\bar{b}$ candidates in the 1HPC1LPC region. Cut 6a requires the 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region to be inside the Higgs mass window. Cut 5b requires 2 identified $a\to b\bar{b}$ candidates in the 2HPC region. Cut 6b requires the 2 $a\to b\bar{b}$ candidates in the 2HPC region to be inside the Higgs mass window.

Event yields for a simulated $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ sample with $m_a = 17.5\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$. Cut 0 corresponds to the initial number of events. Cut 1 requires the single lepton trigger. Cut 2 requires 2 identified leptons. Cut 3 requires the Z-boson mass window. Cut 4 requires 2 reconstructed $a\to b\bar{b}$ candidates. Cut 5a requires 2 identified $a\to b\bar{b}$ candidates in the 1HPC1LPC region. Cut 6a requires the 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region to be inside the Higgs mass window. Cut 5b requires 2 identified $a\to b\bar{b}$ candidates in the 2HPC region. Cut 6b requires the 2 $a\to b\bar{b}$ candidates in the 2HPC region to be inside the Higgs mass window.

Background yield table for Z+jets, $t\bar{t}$, and rare sources. Observed data yield. Signal $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ yield with $m_a = 20\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$. The table includes the yields in two signal regions with leptons consistent with an on-shell Z-boson decay, one with 2 $a\to b\bar{b}$ candidates in the 2HPC region and one with 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region. The table also includes the yields in four control regions, one with leptons consistent with an on-shell Z-boson decay and 2 $a\to b\bar{b}$ candidates in the Low Purity Category (LPC), and three others where the leptons are not consistent an on-shell Z-boson decay.

Background yield table for Z+jets, $t\bar{t}$, and rare sources. Observed data yield. Signal $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ yield with $m_a = 20\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$, with a branching ratio set to 1 for the $H \rightarrow aa$ decay, whereas the ATLAS figure attached to this entry instead uses the upper-limit branching ratio (smaller than 1). The table includes the yields in two signal regions with leptons consistent with an on-shell Z-boson decay, one with 2 $a\to b\bar{b}$ candidates in the 2HPC region and one with 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region. The table also includes the yields in four control regions, one with leptons consistent with an on-shell Z-boson decay and 2 $a\to b\bar{b}$ candidates in the Low Purity Category (LPC), and three others where the leptons are not consistent an on-shell Z-boson decay.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) × B(s&rarr; 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 &sigma;(pp &rarr; s &chi;&chi;) × B(s&rarr; VV) process for m_Z'=0.5 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) × B(s&rarr; 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 &sigma;(pp &rarr; s &chi;&chi;) × B(s&rarr; VV) process for m_Z'=1 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) × B(s&rarr; 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 &sigma;(pp &rarr; s &chi;&chi;) × B(s&rarr; VV) process for m_Z'=1.7 TeV, shown in dashed blue.

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, &Delta; &phi;_{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&gt;60 GeV that is supplemented by up to two additional small-R jets within &Delta;R_cone=2.5 of the TAR jet. Here, &Delta; &phi;_{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 &times; &varepsilon;), 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^&plusmn;W^&#8723;) dark Higgs signal points with dark Higgs boson mass m_s and Z' boson mass m_Z'.

The product of acceptance and efficiency (A &times; &varepsilon;), 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'.