Figure 6. Postfit $m_{j\gamma}$ spectra in the SRZ2. The lower panel shows the pull distributions with respect to the best-fit function. The signals with the largest local significances are shown normalized to the observed cross section upper limits.

Figure 7. The 95% CL upper limits on the product of production cross section and branching fraction $\sigma\mathcal{B}$ for $Z'\rightarrow H\gamma$ 0.01% width. Observed (expected) limits are shown with solid (dashed) lines. The colored bands represent the 68% and 95% CL intervals for the expected limits. The red line represents the theory benchmark model used for $Z'$ signal simulation.

Figure 7. The 95% CL upper limits on the product of production cross section and branching fraction $\sigma\mathcal{B}$ for $S\rightarrow Z\gamma$ 0.014% width. Observed (expected) limits are shown with solid (dashed) lines. The colored bands represent the 68% and 95% CL intervals for the expected limits

Figure 7. The 95% CL upper limits on the product of production cross section and branching fraction $\sigma\mathcal{B}$ for $S\rightarrow Z\gamma$ 5.6% width. Observed (expected) limits are shown with solid (dashed) lines. The colored bands represent the 68% and 95% CL intervals for the expected limits

Figure 7. The 95% CL upper limits on the product of production cross section and branching fraction $\sigma\mathcal{B}$ for $S\rightarrow Z\gamma$ 10% width. Observed (expected) limits are shown with solid (dashed) lines. The colored bands represent the 68% and 95% CL intervals for the expected limits

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\omega}$ meson $c\tau$ for $m_{\tilde{\omega}}=5~\textrm{GeV}$ and $\mathcal{B}(\tilde{\omega}\to\mu\mu)=0.162$ in the vector portal model. It is assumed that $m_{\tilde{\omega}}=\tilde{\Lambda}=m_{\tilde{\eta}}$, where $m_{\tilde{\omega}}$, $\tilde{\Lambda}$, and $m_{\tilde{\eta}}$ are parameters of the dark sector: the mass of the spin-one meson, confinement scale, and the mass of the spin-zero meson, respectively.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\omega}$ meson $c\tau$ for $m_{\tilde{\omega}}=15~\textrm{GeV}$ and $\mathcal{B}(\tilde{\omega}\to\mu\mu)=0.15$ in the vector portal model. It is assumed that $m_{\tilde{\omega}}=\tilde{\Lambda}=m_{\tilde{\eta}}$, where $m_{\tilde{\omega}}$, $\tilde{\Lambda}$, and $m_{\tilde{\eta}}$ are parameters of the dark sector: the mass of the spin-one meson, confinement scale, and the mass of the spin-zero meson, respectively.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\omega}$ meson $c\tau$ for $m_{\tilde{\omega}}=20~\textrm{GeV}$ and $\mathcal{B}(\tilde{\omega}\to\mu\mu)=0.15$ in the vector portal model. It is assumed that $m_{\tilde{\omega}}=\tilde{\Lambda}=m_{\tilde{\eta}}$, where $m_{\tilde{\omega}}$, $\tilde{\Lambda}$, and $m_{\tilde{\eta}}$ are parameters of the dark sector: the mass of the spin-one meson, confinement scale, and the mass of the spin-zero meson, respectively.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\textrm{A}'$ $c\tau$ for $m_{\tilde{\pi}_{3}}=1~\textrm{GeV}$ and $m_{\textrm{A}'}=0.33~\textrm{GeV}$ , and $\mathcal{B}(\textrm{A}'\to\mu\mu)=0.458$ in Scenario A. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\textrm{A}'$ $c\tau$ for $m_{\tilde{\pi}_{3}}=4~\textrm{GeV}$ and $m_{\textrm{A}'}=0.4~\textrm{GeV}$ , and $\mathcal{B}(\textrm{A}'\to\mu\mu)=0.436$ in Scenario A. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\textrm{A}'$ $c\tau$ for $m_{\tilde{\pi}_{3}}=2~\textrm{GeV}$ and $m_{\textrm{A}'}=0.67~\textrm{GeV}$, and $\mathcal{B}(\textrm{A}'\to\mu\mu)=0.17$ in Scenario A. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\textrm{A}'$ $c\tau$ for $m_{\tilde{\pi}_{3}}=4~\textrm{GeV}$ and $m_{\textrm{A}'}=1.33~\textrm{GeV}$, and $\mathcal{B}(\textrm{A}'\to\mu\mu)=0.305$ in Scenario A. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\pi}_{3}$ meson $c\tau$ for $m_{\tilde{\pi}_{3}}=1~\textrm{GeV}$ and $m_{\textrm{A}'}=0.33~\textrm{GeV}$ , and $\mathcal{B}(\textrm{A}'\to\mu\mu)=0.458$ in Scenario B1. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\pi}_{3}$ meson $c\tau$ for $m_{\tilde{\pi}_{3}}=2~\textrm{GeV}$ and $m_{\textrm{A}'}=0.4~\textrm{GeV}$ , and $\mathcal{B}(\textrm{A}'\to\mu\mu)=0.436$ in Scenario B1. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\pi}_{3}$ meson $c\tau$ for $m_{\tilde{\pi}_{3}}=2~\textrm{GeV}$ and $m_{\textrm{A}'}=0.67~\textrm{GeV}$ , and $\mathcal{B}(\textrm{A}'\to\mu\mu)=0.17$ in Scenario B1. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\pi}_{3}$ meson $c\tau$ for $m_{\tilde{\pi}_{3}}=4~\textrm{GeV}$ and $m_{\textrm{A}'}=1.33~\textrm{GeV}$, and $\mathcal{B}(\textrm{A}'\to\mu\mu)=0.305$ in Scenario B1. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\omega}$ meson mass for $c\tau=10~\textrm{mm}$ in the vector portal model. It is assumed that $m_{\tilde{\omega}}=\tilde{\Lambda}=m_{\tilde{\eta}}$, where $m_{\tilde{\omega}}$, $\tilde{\Lambda}$, and $m_{\tilde{\eta}}$ are parameters of the dark sector: the mass of the spin-one meson, confinement scale, and the mass of the spin-zero meson, respectively. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\omega}$ meson mass for $c\tau=100~\textrm{mm}$ in the vector portal model. It is assumed that $m_{\tilde{\omega}}=\tilde{\Lambda}=m_{\tilde{\eta}}$, where $m_{\tilde{\omega}}$, $\tilde{\Lambda}$, and $m_{\tilde{\eta}}$ are parameters of the dark sector: the mass of the spin-one meson, confinement scale, and the mass of the spin-zero meson, respectively. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\to\psi\overline{\psi}\right)$ as a function of the $\mathrm{A}'$ mass in Scenario A, for $\mathrm{A}'$ $c\tau=10~\textrm{mm}$, and $m_{\tilde{\pi}_{3}}=10m_{\textrm{A}'}$. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and ${\sin\theta=0.1}$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\to \mathrm{A}'\mathrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\to\psi\overline{\psi}\right)$ as a function of the $\mathrm{A}'$ mass in Scenario A, for $\mathrm{A}'$ $c\tau=10~\textrm{mm}$, and $m_{\tilde{\pi}_{3}}=3m_{\textrm{A}'}$. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and ${\sin\theta=0.1}$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\to \mathrm{A}'\mathrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\to\psi\overline{\psi}\right)$ as a function of the $\mathrm{A}'$ mass in Scenario A, for $\mathrm{A}'$ $c\tau=100~\textrm{mm}$, and $m_{\tilde{\pi}_{3}}=10m_{\textrm{A}'}$. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and ${\sin\theta=0.1}$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\to \mathrm{A}'\mathrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\to\psi\overline{\psi}\right)$ as a function of the $\mathrm{A}'$ mass in Scenario A, for $\mathrm{A}'$ $c\tau=100~\textrm{mm}$, and $m_{\tilde{\pi}_{3}}=3m_{\textrm{A}'}$. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and ${\sin\theta=0.1}$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\to \mathrm{A}'\mathrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\to\psi\overline{\psi}\right)$ as a function of the $\mathrm{A}'$ mass in Scenario B1, for $\tilde{\pi}_{3}$ $c\tau=10~\textrm{mm}$, and $m_{\tilde{\pi}_{3}}=10m_{\textrm{A}'}$. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and ${\sin\theta=0.1}$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\to \mathrm{A}'\mathrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\to\psi\overline{\psi}\right)$ as a function of the $\mathrm{A}'$ mass in Scenario B1, for $\tilde{\pi}_{3}$ $c\tau=10~\textrm{mm}$, and $m_{\tilde{\pi}_{3}}=3m_{\textrm{A}'}$. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and ${\sin\theta=0.1}$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\to \mathrm{A}'\mathrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\to\psi\overline{\psi}\right)$ as a function of the $\mathrm{A}'$ mass in Scenario B1, for $\tilde{\pi}_{3}$ $c\tau=100~\textrm{mm}$, and $m_{\tilde{\pi}_{3}}=10m_{\textrm{A}'}$. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and ${\sin\theta=0.1}$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\to \mathrm{A}'\mathrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\to\psi\overline{\psi}\right)$ as a function of the $\mathrm{A}'$ mass in Scenario B1, for $\tilde{\pi}_{3}$ $c\tau=100~\textrm{mm}$, and $m_{\tilde{\pi}_{3}}=3m_{\textrm{A}'}$. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and ${\sin\theta=0.1}$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\to \mathrm{A}'\mathrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Observed upper limits at 95% CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as a function of the $\tilde{\omega}$ meson mass and $c\tau$ in the vector portal model. The phase space region that is left out is kinematically forbidden in the model. It is assumed that $m_{\tilde{\omega}}=\tilde{\Lambda}=m_{\tilde{\eta}}$, where $m_{\tilde{\omega}}$, $\tilde{\Lambda}$, and $m_{\tilde{\eta}}$ are parameters of the dark sector: the mass of the spin-one meson, confinement scale, and the mass of the spin-zero meson, respectively. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Observed upper limits at $95\%$ CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as functions of the $\textrm{A}'$ mass and $c\tau$ for $m_{\tilde{\pi}_{3}}=10m_{A'}$ in Scenario A. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Observed upper limits at $95\%$ CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as functions of the $\textrm{A}'$ mass and $c\tau$ for $m_{\tilde{\pi}_{3}}=3m_{A'}$ in Scenario A. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Observed upper limits at $95\%$ CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as functions of the $\textrm{A}'$ mass and $c\tau$ for $m_{\tilde{\pi}_{3}}=10m_{A'}$ in Scenario B1. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Observed upper limits at $95\%$ CL on the branching fraction $\mathcal{B}\left(\textrm{H}\rightarrow\psi\overline{\psi}\right)$ as functions of the $\textrm{A}'$ mass and $c\tau$ for $m_{\tilde{\pi}_{3}}=3m_{A'}$ in Scenario B1. It is assumed that $m_{\tilde{\eta}}=\tilde{\Lambda}=4m_{\tilde{\pi}_{2}}$ and $\sin\theta=0.1$, where $m_{\tilde{\eta}}$ is the mass of the dark-sector pseudoscalar meson and $\theta$ is the mixing angle parametrising the isospin violation. The branching fraction $\mathcal{B}(\tilde{\pi}_{3}\rightarrow \textrm{A}'\textrm{A}')$ is assumed to be 1. The grey bands correspond to mass regions of known SM resonances, which are masked in the search.

