Correlation matrix for the differential branching fraction extraction between different $q^2$ bins in the simultaneous fit.

The product of acceptance and efficiency ($A\epsilon$) of the $\mathcal{B}(\mathrm{B}^+\to\mathrm{K}^+\mu^+\mu^-)$ channel, as a function of the muon pair $q^2$, as measured in simulated signal events, after all the corrections applied. In the figure, regions corresponding to resonances are displayed with red markers.

Same as Table 5, but for the bins containing the J$\psi$ and $\psi$(2S) resonances.

Integrated branching fraction $\mathcal{B}(\mathrm{B}^\pm \to \mathrm{K}^\pm\mu^+\mu^-)$ in the low-$q^2$ region.

Ratio of branching fractions $\mathcal{B}(\text{B}^\pm \to \text{K}^\pm\mu^+\mu^-)$ to $\mathcal{B}(\text{B}^\pm \to \text{K}^\pm\text{e}^+\text{e}^-)$, determined from the measured double ratio $R$(K) of these decays to the respective branching fractions of the decay chains $\text{B}^\pm\to\text{J/}\psi\text{K}^\pm$ with $\text{J/}\psi\to\mu^+\mu^-$ and $\text{e}^+\text{e}^-$.

Negative log likelihood function from the fit profiled as a function of the inverse ratio $R(\mathrm{K})^{-1}$.

Signal yields for two charge scenarios considered in the analysis, as well as their associated uncertainties.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the early 2016 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 6%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text. Comparing with respect to the binomial fit starting at $N_{\text{hits}}^{\text{low dE/dx}} = 2$, and not $N_{\text{hits}}^{\text{low dE/dx}} = 1$, is needed to account for the fact that early 2016 data is more strongly affected low dE/dx by instrumental effects that widen the N hits distribution.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the late 2016 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 78%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the 2017 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 65%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the 2018 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 9%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Exclusion region (hatched) at 95% CL in the FCP charge-mass plane for the considered signal. The expected exclusion is shown with the associated 1 and 2 standard deviations $\sigma$ bands. Signal points at charges 0.9, 0.8, 2/3, 0.5, and 1/3 e are connected by straight lines to guide the eye. This is a conservative interpolation. Previous exclusions from CMS [Phys. Rev. D 87 (2013) 092008), JHEP 07 (2013) 122] as well as OPAL [Phys. Lett. B 572 (2003) 8] are given for comparison.

Example description

Example description

Example description

Unfolded E2C distributions in data compared to MC predictions.

Unfolded E2C distributions in data compared to MC predictions.

Unfolded E2C distributions in data compared to MC predictions.

Unfolded E2C distributions in data compared to MC predictions.

Unfolded E2C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

The fitted slopes of the E3C/E2C data distributions as a function of jet pt are used to illustrate the dependency of alphas on jet pt.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

The chi2 scan result using eq.3 for different alphas.

The correlation matrix is composed by 10 jet pt region, each region is represented by a block in the plot. Inside each block, there are 22 xL bins same as the E2C, E3C and E3C/E2C distributions. Therefore, the x and y bins of the correlation matrix is given by, binNumber = pT_index * 22 + xL_index.

The correlation matrix is composed by 10 jet pt region, each region is represented by a block in the plot. Inside each block, there are 22 xL bins same as the E2C, E3C and E3C/E2C distributions. Therefore, the x and y bins of the correlation matrix is given by, binNumber = pT_index * 22 + xL_index.

The correlation matrix is composed by 10 jet pt region, each region is represented by a block in the plot. Inside each block, there are 22 xL bins same as the E2C, E3C and E3C/E2C distributions. Therefore, the x and y bins of the correlation matrix is given by, binNumber = pT_index * 22 + xL_index.

The energy weight of E2C as defined in eq.1 in the paper is used in the unfolding, the binning is listed in this table.

The energy weight of E3C as defined in eq.1 in the paper is used in the unfolding, the binning is listed in this table.

Transverse dilepton-$p_{T}^{miss}$ mass pre-fit distribution for selected events in SR2 of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton mass pre-fit distribution for selected events in SR3 of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton $p_T$ pre-fit distribution for selected events in SR3 of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dijet $p_T$ pre-fit distribution for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Lepton-dijet $p_T$ pre-fit distribution for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Maximum $p_T$ of the visible particles divided by the missing $p_T$, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between jets, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and missing $p_T$, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and the dijet system, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta R$ distance between the leptons, for selected events in the non-prompt validation region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the di-lepton system and the missing $p_T$, for selected events in the non-prompt validation region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and the di-jet system, for selected events in the non-prompt validation region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton-dijet system and the missing $p_T$, for selected events in the non-prompt validation region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Leading lepton $p_T$ pre-fit distribution for selected events in the top control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Trailing lepton $p_T$ pre-fit distribution for selected events in the top control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton mass pre-fit distribution for selected events in the Drell-Yan control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Missing $p_T$ pre-fit distribution for selected events in the Drell-Yan control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Missing $p_T$ pre-fit distribution for selected events in the $WW$ control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton $p_T$ pre-fit distribution for selected events in the $WW$ control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Minimum $p_T$ of the visible particles divided by the missing $p_T$, for selected events in the top control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between jets, for selected events in the top control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and missing $p_T$, for selected events in the $W$+jets control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and the dijet system, for selected events in the $W$+jets control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Unrolled $m_{ll}-m_{T}^{l min, p_{T}^{miss}}$ post-fit distributions in the di-leptonic channel in SR1 for the full data set. The histogram bins are spaced uniformly. Each group of five bins (from left to right) corresponds to $m_{T}^{l min, p_{T}^{miss}}$ distribution in a $m_{ll}$ region, placed in ascending order. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$.

