Single-differential production cross-sections for nonprompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainties are statistical, the second are correlated systematic uncertainties shared between intervals, and the last are uncorrelated systematic uncertainties.

Single-differential production cross-sections for prompt $J/\psi$ mesons as a function of $y$. The first uncertainties are statistical, the second are correlated systematic uncertainties shared between intervals, and the last are uncorrelated systematic uncertainties.

Single-differential production cross-sections for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainties are statistical, the second are correlated systematic uncertainties shared between intervals, and the last are uncorrelated systematic uncertainties.

Fraction of nonprompt $J/\psi$ mesons (in \%) in ($p_\text{T},y$) intervals. The first uncertainty is statistical and the second is systematic.

Nuclear modification factor $R_{p\text{Pb}}$ for prompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Nuclear modification factor $R_{p\text{Pb}}$ for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for prompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 8 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $p_\text{T}$. The first uncertainty is statistical and the second is systematic.

Cross-section ratios between 13 TeV and 5 TeV measurements for nonprompt $J/\psi$ mesons as a function of $y$. The first uncertainty is statistical and the second is systematic.

Relative changes of cross-sections (in \%), for a polarisation of $\lambda_{\theta}=-0.2$ rather than zero, in ($p_\text{T},y$) intervals.

Relative changes of cross-sections (in \%), for a polarisation of $\lambda_{\theta}=-1$ rather than zero, in ($p_\text{T},y$) intervals.

Relative changes of cross-sections (in \%), for a polarisation of $\lambda_{\theta}=+1$ rather than zero, in ($p_\text{T},y$) intervals.

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.

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 event yields in the signal and control regions obtained from a scan over $\epsilon_{VV}$ exploiting both shape and rate information. The quoted uncertainties include the theoretical and experimental systematic sources and those due to sample statistics. The fit constrains the total expected yield to the observed yield. The diboson background is split into $W W$ and non-$W W$ contributions.

Best-fit values and their uncertainties as obtained from the shape-only and shape-plus-rate likelihood fits to the Asimov dataset and to ATLAS data. Results of both shape-only and shape+rate fits for $a_L$ and $a_T$ are shown. Results of fits to one parameter with the other one fixed or profiled are presented.

Best-fit values and their uncertainties as obtained from the shape-only and shape-plus-rate likelihood fits to the Asimov dataset and to ATLAS data. Results of both shape-only and shape+rate fits for $\epsilon_{VV}$ and $\kappa_{VV}$ are shown. Results of fits to one parameter with the other one fixed or profiled are presented.

The contributions of the leading individual systematic uncertainties together with the data statistical uncertainties, in the one dimensional fit for the pseudo-observables $\kappa_{VV}$ (a) and $\epsilon_{VV}$ (b) for electroweak-boson polarisation in the VBF $H\to WW$ channel. Both shape and rate informations are exploited in the fit. The theoretical and experimental uncertainties are subdivided further into categories.

The contributions of the leading individual systematic uncertainties together with the data statistical uncertainties, in the one dimensional fit for the pseudo-observables $\kappa_{VV}$ (a) and $\epsilon_{VV}$ (b) for electroweak-boson polarisation in the VBF $H\to WW$ channel. Both shape and rate informations are exploited in the fit. The theoretical and experimental uncertainties are subdivided further into categories.

Post-fit distribution of the BDT response observable presented in the four $|\Delta \eta jj|$ categories of the ggF +2 jets signal region, with signal and background yields fixed from the fit for $\mu^{\text{ggF+2jets}}$. The distributions of the ggF + 2 jets and VBF processes are overlaid with their respective contributions multiplied by 50.

The weighted $\Delta \Phi_{jj}$ post-fit distribution in the ggF +2 jets signal region, with signal and background yields fixed from the fit to $\tan \alpha$ using shape and rate information.

Expected and observed likelihood curves for scans over $\tan \alpha$ where only the shape is taken into account in the fit, $\mu_{VBF}$ is fixed.

Expected and observed likelihood curves for scans over $\tan \alpha$ where both shape and normalisation are taken into account in the fit, $\mu_{VBF}$ is fixed.

68% and 95% CL two-dimensional likelihood contours of the CP-even and CP-odd coupling parameters $K_{gg} \cos(\alpha)$ and $K_{gg} \sin(\alpha)$. The minima are represented by black stars, while the SM value is shown as a red star.

The weighted $\Delta \Phi jj$ distribution in the VBF signal region, with signal and background yields fixed from the fit for $\epsilon_{VV}$ using shape and rate information.

Likelihood scans over the transversally polarised couplings. The fit is using shape-only information. All relevant experimental and modelling systematic uncertainties are considered in the fit.

Likelihood scans over the transversally polarised couplings. The fit is using shape + rate information. All relevant experimental and modelling systematic uncertainties are considered in the fit.