Cutflow efficiencies for representative signal model points in the vector portal.

Cutflow efficiencies for representative signal model points in Scenario A.

Cutflow efficiencies for representative signal model points in Scenario B1.

Observed 95% CL upper limits on the $\tilde{t}$ production cross section, as functions of $m_{\tilde{t}}$ and $\Delta m$, for $\mathcal{B}(\tilde{t} \to bf\overline{f}'\tilde{\chi}^{0}_{1})$ of 50%. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the $\tilde{t}$ production cross section, as functions of $m_{\tilde{t}}$ and $\Delta m$, for $\mathcal{B}(\tilde{t} \to bf\overline{f}'\tilde{\chi}^{0}_{1})$ of 50%. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the $\tilde{t}$ production cross section, as functions of $m_{\tilde{t}}$ and $\Delta m$, for $\mathcal{B}(\tilde{t} \to bf\overline{f}'\tilde{\chi}^{0}_{1})$ of 100%. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the $\tilde{t}$ production cross section, as functions of $m_{\tilde{t}}$ and $\Delta m$, for $\mathcal{B}(\tilde{t} \to bf\overline{f}'\tilde{\chi}^{0}_{1})$ of 100%. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the production cross section for the bino-wino NLSP model, as functions of $m_{ ext{LLP}}$ and $c\tau$, for $\Delta m$ of 12 GeV. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the production cross section for the bino-wino NLSP model, as functions of $m_{ ext{LLP}}$ and $c\tau$, for $\Delta m$ of 12 GeV. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the production cross section for the bino-wino NLSP model, as functions of the $m_{ ext{LLP}}$ and $c\tau$, for $\Delta m$ of 15 GeV. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the production cross section for the bino-wino NLSP model, as functions of the $m_{ ext{LLP}}$ and $c\tau$, for $\Delta m$ of 15 GeV. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the production cross section for the bino-wino NLSP model, as functions of the $m_{ ext{LLP}}$ and $c\tau$, for $\Delta m$ of 20 GeV. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the production cross section for the bino-wino NLSP model, as functions of the $m_{ ext{LLP}}$ and $c\tau$, for $\Delta m$ of 20 GeV. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the production cross section for the bino-wino NLSP model, as functions of the $m_{ ext{LLP}}$ and $c\tau$, for $\Delta m$ of 25 GeV. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Observed 95% CL upper limits on the production cross section for the bino-wino NLSP model, as functions of the $m_{ ext{LLP}}$ and $c\tau$, for $\Delta m$ of 25 GeV. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the plots. The search excludes the region to the left of the exclusion curves.