Unrolled $m_{ll}-m_{T}^{l min, p_{T}^{miss}}$ post-fit distributions in the di-leptonic channel in SR2 for the full data set. The histogram bins are spaced uniformly. Each group of five bins (from left to right) corresponds to $m_{T}^{l min, p_{T}^{miss}}$ distribution in a $m_{ll}$ region, placed in ascending order. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$.

Unrolled $m_{ll}-m_{T}^{l min, p_{T}^{miss}}$ post-fit distributions in the di-leptonic channel in SR3 for the full data set. The histogram bins are spaced uniformly. Each group of five bins (from left to right) corresponds to $m_{T}^{l min, p_{T}^{miss}}$ distribution in a $m_{ll}$ region, placed in ascending order. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$.

Post-fit BDT distributions in the semi-leptonic channel for the full data set in the top control region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Post-fit BDT distributions in the semi-leptonic channel for the full data set in the $W$+jets control region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Post-fit BDT distributions in the semi-leptonic channel in the 2016 signal region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Post-fit BDT distributions in the semi-leptonic channel in the 2017+2018 signal region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Observed and expected 95% CL upper limits on the nonresonant HH production via vector boson fusion cross section obtained for both single-lepton and dilepton channels, and from their combination

Observed and expected 95% CL upper limits on the nonresonant HH production cross section for two different benchmark scenarios “JHEP04(2016)01” and “JHEP03(2020)91”

Observed and expected 95% CL upper limits on the nonresonant HH production cross section as a function of the Higgs boson self-coupling strength modifier κλ All Higgs boson couplings other than λ are assumed to have the values predicted by the SM

Observed and expected 95% CL upper limits on the nonresonant HH production via VBF cross section as a function of the effective coupling κ2V. The ggF contribution in this case is set to the SM expectation. All other Higgs boson couplings are assumed to have the values predicted by the SM.

Observed and expected 95% CL upper limits on the nonresonant HH production cross section as a function of the effective coupling c2. All other Higgs boson couplings are assumed to have the values predicted in the SM

Upper limits at 95% CL on $\sigma$B(pp→X→Y(bb)H) for combination as a function of m$_Y$.

Expected upper limits at 95% CL on the X → HH cross section as a function of mX from the combination of the three channels bbγγ, bbττ, and bbbb for an integrated luminosity of 3000 fb−1 with different assumptions on the systematic uncertainties.

Expected upper limits on the X → YH cross section times B(Y → bb) for mX ≤ 1000 GeV at 95 % CL as a function of m$_X$ and m$_Y$ from the combination of all considered channels projected to an integrated luminosity of 3000 fb$^{-1}$. Systematic uncertainties according to the ”S2” scenario are assumed.

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.

Invariant mass distribution of the diphoton pairs for the elastic selection region with events satisfying a < 0.005 (signal contribution is magnified by a factor 5000).

Predicted number of events having an elastic diphoton pair in association with a pair of protons observed within the range where the proton detectors have a radiation inefficiency less than 10%. The yields where the two-photon and two-proton systems mass and rapidity are matching at 2 and 3$\sigma$ are also quoted. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

Predicted number of events having an elastic diphoton pair in association with a pair of protons observed within the range where the proton detectors have a radiation inefficiency less than 10%. The yields where the two-photon and two-proton systems mass and rapidity are matching at 2 and 3$\sigma$ are also quoted. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

95% CL expected and observed frequentist upper limit for the LbL cross section within the $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$ search region.

95% CL expected and observed frequentist upper limit for the LbL cross section within the $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$ search region.

95% expected and observed one-dimensional limits on $\zeta_1$ and $\zeta_2$ anomalous LbyL production parameters, when the other parameter is set to zero. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

95% expected and observed one-dimensional limits on $\zeta_1$ and $\zeta_2$ anomalous LbyL production parameters, when the other parameter is set to zero. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

95% expected and observed limits on $\zeta_1$ and $\zeta_2$ anomalous LbyL production parameters. The parametric elliptic form is assumed: $a_0\zeta_1^2+a_1\zeta_1\zeta_2+a_2\zeta_2^2=1$. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.

95% expected and observed limits on $\zeta_1$ and $\zeta_2$ anomalous LbyL production parameters. The parametric elliptic form is assumed: $a_0\zeta_1^2+a_1\zeta_1\zeta_2+a_2\zeta_2^2=1$. This corresponds to a search region of $m_{\gamma\gamma} > 350$ GeV, $0.070 < \xi^+ < 0.111$, and $0.070 < \xi^- < 0.138$.