Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected (dashed black line) and observed (solid red line) 95% CL exclusion limits on the higgsino simplified model being considered. These are shown with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively. The limits set by the latest ATLAS searches using the soft lepton and disappearing track signatures are illustrated by the blue and green regions, respectively, while the limit imposed by the LEP experiments is shown in gray. The dot-dashed gray line indicates the predicted mass-splitting for the pure higgsino scenario.

Expected and observed CLs values per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Expected and observed CLs values per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Expected and observed cross-section upper-limit per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Expected and observed cross-section upper-limit per signal point represented by the grey numbers. The expected (dashed) and observed (solid) 95% CL exclusion limits are overlaid along with $\pm 1\sigma_{\mathrm{exp}}$ (yellow band) from experimental systematic and statistical uncertainties, and with $\pm 1\sigma_{\mathrm{theory}}^{\mathrm{SUSY}}$ (red dotted lines) from signal cross-section uncertainties, respectively.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Truth-level signal acceptances for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$) in a SR with the $S(d_0)$ requirement removed. The acceptance is defined as the fraction of accepted events divided by the total number of events in the generator-level signal Monte Carlo simulation, where the signal candidate track is identified as the charged particle with the largest distance between the interaction vertex and the secondary vertex of the higgsino decays.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-Low for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Signal efficiencies in SR-High for each production process ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$), defined by the number of events of reconstructed-level signal simulation divided by the number of events obtained at generator level, where the $S(d_0)$ selecton efficiency has the largest impact. The higgsino decay products from $\Delta \mathrm{m}(\tilde{\chi}_1^\pm,\tilde{\chi}_1^0) < 0.4$ GeV signal have $p_{\mathrm{T}}$ too low to be reconstructed as the signal candidate tracks, and therefore the identified signal candidate tracks are typically from pile-up collisions or underlying events similar to the QCD track background, causing a low $S(d_0)$ selection efficiency in these plots.

Event selection cutflows for signal samples with $m(\tilde{\chi}_{1}^0)$ = 150 GeV and $\Delta m(\tilde{\chi}_{1}^\pm, \tilde{\chi}_{1}^0)$ = 1.5, 1.0, and 0.75 GeV, including all six production processes ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$). The cross-section used to obtain the initial number of events ($\sigma(\mathrm{n}_{\mathrm{jets}}) \geq 1$) refers to an emission of at least one gluon or quark with $p_{\mathrm{T}} > 50$ GeV at the parton level.

Event selection cutflows for signal samples with $m(\tilde{\chi}_{1}^0)$ = 150 GeV and $\Delta m(\tilde{\chi}_{1}^\pm, \tilde{\chi}_{1}^0)$ = 0.5, 0.35, and 0.25 GeV, including all six production processes ($\tilde{\chi}_1^\pm \tilde{\chi}_1^0$, $\tilde{\chi}_1^\pm \tilde{\chi}_2^0$, $\tilde{\chi}_1^+ \tilde{\chi}_1^-$, and $\tilde{\chi}_2^0 \tilde{\chi}_1^0$). The cross-section used to obtain the initial number of events ($\sigma(\mathrm{n}_{\mathrm{jets}}) \geq 1$) refers to an emission of at least one gluon or quark with $p_{\mathrm{T}} > 50$ GeV at the parton level.

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.

Impact of systematic uncertainties on measured invisible width of the Z boson.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Distribution of $m_{\mathrm{JJ}}$ in the SR (black points) compared to the background-only expectation, where the shape is derived from the CR and the normalisation is fitted in the SR.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model A (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model A (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model A, the predicted cross-section is shown for $g_q=0.05$ and $g_{q_d} = 0.2$.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model B (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model B (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model B, the predicted cross-section is shown for $g_q=0.15$ and $g_{q_d} = 0.5$.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model C (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model C (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model C, the predicted cross-section is shown for $g_q=0.15$ and $g_{q_d}= 0.5$.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model D (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model D (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model D, the predicted cross-section is shown for $g_q=0.05$ and $g_{q_d} = 0.2$.

Relative efficiency (in %) of each stage of the analysis selection with respect to the previous one for models A through D and $m_{Z^\prime}=2.5$ TeV. Approximately 90000 events were generated per signal mass hypothesis.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63717.4 $\pm$ 252.4, NPowHeg truth =338714.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63855.8 $\pm$ 252.7 , NPowHeg truth =338708.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64809.8 $\pm$ 254.6, NPowHeg truth =634285.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Truth $P_{T}^{Z}$ distribution after reweghting to data.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63717.4 $\pm$ 252.4, NPowHeg truth =338714.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63855.8 $\pm$ 252.7 , NPowHeg truth =338708.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64809.8 $\pm$ 254.6, NPowHeg truth =634285.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63717.4 $\pm$ 252.4, NPowHeg truth =338714.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63855.8 $\pm$ 252.7 , NPowHeg truth =338708.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64809.8 $\pm$ 254.6, NPowHeg truth =634285.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Unfolded $M(l^{+}\gamma\gamma)$ distribution for $Z \to ll\gamma\gamma$ process with dressed leptons.

Unfolded $M(l^{+}\gamma\gamma)$ distribution for $Z \to ll\gamma\gamma$ process with bare leptons.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for leading photon.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with bare leptons for leading photon.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for second photon.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with bare leptons for second photon.

Unfolded dRgg distribution for $Z \to ll\gamma\gamma$ process with dressed leptons.

Unfolded dRgg distribution for $Z \to ll\gamma\gamma$ process with bare leptons.

Unfolded $P_{T}^{\gamma 1}$ distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for leading photon.

Unfolded $P_{T}^{\gamma 1}$ distribution for $Z \to ll\gamma\gamma$ process with bare leptons for leading photon.

Unfolded $P_{T}^{\gamma 2}$ distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for second photon.

Unfolded $P_{T}^{\gamma 2}$ distribution for $Z \to ll\gamma\gamma$ process with bare leptons for second photon.

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

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

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

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

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

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

Measured $p_T$ cross sections in full-lepton phase space for 1.2 < |y| < 1.6.

Measured $p_T$ cross sections in full-lepton phase space for 1.6 < |y| < 2.0.

Measured $p_T$ cross sections in full-lepton phase space for 2.0 < |y| < 2.4.

Measured $p_T$ cross sections in full-lepton phase space for 2.4 < |y| < 2.8.

Measured $p_T$ cross sections in full-lepton phase space for 2.8 < |y| < 3.6.

Measured cross sections in full-lepton phase space as a function of |y|.

Normalised measured $p_T$ cross sections in full-lepton phase space for |y| < 1.6.