The spin-density matrix elements obtained from the low Q**2 BPC data set in two Q**2 intervals.

The spin-density matrix elements obtained from the low Q**2 BPC data set in three T intervals.

The spin-density matrix elements obtained from the low Q**2 BPC data set in two region of the 2PI mass from the RHO0 decay.

Correlation matrix for the 15 combinations of spin-density matrix elements obtained from the BPC data.

Correlation matrix for the 15 combinations of spin-density matrix elements obtained from the DIS data.

Pseusdorapidity densities for the INEL>0 class of data as a function of pseudorapidity for centre-of mass energy 0.9 TeV. Note that this data is not in the paper but has been approved by the collaboration.

Pseusdorapidity densities for the INEL>0 class of data as a function of pseudorapidity for centre-of mass energy 2.36 TeV. Note that this data is not in the paper but has been approved by the collaboration.

Pseusdorapidity densities for the INEL>0 class of data as a function of pseudorapidity for centre-of mass energy 7 TeV. Note that this data is not in the paper but has been approved by the collaboration.

ASYM for ETA mesons measured as a function of PT for ABS(XF) > 0.2. Uncertainties listed are those due to the statistics, the PT uncorrelated uncertainties due to extracting ASYM, and the correlated relative luminosity uncertainty.

The measured production spectra for xi hadrons as a function of pT.

The measured production spectra for phi hadrons as a function of pT.

The ratio of cross sections as a function of PT for LAMBDA/K0S production.

Distribution of the missing transverse energy in CRTop. The Z+HF and $t\bar{t}$ scale factors, described in the text, have been applied to the simulated samples. The distribution labeled "Signal" is for the model with ($m_a$, $m_{\tilde{\chi}_{1}^{0}}$, $m_{\tilde{\chi}_{2}^{0}}$) = (45 GeV, 10 GeV, 80 GeV).

Distribution of the dijet invariant mass in the signal region, shown together with the parameterized background model (labelled "Bkg Model"). For reference, the MC prediction for the SM background is also shown (labelled "SM MC"). The Z+HF and $t\bar{t}$ scale factors, described in the text, have been applied to the simulated samples. The signal region is defined to have dijet invariant mass > 20 GeV. The distribution labeled "Signal" is for the model with ($m_a$, $m_{\tilde{\chi}_{1}^{0}}$, $m_{\tilde{\chi}_{2}^{0}}$) = (45 GeV, 10 GeV, 80 GeV).

Upper limits at 95% CL on the $pp \rightarrow ZH$ cross section times the branching ratio for $Z \rightarrow \ell^{+}\ell^{-}$ (where $\ell = e, \mu \;\mathrm{or}\; \tau$) and $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 65$ GeV in the NMSSM scenario described in the text. All branching ratios in the Higgs boson decay chain after the decay $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0}$ are set to 100%.

Upper limits at 95% CL on the $pp \rightarrow ZH$ cross section times the branching ratio for $Z \rightarrow \ell^{+}\ell^{-}$ (where $\ell = e, \mu \;\mathrm{or}\; \tau$) and $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 80$ GeV in the NMSSM scenario described in the text. All branching ratios in the Higgs boson decay chain after the decay $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0}$ are set to 100%.

Upper limits at 95% CL on the $pp \rightarrow ZH$ cross section times the branching ratio for $Z \rightarrow \ell^{+}\ell^{-}$ (where $\ell = e, \mu\; \mathrm{or}\; \tau$) and $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 95$ GeV in the NMSSM scenario described in the text. All branching ratios in the Higgs boson decay chain after the decay $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0}$ are set to 100%.

Upper limits at 95% CL on the $pp \rightarrow ZH$ cross section times the branching ratio for $Z \rightarrow \ell^{+}\ell^{-}$ (where $\ell = e, \mu\; \mathrm{or}\; \tau$) and $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 110$ GeV in the NMSSM scenario described in the text. All branching ratios in the Higgs boson decay chain after the decay $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0}$ are set to 100%.

Upper limits at 95% CL on the $pp \rightarrow ZH$ cross section times the branching ratio for $Z \rightarrow \ell^{+}\ell^{-}$ (where $\ell = e, \mu\; \mathrm{or}\; \tau$) and $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 20$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 80$ GeV in the NMSSM scenario described in the text. All branching ratios in the Higgs boson decay chain after the decay $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0}$ are set to 100%.

Upper limits at 95% CL on the $pp \rightarrow ZH$ cross section times the branching ratio for $Z \rightarrow \ell^{+}\ell^{-}$ (where $\ell = e, \mu \;\mathrm{or}\; \tau$) and $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 30$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 80$ GeV in the NMSSM scenario described in the text. All branching ratios in the Higgs boson decay chain after the decay $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0}$ are set to 100%.

Upper limits at 95% CL on the branching ratio $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 65$ GeV in the NMSSM scenario described in the text. The SM $ZH$ cross section is assumed.