Likelihood scans over the longitudinally polarised couplings. The fit is using shape + rate information. All relevant experimental and modelling systematic uncertainties are considered in the fit.

Likelihood scans over $\kappa_{VV}$ with the $\epsilon_{VV}$ profiled. The fit is performed using both shape and rate information. All relevant experimental and theoretical systematic uncertainties are considered in the fit.

Likelihood scans over $\epsilon_{VV}$ with the $\kappa_{VV}$ profiled. The fit is performed using both shape and rate information. All relevant experimental and theoretical systematic uncertainties are considered in the fit.

The contributions of the leading individual systematic uncertainties together with the data statistical uncertainties, in the one dimensional fit for electroweak-boson polarisation in the VBF $H\to WW$ channel, using (aL, aT) parametrisation. Both shape and rate informations are exploited in the fit. The theoretical and experimental uncertainties are subdivided further into categories.

The contributions of the leading individual systematic uncertainties together with the data statistical uncertainties, in the one dimensional fit for electroweak-boson polarisation in the VBF $H\to WW$ channel, using (aL, aT) parametrisation.. Both shape and rate informations are exploited in the fit. The theoretical and experimental uncertainties are subdivided further into categories.

Upper limits on $\mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ at 95% CL including the BDT cut as a function of the signal mass.

Upper limits on $\mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ at 95% CL without the BDT cut as a function of the signal mass.

The number of events when applying each cut in the analysis for $m_a$ = 16 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 20 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 30 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 40 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 50 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

The number of events when applying each cut in the analysis for $m_a$ = 60 GeV. The first row shows the approximate total number of events scaled to the theoretical signal cross-section and luminosity of the full Run 2 dataset, assuming $ \mathcal{B}(H\rightarrow aa\rightarrow bb\mu\mu)$ = 0.16%.

Covariances $c_{ij}$ for the highest ranked systematic uncertainties in the energy asymmetry measurement differential in $\theta_j$.

Definition of the fiducial phase space with lepton candidate ($\ell$), hadronic top candidate ($j_h$), leptonic top candidate ($j_l$) and associated jet candidate ($j_a$).

Bounds on individual Wilson coefficients $C($TeV$/\Lambda)^2$ from the energy asymmetry, obtained from a combined fit to its values in all three $\theta_j$ bins.

Forward neutron single spin asymmetries in p+p, p+Al and p+Au collisions as a function of pT in bins of xF correlated with hard collision activity

Forward neutron single spin asymmetries in p+p, p+Al and p+Au collisions as a function of pT in bins of xF without hard collision activity

Forward neutron single spin asymmetries in p+p collisions as a function of xF in bins of pT (anti-)correlated with hard collision activity

Forward neutron single spin asymmetries in p+Al collisions as a function of xF in bins of pT (anti-)correlated with hard collision activity

Forward neutron single spin asymmetries in p+Au collisions as a function of xF in bins of pT (anti-)correlated with hard collision activity

Data from Figure 11a

Model-independent per-channel efficiencies calculated in the fiducial volume of the High Mass channel for the $2e2\mu$ final state.

Model-independent per-channel efficiencies calculated in the fiducial volume of the High Mass channel for the $4e$ final state.

Model-dependent per-channel fiducial region acceptances for the $H\to aa\to 4\mu$ process. The regions of 2 GeV to 4.4 GeV, and 8 GeV to 12 GeV are each excluded by the quarkonia veto. Both, the Low Mass and the High Mass channel, include $m_X = $ 15 GeV at their signal region's edges, hence the dual entry for $m_X = $ 15 GeV.

Model-dependent per-channel fiducial region acceptances for the $H\to Z_dZ_d\to 4\ell$ process. The regions of 2 GeV to 4.4 GeV, and 8 GeV to 12 GeV are each excluded by the quarkonia veto. Both, the Low Mass and the High Mass channel, include $m_X = $ 15 GeV at their signal region's edges, hence the dual entry for $m_X = $ 15 GeV.

Model-dependent per-channel fiducial region acceptances for the $H\to Z_dZ_d\to 2e2\mu$ process.

Model-dependent per-channel fiducial region acceptances for the $H\to Z_dZ_d\to 4e$ process.

Upper limits at 95% CL on fiducial cross-sections for the $H\rightarrow XX \rightarrow 4\mu$ process. The step change at $m_X$ = 15 GeV is due to the change in efficiency caused by the change in fiducial phase-space definition. The regions of 2 GeV to 4.4 GeV, and 8 GeV to 12 GeV are each excluded by the quarkonia veto.

Upper limits at 95% CL on fiducial cross-sections for the $H\rightarrow XX \rightarrow 4e$ process.

