The post-fit $m_{CT}$ distribution is shown in the validation region VR-offLM after all the selection requirements are applied other than the $m_{CT}$ cut. The stacked histograms show the expected SM backgrounds. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. The red line with arrow indicates the $m_{CT}$ cut used in SR selection. The first and the last bin include the underflow and overflow events (where present), respectively.

The post-fit $m_{CT}$ distribution is shown in the validation region VR-offMM after all the selection requirements are applied other than the $m_{CT}$ cut. The stacked histograms show the expected SM backgrounds. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. The red line with arrow indicates the $m_{CT}$ cut used in SR selection. The first and the last bin include the underflow and overflow events (where present), respectively.

The post-fit $m_{CT}$ distribution is shown in the validation region VR-offHM after all the selection requirements are applied other than the $m_{CT}$ cut. The stacked histograms show the expected SM backgrounds. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. The red line with arrow indicates the $m_{CT}$ cut used in SR selection. The first and the last bin include the underflow and overflow events (where present), respectively.

The post-fit $m_{CT}$ distribution for SR-HM. The stacked histograms show the expected SM backgrounds. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distribution of the SUSY reference points are also shown as dashed lines. The overflow events, where present, are included in the last bin.

The post-fit $m_{CT}$ distribution for SR-MM. The stacked histograms show the expected SM backgrounds. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distribution of the SUSY reference points are also shown as dashed lines. The overflow events, where present, are included in the last bin.

The post-fit $m_{CT}$ distribution for SR-LM. The stacked histograms show the expected SM backgrounds. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distribution of the SUSY reference points are also shown as dashed lines. The overflow events, where present, are included in the last bin.

The post-fit $m_{bb}$ distribution is shown in the signal region SR-HM after all the selection requirements are applied other than the $m_{bb}$ cut. The stacked histograms show the expected SM backgrounds. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distribution of the SUSY reference points are also shown as dashed lines. The red line with arrow indicates the $m_{bb}$ cut used in SR selection.The overflow events, where present, are included in the last bin.

The post-fit $m_{bb}$ distribution is shown in the signal region SR-MM after all the selection requirements are applied other than the $m_{bb}$ cut. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distribution of the SUSY reference points are also shown as dashed lines. The red line with arrow indicates the $m_{bb}$ cut used in SR selection. The overflow events, where present, are included in the last bin.

The post-fit $m_{bb}$ distribution is shown in the signal region SR-LM after all the selection requirements are applied other than the $m_{bb}$ cut. The hatched bands represent the sum in quadrature of systematic and statistical uncertainties of the total SM background. For illustration, the distribution of the SUSY reference points are also shown as dashed lines. The red line with arrow indicates the $m_{bb}$ cut used in SR selection. The overflow events, where present, are included in the last bin.

The observed exclusion for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production. Experimental and theoretical systematic uncertainties are applied to background and signal samples and illustrated by the yellow band and the red dotted contour lines, respectively. The red dotted lines indicate the $\pm$ 1 standard-deviation variation on the observed exclusion limit due to theoretical uncertainties in the signal cross-section.

The observed exclusion up limit for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production. The red dotted lines indicate the $\pm 1 \sigma$ on the observed exclusion limit due to the theoretical uncertainties in the signal cross-section.

The observed exclusion down limit for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production. The red dotted lines indicate the $\pm 1 \sigma$ on the observed exclusion limit due to the theoretical uncertainties in the signal cross-section.

The expected exclusion for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production. Experimental and theoretical systematic uncertainties are applied to background and signal samples and illustrated by the yellow band and the red dotted contour lines, respectively. The red dotted lines indicate the $\pm$ 1 standard-deviation variation on the observed exclusion limit due to theoretical uncertainties in the signal cross-section.

Upper limits on the cross sections for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-LM for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-LM low $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-LM med. $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-LM high $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-MM for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-MM low $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-MM med. $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-MM high $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-HM for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production. 1lb\bar{b}$ production

Signal acceptance in SR-HM low $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-HM med. $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal acceptance in SR-HM high $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-LM for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-LM low $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-LM med. $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-LM high $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-MM for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-MM low $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-MM med. $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-MM high $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-HM for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-HM low $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-HM med. $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Signal efficiency in SR-HM high $m_{CT}$ for simplified models with $\tilde\chi^\pm_1 \tilde\chi^0_2 \rightarrow Wh\tilde\chi^0_1\tilde\chi^0_1, W \rightarrow l\nu, h \rightarrow b\bar{b}$ production.

Event selection cutflow for a representative signal sample for the SR-LM low $m_{CT}$. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the SR-LM med. $m_{CT}$. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the SR-LM high $m_{CT}$. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the SR-MM low $m_{CT}$. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the SR-MM med. $m_{CT}$. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the SR-MM high $m_{CT}$. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the SR-HM low $m_{CT}$. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the SR-HM med. $m_{CT}$. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the SR-HM high $m_{CT}$. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the discovery SR-LM. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the discovery SR-MM. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Event selection cutflow for a representative signal sample for the discovery SR-HM. The masses of next-lightest-neutralinos and LSPs are reported. While the first row of the table reports the total raw MC events produced, all subsequent rows show weighted events. Only statistical uncertainties are shown. Samples are produced with generator filters which selects $h\rightarrow b\bar{b}$ and $W\rightarrow\ell\nu$ decays.

Dijet invariant mass distribution for the category with at least one b-tagged jet.

Dijet invariant mass distribution for the category with both jets b-tagged.

The 95% CL upper limits on the cross-section times acceptance times branching ratio into two jets as a function of the mass of q* signal.

The 95% CL upper limits on the cross-section times acceptance times branching ratio into two jets as a function of the mass of QBH signal.

The 95% CL upper limits on the cross-section times acceptance times branching ratio into two jets as a function of the mass of W' signal.

The 95% CL upper limits on the cross-section times acceptance times branching ratio into two jets as a function of the mass of W* signal.

The upper limits on the DM mediator Z' signal at 95% CL from the inclusive category. The 95% CL upper limits are set on the universal quark coupling $g_q$ as a function of the Z' mass.

The 95% CL upper limit on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the mass of the b* signal.

The 95% CL upper limit on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the signal mass in the DM mediator Z' with $g_q$ = 0.25 model.

The 95% CL upper limit on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the signal mass in the SSM Z' model.

The 95% CL upper limit on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the signal mass in the graviton with $k/\overline{M}_{PL}=0.2$ model.

The 95% CL upper limit on the cross-section times kinematic acceptance times branching ratio for resonances with a generic Gaussian shape, as a function of the Gaussian mean mass in the inclusive category. Different widths, from 0% up to 15% of the signal mass, are considered. Gaussian-shape signals with 0% widths correspond to signal widths smaller than the experimental resolution. For a Gaussian-shaped signal with a relative width of 15%, the limits are truncated at high mass when the broad signal starts to overlap the upper end of the mass spectrum.

The 95% CL upper limit on the cross-section times kinematic acceptance times b-tagging efficiency times branching ratio for resonances with a generic Gaussian shape, as a function of the Gaussian mean mass in the 1 b category. Different widths, from 0% up to 15% of the signal mass, are considered. Gaussian-shape signals with 0% widths correspond to signal widths smaller than the experimental resolution. For a Gaussian-shaped signal with a relative width of 15%, the limits are truncated at high mass when the broad signal starts to overlap the upper end of the mass spectrum.

The 95% CL upper limit on the cross-section times kinematic acceptance times b-tagging efficiency times branching ratio for resonances with a generic Gaussian shape, as a function of the Gaussian mean mass in the 2 b category. Different widths, from 0% up to 15% of the signal mass, are considered. Gaussian-shape signals with 0% widths correspond to signal widths smaller than the experimental resolution. For a Gaussian-shaped signal with a relative width of 15%, the limits are truncated at high mass when the broad signal starts to overlap the upper end of the mass spectrum.

The expected 95% CL upper limits on the cross-section times acceptance times b-tagging efficiency times branching ratio as a function of the DM mediator Z' mass for the current and previous iterations of the analysis. The upper limit of the previous result was obtained with the Bayesian method and is also shown scaled to the 139 fb$^{-1}$ integrated luminosity of the current result to illustrate the effect of the analysis improvements.

Acceptance for the QBH, Q* and W' benchmark signal models in the inclusive category, as a function of the simulated mass.

Acceptance for the DM Z' benchmark signal model for various $g_q$ coupling parameters in the inclusive category, as a function of the simulated mass.

Acceptance for the W* benchmark signal model in the inclusive category, as a function of the simulated mass.

Acceptance times b-tagging efficiency for the b* benchmark signal model in the 1b category, as a function of the simulated mass.

