The post-fit BDT score distribution for the direct stau channel, showing the scores for BDT4, before the selections on the BDT score is made. The black arrow depicts the BDT score selection for the SR-BDT. A few example SUSY scenarios targeted by each BDT are overlaid for illustration.

The post-fit kinematic distribution for the Intermediate stau channel, showing the $m_{\mathrm{T}2}$ distribution in SR-C1C1-LM, before the selection on $m_{\mathrm{T}2}$ is made. The black arrow depicts the selection for the signal region. An example SUSY scenario is overlaid for illustration.

The post-fit kinematic distribution for the Intermediate stau channel, showing the $m_{\mathrm{T}2}$ distribution in SR-C1N2OS-LM, before the selection on $m_{\mathrm{T}2}$ is made. The black arrow depicts the selection for the signal region. An example SUSY scenario is overlaid for illustration.

The post-fit kinematic distribution for the Intermediate stau channel, showing the $m_{\mathrm{T}2}$ distribution in SR-C1N2SS-LM, before the selection on $m_{\mathrm{T}2}$ is made. The black arrow depicts the selection for the signal region. An example SUSY scenario is overlaid for illustration.

The post-fit kinematic distribution for the Intermediate $Wh$ channel, showing the $m_{\mathrm{T}2}$ distribution in SR-Wh-LM, before the selection on $m_{\mathrm{T}2}$ is made. The black arrow depicts the selection for the signal region. An example SUSY scenario is overlaid for illustration.

The post-fit kinematic distribution for the Intermediate $Wh$ channel, showing the $m_{\mathrm{Tsum}}$ distribution in SR-Wh-HM, before the selection on $m_{\mathrm{Tsum}}$ is made. The black arrow depicts the selection for the signal region. An example SUSY scenario is overlaid for illustration.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$ production.

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$ production.

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$ production.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$ production.

$+1\sigma$ observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$ production.

$-1\sigma$ observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$ production.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L}\tilde{\tau}_{L}$ production.

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L}\tilde{\tau}_{L}$ production.

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L}\tilde{\tau}_{L}$ production.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L}\tilde{\tau}_{L}$ production.

$+1\sigma$ observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L}\tilde{\tau}_{L}$ production.

$-1\sigma$ observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{L}\tilde{\tau}_{L}$ production.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{R}\tilde{\tau}_{R}$ production.

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{R}\tilde{\tau}_{R}$ production.

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{R}\tilde{\tau}_{R}$ production.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{R}\tilde{\tau}_{R}$ production.

$+1\sigma$ observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{R}\tilde{\tau}_{R}$ production.

$-1\sigma$ observed 95\% CL exclusion limits on the simplified models of $\tilde{\tau}_{R}\tilde{\tau}_{R}$ production.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ production decaying via staus.

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ production decaying via staus.

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ production decaying via staus.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ production decaying via staus.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$).

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$).

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$).

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$).

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$) using the OS signal regions.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$) using the OS signal regions.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$) using the SS signal regions.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via staus (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_1^0)$) using the SS signal regions.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$).

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$).

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$).

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$).

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$) using the OS signal regions.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$) using the OS signal regions.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$) using the SS signal regions.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$ (parameterised as $m(\tilde{\chi}_1^{\pm}/\tilde{\chi}_2^0)$ and $m(\tilde{\chi}_2^0)-m(\tilde{\chi}_1^0)$) using the SS signal regions.

Expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$.

$+1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$.

$-1\sigma$ expected 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$.

Observed 95\% CL exclusion limits on the simplified models of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ production decaying via $Wh$.

Observed upper limit on the signal cross section in fb for the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

Observed upper limit on the signal cross section in fb for the production of $\tilde{\tau}_{L}\tilde{\tau}_{L}$.

Observed upper limit on the signal cross section in fb for the production of $\tilde{\tau}_{R}\tilde{\tau}_{R}$.

The best expected signal region used to set limits in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The best expected signal region used to set limits in models of the production of $\tilde{\tau}_{L}\tilde{\tau}_{L}$.

The best expected signal region used to set limits in models of the production of $\tilde{\tau}_{R}\tilde{\tau}_{R}$.

The acceptance $A$ of the direct stau production signal region SR-BDT1 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The efficiency $\epsilon$ of the direct stau production signal region SR-BDT1 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The acceptance $A$ of the direct stau production signal region SR-BDT2 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The efficiency $\epsilon$ of the direct stau production signal region SR-BDT2 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The acceptance $A$ of the direct stau production signal region SR-BDT3 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The efficiency $\epsilon$ of the direct stau production signal region SR-BDT3 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The acceptance $A$ of the direct stau production signal region SR-BDT4 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

The efficiency $\epsilon$ of the direct stau production signal region SR-BDT4 in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$.

