$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 0.5$ in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $0.5 \leq R_{\mathrm{T}} < 1.5$ in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $1.5 \leq R_{\mathrm{T}} < 2.5$ in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $2.5 \leq R_{\mathrm{T}} < 5$ in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 5$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 0.5$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $0.5 \leq R_{\mathrm{T}} < 1.5$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $1.5 \leq R_{\mathrm{T}} < 2.5$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $2.5 \leq R_{\mathrm{T}} < 5$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 5$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 0.5$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $0.5 \leq R_{\mathrm{T}} < 1.5$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $1.5 \leq R_{\mathrm{T}} < 2.5$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $2.5 \leq R_{\mathrm{T}} < 5$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 5$ in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 0.5$ in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $0.5 \leq R_{\mathrm{T}} < 1.5$ in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $1.5 \leq R_{\mathrm{T}} < 2.5$ in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $2.5 \leq R_{\mathrm{T}} < 5$ in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 5$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 0.5$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $0.5 \leq R_{\mathrm{T}} < 1.5$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $1.5 \leq R_{\mathrm{T}} < 2.5$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $2.5 \leq R_{\mathrm{T}} < 5$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 5$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 0.5$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $0.5 \leq R_{\mathrm{T}} < 1.5$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $1.5 \leq R_{\mathrm{T}} < 2.5$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $2.5 \leq R_{\mathrm{T}} < 5$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 5$ in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 0.5$ in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $0.5 \leq R_{\mathrm{T}} < 1.5$ in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $1.5 \leq R_{\mathrm{T}} < 2.5$ in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $2.5 \leq R_{\mathrm{T}} < 5$ in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 5$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 0.5$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $0.5 \leq R_{\mathrm{T}} < 1.5$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $1.5 \leq R_{\mathrm{T}} < 2.5$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $2.5 \leq R_{\mathrm{T}} < 5$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 5$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 0.5$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $0.5 \leq R_{\mathrm{T}} < 1.5$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $1.5 \leq R_{\mathrm{T}} < 2.5$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $2.5 \leq R_{\mathrm{T}} < 5$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 0.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $0.5 \leq R_{\mathrm{T}} < 1.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $1.5 \leq R_{\mathrm{T}} < 2.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $2.5 \leq R_{\mathrm{T}} < 5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 0.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $0.5 \leq R_{\mathrm{T}} < 1.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $1.5 \leq R_{\mathrm{T}} < 2.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $2.5 \leq R_{\mathrm{T}} < 5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 0.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $0.5 \leq R_{\mathrm{T}} < 1.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $1.5 \leq R_{\mathrm{T}} < 2.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ transverse momentum spectrum for events with $2.5 \leq R_{\mathrm{T}} < 5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 0.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $0.5 \leq R_{\mathrm{T}} < 1.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $1.5 \leq R_{\mathrm{T}} < 2.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $2.5 \leq R_{\mathrm{T}} < 5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 0.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $0.5 \leq R_{\mathrm{T}} < 1.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $1.5 \leq R_{\mathrm{T}} < 2.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $2.5 \leq R_{\mathrm{T}} < 5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 0.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $0.5 \leq R_{\mathrm{T}} < 1.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $1.5 \leq R_{\mathrm{T}} < 2.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ transverse momentum spectrum for events with $2.5 \leq R_{\mathrm{T}} < 5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 0.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $0.5 \leq R_{\mathrm{T}} < 1.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $1.5 \leq R_{\mathrm{T}} < 2.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $2.5 \leq R_{\mathrm{T}} < 5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 0.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $0.5 \leq R_{\mathrm{T}} < 1.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $1.5 \leq R_{\mathrm{T}} < 2.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $2.5 \leq R_{\mathrm{T}} < 5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $0 \leq R_{\mathrm{T}} < 0.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $0.5 \leq R_{\mathrm{T}} < 1.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $1.5 \leq R_{\mathrm{T}} < 2.5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ transverse momentum spectrum for events with $2.5 \leq R_{\mathrm{T}} < 5$ normalised to the $R_{\mathrm{T}}$-integrated spectrum in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{K}^{+}+\mathrm{K}^{-})/(\pi^{+}+\pi^{-})$ for events in the $0 \leq R_{\mathrm{T}} < 5$ class in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{K}^{+}+\mathrm{K}^{-})/(\pi^{+}+\pi^{-})$ for events in the $0 \leq R_{\mathrm{T}} < 0.5$ class in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{K}^{+}+\mathrm{K}^{-})/(\pi^{+}+\pi^{-})$ for events in the $0.5 \leq R_{\mathrm{T}} < 1.5$ class in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{K}^{+}+\mathrm{K}^{-})/(\pi^{+}+\pi^{-})$ for events in the $1.5 \leq R_{\mathrm{T}} < 2.5$ class in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{K}^{+}+\mathrm{K}^{-})/(\pi^{+}+\pi^{-})$ for events in the $2.5 \leq R_{\mathrm{T}} < 5$ class in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{K}^{+}+\mathrm{K}^{-})/(\pi^{+}+\pi^{-})$ for events in the $0 \leq R_{\mathrm{T}} < 5$ class in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{K}^{+}+\mathrm{K}^{-})/(\pi^{+}+\pi^{-})$ for events in the $0 \leq R_{\mathrm{T}} < 0.5$ class in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{K}^{+}+\mathrm{K}^{-})/(\pi^{+}+\pi^{-})$ for events in the $0.5 \leq R_{\mathrm{T}} < 1.5$ class in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{K}^{+}+\mathrm{K}^{-})/(\pi^{+}+\pi^{-})$ for events in the $1.5 \leq R_{\mathrm{T}} < 2.5$ class in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{K}^{+}+\mathrm{K}^{-})/(\pi^{+}+\pi^{-})$ for events in the $2.5 \leq R_{\mathrm{T}} < 5$ class in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{K}^{+}+\mathrm{K}^{-})/(\pi^{+}+\pi^{-})$ for events in the $0 \leq R_{\mathrm{T}} < 5$ class in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{K}^{+}+\mathrm{K}^{-})/(\pi^{+}+\pi^{-})$ for events in the $0 \leq R_{\mathrm{T}} < 0.5$ class in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{K}^{+}+\mathrm{K}^{-})/(\pi^{+}+\pi^{-})$ for events in the $0.5 \leq R_{\mathrm{T}} < 1.5$ class in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{K}^{+}+\mathrm{K}^{-})/(\pi^{+}+\pi^{-})$ for events in the $1.5 \leq R_{\mathrm{T}} < 2.5$ class in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{K}^{+}+\mathrm{K}^{-})/(\pi^{+}+\pi^{-})$ for events in the $2.5 \leq R_{\mathrm{T}} < 5$ class in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{p} + \mathrm {\overline{p}})/(\pi^{+} + \pi^{-})$ for events in the $0 \leq R_{\mathrm{T}} < 5$ class in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{p} + \mathrm {\overline{p}})/(\pi^{+} + \pi^{-})$ for events in the $0 \leq R_{\mathrm{T}} < 0.5$ class in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{p} + \mathrm {\overline{p}})/(\pi^{+} + \pi^{-})$ for events in the $0.5 \leq R_{\mathrm{T}} < 1.5$ class in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{p} + \mathrm {\overline{p}})/(\pi^{+} + \pi^{-})$ for events in the $1.5 \leq R_{\mathrm{T}} < 2.5$ class in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{p} + \mathrm {\overline{p}})/(\pi^{+} + \pi^{-})$ for events in the $2.5 \leq R_{\mathrm{T}} < 5$ class in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{p} + \mathrm {\overline{p}})/(\pi^{+} + \pi^{-})$ for events in the $0 \leq R_{\mathrm{T}} < 5$ class in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{p} + \mathrm {\overline{p}})/(\pi^{+} + \pi^{-})$ for events in the $0 \leq R_{\mathrm{T}} < 0.5$ class in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{p} + \mathrm {\overline{p}})/(\pi^{+} + \pi^{-})$ for events in the $0.5 \leq R_{\mathrm{T}} < 1.5$ class in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{p} + \mathrm {\overline{p}})/(\pi^{+} + \pi^{-})$ for events in the $1.5 \leq R_{\mathrm{T}} < 2.5$ class in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{p} + \mathrm {\overline{p}})/(\pi^{+} + \pi^{-})$ for events in the $2.5 \leq R_{\mathrm{T}} < 5$ class in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{p} + \mathrm {\overline{p}})/(\pi^{+} + \pi^{-})$ for events in the $0 \leq R_{\mathrm{T}} < 5$ class in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{p} + \mathrm {\overline{p}})/(\pi^{+} + \pi^{-})$ for events in the $0 \leq R_{\mathrm{T}} < 0.5$ class in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{p} + \mathrm {\overline{p}})/(\pi^{+} + \pi^{-})$ for events in the $0.5 \leq R_{\mathrm{T}} < 1.5$ class in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{p} + \mathrm {\overline{p}})/(\pi^{+} + \pi^{-})$ for events in the $1.5 \leq R_{\mathrm{T}} < 2.5$ class in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-differential $(\mathrm{p} + \mathrm {\overline{p}})/(\pi^{+} + \pi^{-})$ for events in the $2.5 \leq R_{\mathrm{T}} < 5$ class in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ average transverse momentum as a function of $R_{\mathrm{T}}$ in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ average transverse momentum as a function of $R_{\mathrm{T}}$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\pi^{+}+\pi^{-}$ average transverse momentum as a function of $R_{\mathrm{T}}$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ average transverse momentum as a function of $R_{\mathrm{T}}$ in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ average transverse momentum as a function of $R_{\mathrm{T}}$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{K}^{+}+\mathrm{K}^{-}$ average transverse momentum as a function of $R_{\mathrm{T}}$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ average transverse momentum as a function of $R_{\mathrm{T}}$ in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ average transverse momentum as a function of $R_{\mathrm{T}}$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$\mathrm{p} + \mathrm {\overline{p}}$ average transverse momentum as a function of $R_{\mathrm{T}}$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-integrated $(\mathrm{K}^{+} + \mathrm{K}^{-})/(\pi^{+} + \pi^{-})$ as a function of $R_{\mathrm{T}}$ in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-integrated $(\mathrm{K}^{+} + \mathrm{K}^{-})/(\pi^{+} + \pi^{-})$ as a function of $R_{\mathrm{T}}$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-integrated $(\mathrm{K}^{+} + \mathrm{K}^{-})/(\pi^{+} + \pi^{-})$ as a function of $R_{\mathrm{T}}$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-integrated $(\mathrm{p} + \mathrm {\overline{p}})/(\pi^{+} + \pi^{-})$ as a function of $R_{\mathrm{T}}$ in the Toward region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-integrated $(\mathrm{p} + \mathrm {\overline{p}})/(\pi^{+} + \pi^{-})$ as a function of $R_{\mathrm{T}}$ in the Away region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