Upper limits at 95% CL on fiducial cross-sections for the $H\rightarrow XX \rightarrow 2e2\mu$ process.

Observed and expected upper limits at 95% CL for the cross section of the $H\to Z_dZ_d\to 4\ell$ process, assuming SM Higgs boson production via the gluon-gluon fusion process. All final states are combined.

Observed and expected upper limits at 95% CL for the cross sections of the $H\to Z_dZ_d\to 4\mu$ processes, assuming SM Higgs boson production via the gluon-gluon fusion process. The regions of 2 GeV to 4.4 GeV, and 8 GeV to 12 GeV are each excluded by the quarkonia veto. The step changes at $m_{Z_d}= 15 \,\text{GeV}$ are due to the change in selection from the Low Mass to the High Mass analysis.

Observed and expected upper limits at 95% CL for the cross sections of the $H\to aa\to 4\mu$ processes, assuming SM Higgs boson production via the gluon-gluon fusion process. The regions of 2 GeV to 4.4 GeV, and 8 GeV to 12 GeV are each excluded by the quarkonia veto. The step changes at $m_{Z_d}= 15 \,\text{GeV}$ are due to the change in selection from the Low Mass to the High Mass analysis.

Model-independent efficiencies for the $H\to ZX$ process for the $4\mu$ and $2e2\mu$ final states combined calculated in the $ZX$ fiducial volume.

Model-independent efficiencies for the $H\to ZX$ process for the $4\mu$, $2e2\mu$, $2\mu 2e$, and $4e$ final states combined calculated in the $ZX$ fiducial volume.

Model-independent efficiencies for the $H\to ZX$ process for the $2\mu 2e$, and $4e$ final states combined calculated in the $ZX$ fiducial volume.

Model-dependent per-channel fiducial region acceptances for the $H\to ZZ_d\to 4\ell$ process for the $4\mu$, and $2e2\mu$ final states.

Model-dependent per-channel fiducial region acceptances for the $H\to ZZ_d\to 4\ell$ process for the $4\mu$, $2e2\mu$, $2\mu 2e$, and $4e$ final states.

Model-dependent per-channel fiducial region acceptances for the $H\to ZZ_d\to 4\ell$ process for the $2\mu 2e$, and $4e$ final states.

Per-channel upper limit at 95% CL on the fiducial cross-section for the $H\rightarrow ZX \to 2\ell2e$ process.

Per-channel upper limit at 95% CL on the fiducial cross-section for the $H\rightarrow ZX \to 2\ell2\mu$ process.

Observed and expected upper limits at 95% CL for the cross sections of the $H\to ZZ_d\to 4\ell$ process, assuming SM Higgs~boson production via the gluon-gluon fusion process.

Observed and expected upper limits at 95% CL for the cross sections of the $H\to Za\to 2\ell2\mu$ process, assuming SM Higgs~boson production via the gluon-gluon fusion process.

95% CL upper limits on the cross section times the model-dependent branching ratio divided by the SM Higgs boson production cross section for the $H\rightarrow Z_dZ_d$ process for the benchmark HAHM. The regions of 2 GeV to 4.4 GeV, and 8 GeV to 12 GeV are each excluded by the quarkonia veto. The step changes at $m_{Z_d}= 15 \,\text{GeV}$ are due to the change in selection from the Low Mass to the High Mass analysis.

95% CL upper limits on the cross section times the model-dependent branching ratio divided by the SM Higgs boson production cross section for the $H\rightarrow aa$ process for the benchmark 2HDM+S model. The regions of 2 GeV to 4.4 GeV, and 8 GeV to 12 GeV are each excluded by the quarkonia veto. The step changes at $m_{Z_d}= 15 \,\text{GeV}$ are due to the change in selection from the Low Mass to the High Mass analysis.

Upper limit at 95% CL on the effective Higgs mixing parameter $\kappa^\prime=\kappa \, m_H^2 / |m_H^2-m_S^2|$. The regions of 2 GeV to 4.4 GeV, and 8 GeV to 12 GeV are each excluded by the quarkonia veto. The step changes at $m_{Z_d}= 15 \,\text{GeV}$ are due to the change in selection from the Low Mass to the High Mass analysis.

Upper limit at 95% CL on the $Z_d$ kinetic mixing parameter $\epsilon$, assuming the SM Higgs boson production cross section.

Upper limit at 95% CL on the $Z$-$Z_d$ mass mixing parameter $\delta^2 \times \textrm{BR}(Z_d\to \ell\ell)$, assuming the SM Higgs boson production cross section.

Exclusion limit for BrHZX_BrXee

Exclusion limit for BrHZX_BrXmumu

Exclusion limit for BrHZX_BrXll

Exclusion limit for cah

Exclusion limit for cZh

Exclusion limit for kappa