Acceptance times b-tagging efficiency for the DM Z' with $g_q$=0.25, SSM Z' and graviton with $k/\overline{M}_{PL}=0.2$ benchmark signal models in the 2b category, as a function of the simulated mass.

The expected 95% CL upper limits on the cross-section times branching ratio as a function of the DM mediator Z' mass for the current and previous iterations of the analysis. The upper limit of the previous result was obtained with the Bayesian method and is also shown scaled to the 139 fb$^{-1}$ integrated luminosity of the current result to illustrate the effect of the analysis improvements. The current b-tagging requirement is tighter than the previous one for high-$p_T$ jets, resulting in a data sample with limited size for mass above 4 TeV. The background rejection, instead, has improved significantly across the entire mass spectrum inspected by the analysis.

Normalised differential cross-section in the fiducial region as a function of lepton |eta|. The first column gives the cross-section including contributions from leptonic tau decays, the second without. Systematic uncertainties are given for ttbar modelling (ttmod), lepton calibration (lept), jet and b-tagging calibration (jet), backgrounds (bkg) and integrated luminosity and beam energy (leb). The corresponding correlation matrices are given in Tables 29 and 30.

Absolute differential cross-section in the fiducial region as a function of dilepton pT. The first column gives the cross-section including contributions from leptonic tau decays, the second without. Systematic uncertainties are given for ttbar modelling (ttmod), lepton calibration (lept), jet and b-tagging calibration (jet), backgrounds (bkg) and integrated luminosity and beam energy (leb). The last bin includes overflow beyond the upper bin boundary. The corresponding correlation matrices are given in Tables 31 and 32.

Normalised differential cross-section in the fiducial region as a function of dilepton pT. The first column gives the cross-section including contributions from leptonic tau decays, the second without. Systematic uncertainties are given for ttbar modelling (ttmod), lepton calibration (lept), jet and b-tagging calibration (jet), backgrounds (bkg) and integrated luminosity and beam energy (leb). The last bin includes overflow beyond the upper bin boundary. The corresponding correlation matrices are given in Tables 33 and 34.

Absolute differential cross-section in the fiducial region as a function of dilepton invariant mass. The first column gives the cross-section including contributions from leptonic tau decays, the second without. Systematic uncertainties are given for ttbar modelling (ttmod), lepton calibration (lept), jet and b-tagging calibration (jet), backgrounds (bkg) and integrated luminosity and beam energy (leb). The last bin includes overflow beyond the upper bin boundary. The corresponding correlation matrices are given in Tables 35 and 36.

Normalised differential cross-section in the fiducial region as a function of dilepton invariant mass. The first column gives the cross-section including contributions from leptonic tau decays, the second without. Systematic uncertainties are given for ttbar modelling (ttmod), lepton calibration (lept), jet and b-tagging calibration (jet), backgrounds (bkg) and integrated luminosity and beam energy (leb). The last bin includes overflow beyond the upper bin boundary. The corresponding correlation matrices are given in Tables 37 and 38.

Absolute differential cross-section in the fiducial region as a function of dilepton |rapidity|. The first column gives the cross-section including contributions from leptonic tau decays, the second without. Systematic uncertainties are given for ttbar modelling (ttmod), lepton calibration (lept), jet and b-tagging calibration (jet), backgrounds (bkg) and integrated luminosity and beam energy (leb). The corresponding correlation matrices are given in Tables 39 and 40.

Normalised differential cross-section in the fiducial region as a function of dilepton |rapidity|. The first column gives the cross-section including contributions from leptonic tau decays, the second without. Systematic uncertainties are given for ttbar modelling (ttmod), lepton calibration (lept), jet and b-tagging calibration (jet), backgrounds (bkg) and integrated luminosity and beam energy (leb). The corresponding correlation matrices are given in Tables 41 and 42.

Absolute differential cross-section in the fiducial region as a function of the azimuthal angle between the leptons. The first column gives the cross-section including contributions from leptonic tau decays, the second without. Systematic uncertainties are given for ttbar modelling (ttmod), lepton calibration (lept), jet and b-tagging calibration (jet), backgrounds (bkg) and integrated luminosity and beam energy (leb). The corresponding correlation matrices are given in Tables 43 and 44.

Normalised differential cross-section in the fiducial region as a function of the azimuthal angle between the leptons. The first column gives the cross-section including contributions from leptonic tau decays, the second without. Systematic uncertainties are given for ttbar modelling (ttmod), lepton calibration (lept), jet and b-tagging calibration (jet), backgrounds (bkg) and integrated luminosity and beam energy (leb). The corresponding correlation matrices are given in Tables 45 and 46.

Absolute differential cross-section in the fiducial region as a function of the sum of pT of the two leptons. The first column gives the cross-section including contributions from leptonic tau decays, the second without. Systematic uncertainties are given for ttbar modelling (ttmod), lepton calibration (lept), jet and b-tagging calibration (jet), backgrounds (bkg) and integrated luminosity and beam energy (leb). The last bin includes overflow beyond the upper bin boundary. The corresponding correlation matrices are given in Tables 47 and 48.

Normalised differential cross-section in the fiducial region as a function of the sum of pT of the two leptons. The first column gives the cross-section including contributions from leptonic tau decays, the second without. Systematic uncertainties are given for ttbar modelling (ttmod), lepton calibration (lept), jet and b-tagging calibration (jet), backgrounds (bkg) and integrated luminosity and beam energy (leb). The last bin includes overflow beyond the upper bin boundary. The corresponding correlation matrices are given in Tables 49 and 50.

Absolute differential cross-section in the fiducial region as a function of the sum of the energies of the two leptons. The first column gives the cross-section including contributions from leptonic tau decays, the second without. Systematic uncertainties are given for ttbar modelling (ttmod), lepton calibration (lept), jet and b-tagging calibration (jet), backgrounds (bkg) and integrated luminosity and beam energy (leb). The last bin includes overflow beyond the upper bin boundary. The corresponding correlation matrices are given in Tables 51 and 52.

Normalised differential cross-section in the fiducial region as a function of the sum of the energies of the two leptons. The first column gives the cross-section including contributions from leptonic tau decays, the second without. Systematic uncertainties are given for ttbar modelling (ttmod), lepton calibration (lept), jet and b-tagging calibration (jet), backgrounds (bkg) and integrated luminosity and beam energy (leb). The last bin includes overflow beyond the upper bin boundary. The corresponding correlation matrices are given in Tables 53 and 54.

Absolute differential cross-section in the fiducial region as a function of lepton |eta| and dilepton invariant mass. The first four columns give the cross-section including contributions from leptonic tau decays in the four mass bins, and the last four columns do not include the leptonic tau contributions. Systematic uncertainties are given for ttbar modelling (ttmod), lepton calibration (lept), jet and b-tagging calibration (jet), backgrounds (bkg) and integrated luminosity and beam energy (leb).

Normalised differential cross-section in the fiducial region as a function of lepton |eta| and dilepton invariant mass. The first four columns give the cross-section including contributions from leptonic tau decays in the four mass bins, and the last four columns do not include the leptonic tau contributions. Systematic uncertainties are given for ttbar modelling (ttmod), lepton calibration (lept), jet and b-tagging calibration (jet), backgrounds (bkg) and integrated luminosity and beam energy (leb).

Absolute differential cross-section in the fiducial region as a function of dilepton |rapidity| and dilepton invariant mass. The first four columns give the cross-section including contributions from leptonic tau decays in the four mass bins, and the last four columns do not include the leptonic tau contributions. Systematic uncertainties are given for ttbar modelling (ttmod), lepton calibration (lept), jet and b-tagging calibration (jet), backgrounds (bkg) and integrated luminosity and beam energy (leb).

Normalised differential cross-section in the fiducial region as a function of dilepton |rapidity| and dilepton invariant mass. The first four columns give the cross-section including contributions from leptonic tau decays in the four mass bins, and the last four columns do not include the leptonic tau contributions. Systematic uncertainties are given for ttbar modelling (ttmod), lepton calibration (lept), jet and b-tagging calibration (jet), backgrounds (bkg) and integrated luminosity and beam energy (leb).

Absolute differential cross-section in the fiducial region as a function of the azimuthal angle between the leptons and dilepton invariant mass. The first four columns give the cross-section including contributions from leptonic tau decays in the four mass bins, and the last four columns do not include the leptonic tau contributions. Systematic uncertainties are given for ttbar modelling (ttmod), lepton calibration (lept), jet and b-tagging calibration (jet), backgrounds (bkg) and integrated luminosity and beam energy (leb).

Normalised differential cross-section in the fiducial region as a function of the azimuthal angle between the leptons and dilepton invariant mass. The first four columns give the cross-section including contributions from leptonic tau decays in the four mass bins, and the last four columns do not include the leptonic tau contributions. Systematic uncertainties are given for ttbar modelling (ttmod), lepton calibration (lept), jet and b-tagging calibration (jet), backgrounds (bkg) and integrated luminosity and beam energy (leb).

