Fiducial differential cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $N_{\mathrm{gap}\,\mathrm{jets}}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 14.

Fiducial differential cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $\xi_{\mathrm{j}3}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 15.

Fiducial differential cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $m_{\ell\ell}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 16.

Fiducial differential cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $m_{\mathrm{T}}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 17.

Fiducial differential cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $m_{\mathrm{jj}}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 18.

Fiducial differential cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $N_{\mathrm{gap}\,\mathrm{jets}}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 19.

Fiducial differential cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $\xi_{\mathrm{j}3}$. The correlation of uncertainties of the measured cross section across bins is presented in Table 20.

Observed correlations between the bins of the LH-unfolded cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $m_{\ell\ell}$. See Table 1.

Observed correlations between the bins of the LH-unfolded cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $m_{\mathrm{T}}$. See Table 2.

Observed correlations between the bins of the LH-unfolded cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $m_{\mathrm{jj}}$. See Table 3.

Observed correlations between the bins of the LH-unfolded cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $N_{\mathrm{gap}\,\mathrm{jets}}$. See Table 4.

Observed correlations between the bins of the LH-unfolded cross section of the electroweak $W^\pm W^\pm jj$ production as a function of $\xi_{\mathrm{j}3}$. See Table 5.

Observed correlations between the bins of the LH-unfolded cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $m_{\ell\ell}$. See Table 6.

Observed correlations between the bins of the LH-unfolded cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $m_{\mathrm{T}}$. See Table 7.

Observed correlations between the bins of the LH-unfolded cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $m_{\mathrm{jj}}$. See Table 8.

Observed correlations between the bins of the LH-unfolded cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $N_{\mathrm{gap}\,\mathrm{jets}}$. See Table 9.

Observed correlations between the bins of the LH-unfolded cross section of the inclusive $W^\pm W^\pm jj$ production as a function of $\xi_{\mathrm{j}3}$. See Table 10.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label M0 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label M1 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label M7 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient. The limits on M7 were obtained without taking into account the SM-EFT interference for the EW W$^\pm$Zjj final state.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label S02 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label S1 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label T0 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label T1 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Evolution of the one-dimensional expected and observed limits at 95% CL on the parameters corresponding to the quartic operators with label T2 as a function of the cut-off scale. The unitarity bounds as a function of the cut-off scale are defined for one non-zero Wilson coefficient.

Expected and observed exclusion limits at 95% CL for $\sigma_{\mathrm{VBF}}(\mathrm{H}^{\pm\pm}_5) \times \mathcal{B}(\mathrm{H}^{\pm\pm}_5 \rightarrow W^{\pm}W^{\pm})$ as a function of the doubly-charged Higgs mass.

Expected and observed exclusion limits at 95% CL for $\sin\theta_{\mathrm{H}}$ as a function of the doubly-charged Higgs mass in the Georgi-Machacek model.

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.

Measured unfolded differential cross-section as a function of the dilepton transverse momentum $p^{ll}_{\mathrm{T}}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the the four-body transverse momentum $p^{ll\gamma\gamma}_{\mathrm{T}}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the diphoton invariant mass $m_{\gamma\gamma}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the four-body invariant mass $m_{ll\gamma\gamma}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Expected and observed $95\%$ confidence intervals for the coupling parameters $f_{T,j}/\Lambda^{4}$ of transverse dimension-8 operators. All parameter values outside of the stated range are excluded at the chosen confidence level. No unitarity constraints are applied.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,8}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,0}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,1}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,2}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,5}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,6}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,7}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,9}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

The $m_{jj}$ distributions for the CRs and the BDT classifier response distribution for the SR after the fit in all regions. The dashed line shows the total background distribution before the fit. The vertical error bars on the data points correspond to the data's statistical uncertainty. Overflows are included in the last bin. The uncertainty band corresponds to the combination of the MC statistical uncertainty and systematic uncertainties obtained in the fit.

The $m_{jj}$ distributions for the CRs and the BDT classifier response distribution for the SR after the fit in all regions. The dashed line shows the total background distribution before the fit. The vertical error bars on the data points correspond to the data's statistical uncertainty. Overflows are included in the last bin. The uncertainty band corresponds to the combination of the MC statistical uncertainty and systematic uncertainties obtained in the fit.

