Lower limits at 95$\%$ CL on $\Lambda$ for the dimuon channel for different signal chiralities in the (constructive/destructive interference) SRs of the analysis.

Lower limits at 95$\%$ CL on $\Lambda$ for the dilepton combination for different signal chiralities in the (constructive/destructive interference) SRs of the analysis.

The dielectron data events passing the full selection.

The dimuon data events passing the full selection.

Distribution of the resonance mass in the c-tag Signal Region of the $ ce$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the b-tag Signal Region of the $ ce$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the untagged Signal Region of the $ c\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the c-tag Signal Region of the $ c\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the b-tag Signal Region of the $ c\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 0-tag Signal Region of the $ be$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 1-tag Signal Region of the $ be$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 2-tag Signal Region of the $ be$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 0-tag Signal Region of the $ b\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 1-tag Signal Region of the $ b\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

Distribution of the resonance mass in the 2-tag Signal Region of the $ b\mu$ channel for the post-fit background, the observed data, and the expected signal with $m_{LQ} = 1$ TeV.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into electrons, shown as a function of $m_{LQ}$ for the $qe$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into muons, shown as a function of $m_{LQ}$ for the $q\mu$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into electrons, shown as a function of $m_{LQ}$ for the $ce$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into muons, shown as a function of $m_{LQ}$ for the $c\mu$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into electrons, shown as a function of $m_{LQ}$ for the $be$ channel.

The observed and expected limits on the leptoquark pair production cross-section at 95% CL for $\mathcal{B}=1$ into muons, shown as a function of $m_{LQ}$ for the $b\mu$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $qe$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $q\mu$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $ce$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $c\mu$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $be$ channel.

The observed and expected limits on the leptoquark branching ratio at 95% CL, shown as a function of $m_{LQ}$ for the $b\mu$ channel.

The signal selection efficiency x acceptance summed over all signal regions, for all masses and LQ decay channels considered.

The observed and expected limits for all masses and LQ decay channels considered.

Cutflow Table in the electron channel, considering signal samples with LQ mass of 1 TeV.

Cutflow Table in the muon channel, considering signal samples with LQ mass of 1 TeV.

Event yields for a simulated $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ sample with $m_a = 17.5\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$. Cut 0 corresponds to the initial number of events. Cut 1 requires the single lepton trigger. Cut 2 requires 2 identified leptons. Cut 3 requires the Z-boson mass window. Cut 4 requires 2 reconstructed $a\to b\bar{b}$ candidates. Cut 5a requires 2 identified $a\to b\bar{b}$ candidates in the 1HPC1LPC region. Cut 6a requires the 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region to be inside the Higgs mass window. Cut 5b requires 2 identified $a\to b\bar{b}$ candidates in the 2HPC region. Cut 6b requires the 2 $a\to b\bar{b}$ candidates in the 2HPC region to be inside the Higgs mass window.

Background yield table for Z+jets, $t\bar{t}$, and rare sources. Observed data yield. Signal $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ yield with $m_a = 20\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$, with a branching ratio set to 1 for the $H \rightarrow aa$ decay, whereas the ATLAS figure attached to this entry instead uses the upper-limit branching ratio (smaller than 1). The table includes the yields in two signal regions with leptons consistent with an on-shell Z-boson decay, one with 2 $a\to b\bar{b}$ candidates in the 2HPC region and one with 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region. The table also includes the yields in four control regions, one with leptons consistent with an on-shell Z-boson decay and 2 $a\to b\bar{b}$ candidates in the Low Purity Category (LPC), and three others where the leptons are not consistent an on-shell Z-boson decay.

Observed $m_{t\bar{t}}^{reco}$ distributions in data for SR$2b$, shown together with the result of the fit with the three-shape-parameter function. The error bars indicate the effect of the fit parameter uncertainty on the background prediction. The bin width of the distributions is chosen to be the same as that used in the background parameterization.

Observed and expected upper limits on the cross-section times branching fraction of the $Z'_{TC2} \rightarrow t\bar{t}$ as a function of the $Z'_{TC2}$ mass ($\Gamma/m = 1.2$%).