Correlation matrix for the absolute differential cross-section in the fiducial region as a function of lepton pT as measured in Table 1, including tau decay contributions.

Correlation matrix for the absolute differential cross-section in the fiducial region as a function of lepton pT as measured in Table 1, not including tau decay contributions.

Correlation matrix for the normalised differential cross-section in the fiducial region as a function of lepton pT as measured in Table 2, including tau decay contributions.

Correlation matrix for the normalised differential cross-section in the fiducial region as a function of lepton pT as measured in Table 2, not including tau decay contributions.

Correlation matrix for the absolute differential cross-section in the fiducial region as a function of lepton |eta| as measured in Table 3, including tau decay contributions.

Correlation matrix for the absolute differential cross-section in the fiducial region as a function of lepton |eta| as measured in Table 3, not including tau decay contributions.

Correlation matrix for the normalised differential cross-section in the fiducial region as a function of lepton |eta| as measured in Table 4, including tau decay contributions.

Correlation matrix for the normalised differential cross-section in the fiducial region as a function of lepton |eta| as measured in Table 4, not including tau decay contributions.

Correlation matrix for the absolute differential cross-section in the fiducial region as a function of dilepton pT as measured in Table 5, including tau decay contributions.

Correlation matrix for the absolute differential cross-section in the fiducial region as a function of dilepton pT as measured in Table 5, not including tau decay contributions.

Correlation matrix for the normalised differential cross-section in the fiducial region as a function of dilepton pT as measured in Table 6, including tau decay contributions.

Correlation matrix for the normalised differential cross-section in the fiducial region as a function of dilepton pT as measured in Table 6, not including tau decay contributions.

Correlation matrix for the absolute differential cross-section in the fiducial region as a function of dilepton invariant mass as measured in Table 7, including tau decay contributions.

Correlation matrix for the absolute differential cross-section in the fiducial region as a function of dilepton invariant mass as measured in Table 7, not including tau decay contributions.

Correlation matrix for the normalised differential cross-section in the fiducial region as a function of dilepton invariant mass as measured in Table 8, including tau decay contributions.

Correlation matrix for the normalised differential cross-section in the fiducial region as a function of dilepton invariant mass as measured in Table 8, not including tau decay contributions.

Correlation matrix for the absolute differential cross-section in the fiducial region as a function of dilepton |rapidity| as measured in Table 9, including tau decay contributions.

Correlation matrix for the absolute differential cross-section in the fiducial region as a function of dilepton |rapidity| as measured in Table 9, not including tau decay contributions.

Correlation matrix for the normalised differential cross-section in the fiducial region as a function of dilepton |rapidity| as measured in Table 10, including tau decay contributions.

Correlation matrix for the normalised differential cross-section in the fiducial region as a function of dilepton |rapidity| as measured in Table 10, not including tau decay contributions.

Correlation matrix for the absolute differential cross-section in the fiducial region as a function of the azimuthal angle between the leptons as measured in Table 11, including tau decay contributions.

Correlation matrix for the absolute differential cross-section in the fiducial region as a function of the azimuthal angle between the leptons as measured in Table 11, not including tau decay contributions.

Correlation matrix for the normalised differential cross-section in the fiducial region as a function of the azimuthal angle between the leptons as measured in Table 12, including tau decay contributions.

Correlation matrix for the normalised differential cross-section in the fiducial region as a function of the azimuthal angle between the leptons as measured in Table 12, not including tau decay contributions.

Correlation matrix for the absolute differential cross-section in the fiducial region as a function of the sum of pT of the two leptons as measured in Table 13, including tau decay contributions.

Correlation matrix for the absolute differential cross-section in the fiducial region as a function of the sum of pT of the two leptons as measured in Table 13, not including tau decay contributions.

Correlation matrix for the normalised differential cross-section in the fiducial region as a function of the sum of pT of the two leptons as measured in Table 14, including tau decay contributions.

Correlation matrix for the normalised differential cross-section in the fiducial region as a function of the sum of pT of the two leptons as measured in Table 14, not including tau decay contributions.

Correlation matrix for the absolute differential cross-section in the fiducial region as a function of the sum of the energies of the two leptons as measured in Table 15, including tau decay contributions.

Correlation matrix for the absolute differential cross-section in the fiducial region as a function of the sum of the energies of the two leptons as measured in Table 15, not including tau decay contributions.

Correlation matrix for the normalised differential cross-section in the fiducial region as a function of the sum of the energies of the two leptons as measured in Table 16, including tau decay contributions.

Correlation matrix for the normalised differential cross-section in the fiducial region as a function of the sum of the energies of the two leptons as measured in Table 16, not including tau decay contributions.

Correlation matrix for the combination of all eight one-dimensional absolute fiducial cross-section measurements, including tau decay contributions. The observables are arranged as follows: lepton pT (bins 1-11),lepton |eta| (bins 12-20),dilepton pT (bins 21-29),dilepton invariant mass (bins 30-41),dilepton |rapidity| (bins 42-50),the azimuthal angle between the leptons (bins 51-60),the sum of pT of the two leptons (bins 61-68),the sum of the energies of the two leptons (bins 69-78).

Correlation matrix for the combination of all eight one-dimensional absolute fiducial cross-section measurements, not including tau decay contributions. The observables are arranged as follows: lepton pT (bins 1-11),lepton |eta| (bins 12-20),dilepton pT (bins 21-29),dilepton invariant mass (bins 30-41),dilepton |rapidity| (bins 42-50),the azimuthal angle between the leptons (bins 51-60),the sum of pT of the two leptons (bins 61-68),the sum of the energies of the two leptons (bins 69-78).

Correlation matrix for the combination of all eight one-dimensional normalised fiducial cross-section measurements, including tau decay contributions. The observables are arranged as follows: lepton pT (bins 1-10),lepton |eta| (bins 11-18),dilepton pT (bins 19-26),dilepton invariant mass (bins 27-37),dilepton |rapidity| (bins 38-45),the azimuthal angle between the leptons (bins 46-54),the sum of pT of the two leptons (bins 55-61),the sum of the energies of the two leptons (bins 62-70).

Correlation matrix for the combination of all eight one-dimensional normalised fiducial cross-section measurements, not including tau decay contributions. The observables are arranged as follows: lepton pT (bins 1-10),lepton |eta| (bins 11-18),dilepton pT (bins 19-26),dilepton invariant mass (bins 27-37),dilepton |rapidity| (bins 38-45),the azimuthal angle between the leptons (bins 46-54),the sum of pT of the two leptons (bins 55-61),the sum of the energies of the two leptons (bins 62-70).

Fiducial region definition

Distribution of $v_{3}$ from MBT events plotted as a function of the A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of $v_{3}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of $v_{3}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of $v_{2}$ from MBT events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of $v_{2}$ from MBT events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of $v_{2}$ from MBT events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-HS pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-HS pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-HS pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-HS pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-HS pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-HS pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Observed $CL_S$ limits on $\sigma \times BR$ for $m_{\Phi} = 200-400$ GeV.

Observed $CL_S$ limits on $\sigma \times BR$ for $m_{\Phi} = 600-1000$ GeV.

Combined limits from this analysis (ID) and the CR and MS analyses for $m_{H} = 125$ GeV, $m_s = 15$ GeV.

Combined limits from this analysis (ID) and the CR and MS analyses for $m_{H} = 125$ GeV, $m_s = 25$ GeV.

Combined limits from this analysis (ID) and the CR and MS analyses for $m_{H} = 125$ GeV, $m_s = 40$ GeV.

Combined limits from this analysis (ID) and the CR and MS analyses for $m_{H} = 125$ GeV, $m_s = 55$ GeV.

Combined limits from this analysis (ID) and the CR and MS analyses for $m_{\Phi} = 200$ GeV, $m_s = 25$ GeV.

Combined limits from this analysis (ID) and the CR and MS analyses for $m_{\Phi} = 200$ GeV, $m_s = 50$ GeV.

Comparison of the IDVx reconstruction and selection efficiency for a Higgs with a mass of 125 GeV decaying to an LLP with a mass of 8 GeV, using only standard tracking versus using standard and large radius tracking for all ID vertices in the signal MC sample.

Comparison of the IDVx reconstruction and selection efficiency for a Higgs with a mass of 125 GeV decaying to an LLP with a mass of 8 GeV, using only standard tracking versus using standard and large radius tracking for ID vertices passing the full IDVx selection criteria.

Comparison of the IDVx reconstruction and selection efficiency for a $\Phi$ with a mass of 1000 GeV decaying to an LLP with a mass of 150 GeV, using only standard tracking versus using standard and large radius tracking for all ID vertices in the signal MC sample.