$p_{\mathrm{T}}$-integrated $(\mathrm{p} + \mathrm {\overline{p}})/(\pi^{+} + \pi^{-})$ as a function of $R_{\mathrm{T}}$ in the Transverse region in pp collisions at $\sqrt{s} = 13~\mathrm{TeV}$.

The ratio of data to the sum of the SM backgrounds. The uncertainties of simulated data are only the statistical unvertainty in the simulation predictions.

Distribution of $p_{T}^{miss}$ after the preselection from 2017 data (black points) and simulation (colored lines). The simulated distribution of two signal points are represented by colored lines, while not being stacked on the background distributions: $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (500, 490) and (500, 420) GeV. The last bin includes the overflow events.

The ratio of data to the sum of the SM backgrounds. The uncertainties of simulated data are only the statistical unvertainty in the simulation predictions.

Distribution of $p_{T}^{miss}$ after the preselection from 2018 data (black points) and simulation (colored lines). The simulated distribution of two signal points are represented by colored lines, while not being stacked on the background distributions: $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (500, 490) and (500, 420) GeV. The last bin includes the overflow events.

The ratio of data to the sum of the SM backgrounds. The uncertainties of simulated data are only the statistical unvertainty in the simulation predictions.

Distribution of $N_{Jet}$ after the preselection from 2017 data (black points) and simulation (colored lines). The simulated distribution of two signal points are represented by colored lines, while not being stacked on the background distributions: $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (500, 490) and (500, 420) GeV. The last bin includes the overflow events.