Cutflow event yields of the direct stau production SRs in models of the production of $\tilde{\tau}_{L,R}\tilde{\tau}_{L,R}$. All yields correspond to weighted events, so that effects from reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalised to the integrated luminosity of the data sample, $\int L dt = 139$ fb$^{-1}$. The preliminary event reduction is a centralised stage requiring the event to: have at least two hadronically-decaying taus; be selected and matched by the asymmetric ditau trigger (with plateau cuts); have at least two medium taus of opposite sign and have $m_{\mathrm{T2}} > 30$ GeV. Cuts except the BDT score cuts are applied consecutively, while the BDT score cuts are applied one at a time.

Observed upper limit on the signal cross section in fb for the production of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$.

Observed upper limit on the signal cross section in fb for the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$.

The acceptance $A$ of the Intermediate stau signal region SR-C1C1-LM in models of the production of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ decaying via staus.

The efficiency $\epsilon$ of the Intermediate stau signal region SR-C1C1-LM in models of the production of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ decaying via staus.

The acceptance $A$ of the Intermediate stau signal region SR-C1C1-HM in models of the production of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ decaying via staus.

The efficiency $\epsilon$ of the Intermediate stau signal region SR-C1C1-HM in models of the production of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$ decaying via staus.

The acceptance $A$ of the Intermediate stau signal region SR-C1N2OS-LM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The efficiency $\epsilon$ of the Intermediate stau signal region SR-C1N2OS-LM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The acceptance $A$ of the Intermediate stau signal region SR-C1N2SS-LM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The efficiency $\epsilon$ of the Intermediate stau signal region SR-C1N2SS-LM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The acceptance $A$ of the Intermediate stau signal region SR-C1N2OS-HM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The efficiency $\epsilon$ of the Intermediate stau signal region SR-C1N2OS-HM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The acceptance $A$ of the Intermediate stau signal region SR-C1N2SS-HM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

The efficiency $\epsilon$ of the Intermediate stau signal region SR-C1N2SS-HM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$ decaying via staus.

Cutflow event yields of the Intermediate stau SRs in models of the production of $\tilde{\chi}_1^{+}\tilde{\chi}_1^{-}$. All yields correspond to weighted events, so that effects from reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalised to the integrated luminosity of the data sample, $\int L dt = 139$ fb$^{-1}$. The preliminary event reduction is a centralised stage requiring the event to: have at least two hadronically-decaying medium taus; be selected and matched by the asymmetric ditau trigger or ditau+MET trigger with offline cuts using HLT threshold values. The different yields of preliminary event reduction are caused by different trigger scale factor.

Cutflow event yields of the Intermediate stau SRs in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^{0}$. All yields correspond to weighted events, so that effects from reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalised to the integrated luminosity of the data sample, $\int L dt = 139$ fb$^{-1}$. The preliminary event reduction is a centralised stage requiring the event to: have at least two hadronically-decaying medium taus; be selected and matched by the asymmetric ditau trigger or ditau+MET trigger with offline cuts using HLT threshold values. The different yields of preliminary event reduction are caused by different trigger scale factor.

Observed upper limit on the signal cross section in fb for the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$.

The best expected signal region used to set limits in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$.

The acceptance $A$ of the Intermediate $Wh$ signal region SR-Wh-LM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$.

The efficiency $\epsilon$ of the Intermediate $Wh$ signal region SR-Wh-LM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$.

The acceptance $A$ of the Intermediate $Wh$ signal region SR-Wh-HM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$.

The efficiency $\epsilon$ of the Intermediate $Wh$ signal region SR-Wh-HM in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$.

Cutflow event yields of the Intermediate $Wh$ SRs in models of the production of $\tilde{\chi}_1^{\pm}\tilde{\chi}_2^0$ decaying via $Wh$. All yields correspond to weighted events, so that effects from reconstruction efficiencies, trigger corrections, pileup reweighting, etc., are included. They are normalised to the integrated luminosity of the data sample, $\int L dt = 139$ fb$^{-1}$. The preliminary event reduction is a centralised stage requiring the event to: have one single lepton and two hadronically decaying medium taus; be selected and matched by the single lepton trigger with plateau cuts.

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.

Upper limits on $\mathcal{B}(H\rightarrow aa\rightarrow 4\gamma)$ at 95% CL as a function of the axion mass and for ALP-photon coupling $C_{a\gamma\gamma}=0.01$.

Upper limits on $\mathcal{B}(H\rightarrow aa\rightarrow 4\gamma)$ at 95% CL as a function of the axion mass and for ALP-photon coupling $C_{a\gamma\gamma}=5\times10^{-4}$

Upper limits on $\mathcal{B}(H\rightarrow aa\rightarrow 4\gamma)$ at 95% CL as a function of the axion mass and for ALP-photon coupling $C_{a\gamma\gamma}=10^{-5}$

Upper limits on $B(H\rightarrow 4\gamma)$ at 95% CL as a function of the signal mass hypothesis and for the assumption of promptly decaying ALPs.

Identification efficiency for isolated photons that pass the tight photon ID requirements in dependence of the opening angle in $\Delta R$ of both decay photons of the axion and the decay radius of the originating axion from the primary vertex for axions with masses between 100 MeV and 300 MeV. All four axion decay photons on MC truth level are required to fulfil $|\eta|<2.5$ and $p_T>5$ GeV. Only a mild dependency on the decay radius can be seen, since the distance between the two photons is small and hence the combined decay signature in the calorimeter still points to the primary vertex.

