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

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1400-1500 MeV.

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1500-1600 MeV.

$M_{p\eta}$ invariant mass distribution under the conditions: $E_\gamma$ = 1400-1500 MeV; 1120 MeV $\leq M_{p\pi^0} \leq$ 1220 MeV; $\theta_p^{lab} \leq$ 25°; 25° $\leq \theta_\gamma^{lab} \leq$ 155°.

$M_{p\eta}$ invariant mass distribution under the conditions: $E_\gamma$ = 1400-1500 MeV; 1120 MeV $\leq M_{p\pi^0} \leq$ 1220 MeV; $\theta_p^{lab} \leq$ 25°; 1° $\leq \theta_\gamma^{lab} \leq$ 155°.

$M_{p\eta}$ invariant mass distribution for $E_\gamma$ = 1400-1500 MeV and $M_{p\pi^0} \leq$ 1190 MeV, allowing for all proton and photon laboratory angles.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1170 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1180 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1200 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1210 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1220 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1400-1420 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1460 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1460-1500 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1500-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1540-1600 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

$M_{p\eta}$ invariant mass distribution for the incident photon energy range $E_\gamma$ = 1420-1540 MeV and $M_{p\pi^0} \leq$ 1190 MeV.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{p\eta}$.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{\pi^0\eta}$.

Excess yield relative to the PWA distributions in the invariant mass distributions $M^2_{p\pi^0}$, with a cut on the structure in the $M^2_{p\eta}$ invariant mass distribution.

Excess yield relative to the PWA distributions in the $p-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $p-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $\pi^0-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Excess yield relative to the PWA distributions in the $\pi^0-\eta$ opening angle distribution in the $\gamma p$ centre-of-mass system.

Differential cross section $d\sigma/dM_{\pi^0\eta}$ for $E_\gamma$ = 1540-1600 MeV.

$M_{p\eta}$ invariant mass distribution, taken as difference to the PWA calculations (s. Fig. 11).

$M_{\pi^0\eta}$ invariant mass distribution, taken as difference to the PWA calculations (s. Fig. 11).

$\theta_{p\eta}$ opening angle distribution, taken as difference to the PWA calculations (s. Fig. 12).

$\theta_{\pi^0\eta}$ opening angle distribution, taken as difference to the PWA calculations (s. Fig. 12).

Correlation coefficients $\rho_{i,j}$ for the statistical and systematic uncertainties between the $i$-th and $j$-th bin of the differential $A_E$ measurement as a function of the jet scattering angle $\theta_j$

The effects on the energy asymmetry of $1\sigma$ variations in its influencing nuisance parameters for the three $\theta_j$ bins. These are extracted from the samples of the posterior distribution with $\sigma_i^{(j)} = c_{ij}/\sqrt{c_{jj}}$ being the estimated shift of bin $i$ in conjunction with a shift $\Delta\theta_j$ of nuisance parameter $j$. The data statistical (Data stat.) uncertainty is obtained from running the unfolding with all nuisance parameters being fixed to their post-marginalised values, the MC statistical uncertainty on the response matrix ($t\bar{t}$ response MC stat.) is evaluated using a bootstrapping method from the covariance matrix of the ensemble of repeated unfolding results with varied response matrices. The $\gamma$ variations denote the statistical uncertainties of the background predictions in the corresponding bin of the $\Delta E$ vs $\theta_{j}$ distribution. The numbers appended to the $W$+jets PDF variations denote the corresponding NNPDF3.0 PDF sets.

The effects on the energy asymmetry of $1\sigma$ variations in its influencing nuisance parameters for the three $\theta_j$ bins. These are extracted from the samples of the posterior distribution with $\sigma_i^{(j)} = c_{ij}/\sqrt{c_{jj}}$ being the estimated shift of bin $i$ in conjunction with a shift $\Delta\theta_j$ of nuisance parameter $j$. The data statistical (Data stat.) uncertainty is obtained from running the unfolding with all nuisance parameters being fixed to their post-marginalised values, the MC statistical uncertainty on the response matrix ($t\bar{t}$ response MC stat.) is evaluated using a bootstrapping method from the covariance matrix of the ensemble of repeated unfolding results with varied response matrices. The $\gamma$ variations denote the statistical uncertainties of the background predictions in the corresponding bin of the $\Delta E$ vs $\theta_{j}$ distribution. The numbers appended to the $W$+jets PDF variations denote the corresponding NNPDF3.0 PDF sets.

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

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

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.

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.

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