The ratio of data to the sum of the SM backgrounds. The uncertainties of simulated data are only the statistical unvertainty in the simulation predictions.

Distribution of $N_{Jet}$ after the preselection from 2018 data (black points) and simulation (colored lines). The simulated distribution of two signal points are represented by colored lines, while not being stacked on the background distributions: $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (500, 490) and (500, 420) GeV. The last bin includes the overflow events.

The ratio of data to the sum of the SM backgrounds. The uncertainties of simulated data are only the statistical unvertainty in the simulation predictions.

Simulated distribution of $p_{T}(l)$ after the preselection. The total background distribution and the signal distributions are all normalized to unit area. The signal distributions are given for a top squark mass of 300 GeV and $\Delta {m}=10$, 30,50, and 80 GeV.

Simulated distribution of $p_{T}(l)$ after the preselection. The total background distribution and the signal distributions are all normalized to unit area. The signal distributions are shown for four different $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ values, all corresponding to the same $\Delta {m} = 30$ GeV.

Simulated distribution of $p_{T}^{miss}$ after the preselection. The total background distribution and the signal distributions are all normalized to unit area. The signal distributions are given for a top squark mass of 300 GeV and $\Delta {m}=10$, 30,50, and 80 GeV.

Simulated distribution of $p_{T}^{miss}$ after the preselection. The total background distribution and the signal distributions are all normalized to unit area. The signal distributions are shown for four different $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ values, all corresponding to the same $\Delta {m} = 30$ GeV.

Simulated distribution of $N_{Jet}$ after the preselection. The total background distribution and the signal distributions are all normalized to unit area. The signal distributions are given for a top squark mass of 300 GeV and $\Delta {m}=10$, 30,50, and 80 GeV.

Simulated distribution of $N_{Jet}$ after the preselection. The total background distribution and the signal distributions are all normalized to unit area. The signal distributions are shown for four different $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ values, all corresponding to the same $\Delta {m} = 30$ GeV.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 10$ GeV for 2017 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 20$ GeV for 2017 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 30$ GeV for 2017 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 40$ GeV for 2017 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 50$ GeV for 2017 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 60$ GeV for 2017 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 70$ GeV for 2017 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 80$ GeV for 2017 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 10$ GeV for 2018 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 20$ GeV for 2018 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 30$ GeV for 2018 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 40$ GeV for 2018 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 50$ GeV for 2018 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 60$ GeV for 2018 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 70$ GeV for 2018 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

BDT output distributions from data (black points) and simulation (colored lines) after the preselection for $\Delta {m} = 80$ GeV for 2018 data. The last bin corresponds to the SR. A representative $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ signal point is also presented, while not added to the SM background.

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

Distribution of $p_{T}(l)$ for $\Delta {m} = 10$ in the VR where $200 < p_{T}^{miss} < 280$ GeV from 2017 data (black points) and simulation (colored lines). The predicted signal distribution is shown for $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (475, 465).

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

Distribution of $p_{T}(l)$ for $\Delta {m} = 10$ in the VR where $200 < p_{T}^{miss} < 280$ GeV from 2018 data (black points) and simulation (colored lines). The predicted signal distribution is shown for $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (475, 465).

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

Distribution of BDT for $\Delta {m} = 10$ in the VR where $200 < p_{T}^{miss} < 280$ GeV from 2017 data (black points) and simulation (colored lines). The predicted signal distribution is shown for $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (475, 465).

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

Distribution of BDT for $\Delta {m} = 10$ in the VR where $200 < p_{T}^{miss} < 280$ GeV from 2018 data (black points) and simulation (colored lines). The predicted signal distribution is shown for $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (475, 465).

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