The $m_{jj}$ distributions for the CRs and the BDT classifier response distribution for the SR after the fit in all regions. The dashed line shows the total background distribution before the fit. The vertical error bars on the data points correspond to the data's statistical uncertainty. Overflows are included in the last bin. The uncertainty band corresponds to the combination of the MC statistical uncertainty and systematic uncertainties obtained in the fit.

Observed and expected event yields for the signal and all of the background processes considered in this analysis after the fit to the data in all of the regions. The uncertainty on the expected yield is the combination of statistical and systematic uncertainties obtained in the fit. The individual uncertainties can be correlated and do not necessarily add in quadrature to equal the total background uncertainty.

Summary of the event yield for processes in all regions, after the fit over all regions. The dashed line shows the total background distribution before the fit. The vertical error bars on the data points correspond to the data's statistical uncertainty. The uncertainty band corresponds to the combination of the MC statistical uncertainty and systematic uncertainties obtained in the fit.

Impact of different components of systematic uncertainty on the measured cross section, without taking into account the correlations. The impact is calculated by fixing the value of the corresponding nuisance parameters to the values obtained in the fit used to measure the cross section, performing the fit, estimating the signal strength uncertainty and subtracting it from the nominal uncertainty in quadrature.

Fitted POI values for this analysis, the previous ATLAS analysis, and their combination. The first and second columns present the values obtained in the individual analyses. The third column presents the values obtained in the combination.

Observed and expected one-dimensional limits on dimension 8 aQGC parameters. Limits are obtained by setting all aQGCs parameters except one to zero. Unitarity is not preserved.

Observed and expected one-dimensional limits on dimension 8 aQGC parameters in the region, where unitarity is preserved. Cutoff scales in MC, $E_{c}$, are given for each parameter. Limits are obtained by setting all aQGCs parameters except one to zero.

The $E_{T}^{\gamma}$ distribution in the SR after the fit in the control regions. The red (green) line shows the expected number of events in the case of non-zero EFT coefficient $f_{T0}/\Lambda^4$ ($f_{M0}/\Lambda^4$) with the value shown in the legend. The vertical error bars on the data points correspond to the data statistical uncertainty. Overflows are included in the last bin. The uncertainty band corresponds to the combination of the MC statistical uncertainty and systematic uncertainties obtained in the fit.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{T0}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 600$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{T5}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 600$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{T8}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 600$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{T9}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 600$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{M0}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 400$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{M1}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 400$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Evolution of the expected (red line) and observed (blue line) limits versus $E_{c}$ values for $f_{M2}/\Lambda^4$. The unitarity bound is shown by the black line. The $E_{c}<4$ TeV regime was obtained with $E_{T}^{\gamma} > 400$ GeV. The $E_{c} \geq 4$ TeV regime was obtained with $E_{T}^{\gamma} > 900$ GeV.

Correlation matrix of the parameters for the fit in all regions. Only parameters with absolute value of the correlation coefficient larger than 10 are shown.

Distribution of ${E}_{T}^{miss}$ after the fit of the multijet backgrounds, in the electron channel, in the $t\bar{t}$ 3-jets VR, without applying the cut on ${E}_{T}^{miss}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of ${E}_{T}^{miss}$ after the fit of the multijet backgrounds, in the electron channel, in the $t\bar{t}$ 4-jets VR, without applying the cut on ${E}_{T}^{miss}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the signal region, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $W$+jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $t\bar{t}$ 3-jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $t\bar{t}$ 4-jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Expected distributions of the MEM discriminant $P(S|X)$ in the SR, for the s-channel single-top signal, and for the $t\bar{t}$ and $W$+jets backgrounds, for MEM discriminant values larger than $2.0\times10^{-4}$. Each distribution is normalised to unity. The binning is the same as the optimised binning used in the signal extraction fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $W$+jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty band indicates the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $t\bar{t}$ 3-jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty band indicates the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $t\bar{t}$ 4-jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty bands indicate the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR before the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$. The lower panel of the figure shows the ratio of the data to the prediction, with different vertical axis ranges. The uncertainty band indicates the total uncertainties and their correlations in each bin. The uncertainties in the $t\bar{t}$ and $W$+jets normalisation factors, as well as in the s-channel signal cross-section, are not defined pre-fit and therefore not included. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR after the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$. The lower panel of the figure shows the ratio of the data to the prediction, with different vertical axis ranges. The uncertainty band indicates the total uncertainties and their correlations in each bin. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR after the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$, after subtraction of all backgrounds. The fitted distribution for the simulation of the signal is shown together with the post-fit uncertainty in the backgrounds. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Pre-fit and post-fit event yields in the SR, for MEM discriminant values larger than $2.0\times10^{-4}$. The central value of the event yield for each process is calculated by summing the values of the discriminant bin contents, using the nominal expected yield for the pre-fit value, and the best-fit estimate for the post-fit value. The error includes statistical and systematic uncertainties summed in quadrature. All sources of systematic uncertainties are included, taking into account correlations and anti-correlations in the post-fit case. The uncertainties in the $t\bar{t}$ and $W$+jets normalisation factors, as well as in the s-channel signal cross-section, are not defined pre-fit and therefore only included in the post-fit uncertainties.