Fit parameter values obtained from fits to data in the two signal regions. The fit function is $f(x) = p_{0} (1-x)^{p_{1}} x^{p_{2} + p_{3} \log(x)}$, where $x = m^{reco}_{t\bar{t}}/\sqrt{s}$.

Selection acceptance times efficiency for the 1 lepton signal events from MC simulations as a function of the resonance mass for VBF production.

Selection acceptance times efficiency for the 2 leptons signal events from MC simulations as a function of the resonance mass for ggF/DY production.

Selection acceptance times efficiency for the 2 leptons signal events from MC simulations as a function of the resonance mass for VBF production.

The fractions of signal events passing the VBF requirement on the RNN score as functions of the resonance mass for both VBF and ggF production. Information on the RNN model, such as the architecture, network weights and features scaling files are provided as additional resources.

Upper 95% CLs limits on the ggF production cross section of Radions in their $WW+ZZ$ decays as a function of the resonance mass.

Upper 95% CLs limits on the VBF production cross section of Radions in their $WW+ZZ$ decays as a function of the resonance mass.

Upper 95% CLs limits on the DY production cross section of HVT $W'$ in their $WZ$ decays as a function of the resonance mass.

Upper 95% CLs limits on the VBF production cross section of HVT $W'$ in their $WZ$ decays as a function of the resonance mass.

Upper 95% CLs limits on the DY production cross section of HVT $Z'$ in their $WW$ decays as a function of the resonance mass.

Upper 95% CLs limits on the VBF production cross section of HVT $Z'$ in their $WW$ decays as a function of the resonance mass.

Upper 95% CLs limits on the ggF production cross section of Gravitons in their $WW+ZZ$ decays as a function of the resonance mass.

Upper 95% CLs limits on the VBF production cross section of Gravitons in their $WW+ZZ$ decays as a function of the resonance mass.

Upper 95% CLs limits on the ggF production cross section of Radions in their $WW$ decays as a function of the resonance mass.

Upper 95% CLs limits on the VBF production cross section of Radions in their $WW$ decays as a function of the resonance mass.

Upper 95% CLs limits on the ggF production cross section of Radions in their $ZZ$ decays as a function of the resonance mass.

Upper 95% CLs limits on the VBF production cross section of Radions in their $ZZ$ decays as a function of the resonance mass.

Upper 95% CLs limits on the ggF production cross section of Gravitons in their $WW$ decays as a function of the resonance mass.

Upper 95% CLs limits on the VBF production cross section of Gravitons in their $WW$ decays as a function of the resonance mass.

Upper 95% CLs limits on the ggF production cross section of Gravitons in their $ZZ$ decays as a function of the resonance mass.

Upper 95% CLs limits on the VBF production cross section of Gravitons in their $ZZ$ decays as a function of the resonance mass.

Correlation matrix between the fiducial cross sections for the four individual decay final states and the $ZZ^{*}$ normalisation factor.

Differential fiducial cross section for the transverse momentum $p_{T}^{4l}$ of the Higgs boson. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 . Measured value in the last bin is un upper limit at 95% CL.

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the transverse momentum $p_{T}^{4l}$ of the Higgs boson.

Differential fiducial cross section for the invariant mass $m_{12}$ of the leading Z boson. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the invariant mass $m_{12}$ of the leading Z boson.

Differential fiducial cross section for the invariant mass $m_{34}$ of the subleading Z boson. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the invariant mass $m_{34}$ of the subleading Z boson.

Differential fiducial cross section for the rapidity $|y_{4l}|$ of the Higgs boson. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the rapidity $|y_{4l}|$ of the Higgs boson.

Differential fiducial cross section for the production angle $|\cos\theta^{*}|$ of the leading Z boson. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the production angle $|\cos\theta^{*}|$ of the leading Z boson.

Differential fiducial cross section for the production angle $\cos\theta_{1}$ of the anti-lepton from the leading Z boson. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the production angle $\cos\theta_{1}$ of the anti-lepton from the leading Z boson.

Differential fiducial cross section for the production angle $\cos\theta_{2}$ of the anti-lepton from the subleading Z boson. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the production angle $\cos\theta_{2}$ of the anti-lepton from the subleading Z boson.