Distribution of $p_{T}^{miss}$ for $\Delta {m} = 60$ in the VR where $p_{T}(l) > 30$ GeV from 2017 data (black points) and simulation (colored lines). The predicted signal distribution is shown for $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (576, 516).

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

Distribution of $p_{T}^{miss}$ for $\Delta {m} = 60$ in the VR where $p_{T}(l) > 30$ GeV from 2018 data (black points) and simulation (colored lines). The predicted signal distribution is shown for $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (576, 516).

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

Distribution of BDT for $\Delta {m} = 60$ in the VR where $p_{T}(l) > 30$ GeV from 2017 data (black points) and simulation (colored lines). The predicted signal distribution is shown for $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (576, 516).

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

Distribution of BDT for $\Delta {m} = 60$ in the VR where $p_{T}(l) > 30$ GeV from 2018 data (black points) and simulation (colored lines). The predicted signal distribution is shown for $(m(\mathrm{\widetilde{t}}_{1}),m(\widetilde{\chi}^{0}_{1}))$ = (576, 516).

The ratio of data to the sum of the SM backgrounds. Only the statistical uncertainty in the data and the simulated background are represented.

Data yield and prediction of the total background and its composition in the eight SRs as defined in for the 2017 analysis. The predictions of W + jets, ttbar, and nonprompt lepton processes are based on data, while those of rare backgrounds are based on simulation. The predictions and their associated uncertainties are pre-fit. The uncertainties are the quadratic sum of the statistical and systematic uncertainties for background processes predicted from data, and the uncertainties on the cross section for the backgrounds predicted from simulation. The yields for two signal points, with $\Delta {m} = 10$ GeV and $\Delta {m} = 80$ GeV, are also shown. The bins corresponding to different $\Delta {m}$ trainings are not statistically independent.

The ratio of the number of observed events over the predicted total background per SR.

Data yield and prediction of the total background and its composition in the eight SRs as defined in for the 2018 analysis. The predictions of W + jets, ttbar, and nonprompt lepton processes are based on data, while those of rare backgrounds are based on simulation. The predictions and their associated uncertainties are pre-fit. The uncertainties are the quadratic sum of the statistical and systematic uncertainties for background processes predicted from data, and the uncertainties on the cross section for the backgrounds predicted from simulation. The yields for two signal points, with $\Delta {m} = 10$ GeV and $\Delta {m} = 80$ GeV, are also shown. The bins corresponding to different $\Delta {m}$ trainings are not statistically independent.

The ratio of the number of observed events over the predicted total background per SR.

Exclusion limit at 95% CL for the four-body decay of the top squark as a function of $m(\mathrm{\widetilde{t}}_{1})$ and $\Delta {m}$, for the combined data of the years 2016, 2017, and 2018. The lines represent the observed limits. These limits are derived using the expected top squark pair production cross section. The green line represents the central values, and the other lines the variations due to the theoretical or experimental uncertainties.

Exclusion limit at 95% CL for the four-body decay of the top squark as a function of $m(\mathrm{\widetilde{t}}_{1})$ and $\Delta {m}$, for the combined data of the years 2016, 2017, and 2018. The lines represent the expected limits. These limits are derived using the expected top squark pair production cross section. The green line represents the central values, and the other lines the variations due to the theoretical or experimental uncertainties.

Exclusion limit at 95% CL for the four-body decay of the top squark as a function of $m(\mathrm{\widetilde{t}}_{1})$ and $\Delta {m}$, for the combined data of the years 2016, 2017, and 2018.

Reconstruction level differential cross section spectla, obtained for the central acceptance, $|\eta| < 3$, for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology compared to the to the EPOS-LHC predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

The number of high purity tracksin the first $\eta$ bin after a gap of $4.5 \leq \Delta\eta^F < 5.0$ for events with the Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

The transeverse momentum of high purity tracksin the first $\eta$ bin after a gap of $4.5 \leq \Delta\eta^F < 5.0$ for events with the Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

The total energy of all Particle Flow candidatesin the first $\eta$ bin after a gap of $4.5 \leq \Delta\eta^F < 5.0$ for events with the Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV

Unfolded diffraction enhanced differential cross section for events with Pomeron-Lead ($\mathrm{I\!P}\mathrm{Pb}$) topology, defined on the region $|\eta| < 5.2$, compared to hadron level predictions of the MC generators. The data are corrected for the contribution from events with undetectable energy in the HF calorimeter adjacent to the rapidity gap. The corrections are obtained using the EPOS-LHC MC sample. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffraction enhanced differential cross section for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology, defined on the region $|\eta| < 5.2$, compared to hadron level predictions of the MC generators. The data are corrected for the contribution from events with undetectable energy in the HF calorimeter adjacent to the rapidity gap. The corrections are obtained using the EPOS-LHC MC sample. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Lead ($\mathrm{I\!P}\mathrm{Pb}$) topology, compared to the to the hadron-level EPOS-LHC predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology, compared to the to the hadron-level EPOS-LHC predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Lead ($\mathrm{I\!P}\mathrm{Pb}$) topology, compared to the to the hadron-level QGSJET II predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Unfolded diffractive enhanced differential cross section spectla for events with Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology, compared to the to the hadron-level QGSJET II predictions, broken down into the non-diffractive (ND), central diffractive (CD), single diffractive (SD) and double diffractive (DD) components. Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