Comparison of the IDVx reconstruction and selection efficiency for a $\Phi$ with a mass of 1000 GeV decaying to an LLP with a mass of 150 GeV, using only standard tracking versus using standard and large radius tracking for ID vertices passing the full IDVx selection criteria.

The IDVx selection efficiency as a function of long-lived particle decay $z$ position for MC signal samples with a 125 GeV Higgs boson decaying to LLPs with masses of 8, 25, and 55 GeV. The efficiency for the 55 GeV LLP sample is shown both with and without the material veto applied.

The IDVx selection efficiency as a function of long-lived particle decay $z$ position for MC signal samples with mediators of masses 200, 400, and 600 GeV decaying to LLPs with a mass of 50 GeV. The efficiency for the 200 GeV mediator sample is shown both with and without the material veto applied.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a Higgs with a mass of 125 GeV decaying to an LLP with a mass of 8 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a Higgs with a mass of 125 GeV decaying to an LLP with a mass of 8 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a Higgs with a mass of 125 GeV decaying to an LLP with a mass of 8 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a $\Phi$ with a mass of 400 GeV decaying to an LLP with a mass of 50 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a $\Phi$ with a mass of 400 GeV decaying to an LLP with a mass of 50 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a $\Phi$ with a mass of 400 GeV decaying to an LLP with a mass of 50 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a $\Phi$ with a mass of 1000 GeV decaying to an LLP with a mass of 150 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a $\Phi$ with a mass of 1000 GeV decaying to an LLP with a mass of 150 GeV.

The impact of the IDVx selections on the selection efficiency for the MC signal sample with a $\Phi$ with a mass of 1000 GeV decaying to an LLP with a mass of 150 GeV.

$CL_S$ limits on $B_{H\rightarrow ss}$ for $m_{H} = 125$ GeV, $m_s = 8$ GeV.

$CL_S$ limits on $B_{H\rightarrow ss}$ for $m_{H} = 125$ GeV, $m_s = 15$ GeV.

$CL_S$ limits on $B_{H\rightarrow ss}$ for $m_{H} = 125$ GeV, $m_s = 25$ GeV.

$CL_S$ limits on $B_{H\rightarrow ss}$ for $m_{H} = 125$ GeV, $m_s = 40$ GeV.

$CL_S$ limits on $B_{H\rightarrow ss}$ for $m_{H} = 125$ GeV, $m_s = 55$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 200$ GeV, $m_s = 8$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 200$ GeV, $m_s = 25$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 200$ GeV, $m_s = 50$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 400$ GeV, $m_s = 50$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 400$ GeV, $m_s = 100$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 600$ GeV, $m_s = 50$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 600$ GeV, $m_s = 150$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 1000$ GeV, $m_s = 50$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 1000$ GeV, $m_s = 150$ GeV.

$CL_S$ limits on $\sigma \times B_{\Phi\rightarrow ss}$ for $m_{\Phi} = 1000$ GeV, $m_s = 400$ GeV.

Observed 95% CL exclusion sensitivity for simplified models of direct wino-bino production, assuming $m(\tilde{\chi}_{2}^{0}) \times m(\tilde{\chi}_{1}^{0})<0$.

Expected 95% CL exclusion sensitivity for simplified models of direct wino-bino production, assuming $m(\tilde{\chi}_{2}^{0}) \times m(\tilde{\chi}_{1}^{0})>0$.

Observed 95% CL exclusion sensitivity for simplified models of direct wino-bino production, assuming $m(\tilde{\chi}_{2}^{0}) \times m(\tilde{\chi}_{1}^{0})>0$.

Expected 95% CL exclusion sensitivity for simplified models of VBF wino-bino production.

Observed 95% CL exclusion sensitivity for simplified models of direct wino-bino production.

Expected 95% CL exclusion sensitivity for simplified models of direct slepton production.

Observed 95% CL exclusion sensitivity for simplified models of direct slepton production.

Expected 95% CL exclusion sensitivity for simplified models of direct smuon production.

Observed 95% CL exclusion sensitivity for simplified models of direct smuon production.

Expected 95% CL exclusion sensitivity for simplified models of direct selectron production.

Observed 95% CL exclusion sensitivity for simplified models of direct selectron production.

Expected 95% CL exclusion sensitivity for simplified models of direct LH slepton production.

Observed 95% CL exclusion sensitivity for simplified models of direct LH slepton production.

Expected 95% CL exclusion sensitivity for simplified models of direct RH slepton production.

Observed 95% CL exclusion sensitivity for simplified models of direct RH slepton production.

Expected 95% CL exclusion sensitivity for simplified models of direct LH smuon production.

Observed 95% CL exclusion sensitivity for simplified models of direct LH smuon production.

Expected 95% CL exclusion sensitivity for simplified models of direct RH smuon production.

Observed 95% CL exclusion sensitivity for simplified models of direct RH smuon production.

Expected 95% CL exclusion sensitivity for simplified models of direct LH selectron production.

Observed 95% CL exclusion sensitivity for simplified models of direct LH selectron production.

Expected 95% CL exclusion sensitivity for simplified models of direct RH selectron production.

Observed 95% CL exclusion sensitivity for simplified models of direct RH selectron production.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the C1C1 higgsino process in the SR-E-high region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Efficiency for the C1C1 higgsino process in the SR-E-high region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the C1C1 higgsino process in the SR-E-low region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Efficiency for the C1C1 higgsino process in the SR-E-low region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the C1C1 higgsino process in the SR-E-med region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Efficiency for the C1C1 higgsino process in the SR-E-med region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2N1 higgsino process in the SR-E-high region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Efficiency for the N2N1 higgsino process in the SR-E-high region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2N1 higgsino process in the SR-E-low region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Efficiency for the N2N1 higgsino process in the SR-E-low region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2N1 higgsino process in the SR-E-med region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Efficiency for the N2N1 higgsino process in the SR-E-med region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2N1 higgsino process in the SR-E-1l1T region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Efficiency for the N2N1 higgsino process in the SR-E-1l1T region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2C1p higgsino process in the SR-E-high region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Efficiency for the N2C1p higgsino process in the SR-E-high region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2C1p higgsino process in the SR-E-low region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Efficiency for the N2C1p higgsino process in the SR-E-low region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2C1p higgsino process in the SR-E-med region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Efficiency for the N2C1p higgsino process in the SR-E-med region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2C1p higgsino process in the SR-E-1l1T region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Efficiency for the N2C1p higgsino process in the SR-E-1l1T region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2C1m higgsino process in the SR-E-high region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Efficiency for the N2C1m higgsino process in the SR-E-high region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2C1m higgsino process in the SR-E-low region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Efficiency for the N2C1m higgsino process in the SR-E-low region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2C1m higgsino process in the SR-E-med region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Efficiency for the N2C1m higgsino process in the SR-E-med region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2C1m higgsino process in the SR-E-1l1T region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Efficiency for the N2C1m higgsino process in the SR-E-1l1T region. Truth dilepton invariant mass is constrained to be within the range [0.5,60] GeV.

Acceptance (note the $z$-axis is in units of $10^{-3}$) for the slepton process in the SR-S region. Truth stransverse mass is constrained to be within the range [100,140] GeV.

Efficiency for the slepton process in the SR-S region. Truth stransverse mass is constrained to be within the range [100,140] GeV.

Acceptance (note the $z$-axis is in units of $10^{-3}$) for the slepton process in the SR-S-high region. Truth stransverse mass is constrained to be within the range [100,140] GeV.

Efficiency for the slepton process in the SR-S-high region. Truth stransverse mass is constrained to be within the range [100,140] GeV.

Acceptance (note the $z$-axis is in units of $10^{-3}$) for the slepton process in the SR-S-low region. Truth stransverse mass is constrained to be within the range [100,140] GeV.

Efficiency for the slepton process in the SR-S-low region. Truth stransverse mass is constrained to be within the range [100,140] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the C1C1 VBF higgsino process in the SR-VBF-high region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Efficiency for the C1C1 VBF higgsino process in the SR-VBF-high region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the C1C1 VBF higgsino process in the SR-VBF-low region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Efficiency for the C1C1 VBF higgsino process in the SR-VBF-low region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the C1C1 VBF higgsino process in the SR-VBF region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Efficiency for the C1C1 VBF higgsino process in the SR-VBF region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2N1 VBF higgsino process in the SR-VBF-high region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Efficiency for the N2N1 VBF higgsino process in the SR-VBF-high region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2N1 VBF higgsino process in the SR-VBF-low region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Efficiency for the N2N1 VBF higgsino process in the SR-VBF-low region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2N1 VBF higgsino process in the SR-VBF region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Efficiency for the N2N1 VBF higgsino process in the SR-VBF region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2C1p VBF higgsino process in the SR-VBF-high region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Efficiency for the N2C1p VBF higgsino process in the SR-VBF-high region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2C1p VBF higgsino process in the SR-VBF-low region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Efficiency for the N2C1p VBF higgsino process in the SR-VBF-low region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2C1p VBF higgsino process in the SR-VBF region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Efficiency for the N2C1p VBF higgsino process in the SR-VBF region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2C1m VBF higgsino process in the SR-VBF-high region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Efficiency for the N2C1m VBF higgsino process in the SR-VBF-high region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2C1m VBF higgsino process in the SR-VBF-low region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Efficiency for the N2C1m VBF higgsino process in the SR-VBF-low region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Acceptance (note the $z$-axis is in units of $10^{-4}$) for the N2C1m VBF higgsino process in the SR-VBF region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Efficiency for the N2C1m VBF higgsino process in the SR-VBF region. Truth dilepton invariant mass is constrained to be within the range [1,40] GeV.