Differential fiducial cross section for the azimuthal angle $\phi$ of the decay planes of the two reconstructed Z bosons. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the azimuthal angle $\phi$ of the decay planes of the two reconstructed Z bosons.

Differential fiducial cross section for the azimuthal angle $\phi_{1}$ of the decay plane of the leading Z boson and the plane formed between its four-momentum and the z-axis. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the azimuthal angle $\phi_{1}$ of the decay plane of the leading Z boson and the plane formed between its four-momentum and the z-axis.

Differential fiducial cross section for the jet multiplicity $N_{jets}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the jet multiplicity $N_{jets}$.

Differential fiducial cross section for the inclusive jet multiplicity $N_{jets}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Differential fiducial cross section for the number of b-quark initiated jets $N_{b-jets}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the number of b-quark initiated jets $N_{b-jets}$.

Differential fiducial cross section for the transverse momentum of the leading jet $p_{T}^{lead.jet}$ in events with at least one jet. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the transverse momentum of the leading jet $p_{T}^{lead.jet}$ in events with at least one jet.

Differential fiducial cross section for the transverse momentum of the subleading jet $p_{T}^{sublead.jet}$ in events with at least two jets. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the transverse momentum of the subleading jet $p_{T}^{sublead.jet}$ in events with at least two jets.

Differential fiducial cross section for the invariant mass of the two highest-pT jets $m_{jj}$ in events with at least two jets. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the invariant mass of the two highest-pT jets $m_{jj}$ in events with at least two jets.

Differential fiducial cross section for the distance between the two highest-pT jets in pseudorapidity $\Delta\eta_{jj}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the distance between the two highest-pT jets in pseudorapidity $\Delta\eta_{jj}$.

Differential fiducial cross section for the distance between the two highest-pT jets in $\phi$ $\Delta\phi_{jj}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the distance between the two highest-pT jets in $\phi$ $\Delta\phi_{jj}$.

Differential fiducial cross section for the transverse momentum of the four lepton plus jet system, in events with at least one jet $p_{T}^{4lj}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the transverse momentum of the four lepton plus jet system, in events with at least one jet $p_{T}^{4lj}$.

Differential fiducial cross section for the transverse momentum of the four lepton plus di-jet system, in events with at least two jets $p_{T}^{4ljj}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 . Measured value in the last bin is un upper limit at 95% CL.

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the transverse momentum of the four lepton plus di-jet system, in events with at least two jets $p_{T}^{4ljj}$.

Differential fiducial cross section for the invariant mass of the four lepton plus jet system in events with at least one jet $m_{4lj}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the invariant mass of the four lepton plus jet system in events with at least one jet $m_{4lj}$.

Differential fiducial cross section for the invariant mass of the four lepton plus di-jet system in events with at least two jets $m_{4ljj}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the invariant mass of the four lepton plus di-jet system in events with at least two jets $m_{4ljj}$.

Differential fiducial cross section for the leading vs. subleading Z boson mass $m_{12}$vs.$m_{34}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the leading vs. subleading Z boson mass $m_{12}$vs.$m_{34}$.

Differential fiducial cross section for the leading vs. subleading Z boson mass $m_{12}$vs.$m_{34}$ in $ll\mu\mu$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Differential fiducial cross section for the leading vs. subleading Z boson mass $m_{12}$vs.$m_{34}$ in $llee$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the leading vs. subleading Z boson mass m12 vs. m34 in $ll\mu\mu$ and $llee$ final states.

Differential fiducial cross section of the $p_{T}^{4l}$ distribution in $|y_{4l}|$ bins. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section of the $p_{T}^{4l}$ distribution in $|y_{4l}|$ bins.

Differential fiducial cross section of the $p_{T}^{4l}$ distribution in $N_{jets}$ bins. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section of the $p_{T}^{4l}$ distribution in $N_{jets}$ bins.

Differential fiducial cross section for transverse momentum of the four lepton system vs. the transverse momentum of the four lepton plus jet system $p_{T}^{4l}$vs.$p_{T}^{4lj}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for transverse momentum of the four lepton system vs. the transverse momentum of the four lepton plus jet system $p_{T}^{4l}$vs.$p_{T}^{4lj}$.

