Observed and expected 95 % CL exclusion limits on $\tan\beta$ as a function of $m_{H^{\pm}}$, shown in the context of the mh125 scenario, for $m_{H^{\pm}}>150$ GeV and $(1 \leq \tan\beta \leq 60)$. The surrounding shaded bands correspond to the 1$\sigma$ and 2$\sigma$ confidence intervals around the expected limit.

Observed and expected 95 % CL exclusions on $\tan\beta$ as a function of $m_{H^{\pm}}$ derived from exclusion limits on $\mathrm{\cal{B}}(t\to bH^+)\times \mathrm{\cal{B}}(H^+ \to \tau \nu)$, for Type I 2HDM low mass scenario. The surrounding shaded bands correspond to the 1$\sigma$ and 2$\sigma$ confidence intervals around the expected limit.

Observed and expected 95 % CL exclusions on $\tan\beta$ as a function of $m_{H^{\pm}}$ derived from exclusion limits on $\mathrm{\cal{B}}(t\to bH^+)\times \mathrm{\cal{B}}(H^+ \to \tau \nu)$, for Lepton Specific 2HDM low mass scenario. The surrounding shaded bands correspond to the 1$\sigma$ and 2$\sigma$ confidence intervals around the expected limit.

Breakdown of relative systematic uncertainties in the measured ${t\bar{t}}$ cross-section. The quoted uncertainties are obtained by repeating the fit with a group of nuisance parameters fixed to their fitted values and subtracting in quadrature the resulting total uncertainty from the uncertainty of the complete fit. However, the total uncertainty is not the quadratic sum of the grouped impacts, as this approach neglects the correlation among the different groups.

The figure shows the measurement of the top quark pair production cross-section, $\sigma_{t\bar{t}}$, in lead-lead (Pb+Pb) collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV, based on 1.9 nb$^{-1}$ of data.

Observed and expected 95% CL upper limits on the total cross-section σ($pp$ → $T$ → $Ht/Zt$) as a function of $T$-quark mass in the SU(2) doublet representation assuming $\kappa$=0.5. The expected limits for the individual analyses are shown. The $HtZt$ analysis is only included in the limit calculation for $m_{\mathrm{T}}$ < 2.1 TeV.

Observed 95% CL exclusion limits on the total cross-section σ($pp$ → $T$ → $Ht/Zt$) as a function of the universal coupling constant and the $T$-quark mass in the SU(2) singlet representation. Limits are only presented in the regime $\Gamma/m_{\mathrm{T}}$ < 50%, where the theory calculations are known to be valid.

Expected 95% CL exclusion limits on the total cross-section σ($pp$ → $T$ → $Ht/Zt$) as a function of the universal coupling constant and the $T$-quark mass in the SU(2) singlet representation. Limits are only presented in the regime $\Gamma/m_{\mathrm{T}}$ < 50%, where the theory calculations are known to be valid.

Observed 95% CL exclusion limits on the total cross-section σ($pp$ → $T$ → $Ht/Zt$) as a function of the universal coupling constant and the $T$-quark mass in the SU(2) doublet representation. Limits are only presented in the regime $\Gamma/m_{\mathrm{T}}$ < 50%, where the theory calculations are known to be valid.

Expected 95% CL exclusion limits on the total cross-section σ($pp$ → $T$ → $Ht/Zt$) as a function of the universal coupling constant and the $T$-quark mass in the SU(2) doublet representation. Limits are only presented in the regime $\Gamma/m_{\mathrm{T}}$ < 50%, where the theory calculations are known to be valid.

Obs.Comb 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) singlet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.Comb 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) singlet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.$+1\sigma$ for the 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) singlet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.$-1\sigma$ for the 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) singlet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.$+2\sigma$ for the 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) singlet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.$-2\sigma$ for the 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) singlet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.OSML 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) singlet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.HTZT 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) singlet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.Monotop 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) singlet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Obs.Comb 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) doublet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.Comb 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) doublet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.$+1\sigma$ for the 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) doublet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.$-1\sigma$ for the 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) doublet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.$+2\sigma$ for the 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) doublet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.$-2\sigma$ for the 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) doublet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.OSML 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) doublet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.HTZT 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) doublet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Observed upper limits at 95% CL on the $T$ quark mass as a function of the relative resonance width ($\Gamma_{\mathrm{T}}/m_{\mathrm{T}}$) and the relative coupling parameter $\xi_W$.