Observed and Expected upper cross-section limits for wino-bino scenarios, assuming $m(\tilde{\chi}_{2}^{0}) \times m(\tilde{\chi}_{1}^{0})>0$.

Observed and Expected upper cross-section limits for wino-bino scenarios, assuming $m(\tilde{\chi}_{2}^{0}) \times m(\tilde{\chi}_{1}^{0})<0$.

Observed and Expected upper cross-section limits for higgsino scenarios.

Observed and Expected upper cross-section limits for higgsino scenarios, assuming VBF production..

Observed and Expected upper cross-section limits for direct slepton scenarios.

Observed and Expected upper cross-section limits for direct LH slepton scenarios.

Observed and Expected upper cross-section limits for direct RH slepton scenarios.

Observed and Expected upper cross-section limits for direct smuon scenarios.

Observed and Expected upper cross-section limits for direct LH smuon scenarios.

Observed and Expected upper cross-section limits for direct RH smuon scenarios.

Observed and Expected upper cross-section limits for direct selectron scenarios.

Observed and Expected upper cross-section limits for direct LH selectron scenarios.

Observed and Expected upper cross-section limits for direct RH selectron scenarios.

Number of signal events in SR-E-1L1T for the (m($\tilde{\chi}_{2}^{0}$),m($\tilde{\chi}_{1}^{0}$)) = (155 GeV, 150 GeV) Higgsino signal model at different stages of selection before and after weighting events to correspond to 140/fb.

Number of signal events in SR-E-high for the (m($\tilde{\chi}_{2}^{0}$),m($\tilde{\chi}_{1}^{0}$)) = (155 GeV, 150 GeV) Higgsino signal model at different stages of selection before and after weighting events to correspond to 140/fb.

Number of signal events in SR-E-low for the (m($\tilde{\chi}_{2}^{0}$),m($\tilde{\chi}_{1}^{0}$)) = (155 GeV, 150 GeV) Higgsino signal model at different stages of selection before and after weighting events to correspond to 140/fb.

Number of signal events in SR-E-med for the (m($\tilde{\chi}_{2}^{0}$),m($\tilde{\chi}_{1}^{0}$)) = (155 GeV, 150 GeV) Higgsino signal model at different stages of selection before and after weighting events to correspond to 140/fb.

Number of signal events in SR-S-high for the (m($\tilde{\ell}$),m($\tilde{\chi}_{1}^{0}$)) = (150 GeV, 140 GeV) Slepton signal model at different stages of selection before and after weighting events to correspond to 140/fb.

Number of signal events in SR-S-low for the (m($\tilde{\ell}$),m($\tilde{\chi}_{1}^{0}$)) = (150 GeV, 140 GeV) Slepton signal model at different stages of selection before and after weighting events to correspond to 140/fb.

Number of signal events in SR-VBF for the (m($\tilde{\chi}_{2}^{0}$),m($\tilde{\chi}_{1}^{0}$)) = (100 GeV, 95 GeV) Higgsino signal model at different stages of selection before and after weighting events to correspond to 140/fb.

The measured combined normalized differential cross-sections, divided by the bin-width, in the fiducial volume at Born level as well as a factor $k_{dressed}$ to translate from the Born particle level to the dressed particle level.

The measured combined normalized differential cross-sections, divided by the bin-width, in the fiducial volume at Born level as well as a factor $k_{dressed}$ to translate from the Born particle level to the dressed particle level.

Results of the normalized differential cross-section $1/\sigma_\mathrm{fid}\times \mathrm{d}\sigma_\mathrm{fid}/\mathrm{d}p_{T}^{ll}$ measured on born level for the $Z\rightarrow\mu\mu$ decay channel. The following naming convention is used: Stat.(Data), Stat.(MC) an Eff.(Uncor.), denote the statistical uncertainties due limited data and MC as well as the uncorrelated lepton efficiency uncertainties; Scale and Res. denote the muon momentum scale and resolution uncertainties; Muon Sag. denotes the uncertainty due to the muon sagitta bias; Eff. (Cor.), Isolation, Trigger and TTVA denote the uncertainties of the data/MC scale-factors for the correlated muon reconstruction, isolation, trigger and track-to-vertex matching efficiencies; the uncertainties due to the primary vertex z-distribution and pile-up reweighting are denoted as Z-Pos and Pile-Up, while the model and background uncertainties are summarized under Model and Bkg.. The sign-information is kept to track bin-to-bin changes.

Results of the normalized differential cross-section $1/\sigma_\mathrm{fid}\times \mathrm{d}\sigma_\mathrm{fid}/\mathrm{d}\phi_{\eta}^{*}$ measured on born level for the $Z\rightarrow\mu\mu$ decay channel. The following naming convention is used: Stat.(Data), Stat.(MC) an Eff.(Uncor.), denote the statistical uncertainties due limited data and MC as well as the uncorrelated lepton efficiency uncertainties; Scale and Res. denote the muon momentum scale and resolution uncertainties; Muon Sag. denotes the uncertainty due to the muon sagitta bias; Eff. (Cor.), Isolation, Trigger and TTVA denote the uncertainties of the data/MC scale-factors for the correlated muon reconstruction, isolation, trigger and track-to-vertex matching efficiencies; the uncertainties due to the primary vertex z-distribution and pile-up reweighting are denoted as Z-Pos and Pile-Up, while the model and background uncertainties are summarized under Model and Bkg.. The sign-information is kept to track bin-to-bin changes.

Results of the normalized differential cross-section $1/\sigma_\mathrm{fid}\times \mathrm{d}\sigma_\mathrm{fid}/\mathrm{d}p_{T}^{ll}$ measured on born level for the $Z\rightarrow ee$ decay channel. The following naming convention is used: Stat.(Data), Stat.(MC) and Eff.(Uncor.), denote the statistical uncertainties due limited data and MC as well as the uncorrelated lepton efficiency uncertainties; Scale and Res. for the electron momentum scale and resolution uncertainties; Elec. (Reco), Elec. (ID), Isolation, Trigger and Charge-ID denote the correlated uncertainties of the data/MC scale-factors for the electron reconstruction, identification, isolation, trigger and charge-identification efficiencies; The uncertainties due to the primary vertex z-distribution and pile-up reweighting are denoted as Z-Pos and Pile-Up, while the model and background uncertainties are summarized under Model and Bkg. The sign-information is kept to track bin-to-bin changes.

Results of the normalized differential cross-section $1/\sigma_\mathrm{fid} \times \mathrm{d}\sigma_\mathrm{fid}/\mathrm{d}\phi_{\eta}^{*}$ measured on born level for the $Z\rightarrow ee$ decay channel. The following naming convention is used: Stat.(Data), Stat.(MC) and Eff.(Uncor.), denote the statistical uncertainties due limited data and MC as well as the uncorrelated lepton efficiency uncertainties; Scale and Res. for the electron momentum scale and resolution uncertainties; Elec. (Reco), Elec. (ID), Isolation, Trigger and Charge-ID denote the correlated uncertainties of the data/MC scale-factors for the electron reconstruction, identification, isolation, trigger and charge-identification efficiencies; The uncertainties due to the primary vertex z-distribution and pile-up reweighting are denoted as Z-Pos and Pile-Up, while the model and background uncertainties are summarized under Model and Bkg. The sign-information is kept to track bin-to-bin changes.

Results of the normalized differential cross-section $1/\sigma_\mathrm{fid} \times \mathrm{d}\sigma_\mathrm{fid}/\mathrm{d}p_{T}^{ll}$ measured on bare level for the $Z\rightarrow\mu\mu$ decay channel. The following naming convention is used: Stat.(Data), Stat.(MC) an Eff.(Uncor.), denote the statistical uncertainties due limited data and MC as well as the uncorrelated lepton efficiency uncertainties; Scale and Res. denote the muon momentum scale and resolution uncertainties; Muon Sag. denotes the uncertainty due to the muon sagitta bias; Eff. (Cor.), Isolation, Trigger and TTVA denote the uncertainties of the data/MC scale-factors for the correlated muon reconstruction, isolation, trigger and track-to-vertex matching efficiencies; the uncertainties due to the primary vertex z-distribution and pile-up reweighting are denoted as Z-Pos and Pile-Up, while the model and background uncertainties are summarized under Model and Bkg.. The sign-information is kept to track bin-to-bin changes.

