$M_{\tau}$ distribution in signal region, (OS$\mu$, Belle II)

$M_{\tau}$ distribution in signal region, (SS$e$, Belle)

$M_{\tau}$ distribution in signal region, (SS$e$, Belle II)

$M_{\tau}$ distribution in signal region, (SS$\mu$, Belle)

$M_{\tau}$ distribution in signal region, (SS$\mu$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$ and $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (OS$e$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$ and $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (OS$e$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$ and $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (OS$\mu$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$ and $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (OS$\mu$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$ and $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (SS$e$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$ and $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (SS$e$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$ and $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (SS$\mu$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$ and $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (SS$\mu$, Belle II)

Efficiency as function the helicity angle of the kaon $\cos\theta_{K}$, and the helicity angle of the lepton $\cos\theta_{\ell}$ in bin e the four-momentum-transfer $q^{2}$, (OS$e$, Belle)

Efficiency as function the helicity angle of the kaon $\cos\theta_{K}$, and the helicity angle of the lepton $\cos\theta_{\ell}$ in bin e the four-momentum-transfer $q^{2}$, (OS$e$, Belle II)

Efficiency as function the helicity angle of the kaon $\cos\theta_{K}$, and the helicity angle of the lepton $\cos\theta_{\ell}$ in bin e the four-momentum-transfer $q^{2}$, (OS$\mu$, Belle)

Efficiency as function the helicity angle of the kaon $\cos\theta_{K}$, and the helicity angle of the lepton $\cos\theta_{\ell}$ in bin e the four-momentum-transfer $q^{2}$, (OS$\mu$, Belle II)

Efficiency as function the helicity angle of the kaon $\cos\theta_{K}$, and the helicity angle of the lepton $\cos\theta_{\ell}$ in bin e the four-momentum-transfer $q^{2}$, (SS$e$, Belle)

Efficiency as function the helicity angle of the kaon $\cos\theta_{K}$, and the helicity angle of the lepton $\cos\theta_{\ell}$ in bin e the four-momentum-transfer $q^{2}$, (SS$e$, Belle II)

Efficiency as function the helicity angle of the kaon $\cos\theta_{K}$, and the helicity angle of the lepton $\cos\theta_{\ell}$ in bin e the four-momentum-transfer $q^{2}$, (SS$\mu$, Belle)

Efficiency as function the helicity angle of the kaon $\cos\theta_{K}$, and the helicity angle of the lepton $\cos\theta_{\ell}$ in bin e the four-momentum-transfer $q^{2}$, (SS$\mu$, Belle II)

Efficiency as function of the helicity angle of the kaon $\cos\theta_{K}$, (OS$e$, Belle)

Efficiency as function of the helicity angle of the kaon $\cos\theta_{K}$, (OS$e$, Belle II)

Efficiency as function of the helicity angle of the kaon $\cos\theta_{K}$, (OS$\mu$, Belle)

Efficiency as function of the helicity angle of the kaon $\cos\theta_{K}$, (OS$\mu$, Belle II)

Efficiency as function of the helicity angle of the kaon $\cos\theta_{K}$, (SS$e$, Belle)

Efficiency as function of the helicity angle of the kaon $\cos\theta_{K}$, (SS$e$, Belle II)

Efficiency as function of the helicity angle of the kaon $\cos\theta_{K}$, (SS$\mu$, Belle)

Efficiency as function of the helicity angle of the kaon $\cos\theta_{K}$, (SS$\mu$, Belle II)

Efficiency as function of the helicity angle of the lepton $\cos\theta_{\ell}$, (OS$e$, Belle)

Efficiency as function of the helicity angle of the lepton $\cos\theta_{\ell}$, (OS$e$, Belle II)

Efficiency as function of the helicity angle of the lepton $\cos\theta_{\ell}$, (OS$\mu$, Belle)

Efficiency as function of the helicity angle of the lepton $\cos\theta_{\ell}$, (OS$\mu$, Belle II)

Efficiency as function of the helicity angle of the lepton $\cos\theta_{\ell}$, (SS$e$, Belle)

Efficiency as function of the helicity angle of the lepton $\cos\theta_{\ell}$, (SS$e$, Belle II)

Efficiency as function of the helicity angle of the lepton $\cos\theta_{\ell}$, (SS$\mu$, Belle)

Efficiency as function of the helicity angle of the lepton $\cos\theta_{\ell}$, (SS$\mu$, Belle II)

Efficiency as function of the acoplanarity angle $\phi$, (OS$e$, Belle)

Efficiency as function of the acoplanarity angle $\phi$, (OS$e$, Belle II)

Efficiency as function of the acoplanarity angle $\phi$, (OS$\mu$, Belle)

Efficiency as function of the acoplanarity angle $\phi$, (OS$\mu$, Belle II)

Efficiency as function of the acoplanarity angle $\phi$, (SS$e$, Belle)

Efficiency as function of the acoplanarity angle $\phi$, (SS$e$, Belle II)

Efficiency as function of the acoplanarity angle $\phi$, (SS$\mu$, Belle)

Efficiency as function of the acoplanarity angle $\phi$, (SS$\mu$, Belle II)

Efficiency as function of the four-momentum-transfer $q^{2}$, (OS$e$, Belle)

Efficiency as function of the four-momentum-transfer $q^{2}$, (OS$e$, Belle II)

Efficiency as function of the four-momentum-transfer $q^{2}$, (OS$\mu$, Belle)

Efficiency as function of the four-momentum-transfer $q^{2}$, (OS$\mu$, Belle II)

Efficiency as function of the four-momentum-transfer $q^{2}$, (SS$e$, Belle)

Efficiency as function of the four-momentum-transfer $q^{2}$, (SS$e$, Belle II)

Efficiency as function of the four-momentum-transfer $q^{2}$, (SS$\mu$, Belle)

Efficiency as function of the four-momentum-transfer $q^{2}$, (SS$\mu$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$, (OS$e$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$, (OS$e$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$, (OS$\mu$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$, (OS$\mu$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$, (SS$e$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$, (SS$e$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$, (SS$\mu$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}\ell)$, (SS$\mu$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (OS$e$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (OS$e$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (OS$\mu$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (OS$\mu$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (SS$e$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (SS$e$, Belle II)

Efficiency as function of the invariant mass squared $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (SS$\mu$, Belle)

Efficiency as function of the invariant mass squared $m^2(K^{*0}t_{\tau})$, where $t_{\tau}$ is the track from the $\tau$, (SS$\mu$, Belle II)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Example description

Differential fiducial cross sections for the number of jets

Example description

Differential fiducial cross sections for the leading jet pT

Example description

Observed best fit differential cross section for the $|y_{H}|$ observable

Observed best fit differential cross section for the $m_{jj}$ (GeV) observable

Observed best fit differential cross section for the $\Delta\eta_{jj}$ observable.

Observed best fit differential cross section for the $\tau_{C}^{j}$ (GeV) observable

Bin-to-bin correlation matrix for the $p_{T}^{H}$ spectrum

Bin-to-bin correlation matrix for the $N_{jets}$ spectrum

Bin-to-bin correlation matrix for the $p_{T}^{j0}$ spectrum

Bin-to-bin correlation matrix for the $|y_{H}|$ spectrum

Bin-to-bin correlation matrix for the $m_{jj}$ spectrum

Bin-to-bin correlation matrix for the $\Delta\eta_{jj}$ spectrum

Bin-to-bin correlation matrix for the $\tau_{C}^{j}$ spectrum

Rotation matrix for the PCA of the SMEFT Wilson coefficients

Correlation matrix of the linear combinations of Wilson coefficients obtained from the PCA, obtained by fitting the observed data to the $p_{T}^{H}$ spectra of all decay channels.

Summary of observed and expected confidence intervals for the first eigenvectors obtained from the PCA of the SMEFT Wilson coefficients.

Signal-selection efficiency as a function of the di-tau invariant mass squared $q^2$ for simulated events in the SR. The error bars indicate the statistical uncertainty.

The unfolded dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with leading and subleading jets at forward and midrapidity regions, respectively.

The unfolded dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with leading and subleading jets at backward and midrapidity regions, respectively.

The ratio of unfolded dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with both jets at the midrapidity regions.

The ratio of unfolded dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at midrapidity and forward regions, respectively.

The ratio of unfolded dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at midrapidity and backward regions, respectively.

The ratio of unfolded dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at forward and midrapidity regions, respectively.

The ratio of unfolded dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at backward and midrapidity regions, respectively.

Mean xj ratio, $\langle x_{j} \rangle$, in the range $10-60$ as unfolded for different multiplicities, for different leading and subleading jet combinations for the midrapidity, forward, and backward regions

The reconstructed (raw) dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with both jets at the midrapidity regions.

The reconstructed (raw) dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with leading and subleading jets at midrapidity and forward regions, respectively.

The reconstructed (raw) dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with leading and subleading jets at midrapidity and backward regions, respectively.

The reconstructed (raw) dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with leading and subleading jets at forward and midrapidity regions, respectively.

The reconstructed (raw) dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with leading and subleading jets at backward and midrapidity regions, respectively.

The ratio of reconstructed (raw) dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with both jets at the midrapidity regions.

The ratio of reconstructed (raw) dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at midrapidity and forward regions, respectively.

The ratio of reconstructed (raw) dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at midrapidity and backward regions, respectively.

The ratio of reconstructed (raw) dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at forward and midrapidity regions, respectively.

The ratio of reconstructed (raw) dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at backward and midrapidity regions, respectively.

Mean xj ratio, $\langle x_{j} \rangle$, in the range $10-60$ as reconstructed (raw) for different multiplicities, for different leading and subleading jet combinations for the midrapidity, forward, and backward regions

An example of the measured energy distribution in a single OHCal tower, showing the MIP distribution from a GEANT-4 simulation of cosmic-ray muons from EcoMug.

$\frac{dE_T}{d\eta}$ measurements in Au+Au collisions at $\sqrt{s_\mathrm{_{NN}}} = 200$ GeV over the pseudorapidity range $\left|\eta\right| < 1.1$.

Measurement of $\frac{dE_T}{d\eta}$ in 0-5% central Au+Au events as a function of $\eta$ over the range $\left|\eta\right| < 1.1$.

Measured $\frac{dE_T}{d\eta}$ normalized by the estimated number of participant pairs ($0.5$$N_{part}$), as a function of $N_{part}$, using the full calorimeter system.

Measured $\frac{dE_T}{d\eta}$ normalized by the estimated number of participant pairs ($0.5$$N_{part}$), as a function of $N_{part}$, predictions from AMPT.

Measured $\frac{dE_T}{d\eta}$ normalized by the estimated number of participant pairs ($0.5$$N_{part}$), as a function of $N_{part}$, predictions from EPOS.

Measured $\frac{dE_T}{d\eta}$ normalized by the estimated number of participant pairs ($0.5$$N_{part}$), as a function of $N_{part}$, predictions from HIJING.