Expected upper limits at 95% CL on the $T$ quark mass as a function of the relative resonance width ($\Gamma_{\mathrm{T}}/m_{\mathrm{T}}$) and the relative coupling parameter $\xi_W$.

Observed and expected 95% CL upper limits on the total cross-section σ($pp$ → $T$ → $Ht/Zt$) as a function of $T$-quark mass in the SU(2) singlet representation assuming $\kappa$=0.7. The expected limits for the individual analyses are shown.

Observed and expected 95% CL upper limits on the total cross-section σ($pp$ → $T$ → $Ht/Zt$) as a function of $T$-quark mass in the SU(2) doublet representation assuming $\kappa$=0.7. The expected limits for the individual analyses are shown.

Observed and expected 95% CL upper limits on the total cross-section σ($pp$ → $T$ → $Ht/Zt$) as a function of $T$-quark mass in the SU(2) singlet representation assuming $\kappa$=0.3. The observed limits for the individual analyses are shown. The $HtZt$ analysis is only included in the limit calculation for $m_{\mathrm{T}}$ < 2.1 TeV.

Observed and expected 95% CL upper limits on the total cross-section σ($pp$ → $T$ → $Ht/Zt$) as a function of $T$-quark mass in the SU(2) singlet representation assuming $\kappa$=0.5. The observed limits for the individual analyses are shown. The $HtZt$ analysis is only included in the limit calculation for $m_{\mathrm{T}}$ < 2.1 TeV.

Observed and expected 95% CL upper limits on the total cross-section σ($pp$ → $T$ → $Ht/Zt$) as a function of $T$-quark mass in the SU(2) singlet representation assuming $\kappa$=0.7. The observed limits for the individual analyses are shown.

Observed and expected 95% CL upper limits on the total cross-section σ($pp$ → $T$ → $Ht/Zt$) as a function of $T$-quark mass in the SU(2) doublet representation assuming $\kappa$=0.3. The observed limits for the individual analyses are shown. The $HtZt$ analysis is only included in the limit calculation for $m_{\mathrm{T}}$ < 2.1 TeV.

Observed and expected 95% CL upper limits on the total cross-section σ($pp$ → $T$ → $Ht/Zt$) as a function of $T$-quark mass in the SU(2) doublet representation assuming $\kappa$=0.5. The observed limits for the individual analyses are shown. The $HtZt$ analysis is only included in the limit calculation for $m_{\mathrm{T}}$ < 2.1 TeV.

Observed and expected 95% CL upper limits on the total cross-section σ($pp$ → $T$ → $Ht/Zt$) as a function of $T$-quark mass in the SU(2) doublet representation assuming $\kappa$=0.7. The observed limits for the individual analyses are shown.

Obs.Comb 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) singlet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.Comb 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) singlet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.$+1\sigma$ for the 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) singlet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.$-1\sigma$ for the 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) singlet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.$+2\sigma$ for the 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) singlet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.$-2\sigma$ for the 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) singlet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Obs.OSML 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) singlet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Obs.HTZT 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) singlet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Obs.Monotop 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) singlet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Obs.Comb 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) doublet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.Comb 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) doublet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.$+1\sigma$ for the 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) doublet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.$-1\sigma$ for the 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) doublet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.$+2\sigma$ for the 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) doublet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Exp.$-2\sigma$ for the 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) doublet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Obs.OSML 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) doublet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Obs.HTZT 95% CL exclusion limits on the universal coupling constant $\kappa$ as a function of the $T$ quark mass in the SU(2) doublet representations for the combination. Limits are only presented in the regime $\Gamma_{\mathrm{T}}/m_{\mathrm{T}} < 50\%$, where the theory calculations are known to be valid.

Figure 8(left) of the article. Differential fiducial cross-section of the Z boson $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ b-jets. The thin inner band corresponds to the statistical uncertainty of the data, and the outer band to statistical and systematic uncertainties of the data, added in quadrature.

Figure 8(right) of the article. Differential fiducial cross-section of the leading $b$-jet $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ b-jets. The thin inner band corresponds to the statistical uncertainty of the data, and the outer band to statistical and systematic uncertainties of the data, added in quadrature.

Figure 9 of the article. Differential fiducial cross-section of the $\Delta R$ between the $Z$ boson and the leading $b$-jet in events with $Z (\rightarrow ll) \ge 1 $ b-jets. The thin inner band corresponds to the statistical uncertainty of the data, and the outer band to statistical and systematic uncertainties of the data, added in quadrature.