Differential fiducial cross section for the transverse momentum of the four lepton plus jet system vs the invariant mass of the four lepton plus jet system $p_{T}^{4l}$vs.$m_{4lj}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the transverse momentum of the four lepton plus jet system vs the invariant mass of the four lepton plus jet system $p_{T}^{4l}$vs.$m_{4lj}$.

Differential fiducial cross section for the transverse momentum of the four lepton vs the transverse momentum of the leading jet $p_{T}^{4l}$vs.$p_{T}^{l.jet}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the transverse momentum of the four lepton vs the transverse momentum of the leading jet $p_{T}^{4l}$vs.$p_{T}^{lead.jet}$.

Differential fiducial cross section for the transverse momentum of the leading jet vs the rapidity of the leading jet $p_{T}^{lead.jet}$vs.$|y^{lead.jet}|$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the transverse momentum of the leading jet vs the rapidity of the leading jet $p_{T}^{lead.jet}$vs.$|y^{lead.jet}|$.

Differential fiducial cross section for the transverse momentum of the leading jet vs the transverse momentum of the subleading jet $p_{T}^{lead.jet}$vs.$p_{T}^{sublead.jet}$. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the transverse momentum of the leading jet vs the transverse momentum of the subleading jet $p_{T}^{lead.jet}$vs.$p_{T}^{sublead.jet}$.

Differential fiducial cross section for the leading Z boson mass $m_{12}$ in $4\mu$ and $4e$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Differential fiducial cross section for the leading Z boson mass $m_{12}$ in $2e2\mu$ and $2\mu2e$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the leading Z boson mass $m_{12}$ in $4l$ and $2l2l$ final states.

Differential fiducial cross section for the subleading Z boson mass $m_{34}$ in $4\mu$ and $4e$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Differential fiducial cross section for the subleading Z boson mass $m_{34}$ in $2e2\mu$ and $2\mu2e$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the subleading Z boson mass $m_{34}$ in $4l$ and $2l2l$ final states.

Differential fiducial cross section for the azimuthal angle $\phi$ of the decay planes of the two reconstructed Z bosons in $4\mu$ and $4e$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Differential fiducial cross section for the azimuthal angle $\phi$ of the decay planes of the two reconstructed Z bosons in $2e2\mu$ and $2\mu2e$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the azimuthal angle $\phi$ of the decay planes of the two reconstructed Z bosons in $4l$ and $2l2l$ final states.

Differential fiducial cross section for the leading vs. subleading Z boson mass $m_{12}$vs.$m_{34}$ in $4\mu$ and $4e$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Differential fiducial cross section for the leading vs. subleading Z boson mass $m_{12}$vs.$m_{34}$ in $2\mu2e$ and $2e2\mu$ final states. The measured cross sections are compared to predictions provided by NNLOPS + XH. NNLOPS is normalised to the N3LO total cross section with a K-factor = 1.1 .

