The distributions of the $\mathrm{D}^{*}-\mathrm{D}^0$ mass difference $\Delta m$ for the normalization channel in data.

The distributions of the dimuon invariant mass $m_{\mu\mu}$ for the signal channel in data with the requirement $0.145<\Delta m<0.146$ GeV.

The distributions of the $\mathrm{D}^{*}-\mathrm{D}^0$ mass difference $\Delta m$ for the signal channel in data with the requirement $1.84<m_{\mu\mu}<1.89$ GeV.

The post-fit event yields for the signal, the combinatorial background, the $\mathrm{D}^0\to\pi^+\pi^-$ background, and the $\mathrm{D}^0\to\pi^-\mu^+\nu$ background. The observed numbers of events are given in the Data column. The subrange is in one dimension with a full range in the other dimension.

Likelihood profiles of the $\kappa_\lambda$ parameter.

Likelihood profiles of the $\kappa_{2V}$ parameter.

The 95% CL upper limits on the $HH$ cross-section for different values of $\kappa_\lambda$. The expected limits assume no $HH$ production. The red line shows the theory prediction for the combined ggF $HH$ and VBF $HH$ cross-section as a function of $\kappa_\lambda$ where all parameters and couplings are set to their SM values except $\kappa_\lambda$. The bands surrounding the red cross-section line indicate the theoretical uncertainty on the predicted cross-section. The $\kappa_\lambda$ values in the interval $[-6.7, 15.3]$ cannot be excluded by this analysis.

The 95% CL upper limits on the VBF $HH$ cross section for different values of $\kappa_{2V}$. The expected limits assume no VBF $HH$ production. SM ggF $HH$ production is considered as background. The red line shows the theory prediction for the VBF $HH$ cross-section as a function of $\kappa_{2V}$ where all parameters and couplings are set to their SM values except $\kappa_{2V}$. The uncertainty band is smaller than the width of the plotted line. The $\kappa_{2V}$ values in the interval $[-0.3, 2.4]$ cannot be excluded by this analysis.

Expected CLs values in the High mA SR, mA vs ma signal grid.

Observed CLs values in the High mA SR, mA vs ma signal grid.

Expected CLs values in the Low mA SR, mA vs ma signal grid.

Observed CLs values in the High mA SR, mA vs ma signal grid.

CLs+b values in the Low mA SR, mA vs tanB signal grid.

CLs+b values in the High mA SR, mA vs ma signal grid.

CLs+b values in the Low mA SR, mA vs ma signal grid.

Cut flow of the 2HDM+a signal points, gluon–gluon fusion production, Low mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

Cut flow of the 2HDM+a signal points, gluon–gluon fusion production, High mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

Cut flow of the 2HDM+a signal points, bb annihilation production, Low mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

Cut flow of the 2HDM+a signal points, bb annihilation production, High mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

A comparison of the observed and expected yields in the four bins of the Low mA SR.

A comparison of the observed and expected yields in the two bins of the High mA SR.

Expected exclusion contours at 95% CL as a function of mA and ma.

Observed exclusion contours at 95% CL as a function of mA and ma.

Expected +- 1sigma exclusion contours at 95% CL as a function of mA and ma.

Expected exclusion contours at 95% CL as a function of mA and ma.

Observed exclusion contours at 95% CL as a function of mA and ma.

Expected +- 1sigma exclusion contours at 95% CL as a function of mA and ma.

Acceptance, High mA SR, mA vs tanB grid, 400-750 GeV, bb prod

Acceptance, High mA SR, mA vs tanB grid, >750 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, 100-250 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, 250-400 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, 400-550 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, >550 GeV, bb prod

Acceptance, High mA SR, mA vs ma grid, 400-750 GeV, bb prod

Acceptance, High mA SR, mA vs ma grid, >750 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, 100-250 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, 250-400 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, 400-550 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, >550 GeV, bb prod

Acceptance, High mA SR, mA vs tanB grid, 400-750 GeV, ggF prod

Acceptance, High mA SR, mA vs tanB grid, >750 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, 100-250 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, 250-400 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, 400-550 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, >550 GeV, ggF prod

