Post-fit distributon of BDT discriminants for $\rho_{tc}=1.0$ with $m_A$ = 350 GeV interfered with H.($m_A - m_H$ = 50 GeV)

Observed and expected 95% CL upper limit on the signal strength as a function of $m_A$ for the g2HDM signal model using different coupling assumptions: $\rho_{tu}$ = 0.1, 0.4, 1.0.

Observed and expected 95% CL upper limit on the signal strength as a function of $m_A$ for the g2HDM signal model using different coupling assumptions: $\rho_{tc}$ = 0.1, 0.4, 1.0.

Observed and expected 95% CL upper limit on the signal strength as a function of $m_A$ for the g2HDM signal model using different coupling assumptions: $\rho_{tu}$ = 0.1, 0.4, 1.0 with A–H interference assuming $m_A$ − $m_H$ = 50 GeV.

Observed and expected 95% CL upper limit on the signal strength as a function of $m_A$ for the g2HDM signal model using different coupling assumptions: $\rho_{tc}$ = 0.1, 0.4, 1.0 with A–H interference assuming $m_A$ − $m_H$ = 50 GeV.

Observed 95% CL upper limit on the signal strength as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model.

Expected exclusion as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model

Observed exclusion as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model

Observed 95% CL upper limit on the signal strength as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model.

Expected exclusion as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model

Observed exclusion as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model

Observed 95% CL upper limit on the signal strength as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV.

Expected exclusion as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV

Observed exclusion as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV

Observed 95% CL upper limit on the signal strength as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV.

Expected exclusion as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV

Observed exclusion as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV

Product of detector efficiency and analysis acceptance for signal samples with various $m_a$ values for the muon channel.

