Di-muon mass distribution for the $\mu\mu\mu$ signal region for data and expected background. The expected signal distribution for $m_a = 35$ GeV is shown assuming $\sigma(t\bar{t}a)\times \text{Br}(a\rightarrow\mu\mu) = 4$ fb. Rare background processes include $ZZ+$jets, $WWZ$, $WZZ$, $ZZZ$, and $t\bar{t}t\bar{t}$.

Fit of the Double Crystal Ball width as a function of the $a$-boson mass in the $e\mu\mu$ channel.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $t\bar{t}a$.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $H^{+}\rightarrow aW$, considering a top pair cross section of $\sigma(pp\rightarrow t\bar{t})= 833$ pb for the charged Higgs boson mass, $m_{H^{+}} = 120$ GeV.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $H^{+}\rightarrow aW$, considering a top pair cross section of $\sigma(pp\rightarrow t\bar{t})= 833$ pb for the charged Higgs boson mass, $m_{H^{+}} = 140$ GeV.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $H^{+}\rightarrow aW$, considering a top pair cross section of $\sigma(pp\rightarrow t\bar{t})= 833$ pb for the charged Higgs boson mass, $m_{H^{+}} = 160$ GeV.

Comparison of the expected background with data in a region dominated by non-prompt background. The region is defined exactly like the $e\mu\mu$ signal region but requiring the two selected muons to have the same electrical charge.

Scan of the observed p-value as a function of $m_a$ for the background-only hypothesis.

Scan of the observed p-value as a function of $m_a$ for the background-only hypothesis for $e\mu\mu$ channel.

Scan of the observed p-value as a function of $m_a$ for the background-only hypothesis for $\mu\mu\mu$ channel.

Acceptance times efficiency for the $t\bar{t}a$, $a\rightarrow \mu\mu$ signal as a function of the $m_a$ hypothesis.

Acceptance times efficiency for the $t \rightarrow bH^+$, $H^+ \rightarrow W^+ a$, $a\rightarrow \mu\mu$ signal as a function of the $m_a$ and $m_{H^+}$ hypotheses, for the $e\mu\mu$ signal region. The upper left region of the plot corresponds to signals for which on-shell decays of $H^+ \rightarrow Wa$ are not possible. The CP conjugate process is also included.

Acceptance times efficiency for the $t \rightarrow bH^+ $, $H^+ \rightarrow W^+ a$, $a\rightarrow \mu\mu$ signal as a function of the $m_a$ and $m_{H^+}$ hypotheses, for the $\mu\mu\mu$ signal region. The upper left region of the plot corresponds to signals for which on-shell decays of $H^+ \rightarrow Wa$ are not possible. The CP conjugate process is also included.

Expected lower limits on $\tan\beta$ for the $t\rightarrow H^{+}b, H^{+}\rightarrow aW, a\rightarrow\mu\mu$ signal in the context of the 2HDM as a function of $m_{H^{+}}$ and $m_a$.

Observed lower limits on $\tan\beta$ for the $t\rightarrow H^{+}b, H^{+}\rightarrow aW, a\rightarrow\mu\mu$ signal in the context of the 2HDM as a function of $m_{H^{+}}$ and $m_a$.

Cutflow for two signal mass points for $t\bar{t}a$, $m_{a} = 20$ GeV and $m_{a} = 60$ GeV, and one signal mass point for $H^+ \rightarrow Wa$ for the electron channel $e\mu\mu$. For the signal yields, a cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Cutflow for the dominant backgrounds estimated from simulation ($t\bar{t}Z$, $WZ$, $t\bar{t}H$) and data, for the electron channel $e\mu\mu$. The yields are presented before the profile likelihood fit.