Figure 10(left) of the article. Differential fiducial cross-section of the $\Delta \phi$ between the leading and sub-leading $b$-jets in events with $Z (\rightarrow ll) \ge 2 $ b-jets. The thin inner band corresponds to the statistical uncertainty of the data, and the outer band to statistical and systematic uncertainties of the data, added in quadrature.

Figure 10(right) of the article. Differential fiducial cross-section of the invariant mass of the leading and sub-leading $b$-jets in events with $Z (\rightarrow ll) \ge 2 $ b-jets. The thin inner band corresponds to the statistical uncertainty of the data, and the outer band to statistical and systematic uncertainties of the data, added in quadrature.

Figure 11(left) of the article. Differential fiducial cross-section of the Z boson $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. The thin inner band corresponds to the statistical uncertainty of the data, and the outer band to statistical and systematic uncertainties of the data, added in quadrature.

Figure 11(right) of the article. Differential fiducial cross-section of the leading $c$-jet $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. The thin inner band corresponds to the statistical uncertainty of the data, and the outer band to statistical and systematic uncertainties of the data, added in quadrature.

Figure 12 of the article. Differential fiducial cross-section of the leading $c$-jet $x_F$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. The thin inner band corresponds to the statistical uncertainty of the data, and the outer band to statistical and systematic uncertainties of the data, added in quadrature.

Correction factors from parton level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level with R 0.3 cone labelling (see Tables 18 & 24) for the differential cross-section of the Z boson $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ b-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties for what concerns the passage from hadron level flavour dressing to hadron level cone labelling, plus statistical uncertainty of Pythia8 simulations used for parton to hadron corrections. See Table 4.

Correction factors from parton level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level with R 0.3 cone labelling (see Tables 19 & 25) for the differential cross-section of the Z boson $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties for what concerns the passage from hadron level flavour dressing to hadron level cone labelling, plus statistical uncertainty of Pythia8 simulations used for parton to hadron corrections. See Table 9.

Correction factors from parton level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level with R 0.3 cone labelling (see Tables 20 & 26) for the differential cross-section of the leading $b$-jet $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ b-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties for what concerns the passage from hadron level flavour dressing to hadron level cone labelling, plus statistical uncertainty of Pythia8 simulations used for parton to hadron corrections. See Table 5.

Correction factors from parton level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level with R 0.3 cone labelling (see Tables 21 & 27) for the differential cross-section of the leading $c$-jet $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties for what concerns the passage from hadron level flavour dressing to hadron level cone labelling, plus statistical uncertainty of Pythia8 simulations used for parton to hadron corrections. See Table 10.

Correction factors from parton level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level with R 0.3 cone labelling (see Tables 22 & 28) for the differential cross-section of the leading $c$-jet $x_F$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties for what concerns the passage from hadron level flavour dressing to hadron level cone labelling, plus statistical uncertainty of Pythia8 simulations used for parton to hadron corrections. See Table 11.

Correction factors from parton level done with flavour dressing algorithms as seen in Ref. [2] of the paper of the paper to hadron level with R 0.3 cone labelling (see Tables 23 & 29) for the differential cross-section of $\Delta R$ between the $Z$ boson and the leading $b$-jet in events with $Z (\rightarrow ll) \ge 1 $ b-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties for what concerns the passage from hadron level flavour dressing to hadron level cone labelling, plus statistical uncertainty of Pythia8 simulations used for parton to hadron corrections. See Table 6.

Correction factors from hadron level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level R 0.3 cone labelling for the differential cross-section of the Z boson $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ b-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties. See Table 12.

Correction factors from hadron level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level with R 0.3 cone labelling for the differential cross-section of the Z boson $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties. See Table 13.

Correction factors from hadron level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level with R 0.3 cone labelling for the differential cross-section of the leading $b$-jet $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ b-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties. See Table 14.

Correction factors from hadron level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level with R 0.3 cone labelling for the differential cross-section of the leading $c$-jet $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties. See Table 15.

Correction factors from hadron level done with flavour dressing algorithms as seen in Ref. [2] of the paper to hadron level with R 0.3 cone labelling for the differential cross-section of the leading $c$-jet $x_F$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties. See Table 16.