Acceptance, High mA SR, mA vs ma grid, 400-750 GeV, ggF prod

Acceptance, High mA SR, mA vs ma grid, >750 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, 100-250 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, 250-400 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, 400-550 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, >550 GeV, ggF prod

Efficiency, High mA SR, mA vs tanB grid, 400-750 GeV, bb prod

Efficiency, High mA SR, mA vs tanB grid, >750 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, 100-250 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, 250-400 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, 400-550 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, >550 GeV, bb prod

Efficiency, High mA SR, mA vs ma grid, 400-750 GeV, bb prod

Efficiency, High mA SR, mA vs ma grid, >750 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, 100-250 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, 250-400 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, 400-550 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, >550 GeV, bb prod

Efficiency, High mA SR, mA vs tanB grid, 400-750 GeV, ggF prod

Efficiency, High mA SR, mA vs tanB grid, >750 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, 100-250 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, 250-400 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, 400-550 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, >550 GeV, ggF prod

Efficiency, High mA SR, mA vs ma grid, 400-750 GeV, ggF prod

Efficiency, High mA SR, mA vs ma grid, >750 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, 100-250 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, 250-400 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, 400-550 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, >550 GeV, ggF prod

Acceptance * reconstruction efficiency for the P P --> bbA --> Zh --> llbb signals in the 2-lepton channel.

Acceptance * reconstruction efficiency for the P P --> Wprime --> Zh --> lvbb/cc signals in the 0-lepton channel.

Acceptance * reconstruction efficiency for the P P --> Wprime --> Zh --> lvbb/cc signals in the 1-lepton channel.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the resolved 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the resolved 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the merged 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the merged 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the resolved 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the resolved 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the merged 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the merged 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 3+ b-tag signal region. The background prediction is shown after a background-only maximum-likelihood bbA fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the resolved 3+ b-tag signal region. The background prediction is shown after a background-only maximum-likelihood bbA fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag signal region with additional b-tagged track jets not associated with the large-R jet. The background prediction is shown after a background-only maximum-likelihood bbA fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the merged 1+2 b-tag signal region with additional b-tagged track jets not associated with the large-R jet. The background prediction is shown after a background-only maximum-likelihood bbA fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the resolved 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the resolved 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the merged 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 1-lepton channel in the merged 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{Vh}$ for the 2-lepton channel in the resolved top control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 1 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag sideband control region. The background prediction is shown after a background-only maximum-likelihood Z' fit to the data. In the plot, the last bin contains the overflow.

Upper limits on Zprime to Z h production cross section times branching fraction in pb.

Upper limits on Wprime to W h production cross section times branching fraction in pb.

Upper limits on ggA to Z h production cross section times branching fraction in pb.

Upper limits on bbA to Z h production cross section times branching fraction in pb.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 220 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 260 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 300 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 340 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 380 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 400 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 420 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 440 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 460 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 500 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 600 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 700 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 800 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 900 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 1000 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 1200 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 1400 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 1600 GeV.

Expected and observed two-dimensional likelihood scans of the b-associated production cross section times branching fraction vs the gluon-fusion production cross section times branching fraction at $m_{A}$ = 2000 GeV.

Acceptance * reconstruction efficiency for the P P --> A --> Zh --> vvbb signal in the 0-lepton channel.

Acceptance * reconstruction efficiency for the P P --> A --> Zh --> llbb signal in the 2-lepton channel.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the resolved 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 1 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Event distributions of $m_{T,Vh}$ for the 0-lepton channel in the merged 2 b-tag signal region. The background prediction is shown after a background-only maximum-likelihood W' fit to the data. In the plot, the last bin contains the overflow.

Distributions of expected upper limits at 95% confidence level on the cross section of P P --> A --> Zh as a function of bbA fraction an signal mass.

Distributions of observed upper limits at 95% confidence level on the cross section of P P --> A --> Zh as a function of bbA fraction an signal mass.

Excluded signal cross section at $95\%$ CL for Hut coupling. Each jet category and their combination are shown separately.

Excluded signal cross section at $95\%$ CL for Hut coupling. Each jet category and their combination are shown separately.