Reconstruction level diffraction enhanced differential cross section spectrum for Pomeron-Proton ($\mathrm{I\!P}\mathrm{p} + \gamma \mathrm{p}$) topology corrected for the contribution from events with undetectable energy in the HF calorimeter adjacent to the rapidity gap, compared with the spectrum obtained with events satisfying the ZDC veto requirement $E_{ZDC−} < 1 \mathrm{~TeV}$ which selects only the events without lead nuclear break up (without correction for HF undetectable energy). Forward Rapidity Gap definition: $|\eta| < 2.5$: $p_{T}^{track} < 200$ MeV and $\sum \limits_{bin} E^{PF} < 6$ GeV $|\eta| \in [2.5,3.0]$: $\sum \limits_{bin} E_{neutral}^{PF} < 13.4$ GeV $|\eta| > 3.0$: No particles

$K^{0}_{S} K^{0}_{S}$ correlation meassurement in $10-20\%$ centrality.

$K^{0}_{S} K^{0}_{S}$ correlation meassurement in $20-30\%$ centrality.

$K^{0}_{S} K^{0}_{S}$ correlation meassurement in $30-40\%$ centrality.

$K^{0}_{S} K^{0}_{S}$ correlation meassurement in $40-50\%$ centrality.

$K^{0}_{S} K^{0}_{S}$ correlation meassurement in $50-60\%$ centrality.

$\Lambda K^{0}_{S} \oplus \overline{\Lambda} K^{0}_{S}$ correlation meassurement in $0-80\%$ centrality.

$\Lambda\Lambda\oplus\overline{\Lambda}\overline{\Lambda}$ correlation meassurement in $0-80\%$ centrality.

The extracted $R_{inv}$ from $K^{0}_{S} K^{0}_{S}$ correlations as a function of centrality.

The extracted $\lambda$ parameters from $K^{0}_{S} K^{0}_{S}$ correlations as a function of centrality.

Extracted effective range ($d_{0}$) and real components of the complex scattering length ($\Re f_{0}$) from the $\Lambda\Lambda\oplus\overline{\Lambda}\overline{\Lambda}$ and $\Lambda K^{0}_{S} \oplus \overline{\Lambda} K^{0}_{S}$ correlations.

Extracted imaginary components of the complex scattering length ($\Im f_{0}$) and real components of the complex scattering length ($\Re f_{0}$) from the $\Lambda K^{0}_{S} \oplus \overline{\Lambda} K^{0}_{S}$ correlation.

Electron $p_{T}$ distribution using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

$p_{T}^{miss}$ distribution in the electron channel using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

Muon $p_T$ distribution using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

$p_T^{miss}$ distribution in the muon channel using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

$M_T$ distribution for the E+INVISIBLE system using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

$M_T$ distribution for the MU+INVISIBLE system using 2016-2018 data sets. The complete set of selection criteria were applied. Two examples of SSM WPRIME boson signal samples with masses of 3.8 and 5.6 TeV are shown. The lower panel shows the ratio of data to SM prediction, and the shaded band represents the systematic uncertainties.

Theoretical prediction for the production and decay cross-section of SSM WPRIME bosons (W --> LEPTON + NU) as a function of the WPRIME mass, for the range 200-6400 GeV and inclusive in phase-space (no cuts). The uncertainty shown covers PDF and ALPHAS uncertainties. These theoretical predictions are obtained with PYTHIA (LO) and corrected to NNLO in QCD using FEWZ program.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the product of WPRIME cross-section production and BF(WPRIME --> LEPTON + NU) for an SSM WPRIME boson, as a function of the boson mass, for the electron decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits. The theoretical predictions for the SSM at NNLO precision in QCD are shown with narrow gray bands corresponding to the PDF and ALPHAS uncertainties.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the product of WPRIME cross-section production and BF(WPRIME --> LEPTON + NU) for an SSM WPRIME boson, as a function of the boson mass, for the muon decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits. The theoretical predictions for the SSM at NNLO precision in QCD are shown with narrow gray bands corresponding to the PDF and ALPHAS uncertainties.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the product of WPRIME cross-section production and BF(WPRIME --> LEPTON + NU) for an SSM WPRIME boson, as a function of the boson mass, for the combined electron and muon decay channels. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits. The theoretical predictions for the SSM at NNLO precision in QCD are shown with narrow gray bands corresponding to the PDF and ALPHAS uncertainties.

The 95% CL observed (solid line) and expected (dashed line) model independent cross section upper limits as function of the $MT_{min}$ threshold, for the electron decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits.

The 95% CL observed (solid line) and expected (dashed line) model independent cross section upper limits as function of the $MT_{min}$ threshold, for the muon decay channel. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits.

The 95% CL observed (solid line) and expected (dashed line) model independent cross section upper limits as function of the $MT_{min}$ threshold, for the combined electron and muon decay channels. The shaded bands represent the one (green) and two (yellow) standard deviation uncertainty bands for the expected limits.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the coupling strength ratio, gWprime/gW as a function of the mass of the WPRIME boson, for the electron channel. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. The dotted line represents the case of SSM coupling, gWprime, being equal to gW, the SM coupling.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the coupling strength ratio, gWprime/gW, as a function of the mass of the WPRIME boson, for the muon channel. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. The dotted line represents the case of SSM coupling, gWprime , being equal to gW, the SM coupling.