Correlation matrix between the measured cross sections and the $ZZ^{*}$ background normalization corresponding to the differential fiducial cross section for the leading vs. subleading Z boson mass $m_{12}$vs.$m_{34}$ in $4l$ and $2l2l$ final states.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 0.33 < ln(R/#DeltaR) < 0.67.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 0.67 < ln(R/#DeltaR) < 1.00.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 1.00 < ln(R/#DeltaR) < 1.33.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 1.33 < ln(R/#DeltaR) < 1.67.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 1.67 < ln(R/#DeltaR) < 2.00.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 2.00 < ln(R/#DeltaR) < 2.33.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 2.33 < ln(R/#DeltaR) < 2.67.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 2.67 < ln(R/#DeltaR) < 3.00.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 3.00 < ln(R/#DeltaR) < 3.33.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 3.33 < ln(R/#DeltaR) < 3.67.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 3.67 < ln(R/#DeltaR) < 4.00.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single vertical slice of the Lund jet plane between 4.00 < ln(R/#DeltaR) < 4.33.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 0.69 < ln(1/z) < 0.97.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 0.97 < ln(1/z) < 1.25.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 1.25 < ln(1/z) < 1.52.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 1.52 < ln(1/z) < 1.80.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 1.80 < ln(1/z) < 2.08.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 2.08 < ln(1/z) < 2.36.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 2.36 < ln(1/z) < 2.63.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 2.63 < ln(1/z) < 2.91.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 2.91 < ln(1/z) < 3.19.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 3.19 < ln(1/z) < 3.47.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 3.47 < ln(1/z) < 3.74.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 3.74 < ln(1/z) < 4.02.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 4.02 < ln(1/z) < 4.30.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 4.30 < ln(1/z) < 4.57.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 4.57 < ln(1/z) < 4.85.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 4.85 < ln(1/z) < 5.13.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 5.13 < ln(1/z) < 5.41.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 5.41 < ln(1/z) < 5.68.

Normalized differential cross-section of the Lund jet plane. The first systematic uncertainty is detector systematics, the second is background systematic uncertainties. The data is presented as a 1D distribution, for a single horizontal slice of the Lund jet plane between 5.68 < ln(1/z) < 5.96.

The summed covariance matrix of all systematic and statistical uncertainties associated with the measurement in bins of $\ln{(1/z)} \times \ln{(R/\Delta R)}$.

The summed covariance matrix of all statistical uncertainties associated with the measurement in bins of $\ln{(1/z)} \times \ln{(R/\Delta R)}$.

The expected (pre-fit) and observed numbers of events for an integrated luminosity of 139 fb$^{-1}$ at $\sqrt{s}$=13 TeV in the signal region 115 < $m_{4l}$< 130 GeV and sideband region $m_{4l}$ in 105-115 GeV or 130-160 GeV (350 GeV for $tXX$-enriched) in each reconstructed event category assuming the SM Higgs boson signal with a mass $m_{H}$= 125 GeV. Combined statistical and systematic uncertainties are included for the predictions. Expected contributions that are below 0.2% of the total yield in each reconstructed event category are not shown and replaced by -.

The expected and the observed (post-fit) the jet multiplicity distribution after the inclusive event selection for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed (post-fit) the four-lepton transverse momentum distribution for events with zero jets, after the inclusive event selection for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed (post-fit) the four-lepton transverse momentum distribution for events with one jet, after the inclusive event selection for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed (post-fit) the four-lepton transverse momentum distribution for events with at least two jets, after the inclusive event selection for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed (post-fit) the dijet invariant mass distribution for events with at least two jets, after the inclusive event selection for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{ggF}$ output (post-fit) distribution in 0$j$-$p_{T}^{4l}$-Low category for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{ggF}$ output (post-fit) distribution in 0$j$-$p_{T}^{4l}$-Med category for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{VBF}$ output (post-fit) distribution in 1$j$-$p_{T}^{4l}$-Low category with $NN_{ZZ}$<0.25 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{ZZ}$ output (post-fit) distribution in 1$j$-$p_{T}^{4l}$-Low category with $NN_{ZZ}$>0.25 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{VBF}$ output (post-fit) distribution in 1$j$-$p_{T}^{4l}$-Med category with $NN_{ZZ}$<0.25 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{ZZ}$ output (post-fit) distribution in 1$j$-$p_{T}^{4l}$-Med category with $NN_{ZZ}$>0.25 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{VBF}$ output (post-fit) distribution in 1$j$-$p_{T}^{4l}$-High category for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{VBF}$ output (post-fit) distribution in 2$j$ category with $NN_{VH}$<0.2 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{VH}$ output (post-fit) distribution in 2$j$ category with $NN_{VH}$>0.2 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{VBF}$ output (post-fit) distribution in 2$j$-BSM-like category for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{ttH}$ output (post-fit) distribution in $ttH$-Had-enriched category with $NN_{tXX}$<0.4 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{tXX}$ output (post-fit) distribution in $ttH$-Had-enriched category with $NN_{tXX}$>0.4 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{ttH}$ output (post-fit) distribution in $VH$-Lep-enriched category for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed events in the categories where no NN discriminant is used for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed events in the side bands used to constraint the $ZZ^{*}$ and $tXX$ backgrounds for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected SM cross section $(\sigma \cdot BR)_{SM}$, the observed value of $(\sigma \cdot BR)$, and their ratio $(\sigma \cdot BR)/(\sigma \cdot BR)_{SM}$ for the inclusive production and for each Production Mode and Reduced Stage-1.1 production bin for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13 TeV.

The correlation matrices between the measured cross-sections and the $ZZ$ and $tXX$ normalisation factors for Production Mode Stage.

The correlation matrices between the measured cross-sections and the $ZZ$ and $tXX$ normalisation factors for Reduced Stage 1.1.

Likelihood contours at 68% CLs in the (ggF,VBF) plane. The observed best fit value is the first value of the table.

Likelihood contours at 95% CLs in the (ggF,VBF) plane. The observed best fit value is the first value of the table.

Likelihood contours at 68% CLs in the (ggF,VH) plane. The observed best fit value is the first value of the table.

Likelihood contours at 95% CLs in the (ggF,VH) plane. The observed best fit value is the first value of the table.

Likelihood contours at 68% CLs in the (VBF,VH) plane. The observed best fit value is the first value of the table.

Likelihood contours at 95% CLs in the (VBF,VH) plane. The observed best fit value is the first value of the table.

Likelihood contours at 68% CLs in the ($gg2H-0j-p_T^H-Low$,$gg2H-0j-p_T^H-High$) plane. The observed best fit value is the first value of the table.

Likelihood contours at 95% CLs in the ($gg2H-0j-p_T^H-Low$,$gg2H-0j-p_T^H-High$) plane. The observed best fit value is the first value of the table.

Likelihood contours at 68% CLs in the ($\kappa_{V}$vs.$\kappa_{F}$) plane. The observed best fit value is the first value of the table.

Likelihood contours at 95% CLs in the ($\kappa_{V}$vs.$\kappa_{F}$) plane. The observed best fit value is the first value of the table.

The expected signal yield ratio for chosen CP-even and CP-odd EFT parameter values together with the corresponding cross-section measurement in each production bin of Reduced Stage 1.1. The parameter values correspond approximately to the expected confidence intervals at the 68% CLs obtained from the statistical interpretation of data.

The expected and observed confidence intervals at 68% and 95% CLs of the SMEFT Wilson coefficients for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13 TeV. For the coupling with two best fit values, the two values have been splitted in the two different columns 'Two Best-fit value - first' and 'Two Best-fit value - second', and the symbol '-' has been inserted in the missing fields.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{HW}$ versus $c_{HB}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{HW}$ versus $c_{HB}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{HB}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{HB}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{uH}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{uH}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

The best-fit values and the corresponding deviation from the SM prediction obtained from two-dimensional likelihood scans of the CP-odd BSM coupling parameters with 139fb$^{-1}$ data at $\sqrt{s}$=13 TeV. For the coupling with two best fit values, the two values have been splitted in the two different columns 'Two Observed best-fit first/second_POI - first' and 'Two Observed best-fit first/second_POI - second', and the symbol '-' has been inserted in the missing fields.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}}$ versus $c_{H\tilde{B}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}}$ versus $c_{H\tilde{B}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{B}}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{B}}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{\tilde{u}H}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{\tilde{u}H}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Signal acceptance times efficiency, expressed as a percentage, obtained from the ratio of the number of simulated signal events satisfying the event selection criteria in each reconstructed category to the total number of events generated in the phase space specified by a given Production model bin.