Observed upper limits on the branching fractions $\mathcal{B}(\mathrm{t}\to\mathrm{Hu})$ and $\mathcal{B}(\mathrm{t}\to\mathrm{Hc})$ at $95\%$ CL.

Expected upper limits on the branching fractions $\mathcal{B}(\mathrm{t}\to\mathrm{Hu})$ and $\mathcal{B}(\mathrm{t}\to\mathrm{Hc})$ at $95\%$ CL.

Expected +1 standard deviation upper limits on the branching fractions $\mathcal{B}(\mathrm{t}\to\mathrm{Hu})$ and $\mathcal{B}(\mathrm{t}\to\mathrm{Hc})$ at $95\%$ CL.

Expected -1 standard deviation upper limits on the branching fractions $\mathcal{B}(\mathrm{t}\to\mathrm{Hu})$ and $\mathcal{B}(\mathrm{t}\to\mathrm{Hc})$ at $95\%$ CL.

Expected +2 standard deviation upper limits on the branching fractions $\mathcal{B}(\mathrm{t}\to\mathrm{Hu})$ and $\mathcal{B}(\mathrm{t}\to\mathrm{Hc})$ at $95\%$ CL.

Expected -2 standard deviation upper limits on the branching fractions $\mathcal{B}(\mathrm{t}\to\mathrm{Hu})$ and $\mathcal{B}(\mathrm{t}\to\mathrm{Hc})$ at $95\%$ CL.

Observed upper limits on the anomalous couplings $\kappa_{\mathrm{Hut}}$ and $\kappa_{\mathrm{Hut}}$ at $95\%$ CL.

Expected upper limits on the anomalous couplings $\kappa_{\mathrm{Hut}}$ and $\kappa_{\mathrm{Hut}}$ at $95\%$ CL.

Expected +1 standard deviation upper limits on the anomalous couplings $\kappa_{\mathrm{Hut}}$ and $\kappa_{\mathrm{Hut}}$ at $95\%$ CL.

Expected -1 standard deviation upper limits on the anomalous couplings $\kappa_{\mathrm{Hut}}$ and $\kappa_{\mathrm{Hut}}$ at $95\%$ CL.

Expected +2 standard deviation upper limits on the anomalous couplings $\kappa_{\mathrm{Hut}}$ and $\kappa_{\mathrm{Hut}}$ at $95\%$ CL.

Expected -2 standard deviation upper limits on the anomalous couplings $\kappa_{\mathrm{Hut}}$ and $\kappa_{\mathrm{Hut}}$ at $95\%$ CL.

The +1 $\sigma$ band of the observed exclusion contor for the $Z^{\prime}_{B}$ model in the $m_{\chi}$-$m_{Z^{\prime}_{B}}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the $Z^{\prime}_{B}$ model in the $m_{\chi}$-$m_{Z^{\prime}_{B}}$ plane.

A comparison of the inferred limits to the constraints from direct detection experiments on the spin-independent DM--nucleon cross section in the context of the $Z'_B$ simplified model with vector couplings. Limits are shown at 90% CL.

The observed exclusion contor for the $Z^{\prime}$-2HDM model in the $m_{A}$-$m_{Z^{\prime}}$ plane.

The expected exclusion contor for the $Z^{\prime}$-2HDM model in the $m_{A}$-$m_{Z^{\prime}}$ plane.

The +1 $\sigma$ band of the observed exclusion contor for the $m_{A}$-$Z^{\prime}$-2HDM model in the $m_{Z^{\prime}}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the $m_{A}$-$Z^{\prime}$-2HDM model in the $m_{Z^{\prime}}$ plane.

The observed exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The expected exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The +1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $m_{A}$-$m_{a}$ plane.

The observed exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The expected exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The +1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The -1 $\sigma$ band of the observed exclusion contor for the 2HDM-a model in the $tan\beta$-$m_{a}$ plane.

The exclusion limits at 95% CL for the 2HDM+a model as a function of $\sin \theta$ for $m_{A,H^{\pm},H}$= 600GeV, $m_a$ = 200GeV, $\tan \beta$ = 1.0$.

