$R_\text{CP}$ plotted as a function of approximated $x_p$ for $-1.0 < y_b < 0.0$ and $0.0 < y^* < 1.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_p$ for $-1.0 < y_b < 0.0$ and $1.0 < y^* < 2.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_p$ for $-1.0 < y_b < 0.0$ and $2.0 < y^* < 3.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_p$ for $0.0 < y_b < 1.0$ and $0.0 < y^* < 1.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_p$ for $0.0 < y_b < 1.0$ and $1.0 < y^* < 2.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_p$ for $0.0 < y_b < 1.0$ and $2.0 < y^* < 4.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_p$ for $1.0 < y_b < 2.0$ and $0.0 < y^* < 1.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_p$ for $1.0 < y_b < 2.0$ and $1.0 < y^* < 2.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_p$ for $1.0 < y_b < 2.0$ and $2.0 < y^* < 3.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_p$ for $2.0 < y_b < 3.0$ and $0.0 < y^* < 1.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_p$ for $2.0 < y_b < 3.0$ and $1.0 < y^* < 2.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_p$ for $3.0 < y_b < 4.0$ and $0.0 < y^* < 1.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_{Pb}$ for $-3.0 < y_b < -2.0$ and $0.0 < y^* < 1.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_{Pb}$ for $-2.0 < y_b < -1.0$ and $0.0 < y^* < 1.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_{Pb}$ for $-2.0 < y_b < -1.0$ and $1.0 < y^* < 2.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_{Pb}$ for $-1.0 < y_b < 0.0$ and $0.0 < y^* < 1.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_{Pb}$ for $-1.0 < y_b < 0.0$ and $1.0 < y^* < 2.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_{Pb}$ for $-1.0 < y_b < 0.0$ and $2.0 < y^* < 3.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_{Pb}$ for $0.0 < y_b < 1.0$ and $0.0 < y^* < 1.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_{Pb}$ for $0.0 < y_b < 1.0$ and $1.0 < y^* < 2.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_{Pb}$ for $0.0 < y_b < 1.0$ and $2.0 < y^* < 4.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_{Pb}$ for $1.0 < y_b < 2.0$ and $0.0 < y^* < 1.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_{Pb}$ for $1.0 < y_b < 2.0$ and $1.0 < y^* < 2.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_{Pb}$ for $1.0 < y_b < 2.0$ and $2.0 < y^* < 3.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_{Pb}$ for $2.0 < y_b < 3.0$ and $0.0 < y^* < 1.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_{Pb}$ for $2.0 < y_b < 3.0$ and $1.0 < y^* < 2.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_{Pb}$ for $3.0 < y_b < 4.0$ and $0.0 < y^* < 1.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_\mathrm{F}$ for $-3.0 < y_b < -2.0$ and $0.0 < y^* < 1.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_\mathrm{F}$ for $-2.0 < y_b < -1.0$ and $0.0 < y^* < 1.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_\mathrm{F}$ for $-2.0 < y_b < -1.0$ and $1.0 < y^* < 2.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_\mathrm{F}$ for $-1.0 < y_b < 0.0$ and $0.0 < y^* < 1.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_\mathrm{F}$ for $-1.0 < y_b < 0.0$ and $1.0 < y^* < 2.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_\mathrm{F}$ for $-1.0 < y_b < 0.0$ and $2.0 < y^* < 3.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_\mathrm{F}$ for $0.0 < y_b < 1.0$ and $0.0 < y^* < 1.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_\mathrm{F}$ for $0.0 < y_b < 1.0$ and $1.0 < y^* < 2.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_\mathrm{F}$ for $0.0 < y_b < 1.0$ and $2.0 < y^* < 4.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_\mathrm{F}$ for $1.0 < y_b < 2.0$ and $0.0 < y^* < 1.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_\mathrm{F}$ for $1.0 < y_b < 2.0$ and $1.0 < y^* < 2.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_\mathrm{F}$ for $1.0 < y_b < 2.0$ and $2.0 < y^* < 3.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_\mathrm{F}$ for $2.0 < y_b < 3.0$ and $0.0 < y^* < 1.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_\mathrm{F}$ for $2.0 < y_b < 3.0$ and $1.0 < y^* < 2.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$R_\text{CP}$ plotted as a function of approximated $x_\mathrm{F}$ for $3.0 < y_b < 4.0$ and $0.0 < y^* < 1.0$, constructed using $\langle y_{\text{b}} \rangle$ and $\langle y^{*} \rangle$. The proton-going direction is defined by $y_{\text{b}} > 0$.

$\langle T_{\text{AB}} \rangle$ normalized per-event dijet yields in $0-10$% collisions as a function of $\langle y_{\text{b}} \rangle$ for $40.0 < p_{\text{T,Avg}} \text{ [GeV]} < 49.6$ and $0.0 < y^{*} < 1.0$.

$\langle T_{\text{AB}} \rangle$ normalized per-event dijet yields in $60-90$% collisions as a function of $\langle y_{\text{b}} \rangle$ for $40.0 < p_{\text{T,Avg}} \text{ [GeV]} < 49.6$ and $0.0 < y^{*} < 1.0$.

$\langle T_{\text{AB}} \rangle$ normalized per-event dijet yields in $0-10$% collisions as a function of $\langle y_{\text{b}} \rangle$ for $49.6 < p_{\text{T,Avg}} \text{ [GeV]} < 61.4$ and $0.0 < y^{*} < 1.0$.

$\langle T_{\text{AB}} \rangle$ normalized per-event dijet yields in $60-90$% collisions as a function of $\langle y_{\text{b}} \rangle$ for $49.6 < p_{\text{T,Avg}} \text{ [GeV]} < 61.4$ and $0.0 < y^{*} < 1.0$.

$\langle T_{\text{AB}} \rangle$ normalized per-event dijet yields in $0-10$% collisions as a function of $\langle y_{\text{b}} \rangle$ for $61.4 < p_{\text{T,Avg}} \text{ [GeV]} < 76.1$ and $0.0 < y^{*} < 1.0$.