Exclusion limits in the 2D plane (1/R, $\mu$) for the split-UED interpretation for the n = 2 case, for the electron channel for the 2016-2018 data sets. R is the radius of the additional spatial dimension and $\mu$ the bulk mass for the fermion field in five dimensions. The expected limit is depicted as a black dashed line. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The experimentally excluded region is the entire area to the left of the solid black line.

Exclusion limits in the 2D plane (1/R, $\mu$) for the split-UED interpretation for the n = 2 case, for the muon channel for the 2016-2018 data sets. R is the radius of the additional spatial dimension and $\mu$ the bulk mass for the fermion field in five dimensions. The expected limit is depicted as a black dashed line. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The experimentally excluded region is the entire area to the left of the solid black line.

Exclusion limits in the 2D plane (1/R, $\mu$) for the split-UED interpretation for the n = 2 case, for the combined electron and muon channels for the 2016-2018 data sets. R is the radius of the additional spatial dimension and $\mu$ the bulk mass for the fermion field in five dimensions. The expected limit is depicted as a black dashed line. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The experimentally excluded region is the entire area to the left of the solid black line.

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the ${\lambda}_{231}$ coupling in the RPV SUSY model with a STAU mediator, as a function of the mass of STAU, for the electron channel and for the production coupling, ${\lambda}^{'}_{3ij}$=0.5 . The couplings ${\lambda}^{'}_{3ij}$, ${\lambda}_{231}$, and ${\lambda}_{132}$ are defined in Fig. 1. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. Limits for other values of ${\lambda}^{'}_{3ij}$ are obtained multiplying the ones shown by the ratio ${\lambda}^{'}_{3ij}$/0.5

Observed (solid line) and expected (dashed line) upper limits at 95% CL on the ${\lambda}_{132}$ coupling in the RPV SUSY model with a stau mediator, as a function of the mass of the stau, for the muon channel and for the production coupling, ${\lambda}^{'}_{3ij}$ = 0.5 . The couplings ${\lambda}^{'}_{3ij}$, ${\lambda}_{231}$, and ${\lambda}_{132}$ are defined in Fig. 1. The one (green) and two (yellow) standard deviation uncertainty bands for the expected limits are shown. The area above the limit curve is excluded. Limits for other values of ${\lambda}^{'}_{3ij}$ are obtained multiplying these ones by the ratio ${\lambda}^{'}_{3ij}$/0.5

Region in the oblique (W,Y) parameter phase space allowed by the current analysis at 95% CL. The combination of the electron and muon channel distributions from the 2017 and 2018 data sets at sqrt(s) = 13 TeV is used, having 101 fb-1 of luminosity.

Negative Log-Likelihood for the determination of the oblique W parameter, for the combination of the electron and muon channel and for 2017 and 2018 data sets (L=101 fb-1) at sqrt(s) = 13 TeV.

Total background estimate and observed data in all bins used in the limit calculation for the SSM WPRIME limits. The total event sample is split into the three years of data taking and the two lepton flavors.The distributions being shown are $M_T$ ones. The associated integrated luminosities for each year are of 35.9, 41.5, and 59.4 fb$̣^{-1}. The mass bins are numbered in ascending order.

Differential cross sections for the EW W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW + QCD W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW +QCD W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW + QCD W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW + QCD W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW + QCD W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW + QCD W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

The $m_{W\gamma}$ distribution for muon events (electron events are also used and the muon events are shown here as an example.) satisfying the aQGC region selection and used to set constraints on the anomalous gauge coupling parameters (left). The gray line represents a nonzero $f_{M,2}/Λ^4$ setting. Events with $m_{W\gamma} > 1500 GeV$ are included in the last bin.

Exclusion limits at the 95% CL for each aQGC coefficient, derived from the mW$\gamma$ distribution, assuming all other coefficients are set to zero. Unitarity bounds corresponding to each operator are also listed. All coupling parameter limits are in TeV^{−4}, while Ubound values are in TeV.

Results summary.

Jet-particle $v_{2}$ as a function of $p_{\mathrm{T}}$ in high-multiplicity (0$-$10%) p$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV, the $p_{\mathrm{T}}$ of associated particles is larger than 1.5 GeV/$\it{c}$

Inclusive charged-particle $v_{2}$ as a function of $p_{\mathrm{T}}$ in 20$-$60% semi-central Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV

Jet-particle $v_{2}$ as a function of $p_{\mathrm{T}}$ in high-multiplicity (20$-$60%) Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV, the $p_{\mathrm{T}}$ of associated particles is larger than 2 GeV/$\it{c}$

Jet-particle $v_{2}$ as a function of $p_{\mathrm{T}}$ in high-multiplicity (20$-$60%) Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV, the $p_{\mathrm{T}}$ of associated particles is larger than 3 GeV/$\it{c}$

Jet-particle $v_{2}$ as a function of $p_{\mathrm{T}}$ in semi-central (20$-$60%) Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV, the $p_{\mathrm{T}}$ of associated particles is larger than 5 GeV/$\it{c}$

The profile likelihood scan as a function of Bs to mu+mu- and B0 to mu+mu- decay branching fractions in 2D.