Post-fit distributions of $M_{\ell v jj}$ in the $R2B_A$ subregion for electrons. All process yields and nuisance parameters are set to the values obtained from the background plus signal fit. The signal considered for the fit corresponds to the purely right-handed production of a W' with $m_{W'}$ of 3.6 TeV and a relative width of 1$\%$ of the $m_{W'}$, and is represented by the solid red line. The lower panels show the data minus the expected number of events, normalized to the statistical uncertainty of the data. The orange band represents the systematic uncertainties, also normalized to the statistical uncertainty of the data.

Post-fit distributions of $M_{\ell v jj}$ in the $RW'_A$ subregion for muons. All process yields and nuisance parameters are set to the values obtained from the background plus signal fit. The signal considered for the fit corresponds to the purely right-handed production of a W' with $m_{W'}$ of 3.6 TeV and a relative width of 1$\%$ of the $m_{W'}$, and is represented by the solid red line. The lower panels show the data minus the expected number of events, normalized to the statistical uncertainty of the data. The orange band represents the systematic uncertainties, also normalized to the statistical uncertainty of the data.

Post-fit distributions of $M_{\ell v jj}$ in the $RW'_A$ subregion for electrons. All process yields and nuisance parameters are set to the values obtained from the background plus signal fit. The signal considered for the fit corresponds to the purely right-handed production of a W' with $m_{W'}$ of 3.6 TeV and a relative width of 1$\%$ of the $m_{W'}$, and is represented by the solid red line. The lower panels show the data minus the expected number of events, normalized to the statistical uncertainty of the data. The orange band represents the systematic uncertainties, also normalized to the statistical uncertainty of the data.

Post-fit distributions of $M_{\ell v jj}$ in the $RT_A$ subregion for muons. All process yields and nuisance parameters are set to the values obtained from the background plus signal fit. The signal considered for the fit corresponds to the purely right-handed production of a W' with $m_{W'}$ of 3.6 TeV and a relative width of 1$\%$ of the $m_{W'}$, and is represented by the solid red line. The lower panels show the data minus the expected number of events, normalized to the statistical uncertainty of the data. The orange band represents the systematic uncertainties, also normalized to the statistical uncertainty of the data.

Post-fit distributions of $M_{\ell v jj}$ in the $RT_A$ subregion for electrons. All process yields and nuisance parameters are set to the values obtained from the background plus signal fit. The signal considered for the fit corresponds to the purely right-handed production of a W' with $m_{W'}$ of 3.6 TeV and a relative width of 1$\%$ of the $m_{W'}$, and is represented by the solid red line. The lower panels show the data minus the expected number of events, normalized to the statistical uncertainty of the data. The orange band represents the systematic uncertainties, also normalized to the statistical uncertainty of the data.

Observed and expected 95% CL upper limits on the product of the production cross section for a left-handed production of a tb quark pair in the $s$-channel, mediated by a W or W' bosons and including interference terms, given as functions of $m_{W'}$ for a relative width of 1%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a left-handed production of a tb quark pair in the $s$-channel, mediated by a W or W' bosons and including interference terms, given as functions of $m_{W'}$ for a relative width of 10%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a left-handed production of a tb quark pair in the $s$-channel, mediated by a W or W' bosons and including interference terms, given as functions of $m_{W'}$ for a relative width of 20%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a left-handed production of a tb quark pair in the $s$-channel, mediated by a W or W' bosons and including interference terms, given as functions of $m_{W'}$ for a relative width of 30%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of $m_{W'}$ for a relative width of 1%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of $m_{W'}$ for a relative width of 10%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of $m_{W'}$ for a relative width of 20%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of $m_{W'}$ for a relative width of 30%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed 95% CL upper limit on the production cross section for a left-handed W' boson in the tb final state, as functions of $m_{W'}$ and relative width $\Gamma/m_{W'}$. Numbers in red, written diagonally, represent values of the excluded cross sections that are lower than the theoretical ones for the analyzed model.

Observed 95% CL upper limit on the production cross section for a right-handed W' boson in the tb final state, as functions of $m_{W'}$ and relative width $\Gamma/m_{W'}$. Numbers in red, written diagonally, represent values of the excluded cross sections that are lower than the theoretical ones for the analyzed model.

Observed 95% CL upper limit on the production cross section for a generalized left-right coupling of the W' boson to a t and a b quark for a mass of the W' boson of 2 TeV.

Observed 95% CL upper limit on the production cross section for a generalized left-right coupling of the W' boson to a t and a b quark for a mass of the W' boson of 2.8 TeV.

Observed 95% CL upper limit on the production cross section for a generalized left-right coupling of the W' boson to a t and a b quark for a mass of the W' boson of 3.6 TeV.

Observed 95% CL upper limit on the production cross section for a generalized left-right coupling of the W' boson to a t and a b quark for a mass of the W' boson of 4.4 TeV.

Observed 95% CL upper limit on the production cross section for a generalized left-right coupling of the W' boson to a t and a b quark for a mass of the W' boson of 5.2 TeV.

Observed 95% CL upper limit on the production cross section for a generalized left-right coupling of the W' boson to a t and a b quark for a mass of the W' boson of 6 TeV.

Expected 95% CL lower limit on $m_{W'}$ for a generalized left-right coupling of the W' boson to a t and a b quark.

Observed 95% CL lower limit on $m_{W'}$ for a generalized left-right coupling of the W' boson to a t and a b quark.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of the $m_{W'}$ for a relative width of 10%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid curves show the theoretical expectation at LO in the case that the couplings, and thus the partial widths, are varied together with the total width. In this interpretation the branching fraction of the W' --> tb decays is the same for each value of the width of the W' boson.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of the $m_{W'}$ for a relative width of 20%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid curves show the theoretical expectation at LO in the case that the couplings, and thus the partial widths, are varied together with the total width. In this interpretation the branching fraction of the W' --> tb decays is the same for each value of the width of the W' boson.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of the $m_{W'}$ for a relative width of 30%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid curves show the theoretical expectation at LO in the case that the couplings, and thus the partial widths, are varied together with the total width. In this interpretation the branching fraction of the W' --> tb decays is the same for each value of the width of the W' boson.

Observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance. The acceptance $\mathcal{A}$ is defined as $\mathcal{A} = N$(events with $m_{X}^{GEN} > 85\% m_{X}^{input}$) / $N$(events generated in the full phase space defined by the CMS default generator settings). Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Mass exclusion ranges of the benchmark signal scenarios are depicted with hatched areas inside the black contours.

Observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance. The acceptance $\mathcal{A}$ is defined as $\mathcal{A} = N$(events with $m_{X}^{GEN} > 85\% m_{X}^{input}$) / $N$(events generated in the full phase space defined by the CMS default generator settings). Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Mass exclusion ranges of the benchmark signal scenarios are depicted with hatched areas inside the black contours.

Efficiencies of the selection requirements on the benchmark signal processes: $Z_{R} \to ggg$ with nominal width ($\Gamma_{X}/m_{X}\sim 3\%$). A value of -1 means that the corresponding efficiency is not calculated for this year.

Efficiencies of the selection requirements on the benchmark signal processes: $Z_{R} \to ggg$ with narrow width ($\Gamma_{X}/m_{X}\sim 0.01\%$). A value of -1 means that the corresponding efficiency is not calculated for this year.

Efficiencies of the selection requirements on the benchmark signal processes: $G_{KK} \to {\varphi}(gg)g$, where $\varphi$ is the radion. A value of -1 means that the corresponding efficiency is not calculated for this year.

Efficiencies of the selection requirements on the benchmark signal processes: $q^{*} \to V(qq)q$, where $V$ is a beyond-the-SM vector boson. A value of -1 means that the corresponding efficiency is not calculated for this year.

Acceptance of the signal selection requirement $m_{X}^{\text{GEN}}/m_{X}^{\text{input}} > 85\%$ on the benchmark signal process $Z_{R} \to ggg$. The acceptance is defined as $\mathcal{A} = N$(events with $m_{X}^{\text{GEN}}/m_{X}^{\text{input}} > 85\%$)/$N$(events generated in the full phase space defined by the CMS default generator settings).

Acceptance of the signal selection requirement $m_{X}^{\text{GEN}}/m_{X}^{\text{input}} > 85\%$ on the benchmark signal process $G_{KK} \to {\varphi}(gg)g$.

Acceptance of the signal selection requirement $m_{X}^{\text{GEN}}/m_{X}^{\text{input}} > 85\%$ on the benchmark signal process $q^{*} \to V(qq)q$.

Observed local significance for a 3-body decay $ggg$ resonance, shown for resonances with nominal width (blue solid line) and narrow width (red dashed line).

Observed local significance for a cascade decay $ggg$ resonance.

Observed local significance for a cascade decay $qqq$ resonance.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.2$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.3$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.4$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.5$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.6$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.7$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.8$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance $m_{Y} / m_{X} = 0.2$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.3$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.4$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.5$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.6$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.7$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.8$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal

Cut flow table of the selection requirments for the $Z_{R}$ model scenarios. Values shown are absolute efficiencies, e.g., values shown in the column of 'Efficiency ($\Delta_{R}^{max} < 3.0$)' represent the cumulative efficiencies achieved through all selection requirements applied up to and including the current selection criterion.

Cut flow table of the selection requirments for the $G_{KK}$ model scenarios. Values shown are absolute efficiencies, e.g., values shown in the column of 'Efficiency ($\Delta_{R}^{max} < 3.0$)' represent the cumulative efficiencies achieved through all selection requirements applied up to and including the current selection criterion.

Cut flow table of the selection requirments for the $q^{*}$ model scenarios. Values shown are absolute efficiencies, e.g., values shown in the column of 'Efficiency ($\Delta_{R}^{max} < 3.0$)' represent the cumulative efficiencies achieved through all selection requirements applied up to and including the current selection criterion.

Predicted production cross section of the benchmark signal process $pp \to Z_{R} \to ggg$.

Predicted production cross section of the benchmark signal process $pp \to G_{KK} \to \varphi(gg)g$.

Predicted production cross section of the benchmark signal process $pp \to q^{*} \to V(qq)q$.

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

The number of events in data and simulations after the fit to data in the Top$\gamma$ CR, SS$WW\gamma$ CR, and signal region.

The distributions of the invariant mass of the dilepton-photon system after the signal region selection before the fit to data.The shaded band represents the Monte Carlo statistical unc

Upper limits on the cross section and derived limits in terms of Yukawa coupling at 95% CL for Hγ production initiated by light quarks.

The distributions of the invariant mass of the transverse mass of the WW system after the signal region selection before the fit to data.The shaded band represents the Monte Carlo statistical unc

The distributions of the invariant mass of the dilepton-photon system after the signal region selection before the fit to data.The shaded band represents the Monte Carlo statistical unc

The unrolled two-dimensional mTWW-mll$\gamma$ distributions with category 0 jet after the fit to data. The data are compared to the sum of the signal and expected background. The black points with error bars represent the data and their statistical uncertainties, whereas the hatched bands represent the total uncertainties in the predictions.

The distributions of the invariant mass of the transverse mass of the WW system after the signal region selection before the fit to data.The shaded band represents the Monte Carlo statistical unc

The unrolled two-dimensional mTWW-mll$\gamma$ distributions with category 1more jets after the fit to data. The data are compared to the sum of the signal and expected background. The black points with error bars represent the data and their statistical uncertainties, whereas the hatched bands represent the total uncertainties in the predictions.

The unrolled two-dimensional mTWW-mll$\gamma$ distributions with category 0 jet after the fit to data. The data are compared to the sum of the signal and expected background. The black points with error bars represent the data and their statistical uncertainties, whereas the hatched bands represent the total uncertainties in the predictions.

The unrolled two-dimensional mTWW-mll$\gamma$ distributions with category 1more jets after the fit to data. The data are compared to the sum of the signal and expected background. The black points with error bars represent the data and their statistical uncertainties, whereas the hatched bands represent the total uncertainties in the predictions.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 350 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 400 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 450 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 500 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 550 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 600 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 650 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 700 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 750 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 800 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 850 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 900 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 950 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 1000 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.