Cutflow for the signal mass point $t\bar{t}a$, $m_{a} = 20$ GeV, as well as the dominant backgrounds estimated from simulation ($t\bar{t}Z$, $WZ$, $t\bar{t}H$) and data, for the corresponding signal mass hypothesis, for the muon channel $\mu\mu\mu$. For the signal yields, a cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Cutflow for the signal mass point $t\bar{t}a$, $m_{a} = 60$ GeV, as well as the dominant backgrounds estimated from simulation ($t\bar{t}Z$, $WZ$, $t\bar{t}H$) and data, for the corresponding signal mass hypothesis, for the muon channel $\mu\mu\mu$. For the signal yields, a cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Cutflow for $H^+ \rightarrow W a$ with $m_{H^{+}} = 120$ GeV, $m_{a} = 30$ GeV for the muon channel $\mu\mu\mu$. A cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Total yields in the $e\mu\mu$ and $\mu\mu\mu$ signal regions for the expected background and data after the profile likelihood fit for a $m_a = 30$ GeV signal mass hypothesis.

Average timing distributions for SR data and the estimated background as determined by the background-only fit, in each of the five exclusive $\rho$ categories. For comparison, the expected timing shapes for a few different signal models are superimposed, with each model labeled by the values of the $\tilde\chi^0_1$ mass and lifetime, as well as decay mode. To provide some indication of the variations in signal yield and shape, three signal models are shown for each of the $\tilde\chi^0_1$ decay modes, namely $\tilde\chi^0_1$ $\rightarrow$ $H \tilde G$ and $\tilde\chi^0_1$ $\rightarrow$ $Z \tilde G$. The models shown include a rather low $\tilde\chi^0_1$ mass value of 135 GeV for lifetimes of either 2 ns or 10 ns, and a higher $\tilde\chi^0_1$ mass value which is near the 95% CL exclusion limit for each decay mode for a lifetime of 2 ns. Each signal model is shown with the signal normalization corresponding to a BR value of unity for the decay mode in question.

Average timing distributions for SR data and the estimated background as determined by the background-only fit, in each of the five exclusive $\rho$ categories. For comparison, the expected timing shapes for a few different signal models are superimposed, with each model labeled by the values of the $\tilde\chi^0_1$ mass and lifetime, as well as decay mode. To provide some indication of the variations in signal yield and shape, three signal models are shown for each of the $\tilde\chi^0_1$ decay modes, namely $\tilde\chi^0_1$ $\rightarrow$ $H \tilde G$ and $\tilde\chi^0_1$ $\rightarrow$ $Z \tilde G$. The models shown include a rather low $\tilde\chi^0_1$ mass value of 135 GeV for lifetimes of either 2 ns or 10 ns, and a higher $\tilde\chi^0_1$ mass value which is near the 95% CL exclusion limit for each decay mode for a lifetime of 2 ns. Each signal model is shown with the signal normalization corresponding to a BR value of unity for the decay mode in question.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ mass (left) and $\tilde\chi^0_1$ lifetime (right), for the different decay modes of $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$) = 1 (top) and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 (bottom). For the limits as a function of mass (lifetime), several signal models with varying lifetime (mass) are overlaid for comparison. Included are the theoretical expectations from higgsino production for each mass hypothesis, calculated from a GMSB SUSY model that assumes nearly degenerate $\tilde\chi^0_1$, $\tilde\chi^\pm_1$, and $\tilde\chi^0_2$.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL limits on $\sigma(pp \rightarrow \tilde\chi^0_1 \tilde\chi^0_1$) in fb as a function of $\tilde\chi^0_1$ branching ratio to the SM Higgs boson, where the assumed cross-section is for higgsino production, and $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow Z +\tilde{G}$) = 1 - $\mathcal{B}$($\tilde\chi^0_1$ $\rightarrow H + \tilde{G}$). Several signal hypotheses are overlaid that are labelled by the $\tilde\chi^0_1$ mass, all with a fixed $\tilde\chi^0_1$ lifetime of 2 ns.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