$\langle T_{\text{AB}} \rangle$ normalized per-event dijet yields in $60-90$% collisions as a function of $\langle y_{\text{b}} \rangle$ for $61.4 < p_{\text{T,Avg}} \text{ [GeV]} < 76.1$ and $0.0 < y^{*} < 1.0$.

$\langle T_{\text{AB}} \rangle$ normalized per-event dijet yields in $0-10$% collisions as a function of $\langle y_{\text{b}} \rangle$ for $40.0 < p_{\text{T,Avg}} \text{ [GeV]} < 49.6$ and $1.0 < y^{*} < 2.0$.

$\langle T_{\text{AB}} \rangle$ normalized per-event dijet yields in $60-90$% collisions as a function of $\langle y_{\text{b}} \rangle$ for $40.0 < p_{\text{T,Avg}} \text{ [GeV]} < 49.6$ and $1.0 < y^{*} < 2.0$.

$\langle T_{\text{AB}} \rangle$ normalized per-event dijet yields in $0-10$% collisions as a function of $\langle y_{\text{b}} \rangle$ for $49.6 < p_{\text{T,Avg}} \text{ [GeV]} < 61.4$ and $1.0 < y^{*} < 2.0$.

$\langle T_{\text{AB}} \rangle$ normalized per-event dijet yields in $60-90$% collisions as a function of $\langle y_{\text{b}} \rangle$ for $49.6 < p_{\text{T,Avg}} \text{ [GeV]} < 61.4$ and $1.0 < y^{*} < 2.0$.

$\langle T_{\text{AB}} \rangle$ normalized per-event dijet yields in $0-10$% collisions as a function of $\langle y_{\text{b}} \rangle$ for $61.4 < p_{\text{T,Avg}} \text{ [GeV]} < 76.1$ and $1.0 < y^{*} < 2.0$.

$\langle T_{\text{AB}} \rangle$ normalized per-event dijet yields in $60-90$% collisions as a function of $\langle y_{\text{b}} \rangle$ for $61.4 < p_{\text{T,Avg}} \text{ [GeV]} < 76.1$ and $1.0 < y^{*} < 2.0$.

$\langle T_{\text{AB}} \rangle$ normalized per-event dijet yields in $0-10$% collisions as a function of $\langle y_{\text{b}} \rangle$ for $40.0 < p_{\text{T,Avg}} \text{ [GeV]} < 49.6$ and $2.0 < y^{*} < 4.0$.

$\langle T_{\text{AB}} \rangle$ normalized per-event dijet yields in $60-90$% collisions as a function of $\langle y_{\text{b}} \rangle$ for $40.0 < p_{\text{T,Avg}} \text{ [GeV]} < 49.6$ and $2.0 < y^{*} < 4.0$.

$\langle T_{\text{AB}} \rangle$ normalized per-event dijet yields in $0-10$% collisions as a function of $\langle y_{\text{b}} \rangle$ for $49.6 < p_{\text{T,Avg}} \text{ [GeV]} < 61.4$ and $2.0 < y^{*} < 4.0$.

$\langle T_{\text{AB}} \rangle$ normalized per-event dijet yields in $60-90$% collisions as a function of $\langle y_{\text{b}} \rangle$ for $49.6 < p_{\text{T,Avg}} \text{ [GeV]} < 61.4$ and $2.0 < y^{*} < 4.0$.

$\langle T_{\text{AB}} \rangle$ normalized per-event dijet yields in $0-10$% collisions as a function of $\langle y_{\text{b}} \rangle$ for $61.4 < p_{\text{T,Avg}} \text{ [GeV]} < 76.1$ and $2.0 < y^{*} < 4.0$.