Feynman diagrams of the effective neutralino pair production in the GMSB simplified model in which the two neutralinos decay into two gravitinos ($\tilde{G}$) and two $Z$ bosons (left), a $Z$ and a Higgs boson ($H$) (center), or two Higgs bosons (right).

The distributions of the most impactful input variables to the TD jet tagger for signal (red, lighter) and collision background (blue, darker). They include the charged (upper left) and neutral (upper right) hadron energy fractions, the number of track constituents in the jet (middle left), the $\Delta R$ between the jet axis and the closest track associated with the PV (middle right), and the jet time (lower).

The distributions of the most impactful input variables to the TD jet tagger for signal (red, lighter) and collision background (blue, darker). They include the charged (upper left) and neutral (upper right) hadron energy fractions, the number of track constituents in the jet (middle left), the $\Delta R$ between the jet axis and the closest track associated with the PV (middle right), and the jet time (lower).

The distributions of the most impactful input variables to the TD jet tagger for signal (red, lighter) and collision background (blue, darker). They include the charged (upper left) and neutral (upper right) hadron energy fractions, the number of track constituents in the jet (middle left), the $\Delta R$ between the jet axis and the closest track associated with the PV (middle right), and the jet time (lower).

The distributions of the most impactful input variables to the TD jet tagger for signal (red, lighter) and collision background (blue, darker). They include the charged (upper left) and neutral (upper right) hadron energy fractions, the number of track constituents in the jet (middle left), the $\Delta R$ between the jet axis and the closest track associated with the PV (middle right), and the jet time (lower).

The distributions of the most impactful input variables to the TD jet tagger for signal (red, lighter) and collision background (blue, darker). They include the charged (upper left) and neutral (upper right) hadron energy fractions, the number of track constituents in the jet (middle left), the $\Delta R$ between the jet axis and the closest track associated with the PV (middle right), and the jet time (lower).

TD jet tagger score distributions (left) for signal (red, lighter) and collision background (blue, darker). Identification probability for the signal versus the misidentification probability for the background (right) with the tagger working point (w.~p.) used in the analysis shown as a blue marker. Both are evaluated using an independent sample of testing events.

TD jet tagger score distributions (left) for signal (red, lighter) and collision background (blue, darker). Identification probability for the signal versus the misidentification probability for the background (right) with the tagger working point (w.~p.) used in the analysis shown as a blue marker. Both are evaluated using an independent sample of testing events.

The efficiency of the TD jet tagger working point used in the analysis is shown as a function of the lab frame transverse decay length for simulated signals with $\chi$ mass of 400 GeV. The uncertainties shown account for lifetime dependence and statistical uncertainty.

The TD jet tagger score distributions for simulation (shaded histogram) and data (black markers) when using electrons from $W\to e\nu_e$ events as proxy objects for signal jets. The histograms and data points have been normalized to unit area. The last bin contains jets with tagger scores greater than 0.996, the threshold used to tag signal jets. Similar levels of agreement are observed for photon proxy objects from the $Z\to\ell^+\ell^-\gamma$ sample.

The TD jet tagger misidentification probability measured using the nominal $W$+jets MR (black round markers) is shown along with the systematic uncertainty (gray band), quantifying the degree of process dependence measured from alternative MRs. The measurements in the alternative MRs are displayed as well ($Z$+jets MR as green round markers, $t\bar{t}$ MR as red squared markers, QCD MR as blue triangular markers) along with their respective statistical uncertainty. On the left, this probability is shown for the first 19.9 fb$^{-1}$ of data collected in 2016, while on the right it is shown for the last 16.4 fb$^{-1}$ of data collected in 2016combined with data collected in 2017-2018.

The TD jet tagger misidentification probability measured using the nominal $W$+jets MR (black round markers) is shown along with the systematic uncertainty (gray band), quantifying the degree of process dependence measured from alternative MRs. The measurements in the alternative MRs are displayed as well ($Z$+jets MR as green round markers, $t\bar{t}$ MR as red squared markers, QCD MR as blue triangular markers) along with their respective statistical uncertainty. On the left, this probability is shown for the first 19.9 fb$^{-1}$ of data collected in 2016, while on the right it is shown for the last 16.4 fb$^{-1}$ of data collected in 2016combined with data collected in 2017-2018.

Distribution of the number of TD tagged jets for the $m_{\chi} = 400$ GeVsimulated signal samples with $c\tau_{\chi} = 0.5$ m (solid red line) and $c\tau_{\chi} = 3.0$ m (dotted green line), estimated background (blue square markers), and data (black round markers). The signal distributions are normalized to the expected cross section limit. The blue shaded region indicates the systematic uncertainty in the background prediction. No background prediction is shown for the bin with zero TD tagged jets as it is the main control region used to predict the background for the other two bins. There are zero observed events in the bin with two or more TD tagged jets.

Expected and observed 95% CL upper limits on $\sigma_{\chi\chi}$ as functions of $m_\chi$ in a scenario with $\mathcal{B}(\chi\to HG) = 0.5$ and $c\tau = 0.5$ m (left) or 3 m (right).

Expected and observed 95% CL upper limits on $\sigma_{\chi\chi}$ as functions of $m_\chi$ in a scenario with $\mathcal{B}(\chi\to HG) = 0.5$ and $c\tau = 0.5$ m (left) or 3 m (right).

Expected and observed 95% CL upper limits on $\sigma_{\chi\chi}$ as functions of $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to H\tilde{G}) = 0.5$ and $m_{\chi} = 400$ GeV (left) or 1000 GeV (right).

