Best fit values of $\sigma(gg\to H\to ZZ)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{ZZ}$ relative to their SM prediction from the combined analysis of the $\sqrt{s}$=7 and 8 TeV data. The results are shown for the combination of ATLAS and CMS, and also separately for each experiment, together with their total uncertainties and their breakdown into the four components described in the text. The expected uncertainties in the measurements are also shown.

Best fit values of $\sigma(gg\to H\to WW)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{WW}$ from the combined analysis of the $\sqrt{s}$=7 and 8 TeV data. The values involving cross sections are given for $\sqrt{s}$=8 TeV, assuming the SM values for $\sigma_i(7 \mathrm{TeV})/\sigma_i(8 \mathrm{TeV})$. The results are shown for the combination of ATLAS and CMS, and also separately for each experiment, together with their total uncertainties and their breakdown into the four components described in the text. The expected uncertainties in the measurements are also shown.

Best fit values of $\sigma(gg\to H\to WW)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{WW}$ relative to their SM prediction from the combined analysis of the $\sqrt{s}$=7 and 8 TeV data. The results are shown for the combination of ATLAS and CMS, and also separately for each experiment, together with their total uncertainties and their breakdown into the four components described in the text. The expected uncertainties in the measurements are also shown.

Best fit values of $\kappa_{gZ}= \kappa_{g}\cdot\kappa_{Z} / \kappa_{H}$ and of the ratios of coupling modifiers, as defined in the most generic parameterisation described in the context of the $\kappa$ framework, from the combined analysis of the $\sqrt{s}$=7 and 8 TeV data. The results are shown for the combination of ATLAS and CMS and also separately for each experiment, together with their total uncertainties and their breakdown into the four components described in the text. The uncertainties in $\lambda_{tg}$ and $\lambda_{WZ}$, for which a negative solution is allowed, are calculated around the overall best fit value. The expected total uncertainties in the measurements are also shown. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Measured global signal strength $\mu$ and its total uncertainty, together with the breakdown of the uncertainty into its four components as defined in the text. The results are shown for the combination of ATLAS and CMS, and separately for each experiment. The expected uncertainty, with its breakdown, is also shown.

Measured signal strengths $\mu$ and their total uncertainties for different Higgs boson production processes. The results are shown for the combination of ATLAS and CMS, and separately for each experiment, for the combined $\sqrt{s}$=7 and 8 TeV data. The expected uncertainties in the measurements are also displayed. These results are obtained assuming that the Higgs boson branching fractions are the same as in the SM.

Measured signal strengths $\mu$ and their total uncertainties for different Higgs boson decay channels. The results are shown for the combination of ATLAS and CMS, and separately for each experiment, for the combined $\sqrt{s}$=7 and 8 TeV data. The expected uncertainties in the measurements are also displayed. These results are obtained assuming that the Higgs boson production process cross sections at $\sqrt{s}$ = 7 and 8 TeV are the same as in the SM.

Results of the ten-parameter fit of $\mu_F^f = \mu_{gg\mathrm{F}+ttH}^f$ and $\mu_V^f = \mu_{\mathrm{VBF}+VH}^f$ for each of the five decay channels. The results are shown for the combination of ATLAS and CMS, together with their measured and expected uncertainties. The measured results are also shown separately for each experiment.

Results of the six-parameter fit of the global ratio $\mu_V/\mu_F = \mu_{\mathrm{VBF}+VH}/\mu_{gg\mathrm{F}+ttH}$ together with $\mu_F^f$ for each of the five decay channels. The results are shown for the combination of ATLAS and CMS, together with their measured and expected uncertainties. The measured results are also shown separately for each experiment.