Signal acceptance times efficiency, expressed as a percentage, obtained from the ratio of the number of simulated signal events satisfying the event selection criteria in each reconstructed category to the total number of events generated in the phase space specified by a given Reduced Stage-1.1 production bin.

Observed covariance matrix for the different parameters of interest when fitting the for Production Mode Stage.

Observed covariance matrix for the different parameters of interest when fitting the for Reduced Stage 1.1.

The expected SM cross section $(\sigma\cdot BR)_{SM}$, the observed value of $(\sigma\cdot BR)$, and their ratio $(\sigma\cdot BR)/(\sigma\cdot BR)_{SM}$ for the alternate Reduced Stage-1.1 scheme for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13 TeV.

The correlation matrices between the measured cross-sections and the ZZ and $tXX$ normalisation factors for alternative Reduced Stage 1.1.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{HW}$ versus $c_{HWB}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{HW}$ versus $c_{HWB}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{HB}$ versus $c_{HWB}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{HB}$ versus $c_{HWB}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{HW}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{HW}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{HWB}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{HWB}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}}$ versus $c_{H\tilde{W}B}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}}$ versus $c_{H\tilde{W}B}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{B}}$ versus $c_{H\tilde{W}B}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{B}}$ versus $c_{H\tilde{W}B}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}B}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}B}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected and observed 95% CL upper limits on $\sigma(pp\to H)B(H\to Z(\eta_c~\text{or}~J/\psi))/$pb. The smaller (larger) quoted ranges around the expected limits represent $\pm 1\sigma$ ($\pm 2\sigma$) fluctuations.