Expected and observed 95% CL upper limits on $\sigma_{\chi\chi}$ as functions of $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to H\tilde{G}) = 0.5$ and $m_{\chi} = 400$ GeV (left) or 1000 GeV (right).

The observed 95% CL upper limit on $\sigma_{\chi\chi}$ as a function of $m_{\chi}$ and $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to H\tilde{G}) = 0.5$. The area enclosed by the dotted black line corresponds to the observed excluded region.

The distribution of the jet charged hadron energy fraction, a variable used as input to the TD jet tagger score, for simulation (shaded histogram) and data (black markers) when using electrons from $W\to e\nu_e$ events as proxy objects for signal jets. The histograms and data points have been normalized to unit area. Similar levels of agreement are observed for photon proxy objects from the $Z\to\ell^+\ell^-\gamma$ sample.

The distribution of the jet neutral hadron energy fraction, a variable used as input to the TD jet tagger score, for simulation (shaded histogram) and data (black markers) when using electrons from $W\to e\nu_e$ events as proxy objects for signal jets. The histograms and data points have been normalized to unit area. Similar levels of agreement are observed for photon proxy objects from the $Z\to\ell^+\ell^-\gamma$ sample.

The distribution of the number of track constituents in the jet, a variable used as input to the TD jet tagger score, for simulation (shaded histogram) and data (black markers) when using electrons from $W\to e\nu_e$ events as proxy objects for signal jets. The histograms and data points have been normalized to unit area. Similar levels of agreement are observed for photon proxy objects from the $Z\to\ell^+\ell^-\gamma$ sample.

The $\eta$ distribution of TD-tagged jets in a background-enriched control region that comprises events with exactly one TD-tagged jet. Observed data (black round markers) and the corresponding prediction based on control samples in data (empty squared markers), measured using the nominal $W$+jets MR, are compared. The prediction uncertainty (gray band) includes the systematic uncertainty quantifying the degree of process dependence measured from alternative MRs. The predictions for the shape and the normalization of the $\eta$ distribution are consistent with the data.

Jet time distribution in a sample of b-tagged jets from dilepton $t \bar{t}$ events in 2017 data-taking period (black round markers) and simulation (filled histogram). A Gaussian smearing procedure is applied to the jet time in the $t \bar{t}$ sample (green line) to correct for effects that are difficult to simulate (timing drift caused by crystals transparency loss due to detector aging, electronics jitter).

The observed 95% CL upper limit on $\sigma_{\chi\chi}$ as a function of $m_{\chi}$ and $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to H\tilde{G}) = 1$. The area enclosed by the dotted black line corresponds to the observed excluded region.

The observed 95% CL upper limit on $\sigma_{\chi\chi}$ as a function of $m_{\chi}$ and $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to H\tilde{G}) = 0.75$, $\mathcal{B}(\chi\to Z\tilde{G}) = 0.25$. The area enclosed by the dotted black line corresponds to the observed excluded region.

The observed 95% CL upper limit on $\sigma_{\chi\chi}$ as a function of $m_{\chi}$ and $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to H\tilde{G}) = 0.25$, $\mathcal{B}(\chi\to Z\tilde{G}) = 0.75$. The area enclosed by the dotted black line corresponds to the observed excluded region.

The observed 95% CL upper limit on $\sigma_{\chi\chi}$ as a function of $m_{\chi}$ and $c\tau_{\chi}$ in a scenario with $\mathcal{B}(\chi\to Z\tilde{G}) = 1$. The area enclosed by the dotted black line corresponds to the observed excluded region.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 127 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 127 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 127 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 150 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 150 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 150 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 175 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 175 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 175 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 200 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 200 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 200 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 250 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 250 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 250 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 300 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 300 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 300 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 400 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 400 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 400 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 600 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 600 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 600 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 800 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 800 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 800 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1000 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1000 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1000 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1250 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1250 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1250 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1500 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1500 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1500 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1800 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the merged topology, namely, the H (or Z) decay products are produced with an angular separation $\Delta R < 0.8$, and the H (or Z) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 3\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1800 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with exactly one quark in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and only one b-quark (or light quark) has $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 5\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

The efficiency of identifying a LLP decay as a TD-tagged jet in bins of the LLP transverse and longitudinal decay position. The sample used to compute the efficiency contains events with pair production of $\chi$ with a mass of 1800 GeV and a lifetime of 0.5 and 3 m, and considering the combinations of the $\chi$ decay modes considered in this search ($H \tilde{G} \rightarrow b\bar{b} \tilde{G}$ or $Z\tilde{G} \rightarrow q\bar{q} \tilde{G}$). The efficiency is calculated on top of the acceptance definition for the resolved topology with two quarks in acceptance, namely, the H (or Z) decay products are produced with an angular separation $\Delta R \geq 0.8$, and both b-quarks (or light quarks) have $p_T > 30$ GeV and $|\eta|<1$. The full simulation signal yield prediction can be reproduced within 7\%. This nonclosure uncertainty is added in quadrature to the statistical uncertainty of each bin.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 127 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 127 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 150 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 150 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 175 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 175 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 200 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 200 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 250 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 250 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 300 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 300 GeV.

Cutflow table for a $\tilde{\chi}_{1}^{0}$ signal sample with a mass of 400 GeV.