Data/MC ratio as a function of reconstructed $\pi^0$ momentum for the 2$\gamma$0p selection

Energy spectra for the 1$\gamma$1p selected events. The figure shows the unconstrained background prediction and breakdowns as a function of reconstructed shower energy.

Energy spectra for the 1$\gamma$1p selected events. This figure shows the total background prediction with systematic uncertainty both before and after the 2$\gamma$ constraint.

Energy spectra for the 1$\gamma$0p selected events. The figure shows the unconstrained background prediction and breakdowns as a function of reconstructed shower energy.

Energy spectra for the 1$\gamma$0p selected events. This figure shows the total background prediction with systematic uncertainty both before and after the 2$\gamma$ constraint.

The observed event rate for the 1$\gamma$1p event sample, and comparisons to unconstrained (left) and constrained (right) background and LEE model predictions.

The observed event rate for the 1$\gamma$0p event sample, and comparisons to unconstrained (left) and constrained (right) background and LEE model predictions.

The fractional systematic covariance matrix for the final selected 1$\gamma$1p, 1$\gamma$0p, 2$\gamma$1p, and 2$\gamma$0p event samples, in bins of reconstructed shower energy for the 1$\gamma$ distributions, and reconstructed $\pi^{0}$ momentum for the 2$\gamma$ distributions. The bin boundaries are defined in Sec. II in supplementary material. This covariance matrix is evaluated at central value prediction and all variances have been rounded to three significant digits.

The fractional systematic covariance matrix for the final selected 1$\gamma$1p, 1$\gamma$0p, 2$\gamma$1p, and 2$\gamma$0p event samples, in bins of reconstructed shower energy for the 1$\gamma$1p signal, 1$\gamma$1p background, 1$\gamma$0p signal and 1$\gamma$0p background distributions, and reconstructed $\pi^{0}$ momentum for the 2$\gamma$ distributions. The bin boundaries are defined in Sec. II in supplementary material. This covariance matrix is evaluated at central value prediction and variances have been rounded to three significant digits.

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.

$z_{g}$ for Matched Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{g}$ for HardCore Recoil jets in AuAu Data anti-kT R$=$0.4

$z_{g}$ for HardCore Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{g}$ for Matched Recoil jets in AuAu Data anti-kT R$=$0.4

$z_{g}$ for Matched Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$R_{g}$ for HardCore Trigger jets in AuAu Data anti-kT R$=$0.4

$R_{g}$ for HardCore Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$R_{g}$ for Matched Trigger jets in AuAu Data anti-kT R$=$0.4

$R_{g}$ for Matched Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$R_{g}$ for HardCore Recoil jets in AuAu Data anti-kT R$=$0.4

$R_{g}$ for HardCore Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$R_{g}$ for Matched Recoil jets in AuAu Data anti-kT R$=$0.4

$R_{g}$ for Matched Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

Fully corrected two subjet $z_{sj}$ for $15 < p_{\rm{T, jet}} < 20$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $z_{sj}$ for $20 < p_{\rm{T, jet}} < 25$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $z_{sj}$ for $25 < p_{\rm{T, jet}} < 30$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $z_{sj}$ for $30 < p_{\rm{T, jet}} < 40$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $\theta_{sj}$ for $15 < p_{\rm{T, jet}} < 20$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $\theta_{sj}$ for $20 < p_{\rm{T, jet}} < 25$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $\theta_{sj}$ for $25 < p_{\rm{T, jet}} < 30$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $\theta_{sj}$ for $30 < p_{\rm{T, jet}} < 40$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

$z_{sj}$ for HardCore Trigger jets in AuAu Data anti-kT R$=$0.4