Observed impact of the different sources of uncertainty on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. The 'Systematic uncertainties' category combines all sources of systematic uncertainties. The statistical uncertainty is obtained by repeating the fit after having fixed all nuisance parameters, including the $t\bar{t}$ and $W$+jets normalisation factors. 'Total' gives the total uncertainty on the measurement.

Observed impact of the different sources of $t\bar{t}$ modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, and 'ME/PS matching' to the matching of the ME to the parton shower.

Observed impact of the different sources of s-channel modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, as described in Section 7 in the paper.

Observed impact of the different sources of t-channel modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, as described in Section 7 in the paper.

Observed impact of the different sources of $tW$ modelling uncertainty on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, and '$t\bar{t}$ overlap' to the algorithm removing the overlap between $tW$ and $t\bar{t}$ production at NLO, as described in Section 7 in the paper.

Observed impact of the different sources of PDF uncertainties on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root.

Comparison between data and prediction after the fit to data in the signal region for the leading-jet $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the leading-jet $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the subleading-jet $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the subleading-jet $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the lepton $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the lepton $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the ${E}_{T}^{miss}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the $m_{T}^{W}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Nuisance parameters ranked according to their post-fit impacts on the best-fit value of the ratio $\mu$ of the measured cross-section to the predicted cross-section. In the figure, only the 20 nuisance parameters with the largest post-fit impacts are shown. The empty (solid) blue rectangles illustrate the pre-fit (post-fit) impact on $\mu$, corresponding to the upper axis. The pre-fit (post-fit) impact of each nuisance parameter, $\Delta\mu$, is calculated as the difference in the fitted value of $\mu$ between the nominal fit and the fit when fixing the corresponding nuisance parameter to $\hat{\theta}\pm\Delta\theta$ ($\hat{\theta}\pm\Delta\hat{\theta}$), where $\hat{\theta}$ is the best-fit value of the nuisance parameter and $\Delta\theta$ ($\Delta\hat{\theta}$) is its pre-fit (post-fit) uncertainty. Several systematic uncertainties are split into different nuisance parameters, which are indicated by NP. JES (JER) indicates jet energy scale (resolution), and $\gamma$ indicates a nuisance parameter associated to the MC statistics in one of the 18 bins numbered from 0 to 17. The black points show the best-fit values of the nuisance parameters, with the error bars representing the post-fit uncertainties. Each nuisance parameter is shown wrt. its nominal value, $\theta_0$, and in units of its pre-fit uncertainty, except the free-floating normalisation factors of the $t\bar{t}$ and $W$+jets backgrounds, and the parameters associated to the MC statistics in each bin, for which the post-fit values and uncertainties are shown.

Correlation matrix of the nuisance parameters and of the ratio $\mu$ of the measured cross-section to the predicted cross-section. The correlations are given after the fit to data. In the figure, only the parameters which have a correlation of at least 0.2 with any other parameter are shown.

Distribution of the MEM discriminant $P(S|X)$ in the SR for MEM discriminant values larger than $2.0\times10^{-4}$, for the collision data used for the measurement, and for 1000 pseudo-data replicas, generated using a bootstrapping technique, in order to assess the statistical correlations between this measurement and others, for the purpose of combinations. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the <a href="https://zenodo.org/record/5361038">BootstrapGenerator</a> software package , which implements a technique described in <a href="https://cds.cern.ch/record/2759945/">ATL-PHYS-PUB-2021-011</a>. The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. Each pseudo-data replica is assigned an index, ranging from 0 to 999, corresponding to the random number index used consistently for each observed data event.