Example description

Example description

Normalized differential distributions of the $\Delta R (\gamma,\mathrm{tt})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Absolute differential distributions of the min. $\Delta R (\gamma,\mathrm{t})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Normalized differential distributions of the min. $\Delta R (\gamma,\mathrm{t})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Absolute differential distributions of the $m(\mathrm{tt})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Normalized differential distributions of the $m(\mathrm{tt})$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Absolute differential distributions of the leading lepton $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Normalized differential distributions of the leading lepton $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Absolute differential distributions of the photon $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Normalized differential distributions of the photon $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Absolute differential distributions of the $\Delta \phi (\ell,\ell)$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Normalized differential distributions of the $\Delta \phi (\ell,\ell)$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Absolute differential distributions of the ratio between the cross sections of $\mathrm{tt\gamma}$ and $\mathrm{tt}$, as a function of the leading top quark $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process.

Absolute differential distributions of the ratio between the cross sections of $\mathrm{tt\gamma}$ and $\mathrm{tt}$, as a function of the leading lepton $p_{\mathrm{T}}$. The nominal MC prediction used to compare the experimental results to is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and Madgraph5 at LO in QCD for photons from the decay part of the process. The alternative prediction is obtained with Madgraph5 at NLO in QCD for photons from the production part of the process and a POWHEG+Pythia $\mathrm{tt}$ simulation at NLO in QCD for photons from the decay part of the process. The theory predictions are computed at fixed order, at NLO in QCD, as described in Section 9.2 of the paper.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the leading top quark $p_{\mathrm{T}}$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the leading top quark $p_{\mathrm{T}}$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the $\Delta R (\gamma,\mathrm{tt})$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the $\Delta R (\gamma,\mathrm{tt})$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the min. $\Delta R (\gamma,\mathrm{t})$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the min. $\Delta R (\gamma,\mathrm{t})$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the $m(\mathrm{tt})$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the $m(\mathrm{tt})$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the leading lepton $p_{\mathrm{T}}$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the leading lepton $p_{\mathrm{T}}$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the photon $p_{\mathrm{T}}$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the photon $p_{\mathrm{T}}$.

