Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{1})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}^{\mathrm{add.}}_{2})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{2})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{add.}})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\Delta\mathrm{R}_{\mathrm{b}\mathrm{b}}^{\mathrm{avg}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}_{3})|$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}_{3})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}_{3})$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}_{3})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}_{4})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}_{4})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{b}}_{\mathrm{T}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{b}}_{\mathrm{T}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $|\eta(\mathrm{b}^{\mathrm{extra}}_{1})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{1})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}^{\mathrm{extra}}_{2})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{2})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{extra}})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $H^{\mathrm{j}}_{\mathrm{T}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{j}}_{\mathrm{T}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\mathrm{m}_{\mathrm{b}\mathrm{b}}^{\mathrm{max}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{light}}_{\mathrm{T}}$ in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $H^{\mathrm{light}}_{\mathrm{T}}$ in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $N_{\mathrm{jets}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $N_{\mathrm{jets}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $N_{\mathrm{b}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{add.}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\Delta\mathrm{R}_{\mathrm{b}\mathrm{b}}^{\mathrm{avg}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}_{4})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}_{4})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{extra}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\mathrm{m}_{\mathrm{b}\mathrm{b}}^{\mathrm{max}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $N_{\mathrm{jets}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $N_{\mathrm{jets}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $N_{\mathrm{b}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{add.}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\Delta\mathrm{R}_{\mathrm{b}\mathrm{b}}^{\mathrm{avg}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}_{4})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}_{4})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{extra}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\mathrm{m}_{\mathrm{b}\mathrm{b}}^{\mathrm{max}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $N_{\mathrm{jets}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $N_{\mathrm{jets}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $N_{\mathrm{b}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Fiducial cross sections from the measurements of all observables, compared to predictions from different ttbb simulation approaches. For each of the normalized differential measurements the fiducial cross section in the respective phase space is also determined. In the paper only one representative observable is quoted for each fiducial phase space, while here the measured cross section with the uncertainties from the fit to the respective observable is summarized.

Postfit distributions of $S_\mathrm{T}^\mathrm{MET}$ in the $\mathrm{e}\mu$ channel of the $\geq$1b category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width. The last bin includes the overflow.

Postfit distributions of $S_\mathrm{T}^\mathrm{MET}$ in the $\ell\tau_\mathrm{h}$ channel of the 0b category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width. The last bin includes the overflow.

Postfit distributions of $S_\mathrm{T}^\mathrm{MET}$ in the $\ell\tau_\mathrm{h}$ channel of the $\geq$1b category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width. The last bin includes the overflow.

Postfit distributions of $S_\mathrm{T}^\mathrm{MET}$ in the $\tau_\mathrm{h}\tau_\mathrm{h}$ channel of the 0b category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width. The last bin includes the overflow.

Postfit distributions of $S_\mathrm{T}^\mathrm{MET}$ in the $\tau_\mathrm{h}\tau_\mathrm{h}$ channel of the $\geq$1b category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width. The last bin includes the overflow.

Postfit distributions of $\chi$ in the $\mathrm{e}\mu$ channel of the 0j, $400 < m_\mathrm{vis} < 600\,\mathrm{GeV}$ category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width.

Postfit distributions of $\chi$ in the $\mathrm{e}\mu$ channel of the 0j, $m_\mathrm{vis} > 600\,\mathrm{GeV}$ category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width.

Postfit distributions of $\chi$ in the $\ell\tau_\mathrm{h}$ channel of the 0j, $400 < m_\mathrm{vis} < 600\,\mathrm{GeV}$ category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width.

Postfit distributions of $\chi$ in the $\ell\tau_\mathrm{h}$ channel of the 0j, $m_\mathrm{vis} > 600\,\mathrm{GeV}$ category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width.

Postfit distributions of $\chi$ in the $\tau_\mathrm{h}\tau_\mathrm{h}$ channel of the 0j, $400 < m_\mathrm{vis} < 600\,\mathrm{GeV}$ category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width.

Postfit distributions of $\chi$ in the $\tau_\mathrm{h}\tau_\mathrm{h}$ channel of the 0j, $m_\mathrm{vis} > 600\,\mathrm{GeV}$ category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width.

Exclusion limits on the total cross section of a scalar LQ as a function of LQ mass under the assumption of exclusive LQ couplings to b quarks and $\tau$ leptons.

Exclusion limits on the total cross section of a vector LQ as a function of LQ mass under the assumption of exclusive LQ couplings to b quarks and $\tau$ leptons.

Exclusion limits on the mass of a scalar LQ as a function of the coupling strength $\lambda$ under the assumption of exclusive LQ couplings to b quarks and $\tau$ leptons.

Exclusion limits on the mass of a vector LQ as a function of the coupling strength $\lambda$ under the assumption of exclusive LQ couplings to b quarks and $\tau$ leptons.

Exclusion limits on the coupling strength $\lambda$ of a scalar LQ as a function of the LQ mass under the assumption of exclusive LQ couplings to b quarks and $\tau$ leptons.

Exclusion limits on the coupling strength $\lambda$ of a vector LQ as a function of the LQ mass under the assumption of exclusive LQ couplings to b quarks and $\tau$ leptons.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet no-btag category with the BDT requirements selecting the most signal-like events. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet no-btag category with the BDT requirements selecting the most signal-like events. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet no-btag category with the BDT requirements selecting the most signal-like events. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to b quarks and $\tau$ leptons, using $\lambda = 1.5$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to light-flavor quarks and $\tau$ leptons, using $\lambda = 1.5$.

Expected upper limit at 95% CL on the coupling strength $\lambda$ of a scalar LQ to b quarks and $\tau$ leptons.

Expected upper limit at 95% CL on the coupling strength $\lambda$ of a scalar LQ to light-flavor quarks and $\tau$ leptons.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet btag category with the intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet no-btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet no-btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet no-btag category with low BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet no-btag category with low BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet no-btag category with low BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet no-btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet no-btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to b quarks and $\tau$ leptons, using $\lambda = 1.0$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to b quarks and $\tau$ leptons, using $\lambda = 2.0$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to b quarks and $\tau$ leptons, using $\lambda = 3.0$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to light-flavor quarks and $\tau$ leptons, using $\lambda = 0.5$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to light-flavor quarks and $\tau$ leptons, using $\lambda = 1.0$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to light-flavor quarks and $\tau$ leptons, using $\lambda = 2.0$.

Observed and expected distributions of the jet multiplicity in the e+jet btag control region.

Observed and expected distributions of the jet multiplicity in the $\mu$+jet btag control region.

Inclusive cross section predictions up to QCD NNLO accuracy using the NNPDF3.1 PDF set

Measured production cross sections ratio $\sigma(W^+ + \overline{c})$ / $\sigma(W^- + c)$

Predictions for the cross sections ratio $\sigma(W^+ + \overline{c})$ / $\sigma(W^- + c)$ at QCD NLO accuracy from MCFM using different PDF sets

Predictions for the cross sections ratio $\sigma(W^+ + \overline{c})$ / $\sigma(W^- + c)$ up to QCD NNLO accuracy using the NNPDF3.1 PDF set

Measured diferential cross sections $\sigma(W^- + c) + \sigma(W^+ + \overline{c})$ as a function of the absolute value of the pseudorapidity of the lepton from the W decay, unfolded to the particle level.

Measured diferential cross sections as a function of the transverse momentum of the lepton from the W decay, unfolded to the particle level.

Measured diferential cross sections $\sigma(W^- + c) + \sigma(W^+ + \overline{c})$ as a function of the absolute value of the pseudorapidity of the lepton from the W decay, unfolded to the parton level.

Measured diferential cross sections as a function of the transverse momentum of the lepton from the W decay, unfolded to the parton level.

Measured diferential cross sections ratio $R=\sigma(W^+ + \overline{c}) / \sigma(W^- + c)$ as a function of the absolute value of the pseudorapidity of the lepton from the W decay.

Measured diferential cross sections ratio $R=\sigma(W^+ + \overline{c}) / \sigma(W^- + c)$ as a function of the transverse momentum of the lepton from the W decay.

Two dimensional observed likelihood profile scan as a function of mu_ggH and mu_qqH

Observed $2 \Delta NLL$ of a 2D scan for $c^{-}_{\varphi Q}$ and $c_{\varphi t}$ with the other WCs profiled.

Observed $2 \Delta NLL$ of a 2D scan for $c_{tG}$ and $c_{t\varphi}$ with the other WCs profiled.

Observed $2 \Delta NLL$ of a 2D scan for $c^{1}_{Qt}$ and $c^{8}_{Qt}$ with the other WCs profiled.

Observed $2 \Delta NLL$ of a 2D scan for $c^{1}_{QQ}$ and $c^{8}_{Qt}$ with the other WCs profiled.

Observed $2 \Delta NLL$ of a 2D scan for $c^{1}_{Qt}$ and $c^{1}_{tt}$ with the other WCs profiled.

Observed $2 \Delta NLL$ of a 2D scan for $c^{1}_{QQ}$ and $c^{1}_{Qt}$ with the other WCs profiled.

The nuclear modification factor $\mathrm{R_{AA}}$ versus $p_{\mathrm{T}}$ for prompt $\Lambda^+_c$ production in centrality regions of 0-90, 0-10, 10-30, 30-50 and 50-90% in PbPb collisions. The global uncertainty includes the uncertainties for the luminosity of pp collisions, number of MB events in PbPb collisions, and tracking efficiency. The global uncertainty for $\mathrm{R_{AA}}$ is 16.5%.

The ratio of the production cross sections of prompt $\Lambda^+_c$ to prompt $\mathrm{D_0}$ versus $p_{\mathrm{T}}$ from pp collisions. The global normalization uncertainty is 6.6%.

The ratio of the $\mathrm{T_{AA}}$-scaled prompt $\Lambda^+_c$ yields to prompt $\mathrm{D_0}$ versus $p_{\mathrm{T}}$ from PbPb collisions for 0-90 and 0-10% centrality classes. The global normalization uncertainty is 7.3%.

$\Delta\phi$ projection of balance function in low $p_{T}$ in 0-10% centrality

$\Delta\phi$ projection of balance function in low $p_{T}$ in 30-40% centrality

1D $\Delta\eta$ projection of balance function in low $p_{T}$ in 70-80% centrality

$\Delta\eta$ projection of balance function in low $p_{T}$ in 0-40 $N_{trk}$

$\Delta\eta$ projection of balance function in low $p_{T}$ in 120-150 $N_{trk}$

$\Delta\eta$ projection of balance function in low $p_{T}$ in 270-300 $N_{trk}$

$\Delta\phi$ projection of balance function in low $p_{T}$ in 0-40 $N_{trk}$

$\Delta\phi$ projection of balance function in low $p_{T}$ in 120-150 $N_{trk}$

$\Delta\phi$ projection of balance function in low $p_{T}$ in 270-300 $N_{trk}$

The width of balance function in $<|\Delta\eta|>$ in low $p_{T}$ for PbPb collisions

Balance function width in low $p_{T}$ PbPb relative $<|\Delta\eta|>/<|\Delta\eta|>_{N_{ch}<65}$

The width of balance function in $<|\Delta\eta|>$ in low $p_{T}$ for pPb collisions

Balance function width in low $p_{T}$ in pPb relative $<|\Delta\eta|>/<|\Delta\eta|>_{N_{ch}<24}$

The width of the balance function in $<|\Delta\phi|>$ in low $p_{T}$ for PbPb collisions

Balance function width in low $p_{T}$ in PbPb for relative $<|\Delta\phi|>/<|\Delta\phi|>_{N_{ch}<65}$

The width of the balance function in $<|\Delta\phi|>$ in low $p_{T}$ for pPb collisions

Balance function width in low $p_{T}$ in pPb for relative $|\Delta\phi|/<|\Delta\phi|>_{N_{ch}<24}$

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in pPb for 0-40 $N_{trk}$

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in pPb for 120-150 $N_{trk}$

Balance function projection as a function of $\Delta\eta$ in intermediate $p_{T}$ in pPb for 270-300 $N_{trk}$

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in pPb for $N_{trk}$ 0-40

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in pPb for $N_{trk}$ 120-150

Balance function projection as a function of $\Delta\phi$ in intermediate $p_{T}$ in pPb for $N_{trk}$ 270-300

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for 0-10% centrality

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for 30-40% centrality

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for 70-80% centrality

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in pPb for $N_{trk}$ 0-40

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in pPb for $N_{trk}$ 120-150

Balance function projection as a function of $\Delta\eta$ in high $p_{T}$ in pPb for $N_{trk}$ 270-300

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for $N_{trk}$ 0-40

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for $N_{trk}$ 120-150

Balance function projection as a function of $\Delta\phi$ in high $p_{T}$ in PbPb for $N_{trk}$ 270-300

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for low $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for high $p_{T}$

Average width of the balance function in $\Delta \phi$ with $p_{T}$ for low $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for high $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for low $p_{T}$

Average width of the balance function in $\Delta \eta$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for high $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for low $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for intermediate $p_{T}$

Average width of the balance function in $\Delta \phi$ with $N_{ch}$ for high $p_{T}$

Feynman diagrams of $\mathrm{Z'}\to\mu^{-}\mu^{+}$ with a $\mathrm{Z'}$ boson produced via $\mathrm{s}\overline{\mathrm{b}}\to\mathrm{Z'}$, with one $\mathrm{b}$ quark and one $\overline{\mathrm{s}}$ quark in the final state.

Distribution of $\min(m_{\mu\mathrm{b}})$ as obtained from simulation in events with $N_{\mathrm{b}} \geq 1$ passing all the other selection requirements. In this search, we require $\min(m_{\mu\mathrm{b}})>175~\mathrm{GeV}$. The stacked histogram displays the expected distribution from the simulation of the SM backgrounds, while the overlaid open histograms illustrate the size and shape of the $\mathrm{Z'}$ contribution from the LFU model described in Eq. (1), for several $\mathrm{Z'}$ mass hypotheses. For illustrative purposes, we choose couplings $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$. The contribution of background processes other than DY and $\mathrm{t}\overline{\mathrm{t}}$ is so small that it is only barely visible at the bottom of the stacked histogram. The hatched region indicates the statistical uncertainty arising from the limited size of the SM simulated samples. Histograms are normalized to unit area.

Distributions of $m_{\mu\mu}$ in the $N_{\mathrm{b}}=1$ event category. The stacked histogram displays the expected distribution from the SM background simulation. The overlaid open distributions illustrate the $\mathrm{Z'}$ contribution from the LFU model at $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$ for a variety of $\mathrm{Z'}$ mass hypotheses. The observed data are shown as black points with statistical error bars. The hatched region indicates the statistical uncertainty arising from the limited size of the SM simulated samples. The size of the bins increases as a function of $m_{\mu\mu}$. In extracting the results of the search, the background is estimated directly from data, so the SM background simulation is only illustrative in these distributions.

Distributions of $m_{\mu\mu}$ in the $N_{\mathrm{b}}\geq 2$ event category. The stacked histogram displays the expected distribution from the SM background simulation. The overlaid open distributions illustrate the $\mathrm{Z'}$ contribution from the LFU model at $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$ for a variety of $\mathrm{Z'}$ mass hypotheses. The observed data are shown as black points with statistical error bars. The hatched region indicates the statistical uncertainty arising from the limited size of the SM simulated samples. The size of the bins increases as a function of $m_{\mu\mu}$. In extracting the results of the search, the background is estimated directly from data, so the SM background simulation is only illustrative in these distributions.

Invariant mass $m_{\mu\mu}$ distributions in the $N_{\mathrm{b}} = 1$ category, shown together with the corresponding selected background functional forms used as input to the $discrete~profiling$ method when probing the $m_{\mathrm{Z'}}=500~\mathrm{GeV}$ hypothesis. The expected signal distribution for the LFU model described in Eq. (1), with couplings $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$, is overlaid. The displayed mass range corresponds to the fit window used for this $m_{\mathrm{Z'}}$ hypothesis, which is $\pm 10\, \sigma_{\mathrm{mass}}$ around the probed $m_{\mathrm{Z'}}$ value. While the likelihood fits are performed on unbinned data, here we present the data in binned histograms with binning chosen to reflect the size of $\sigma_{\mathrm{mass}}$.

Invariant mass $m_{\mu\mu}$ distributions in the $N_{\mathrm{b}} \geq 2$ category, shown together with the corresponding selected background functional forms used as input to the $discrete~profiling$ method when probing the $m_{\mathrm{Z'}}=500~\mathrm{GeV}$ hypothesis. The expected signal distribution for the LFU model described in Eq. (1), with couplings $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$, is overlaid. The displayed mass range corresponds to the fit window used for this $m_{\mathrm{Z'}}$ hypothesis, which is $\pm 10\, \sigma_{\mathrm{mass}}$ around the probed $m_{\mathrm{Z'}}$ value. While the likelihood fits are performed on unbinned data, here we present the data in binned histograms with binning chosen to reflect the size of $\sigma_{\mathrm{mass}}$.

Exclusion limits at $95\%$ CL on the number of selected BSM events with $N_{\mathrm{b}} \geq 1$ as functions of $m_{\mathrm{Z'}}$ for $f_{2\mathrm{b}}=0$. The quantity $f_{2\mathrm{b}}$ is the fraction of BSM events passing the analysis selection that have at least two $\mathrm{b}$ quark jets. The solid black (dashed red) curve represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing $68~(95)\%$ of the distribution of limits expected under the background-only hypothesis.

Exclusion limits at $95\%$ CL on the number of selected BSM events with $N_{\mathrm{b}} \geq 1$ as functions of $m_{\mathrm{Z'}}$ for $f_{2\mathrm{b}}=0.25$. The quantity $f_{2\mathrm{b}}$ is the fraction of BSM events passing the analysis selection that have at least two $\mathrm{b}$ quark jets. The solid black (dashed red) curve represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing $68~(95)\%$ of the distribution of limits expected under the background-only hypothesis.

Exclusion limits at $95\%$ CL on the number of selected BSM events with $N_{\mathrm{b}} \geq 1$ as functions of $m_{\mathrm{Z'}}$ for $f_{2\mathrm{b}}=0.75$. The quantity $f_{2\mathrm{b}}$ is the fraction of BSM events passing the analysis selection that have at least two $\mathrm{b}$ quark jets. The solid black (dashed red) curve represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing $68~(95)\%$ of the distribution of limits expected under the background-only hypothesis.

Exclusion limits at $95\%$ CL on the number of selected BSM events with $N_{\mathrm{b}} \geq 1$ as functions of $m_{\mathrm{Z'}}$ for $f_{2\mathrm{b}}=1$. The quantity $f_{2\mathrm{b}}$ is the fraction of BSM events passing the analysis selection that have at least two $\mathrm{b}$ quark jets. The solid black (dashed red) curve represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing $68~(95)\%$ of the distribution of limits expected under the background-only hypothesis.

Observed (solid) and median expected (dashed) exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution, marked by the dotted curves. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed (solid) and median expected (dashed) exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution, marked by the dotted curves. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The solid black (dashed red) curves represent the observed (median expected) exclusions. The dotted curves denote the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and dotted curves is (expected to be) excluded. The dotted curve for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$ lies beyond the displayed $|g_{\mathrm{Z'}}|$ range and is, therefore, not shown. The shaded blue area represents the region preferred from the global fit in JHEP 04 (2023) 033 at $95\%$ CL. The region above the green dash-dotted curve is incompatible at $95\%$ CL with the measurement of the mass difference between the mass eigenstates of the neutral $\mathrm{B}_\mathrm{s}$ mesons.

Observed exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$.

Expected exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The solid black (dashed red) curves represent the observed (median expected) exclusions. The dotted curves denote the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and dotted curves is (expected to be) excluded. The shaded blue area represents the region preferred from the global fit in JHEP 04 (2023) 033 at $95\%$ CL. The region above the green dash-dotted curve is incompatible at $95\%$ CL with the measurement of the mass difference between the mass eigenstates of the neutral $\mathrm{B}_\mathrm{s}$ mesons.

Observed exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$.

Expected exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The solid black (dashed red) curves represent the observed (median expected) exclusions. The dotted curves denote the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and dotted curves is (expected to be) excluded. The shaded blue area represents the region preferred from the global fit in JHEP 04 (2023) 033 at $95\%$ CL. The region above the green dash-dotted curve is incompatible at $95\%$ CL with the measurement of the mass difference between the mass eigenstates of the neutral $\mathrm{B}_\mathrm{s}$ mesons.

Observed exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$.

Expected exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The solid black (dashed red) curves represent the observed (median expected) exclusions. The dotted curves denote the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and dotted curves is (expected to be) excluded. The shaded blue area represents the region preferred from the global fit in JHEP 04 (2023) 033 at $95\%$ CL. The region above the green dash-dotted curve is incompatible at $95\%$ CL with the measurement of the mass difference between the mass eigenstates of the neutral $\mathrm{B}_\mathrm{s}$ mesons.

Observed exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$.

Expected exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{Z'}}|$--$m_{\mathrm{Z'}}$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for a fixed value of $\theta_{23}=0$. The solid black (dashed red) curve represents the observed (median expected) exclusion. The dotted curve denotes the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and the dotted curves is (expected to be) excluded. The shaded blue area represents the $|g_{\mathrm{Z'}}|$ range preferred from a global fit in JHEP 04 (2023) 033 at $95\%$ CL.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{Z'}}|$--$m_{\mathrm{Z'}}$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for a fixed value of $\theta_{23}=0$.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{Z'}}|$--$m_{\mathrm{Z'}}$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for a fixed value of $\theta_{23}=0$.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{Z'}}|$--$m_{\mathrm{Z'}}$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for a fixed value of $\theta_{23}=0$. The coupling values are shown at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid.

Reconstructed mass of the Z' candidate in CR1, which consists of flavor sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Pre-fit (post-fit) results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR1, which consists of flavor sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Pre-fit (post-fit) results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR1, which consists of flavor sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Pre-fit (post-fit) results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR2, which consists of the dilepton mass sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Post-fit dielectron (dimuon) channel results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR2, which consists of the dilepton mass sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Post-fit dielectron (dimuon) channel results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR2, which consists of the dilepton mass sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Post-fit dielectron (dimuon) channel results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR2, which consists of the dilepton mass sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Post-fit dielectron (dimuon) channel results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR2, which consists of the dilepton mass sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Post-fit dielectron (dimuon) channel results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR2, which consists of the dilepton mass sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Post-fit dielectron (dimuon) channel results are shown on the left (right). The fitting procedure is described in Section 7.

Distributions of the reconstructed Z' candidate mass in SR1 (upper row), SR2 (middle row) and SR3 (lower row). Left (right) column corresponds to the dielectron (dimuon) channel. Signal samples with ($m_{Z'}$, $m_{N}$) = (4000 GeV, 200 GeV)$ and ($m_{Z'}$, $m_{N}$) = (4000 GeV, 1200 GeV)$ are plotted together with a cross section times branching fraction of 1 fb

Distributions of the reconstructed Z' candidate mass in SR1 (upper row), SR2 (middle row) and SR3 (lower row). Left (right) column corresponds to the dielectron (dimuon) channel. Signal samples with ($m_{Z'}$, $m_{N}$) = (4000 GeV, 200 GeV)$ and ($m_{Z'}$, $m_{N}$) = (4000 GeV, 1200 GeV)$ are plotted together with a cross section times branching fraction of 1 fb

Distributions of the reconstructed Z' candidate mass in SR1 (upper row), SR2 (middle row) and SR3 (lower row). Left (right) column corresponds to the dielectron (dimuon) channel. Signal samples with ($m_{Z'}$, $m_{N}$) = (4000 GeV, 200 GeV)$ and ($m_{Z'}$, $m_{N}$) = (4000 GeV, 1200 GeV)$ are plotted together with a cross section times branching fraction of 1 fb

Distributions of the reconstructed Z' candidate mass in SR1 (upper row), SR2 (middle row) and SR3 (lower row). Left (right) column corresponds to the dielectron (dimuon) channel. Signal samples with ($m_{Z'}$, $m_{N}$) = (4000 GeV, 200 GeV)$ and ($m_{Z'}$, $m_{N}$) = (4000 GeV, 1200 GeV)$ are plotted together with a cross section times branching fraction of 1 fb

Distributions of the reconstructed Z' candidate mass in SR1 (upper row), SR2 (middle row) and SR3 (lower row). Left (right) column corresponds to the dielectron (dimuon) channel. Signal samples with ($m_{Z'}$, $m_{N}$) = (4000 GeV, 200 GeV)$ and ($m_{Z'}$, $m_{N}$) = (4000 GeV, 1200 GeV)$ are plotted together with a cross section times branching fraction of 1 fb

Distributions of the reconstructed Z' candidate mass in SR1 (upper row), SR2 (middle row) and SR3 (lower row). Left (right) column corresponds to the dielectron (dimuon) channel. Signal samples with ($m_{Z'}$, $m_{N}$) = (4000 GeV, 200 GeV)$ and ($m_{Z'}$, $m_{N}$) = (4000 GeV, 1200 GeV)$ are plotted together with a cross section times branching fraction of 1 fb

The observed and expected 95\% CL upper limits on the product of \Zp boson cross section times branching fraction are shown for the case of ${m}_{N}$ = ${m}_{Z'}/4$ (upper row) and ${m}_{N}$ = 100 GeV (lower row) for dielectron (left column) and dimuon (right column) channels. The green and yellow bands indicate the 68\% and 95\% CL regions around the expected limit.

The observed and expected 95\% CL upper limits on the product of \Zp boson cross section times branching fraction are shown for the case of ${m}_{N}$ = ${m}_{Z'}/4$ (upper row) and ${m}_{N}$ = 100 GeV (lower row) for dielectron (left column) and dimuon (right column) channels. The green and yellow bands indicate the 68\% and 95\% CL regions around the expected limit.

The observed and expected 95\% CL upper limits on the product of \Zp boson cross section times branching fraction are shown for the case of ${m}_{N}$ = ${m}_{Z'}/4$ (upper row) and ${m}_{N}$ = 100 GeV (lower row) for dielectron (left column) and dimuon (right column) channels. The green and yellow bands indicate the 68\% and 95\% CL regions around the expected limit.

The observed and expected 95\% CL upper limits on the product of \Zp boson cross section times branching fraction are shown for the case of ${m}_{N}$ = ${m}_{Z'}/4$ (upper row) and ${m}_{N}$ = 100 GeV (lower row) for dielectron (left column) and dimuon (right column) channels. The green and yellow bands indicate the 68\% and 95\% CL regions around the expected limit.

Observed and expected exclusion regions at 95% CL in the 2D phase space of ${m}_{Z'}$ vs. ${m}_{N}$ for dielectron (left) and dimuon (right) channels. One standard deviation (s.d.) limits are also shown. Both opposite- and same-sign (OS and SS) dileptons from decays of heavy neutrino pairs are considered.

Observed and expected exclusion regions at 95% CL in the 2D phase space of ${m}_{Z'}$ vs. ${m}_{N}$ for dielectron (left) and dimuon (right) channels. One standard deviation (s.d.) limits are also shown. Both opposite- and same-sign (OS and SS) dileptons from decays of heavy neutrino pairs are considered.

The Levy scale parameter $R$ versus the transverse mass $m_{\mathrm{T}}$ in different centrality classes for negatively and positively charged hadron pairs.

The $1/R^2$ distribution vs. transverse mass $m_{\mathrm{T}}$ in different centrality classes for negatively and positively charged hadron pairs.

The linear-fit parameter $A$ versus $\langle N_{\text{part}}\rangle$ for positively and negatively charged hadron pairs. The bin range on $\langle N_{\text{part}}\rangle$ represents the 68% confidence interval associated to its uncertainty.

The linear-fit parameter $B$ versus $\langle N_{\text{part}}\rangle$ for positively and negatively charged hadron pairs. The bin range on $\langle N_{\text{part}}\rangle$ represents the 68% confidence interval associated to its uncertainty.

The Levy scale parameter $R$ versus $\langle N_{\text{part}}\rangle^{1/3}$ in different $m_{\mathrm{T}}$ classes for negatively and positively charged hadron pairs. The bin range on $\langle N_{\text{part}}\rangle^{1/3}$ represents the 68% confidence interval associated to its uncertainty.

The Levy stability index $\alpha$ versus the transverse mass $m_{\mathrm{T}}$ in different centrality classes for negatively and positively charged hadron pairs.

The average Levy stability index $\langle \alpha \rangle$ versus $\langle N_{\text{part}}\rangle$ for both positively and negatively charged hadron pairs. The bin range on $\langle N_{\text{part}}\rangle$ represents the 68% confidence interval associated to its uncertainty.

The correlation strength $\lambda$ versus the transverse mass $m_{\mathrm{T}}$ in different centrality classes for negatively and positively charged hadron pairs.

The correlation strength (re-scaled with the square of the pion fraction) $\lambda ^*$ versus the transverse mass $m_{\mathrm{T}}$ in different centrality classes for positively and negatively charged hadron pairs.

The $m_{e\mu}$ distribution of the observed data is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in red solid line, and the background component of the fit in red dotted line. Events and fit in each category are weighted by $S/(S+B)$. The one and two standard deviation uncertainty bands of the background component are shown in green and yellow. The lower panel shows the residuals after subtracting the background component of the fit from data.

Observed (expected) 95% confidence level upper limits on $\sigma(p p \to X(146) \to e \mu)$ for each individual analysis category (as shown in the left axis label) and for the combination of all analysis categories assuming the relative SM-like production cross sections of the ggH and VBF production modes.

The $m_{e\mu}$ distribution of the observed data in the category ggH cat 0 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category ggH cat 1 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category ggH cat 2 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category ggH cat 3 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category VBF cat 0 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category VBF cat 1 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

Best fit of $\sigma(p p \to X(146) \to e \mu)$ for each analysis category (black point) compared to the best fit for the combination of all analysis categories (blue line). The one standard deviation uncertainty of the per-category best fit is shown in red line and the one standard deviation uncertainty of the combined best fit is shown in green band.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.0.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.1.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.2.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.3.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.4.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.5.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.6.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.7.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.8.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.9.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 1.0.

The $v_{3}$ values as a function of centrality for prompt and nonprompt J/$\psi$ mesons.The kinematic range is 6.5 $< p_{\text{T}} <$ 50 GeV/c and $|y| < 2.4$.

The $v_{2}$ values of prompt J/$\psi$ mesons as a function of $p_{\text{T}}$ in the 10–60% centrality range. The results for 3 $< p_{\text{T}} <$ 6.5 and 6.5 $< p_{\text{T}} <$ 50 GeV/c are studied in the rapidity range of 1.6 $< |y| <$ 2.4 and $|y| <$ 2.4, respectively.

The $v_{2}$ values as functions of $p_{\text{T}}$ for prompt J/$\psi$ and prompt $\psi$(2S) mesons in the 10–60% centrality range. The results for $p_{\text{T}}$ < 6.5 and $p_{\text{T}}$ > 6.5 GeV/c are studied in the rapidity range of 1.6 $< |y| <$ 2.4 and $|y| <$ 2.4, respectively.

The $v_{2}$ values as functions of centrality for prompt J/$\psi$ and prompt $\psi$(2S) mesons. The kinematic range is 6.5 $< p_{\text{T}} <$ 50 GeV/c and $|y| <$ 2.4.

The $v_{3}$ values as functions of $p_{\text{T}}$ for prompt J/$\psi$ and prompt $\psi$(2S) mesons in the 10–60% centrality range. The results for $p_{\text{T}}$ < 6.5 and $p_{\text{T}}$ > 6.5 GeV/c are studied in the rapidity range of 1.6 $< |y| <$ 2.4 and $|y| <$ 2.4, respectively.

The $v_{3}$ values as functions of centrality for prompt J/$\psi$ and prompt $\psi$(2S) mesons. The kinematic range is 6.5 $< p_{\text{T}} <$ 50 GeV/c and $|y| <$ 2.4.

Measured min-$d_{xy}$ distribution in the 2-match category, after requiring the $d_{xy}$ muon to pass the isolation requirement $I_{\mathrm{rel}}^{\mathrm{PF}} <0.25$ (i. e., the B and D bins of the ABCD plane). Overlaid with a red histogram is the background predicted from the region of the ABCD plane failing the same requirement (the A and C bins), as well as three signal benchmark hypotheses (as defined in the legends), assuming $\alpha_D = \alpha_{\mathrm{EM}}$ (the fine-structure constant). The red hatched bands correspond to the background prediction uncertainty. The last bin includes the overflow.

Two-dimensional exclusion surfaces for $\Delta = 0.1 \, m_1$ as a function of the DM mass $m_1$ and the signal strength $y$, with $m_{A'} = 3 \, m_1$. Filled histograms denote observed limits on $\sigma(\mathrm{pp} \to A' \to \chi_2 \chi_1) \, \mathcal{B}(\chi_2 \to \chi_1 \mu^+ \mu^-)$. Solid (dashed) curves denote the observed (expected) exclusion limits at 95% CL, with 68% CL uncertainty bands around the expectation. Regions above the curves are excluded, depending on the $\alpha_D$ hypothesis: $\alpha_{\mathrm{D}} = \alpha_{\mathrm{EM}}$ (dark blue) or 0.1 (light magenta).

Two-dimensional exclusion surfaces for $\Delta = 0.4 \, m_1$ as a function of the DM mass $m_1$ and the signal strength $y$, with $m_{A'} = 3 \, m_1$. Filled histograms denote observed limits on $\sigma(\mathrm{pp} \to A' \to \chi_2 \chi_1) \, \mathcal{B}(\chi_2 \to \chi_1 \mu^+ \mu^-)$. Solid (dashed) curves denote the observed (expected) exclusion limits at 95% CL, with 68% CL uncertainty bands around the expectation. Regions above the curves are excluded, depending on the $\alpha_D$ hypothesis: $\alpha_{\mathrm{D}} = \alpha_{\mathrm{EM}}$ (dark blue) or 0.1 (light magenta).

The signal cross-sections values of the signal QBH RS1 with n=1 extra dimensions, as a function of threshold mass

The product of acceptance and selection efficiency for QBH ADD model, as a function of threshold mass of the signal.

The product of acceptance and selection efficiency for QBH RS1 models, as a function of threshold mass of the signal.

The product of acceptance and selection efficiency for q* signal, as a function of resonance mass

The product of acceptance and selection efficiency for b* signal, as a function of resonance mass of the signal.

Distributions of the γ+jet invariant mass generated from various signal samples for reconstructed mass 4 TeV

Fitting distribution of data with expected q*, and QBH signal events, and fit parameters p0= 3.05674e+00, p1 = 10.5343, p3 = 4.15163 and p4 = 0.107405

Fitting distribution of data with expected b* events, and fit parameters p0= 1.45055e+03, p1 = 2.04830e+01, p3 = -4.16760e-01 and p4 = -5.07768e-01

Fitting distribution of data with expected b* events for 0btag category, and fit parameters p0= 2.44120e+00, p1 = 1.08886e+01, p3 = 4.23198e+00 and p4 = 1.16300e-01

Exclusion limits on the product of cross section and branching fraction for q* signal coupling f = 1.0, as a function of the resonance mass hypothesis.

Theoretical cross-section of q* signal for coupling multiplier f =1.0

Exclusion limits on the product of cross section and branching fraction for the q* signal coupling f = 0.5, as a function of the resonance mass hypothesis.

Theoretical cross-section of q* signal for coupling multiplier f =0.5

Exclusion limits on the product of cross section and branching fraction for q* coupling f = 0.1, as a function of the resonance mass hypothesis.

Theoretical cross-section of q* signal for coupling multiplier f =0.1

Exclusion limits on the product of cross section and branching fraction for b* signal coupling f = 1.0, as a function of the resonance mass hypothesis.

Theoretical cross-section of b* signal for coupling multiplier f =1.0

Exclusion limits on the product of cross section and branching fraction for b* signal coupling f = 0.5, as a function of the resonance mass hypothesis.

Theoretical cross-section of b* signal for coupling multiplier f =0.5

Exclusion limits on the product of cross section and branching fraction for b* signal coupling f = 0.1, as a function of the resonance mass hypothesis.

Theoretical cross-section of b* signal for coupling multiplier f =0.1

Exclusion limits on the product of cross section and branching fraction for QBH ADD signal, as a function of the mass hypothesis.

Theoretical cross-section of ADD signal

Exclusion limits on the product of cross section and branching fraction for QBH RS1 signal, as a function of the mass hypothesis.

Theoretical cross-section of RS1 signal

The expected and observed exclusion mass limits for q* and b* signals, as a function coupling multiplier

The observed exclusion mass limits for b* signal, as a function of threshold mass of the signal.

Differential cross section measurements in bins of pT4l (v3)

Differential cross section measurements in bins of rapidity4l (v3)

Differential cross section measurements in bins of costheta1 (v3)

Differential cross section measurements in bins of costheta1 (v4 (2e2mu))

Differential cross section measurements in bins of costheta1 (v4 (4e+4mu))

Differential cross section measurements in bins of costheta2 (v3)

Differential cross section measurements in bins of costheta2 (v4 (2e2mu))

Differential cross section measurements in bins of costheta2 (v4 (4e+4mu))

Differential cross section measurements in bins of phi (v3)

Differential cross section measurements in bins of phi (v4 (2e2mu))

Differential cross section measurements in bins of phi (v4 (4e+4mu))

Differential cross section measurements in bins of phi1 (v3)

Differential cross section measurements in bins of phi1 (v4 (2e2mu))

Differential cross section measurements in bins of phi1 (v4 (4e+4mu))

Differential cross section measurements in bins of costhetastar (v3)

Differential cross section measurements in bins of costhetastar (v4 (2e2mu))

Differential cross section measurements in bins of costhetastar (v4 (4e+4mu))

Differential cross section measurements in bins of massZ1 (v3)

Differential cross section measurements in bins of massZ1 (v4 (2e2mu))

Differential cross section measurements in bins of massZ1 (v4 (4e+4mu))

Differential cross section measurements in bins of massZ2 (v3)

Differential cross section measurements in bins of massZ2 (v4 (2e2mu))

Differential cross section measurements in bins of massZ2 (v4 (4e+4mu))

Differential cross section measurements in bins of pTj1 (v3)

Differential cross section measurements in bins of pTHj (v3)

Differential cross section measurements in bins of mHj (v3)

Differential cross section measurements in bins of pTj2 (v3)

Differential cross section measurements in bins of mjj (v3)

Differential cross section measurements in bins of absdetajj (v3)

Differential cross section measurements in bins of dphijj (v3)

Differential cross section measurements in bins of pTHjj (v3)

Differential cross section measurements in bins of TCjmax (v3)

Differential cross section measurements in bins of TBjmax (v3)

Differential cross section measurements in bins of D0m (v3)

Differential cross section measurements in bins of D0m (v4 (2e2mu))

Differential cross section measurements in bins of D0m (v4 (4e+4mu))

Differential cross section measurements in bins of Dcp (v3)

Differential cross section measurements in bins of Dcp (v4 (2e2mu))

Differential cross section measurements in bins of Dcp (v4 (4e+4mu))

Differential cross section measurements in bins of D0hp (v3)

Differential cross section measurements in bins of D0hp (v4 (2e2mu))

Differential cross section measurements in bins of D0hp (v4 (4e+4mu))

Differential cross section measurements in bins of Dint (v3)

Differential cross section measurements in bins of Dint (v4 (2e2mu))

Differential cross section measurements in bins of Dint (v4 (4e+4mu))

Differential cross section measurements in bins of DL1 (v3)

Differential cross section measurements in bins of DL1 (v4 (2e2mu))

Differential cross section measurements in bins of DL1 (v4 (4e+4mu))

Differential cross section measurements in bins of DL1Zg (v3)

Differential cross section measurements in bins of DL1Zg (v4 (2e2mu))

Differential cross section measurements in bins of DL1Zg (v4 (4e+4mu))

Differential cross section measurements in bins of rapidity4l_pT4l (v3, bin 0)

Differential cross section measurements in bins of rapidity4l_pT4l (v3, bin 1)

Differential cross section measurements in bins of rapidity4l_pT4l (v3, bin 2)

Differential cross section measurements in bins of njets_pt30_eta4p7_pT4l (v3, bin 0)

Differential cross section measurements in bins of njets_pt30_eta4p7_pT4l (v3, bin 1)

Differential cross section measurements in bins of njets_pt30_eta4p7_pT4l (v3, bin 2)

Differential cross section measurements in bins of pTj1_pTj2 (v3)

Differential cross section measurements in bins of pT4l_pTHj (v3)

Differential cross section measurements in bins of massZ1_massZ2 (v3)

Differential cross section measurements in bins of TCjmax_pT4l (v3, bin 0)

Differential cross section measurements in bins of TCjmax_pT4l (v3, bin 1)

Differential cross section measurements in bins of TCjmax_pT4l (v3, bin 2)

Differential cross section measurements in bins of TCjmax_pT4l (v3, bin 3)

Correlations of higgs fiducial cross section in bins of absdetajj in v3.

Correlations of higgs fiducial cross section in bins of TBjmax in v3.

Correlations of higgs fiducial cross section in bins of massZ1 in v3.

Correlations of higgs fiducial cross section in bins of D0m in v3.

Correlations of higgs fiducial cross section in bins of DL1 in v3.

Correlations of higgs fiducial cross section in bins of mass4l in v3.

Correlations of higgs fiducial cross section in bins of pT4l in v3.

Correlations of higgs fiducial cross section in bins of phi in v3.

Correlations of higgs fiducial cross section in bins of pTHj in v3.

Correlations of higgs fiducial cross section in bins of pTj1 in v3.

Correlations of higgs fiducial cross section in bins of mHj in v3.

Correlations of higgs fiducial cross section in bins of pTHjj in v3.

Correlations of higgs fiducial cross section in bins of DL1Zg in v3.

Correlations of higgs fiducial cross section in bins of njets.

Correlations of higgs fiducial cross section in bins of pTj2 in v3.

Correlations of higgs fiducial cross section in bins of TCjmax in v3.

Correlations of higgs fiducial cross section in bins of rapidity4l in v3.

Correlations of higgs fiducial cross section in bins of costheta2 in v3.

Correlations of higgs fiducial cross section in bins of Dint in v3.

Correlations of higgs fiducial cross section in bins of Dcp in v3.

Correlations of higgs fiducial cross section in bins of phi1 in v3.

Correlations of higgs fiducial cross section in bins of costhetastar in v3.

Correlations of higgs fiducial cross section in bins of D0hp in v3.

Correlations of higgs fiducial cross section in bins of mjj in v3.

Correlations of higgs fiducial cross section in bins of massZ2 in v3.

Correlations of higgs fiducial cross section in bins of costheta1 in v3.

Correlations of higgs fiducial cross section in bins of dphijj in v3.

Correlations of higgs fiducial cross section in bins of massZ1 in v4.

Correlations of higgs fiducial cross section in bins of Dint in v4.

Correlations of higgs fiducial cross section in bins of DL1 in v4.

Correlations of higgs fiducial cross section in bins of D0m in v4.

Correlations of higgs fiducial cross section in bins of massZ2 in v4.

Correlations of higgs fiducial cross section in bins of DL1Zg in v4.

Correlations of higgs fiducial cross section in bins of Dcp in v4.

Correlations of higgs fiducial cross section in bins of D0hp in v4.

Correlations of higgs fiducial cross section in bins of phi in v4.

Correlations of higgs fiducial cross section in bins of costheta2 in v4.

Correlations of higgs fiducial cross section in bins of phi1 in v4.

Correlations of higgs fiducial cross section in bins of costheta1 in v4.

Correlations of higgs fiducial cross section in bins of costhetastar in v4.

Correlations of higgs fiducial cross section in bins of njets vs pt.

Correlations of higgs fiducial cross section in bins of rapidity4l vs pT4l in v3.

Correlations of higgs fiducial cross section in bins of pTj1 vs pTj2 in v3.

Correlations of higgs fiducial cross section in bins of pT4l vs pTHj in v3.

Correlations of higgs fiducial cross section in bins of massZ1 vs massZ2 in v3.

Correlations of higgs fiducial cross section in bins of TCjmax vs pT4l in v3.

Measured ratio of $\mathcal{B}_{4\mu}/\mathcal{B}_{2\mu}$

Measured branching fraction $\mathcal{B}_{4\mu}$

Measured branching fraction $\mathcal{B}_{4\mu}$

Measured $R_\text{AA}$ for the $\Upsilon(3S)$ state in the 0–90% centrality interval, integrated over the full kinematic range $p_\text{T}$ < 30 GeV/c and |y| < 2.4. The global uncertainty "PP MB" represents the pp luminosity and PbPb $N_\text{MB}$ combined uncertainties, whereas the global uncertainty "PP 3S" corresponds to the uncertainty on the$\Upsilon(3S)$ pp yields.

Measured $R_\text{AA}$ for the $\Upsilon(2S)$ state as functions of transverse momentum $p_\text{T}$, integrated over the rapidity range of |y| < 2.4 and the 0–90% centrality interval. The global uncertainty combines the uncertainties of $T_\text{AA}$ , pp luminosity, and PbPb $N_\text{MB}$.

Measured $R_\text{AA}$ for the $\Upsilon(3S)$ state as functions of transverse momentum $p_\text{T}$, integrated over the rapidity range of |y| < 2.4 and the 0–90% centrality interval. The global uncertainty combines the uncertainties of $T_\text{AA}$ , pp luminosity, and PbPb $N_\text{MB}$.

The double ratios of $\Upsilon(3S)$ /$\Upsilon(2S)$ as functions of PbPb collision centrality, integrated over the full kinematic range $p_\text{T}$ < 30 GeV/c and |y| < 2.4. The global uncertainty represents the combined systematic and statistical uncertainties from pp data.

The double ratios of $\Upsilon(3S)$ /$\Upsilon(2S)$ in the 0–90% centrality interval, integrated over the full kinematic range $p_\text{T}$ < 30 GeV/c and |y| < 2.4. The global uncertainty represents the combined systematic and statistical uncertainties from pp data.

The double ratios of $\Upsilon(3S)$ /$\Upsilon(2S)$ as functions of transverse momentum $p_\text{T}$, integrated over the rapidity range of |y| < 2.4 and the 0–90% centrality interval.

The total coherent $\mathrm{J}/\psi$ photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$, measured in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$ values used correspond to the center of each rapidity range. The theoretical uncertainties is due to the uncertainties in the photon flux.

The nuclear gluon suppression factor $R_{\mathrm{g}}^{\mathrm{Pb}}$ as a function of Bjorken $x$ extracted from the CMS measurement of the coherent $\mathrm{J}/\psi$ photoproduction in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $x$ values are evaluated at the centers of their corresponding rapidity ranges. The theoretical uncertainties are due to the uncertainties in the photon flux and the impulse approximation model.

The nuclear gluon suppression factor $R_{\mathrm{g}}^{\mathrm{Pb}}$ as a function of Bjorken $x$ extracted from the CMS measurement of the coherent $\mathrm{J}/\psi$ photoproduction in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $x$ values are evaluated at the centers of their corresponding rapidity ranges. The theoretical uncertainties are due to the uncertainties in the photon flux and the impulse approximation model.

The total covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The covariance matrix includes both the experimental and theoretical (photon flux) uncertainties. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The photon flux values from STARLight as a function of rapidity, in different neutron multiplicity classes: 0n0n, 0nXn, and XnXn.

The experimental covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The covariance matrix for the flux in the 0n0n neutron multiplicity class.

The theoretical (photon flux) covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The covariance matrix for the flux in the 0nXn neutron multiplicity class.

The covariance matrix for the flux in the XnXn neutron multiplicity class.

The total covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The covariance matrix includes both the experimental and theoretical (photon flux) uncertainties. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The experimental covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The theoretical (photon flux) covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

Comparison of the $m_{miss}$ shapes for the simulated signal events within the fiducial region and those outside it, after including the effect of PU protons as describe in the text, for a generated $m_{X}$ mass of 1000 GeV. The distributions are shown for single(+z)-single(−z) proton reconstruction categories.

Product of the acceptance times the combined reconstruction and identification efficiency, as a function of $m_{X}$ , for events generated inside the fiducial volume defined in Table 1. The curves shown panel display the different final states. The definition of the fiducial region and signal model used to estimate the acceptance is provided in the text.

Product of the acceptance times the combined reconstruction and identification efficiency, as a function of $m_{X}$ , for different proton reconstruction categories in the $Z \rightarrow \mu\mu $ analysis, and for events generated inside the fiducial volume defined in Table 1. The curves shown panel display the different final states. The definition of the fiducial region and signal model used to estimate the acceptance is provided in the text.

Distributions of the reconstructed proton $\xi$ in the negative arm from the proton mixing procedure with simulated MC events, are compared to data. The ee final states are shown.The ee events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the positive arm from the proton mixing procedure with simulated MC events, are compared to data. The ee final states are shown.The ee events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed di-proton rapidity from the proton mixing procedure with simulated MC events, are compared to data. The ee final states are shown.The ee events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the negative arm from the proton mixing procedure with simulated MC events, are compared to data. The $\mu\mu$ final states are shown.The $\mu\mu$ events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the positive arm from the proton mixing procedure with simulated MC events, are compared to data. The $\mu\mu$ final states are shown.The $\mu\mu$ events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed di-proton rapidity from the proton mixing procedure with simulated MC events, are compared to data. The $\mu\mu$ final states are shown.The $\mu\mu$ events are displayed without the Z boson $p_T$ requirement.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the negative arm from the proton mixing procedure with simulated MC events, are compared to data. The $\gamma$ final states are shown.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed proton $\xi$ in the positive arm from the proton mixing procedure with simulated MC events, are compared to data. The $\gamma$ final states are shown.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Distributions of the reconstructed di-proton rapidity from the proton mixing procedure with simulated MC events, are compared to data. The $\gamma$ final states are shown.For illustration, the simulated signal distributions are superimposed for various choices of $m_{X}$ , normalised to a generated fiducial cross section of 100 pb.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The proton $\xi$ distribution for the positive CT-PPS arm is shown.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The proton $\xi$ distribution for the negative CT-PPS arm is shown.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The di-proton invariant mass is shown.

Validation of the background modelling method, using the $e\mu$ control sample. Selected $e\mu$ events are mixed with protons from $Z\rightarrow \mu \mu$ events with $p_{T}(Z)<10$ GeV to simulate the combinatorial background shape, while the data points are unaltered $e\mu$ events. The $m_{miss}$ is shown.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the multi(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the multi(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the single(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZeeX$, with the protons reconstructed with the single(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the multi(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the multi(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the single(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow ppZ\mu\mu X$, with the protons reconstructed with the single(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 1 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the multi(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the multi(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the single(+z)-multi(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

$m_{miss}$ distributions in the $pp\rightarrow pp\gamma X$, with the protons reconstructed with the single(+z)-single(-z)method. The background distributions are shown after the fit. The expectations for a signal with $m_{X} = 1000$ GeV are superimposed, where the fiducial production cross section is normalised to 10 pb.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $Z\rightarrow ee$ analysis.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $Z\rightarrow \mu\mu$ analysis.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $Z$ analysis.

Upper limits on the $pp\rightarrow ppZ/\gamma X$ cross section at 95% CL, as a function of $m_X$. The 68 and 95% central intervals of the expected limits are represented by the green and yellow bands, respectively, while the observed limit is superimposed as a curve. The plot corresponds to the $\gamma$ analysis.

The distribution of the BDT discriminants for the 2016--2018 data set for three different categories in the combined single-electron and single-muon channels. The three categories are defined by the number of resolved t tags ($N_\textrm{RT}$), b tags ($N_\textrm{b}$), and jets ($N_\textrm{j}$), selected as representative based on their sensitivity to signal. Here, $\textrm{t}\bar{\textrm{t}} + \geq 1 \textrm{b}$ refers to $\textrm{t}\bar{\textrm{t}}$ events with at least one additional b jet, while $\textrm{t}\bar{\textrm{t}} + 0 \textrm{b}$ includes all other $\textrm{t}\bar{\textrm{t}}$ events not produced in association with a boson. The TOP grouping contains single top quark production along with the other $\textrm{t}\bar{\textrm{t}}$ processes not explicitly shown, and EW refers to events that contain W and Z bosons but no top quarks.

Measured signal strength ($\mu = \sigma_{\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}}/ \sigma_{\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}} ^{\textrm{SM}}$), corresponding cross section (in fb), and the expected and observed significance (in standard deviations) for $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ production from all analysis channels. Uncertainties shown in this table are statistical and systematic uncertainties respectively. In this table, SSDL\&ML and OSDL (2016) represent results from EPJC 80 (2020) 75 and JHEP 11 (2019) 082 respectively.

The distribution of the BDT discriminants for the 2016--2018 data set for three different categories in the combined single-electron and single-muon channels. The three categories are defined by the number of resolved t tags ($N_\textrm{RT}$), b tags ($N_\textrm{b}$), and jets ($N_\textrm{j}$), selected as representative based on their sensitivity to signal. Here, $\textrm{t}\bar{\textrm{t}} + \geq 1 \textrm{b}$ refers to $\textrm{t}\bar{\textrm{t}}$ events with at least one additional b jet, while $\textrm{t}\bar{\textrm{t}} + 0 \textrm{b}$ includes all other $\textrm{t}\bar{\textrm{t}}$ events not produced in association with a boson. The TOP grouping contains single top quark production along with the other $\textrm{t}\bar{\textrm{t}}$ processes not explicitly shown, and EW refers to events that contain W and Z bosons but no top quarks. The backgrounds and $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ signal (derived from the fit) are shown as a stacked histogram. The hatched bands correspond to the estimated total uncertainty after the fit. While the bins are shown to be equal width, they do not correspond to equal width in BDT value.

The distribution of the BDT discriminants for the full 2016--2018 data set in the all-hadronic channel. The SR category shown in this figure is defined by $N_\textrm{RT}=1$, $N_\textrm{BT} \geq 1$, HT > 1400 GeV. The background from QCD multijet and \ttbar production is derived from control regions in the data. Estimates for the \tttt signal and other backgrounds are shown using simulated samples.

The distribution of the BDT discriminants for the full 2016--2018 data set in the all-hadronic channel. The sample VR category shown is defined by $N_\textrm{RT}=1$, $N_\textrm{BT} \geq 1$, $H_T > 1400$ GeV. The background from QCD multijet and $\textrm{t}\bar{\textrm{t}}$ production is derived from control regions in the data. Estimates for the $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ signal and other backgrounds are shown using simulated samples. The hatched bands correspond to the estimated total uncertainty.

The distribution of the BDT discriminants for the full 2016--2018 data set in the all-hadronic channel. The SR category shown in this figure is defined by $N_\textrm{BT} \geq 2$, HT > 1100 GeV. The background from QCD multijet and \ttbar production is derived from control regions in the data. Estimates for the \tttt signal and other backgrounds are shown using simulated samples.

The distribution of the BDT discriminants for the full 2016--2018 data set in the all-hadronic channel. The sample VR category shown is defined by $N_\textrm{RT} \geq 2$, $H_T > 1100$ GeV. The background from QCD multijet and $\textrm{t}\bar{\textrm{t}}$ production is derived from control regions in the data. Estimates for the $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ signal and other backgrounds are shown using simulated samples. The hatched bands correspond to the estimated total uncertainty.

Expected and observed significance (in standard deviations) for $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ production from each final state and the combination with previous CMS results (EPJC 80 (2020) 75 and JHEP 11 (2019) 082). The same-sign dilepton and multilepton (SSDL\&ML) final state results are from EPJC 80 (2020) 75.

The distribution of the BDT discriminants for the full 2016--2018 data set in the all-hadronic channel. The SR category shown in this figure is defined by $N_\textrm{RT}=1$, $N_\textrm{BT} \geq 1$, $H_T > 1400$ GeV. The background from QCD multijet and $\textrm{t}\bar{\textrm{t}}$ production is derived from control regions in the data. Estimates for the $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ signal and other backgrounds are shown using simulated samples. The hatched bands correspond to the estimated total uncertainty after the fit.

Measured signal strength ($\mu = \sigma_{\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}}/ \sigma_{\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}} ^{\textrm{SM}}$), corresponding cross section (in fb), and the expected and observed significance (in standard deviations) for $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ production from all analysis channels. Uncertainties shown in this table are statistical and systematic uncertainties respectively. In this table, SSDL\&ML and OSDL (2016) represent results from EPJC 80 (2020) 75 and JHEP 11 (2019) 082 respectively.

The distribution of the BDT discriminants for the full 2016--2018 data set in the all-hadronic channel. The SR category shown in this figure is defined by $N_\textrm{RT} = 2$, $H_T > 1100$ GeV. The background from QCD multijet and $\textrm{t}\bar{\textrm{t}}$ production is derived from control regions in the data. Estimates for the $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ signal and other backgrounds are shown using simulated samples. The hatched bands correspond to the estimated total uncertainty after the fit.

Expected and observed significance (in standard deviations) for $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ production from each final state and the combination with previous CMS results (EPJC 80 (2020) 75 and JHEP 11 (2019) 082). The same-sign dilepton and multilepton (SSDL\&ML) final state results are from EPJC 80 (2020) 75.

Measured signal strength ($\mu = \sigma_{\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}}/ \sigma_{\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}} ^{\textrm{SM}}$), corresponding cross section (in fb), and the expected and observed significance (in standard deviations) for $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ production from all analysis channels. This table shows production from each analysis channel in this Letter, the combination of those channels, the results from previously published results, and the full combination of all CMS 2016-2018 results.

Observed and expected 95% CL upper limits on ${{(\sigma_{\text{H}}/\sigma_{\text{H}}^{\mathrm{SM}}) \times {{\mathcal{B}(\text{H} \to \text{inv})}}}}$ for the combined VBF, ttH, VH and ggH channels using all available CMS data for 2016--2018, Run1, 2015, and their combination, assuming a SM Higgs boson with a mass of 125 GeV.

Observed and expected 95% CL upper limits on ${{(\sigma_{\text{H}}/\sigma_{\text{H}}^{\mathrm{SM}}) \times {{\mathcal{B}(\text{H} \to \text{inv})}}}}$ for the VBF, ttH(lep), ttH(had), V(qq)H, Z(ll)H, Z(bb)H and ggH channels using all available CMS data, and their combination, assuming a SM Higgs boson with a mass of 125 GeV.

The 90% CL upper limits on the spin-independent DM-nucleon scattering cross section in Higgs-portal models, assuming a scalar, fermionic or vectorial DM candidate (with a UV-complete vector model), compared to direct detection experimets CRESST-III, DarkSide-50, LUX_ZEPLIN and PandaX-4T (2022 results).

Postfit event yields in each category and hadronic recoil bin of the ttH category SR for the different background and signal processes summed across 2016, 2017 and 2018 samples. The background and signal yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only. The Postfit signal contributions for a Higgs boson, produced in ttH, VH (split into WH and ZH), ggH and VBF, decaying to invisible particles with a branching fraction of $\mathcal{B}(\text{H} \to \text{inv}) = 1$, and the observed event yields are reported.

Postfit event yields in each category and hadronic recoil bin of the VH category SR for the different background and signal processes summed across 2016, 2017 and 2018 samples. The background and signal yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only. The Postfit signal contributions for a Higgs boson, produced in ttH, VH (split into WH and ZH), ggH and VBF, decaying to invisible particles with a branching fraction of $\mathcal{B}(\text{H} \to \text{inv}) = 1$, and the observed event yields are reported.

Postfit event yields in each category and $p\mathrm{_{T}^{miss}}$ bin of the ttH category e$+$jets CR for the different background processes summed across 2016, 2017 and 2018 samples. The background yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only.

Postfit event yields in each category and $p\mathrm{_{T}^{miss}}$ bin of the VH category e$+$jets CR for the different background processes summed across 2016, 2017 and 2018 samples. The background yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only.

Postfit event yields in each category and $p\mathrm{_{T}^{miss}}$ bin of the ttH category $\mu+$jets CR for the different background processes summed across 2016, 2017 and 2018 samples. The background yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only.

Postfit event yields in each category and $p\mathrm{_{T}^{miss}}$ bin of the VH category $\mu+$jets CR for the different background processes summed across 2016, 2017 and 2018 samples. The background yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only.

Postfit event yields in each category and $p\mathrm{_{T}^{miss}}$ bin of the ttH category ee$+$jets, $\mu\mu+$jets, and $\ell\ell+$jets CRs for the different background processes summed across 2016, 2017 and 2018 samples. The background yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only.

Postfit event yields in each category and $p\mathrm{_{T}^{miss}}$ bin of the VH category ee$+$jets, $\mu\mu+$jets, and $\gamma+$jets CRs for the different background processes summed across 2016, 2017 and 2018 samples. The background yields and the corresponding uncertainties are obtained after performing a combined fit across all of the CRs and SR or CR-only.

The expected (dotted black) and observed (solid black) upper limits at 95% CL on $\sigma_{T^\prime bq(T^\prime \rightarrow tH)}$ in the hadronic channel displayed as a function of $\rm{M}_{T^\prime}$. The green (yellow) band represents the 68% (95%) of the limit values expected under the background-only hypothesis. The theoretical cross sections for the singlet $T^\prime$ production with representative $\kappa_T$-values fixed at 0.1, 0.15, 0.2 and 0.25 (for $\Gamma/\rm{M}_{T^\prime} < 5\%$) are shown as solid red lines.

Selection of the hadronic signal region for the three BDTs and $\rm{M}_{T^\prime}$ windows.

Selection of the leptonic signal region for the three BDTs and $\rm{M}_{T^\prime}$ windows.

The expected yields of different processes in each signal window for events with a $T^{\prime}$ with mass in the range $\rm{M}_{T^\prime} \in [600, 1200]$ GeV, and the observed number of events in the signal region $m_{\gamma\gamma} \in [115, 135]$. Here, the yields for the $T^\prime$ are for $\kappa_{T}$ fixed at 0.2.

Signal efficiency for the hadronic channel (left) and the leptonic channel (right) in the signal regions as optimized for each of the three different $\rm{M}_{T^\prime}$ ranges, [600, 700], [700, 1000] and [1000, 1200] GeV, represented as red, green and blue curves respectively. The efficiency is defined as the ratio of the events after the final selection to the total expected events.

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

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.

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).

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).

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).

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).

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).

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.

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.

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 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 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 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 misidentification probability measured using the nominal $W$+jets MR is shown along with the systematic uncertainty, quantifying the degree of process dependence measured from alternative MRs. 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 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 is shown along with the systematic uncertainty, quantifying the degree of process dependence measured from alternative MRs. 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 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.

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.

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).

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 $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.

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.

Left pannel. The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 2.3%) is not graphically represented. The blue line represents the average for $p_T > 18$ GeV. For comparison, the LHCb measurement [10.1103/PhysRevLett.124.122002] is also shown. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (2.3%)

Left pannel. The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 2.3%) is not graphically represented. The blue line represents the average for $p_T > 18$ GeV. For comparison, the LHCb measurement [10.1103/PhysRevLett.124.122002] is also shown. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (2.3%)

Left pannel. The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 2.3%) is not graphically represented. The blue line represents the average for $p_T > 18$ GeV. For comparison, the LHCb measurement [10.1103/PhysRevLett.124.122002] is also shown. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (2.3%)

Right pannel. The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 2.3%) is not graphically represented. The blue line represents the average for $p_T > 18$ GeV. For comparison, the LHCb measurement [10.1103/PhysRevLett.124.122002] is also shown. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (2.3%)

Right pannel. The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 2.3%) is not graphically represented. The blue line represents the average for $p_T > 18$ GeV. For comparison, the LHCb measurement [10.1103/PhysRevLett.124.122002] is also shown. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (2.3%)

Right pannel. The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 2.3%) is not graphically represented. The blue line represents the average for $p_T > 18$ GeV. For comparison, the LHCb measurement [10.1103/PhysRevLett.124.122002] is also shown. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (2.3%)

Left panel.The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 4.9\%) is not graphically represented. The blue line represents the average of all the points. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (4.9%)

Left panel.The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 5.7%) is not graphically represented. The blue line represents the average of all the points. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (5.7%)

Left panel.The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 5.7%) is not graphically represented. The blue line represents the average of all the points. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (5.7%)

Right panel.The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 4.9\%) is not graphically represented. The blue line represents the average of all the points. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (4.9%)

Right panel.The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 5.7%) is not graphically represented. The blue line represents the average of all the points. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (5.7%)

Right panel.The vertical bars (boxes) represent the statistical (bin-to-bin systematic) uncertainties, while the horizontal bars give the bin widths. The global uncertainty (of 5.7%) is not graphically represented. The blue line represents the average of all the points. $ 12 < \mathrm{B} \, p_T < 70$ GeV and $ 0 < |y| < 2.4 $. Global uncertanties are not included in the table (5.7%)

Expected and observed upper limits on the AQGC operators $a^Z_C/\Lambda^2$, with no unitarization. The $y$ axis shows the limit on the ratio of the observed cross section to the cross section predicted for each anomalous coupling value ($\sigma_\mathrm{AQGC}$).

Limits on LEP-like dimension-6 anomalous quartic gauge coupling parameters, with and without unitarization via a clipping procedure.

Conversion of limits on $a^W_0$ to dimension-8 $f_{M,i}$ operators, using the assumption of vanishing $WWZ\gamma$ couplings to eliminate some parameters. When quoting limits on one of the operators, the other is fixed to zero. The results for $|f_{M,0}/\Lambda^{4}|$ and $|f_{M,4}/\Lambda^{4}|$ are shown with and without clipping of the signal model at 1.4 TeV, when the other parameter is fixed to the SM value of zero.

Conversion of limits on $a^W_0$ and $a^W_C$ to dimension-8 $f_{M,i}$ operators, using the assumption that all $f_{M,i}$ except one are equal to zero. The results are shown with and without clipping of the signal model at 1.4 TeV.

{Expected and observed limits in the two-dimensional plane of $a^W_0/\Lambda^2$ vs. $a^W_C/\Lambda^2$. The limits are described by analytical ellipses of equation $(x-x0)^2/a^2 + (y-y0)^2/b^2 = 1$ and rotated counter-clockwise by $\theta$ degrees, where $x$ and $y$ in the equation correspond to the $a_0^W$ and $a_C^W$ couplings, respectively.

{Expected and observed limits in the two-dimensional plane of $a^W_0/\Lambda^2$ vs. $a^W_C/\Lambda^2$ with unitarization imposed by clipping the signal model at 1.4 TeV. The limits are described by analytical ellipses of equation $(x-x0)^2/a^2 + (y-y0)^2/b^2 = 1$ and rotated counter-clockwise by $\theta$ degrees, where $x$ and $y$ in the equation correspond to the $a_0^W$ and $a_C^W$ couplings, respectively.

{Expected and observed limits in the two-dimensional plane of $a^Z_0/\Lambda^2$ vs. $a^Z_C/\Lambda^2$. The limits are described by analytical ellipses of equation $(x-x0)^2/a^2 + (y-y0)^2/b^2 = 1$ and rotated counter-clockwise by $\theta$ degrees, where $x$ and $y$ in the equation correspond to the $a_0^Z$ and $a_C^Z$ couplings, respectively.

Distributions of $\Delta\phi$ as obtained from simulation, requiring various $\textrm{t}$ tag multiplicities for the signal in two representative combinations of (gluino, neutralino) masses with large (2.2, 0.1)$\mathrm{TeV}$ and small (1.8, 1.3)$\mathrm{TeV}$ mass difference.

Results of fits to the $n_{\textrm{b}}$ multiplicity for control regions for the muon channel and with the requirements $3\leq n_{\textrm{jet}}\leq4$, $250<L_T<350~\mathrm{GeV}$, $500<H_T<750~\mathrm{GeV}$, $n_{\textrm{W}}\geq1$, $\Delta \phi<1$. The shaded area shows the fit uncertainty of the total background.

Results of fits to the $n_{\textrm{b}}$ multiplicity for control regions for the muon channel and with the requirements $3\leq n_{\textrm{jet}}\leq4$, $350<L_T<450~\mathrm{GeV}$, $H_T>1000~\mathrm{GeV}$, $n_{\textrm{W}}\geq0$, $\Delta \phi<1$. The shaded area shows the fit uncertainty of the total background.

Jet multiplicity distribution after the single-lepton baseline selection excluding the SRs for the multi-b analysis. The simulation is normalized to data with the SF mentioned in the plot.

Jet multiplicity distribution after the single-lepton baseline selection excluding the SRs for the the zero-b analysis (right). The simulation is normalized to data with the SF mentioned in the plot.

Jet multiplicity distribution in the dilepton CRs for the multi-b analysis. The simulation is normalized to data with the SF mentioned in the plot.

Jet multiplicity distribution in the dilepton CRs for the zero-b analysis. The simulation is normalized to data with the SF mentioned in the plot.

The double ratio of the single-lepton and dilepton ratio between data and simulation together with fit results and their uncertainties is shown for the multi-b analysis. The fits are performed for each data-taking year; 2018 is shown as an example.

The double ratio of the single-lepton and dilepton ratio between data and simulation together with fit results and their uncertainties is shown for the zero-b analysis. The fits are performed for each data-taking year; 2018 is shown as an example.

The prefit $L_{\mathrm{P}}$ distribution for selected electron candidates in the baseline QCD selection, with modified requirements of $n_{\text{jet}}\in[3,4]$ and $n_{\mathrm{b}}=0$.

The prefit $L_{\mathrm{P}}$ distribution for anti-selected electron candidates in the baseline QCD selection, with modified requirements of $n_{\text{jet}}\in[3,4]$ and $n_{\mathrm{b}}=0$.

Observed event yields in the MB SRs of the multi-b analysis compared to signal and background predictions. The relative fraction of the different SM EW background contributions determined in simulation is shown by the stacked, colored histograms, normalized so that their sum is equal to the background estimated using data control regions. The QCD background is predicted using the $L_p$ method. The signal is shown for two representative combinations of (gluino, neutralino) masses with large (2.2, 0.1) $\mathrm{TeV}$ and small (1.8, 1.3) $\mathrm{TeV}$ mass differences.

Observed event yields in the MB SRs of the zero-b analysis compared to signal and background predictions. The $\textrm{W}$+jets, $\textrm{t}\bar{\textrm{t}}$+jets, and QCD predictions are extracted from data control samples, while the other background contributions are estimated from simulation. The signal is shown for two representative combinations of (gluino, neutralino) masses with large (2.2, 0.1) $\mathrm{TeV}$ and small (1.8, 1.3) $\mathrm{TeV}$ mass differences.

Cross section limits at 95% CL for the T1tttt (left) model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Observed number of events in the MB SR bins of the multi-b analysis. All bins are defined with $\Delta\phi > 0.75$.

Observed number of events in the MB SR bins of the zero-b analysis.

The simplified likelihood can be constructed using this covariance matrix of nuisance parameters in the signal region (MB SR). The bin contents have been aggregated accross the three years to simplify the likelihood.

Selection efficiency for the multi-b signal for the SRs.

The simplified likelihood can be constructed using this covariance matrix of nuisance parameters in the signal region (MB SR). The bins were aggregated across HT to stabilize the simplified likelihood.

Selection efficiency for the zero-b signal for the SRs. The bins correspond to the aggregated bins that were used for the covariance matrix.

The covariance matrix containing the experimental uncertainties of the particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section as a function of the XCone-jet mass.

The covariance matrix containing the model uncertainties of the particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section as a function of the XCone-jet mass.

The normalized particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section in the fiducial region as a function of the XCone-jet mass.

Correlations between bins in the normalized particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section as a function of the XCone-jet mass.

The covariance matrix containing the statistical uncertainties of the normalized particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section as a function of the XCone-jet mass.

The covariance matrix containing the experimental uncertainties of the normalized particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section as a function of the XCone-jet mass.

The covariance matrix containing the model uncertainties of the normalized particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section as a function of the XCone-jet mass.

Relative experimental uncertainties of the particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section in the fiducial region as a function of the XCone-jet mass.

Relative model uncertainties of the particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section in the fiducial region as a function of the XCone-jet mass.

Relative experimental uncertainties of the normalized particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section in the fiducial region as a function of the XCone-jet mass.

Relative model uncertainties of the normalized particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section in the fiducial region as a function of the XCone-jet mass.

The measured cross section as a function of exclusive jet multiplicity, $N_{\text{jets}}$, when $50<p_T<100$ GeV

The measured cross section as a function of exclusive jet multiplicity, $N_{\text{jets}}$, when $p_T>100$ GeV

The measured cross section as a function of $\Delta\phi_{Z,jet1}$, when $p_T<10$ GeV

The measured cross section as a function of $\Delta\phi_{Z,jet1}$, when $10<p_T<30$ GeV

The measured cross section as a function of $\Delta\phi_{Z,jet1}$, when $30<p_T<50$ GeV

The measured cross section as a function of $\Delta\phi_{Z,jet1}$, when $50<p_T<100$ GeV

The measured cross section as a function of $\Delta\phi_{Z,jet1}$, when $p_T>100$ GeV

The measured cross section as a function of $\Delta\phi_{jet1,jet2}$, when $p_T<10$ GeV

The measured cross section as a function of $\Delta\phi_{jet1,jet2}$, when $10<p_T<30$ GeV

The measured cross section as a function of $\Delta\phi_{jet1,jet2}$, when $30<p_T<50$ GeV

The measured cross section as a function of $\Delta\phi_{jet1,jet2}$, when $50<p_T<100$ GeV

The measured cross section as a function of $\Delta\phi_{jet1,jet2}$, when $p_T>100$ GeV

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with 400 < $p_{T1}$ < 800 GeV and for an azimuthal separation between the two leading jets of $0 < \Delta\Phi_{1,2} < 150^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with 400 < $p_{T1}$ < 800 GeV and for an azimuthal separation between the two leading jets of $150 < \Delta\Phi_{1,2} < 170^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with 400 < $p_{T1}$ < 800 GeV and for an azimuthal separation between the two leading jets of $170 < \Delta\Phi_{1,2} < 180^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with $p_{T1}$ > 800 and for an azimuthal separation between the two leading jets of $0 < \Delta\Phi_{1,2} < 150^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with $p_{T1}$ > 800 and for an azimuthal separation between the two leading jets of $150 < \Delta\Phi_{1,2} < 170^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Jet multiplicity measured for a leading-pT jet ($p_{T1}$) with $p_{T1}$ > 800 and for an azimuthal separation between the two leading jets of $170 < \Delta\Phi_{1,2} < 180^{\circ}$. The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Measured transverse momentum of the leading $p_{T}$ jet ($p_{T1}$). The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Measured transverse momentum of the subleading $p_{T}$ jet ($p_{T2}$). The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.

Measured transverse momentum of the third leading $p_{T}$ jet ($p_{T3}$). The full breakdown of the uncertainties is displayed, with PU corresponding to Pileup, PREF to Trigger Prefering, PTHAT to the hard-scale (renormalization and factorization scales), MISS and FAKE to the inefficienties and background, LUMI to integrated luminosity. With JES, JER and stat. unc. following the notation in the paper.