The 95% CL exclusion limits on the target signal hypothesis, for $\tilde\chi^0_1$ lifetime in ns as a function of $\tilde\chi^0_1$ mass in GeV. The overlaid curves correspond to different decay hypotheses, where the assumed cross-section is for higgsino production, and the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ or $Z + \tilde{G}$ such that $\mathcal{B}(H + \tilde{G}) + \mathcal{B}(Z + \tilde{G})$ = 100%. The curve shown in red represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $Z + \tilde{G}$ with 100% branching ratio. The curve shown in blue represents the decay hypothesis where the $\tilde\chi^0_1$ decays to $H + \tilde{G}$ with 100% branching ratio.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 10 ns, in the H decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 10 ns, in the Z decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 2 ns, in the H decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Cutflow for an example higgsino signal with mass 225 GeV and lifetime 2 ns, in the Z decay mode. Acceptance is defined at truth level, and efficiency compares the events passing at reconstruction level with respect to truth.

Acceptance across the H decay mode signal grid, calculated using truth information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns).

Acceptance across the Z decay mode signal grid, calculated using truth information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns).

Efficiency across the H decay mode signal grid, calculated using reco information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns). Here, the numerator is the signal yield passing the reco selection and the denominator is the signal yield passing the truth selection.