Correction factors from hadron level done with flavour dressing algorithms as seen in Ref. [2] of the paper of the paper to hadron level with R 0.3 cone labelling for the differential cross-section of $\Delta R$ between the $Z$ boson and the leading $b$-jet in events with $Z (\rightarrow ll) \ge 1 $ b-jets. Systematic errors include Sherpa 2.2.11 and MGaMC+Py8 FxFx uncertainties. See Table 17.

Correction factors from parton level to hadron level for the differential cross-section of the Z boson $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ b-jets. Systematic errors statistical uncertainty of Pythia8 simulations. See Table 12.

Correction factors from parton level to hadron level for the differential cross-section of the Z boson $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. Systematic errors include statistical uncertainty of Pythia8 simulations. See Table 13.

Correction factors from parton level to hadron level for the differential cross-section of the leading $b$-jet $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ b-jets. Systematic errors include statistical uncertainty of Pythia8 simulations. See Table 14.

Correction factors from parton level to hadron level for the differential cross-section of the leading $c$-jet $p_T$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. Systematic errors include statistical uncertainty of Pythia8 simulations. See Table 15.

Correction factors from parton level to hadron level for the differential cross-section of the leading $c$-jet $x_F$ in events with $Z (\rightarrow ll) \ge 1 $ c-jets. Systematic errors include statistical uncertainty of Pythia8 simulations. See Table 16.

Correction factors from parton level to hadron level for the differential cross-section of $\Delta R$ between the $Z$ boson and the leading $b$-jet in events with $Z (\rightarrow ll) \ge 1 $ b-jets. Systematic errors include statistical uncertainty of Pythia8 simulations. See Table 17.

Observed event yields in $\textrm{SR}$ compared with post-fit expectations from Monte Carlo simulations, as a function of the scalar sum of lepton and jet transverse momenta, $H_{\mathrm{T}}$. The last bin includes overflow events. `Signal (prod.)' and `Signal (dec.)' refer to the single-top-quark production and top-quark pair decay signal contributions, respectively.

Expected and observed 95$\%$ CL upper limits on the branching ratio corresponding to the decay of a top quark to a muon and a tau lepton through a cLFV process for scalar, vector and tensor couplings. In the vector case, all vector operators are assumed to contribute simultaneously with the same effective coupling strength. The table shows the endpoint values of $\mathcal{B}(\mathrm{t}\to\mu\tau\mathrm{u})$ and $\mathcal{B}(\mathrm{t}\to\mu\tau\mathrm{c})$, which are interpolated linearly for each limit.

Expected and observed 95$\%$ CL upper limits on the Wilson coefficients for scalar, vector and tensor couplings of a top quark to a muon and a tau lepton through a cLFV process. In the vector case, all vector operators are assumed to contribute simultaneously with the same effective coupling strength. The table shows the endpoint values of the $|c^{\mu\tau\mathrm{ut}}|$ and $|c^{\mu\tau\mathrm{ct}}|$ Wilson coefficients, which are interpolated assuming a linear relationship between $\mathcal{B}(\mathrm{t}\to\mu\tau\mathrm{u})$ and $\mathcal{B}(\mathrm{t}\to\mu\tau\mathrm{c})$.

Observed event yields in $\textrm{SR}$ compared with pre-fit expectations from Monte Carlo simulations, as a function of the scalar sum of lepton and jet transverse momenta, $H_{\mathrm{T}}$. The last bin includes overflow events. The signal yields represent a leptoquark mass of $m_{S_{1}} = 1$ TeV and a coupling strength of $\lambda^{\textrm{LQ}} = 2.0$.

Observed event yields in $\textrm{SR}$ compared with post-fit expectations from Monte Carlo simulations, as a function of the scalar sum of lepton and jet transverse momenta, $H_{\mathrm{T}}$. The last bin includes overflow events. The signal yields represent a leptoquark mass of $m_{S_{1}} = 1$ TeV and a coupling strength of $\lambda^{\textrm{LQ}} = 2.0$.

Observed and expected exclusion 95$\%$ CL upper limits on the $S_{1}$ leptoquark coupling strength $\lambda^{\textrm{LQ}}$ as a function of LQ mass, $m_{S_{1}}$.