Results of the normalized differential cross-section $1/\sigma_\mathrm{fid} \times \mathrm{d}\sigma_\mathrm{fid}/\mathrm{d}\phi_{\eta}^{*}$ measured on bare level for the $Z\rightarrow\mu\mu$ decay channel. The following naming convention is used: Stat.(Data), Stat.(MC) an Eff.(Uncor.), denote the statistical uncertainties due limited data and MC as well as the uncorrelated lepton efficiency uncertainties; Scale and Res. denote the muon momentum scale and resolution uncertainties; Muon Sag. denotes the uncertainty due to the muon sagitta bias; Eff. (Cor.), Isolation, Trigger and TTVA denote the uncertainties of the data/MC scale-factors for the correlated muon reconstruction, isolation, trigger and track-to-vertex matching efficiencies; the uncertainties due to the primary vertex z-distribution and pile-up reweighting are denoted as Z-Pos and Pile-Up, while the model and background uncertainties are summarized under Model and Bkg.. The sign-information is kept to track bin-to-bin changes.

Measured combined normalized differential cross-section in the fiducial volume at Born level as well as a factor $k_{dressed}$ to translate from the Born particle level to the dressed particle.

Measured combined normalized differential cross-section in the fiducial volume at Born level as well as a factor $k_{dressed}$ to translate from the Born particle level to the dressed particle.

The distribution of events passing the selection requirements in the electron channel as a function of dilepton invariant mass $m_{ll}$ , the latter with one entry for each lepton per event. The MC signal sample is simulated using Powheg+Pythia8. The predictions of the MC signal sample together with the MC background samples are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the muon channel as a function of dilepton invariant mass $m_{ll}$, the latter with one entry for each lepton per event. The MC signal sample is simulated using Powheg+Pythia8. The predictions of the MC signal sample together with the MC background samples are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the electron channel as a function of lepton pseudorapidity $\eta$, the latter with one entry for each lepton per event. The MC signal sample is simulated using Powheg+Pythia8. The predictions of the MC signal sample together with the MC background samples are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the muon channel as a function of lepton pseudorapidity $\eta$, the latter with one entry for each lepton per event. The MC signal sample is simulated using Powheg+Pythia8. The predictions of the MC signal sample together with the MC background samples are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the electron channel as a function of dilepton transverse momentum. The MC signal sample is simulated using Powheg+Pythia8. The predictions are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the muon channel as a function of dilepton transverse momentum. The MC signal sample is simulated using Powheg+Pythia8. The predictions are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the electron channel as a function of $\phi_{\eta}^{*}$. The MC signal sample is simulated using Powheg+Pythia8. The predictions are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The distribution of events passing the selection requirements in the muon channel as a function of $\phi_{\eta}^{*}$. The MC signal sample is simulated using Powheg+Pythia8. The predictions are normalized to the integral of the data and the total experimental uncertainty of the predicted values is shown as a grey band in the ratio of the prediction to data.

The measured normalized cross section as a function of $p_{ll}$ for the electron and muon channels and the combined result as well as their ratio together with the total uncertainties, shown as a blue band. The pull distribution between the electron and muon channels, defined as the difference between the two channels divided by the combined uncorrelated uncertainty, is also shown. The $p_{ll}$ distribution is split into linear and logarithmic scales at 30 GeV.

The measured normalized cross section as a function of $\phi_{\eta}^{*}$ for the electron and muon channels and the combined result as well as their ratio together with the total uncertainties, shown as a blue band. The pull distribution between the electron and muon channels, defined as the difference between the two channels divided by the combined uncorrelated uncertainty, is also shown.

Comparison of the normalized $p_{ll}$ distributions predicted by different computations: Pythia8 with the AZ tune, Powheg+Pythia8 with the AZNLO tune, Sherpa v2.2.1 and RadISH with the Born level combined measurement. The uncertainties of the measurement are shown as vertical bars and uncertainties of the Sherpa and RadISH predictions are indicated by the coloured bands.

Comparison of the normalized $\phi_{\eta}^{*}$ distributions predicted by different computations: Pythia8 with the AZ tune, Powheg+Pythia8 with the AZNLO tune, Sherpa v2.2.1 and RadISH with the Born level combined measurement. The uncertainties of the measurement are shown as vertical bars and uncertainties of the Sherpa and RadISH predictions are indicated by the coloured bands.

Comparison of the normalized $p_{ll}$ distribution in the range $p_{ll}$ > 10 GeV. The Born level combined measurement is compared with predictions by Sherpa v2.2.1, fixed-order NNLOjet and NNLOjet supplied with NLO electroweak corrections. The uncertainties in the measurement are shown as vertical bars and the uncertainties in the predictions are indicated by the coloured bands.

The measured combined normalized differential cross-sections, divided by the bin-width, in the fiducial volume at dressed level.

The measured combined normalized differential cross-sections, divided by the bin-width, in the fiducial volume at dressed level.