The exclusion limits at 95% CL for the 2HDM+a model as a function of $\sin \theta$ for $m_{A,H^{\pm},H}$= 1000GeV, $m_a$ = 350GeV, $\tan \beta$ = 1.0$.

Breakdown of the dominant systematic uncertainties.

Event yields in the range of 120 $<m_{\gamma\gamma}<$ 130 GeV for data, signal models, the SM Higgs boson background and non-resonant background in each analysis category, for an integrated luminosity of $139$fb$^{-1}$.

Detailed background contributions from the SM Higgs boson and continuum background for each cut

Detailed contributions from the signals for each cut.

Acceptance times efficiency for several signals in each category.

Upper limits on A to Z h production cross section x branching fraction in pb (gluon fusion production)

Upper limits on A to Z h production cross section x branching fraction in pb ( production with associated b-quarks)

Acceptance * Reconstruction efficiency for pp-> Zprime

Acceptance * Reconstruction efficiency for pp-> Zprime

Acceptance * Reconstruction efficiency for pp-> Wprime

Acceptance * Reconstruction efficiency for pp-> A (gluon fusion)

Acceptance * Reconstruction efficiency for pp-> A (gluon fusion)

Acceptance * Reconstruction efficiency for pp-> A (b-quark associated)

Acceptance * Reconstruction efficiency for pp-> A (b-quark associated)

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 0% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 10%.

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 0% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 20%.

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 0% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 30%.

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 0% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 40%.

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 0% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 50%.

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 0% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 60%.

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 0% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 70%.

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 0% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 80%.

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 0% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 90%.

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 1% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 0%.

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 2% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 0%.

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 3% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 0%.

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 4% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 0%.

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 5% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 0%.

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 6% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 0%.

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 7% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 0%.

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 8% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 0%.

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 9% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 0%.

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 10% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 0%.

Upper limits at the 95% CL on the product of the production cross-section for pp->A and the branching fractions for A->Zh and h->bb evaluated by combining the 0-lepton and 2-lepton channels. The signal is smeared by a Breit-Wigner function with A boson width of 11% , assuming a combination of the gluon--gluon fusion and b-quark associated production modes with a bbA fraction of 0%.

Event distributions of mT,Vh for the 0-lepton channel in the resolved 1-btag category. The background prediction is shown after a background-only maximum-likelihood fit to the data.

Event distributions of mT,Vh for the 0-lepton channel in the resolved 2-btag category. The background prediction is shown after a background-only maximum-likelihood fit to the data.

Event distributions of m,Vh for the 1-lepton channel in the resolved 1-btag category. The background prediction is shown after a background-only maximum-likelihood fit to the data.

Event distributions of m,Vh for the 1-lepton channel in the resolved 2-btag category. The background prediction is shown after a background-only maximum-likelihood fit to the data.

Event distributions of mT,Vh for the 2-lepton channel in the resolved 1-btag category. The background prediction is shown after a background-only maximum-likelihood fit to the data.

Event distributions of m,Vh for the 2-lepton channel in the resolved 2-btag category. The background prediction is shown after a background-only maximum-likelihood fit to the data.

Event distributions of mT,Vh for the 0-lepton channel in the boosted 1-btag category. The background prediction is shown after a background-only maximum-likelihood fit to the data.

Event distributions of mT,Vh for the 0-lepton channel in the boosted 2-btag category. The background prediction is shown after a background-only maximum-likelihood fit to the data.

Event distributions of m,Vh for the 1-lepton channel in the boosted 1-btag category. The background prediction is shown after a background-only maximum-likelihood fit to the data.

Event distributions of m,Vh for the 1-lepton channel in the boosted 2-btag category. The background prediction is shown after a background-only maximum-likelihood fit to the data.

Event distributions of m,Vh for the 2-lepton channel in the boosted 1-btag category. The background prediction is shown after a background-only maximum-likelihood fit to the data.

Event distributions of m,Vh for the 2-lepton channel in the boosted 2-btag category. The background prediction is shown after a background-only maximum-likelihood fit to the data.