Theoretical cross-sections for single-top-quark production and top-quark decays through cLFV interactions involving vector, scalar and tensor EFT Wilson coefficients. The column titled as $c^{(ijk3)}_{\textrm{vector}}$ represents the individual cross-section contributions from each of $c_{\mathsf{lq}}^{-(ijk3)}$, $c_{\mathsf{eq}}^{(ijk3)}$, $c_{\mathsf{lu}}^{(ijk3)}$ and $c_{\mathsf{eu}}^{(ijk3)}$. The coefficient indices represent the lepton flavour generations ($i,j = 1,2,3$ where $i \neq j$) and light quark flavour generations ($k = 1,2$). The single-top-quark production cross-sections are quoted for $u$- and $c$-quark couplings separately, while they are combined for the $t\bar{t}$ decay process ($q_k = u,c$). The scale and PDF uncertainties are given. The value of each Wilson coefficient is set to 1.0 for the calculation of the cross-section.

Requirements for each analysis region. The symbol $\ell_{3}$ denotes the lowest $p_{\mathrm{T}}$ lepton. In $\mathrm{CR}t\bar{t}\mu$ an additional requirement is placed on the leading muon $p_{\mathrm{T}}$ ($p^{\mu_1}_\mathrm{T}$) and dilepton invariant masses in order to reject signal events.

Post-fit event yields for each analysis region entering the fit, with correlations on the full systematic uncertainties taken into account as determined in the fit under a signal+background hypothesis. The 'fake electron' category in the CR$t\bar{t}\mu$ control region collects small contributions primarily from $t\bar{t}\gamma$ and $Z\gamma$ which enter the event selection due to photon conversions.

Expected and observed 95$\%$ CL upper limits on Wilson coefficients corresponding to 2Q2L EFT operators which could introduce a cLFV top decay in the $\mu\tau$ channel, and existing limits from JHEP04 (2019) 014 (previous). The previous limits shown for $c^{1(ijk3)}_{\mathsf{lequ}}$ and $c^{3(ijk3)}_{\mathsf{lequ}}$ are tightened by a factor of $\sqrt{2}$, as JHEP04 (2019) 014 does not assume that these operators are Hermitian. Results are shown separately for the $\mu\tau ut$ and $\mu\tau ct$ interactions. The lepton generations are denoted by $i,j=2,3$ for $\mu$ and $\tau$ (where $i \neq j$) and the quark generations are denoted by $k=1,2$ for $u$ and $c$, respectively.

Expected and observed 95$\%$ CL upper limits on the branching ratio ($\mathcal{B}$) corresponding to the decay of a top quark to a muon and a $\tau$ lepton through a cLFV process using specific Wilson coefficients corresponding to 2Q2L EFT operators. Results are shown separately for $\mu\tau ut$ and $\mu\tau ct$ interactions. The lepton generations are denoted by $i,j=2,3$ for $\mu$ and $\tau$ (where $i \neq j$) and the quark generations are denoted by $k=1,2$ for $u$ and $c$, respectively.

Expected and observed 95% CL upper limits on the inclusive branching ratio ($\mathcal{B}$) corresponding to the decay of a top quark to a muon and a $\tau$-lepton through a cLFV process. Limits are shown for the statistical uncertainty only and for the full set of statistical and systematic uncertainties.

Expected and observed 95$\%$ CL upper limits on the Wilson coefficients corresponding to EFT operators inducing the decay of a top quark to a muon, a tau lepton and an up quark through a cLFV process ($\mu\tau ut$ interaction).

Expected and observed 95$\%$ CL upper limits on the Wilson coefficients corresponding to EFT operators inducing the decay of a top quark to a muon, a tau lepton and a charm quark through a cLFV process ($\mu\tau ct$ interaction).

Ranking of nuisance parameters by post-fit impact on cLFV signal strength $\mu$ (the parameter of interest, POI). The pre-fit (post-fit) impact of each nuisance parameter is calculated as the difference to the fitted value of cLFV signal strength between the nominal fit and a fit when fixing the corresponding nuisance parameter to $\hat{\theta}\pm\Delta\theta$ ($\hat{\theta}\pm\Delta\hat{\theta}$), where $\hat{\theta}$ is the best-fit value of the nuisance parameter and $\Delta\theta$ ($\Delta\hat{\theta}$) is its pre-fit (post-fit) uncertainty. The pulls of the nuisance parameters are calculated relative to zero. These pulls and their post-fit errors, $\Delta\hat{\theta} / \Delta\theta$, are given in the second column of the table.