Data from Fig 6d. The unfolded charged-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from Fig 6e. The unfolded $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from Fig 6f. The unfolded charged-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from Fig 7a. The unfolded all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in (300, 400, 600, 800, 1000, infinity) and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 7b. The unfolded all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in (300, 400, 600, 800, 1000, infinity) and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 7c. The unfolded all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in (300, 400, 600, 800, 1000, infinity) and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 7d. The unfolded all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in (300, 400, 600, 800, 1000, infinity) and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 7e. The unfolded all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in (300, 400, 600, 800, 1000, infinity) and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 7f. The unfolded all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in (300, 400, 600, 800, 1000, infinity) and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 8a. The unfolded all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 8b. The unfolded all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 8c. The unfolded all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 8d. The unfolded all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 8e. The unfolded all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 8f. The unfolded all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 14a. The unfolded all-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from Fig 14b. The unfolded charged-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 4b. The unfolded all-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 21b. The unfolded charged-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 5a. The unfolded $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 5b. The unfolded charged-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from Fig 14c. The unfolded all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 14d. The unfolded all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 4c. The unfolded all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 4d. The unfolded all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 5c. The unfolded all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 5d. The unfolded all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 14e. The unfolded all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 14f. The unfolded all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 4e. The unfolded all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 4f. The unfolded all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 5e. The unfolded all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 5f. The unfolded all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 14a. The unfolded all-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from Fig 14b. The unfolded charged-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 4a. The unfolded all-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 4b. The unfolded charged-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 5a. The unfolded $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 5b. The unfolded charged-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from Fig 14c. The unfolded all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 14d. The unfolded all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 4c. The unfolded all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 4d. The unfolded all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 5c. The unfolded all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 5d. The unfolded all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 14e. The unfolded all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 14f. The unfolded all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 4e. The unfolded all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 4f. The unfolded all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 5e. The unfolded all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 5f. The unfolded all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 36-40a. The unfolded all-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in (300, 400, 600, 800, 1000, infinity) and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from Fig 81-85a. The unfolded charged-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from Fig 36-40b. The unfolded all-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from Fig 81-85b. The unfolded charged-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from Fig 36-40c. The unfolded $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from Fig 81-85c. The unfolded charged-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from Fig 51-55a. The unfolded all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 101-105a. The unfolded all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 51-55b. The unfolded all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 101-105b. The unfolded all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 51-55c. The unfolded all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 101-105c. The unfolded all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 66-70a. The unfolded all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from Fig 106-110a. The unfolded all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from Fig 66-70b. The unfolded all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from Fig 106-110b. The unfolded all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from Fig 66-70c. The unfolded all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from Fig 106-110c. The unfolded all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from Fig 26-30a. The unfolded all-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from Fig 71-75a. The unfolded charged-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from Fig 26-30b. The unfolded all-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from Fig 71-75b. The unfolded charged-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from Fig 26-30c. The unfolded $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from Fig 71-75c. The unfolded charged-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from Fig 41-45a. The unfolded all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 86-90a. The unfolded all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 41-45b. The unfolded all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 86-90b. The unfolded all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 41-45c. The unfolded all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 86-90c. The unfolded all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 56-60a. The unfolded all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from Fig 101-105a. The unfolded all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from Fig 56-60b. The unfolded all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from Fig 101-105b. The unfolded all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from Fig 56-60c. The unfolded all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from Fig 101-105c. The unfolded all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from Fig 31-35a. The unfolded all-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from Fig 76-80a. The unfolded charged-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from Fig 31-35b. The unfolded all-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from Fig 76-80b. The unfolded charged-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from Fig 31-35c. The unfolded $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from Fig 76-80c. The unfolded charged-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from Fig 46-50a. The unfolded all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 91-95a. The unfolded all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 46-50b. The unfolded all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 91-95b. The unfolded all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 46-50c. The unfolded all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 91-95c. The unfolded all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from Fig 61-65a. The unfolded all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from Fig 106-110a. The unfolded all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from Fig 61-65b. The unfolded all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from Fig 106-110b. The unfolded all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from Fig 61-65c. The unfolded all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from Fig 106-110c. The unfolded all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from FigAux 6a. The extracted quark-distribution from the unfolded all-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from Fig 15a. Theextracted quark-distribution from the unfolded charged-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 6b. The extracted quark-distribution from the unfolded all-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from Fig 15b. The extracted quark-distribution from the unfolded charged-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 6c. The extracted quark-distribution from the unfolded all-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from Fig 15c. The extracted quark-distribution from the unfolded charged-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 7a. The extracted quark-distribution from the unfolded all-particle $z_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 16a. The extracted quark-distribution from the unfolded charged-particle $z_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 7b. The extracted quark-distribution from the unfolded all-particle $z_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 16b. The extracted quark-distribution from the unfolded charged-particle $z_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 7c. The extracted quark-distribution from the unfolded all-particle $z_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 16c. The extracted quark-distribution from the unfolded charged-particle $z_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 8a. The extracted quark-distribution from the unfolded all-particle $R_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 17a. The extracted quark-distribution from the unfolded charged-particle $R_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 8b. The extracted quark-distribution from the unfolded all-particle $R_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 17b. The extracted quark-distribution from the unfolded charged-particle $R_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 8c. The extracted quark-distribution from the unfolded all-particle $R_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 17c. The extracted quark-distribution from the unfolded charged-particle $R_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 6a. The extracted gluon-distribution from the unfolded all-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from Fig 15a. Theextracted gluon-distribution from the unfolded charged-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 6b. The extracted gluon-distribution from the unfolded all-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from Fig 15b. The extracted gluon-distribution from the unfolded charged-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 6c. The extracted gluon-distribution from the unfolded all-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from Fig 15c. The extracted gluon-distribution from the unfolded charged-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 7a. The extracted gluon-distribution from the unfolded all-particle $z_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 16a. The extracted gluon-distribution from the unfolded charged-particle $z_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 7b. The extracted gluon-distribution from the unfolded all-particle $z_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 16b. The extracted gluon-distribution from the unfolded charged-particle $z_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 7c. The extracted gluon-distribution from the unfolded all-particle $z_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 16c. The extracted gluon-distribution from the unfolded charged-particle $z_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 8a. The extracted gluon-distribution from the unfolded all-particle $R_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 17a. The extracted gluon-distribution from the unfolded charged-particle $R_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 8b. The extracted gluon-distribution from the unfolded all-particle $R_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 17b. The extracted gluon-distribution from the unfolded charged-particle $R_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 8c. The extracted gluon-distribution from the unfolded all-particle $R_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from Fig 17c. The extracted gluon-distribution from the unfolded charged-particle $R_g$ distribution for anti-kt R=0.8 jets with 600 < $p_T$ < 800 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. All uncertainties described in the text are shown on the data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 99a. The full covariance matrices for the all-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 100a. The full covariance matrices for the charged-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 99b. The full covariance matrices for the all-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 100b. The full covariance matrices for the charged-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 99c. The full covariance matrices for the $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 100c. The full covariance matrices for the charged-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 101a. The full covariance matrices for the all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 102a. The full covariance matrices for the all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 101b. The full covariance matrices for the all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 102b. The full covariance matrices for the all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 101c. The full covariance matrices for the all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 102c. The full covariance matrices for the all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 103a. The full covariance matrices for the all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 104a. The full covariance matrices for the all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 103b. The full covariance matrices for the all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 104b. The full covariance matrices for the all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 103c. The full covariance matrices for the all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 104c. The full covariance matrices for the all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 105a. The full covariance matrices for the all-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 106a. The full covariance matrices for the charged-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 105b. The full covariance matrices for the all-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 106b. The full covariance matrices for the charged-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 105c. The full covariance matrices for the $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 106c. The full covariance matrices for the charged-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 107a. The full covariance matrices for the all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 108a. The full covariance matrices for the all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 107b. The full covariance matrices for the all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 108b. The full covariance matrices for the all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 107c. The full covariance matrices for the all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 108c. The full covariance matrices for the all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 109a. The full covariance matrices for the all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 110a. The full covariance matrices for the all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 109b. The full covariance matrices for the all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 110b. The full covariance matrices for the all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 109c. The full covariance matrices for the all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 110c. The full covariance matrices for the all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 111a. The full covariance matrices for the all-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 112a. The full covariance matrices for the charged-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 111b. The full covariance matrices for the all-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 112b. The full covariance matrices for the charged-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 111c. The full covariance matrices for the $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 112c. The full covariance matrices for the charged-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 113a. The full covariance matrices for the all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 114a. The full covariance matrices for the all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 113b. The full covariance matrices for the all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 114b. The full covariance matrices for the all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 113c. The full covariance matrices for the all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 114c. The full covariance matrices for the all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 115a. The full covariance matrices for the all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 116a. The full covariance matrices for the all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 115b. The full covariance matrices for the all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 116b. The full covariance matrices for the all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 115c. The full covariance matrices for the all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 116c. The full covariance matrices for the all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$.

Data from FigAux 99d. The full covariance matrices for the all-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from FigAux 100d. The full covariance matrices for the charged-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from FigAux 99e. The full covariance matrices for the all-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from FigAux 100e. The full covariance matrices for the charged-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from FigAux 99f. The full covariance matrices for the $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from FigAux 100f. The full covariance matrices for the charged-particle $log_{10}(\rho^2)$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from FigAux 101d. The full covariance matrices for the all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from FigAux 102d. The full covariance matrices for the all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from FigAux 101e. The full covariance matrices for the all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from FigAux 102e. The full covariance matrices for the all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from FigAux 101f. The full covariance matrices for the all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from FigAux 102f. The full covariance matrices for the all-particle $z_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from FigAux 103d. The full covariance matrices for the all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from FigAux 104d. The full covariance matrices for the all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from FigAux 103e. The full covariance matrices for the all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from FigAux 104e. The full covariance matrices for the all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from FigAux 103f. The full covariance matrices for the all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from FigAux 104f. The full covariance matrices for the all-particle $R_g$ distribution for anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from FigAux 105d. The full covariance matrices for the all-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from FigAux 106d. The full covariance matrices for the charged-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from FigAux 105e. The full covariance matrices for the all-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from FigAux 106e. The full covariance matrices for the charged-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from FigAux 105f. The full covariance matrices for the $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from FigAux 106f. The full covariance matrices for the charged-particle $log_{10}(\rho^2)$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from FigAux 107d. The full covariance matrices for the all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from FigAux 108d. The full covariance matrices for the all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from FigAux 107e. The full covariance matrices for the all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from FigAux 108e. The full covariance matrices for the all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from FigAux 107f. The full covariance matrices for the all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from FigAux 108f. The full covariance matrices for the all-particle $z_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from FigAux 109d. The full covariance matrices for the all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from FigAux 110d. The full covariance matrices for the all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from FigAux 109e. The full covariance matrices for the all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from FigAux 110e. The full covariance matrices for the all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from FigAux 109f. The full covariance matrices for the all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from FigAux 110f. The full covariance matrices for the all-particle $R_g$ distribution for the more central of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from FigAux 111d. The full covariance matrices for the all-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 112d. The full covariance matrices for the charged-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 111e. The full covariance matrices for the all-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$.

Data from FigAux 112e. The full covariance matrices for the charged-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from FigAux 111f. The full covariance matrices for the $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from FigAux 112f. The full covariance matrices for the charged-particle $log_{10}(\rho^2)$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$(resum), measured in the resummation region, $-3.7 < log_{10}(\rho^2) < -1.7$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $log_{10}(\rho^2)$ from -4.5 to -0.5.

Data from FigAux 113d. The full covariance matrices for the all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from FigAux 114d. The full covariance matrices for the all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from FigAux 113e. The full covariance matrices for the all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from FigAux 114e. The full covariance matrices for the all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from FigAux 113f. The full covariance matrices for the all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from FigAux 114f. The full covariance matrices for the all-particle $z_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 10 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 10 evenly spaced bins in $z_g$ from 0.0 to 0.5.