Fit results allowing BSM loop couplings, assuming that $|\kappa_{V}| \le$ 1, where $\kappa_{V}$ denotes $\kappa_{Z}$ or $\kappa_{W}$, and that $\mathrm{B_{BSM}} \ge$ 0. The results for the combination of ATLAS and CMS are reported with their measured and expected uncertainties. Also shown are the results from each experiment. For the parameters with both signs allowed, the other 1$\sigma$ interval is shown on a second line. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Fit results allowing BSM loop couplings, assuming that there are no additional BSM contributions to the Higgs boson width, i.e. $\mathrm{B_{BSM}}=0$. The results for the combination of ATLAS and CMS are reported with their measured and expected uncertainties. Also shown are the results from each experiment. For the parameters with both signs allowed, the other 1$\sigma$ interval is shown on a second line. When a parameter is constrained and reaches a boundary, namely $\mathrm{B_{BSM}}$ = 0, the uncertainty is not indicated. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Fit results for the parameterisation assuming the absence of BSM particles in the loops ($\mathrm{B_{BSM}}$=0). The results with their measured and expected uncertainties are reported for the combination of ATLAS and CMS, together with the individual results from each experiment. For the parameters with both signs allowed, the other 1$\sigma$ CL interval is shown on a second line. When a parameter is constrained and reaches a boundary, namely $|\kappa_\mu|$ = 0, the uncertainty is not indicated. For those parameters with no sensitivity to the sign, only the absolute values are shown.

Summary of fit results for a parameterisation probing the ratios of coupling modifiers for up-type versus down-type fermions. The results for the combination of ATLAS and CMS are reported together with their measured and expected uncertainties. Also shown are the results from each experiment. The parameter $\kappa_{uu}$ is positive definite since $\kappa_H$ is always assumed to be positive. For the parameter $\lambda_{du}$, for which both signs are allowed, the other 1$\sigma$ CL interval is shown on a second line. Negative values for the parameter $\lambda_{Vu}$ are excluded by more than 4$\sigma$.

Summary of fit results for a parameterisation probing the ratios of coupling modifiers for leptons versus quarks. The results for the combination of ATLAS and CMS are reported together with their measured and expected uncertainties. Also shown are the results from each experiment. The parameter $\kappa_{qq}$ is positive definite since $\kappa_H$ is always assumed to be positive. For the parameter $\lambda_{lq}$, for which there is no sensitivity to the sign, only the absolute values are shown. Negative values for the parameter $\lambda_{Vq}$ are excluded by more than 4$\sigma$.

Correlation matrix obtained from the fit combining the ATLAS and CMS data using the generic parameterisation with 23 parameters, $\sigma_i\cdot\mathrm{B}^f$ for each specific channel $i\to H\to f$. Only 20 parameters are shown because the other three, corresponding to the $H\to ZZ$ decay channel for the $WH$, $ZH$, and $ttH$ production processes, are not measured with a meaningful precision.

Correlation matrix obtained from the fit combining the ATLAS and CMS data using the generic parameterisation with nine parameters, $\sigma(gg\to H\to ZZ)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{ZZ}$.

Correlation matrix obtained from the fit combining the ATLAS and CMS data using the generic parameterisation with seven parameters, $\kappa_{gZ}=\kappa_{g}\cdot\kappa_{Z} / \kappa_{H}$ and the ratios of coupling modifiers.

Expected correlation matrix obtained from the fit combining ATLAS and CMS pre-fit Asimov data sets using the generic parameterisation with 23 parameters, $\sigma_i\cdot\mathrm{B}^f$ for each specific channel $i\to H\to f$. Only 20 parameters are shown because the other three, corresponding to the $H\to ZZ$ decay channel for the $WH$, $ZH$, and $ttH$ production processes, are not measured with a meaningful precision.

Expected correlation matrix obtained from the fit combining ATLAS and CMS pre-fit Asimov data sets using the generic parameterisation with nine parameters, $\sigma(gg\to H\to ZZ)$, $\sigma_i/\sigma_{gg\mathrm{F}}$, and $\mathrm{B}^f/\mathrm{B}^{ZZ}$.

Expected correlation matrix obtained from the fit combining ATLAS and CMS pre-fit Asimov data sets using the generic parameterisation with seven parameters, $\kappa_{gZ}=\kappa_{g}\cdot\kappa_{Z} / \kappa_{H}$ and the ratios of coupling modifiers.