$\langle T_{\text{AB}} \rangle$ normalized per-event dijet yields in $60-90$% collisions as a function of $\langle y_{\text{b}} \rangle$ for $61.4 < p_{\text{T,Avg}} \text{ [GeV]} < 76.1$ and $2.0 < y^{*} < 4.0$.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of the number of $b$-jets compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of the number of $l/c$-jets compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $H_{\text{T}}^{\text{had}}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $\Delta R_{\text{avg}}^{bb}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{1})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{1}^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{2}^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{3})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{1}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $m(b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $m(bb^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $p_{\text{T}}(bb^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $\Delta R(e\mu bb^{\text{top}}, b_{1}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets and at least one $l/c$-jet as a function of $\Delta R(e\mu bb^{\text{top}}, l/c\text{-jet}_{1})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets and at least one $l/c$-jet as a function of $p_{\text{T}}(l/c\text{-jet}_{1})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets and at least one $l/c$-jet as a function of $p_{\text{T}}(l/c\text{-jet}_{1}) - p_{\text{T}}(b_{1}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $m(bb^{\text{min}\Delta R})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(bb^{\text{min}\Delta R})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $m(bb^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(bb^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(b_{3})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(b_{1}^{\text{add}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $\Delta R(b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $m(e\mu bb^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(l/c\text{-jet}_{1})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $\Delta\eta_{\text{max}}^{jj}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $H_{\text{T}}^{\text{all}}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $m(e\mu b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(b_{1})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(b_{2})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(b_{1}^{\text{top}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least three $b$-jets as a function of $|\eta(b_{2}^{\text{top}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{1})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{1}^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{2}^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{3})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{4})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{1}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{2}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $m(b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $m(bb^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $p_{\text{T}}(bb^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $H_{\text{T}}^{\text{all}}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $m(e\mu b_{1}b_{2})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $m(e\mu bb^{\text{top}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $H_{\text{T}}^{\text{had}}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $\text{min}\Delta R(bb)$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $\Delta R(e\mu bb^{\text{top}}, b_{1}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $\Delta R_{\text{avg}}^{bb}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $\Delta\eta_{\text{max}}^{jj}$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of the number of $l/c$-jets compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets and at least one $l/c$-jet as a function of $p_{\text{T}}(l/c\text{-jet}_{1})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets and at least one $l/c$-jet as a function of $|\eta(l/c\text{-jet}_{1})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets and at least one $l/c$-jet as a function of $\Delta R(e\mu bb^{\text{top}}, l/c\text{-jet}_{1})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets and at least one $l/c$-jet as a function of $p_{\text{T}}(l/c\text{-jet}_{1}) - p_{\text{T}}(b_{1}^{\text{add}})$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{1})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{2})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{1}^{\text{top}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at leastfour $b$-jets as a function of $|\eta(b_{2}^{\text{top}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{3})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{4})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{1}^{\text{add}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

Data bootstraps post unfolding for the normalised differential cross-section in the phase space with at least four $b$-jets as a function of $|\eta(b_{2}^{\text{add}})|$ compared with predictions. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. The last bin contains the overflow.

The measured normalised differential cross-section as a function of $N_{b-\text{jets}}$ in the $e\mu+\geq2b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $H_{\text{T}}^{\text{had}}$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $H_{\text{T}}^{\text{all}}$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R_{\text{avg}}^{bb}$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta\eta_{\text{max}}^{jj}$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1}^{\text{top}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{2})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{2}^{\text{top}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{3})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1}^{\text{add}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{1})|$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{1}^{\text{top}})|$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{2})|$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{2}^{\text{top}})|$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{3})|$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{1}^{\text{add}})|$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(b_{1}b_{2})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1}b_{2})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(bb^{\text{top}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(bb^{\text{top}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(e\mu b_{1}b_{2})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(e\mu bb^{\text{top}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(b_{1}b_{2})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $N_{l/c-\text{jets}}$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(e\mu b_{1}b_{2},b_{3})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(e\mu bb^{\text{top}}, b_{1}^{\text{add}})$ in the $e\mu+\geq3b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(e\mu bb^{\text{top}},l/c-\text{jet})$ in the $e\mu+\geq3b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(l/c\text{-jet}_{1}) - p_{\text{T}}(b_{1}^{\text{add}})$ in the $e\mu+\geq3b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(l/c\text{-jet}_{1})|$ in the $e\mu+\geq3b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(l/c\text{-jet}_{1})$ in the $e\mu+\geq3b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $H_{\text{T}}^{\text{had}}$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $H_{\text{T}}^{\text{all}}$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R_{\text{avg}}^{bb}$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta\eta_{\text{max}}^{jj}$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1}^{\text{top}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{2})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{2}^{\text{top}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{3})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1}^{\text{add}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{4})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{2}^{\text{add}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{1})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{1}^{\text{top}})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{2})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{2}^{\text{top}})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{3})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{1}^{\text{add}})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{4})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(b_{2}^{\text{add}})|$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(b_{1}b_{2})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(b_{1}b_{2})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(bb^{\text{top}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(bb^{\text{top}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(e\mu b_{1}b_{2})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(e\mu bb^{\text{top}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(bb^{\text{min}\Delta R})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(bb^{\text{min}\Delta R})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $m(bb^{\text{add}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(bb^{\text{add}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\text{min}\Delta R(bb)$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(b_{1}b_{2})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $N_{l/c-\text{jets}}$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(e\mu b_{1}b_{2},b_{3})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(e\mu bb^{\text{top}}, b_{1}^{\text{add}})$ in the $e\mu+\geq4b$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $\Delta R(e\mu bb^{\text{top}}, l/c\text{-jet}_{1})$ in the $e\mu+\geq4b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(l/c\text{-jet}_{1}) - p_{\text{T}}(b_{1}^{\text{add}})$ in the $e\mu+\geq4b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $|\eta(l/c\text{-jet}_{1})|$ in the $e\mu+\geq4b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The measured normalised differential cross-section as a function of $p_{\text{T}}(l/c\text{-jet}_{1})$ in the $e\mu+\geq4b+\geq1l/c-\text{jet}$ phase space. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $N_{b-\text{jets}}$ in the phase space with at least two b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $N_{b-\text{jets}}$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $H_{\text{T}}^{\text{had}}$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $H_{\text{T}}^{\text{all}}$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R_{\text{avg}}^{bb}$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta\eta_{\text{max}}^{jj}$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1}^{\text{top}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{2})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{2}^{\text{top}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{3})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1}^{\text{add}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{1})|$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{1}^{\text{top}})|$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{2}^{\text{top}})|$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{2}^{\text{top}})|$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{3})|$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{1}^{\text{add}})|$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(b_{1}b_{2})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1}b_{2})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(bb^{\text{top}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(bb^{\text{top}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(e\mu b_{1}b_{2})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(e\mu bb^{\text{top}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(b_{1}b_{2})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $N_{l/c-\text{jets}}$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(e\mu b_{1}b_{2},b_{3})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(e\mu bb^{\text{top}}, b_{1}^{\text{add}})$ in the phase space with at least three b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(e\mu bb^{\text{top}},l/c-\text{jet})$ in the phase space with at least three b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(l/c\text{-jet}_{1}) - p_{\text{T}}(b_{1}^{\text{add}})$ in the phase space with at least three b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(l/c\text{-jet}_{1})|$ in the phase space with at least three b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(l/c\text{-jet}_{1})$ in the phase space with at least three b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $H_{\text{T}}^{\text{had}}$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $H_{\text{T}}^{\text{all}}$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R_{\text{avg}}^{bb}$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta\eta_{\text{max}}^{jj}$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1}^{\text{top}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{2})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{2}^{\text{top}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{3})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1}^{\text{add}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{4})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{2}^{\text{add}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{1})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{1}^{\text{top}})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{2})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{2}^{\text{top}})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{3})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{1}^{\text{add}})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{4})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(b_{2}^{\text{add}})|$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(b_{1}b_{2})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(b_{1}b_{2})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(bb^{\text{top}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(bb^{\text{top}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(e\mu b_{1}b_{2})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(e\mu bb^{\text{top}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(bb^{\text{min}\Delta R})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(bb^{\text{min}\Delta R})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $m(bb^{\text{add}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(bb^{\text{add}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\text{min}\Delta R(bb)$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(b_{1}b_{2})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $N_{l/c-\text{jets}}$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(e\mu b_{1}b_{2},b_{3})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(e\mu bb^{\text{top}}, b_{1}^{\text{add}})$ in the phase space with at least four b-jets. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $\Delta R(e\mu bb^{\text{top}}, l/c\text{-jet}_{1})$ in the phase space with at least four b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(l/c\text{-jet}_{1}) - p_{\text{T}}(b_{1}^{\text{add}})$ in the phase space with at least four b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $|\eta(l/c\text{-jet}_{1})|$ in the phase space with at least four b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

The correlation matrix for the measured normalised differential cross-section in terms of $p_{\text{T}}(l/c\text{-jet}_{1})$ in the phase space with at least four b-jets and at least one $l/c$-jet. The overflow is included in the last bin.

Distribution of the BDT output score in the the 3l channel signal region.

Distribution of the BDT output score in the 2lSC channel signal region.

Distribution of the BDT output score in the 2lSC+$\tau_{had}$ channel signal region.

Distribution of the BDT output score in the the 2l+2$\tau_{had}$ channel signal region.

Distribution of the BDT output score in the l+2$\tau_{had}$ channel signal region.

Invariant mass of the diphoton system in the $\gamma\gamma$+2(l,$\tau_{had}$) channel signal region.

Invariant mass of the diphoton system in the $\gamma\gamma$+$e/\mu$ channel in the loose BDT signal region.

Invariant mass of the diphoton system in the $\gamma\gamma$+$e/\mu$ channel in the medium BDT signal region.

Invariant mass of the diphoton system in the $\gamma\gamma$+$e/\mu$ channel in the tight BDT signal region.

Invariant mass of the diphoton system in the $\gamma\gamma$+$\tau$ channel in the loose BDT signal region.

Invariant mass of the diphoton system in the $\gamma\gamma$+$\tau$ channel in the medium BDT signal region.

Invariant mass of the diphoton system in the $\gamma\gamma$+$\tau$ channel in the tight BDT signal region.

Observed and expected 95% CL upper limits on the signal strength for the ML channels and the $\gamma\gamma$+ML channels.

Expected values of -2ln$\Lambda$ as a function of $\kappa_\lambda$ for the combination. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Observed values of -2ln$\Lambda$ as a function of $\kappa_\lambda$ for the combination. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Expected values of -2ln$\Lambda$ as a function of $\kappa_{2V}$ for the combination. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Observed values of -2ln$\Lambda$ as a function of $\kappa_{2V}$ for the combination. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Expected values of -2ln$\Lambda$ as a function of $\kappa_\lambda$ for the ML channels. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Observed values of -2ln$\Lambda$ as a function of $\kappa_\lambda$ for the ML channels. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Expected values of -2ln$\Lambda$ as a function of $\kappa_\lambda$ for the $\gamma\gamma$+ML channels. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Observed values of -2ln$\Lambda$ as a function of $\kappa_\lambda$ for the $\gamma\gamma$+ML channels. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Expected values of -2ln$\Lambda$ as a function of $\kappa_{2V}$ for the ML channels. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Observed values of -2ln$\Lambda$ as a function of $\kappa_{2V}$ for the ML channels. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Expected values of -2ln$\Lambda$ as a function of $\kappa_{2V}$ for the $\gamma\gamma$+ML channels. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Observed values of -2ln$\Lambda$ as a function of $\kappa_{2V}$ for the $\gamma\gamma$+ML channels. All other coupling modifiers are fixed to their SM predictions and the expected limits are computed assuming the SM.

Absolute differential cross-section as a function of RC large-R jet $\tau_{21}$ at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of RC large-R jet $\tau_{21}$ at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative differential cross-section as a function of RC large-R jet $\tau_{21}$ at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of RC large-R jet $\tau_{3}$ at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of RC large-R jet $\tau_{3}$ at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative differential cross-section as a function of RC large-R jet $\tau_{3}$ at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of RC large-R jet ECF2 at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of RC large-R jet ECF2 at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative differential cross-section as a function of RC large-R jet ECF2 at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of RC large-R jet $D_{2}$ at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of RC large-R jet $D_{2}$ at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative differential cross-section as a function of RC large-R jet $D_{2}$ at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of RC large-R jet $C_{3}$ at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of RC large-R jet $C_{3}$ at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative differential cross-section as a function of RC large-R jet $C_{3}$ at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of RC large-R jet $p_\mathrm{T}^{d,*}$ (jet $p_T$ dispersion) at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of RC large-R jet $p_\mathrm{T}^{d,*}$ (jet $p_T$ dispersion) at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative differential cross-section as a function of RC large-R jet $p_\mathrm{T}^{d,*}$ (jet $p_T$ dispersion) at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute differential cross-section as a function of RC large-R jet LHA at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix of the Absolute differential cross-section as function of RC large-R jet LHA at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative differential cross-section as a function of RC large-R jet LHA at particle level in the $\ell$+jets channel. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 122.5 GeV < $m^{top}$ < 157.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 157.5 GeV < $m^{top}$ < 187.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 187.5 GeV < $m^{top}$ < 222.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 122.5 GeV < $m^{top}$ < 157.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 122.5 GeV < $m^{top}$ < 157.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 157.5 GeV < $m^{top}$ < 187.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 122.5 GeV < $m^{top}$ < 157.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 157.5 GeV < $m^{top}$ < 187.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 157.5 GeV < $m^{top}$ < 187.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 187.5 GeV < $m^{top}$ < 222.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 122.5 GeV < $m^{top}$ < 157.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 187.5 GeV < $m^{top}$ < 222.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 157.5 GeV < $m^{top}$ < 187.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 187.5 GeV < $m^{top}$ < 222.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $m^{top}$ in 187.5 GeV < $m^{top}$ < 222.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 122.5 GeV < $m^{top}$ < 157.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 157.5 GeV < $m^{top}$ < 187.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 187.5 GeV < $m^{top}$ < 222.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 122.5 GeV < $m^{top}$ < 157.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 157.5 GeV < $m^{top}$ < 187.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 187.5 GeV < $m^{top}$ < 222.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 122.5 GeV < $m^{top}$ < 157.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 122.5 GeV < $m^{top}$ < 157.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 157.5 GeV < $m^{top}$ < 187.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 122.5 GeV < $m^{top}$ < 157.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 157.5 GeV < $m^{top}$ < 187.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 157.5 GeV < $m^{top}$ < 187.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 187.5 GeV < $m^{top}$ < 222.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 122.5 GeV < $m^{top}$ < 157.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 187.5 GeV < $m^{top}$ < 222.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 157.5 GeV < $m^{top}$ < 187.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 187.5 GeV < $m^{top}$ < 222.5 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $m^{top}$ in 187.5 GeV < $m^{top}$ < 222.5 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 122.5 GeV < $m^{top}$ < 157.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 157.5 GeV < $m^{top}$ < 187.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $m^{top}$ at particle level in the $\ell$+jets channel in 187.5 GeV < $m^{top}$ < 222.5 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 350.0 GeV < $p_{T}$ < 500.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 500.0 GeV < $p_{T}$ < 550.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 550.0 GeV < $p_{T}$ < 650.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 650.0 GeV < $p_{T}$ < 750.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 750.0 GeV < $p_{T}$ < 2000.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $\tau_{32}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 350.0 GeV < $p_{T}$ < 500.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 500.0 GeV < $p_{T}$ < 550.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 550.0 GeV < $p_{T}$ < 650.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 650.0 GeV < $p_{T}$ < 750.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $\tau_{32}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 750.0 GeV < $p_{T}$ < 2000.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 350.0 GeV < $p_{T}$ < 500.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 500.0 GeV < $p_{T}$ < 550.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 550.0 GeV < $p_{T}$ < 650.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 650.0 GeV < $p_{T}$ < 750.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Absolute double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 750.0 GeV < $p_{T}$ < 2000.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 350.0 GeV < $p_{T}$ < 500.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 500.0 GeV < $p_{T}$ < 550.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 550.0 GeV < $p_{T}$ < 650.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 650.0 GeV < $p_{T}$ < 750.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Covariance matrix between the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV and the Absolute double-differential cross-section as function of RC large-R jet $D_{2}$ vs $p_{T}$ in 750.0 GeV < $p_{T}$ < 2000.0 GeV at particle level in the $\ell$+jets channel, accounting for the statistical uncertainty.

Relative double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 350.0 GeV < $p_{T}$ < 500.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 500.0 GeV < $p_{T}$ < 550.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 550.0 GeV < $p_{T}$ < 650.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 650.0 GeV < $p_{T}$ < 750.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Relative double-differential cross-section as a function of RC large-R jet $D_{2}$ vs $p_{T}$ at particle level in the $\ell$+jets channel in 750.0 GeV < $p_{T}$ < 2000.0 GeV. The measured differential cross-section is compared with the prediction obtained with the Powheg+Pythia8 Monte Carlo generator.

Local observed significance for signal discovery at different (mX, mS). The points show where the significance was evaluated. The band at mS = 125 GeV is not shown as those points are equivalent to those already probed in Phys. Rev. D 106 (2022) 052001.

Signal region yield for each of the simulated signal samples and interpolated signal points that were analyzed in the 1 b-jet Signal Region, normalized to a cross section of 1 fb. The signal region yields for one of the signals is found in Table 2 of the paper.

Signal region yield for each of the simulated signal samples and interpolated signal points that were analyzed in the 2 b-jet Signal Region, normalized to a cross section of 1 fb. The signal region yields for one of the signals is found in Table 2 of the paper.

Invariant mass of the lepton pair plus jet system, for data, background (pre- and post-reweighting) and three $H\rightarrow Za$ signal hypotheses (for $a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection but not the classification NN requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Regression NN output variable, for data, reweighted background, and three $H\rightarrow Za$ signal hypotheses ($a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection, including a $120<m_{\ell\ell j}<140$ GeV requirement, but not the classification NN output variable requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Classification NN output variable, for data, reweighted background, and three $H\rightarrow Za$ signal hypotheses ($a\rightarrow q\bar{q}/gg$ inclusively). Events are required to pass the complete event selection, including a $120<m_{\ell\ell j}<140$ GeV requirement, but not the classification NN output variable requirement. The background normalization is set equal to that of the data for events passing the preselection and being in the $m_{\ell\ell j}$ 100-180 GeV region. The signal normalization assumes the SM Higgs boson inclusive production cross-section, $\mathcal{B}(H\to Za)=100\%$, and it is scaled up by a factor of 100. The error bars (hatched regions) represent the data (MC) sample's statistical uncertainty in the histograms and the ratio plots. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Invariant mass of the lepton pair plus jet system for data, reweighted background, and various signal hypotheses. Events are required to pass the complete event selection (preselection and classification NN requirement). The pre-fit background distribution is shown along with three $H\rightarrow Za$ signal hypotheses coming from the pure MC sample ($a\rightarrow q\bar{q}/gg$ inclusively). The background normalization is set equal to that of the data for events passing the preselection and being in the region 100-180 GeV. The error bars represent the data sample's statistical uncertainty and the grey bands represent the full background model uncertainty before the fit (assuming all the uncertainties are uncorrelated). The lower panel shows the ratio of the data to the pre-fit background prediction.

Invariant mass of the lepton pair plus jet system for data, reweighted background, and 2.0 GeV signal hypotheses. Events are required to pass the complete event selection (preselection and classification NN requirement). The background shape and normalization come from the fit to the data. The signal includes only $gg$ decays, its shape comes from the fit, and it is normalized to $\mathcal{B}(H\rightarrow Za)=100\%$. The error bars represent the data sample's statistical uncertainty and the grey bands represent the full background model uncertainty after the fit. The lower panel shows the per-bin signal significance, $(N_\text{data}-N_\text{MC})\,\text{\Large{/}}\sqrt{\sigma_\text{data}^2+\sigma_\text{MC}^2}$. The post-fit background uncertainty is calculated from a fit under the background-only hypothesis. Vertical arrows indicate data points that fall outside the displayed $y$-axis range.

Observed 95$\%$ CL upper limits on $\mathcal{B}(H \to Za)=100\%$ for $\mathcal{B}(a \to gg)=100\%$, and the expectation under the background-only hypothesis together with its $\pm1\sigma$ and $\pm2\sigma$ intervals. The weaker limits from the previous version of the analysis (PRL 125(2020) 221802) are also shown. Linear interpolation is used to set limits between the mass points for which MC signal samples were generated.

Observed 95$\%$ confidence level upper limits on $\mathcal{B}(H \to Za)=100\%$ for $\mathcal{B}(a\rightarrow q\bar{q})=100\%$, and the expectation under the background-only hypothesis together with its $\pm1\sigma$ and $\pm2\sigma$ intervals. The weaker limits from the previous version of the analysis (PRL 125(2020) 221802) are also shown. Linear interpolation is used to set limits between the mass points for which MC signal samples were generated.

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{loose}}^{2\ell3j}$ signal region, which uses the invariant mass of the c-tagged jet and the geometrically closest b-tagged jet, $m_{\mathit{cb}}^{\mathrm{min}\Delta R}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{loose}}^{1\ell6ji}$ signal region, which uses the jet multiplicity, $N_{\mathrm{jets}}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{tight}}^{1\ell6ji}$ signal region, which uses the jet multiplicity, $N_{\mathrm{jets}}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{loose}}^{2\ell4ji}$ signal region, which uses the jet multiplicity, $N_{\mathrm{jets}}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Post-fit agreement between data and MC prediction for the $SR_{\mathrm{tight}}^{2\ell4ji}$ signal region, which uses the jet multiplicity, $N_{\mathrm{jets}}$, as an observable. The hatched uncertainty bands include all uncertainties and their correlations. The last bins contain overflow events. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes.

Jet multiplicity and tagging criteria for the SRs and CRs in the single-lepton and dilepton channels. All single-lepton regions are defined in the 5-jet-exclusive and 6-jet-inclusive jet selections, all dilepton regions in the 3-jet-exclusive and 4-jet-inclusive jet selections. By construction, $\mathrm{SR}^{2\ell}_{\mathrm{tight}}$ only exists for the 4-jet-inclusive selection. Identical requirements on inclusive and exclusive working points, e.g., one c@22% and one c@11% in the $\mathrm{CR}_1^{1\ell}$, indicate a veto of additional tags at the inclusive working point.

Breakdown of the fiducial $t\bar{t}+{\geq}2c$ and $t\bar{t}+1c$ cross-section uncertainties. Systematic uncertainties are grouped into signal and background modeling, instrumental, and MC statistics categories. The table lists the fractional uncertainty in the measured cross-sections in percent and is estimated from the covariance matrix of the profile likelihood fit, following EPJC 84 (2024) 593. Individual groups of uncertainties can be added in quadrature to obtain the total uncertainty. JES and JER denote the jet energy scale and resolution uncertainties, respectively. The fractional uncertainties in the fitted $t\bar{t}+{\geq}1c$ and $t\bar{t}+\text{light}$ normalization factors are included in the data statistics category. For presentation purposes, all uncertainties are symmetrized.

Measured fiducial cross-section values in comparison with various NLO+PS predictions. The uncertainties in the predictions include independent and simultaneous variations of the $\mu_R$ and $\mu_F$ scales and uncertainties in the choice of the PDF set, but no uncertainties in the parton shower or hadronization. The uncertainties in the measured values include all statistical and systematic uncertainties. For presentation purposes, the quoted measurement and prediction uncertainties are symmetrized.

Measured and predicted values for the $t\bar{t}+{\geq}1b$, $t\bar{t}+{\geq}2c$, and $t\bar{t}+1c$ cross-sections and for the cross-section ratios of these processes to total $t\bar{t}+\text{jets}$ production. The quoted uncertainties in the measurements include statistical and systematic uncertainties. The predictions are taken from the nominal setup using inclusive $t\bar{t}$ Powheg+Pythia8 predictions for the $t\bar{t}+{\geq}2c$, $t\bar{t}+1c$, and $t\bar{t}+\text{light}$ processes, and $t\bar{t}+b\bar{b}$ Powheg+Pythia8 predictions for the $t\bar{t}+{\geq}1b$ process. The first use the 5FS, the latter the 4FS, where the $b$-quarks are treated as massive particles and are included in the matrix element calculation. The uncertainties in the predictions include independent and simultaneous variations of the $\mu_R$ and $\mu_F$ scales and uncertainties in the choice of the PDF set.

Nuisance parameters with a large impact on the fiducial measurement: ranking of the ten most impactful nuisance parameters for the $t\bar{t}+{\geq}2c$ cross-section measurement. The colored boxes indicate the impact of the nuisance parameters, $\Delta\mu$ (top axis), on the corresponding parameter of interest, $\mu$. They are calculated from the covariance matrix using the post-fit up (down) uncertainties of the parameter of interest and the nuisance parameter, $\sigma_{\mathrm{POI}}^{\mathrm{up}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{up}}$ ($\sigma_{\mathrm{POI}}^{\mathrm{down}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{down}}$), and their linear correlation coefficient, $\rho$, following EPJC 84 (2024) 593. The data points and error bars indicate the post-fit values and uncertainties of the nuisance parameters relative to their pre-fit values, $\theta_0$, in units of their pre-fit uncertainties, $\Delta \theta$ (bottom axis). A hollow circle indicates that the nuisance parameter has a Poisson constraint.

Nuisance parameters with a large impact on the fiducial measurement: ranking of the ten most impactful nuisance parameters for the $t\bar{t}+1c$ cross-section measurement. The colored boxes indicate the impact of the nuisance parameters, $\Delta\mu$ (top axis), on the corresponding parameter of interest, $\mu$. They are calculated from the covariance matrix using the post-fit up (down) uncertainties of the parameter of interest and the nuisance parameter, $\sigma_{\mathrm{POI}}^{\mathrm{up}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{up}}$ ($\sigma_{\mathrm{POI}}^{\mathrm{down}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{down}}$), and their linear correlation coefficient, $\rho$, following EPJC 84 (2024) 593. The data points and error bars indicate the post-fit values and uncertainties of the nuisance parameters relative to their pre-fit values, $\theta_0$, in units of their pre-fit uncertainties, $\Delta \theta$ (bottom axis). A hollow circle indicates that the nuisance parameter has a Poisson constraint.

Nuisance parameters with a large impact on the fiducial measurement: ranking of the ten most impactful nuisance parameters for the $t\bar{t}+{\geq}2c$ ratio to total $t\bar{t}+\text{jets}$ production. The colored boxes indicate the impact of the nuisance parameters, $\Delta\mu$ (top axis), on the corresponding parameter of interest, $\mu$. They are calculated from the covariance matrix using the post-fit up (down) uncertainties of the parameter of interest and the nuisance parameter, $\sigma_{\mathrm{POI}}^{\mathrm{up}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{up}}$ ($\sigma_{\mathrm{POI}}^{\mathrm{down}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{down}}$), and their linear correlation coefficient, $\rho$, following EPJC 84 (2024) 593. The data points and error bars indicate the post-fit values and uncertainties of the nuisance parameters relative to their pre-fit values, $\theta_0$, in units of their pre-fit uncertainties, $\Delta \theta$ (bottom axis). A hollow circle indicates that the nuisance parameter has a Poisson constraint.

Nuisance parameters with a large impact on the fiducial measurement: ranking of the ten most impactful nuisance parameters for the $t\bar{t}+1c$ ratio to total $t\bar{t}+\text{jets}$ production. The colored boxes indicate the impact of the nuisance parameters, $\Delta\mu$ (top axis), on the corresponding parameter of interest, $\mu$. They are calculated from the covariance matrix using the post-fit up (down) uncertainties of the parameter of interest and the nuisance parameter, $\sigma_{\mathrm{POI}}^{\mathrm{up}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{up}}$ ($\sigma_{\mathrm{POI}}^{\mathrm{down}}$ and $\sigma_{\mathrm{NP}}^{\mathrm{down}}$), and their linear correlation coefficient, $\rho$, following EPJC 84 (2024) 593. The data points and error bars indicate the post-fit values and uncertainties of the nuisance parameters relative to their pre-fit values, $\theta_0$, in units of their pre-fit uncertainties, $\Delta \theta$ (bottom axis). A hollow circle indicates that the nuisance parameter has a Poisson constraint.

Efficiencies and rejection rates for the $b$- and $c$-jet tagging working points of the $b/c$-tagger, as well as the selection cuts on the discriminants, $\mathcal{D}_b'$ and $\mathcal{D}_c'$ that define the tagging bins. The primes indicate that a standard logistic function is applied to the discriminants. The b@70% and c@22% working points are inclusive of their respective tighter working points. The values were estimated from single-lepton and dilepton $t\bar{t}$ events simulated with Powheg + Pythia 8.

Efficiencies and rejection rates for the $b$- and $c$-jet tagging working points of the $b/c$-tagger, as well as the selection cuts on the discriminants, $\mathcal{D}_b'$ and $\mathcal{D}_c'$ that define the tagging bins. The primes indicate that a standard logistic function is applied to the discriminants. The b@70% and c@22% working points are inclusive of their respective tighter working points. The values were estimated from single-lepton and dilepton $t\bar{t}$ events simulated with Powheg + Pythia 8.

Pre-fit event yields in the twelve single-lepton and dilepton control regions (CRs) of the analysis. The quoted uncertainties include both statistical and systematic components, but all normalization factors are fixed to their pre-fit values with no uncertainties. All uncertainties are assumed to be uncorrelated. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes. The "Fakes" category includes events from misidentified leptons and non-prompt leptons, estimated through the data-driven matrix method in the single-lepton channel and through MC predictions in the dilepton channel.

Pre-fit event yields in the seven signal regions (SRs) of the analysis. The quoted uncertainties include both statistical and systematic components, but all normalization factors are fixed to their pre-fit values with no uncertainties. All uncertainties are assumed to be uncorrelated. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes. The "Fakes" category includes events from misidentified leptons and non-prompt leptons, estimated through the data-driven matrix method in the single-lepton channel and through MC predictions in the dilepton channel.

Post-fit event yields in the twelve single-lepton and dilepton control regions (CRs) of the analysis. The quoted uncertainties include both statistical and systematic components, and post-fit correlations between uncertainties are taken into account. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes. The "Fakes" category includes events from misidentified leptons and non-prompt leptons, estimated through the data-driven matrix method in the single-lepton channel and through MC predictions in the dilepton channel.

Post-fit event yields in the seven signal regions (SRs) of the analysis. The quoted uncertainties include both statistical and systematic components, and post-fit correlations between uncertainties are taken into account. "Other Top" includes single-top-quark production and associated production of $t\bar{t}$ and single top quarks with bosons. "Non-Top" includes W+jets, Z+jets, and diboson processes. The "Fakes" category includes events from misidentified leptons and non-prompt leptons, estimated through the data-driven matrix method in the single-lepton channel and through MC predictions in the dilepton channel.

The cutflow table for $N1$ and $N2$ signal samples for different mass points. The samples consist of both $ee$ and $\mu\mu$ final states in equal ratio. Each sample is normalized to a total cross-section of 0.1 pb. The abbreviations used in the table are, SC = Same Charge, SF = Same Flavor.

The cutflow table for $N3$ signal samples for different mass points. The samples consist of both $ee$ and $\mu\mu$ final states in equal ratio. Each sample is normalized to a total cross-section of 10 pb. The abbreviations used in the table are, SS = Same Charge, SF = Same Flavor.

Expected and observed upper limits on the strength of HNL mixing with electron neutrinos

Expected and observed upper limits on the strength of HNL mixing with muon neutrinos

Expected and observed upper limits on the strength of HNL mixing with tau neutrinos

Distributions of the $\mathrm{NN^{SvsB}}$ output for data and the expected background after the likelihood fit in the $SR_{cQu --}$ signal region. The post-fit background expectations are shown as filled histograms, the combined pre-fit background expectations are shown as dashed lines. The signal distribution using the Wilson coefficient values $c_{tu}^{(1)}=0.04$, $c_{Qu}^{(1)}=0.1$, $c_{Qu}^{(8)}=0.1$ is shown with a dotted line, normalized to the same number of events as the background.

Summary of the observed and predicted number of events in the five 2$\ell$ control regions. The background prediction is shown after the combined likelihood fit to data under the signal-plus-background hypothesis across all control regions and signal regions. The uncertainties in the total yields are smaller than the sum in quadrature of the uncertainties in the individual contributions due to anti-correlations resulting from the likelihood fit. Int Conv refers to events with photon conversions in the inner detector, while Mat Conv refers to events with photon conversions in the detector material. QMisID refers to events with charge-misidentified lepton.

Summary of the observed and predicted number of events in the four 3$\ell$ control regions. The background prediction is shown after the combined likelihood fit to data under the signal-plus-background hypothesis across all control regions and signal regions. The uncertainties in the total yields are smaller than the sum in quadrature of the uncertainties in the individual contributions due to anti-correlations resulting from the likelihood fit. Int Conv refers to events with photon conversions in the inner detector, while Mat Conv refers to events with photon conversions in the detector material. QMisID refers to events with charge-misidentified lepton.

Comparison of the observed limits on the three WCs to the limit obtained by the ATLAS Run~1 analysis of events with $b$-tagged jets and a same-sign lepton pair [JHEP10(2015)150]. The cross-section limits from that analysis are converted to the limits on the WCs by using the fitted EFT parametrization. The bounds on the WCs are shown at the 68$\%$ CL (solid) and/or 95$\%$ CL (dashed) levels. For the ATLAS Run~1 analysis only the bounds at 95$\%$ CL are available. The vertical bar represents the SM prediction. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$.

Observed lower limits at 95$\%$ confidence level on the scale of new physics $\Lambda$ for Wilson coefficient values of 0.01, 1 and $4\pi^2$. The limits are compared with the limits obtained by the ATLAS Run~1 analysis of events with $b$-tagged jets and a same-sign lepton pair.

Observed and expected limits at 95$\%$ confidence level on the cross-section of $\sigma(pp \rightarrow tt)$.

2D likelihood scans of the Wilson coefficients $c_{tu}^{(1)}$ versus $c_{Qu}^{(1)}$. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$. In the table the $2\Delta$NLL values are shown which correspond to the minimized negative log likelihood for a given Wilson coefficient pair.

2D likelihood scans of the Wilson coefficients $c_{tu}^{(1)}$ versus $c_{Qu}^{(8)}$. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$. In the table the $2\Delta$NLL values are shown which correspond to the minimized negative log likelihood for a given Wilson coefficient pair.

2D likelihood scans of the Wilson coefficients $c_{Qu}^{(1)}$ versus $c_{Qu}^{(8)}$. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$. In the table the $2\Delta$NLL values are shown which correspond to the minimized negative log likelihood for a given Wilson coefficient pair.

1D-likelihood scan for $c_{tu}^{(1)}$. The observed and expected limits at 68\% and 95\% confidence level (CL) can be read off by finding the intersection of the likelihood scan curve with the 1$\sigma$ and 2$\sigma$ lines. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$. In the table the $2\Delta$NLL values are shown which correspond to the minimized negative log likelihood for a given Wilson coefficient value.

1D-likelihood scan for $c_{Qu}^{(1)}$. The observed and expected limits at 68\% and 95\% confidence level (CL) can be read off by finding the intersection of the likelihood scan curve with the 1$\sigma$ and 2$\sigma$ lines. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$. In the table the $2\Delta$NLL values are shown which correspond to the minimized negative log likelihood for a given Wilson coefficient value.

1D-likelihood scan for $c_{Qu}^{(8)}$. The observed and expected limits at 68\% and 95\% confidence level (CL) can be read off by finding the intersection of the likelihood scan curve with the 1$\sigma$ and 2$\sigma$ lines. The new-physics scale is set to $\Lambda = 1 \, \text{TeV}$. In the table the $2\Delta$NLL values are shown which correspond to the minimized negative log likelihood for a given Wilson coefficient value.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the combined fit, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the combined fit, when all other coupling modifiers are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the combined fit, when all other coupling modifiers are fixed to their SM predictions.

Likelihood observed contours at 68% CL (solid line) and 95% CL (dashed line) in the ($\kappa_{\lambda}$, $\kappa_{2V}$) parameter space, when all other coupling modifiers are fixed to their SM predictions. The best-fit value, denoted by a cross in the plot, corresponds to the very first entry in the table.

Likelihood expected contours at 68% CL (teal shaded region) and 95% CL (yellow shaded region) in the ($\kappa_{\lambda}$, $\kappa_{2V}$) parameter space, when all other coupling modifiers are fixed to their SM predictions. In the plot, the SM prediction is indicated by the star.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{gghh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross. The best-fit value corresponds to the very first entry in the table.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{gghh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{tthh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross. The best-fit value corresponds to the very first entry in the table.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{tthh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{H}$, $c_{H◻}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross. The best-fit value corresponds to the very first entry in the table.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{H}$, $c_{H◻}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross.

The acceptance times efficiency [in %] for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ SLT channel. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ LTT channel. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ SLT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ LTT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ SLT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ LTT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Summary of the expected (blue) and observed (orange) one-dimensional constraints at 68% (solid error bars) and 95% (dashed error bars) CL on the coefficients of selected HEFT (top) and SMEFT (bottom) operators describing Higgs boson interactions affecting $HH$ production through gluon-gluon fusion. The results are obtained from one dimensional scans of the profile log-likelihood assuming that all other Wilson coefficients are fixed to their SM values. The contribution from VBF $HH$ production is neglected.