Fraction of events where both bosons are longitudinally polarized in the region with $p_T^Z>200$ GeV using two unconstrained parameters.

Numbers of observed and expected events in the 00-enhanced signal regions, before the fit. The total uncertainties are quoted.

Summary of the relative uncertainties in the measured longitudinal-longitudinal fractions $f_{00}$. The uncertainties are reported as percentages.

Numbers of observed and expected events in the 00-enhanced signal regions, after the fit. The total uncertainties are quoted.

$|\Delta Y(l_{W}Z)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 20$ GeV of the TT state.

$|\Delta Y(l_{W}Z)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 40$ GeV of the TT state.

$|\Delta Y(l_{W}Z)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 70$ GeV of the TT state.

$|\Delta Y(l_{W}Z)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 20$ GeV of of the sum of all polarization.

$|\Delta Y(l_{W}Z)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 40$ GeV of of the sum of all polarization.

$|\Delta Y(l_{W}Z)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 70$ GeV of of the sum of all polarization.

$|\Delta Y(WZ)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 20$ GeV of the TT state.

$|\Delta Y(WZ)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 40$ GeV of the TT state.

$|\Delta Y(WZ)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 70$ GeV of the TT state.

$|\Delta Y(WZ)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 20$ GeV of of the sum of all polarization.

$|\Delta Y(WZ)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 40$ GeV of of the sum of all polarization.

$|\Delta Y(WZ)|$ Migration matrix RAZ Region $p_{T}^{WZ}< 70$ GeV of of the sum of all polarization.

The $\mathcal{D}$ value for the unfolded $|\Delta Y(l_{W}Z)|$ distributions of the TT polarization state as a function of the $p_T^{WZ}$ cut value.

The $\mathcal{D}$ value for the unfolded $|\Delta Y(WZ)|$ distributions of the TT polarization state as a function of the $p_T^{WZ}$ cut value.

Measured normalized differential $|\Delta Y(l_{W}Z)|$ cross-section of the TT state in the RAZ Region $p_{T}^{WZ}< 20$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(l_{W}Z)|$ cross-section of the TT state in the RAZ Region $p_{T}^{WZ}< 40$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(l_{W}Z)|$ cross-section of the TT state in the RAZ Region $p_{T}^{WZ}< 70$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(WZ)|$ cross-section of the TT state in the RAZ Region $p_{T}^{WZ}< 20$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(WZ)|$ cross-section of the TT state in the RAZ Region $p_{T}^{WZ}< 40$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(WZ)|$ cross-section of the TT state in the RAZ Region $p_{T}^{WZ}< 70$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(l_{W}Z)|$ cross-section of the sum of polarization states in the RAZ Region $p_{T}^{WZ}< 20$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(l_{W}Z)|$ cross-section of the sum of polarization states in the RAZ Region $p_{T}^{WZ}< 40$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(l_{W}Z)|$ cross-section of the sum of polarization states in the RAZ Region $p_{T}^{WZ}< 70$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(WZ)|$ cross-section of the sum of polarization states in the RAZ Region $p_{T}^{WZ}< 20$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(WZ)|$ cross-section of the sum of polarization states in the RAZ Region $p_{T}^{WZ}< 40$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Measured normalized differential $|\Delta Y(WZ)|$ cross-section of the sum of polarization states in the RAZ Region $p_{T}^{WZ}< 70$ GeV. The total uncertainties are quoted. The last bin covers all events above the lower end of the bin.

Summary of the relative uncertainties in the measured depth of the TT state in the RAZ Region $p_{T}^{WZ}<20$ GeV using $|\Delta Y(l_{W}Z)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the TT state in the RAZ Region $p_{T}^{WZ}<40$ GeV using $|\Delta Y(l_{W}Z)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the TT state in the RAZ Region $p_{T}^{WZ}<70$ GeV using $|\Delta Y(l_{W}Z)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the TT state in the RAZ Region $p_{T}^{WZ}<20$ GeV using $|\Delta Y(WZ)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the TT state in the RAZ Region $p_{T}^{WZ}<40$ GeV using $|\Delta Y(WZ)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the TT state in the RAZ Region $p_{T}^{WZ}<70$ GeV using $|\Delta Y(WZ)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the sum of polarization states in the RAZ Region $p_{T}^{WZ}<20$ GeV using $|\Delta Y(l_{W}Z)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the sum of polarization states in the RAZ Region $p_{T}^{WZ}<40$ GeV using $|\Delta Y(l_{W}Z)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the sum of polarization states in the RAZ Region $p_{T}^{WZ}<70$ GeV using $|\Delta Y(l_{W}Z)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the sum of polarization states in the RAZ Region $p_{T}^{WZ}<20$ GeV using $|\Delta Y(WZ)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the sum of polarization states in the RAZ Region $p_{T}^{WZ}<40$ GeV using $|\Delta Y(WZ)|$. The uncertainties are reported as percentages.