ATLAS+CMS combined best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{g}$,$\kappa_{\gamma}$) plane.

ATLAS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{g}$,$\kappa_{\gamma}$) plane.

CMS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{g}$,$\kappa_{\gamma}$) plane.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{F}$,$\kappa_{V}$) plane.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\gamma\gamma$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow ZZ$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow WW$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\tau\tau$ channel.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow bb$ channel.

ATLAS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{F}$,$\kappa_{V}$) plane.

CMS best-fit and 68% and 95% CL negative log-likelihood contours in the ($\kappa_{F}$,$\kappa_{V}$) plane.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\gamma\gamma$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow ZZ$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow WW$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow\tau\tau$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% CL negative log-likelihood contour in the ($\kappa_{F}$,$\kappa_{V}$) plane for the $H\rightarrow bb$ channel, assuming $\kappa_{F}>0$.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow\gamma\gamma$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow ZZ$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow WW$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow\tau\tau$ channel.

Best-fit and 68% and 95% CL negative log-likelihood contours in the ($\mu^{f}_{gg\mathrm{F}+ttH}$,$\mu^{f}_{\mathrm{VBF}+VH}$) plane for the $H\rightarrow bb$ channel.

Upper limits on the production cross section times branching fraction of the b* RH hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Upper limits on the production cross section times branching fraction of the b* VL hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Upper limits on the production cross section times branching fraction of the B+b hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

Upper limits on the production cross section times branching fraction of the B+t hypothesis at a 95% CL. Dashed colored lines show the expected limits from the l+jets and all-hadronic channels, where the latter start at resonance masses of 1.4 TeV. The observed and expected limits from the combination are shown as solid and dashed black lines, respectively. The green and yellow bands show the 68 and 95% confidence intervals on the combined expected limits.

The correlation matrix for the Higgs pT measurements, both for the unregularized and regularized fits. The last bin is inclusive.

The correlation matrix for the jet multiplicity measurements, both for the unregularized and regularized fits. The last bin is inclusive.

The correlation matrix for the leading jet pT measurements, both for the unregularized and regularized fits. The last bin is inclusive.

The fiducial integrated signal strength and cross section.

Summary of the sources of systematic and statistical uncertainties and their impact on the measurement of the ttH and tH signal rates, and the measured value of the unconstrained nuisance parameters. The quantity $\Delta r_{x}/r_{x}$ corresponds to the change in uncertainty when fixing the nuisance parameters associated with that uncertainty in the fit. Under the label "MC and sideband statistical uncertainty" are the uncertainties associated with the limited number of simulated MC events and the amount of data events in the application region of the MP method.

Number of events selected in the 2$\ell$ss+0$\tau_{\mathrm{h}}$, 3$\ell$+0$\tau_{\mathrm{h}}$, and 2$\ell$ss+1$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 2$\ell$ss+0$\tau_{\mathrm{h}}$, 3$\ell$+0$\tau_{\mathrm{h}}$, and 2$\ell$ss+1$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 1$\ell$+1$\tau_{\mathrm{h}}$, 0$\ell$+2$\tau_{\mathrm{h}}$, 2$\ell$os+1$\tau_{\mathrm{h}}$, and 1$\ell$+2$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 1$\ell$+1$\tau_{\mathrm{h}}$, 0$\ell$+2$\tau_{\mathrm{h}}$, 2$\ell$os+1$\tau_{\mathrm{h}}$, and 1$\ell$+2$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 4$\ell$+0$\tau_{\mathrm{h}}$, 3$\ell$+1$\tau_{\mathrm{h}}$, and 2$\ell$+2$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Number of events selected in the 4$\ell$+0$\tau_{\mathrm{h}}$, 3$\ell$+1$\tau_{\mathrm{h}}$, and 2$\ell$+2$\tau_{\mathrm{h}}$ analysis channels compared to the event yields expected from the $ tH$ and $ H$ signals and from background processes. The expected event yields are computed for the values of nuisance parameters and of the POI obtained from the ML fit. The best fit values of the POI amount to $\hat{r_{ttH}} = 0.92$ and $\hat{r_{tH}} = 5.7$. Quoted uncertainties represent the sum of statistical and systematic components.