Correlation matrix of the total uncertainty in the absolute differential cross section as a function of the $\Delta \phi (\ell,\ell)$.

Correlation matrix of the statistical uncertainty in the absolute differential cross section as a function of the $\Delta \phi (\ell,\ell)$.

Correlation matrix of the total uncertainty in the absolute differential $\mathrm{tt\gamma}$ / $\mathrm{tt}$ cross section ratio as a function of the leading top $p_{\mathrm{T}}$.

Correlation matrix of the statistical uncertainty in the absolute differential $\mathrm{tt\gamma}$ / $\mathrm{tt}$ cross section ratio as a function of the leading top $p_{\mathrm{T}}$.

Correlation matrix of the total uncertainty in the absolute differential $\mathrm{tt\gamma}$ / $\mathrm{tt}$ cross section ratio as a function of the leading lepton $p_{\mathrm{T}}$.

Correlation matrix of the statistical uncertainty in the absolute differential $\mathrm{tt\gamma}$ / $\mathrm{tt}$ cross section ratio as a function of the leading lepton $p_{\mathrm{T}}$.

Distributions in the trailing lepton pt for the DY CR after a fit to the data.

Distributions in the trailing lepton pt for the tt CR after a fit to the data.

Likelihood scan of the signal strengths for tWZ and ttZ production. The black cross shows the best fit value, while the black diamond indicates the SM expected value. The con- tours in different shades of blue correspond to the 68, 95, and 99% confidence level.

The $v_{3}\{2,\lvert\Delta\eta\rvert>2\}$ ratios for charged particles as functions of centrality in NeNe to OO collisions at 5.36 TeV.

Signal strength and significance for $\text{t}\overline{\text{t}}\text{Z}(\text{Z}\to\text{b}\overline{\text{b}})$ decays with respect to the standard model expectation of unity.