Summary of the relative uncertainties in the measured depth of the sum of polarization states in the RAZ Region $p_{T}^{WZ}<70$ GeV using $|\Delta Y(WZ)|$. The uncertainties are reported as percentages.

Observed and expected 95% CL upper limits on the visible $\tau\nu$ production cross section, $\sigma_{\rm vis} = \sigma(pp \to \tau\nu +X) \cdot \mathcal{A} \cdot \varepsilon$, for $m_{\rm T}$ thresholds ranging from 200 to 2950 GeV. See HepData abstract for details on how to use this data for reinterpretation.

Observed 95% CL upper limits on the cross-section times branching ratio as a function of the W′ mass in several NUGIM models defined by the value of the parameter cot$\theta$.

Expected 95% CL upper limits on the cross-section times branching ratio as a function of the W′ mass in several NUGIM models defined by the value of the parameter cot$\theta$.

Acceptance for $W^{\prime}_{\rm SSM}$ as a function of the $W^{\prime}_{\rm SSM}$ mass, shown after successively applying selection at generator-level. The acceptance times efficiency is calculated with respect to all $W^{\prime}_{\rm SSM} \to \tau\nu$ events with a generated $\tau\nu$ mass above 25 GeV. The "selected tau" criteria include the requirement of a $\tau_{\rm had-vis}$ with $p_{\rm T}$ > 30 GeV and $|\eta|$ < 2.4.

Acceptance times efficiency for $W^{\prime}_{\rm SSM}$ as a function of the $W^{\prime}_{\rm SSM}$ mass, shown after successively applying selection at reconstruction-level. The acceptance times efficiency is calculated with respect to all $W^{\prime}_{\rm SSM} \to \tau\nu$ events with a generated $\tau\nu$ mass above 25 GeV. "Preselection" includes all criteria prior to those shown.

Observed 95% CL exclusion limits on the simplified models of chargino-pair production decaying via W bosons.

$+1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino-pair production decaying via W bosons.

$-1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino-pair production decaying via W bosons.

Expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$+1\sigma$ expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$-1\sigma$ expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

Observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$+1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$-1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

Expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and h bosons.

$+1\sigma$ expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and h bosons.

$-1\sigma$ expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and h bosons.

Observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and h bosons.

$+1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and h bosons.

$-1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and h bosons.

Expected 95% CL exclusion limits on the simplified models of higgsino GGM scenarios.

$+1\sigma$ expected 95% CL exclusion limits on the simplified models of higgsino GGM scenarios.

$-1\sigma$ expected 95% CL exclusion limits on the simplified models of higgsino GGM scenarios.

Observed 95% CL exclusion limits on the simplified models of higgsino GGM scenarios.

$+1\sigma$ observed 95% CL exclusion limits on the simplified models of higgsino GGM scenarios.

$-1\sigma$ observed 95% CL exclusion limits on the simplified models of higgsino GGM scenarios.

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

The analyses used in combination for each scenario to set limits in models of the production of $\tilde{\chi}_1^{+}\tilde{\chi}_{1}^{-}$.

Observed upper limit on the signal cross section in fb for chargino--neutralino production decaying via W and Z bosons.

The analyses used in combination for each scenario to set limits in models of chargino--neutralino production decaying via W and Z bosons.

Expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$+1\sigma$ expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$-1\sigma$ expected 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

Observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$+1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

$-1\sigma$ observed 95% CL exclusion limits on the simplified models of chargino--neutralino production decaying via W and Z bosons.

Observed upper limit on the signal cross section in fb for chargino--neutralino production decaying via W and h bosons.

The analyses used in combination for each scenario to set limits in models of chargino--neutralino production decaying via W and h bosons.

Observed upper limit on the signal cross section in fb for higgsino GGM scenarios.

The analyses used in combination for each scenario to set limits in higgsino GGM scenarios.

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.

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.