Probability for ttH signal events to pass the event selection in each of the analysis categories employed.

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.

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

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.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the dijet-1 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the dijet-2 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the dijet-3 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the untagged 1 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the untagged 2 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the untagged 3 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Fits to the $m_{\ell^+\ell^-\gamma}$ data distribution in the untagged 4 categories. In the upper panel, the red solid line shows the result of a signal-plus-background fit to the given category. The red dashed line shows the background component of the fit. The green and yellow bands represent the 68 and 95% CL uncertainties in the fit. Also plotted is the expected SM signal, scaled by a factor of 10. In the lower panel, the data minus the background component of the fit is shown.

Observed signal strength ($\mu$) for a SM Higgs boson at 125.38 GeV. The labels untagged combined,dijet combined and combined represent the results obtained from simultaneous fits of the untagged categories, dijet categories, and full set of categories, respectively. The black solid line shows $\mu=1$, and the red dashed line shows the best fit value $\hat{\mu}=2.4$ of all categories combined.

Upper limit (95% CL) on the signal strength ($\mu$) relative to the SM prediction, as a function of the assumed value of the Higgs boson mass used in the fit.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow d\bar{d}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 7 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 7 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow \tau^{+} \tau^{-}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 15 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 40 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The 95% CL observed and expected limits on the branching fraction B(H $\rightarrow$ SS) for 55 GeV mass and $ S \rightarrow b\bar{b}$ decay mode.

The cluster efficiency in bins of hadronic and EM energy in region A. Region A is defined as 391 cm $< r <$ 695.5 cm and 400 cm $< |z| <$ 671 cm. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The cluster efficiency in bins of hadronic and EM energy in region A. Region A is defined as 391 cm $< r <$ 695.5 cm and 400 cm $< |z| <$ 671 cm. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The cluster efficiency in bins of hadronic and EM energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The cluster efficiency in bins of hadronic and EM energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bins correspond to LLPs that decayed leptonically with 0 hadronic energy. The cluster efficiency includes all cluster-level selections described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 35% and 20% for region A and B, respectively.

The efficiency of $N_{station} > 1$ requirement in bins of hadronic energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bin corresponds to LLPs that decayed leptonically with 0 hadronic energy. The efficiency is calculated with respect to clusters that pass all cluster-level cuts described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 10%.

The efficiency of $N_{station} > 1$ requirement in bins of hadronic energy in region B. Region B is defined as 671 cm $< |z| <$ 1100 cm, $r <$ 695.5 cm, and $|\eta| <$ 2. The cluster efficiency is estimated with LLPs decaying to $\tau^{+} \tau^{-}$. The sample contains equal fractions of events with LLP mass of 7, 15, 40, and 55 GeV and LLP lifetime of 0.1, 1, 10, and 100m. The first hadronic energy bin corresponds to LLPs that decayed leptonically with 0 hadronic energy. The efficiency is calculated with respect to clusters that pass all cluster-level cuts described in the paper, except for the jet veto, time cut, and $\Delta\phi$ cut. The full simulation signal yield prediction for samples with various LLP mass between 7 - 55 GeV, lifetime between 0.1 - 100 m, and decay mode to $d\bar{d}$ and $\tau^{+} \tau^{-}$ can be reproduced using this parameterization to within 10%.

The geometric signal acceptance as a function of $c\tau$. The acceptance is shown for four different LLP mass hypotheses: 7, 15, 40, and 55 GeV. The acceptance is defined by requiring at least 1 LLP to decay in the region defined as 400 cm $< |z| <$ 1100 cm, $r < $ 695.5 cm, and $|\eta| <$ 2.4.