Efficiency across the Z decay mode signal grid, calculated using reco information. The selection applied corresponds to the model-independent signal region (i.e. the standard SR with $t_{\text{avg}$ > 0.9 ns). Here, the numerator is the signal yield passing the reco selection and the denominator is the signal yield passing the truth selection.

Total efficiencies for the selection of signal events after categorization and application of the tighter photon $E_{\mathrm{T}}^{\gamma}$ selection used to optimize the signal significance spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$. In addition to the total efficiency, contributions to the signal selection from each of the separate event categories are shown. The efficiencies calculated from MC samples with $W/Z$ hadronic decays are shown as the points on each curve. The line presents interpolated results.

Total efficiencies for the selection of signal events after categorization and application of the tighter photon $E_{\mathrm{T}}^{\gamma}$ selection used to optimize the signal significance spin-2 $q\bar{q}\to X^0 \to Z\gamma$. In addition to the total efficiency, contributions to the signal selection from each of the separate event categories are shown. The efficiencies calculated from MC samples with $W/Z$ hadronic decays are shown as the points on each curve. The line presents interpolated results.

Total efficiencies for the selection of signal events after categorization and application of the tighter photon $E_{\mathrm{T}}^{\gamma}$ selection used to optimize the signal significance spin-2 $gg\to X^0 \to Z\gamma$. In addition to the total efficiency, contributions to the signal selection from each of the separate event categories are shown. The efficiencies calculated from MC samples with $W/Z$ hadronic decays are shown as the points on each curve. The line presents interpolated results.

The $m_{J\gamma}$ distributions of data events selected for the spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$ search in the WMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-0 $gg \to X^{0} \to Z\gamma$ search in the BTAG category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-0 $gg \to X^{0} \to Z\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-0 $gg \to X^{0} \to Z\gamma$ search in the ZMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $gg \to X^{0} \to Z\gamma$ search in the BTAG category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $gg \to X^{0} \to Z\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $gg \to X^{0} \to Z\gamma$ search in the ZMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $q\bar{q} \to X^{0} \to Z\gamma$ search in the BTAG category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $q\bar{q} \to X^{0} \to Z\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $q\bar{q} \to X^{0} \to Z\gamma$ search in the ZMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-0 $gg \to X^{0} \to Z\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-2 $gg \to X^{0} \to Z\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-2 $q\bar{q} \to X^{0} \to Z\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The jet mass distribution of large-$R$ jets originating from the hadronic decay of $W$ and $Z$ bosons produced from the decay of BSM bosons with mass $m_X = 2000$ GeV. The decays simulated are for the production models $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$ with a spin-1 resonance $X^{\pm}$ and $gg\to X^0 \to Z\gamma$ with a spin-0 resonance $X^{0}$.

The photon $E_{T}$ selection criteria as a function of resonance mass for D2 categories. Selections for four different signal channels are shown together.

The photon $E_{T}$ selection criteria as a function of resonance mass for ZMASS or WMASS categories. Selections for four different signal channels are shown together.

Background-only fit to the observed numbers of events (black points) as a function of the diphoton invariant mass $m_{\gamma\gamma}$.

Expected and observed 95% CL upper limits on the signal fiducial cross section $\sigma_{\textrm{fid}}$, assuming 100% branching ratio for ALP decay into two photons, as a function of the hypothetical ALP mass $m_{\textrm{X}}$. The 1$\sigma$ and 2$\sigma$ confidence intervals are shown by the coloured bands.

Expected and observed 95% CL upper limits on the the ALP coupling constant, assuming 100% branching ratio for ALP decay into two photons, as a function of the hypothetical ALP mass $m_{\textrm{X}}$. The 1$\sigma$ and 2$\sigma$ confidence intervals are shown by the coloured bands. Contours of the ALP natural width $\Gamma$ are illustrated by the smooth blue solid lines.

Fiducial region acceptance of signal events for coupling constant f$^{-1}=0.05$ TeV$^{-1}$. EL, SD, and DD stand for the exclusive, single-dissociative, and double-dissociative processes. The acceptances are parameterised as $p_0 + p^{i}_{1} e^{p_{2}^{i} m_{\textrm{X}}} + p_{3}^{i} e^{p_{4}^{i} m_{\textrm{X}}} $ , where $i$ specifies the process type, exclusive, single-dissociative, or double-dissociative.

Correction factors $A_{\textrm{X}}/A_{0}$ and $C_{\textrm{X}}$ for the exclusive signal events with ALP coupling constant f$^{-1}$=0.05 TeV$^{-1}$, where acceptance $A_{\textrm{X}}$ and efficiency $C_{\textrm{X}}$ are defined in arXiv:2102.13405. $A_{\textrm{X}}$/$A_{0}$ and $C_{\textrm{X}}$ are parameterised as $p_0 + p^{i}_{1} e^{p_{2}^{i} m_{\textrm{X}}} + p_{3}^{i} e^{p_{4}^{i} m_{\textrm{X}}}$, where $i$ specifies the process type.

Simulated ALP signal cutflow for $m_{\rm{X}}$=300 GeV and $m_{\rm{X}}$=1200 GeV. The yields are normalized to a luminosity of 14.6 fb$^{-1}$. Selections are applied for each side after the selection on acoplanarity. The union of the sets of selected events is taken finally.

Cutflow for two representative signal samples used in this analysis. The excited tau mass $m_{\tau^\ast} = 2.75$ TeV and the contact interaction scale $\Lambda=10$ TeV. The LQ mass $m_\mathrm{LQ} = 1.3$ TeV. The event yields include all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$.

Acceptance x efficiency of the $\tau^\ast$ signal SR selection

Acceptance x efficiency of the LQ signal SR selection

Signal selection efficiency ($\epsilon$) times acceptance ($A$) as a function of $H^{\pm}$. The estimated $\epsilon\times A$ arises from the lepton selection and triggering (∼30%) as well as jet selection and flavour tagging (∼10% or lower). The decrease of $\epsilon\times A$ for $m_{H^{\pm}}$ = 120 GeV and higher is expected from the kinematic constraint on the $H^{\pm}$ decay products due to the top-quark mass.

Fit results of the simultaneous measurements of the $H\to e\tau$ and $H\to \mu\tau$ signals (2POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to \mu\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to e\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to e\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to \mu\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to \mu\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $e\tau_{\mu}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $\mu\tau_{e}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $e\tau_{\mu}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $\mu\tau_{e}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $e\tau_{had}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $\mu\tau_{had}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $e\tau_{had}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $\mu\tau_{had}$ final state in the VBF category.

NN score distribution, after an independent fit of the $H\to \mu\tau$ signal, obtained by fitting the data of the Symmetry $\ell\tau_{\ell'}$ channel, for the $\mu\tau_{e}$ final state in the non-VBF category.

NN score distribution, after an independent fit of the $H\to \mu\tau$ signal, obtained by fitting the data of the Symmetry $\ell\tau_{\ell'}$ channel, for the $\mu\tau_{e}$ final state in the VBF category.

Best-fit $\mathcal{B}(H\to \mu\tau)- \mathcal{B}(H\to e\tau)$ values, given in %, obtained from the standalone fits of the $\ell\tau_\ell'$ final state with either background estimation method. For the MC-template method, the correlated uncertainties are fixed to their best-fit values and the uncertainties are re-derived considering only the uncorrelated sources. For the Symmetry method, the full uncertainty is considered.

Best-fit value of the branching ratios $\mathcal{B}(H\to e\tau)$ and $\mathcal{B}(H\to \mu\tau)$, given in %, and likelihood contours at 68% and 95% C.L. Obtained from the simultaneous fit of $H\to e\tau$ and $H\to \mu\tau$ signals based on the MC-template method, compared to the SM expectation.

The performance of jet energy resolution (JER) for jets with |y| < 2.1 evaluated as a function of pT_truth in different centrality bins. Simulated hard scatter events were overlaid onto events from a dedicated sample of minimum-bias Xe+Xe data. The fit parameters are listed in a sperate table (Extras 1)

The relative magnitude of systematic uncertainties for per-pair normalized xJ distributions in 0-10% Xe+Xe centrality

The relative magnitude of systematic uncertainties for absolutely normalized xJ distributions in 0-10% Xe+Xe centrality

The relative magnitude of systematic uncertainties for rho distributions for leading jets in 0-10% Xe+Xe centrality

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

Parameter a,b, and c from JER fits in Figure 1b.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Distribution of the NN output in the $\geq$1fj$\,$SR in data and the expected contribution of the signal and background processes after the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions considering the correlations of the uncertainties as obtained by the fit.

Distribution of the NN output in the $t\bar{t}\gamma\,$CR in data and the expected contribution of the signal and background processes after the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions considering the correlations of the uncertainties as obtained by the fit.

Total event yield in the $W\gamma\,$CR in data and the expected contribution of the signal and background processes after the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions considering the correlations of the uncertainties as obtained by the fit.

Distribution of the scalar sum of the jet transverse momenta in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions. The first and last bins include the underflow and overflow, respectively.

Distribution of the $\eta$ of the $b$-tagged jet in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions.

Distribution of the reconstructed top-quark mass in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions. The first and last bins include the underflow and overflow, respectively.

Distribution of the $p_{\mathrm{T}}$ of the top-quark+photon system in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the photon $p_{\mathrm{T}}$ in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the photon $\eta$ in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions.

Distribution of the scalar sum of the jet transverse momenta in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the invariant mass of the $b$-tagged jet and the highest-$p_{\mathrm{T}}$ forward jet in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of $p_{\mathrm{T}}$ of the highest-$p_{\mathrm{T}}$ forward jet in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the difference in $\eta$ between the highest-$p_{\mathrm{T}}$ forward jet and the photon in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the energy of the system formed by the highest-$p_{\mathrm{T}}$ forward jet and the photon in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions. The first and last bins include the underflow and overflow, respectively.

Distribution of the $\eta$ of the $b$-tagged jet in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions.

Distribution of the reconstructed top-quark mass in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions. The first and last bins include the underflow and overflow, respectively.

Distribution of the $p_{\mathrm{T}}$ of the top-quark and photon system in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the photon $p_{\mathrm{T}}$ in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the photon $\eta$ in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions.

Ordered list of the 30 systematic uncertainties with the largest impact on the measured signal normalisation in the fit to data in the parton-level measurement considered as nuisance parameters (NPs) in the profile-likelihood fit. The column "NP value, error" corresponds to the nominal best-fit values and the corresponding uncertainties. The impact of each NP, $\Delta\sigma$/$\sigma_{\mathrm{pred}}$, is computed by comparing the nominal best-fit value of the POI ($\sigma$/$\sigma_{\mathrm{pred}}$) with the result of the fits when fixing the considered NP to its best-fit value shifted by its pre-fit and post-fit uncertainties. The corresponding impacts are listed in the "POI impact prefit high/low" and "POI impact high/low" columns, respectively. The "MC stat." NPs represent the MC statistical uncertainty and they enter the likelihood with a Poisson term, while all the other NPs enter the likelihood via a Gaussian term.

Ordered list of the 30 systematic uncertainties with the largest impact on the measured signal normalisation in the fit to data in the particle-level measurement considered as nuisance parameters (NPs) in the profile-likelihood fit. The column "NP value, error" corresponds to the nominal best-fit values and the corresponding uncertainties. The impact of each NP, $\Delta\sigma$/$\sigma_{\mathrm{pred}}$, is computed by comparing the nominal best-fit value of the POI ($\sigma$/$\sigma_{\mathrm{pred}}$) with the result of the fits when fixing the considered NP to its best-fit value shifted by its pre-fit and post-fit uncertainties. The corresponding impacts are listed in the "POI impact prefit high/low" and "POI impact high/low" columns, respectively. The "MC stat." NPs represent the MC statistical uncertainty and they enter the likelihood with a Poisson term, while all the other NPs enter the likelihood via a Gaussian term.

This table lists the kinematic requirements on parton-level objects used to define of the fiducial phase space for the parton-level measurement. Frixione isolation ($\href{https://arxiv.org/abs/hep-ph/9801442}{\text{hep-ph/9801442}}$) with a chosen radius of $\Delta R = 0.2$ is applied to photons ($\gamma$). The measured fiducial parton-level cross section is $\sigma_{tq\gamma}\times\mathcal{B}(t\rightarrow l\nu b) = 688\pm 23(\text{stat.})^{+75}_{-71}(\text{syst.})\,$fb.

This table lists the kinematic requirements on particle-level objects used to define of the fiducial phase space for the particle-level measurement. The particle level objects are photons ($\gamma$) not from a hadron decay, neutrinos not from a hadron decay ($\nu$), prompt electrons and muons ($\ell$) "dressed" by adding close-by ($\Delta R < 0.1$) photons, and anti-$k_t$ $R = 0.4$ jets built from stable particles ($\tau > 30\,$ps) and tau leptons excluding neutrinos and prompt dressed muons. Jets are $b$-tagged ($b$-jet) using ghost-matched $b$-hadrons with $p_{\text{T}} > 5\,$GeV. Apart from the kinematic requirements, isolation and overlap removal criteria are applied. Jets within $\Delta R = 0.4$ of a photon are removed if the $p_{\text{T}}$ of charged particles within $\Delta R = 0.3$ of the photon is smaller than $10\,\%$ of its $p_{\text{T}}$. Jets within $\Delta R = 0.4$ of a lepton are removed. Events are removed where a photon is close ($\Delta R < 0.4$) to a lepton or a surviving jet. The measured fiducial particle-level cross section is $\sigma_{tq\gamma}\times\mathcal{B}(t\rightarrow l\nu b)+\sigma_{t(\rightarrow l\nu b\gamma)q} = 303\pm 9(\text{stat.})^{+33}_{-32}(\text{syst.})\,$fb.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and photon isolation cone radius $R=0.4$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and photon isolation cone radius $R=0.4$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and photon isolation cone radius $R=0.4$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and photon isolation cone radius $R=0.2$.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and isolation cone radius $0.2$ at NNLO QCD.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ at NNLO QCD.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $|\eta^{\gamma}|<0.6$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.6<|\eta^{\gamma}|<0.8$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.8<|\eta^{\gamma}|<1.37$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.56<|\eta^{\gamma}|<1.81$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.81<|\eta^{\gamma}|<2.01$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $2.01<|\eta^{\gamma}|<2.37$.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $|\eta^{\gamma}|<0.6$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.6<|\eta^{\gamma}|<0.8$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.8<|\eta^{\gamma}|<1.37$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.56<|\eta^{\gamma}|<1.81$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.81<|\eta^{\gamma}|<2.01$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $2.01<|\eta^{\gamma}|<2.37$ at NNLO QCD.