Track multiplicity $n_{Tracks}$ for preselected DVs in muon-triggered events with at least one muon passing the full selection. Along with the data shown with black markers, the stacked filled histograms represent the background estimates, and predictions for signal scenarios are overlaid with dashed lines. The errors include statistical and systematic uncertainties and are indicated by hatched bands. The DV full selection requirements, $n_{Tracks} \geq 3$ and $m_{DV} > 20$ GeV are visualized with a black arrow.

Invariant mass $m_{DV}$ for the highest-mass preselected DV with at least three associated tracks in MET-triggered events with at least one muon passing the full selection. Along with the data shown with black markers, the stacked filled histograms represent the background estimates, and predictions for signal scenarios are overlaid with dashed lines. The errors include statistical and systematic uncertainties and are indicated by hatched bands. The DV full selection requirements, $n_{Tracks} \geq 3$ and $m_{DV} > 20$ GeV are visualized with a black arrow.

Invariant mass $m_{DV}$ for the highest-mass preselected DV with at least three associated tracks in muon-triggered events with at least one muon passing the full selection. Along with the data shown with black markers, the stacked filled histograms represent the background estimates, and predictions for signal scenarios are overlaid with dashed lines. The errors include statistical and systematic uncertainties and are indicated by hatched bands. The DV full selection requirements, $n_{Tracks} \geq 3$ and $m_{DV} > 20$ GeV are visualized with a black arrow.

The observed event yields in the control, validation and signal regions are shown for the MET Trigger selections, along with the predicted background yields. The bottom panel shows the ratio of observed events to the total background yields. The errors represent the total uncertainty of the backgrounds prediction, including the statistical and systematic uncertainties added in quadrature.

The observed event yields in the control, validation and signal regions are shown for the Muon Trigger selections, along with the predicted background yields. The bottom panel shows the ratio of observed events to the total background yields. The errors represent the total uncertainty of the backgrounds prediction, including the statistical and systematic uncertainties added in quadrature.

Expected exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Expected (1 sigma band) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Expected (2 sigma band) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Observed exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Observed (+1 sigma) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Observed (-1 sigma) exclusion limits at 95% CL on m($\tilde{t}$) as a function of $\tau(\tilde{t})$.

Exclusion limits on the production cross section as a function of m($\tilde{t}$) are shown for several values of $\tau(\tilde{t})$ along with the nominal signal production cross section and its theoretical uncertainty.

Parameterized event selection efficiencies for the $E_{T}^{miss}$ Trigger SR. The event-level efficiencies for each SR are extracted for all events passing the acceptance of the corresponding SR.

Parameterized event selection efficiencies for the Muon Trigger SR. The event-level efficiencies for each SR are extracted for all events passing the acceptance of the corresponding SR.

Parameterized muon-level reconstruction efficiencies as a function of the muon $p_{T}$ and $d_{0}$. The muon-level efficiencies are extracted using muons passing the muon acceptance criteria.

Parameterized vertex-level reconstruction efficiencies as a function of the radial position of the truth vertex. The efficiency is calculated independent of the muons originating from this truth vertex.

Parameterized vertex-level reconstruction efficiencies as a function of the radial position of the truth vertex. The efficiency is calculated only for truth vertices which have a muon originating from them which is matched to a reconstructed muon.

The $p_{T}$ distribution of all muons originating from LLP decays in the samples used to calculate and validate the efficiencies.

The invariant mass and multiplicity of selected decay products of all truth vertices used in the calculation and validation of the reconstructed efficiencies.