Measured values of the signal cross-section and of the $t\bar{t}$ and $W$+jets normalisation factors, obtained by statistical-only fits to the collision data used for the measurement, and to 1000 pseudo-data replicas, generated using a bootstrapping technique, in order to assess the statistical correlations between this measurement and others, for the purpose of combinations. The central values and their statistical uncertainties are obtained by repeating the fit after having fixed all nuisance parameters, except the $t\bar{t}$ and $W$+jets normalisation factors, which are let free-floating (unlike for the statistical uncertainty on the cross-section quoted in the paper). The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the <a href="https://zenodo.org/record/5361038">BootstrapGenerator</a> software package , which implements a technique described in <a href="https://cds.cern.ch/record/2759945/">ATL-PHYS-PUB-2021-011</a>. The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. Each pseudo-data replica is assigned an index, ranging from 0 to 999, corresponding to the random number index used consistently for each observed data event.

Expected 2-D likelihood scan of the signal strength $\mu_{\tau\tau}$ versus the $CP$-mixing angle $\phi_{\tau}$.

Free-floating parameters in the measurement.

Summary of results for cross-section of non-prompt $\psi(2S)$ decaying to a muon pair for 13 TeV data in fb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $J/\psi$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $\psi(2S)$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Correction factors for an assumed angular dependence $\propto 1+\lambda_{\theta} \cos^2\theta $ in the helicity frame, for several values of the parameter $\lambda_{\theta}$. The correction factors were found to be the same (within 1% - 2%) for $J\psi$ and $\psi(2S)$ mesons, for prompt and non-prompt production mechanisms, and for the three rapidity intervals considered in this paper.

Correlation matrix for the measurement of STXS parameters in each of the 28 regions considered.

Fitted values for kappa_g and kappa_y.

Correlation matrix for the measurement of kappa_g and kappa_y.

Summary of the 68% CL confidence intervals for individual measurements of SMEFT parameters observed in data. In each case, SMEFT parameters other than the one measured are fixed to 0.

Results of the EV parameter measurement in data, in the linear and linear+quadratic parameterizations of the SMEFT. All the EVs parameters are free to vary in the fits. The ranges correspond to 68% CL confidence intervals.

Observed linear correlation coefficients of the EVs parameters in the linear parameterization.

Observed linear correlation coefficients of the EVs parameters in the linear+quadratic parameterization.

Cross-sections times H->yy branching ratio for ggF +bbH, VBF, VH, ttH, and tH production. The values are obtained from a simultaneous fit to all categories. The theory uncertainties in the predictions include uncertainties due to missing higher-order terms in the perturbative QCD calculations and choices of parton distribution functions and value of alpha_s, as well as the H->yy branching ratio uncertainty.

Best-fit values and uncertainties for STXS parameters in each of the 28 regions considered. The values for the gg->H process also include the contributions from bbH production.

Correlation matrix for the measurement of STXS parameters in each of the 33 regions considered.

The rho-parameter, i.e. the ratio of the real to imaginary part of the elastic scattering amplitude extrapolated to t=0. The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter B from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter B from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter C from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter C from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter D from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The nuclear slope parameter D from a fit of the form exp(-Bt-Ct^2-Dt^3). The systematic uncertainty includes experimental and theoretical uncerainties.

The total elastic cross section measured inside the fiducial volume. The systematic uncertainty includes experimental uncertainties.

The total elastic cross section measured inside the fiducial volume. The systematic uncertainty includes experimental uncertainties.

The total elastic cross section obtained from the fitted parameters, extrapolated to full phase space using only the nuclear amplitude.

The total elastic cross section obtained from the fitted parameters, extrapolated to full phase space using only the nuclear amplitude.

The total inelastic cross section.

The total inelastic cross section.

The ratio of elastic to total cross section.

The ratio of elastic to total cross section.

The differential elastic cross section as function of t with statistical and systematic uncertainties. The systematic uncertainties are given as signed relative change for 20 sources of experimental uncertainty associated to nuisance parameters used in the fit for the extraction of physics parameters.

The differential elastic cross section as function of t with statistical and systematic uncertainties. The systematic uncertainties are given as signed relative change for 20 sources of experimental uncertainty associated to nuisance parameters used in the fit for the extraction of physics parameters.

Statistical covariance matrix for the measurement of the differential elastic cross section as function of t.

Statistical covariance matrix for the measurement of the differential elastic cross section as function of t.

Signal region detector-level distribution for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$.

Signal region detector-level distribution for the observable $m_{e\mu}$.

Signal region detector-level distribution for the observable $p_{\mathrm{T}}^{e\mu}$.

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$