Efficiency for classified merged isolated photons in dependence of the opening angle in $\Delta R$ of both decay photons of the axion and the decay radius of the originating axion from the primary vertex for axions with masses between 100 MeV and 300 MeV. All four axion decay photons on MC truth level are required to fulfil $|\eta|<2.5$ and $p_T>5$ GeV. Only a mild dependency on the decay radius can be seen, since the distance between the two photons is small and hence the combined decay signature in the calorimeter still points to the primary vertex.

Limits on the ALP mass and coupling to photons at 95\%~CL, assuming $\mathcal{B}(a\rightarrow\gamma\gamma) = 1$, $\Lambda = 1\,TeV$ with $|C_{aH}^{\text{eff}}|$ = $1$ (solid line) and $|C_{aH}^{\text{eff}}|$ = $0.1$.

The number of data and estimated background events in the signal region of the most sensitive categories..

Expected and observed 95% CL upper limits on the branching ratio of the Higgs boson decay to $Za$ times the branching ratio $a $->$ \gamma \gamma$ as a function of the $a$ particle mass in the merged ($m_a \le 2$ GeV) and the resolved ($m_a > 2$ GeV) categories. The figure uses $pp$ collision data at $\sqrt{s}=13$ TeV corresponding to 139 fb$^{-1}$.

ATLAS observed 95% CL exclusion contours limits in terms of the ALP mass and its effective coupling to photons, $|C_{\gamma\gamma}|/\Lambda$, for different values of the Higgs coupling to $Za$, $|C_{Zh}|/\Lambda$. The figure uses $pp$ collision data at $\sqrt{s}=13$ TeV corresponding to 139 fb$^{-1}$. The value 100 in the tables is presentational, denoting the upper edge of the plot.

Isoscalar WIMP-nucleon elastic coupling limit for Operator 9

Isoscalar WIMP-nucleon elastic coupling limit for Operator 12

Isoscalar WIMP-nucleon elastic coupling limit for Operator 7

Isoscalar WIMP-nucleon elastic coupling limit for Operator 13

Isoscalar WIMP-nucleon elastic coupling limit for Operator 3

Isoscalar WIMP-nucleon elastic coupling limit for Operator 5

Isoscalar WIMP-nucleon elastic coupling limit for Operator 11

Isoscalar WIMP-nucleon elastic coupling limit for Operator 14

Isoscalar WIMP-nucleon elastic coupling limit for Operator 4

Isoscalar WIMP-nucleon elastic coupling limit for Operator 10

Isoscalar WIMP-nucleon elastic coupling limit for Operator 1

Isoscalar WIMP-nucleon elastic coupling limit for Operator 6

Isoscalar WIMP-nucleon elastic coupling limit for Operator 15

Isovector WIMP-nucleon elastic coupling limit for Operator 3

Isovector WIMP-nucleon elastic coupling limit for Operator 10

Isovector WIMP-nucleon elastic coupling limit for Operator 13

Isovector WIMP-nucleon elastic coupling limit for Operator 12

Isovector WIMP-nucleon elastic coupling limit for Operator 8

Isovector WIMP-nucleon elastic coupling limit for Operator 1

Isovector WIMP-nucleon elastic coupling limit for Operator 4

Isovector WIMP-nucleon elastic coupling limit for Operator 5

Isovector WIMP-nucleon elastic coupling limit for Operator 7

Isovector WIMP-nucleon elastic coupling limit for Operator 11

Isovector WIMP-nucleon elastic coupling limit for Operator 15

Isovector WIMP-nucleon elastic coupling limit for Operator 9

Isovector WIMP-nucleon elastic coupling limit for Operator 14

Isovector WIMP-nucleon elastic coupling limit for Operator 6

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 12

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 3

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 9

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 7

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 10

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 15

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 8

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 4

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 5

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 6

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 13

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 11

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 14

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 1

Isovector WIMP-nucleon inelastic coupling limit for Operator 13

Isovector WIMP-nucleon inelastic coupling limit for Operator 1

Isovector WIMP-nucleon inelastic coupling limit for Operator 4

Isovector WIMP-nucleon inelastic coupling limit for Operator 7

Isovector WIMP-nucleon inelastic coupling limit for Operator 9

Isovector WIMP-nucleon inelastic coupling limit for Operator 3

Isovector WIMP-nucleon inelastic coupling limit for Operator 11

Isovector WIMP-nucleon inelastic coupling limit for Operator 5

Isovector WIMP-nucleon inelastic coupling limit for Operator 10

Isovector WIMP-nucleon inelastic coupling limit for Operator 8

Isovector WIMP-nucleon inelastic coupling limit for Operator 14

Isovector WIMP-nucleon inelastic coupling limit for Operator 15

Isovector WIMP-nucleon inelastic coupling limit for Operator 6

Isovector WIMP-nucleon inelastic coupling limit for Operator 12

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.

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.