Data from FigAux 115d. The full covariance matrices for the all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from FigAux 116d. The full covariance matrices for the all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 0, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from FigAux 115e. The full covariance matrices for the all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from FigAux 116e. The full covariance matrices for the all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 1, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from FigAux 115f. The full covariance matrices for the all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Data from FigAux 116f. The full covariance matrices for the all-particle $R_g$ distribution for the more forward of the two anti-kt R=0.8 jets with $p_T$ > 300 GeV, after the soft drop algorithm is applied for $\beta$ = 2, in data. The distributions are normalized to the integrated cross section, $\sigma$. Each set of 6 bins corresponds to one $p_T$ bin in {300, 400, 600, 800, 1000, infinity } and 6 bins in $r_g$ (0.06310, 0.10000, 0.15849, 0.25119, 0.39811, 0.63096, 0.80000).

Post-fit BDT distributions after the VBF event selection for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. The VBF signal is shown for $\mu = 0.73$ and $\tilde d = -0.01$. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit $m_{\tau\tau}^{\mathrm{MMC}}$ distributions in the low BDT score CR for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit Optimal Observable distributions in the low BDT score CR for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. ''Other bkg'' denotes all background contributions not listed explicitly in the legend. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit distributions of the event yields as a function of the Optimal Observable in the SR for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ analysis channel. The values of $\tilde d$, the signal strength $\mu$, the normalization of background processes, and nuisance parameters for the event yield prediction are set to those which minimize the NLL. The size of the combined statistical, experimental and theoretical uncertainties is given.

The observed $\Delta\mathrm{NLL}$ curve as a function of $\tilde d$ values. For comparison, expected $\Delta\mathrm{NLL}$ curves are also shown. The constraints on the nuisance parameters and normalization factors are first determined in a CR-only fit, and then a fit to pseudo-data corresponding to these nuisance parameters, normalization factors, and to $\tilde d=0, \mu = 1$ or $\tilde d =0, \mu = 0.73$ is performed to obtain these $\Delta\mathrm{NLL}$ curves. Moreover, a pre-fit expected $\Delta\mathrm{NLL}$ is shown, using pseudo-data corresponding to $\tilde d =0$ and $\mu = 1$ in the signal and control regions.

The expected $\Delta\mathrm{NLL}$ curves comparing the sensitivity of the fit with and without systematic uncertainties. For comparison, other curves are shown which remove the effect of jet-based systematic uncertainties, $\tau$-based systematic uncertainties, and MC statistical uncertainties.

The observed $\Delta\mathrm{NLL}$ curves for each analysis channel as a function of $\tilde d$, compared to the combined result. For the individual analysis channel $\Delta\mathrm{NLL}$ curves, only event yield information in the other SRs is used, ensuring that the Optimal Observable distributions in the other SRs do not influence the preferred value of $\tilde d$. The signal strength is constrained to be positive in these individual channel $\Delta\mathrm{NLL}$ curves. The exact value of the $p_{\mathrm{T}}$ cut on the leading lepton depends on the trigger.

Post-fit BDT distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit BDT distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit BDT distributions in the $Z\to \ell\ell$ CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit Optimal Observable distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit Optimal Observable distributions in the top-quark CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ DF channel. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Post-fit Optimal Observable distributions in the $Z\to \ell\ell$ CR for the $\tau_{\mathrm{lep}}\tau_{\mathrm{lep}}$ SF analysis channel. The size of the combined statistical, experimental, and theoretical uncertainties is given. The exact value of the $p_{\mathrm{T}}$ cut on the leptons depends on the trigger.

Post-fit distribution of weighted event yields as a function of the Optimal Observable for all four SRs combined. The contributions of the different SRs are weighted by a factor of ln(1 + S/B), where S and B are the post-fit expected numbers of signal and background events in that region, respectively. The size of the combined statistical, experimental, and theoretical uncertainties is given.

Observed and expected 95% credibility-level upper limits on the cross-section times branching ratio for the simplified dark-matter model. The table also shows the corresponding $1\sigma$ and $2\sigma$ bands for the expected limits. The limits are calculated using jets in events with at least one isolated lepton ($e$ or $\mu$) with $p_\text{T}^\ell \ge 60$ GeV.

The 95% CL upper limits on the cross-section times acceptance times efficiency times branching ratio, as a function of the resonance mass, $m_X$. The limits are calculated using jets in events with at least one isolated lepton ($e$ or $\mu$) with $p_\text{T}^\ell \ge 60$ GeV. The limits are shown for a generic Gaussian signal with the width $\sigma_{X}/m_{X} =0$. The table also shows the expected limits and the corresponding $1\sigma$ and $2\sigma$ bands. The tabulated data has been used in Figure 04 of the publication.

The 95% CL upper limits on the cross-section times acceptance times efficiency times branching ratio, as a function of the resonance mass, $m_X$. The limits are calculated using jets in events with at least one isolated lepton ($e$ or $\mu$) with $p_\text{T}^\ell \ge 60$ GeV. The limits are shown for a generic Gaussian signal with the width $\sigma_{X}/m_{X} =0.05$. The table also shows the expected limits and the corresponding $1\sigma$ and $2\sigma$ bands. The tabulated data has been used in Figure 04 of the publication.

The 95% CL upper limits on the cross-section times acceptance times efficiency times branching ratio, as a function of the resonance mass, $m_X$. The limits are calculated using jets in events with at least one isolated lepton ($e$ or $\mu$) with $p_\text{T}^\ell \ge 60$ GeV. The limits are shown for a generic Gaussian signal with the width $\sigma_{X}/m_{X} =0.10$. The table also shows the expected limits and the corresponding $1\sigma$ and $2\sigma$ bands. The tabulated data has been used in Figure 04 of the publication.

The 95% CL upper limits on the cross-section times acceptance times efficiency times branching ratio, as a function of the resonance mass, $m_X$. The limits are calculated using jets in events with at least one isolated lepton ($e$ or $\mu$) with $p_\text{T}^\ell \ge 60$ GeV. The limits are shown for a generic Gaussian signal with the width $\sigma_{X}/m_{X} =0.15$. The table also shows the expected limits and the corresponding $1\sigma$ and $2\sigma$ bands. The tabulated data has been used in Figure 04 of the publication.

Truth-level cross-sections ($\sigma$), acceptances ($A$), detector efficiencies ($\epsilon$) and the detector-level visible cross-sections ($\sigma^{vis}$) for different masses of the $\rho_T$ model. Only statistical uncertainties on the acceptances and efficiencies are shown. The acceptance is defined by requiring jets with $m_{jj} >216$ GeV and $|\eta|<2.4$, as well as isolated leptons ($e$ or $\mu$) with $p_\text{T}^\ell \ge 60$ GeV and $|\eta^\mu|<2.5$ ($|\eta^e|<2.47$, excluding of the gap region 1.37--1.52) for muons (electrons). The tabulated data has been used in Figure 05a of the publication.

Truth-level cross-sections ($\sigma$), acceptances ($A$), detector efficiencies ($\epsilon$) and the detector-level visible cross-sections ($\sigma^{vis}$) for different masses of the $W'/Z'$ SSM model. Only statistical uncertainties on the acceptances and efficiencies are shown. The acceptance is defined by requiring jets with $m_{jj} >216$ GeV and $|\eta|<2.4$, as well as leptons with $p_\text{T}^\ell \ge 60$ GeV and $|\eta^\mu|<2.5$ ($|\eta^e|<2.47$ with exclusion of the gap region 1.37-1.52) for muons (electrons). The tabulated data has been used in Figure 05b of the publication.

Truth-level cross-sections ($\sigma$), acceptances ($A$), detector efficiencies ($\epsilon$) and the detector-level visible cross-sections ($\sigma^{vis}$) for different masses of the $H^+$ model for $\tan \beta=1$. Only statistical uncertainties on the acceptances and efficiencies are shown. The acceptance correction was calculated assuming $\tan \beta=1$ and the narrow-width approximation. The acceptance is defined by requiring jets with $m_{jj} >216$ GeV and $|\eta|<2.4$, as well as leptons with $p_\text{T}^\ell \ge 60$ GeV and $|\eta^\mu|<2.5$ ($|\eta^e|<2.47$ with exclusion of the gap region 1.37-1.52) for muons (electrons). The tabulated data has been used in Figure 05c of the publication.

Truth-level cross-sections ($\sigma$), acceptances ($A$), detector efficiencies ($\epsilon$) and the detector-level visible cross-sections ($\sigma^{vis}$) for different masses of the $Z'$ (DM) model. Only statistical uncertainties on the acceptances and efficiencies are shown. The acceptance is defined by requiring jets with $m_{jj} > 216$ GeV and $|\eta|<2.4$, as well as leptons with $p_\text{T}^\ell \ge 60$ GeV and $|\eta^\mu|<2.5$ ($|\eta^e|<2.47$ with exclusion of the gap region 1.37-1.52) for muons (electrons). The tabulated data has been used in Figure 05d of the publication.