Signal strength and significance for $\text{t}\overline{\text{t}}\text{Z}(\text{Z}\to\text{c}\overline{\text{c}})$ decays with respect to the standard model expectation of unity.

2D likelihood scans for $\kappa_{b}$ and $\kappa_{c}$ for $\text{t}\overline{\text{t}}\text{H}(\text{H}\to\text{c}\overline{\text{c}})$ and $\text{t}\overline{\text{t}}\text{H}(\text{H}\to\text{b}\overline{\text{b}})$.

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

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

The $m_{jj}$ and $m_{J}$ projections for the number of observed events (black markers) compared with the backgrounds estimated in the fit to the data (filled histograms) in the CR. Pass and Fail categories are shown. The high level of agreement between the model and the data in the Fail region is due to the nature of the background estimate. The lower panels show the ``Pull'' defined as $(\text{observed events}{-}\text{expected events})/\sqrt{\smash[b]{\sigma_\text{obs}^{2} + \sigma_\text{exp}^{2}}}$, where $\sigma_\text{obs}$ and $\sigma_\text{exp}$ are the total uncertainties in the observation and the background estimation, respectively.

The $m_{jj}$ and $m_{J}$ projections for the number of observed events (black markers) compared with the backgrounds estimated in the fit to the data (filled histograms) in the MR. Pass and Fail categories are shown. The expected contribution from the signal benchmark with $M_{X} = 2200\,\mathrm{GeV}$ and $M_{Y} = 250\,\mathrm{GeV}$ is overlaid in the MR Pass, assuming a production cross section of 5 fb. The high level of agreement between the model and the data in the Fail region is due to the nature of the background estimate.

The $m_{jj}$ and $m_{J}$ projections for the number of observed events (black markers) compared with the backgrounds estimated in the fit to the data (filled histograms) in the MR. Pass and Fail categories are shown. The expected contribution from the signal benchmark with $M_{X} = 2200\,\mathrm{GeV}$ and $M_{Y} = 250\,\mathrm{GeV}$ is overlaid in the MR Pass, assuming a production cross section of 5 fb. The high level of agreement between the model and the data in the Fail region is due to the nature of the background estimate.

The $m_{jj}$ and $m_{J}$ projections for the number of observed events (black markers) compared with the backgrounds estimated in the fit to the data (filled histograms) in the MR. Pass and Fail categories are shown.

The $m_{jj}$ and $m_{J}$ projections for the number of observed events (black markers) compared with the backgrounds estimated in the fit to the data (filled histograms) in the MR. Pass and Fail categories are shown.

The expected and observed 95% confidence level upper limits on the cross section $\sigma(\mathrm{pp} \to \mathrm{X} \to (\mathrm{H} \to b\bar{b})(\mathrm{Y} \to WW \to 4q))$ are shown for different values of $M_X$ and $M_Y$. The limits have been evaluated in discrete steps corresponding to the centers of the boxes. The numbers in the boxes are given in fb. In addition to the observed and median expected limits, the expected limits at $\pm1\sigma$ and $\pm2\sigma$ quantiles are also provided. These represent the regions containing 68% and 95%, respectively, of the distribution of limits expected under the background-only hypothesis.

The median expected (dashed line) and observed (solid line) 95% confidence level upper limits on the main and three alternative signal scenarios as a function of $M_{X}$. The inner (green) band and outer (yellow) band represent the regions containing 68% and 95%, respectively, of the distribution of limits expected under the background-only hypothesis.

The median expected (dashed line) and observed (solid line) 95% confidence level upper limits on the main and three alternative signal scenarios as a function of $M_{X}$. The inner (green) band and outer (yellow) band represent the regions containing 68% and 95%, respectively, of the distribution of limits expected under the background-only hypothesis.

The median expected (dashed line) and observed (solid line) 95% confidence level upper limits on the main and three alternative signal scenarios as a function of $M_{X}$. The inner (green) band and outer (yellow) band represent the regions containing 68% and 95%, respectively, of the distribution of limits expected under the background-only hypothesis.

The median expected (dashed line) and observed (solid line) 95% confidence level upper limits on the main and three alternative signal scenarios as a function of $M_{X}$. The inner (green) band and outer (yellow) band represent the regions containing 68% and 95%, respectively, of the distribution of limits expected under the background-only hypothesis.

Selection efficiency as a function of MX and MY

Selection efficiency as a function of MX for fixed MY

Selection efficiency as a function of MX for fixed MY

Selection efficiency as a function of MX for fixed MY