Observed and predicted mTtot distribution in the b-tag category of the 2tau_h channel. Please note that the bin content is divided by the bin width in the paper figure, but not in the HepData table. The last bin includes overflows. The combined prediction for A and H bosons with masses of 400, 1000 and 1500 GeV and $\tan\beta$ = 6, 12 and 25 respectively in the mh125 scenario are also provided. The combined prediction for A and H bosons with masses of 1000 and 1500 GeV is scaled by 100 in the paper figure, but not in the HepData table.

Observed and expected 95% CL upper limits on the gluon-gluon fusion Higgs boson production cross section times ditau branching fraction as a function of the Higgs boson mass.

Observed and expected 95% CL upper limits on the b-associated Higgs boson production cross section times ditau branching fraction as a function of the boson mass.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the $M_{h}^{125}$ scenario. The lowest value of $\tan\beta$ considered for the $M_{h}^{125}$ scenario is 0.5. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The observed 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus one sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus one sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with plus two sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

The expected 95% CL upper limits with minus two sigma on $\tan\beta$ as a function of $m_{A}$ in the hMSSM scenario. The lowest value of $\tan\beta$ considered by the hMSSM scenario is 0.8 and the highest value of mass is 2 TeV. The points in the region which is called "Not applicable" in the paper figure are kept in the HEPData table. Linear connection is applied in the range of signal mass points from 400 to 1000 GeV in the paper figure. The range of $\tan\beta$ shown in the paper figure and the HEPData is from 1 to 60.

Acceptance times efficiency for a scalar boson produced by gluon-gluon fusion as a function of the scalar boson mass.

Acceptance times efficiency for a scalar boson produced by b-associated production as a function of the scalar boson mass.

Observed 95% CL upper limits on the scalar boson production cross section times ditau branching fraction as a function of the scalar boson mass and the fraction of the b-associated production. The limits are calculated from a statistical combination of the 1l1tau_h and 2tau_h channels.

Expected 95% CL upper limits on the scalar boson production cross section times ditau branching fraction as a function of the scalar boson mass and the fraction of the b-associated production. The limits are calculated from a statistical combination of the 1l1tau_h and 2tau_h channels.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 200 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 250 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 300 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 350 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 400 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 500 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 600 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 700 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 800 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1000 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1200 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1500 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2000 GeV signal mass point is shown in the HEPData table.

Two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2500 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 200 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 250 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 300 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 350 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 400 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 500 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 600 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 700 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 800 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1000 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1200 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 1500 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2000 GeV signal mass point is shown in the HEPData table.

Expected two dimensional likelihood scan of the gluon-gluon fusion cross section times branching fraction, $\sigma(gg\phi)\times B(\phi\to\tau\tau)$, vs the b-associated production times branching fraction, $\sigma(bb\phi)\times B(\phi\to\tau\tau)$ for the scalar boson mass ($m_\phi$) indicated in the table. For each mass, 10000 points are scanned. At each point $\Delta(\mathrm{NLL})$ is calculated, defined as the negative-log-likelihood (NLL) of the conditional fit with $\sigma(gg\phi)$ and $\sigma(bb\phi)$ fixed to their values at the point and with the minimum NLL value at any point subtracted. The best-fit point and the preferred 68% and 95% boundaries are found at $2\Delta(\mathrm{NLL})$ values of 0.0, 2.30 and 5.90, respectively. The value of $2\Delta(\mathrm{NLL})$ for 2500 GeV signal mass point is shown in the HEPData table.

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 limits at 95% CL on the cross-sections of narrow resonance $X$ to $HH$ as a function of the mass $m_{X}$.

Observed and expected limits at 95% CL on the cross-sections of broad resonance $X$ to $HH$ as a function of the mass $m_{X}$.

Observed and expected limits at 95% CL on the cross-sections of non-resonant $HH$ production via VBF as a function of the di-vector-boson-di-Higgs-boson coupling modifier $\kappa_{2V}$. The theory prediction of the cross-section as a function of $\kappa_{2V}$ is also shown. The minimum theory cross-section is 1.37 fb at $\kappa_{2V}=1.17$.

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