Upper limits at 95% CL on the branching ratio $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 80$ GeV in the NMSSM scenario described in the text. The SM $ZH$ cross section is assumed.

Upper limits at 95% CL on the branching ratio $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 95$ GeV in the NMSSM scenario described in the text. The SM $ZH$ cross section is assumed.

Upper limits at 95% CL on the branching ratio $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 10$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 110$ GeV in the NMSSM scenario described in the text. The SM $ZH$ cross section is assumed.

Upper limits at 95% CL on the branching ratio $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 20$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 80$ GeV in the NMSSM scenario described in the text. The SM $ZH$ cross section is assumed.

Upper limits at 95% CL on the branching ratio $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ as a function of $m_{a}$ for $m_{\tilde{\chi}_{1}^{0}} = 30$ GeV and $m_{\tilde{\chi}_{2}^{0}} = 80$ GeV in the NMSSM scenario described in the text. The SM $ZH$ cross section is assumed.

Unweighted and weighted number of events after each stage of selection for the NMSSM scenario with $pp \rightarrow ZH$, $Z \rightarrow \ell^{+}\ell^{-}$, $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ where $(m_{a}, m_{\tilde{\chi}_{1}^{0}}, m_{\tilde{\chi}_{2}^{0}}) = (45,10,80)$ GeV and the Z boson decaying to $e^{+}e^{-}, \;\mu^{+}\mu^{-} \;\mathrm{or}\; \tau^{+}\tau^{-}$. All branching ratios in the Higgs boson decay chain are set to 100%. The weighted number of events corresponds to an integrated luminosity of $139 \;\mathrm{fb}^{-1}$. The "skimming selection" required either a single electron or muon with $p_{T} > 100\;$ GeV or a pair of electrons or muons each with $p_{T} > 20\;$ GeV. The preselection requirements include the trigger, absence of a bad jet or bad muon, and two leptons, where a lepton is either an electron or a muon.

Acceptance and efficiency of this analysis for the signal models considered in this paper. The signal is an NMSSM scenario with $pp \rightarrow ZH$, $Z \rightarrow \ell^{+}\ell^{-}$, $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ where the values of $(m_{a}, m_{\tilde{\chi}_{1}^{0}}, m_{\tilde{\chi}_{2}^{0}})$ are varied. The Z boson decays to $e^{+}e^{-}, \mu^{+}\mu^{-} \;\mathrm{or}\; \tau^{+}\tau^{-}$. All branching ratios in the Higgs boson decay chain are set to 100%. The product of acceptance times efficiency is defined by the fraction of simulated events that pass all the selection criteria of this analysis. The acceptance is defined by the fraction of simulated events that pass all the selection criteria as applied to Monte Carlo truth-level quantities. The efficiency is then defined as the acceptance times efficiency, divided by the acceptance.

Acceptance and efficiency of this analysis for the signal models considered in this paper. The signal is an NMSSM scenario with $pp \rightarrow ZH$, $Z \rightarrow \ell^{+}\ell^{-}$, $H \rightarrow \tilde{\chi}_{2}^{0}\tilde{\chi}_{1}^{0} \rightarrow a \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0} \rightarrow b\bar{b} \tilde{\chi}_{1}^{0}\tilde{\chi}_{1}^{0}$ where the values of $(m_{a}, m_{\tilde{\chi}_{1}^{0}}, m_{\tilde{\chi}_{2}^{0}})$ are varied. The Z boson decays to $e^{+}e^{-}, \mu^{+}\mu^{-} \;\mathrm{or}\; \tau^{+}\tau^{-}$. All branching ratios in the Higgs boson decay chain are set to 100%. The product of acceptance times efficiency is defined by the fraction of simulated events that pass all the selection criteria of this analysis. The acceptance is defined by the fraction of simulated events that pass all the selection criteria as applied to Monte Carlo truth-level quantities. The efficiency is then defined as the acceptance times efficiency, divided by the acceptance.

Post-fit $m_\text{T}(\gamma, E_\text{T}^\text{miss})$ distribution in the inclusive signal region for the dark-photon search with the 125 GeV mass $\mathcal{B}(H \to \gamma \gamma_\text{d})$ signal normalization set to zero. A $H \to \gamma \gamma_\text{d}$ decay signal is shown for two different mass hypotheses, 125 GeV and 500 GeV, and scaled to a $\mathcal{B}(H \to \gamma \gamma_\text{d})$ of 2% and 1%, respectively. Events with $m_\text{T}(\gamma, E_\text{T}^\text{miss})$ larger than the rightmost bin boundary are added to that bin.

The 95% CL upper limit on the Higgs boson production cross-section times branching ratio to $\gamma \gamma_\text{d}$ is shown for different VBF-produced scalar-mediator-mass hypotheses in the NWA. The theoretically predicted cross-section of a Higgs boson produced via VBF and with the $\mathcal{B}(H \to \gamma \gamma_\text{d}) =$ 5% is superimposed on the $\pm 1\sigma$ and $\pm 2\sigma$ NNLO QCD + NLO EW uncertainty band of the expected production cross-section limit.