$z_{sj}$ for HardCore Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{sj}$ for Matched Trigger jets in AuAu Data anti-kT R$=$0.4

$z_{sj}$ for Matched Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{sj}$ for HardCore Recoil jets in AuAu Data anti-kT R$=$0.4

$z_{sj}$ for HardCore Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{sj}$ for Matched Recoil jets in AuAu Data anti-kT R$=$0.4

$z_{sj}$ for Matched Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for HardCore Trigger jets in AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for HardCore Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for Matched Trigger jets in AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for Matched Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for HardCore Recoil jets in AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for HardCore Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for Matched Recoil jets in AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for Matched Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$A_{J}$ for HardCore R=0.4 anti-kT jets with $\theta_{sj} > 0.1$ in Au$+$Au collisions

$A_{J}$ for HardCore R=0.4 anti-kT jets with $\theta_{sj} > 0.1$ in pp$+$AuAu reference

$A_{J}$ for HardCore R=0.4 anti-kT jets with $0.1 < \theta_{sj} < 0.2$ in Au$+$Au collisions

$A_{J}$ for HardCore R=0.4 anti-kT jets with $0.1 < \theta_{sj} < 0.2$ in pp$+$AuAu reference

$A_{J}$ for HardCore R=0.4 anti-kT jets with $0.2 < \theta_{sj} < 0.3$ in Au$+$Au collisions

$A_{J}$ for HardCore R=0.4 anti-kT jets with $0.2 < \theta_{sj} < 0.3$ in pp$+$AuAu reference

$A_{J}$ for Matched R=0.4 anti-kT jets with $\theta_{sj} > 0.1$ in Au$+$Au collisions

$A_{J}$ for Matched R=0.4 anti-kT jets with $\theta_{sj} > 0.1$ in pp$+$AuAu reference

$A_{J}$ for Matched R=0.4 anti-kT jets with $0.1 < \theta_{sj} < 0.2$ in Au$+$Au collisions

$A_{J}$ for Matched R=0.4 anti-kT jets with $0.1 < \theta_{sj} < 0.2$ in pp$+$AuAu reference

$A_{J}$ for Matched R=0.4 anti-kT jets with $0.2 < \theta_{sj} < 0.3$ in Au$+$Au collisions

$A_{J}$ for Matched R=0.4 anti-kT jets with $0.2 < \theta_{sj} < 0.3$ in pp$+$AuAu reference

tracking efficiency in pp collisions as a function of track momentum

tracking efficiency in AuAu collisions as a function of track momentum

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.

Fourth $q^2$ moment in GeV$^8$ for the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

First $q^2$ moment in GeV$^2$ for the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Second $q^2$ moment in GeV$^4$ for the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Third $q^2$ moment in GeV$^6$ for the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Fourth $q^2$ moment in GeV$^8$ for the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Second central $q^2$ moment in GeV$^4$ for the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Third central $q^2$ moment in GeV$^6$ for the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Fourth central $q^2$ moment in GeV$^8$ for the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Second central $q^2$ moment in GeV$^4$ for the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Third central $q^2$ moment in GeV$^6$ for the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Fourth central $q^2$ moment in GeV$^8$ for the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the first $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the second $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the third $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the fourth $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the second $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the third $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the fourth $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the third and the third $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the third and the fourth $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the fourth and the fourth $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the first $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the second $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the third $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the fourth $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the second $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the third $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the fourth $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the third and the third $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the third and the fourth $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the fourth and the fourth $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the second central central $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the third central $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the fourth central $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the second central $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the third central $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the fourth central $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the third and the third central $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the third and the fourth central $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the fourth and the fourth central $q^2$ moment in the electron channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the second central central $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the third central $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the first and the fourth central $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the second central $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the third central $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the second and the fourth central $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the third and the third central $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the third and the fourth central $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

Correlation coefficients for the different lower thresholds between the fourth and the fourth central $q^2$ moment in the muon channel with lower $q^2$ thresholds ranging from $3.0$ GeV$^2$ to $10.0$ GeV$^2$ in steps of $0.5$ GeV$^2$.

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.