Post-fit $m_\text{jj}$ distribution in the inclusive signal region. The Higgs boson invisible decay signal is scaled to a $\mathcal{B}_\text{inv}$ of 37%. Events with $m_\text{jj}$ larger than the rightmost bin boundary are added to that bin.

Post-fit $m_\text{jj}$ distribution in the one-lepton control region $W_{\ell \nu}^\gamma$ CR. Events with $m_\text{jj}$ larger than the rightmost bin boundary are added to that bin.

Post-fit $m_\text{T}$ distribution in the one lepton control region. Events with $m_\text{T}$ larger than the rightmost bin boundary are added to that bin.

Post-fit photon centrality distribution in the zero lepton signal plus control region with the $\mathcal{B}_\text{inv}$ signal normalization set to zero in the fit.

Post-fit photon $E_\text{T}$ distribution in the zero lepton signal region with the $\mathcal{B}_\text{inv}$ signal normalization set to zero in the fit.

Post-fit photon centrality distribution in the zero lepton signal plus control region resulting from the fit to the $m_\text{jj}$ distribution for EW $Z \gamma + \text{jets}$. The post-fit uncertainties include statistical, experimental, and theory contributions.

Post-fit photon $E_\text{T}$ distribution in the zero lepton signal region resulting from the fit to the $m_\text{jj}$ distribution for EW $Z \gamma + \text{jets}$. The post-fit uncertainties include statistical, experimental, and theory contributions.

Post-fit DNN output score distribution in the one lepton control region.

Yields for the EW $Z \gamma + \text{jets}$ process are shown after each selection along with relative and absolute signal acceptance efficiencies.

Yields for the 125 GeV Higgs boson with $\mathcal{B}_\text{inv.} =$ 1 signal produced by the vector boson fusion process in association with a final state photon are shown after each selection along with relative and absolute signal acceptance efficiencies.

Yields for the 125 GeV Higgs boson with $\mathcal{B}(H \to \gamma \gamma_\text{d}) =$ 1 signal produced by the vector boson fusion process are shown after each selection along with relative and absolute signal acceptance efficiencies.

Expected +- 2 sigma 95% CL exclusion limit for the Zprime-2HDM model.

Observed 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected +- 1 sigma 95% CL exclusion limit for the 2HDM+a model ggF production.

Expected +- 2 sigma 95% CL exclusion limit for the 2HDM+a model ggF production.

Observed 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected +- 1 sigma 95% CL exclusion limit for the 2HDM+a model bbA production.

Expected +- 2 sigma 95% CL exclusion limit for the 2HDM+a model bbA production.

Observed 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected +- 1 sigma 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected +- 2 sigma 95% CL exclusion limit for the Zprime-2HDM model with the benchmark used in arXiv:1707.01302.

Expected and observed upper limits at 95% CL on cross-section for Zprime-2HDM model.

Expected and observed upper limits at 95% CL on cross-section for ggF producton in the 2HDM+a model.

Expected and observed upper limits at 95% CL on cross-section for bbA producton in the 2HDM+a model.

Model-independent upper limits on the visible cross-section $σ_{vis, $h(\bar{b})+DM} ≡ σ_{h+DM} \times B(h \to b\bar{b}) \times \mathcal{A} \times \epsilon$ in the different signal regions.

Theory cross-section for Zprime-2HDM model.

Theory cross-section for bbA production in the 2HDM+a model.

Theory cross-section for ggF production in the 2HDM+a model.

Distribution of Higgs boson candidate mass in 2b region with MET=150-200 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=200-350 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=350-500 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET=500-750 GeV.

Distribution of Higgs boson candidate mass in 2b region with MET > 750 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=150-200 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=200-350 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET=350-500 GeV.

Distribution of Higgs boson candidate mass in 3b region with MET > 500 GeV.

Yields in 1-lepton control region.

Yields in 2-lepton control region.

MET distribution in 2b region of the 0-lepton channel.

MET distribution in 3b region of the 0-lepton channel.

Expected signal yields after certain selection cuts in 2b region with MET=150-200 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=200-350 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=350-500 GeV.

Expected signal yields after certain selection cuts in 2b region with MET=500-750 GeV.

Expected signal yields after certain selection cuts in 2b region with MET > 750 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=150-200 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=200-350 GeV.

Expected signal yields after certain selection cuts in 3b region with MET=350-500 GeV.

Expected signal yields after certain selection cuts in 3b region with MET > 500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 2b region with MET > 750 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for bbA production in the 2HDM+a model - 3b region with MET>500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 2b region with MET > 750 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the 2HDM+a model - 3b region with MET > 500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET=500-750 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 2b region with MET > 750 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=150-200 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=200-350 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET=350-500 GeV.

Acceptance times efficiency for ggF production in the Zprime-2HDM model - 3b region with MET > 500 GeV.