Data from Fig. 3 lower right panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 13$ TeV as a function of the jet multiplicity with $|\eta_{j}| < 2.4$.

Data from Fig. 4 upper right panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 13$ TeV as a function of the p$_T$-leading jet transverse momentum with $|\eta_{j}| < 4.7$.

Data from Fig. 5 left panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 13$ TeV as a function of the p$_T$-leading jet transverse momentum with $|\eta_{j}| < 4.7$.

Data from Fig. 4 lower right panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 13$ TeV as a function of the p$_T$-leading jet $|\eta_{j}|$ with $|\eta_{j}| < 4.7$.

Data from Fig. 5 right panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 13$ TeV as a function of the p$_T$-subleading jet $|\eta_{j}|$ with $|\eta_{j}| < 4.7$.

Data from Fig. 6 left panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 13$ TeV as a function of the invariant mass of the two p$_T$ leading jets with $|\eta_{j}| < 4.7$.

Data from Fig. 6 right panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 13$ TeV as a function of the invariant of pseudorapidity separation of the two $p$_T$ leading jets with $|\eta_{j}| < 4.7$.

Data from Fig. 2 upper left. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ differential cross section at $\sqrt{s} = 8$ TeV as a function of the jet multiplicity with $|\eta| < 4.7$.

Data from Fig. 3 upper left panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 8$ TeV as a function of the jet multiplicity with $|\eta| < 4.7$.

Data from Fig. 2 lower left panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ differential cross section at $\sqrt{s} = 8$ TeV as a function of the jet multiplicity with $|\eta| < 2.4$.

Data from Fig. 3 lower left panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 8$ TeV as a function of the jet multiplicity with $|\eta| < 2.4$.

Data from Fig. 4 upper left panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 8$ TeV as a function of the p$_T$-leading jet transverse momentum with $|\eta| < 4.7$.

Data from Fig. 4 lower left panel. The $\textrm{pp} \to \textrm{ZZ}\to \ell\ell\ell^{\prime}\ell^{\prime}$ normalized differential cross section at $\sqrt{s} = 8$ TeV as a function of the p$_T$-leading jet transverse $\eta$ with $|\eta| < 4.7$.

Comparison of the distribution of the scalar sum of small-$R$ jet and lepton transverse momenta, $S_T$, between data and the background prediction in the signal region of the pair-production (PP) $\geq 3\ell$ channel. The background prediction is shown post-fit, i.e. after the fit to the data $S_T$ distributions under the background-only hypothesis. The last bin contains the overflow. An example distribution for a $B\bar B$ signal in the singlet model with $m_B$ = 900 GeV is overlaid.

Comparison of the distribution of the invariant mass of the $Z$ boson candidate and the highest-$p_T$ top-tagged large-$R$ jet, $m_{Zt}$, between data and the background prediction in the signal region of the single-production (SP) 2$\ell$ channel. The background prediction is shown post-fit, i.e. after the fit to the data $m_{Zt}$ distributions under the background-only hypothesis. The last bin contains the overflow. An example distribution for a single-$T$-quark signal with $m_T$ = 900 GeV and $\kappa_T = 0.5$ is overlaid.

Comparison of the distribution of the scalar sum of small-$R$ jet and lepton transverse momenta, $S_T$, between data and the background prediction in the signal region of the single-production (SP) $\geq 3\ell$ channel. The background prediction is shown post-fit, i.e. after the fit to the data $S_T$ distributions under the background-only hypothesis. The last bin contains the overflow. An example distribution for a single-$T$-quark signal with $m_T$ = 900 GeV and $\kappa_T$ = 0.5 is overlaid.

Upper limits at 95$\%$ CL on the cross section of vector-like quark pair production (PP) for $T\bar T$ in the singlet model. The expected limits are shown for the combination of the channels, as are the observed limits.

Upper limits at 95$\%$ CL on the cross section of vector-like quark pair production (PP) for $B\bar B$ in the singlet model. The expected limits are shown for the combination of the channels, as are the observed limits.

Upper limits at 95$\%$ CL on the cross section of vector-like quark pair production (PP) for $T\bar T$ in the doublet model. The expected limits are shown for the combination of the channels, as are the observed limits.

Upper limits at 95$\%$ CL on the cross section of vector-like quark pair production (PP) for $B\bar B$ in the doublet model. The expected limits are shown for the combination of the channels, as are the observed limits.

Upper limits at 95$\%$ CL on the cross section of vector-like quark pair production (PP) for $T\bar T$ with a BR of 100$\%$ to $Zt$. The expected limits are shown for the combination of the channels, as are the observed limits.

Upper limits at 95$\%$ CL on the cross section of vector-like quark pair production (PP) for $B\bar B$ with a BR of 100$\%$ to $Zb$. The expected limits are shown for the combination of the channels, as are the observed limits.

Expected 95$\%$ CL lower limits from the combination of the pair-production channels on the mass of vector-like quarks for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$, adding up to unity.

Observed 95$\%$ CL lower limits from the combination of the pair-production channels on the mass of vector-like quarks for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$, adding up to unity.

Expected 95$\%$ CL lower limits from the combination of the pair-production channels on the mass of vector-like quarks for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$, adding up to unity.

Observed 95$\%$ CL lower limits from the combination of the pair-production channels on the mass of vector-like quarks for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$, adding up to unity.

Upper limits at 95$\%$ CL on the cross section times BR to $Zt$ of single production (SP) of a $T$-quark. The expected limits are shown for the combination of the channels, as are the observed limits. The expected cross section times BR to $Zt$ for single-$T$-quark production is also shown for a coupling $\kappa_T = 0.5$, which corresponds to a coupling of $c_W = \sqrt{c^2_{W,L} + c^2_{W,R}} = 0.45$ from Ref. [14]. The BR assumed here corresponds to the singlet benchmark model, i.e. $\approx 25\%$.

Expected and observed 95$\%$ CL limits from the combination of the single-production channels on the coupling of the $T$ quark to SM particles, $c_W = \sqrt c^2_L + c^2_R$, assuming the singlet-model BR of $\approx 25\%$, as a function of the mass of the $T$ quark, $m_T$. Values of $c_W$ larger than the observed limit are excluded.

Expected and observed 95$\%$ CL limits from the combination of the single-production channels on the mixing angle in the singlet model between the $T$ quark and the top quark, $|\sin\theta_L|$, as a function of the mass of the $T$ quark, $m_T$. Values of $|\sin\theta_L|$ enclosed by the observed limit are excluded, i.e. for $m_T$ larger than $\approx 1200$ GeV, no value of $|\sin\theta_L|$ is excluded.

Expected lower limit from the combination of the single-production channels on the mass of the $T$ quark as a function of the couplings of the $T$ quark to the $W$ boson, $\sqrt c^2_L + c^2_R$, and to the $Z$ boson, $\sqrt c^2_L + c^2_R$ with the assumption of equal BRs for $T\rightarrow Zt$ and $T\rightarrow Ht$ in the limit of large $T$-quark masses.

Observed lower limit from the combination of the single-production channels on the mass of the $T$ quark as a function of the couplings of the $T$ quark to the $W$ boson, $\sqrt c^2_L + c^2_R$, and to the $Z$ boson, $\sqrt c^2_L + c^2_R$ with the assumption of equal BRs for $T\rightarrow Zt$ and $T\rightarrow Ht$ in the limit of large $T$-quark masses.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 500 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 600 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 700 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 750 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 800 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 850 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 900 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 950 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 1000 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 1050 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 1100 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 1150 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 1200 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 1300 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow T\bar T)$ for the combination of the pair-production channels for all combinations of BRs for $T\rightarrow Zt$, $T\rightarrow Ht$, $T\rightarrow Wb$ adding up to unity for a $T$-quark mass of 1400 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 500 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 600 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 700 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 750 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 800 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 850 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 900 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 950 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 1000 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 1050 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 1100 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 1150 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 1200 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 1300 GeV.

Observed upper limits at 95$\%$ CL on $\sigma(pp\rightarrow B\bar B)$ for the combination of the pair-production channels for all combinations of BRs for $B\rightarrow Zb$, $B\rightarrow Hb$, $B\rightarrow Wt$ adding up to unity for a $B$-quark mass of 1400 GeV.

Observed upper limits at 95$\%$ CL on the cross section times BR to $Zt$ as a function of the vector-like $T$-quark mass and coupling $\kappa_T$ for the combination of the SP $2\ell$ and SP $\geq 3\ell$ channels.

Signal cutflows for the $B\bar B$ process in the singlet model in the PP $2\ell$ $0-1$J channel in the 0-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the singlet model in the PP $2\ell$ $0-1$J channel in the 1-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the doublet model in the PP $2\ell$ $0-1$J channel in the 0-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the doublet model in the PP $2\ell$ $0-1$J channel in the 1-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the case of 100$\%$ BR for $B \rightarrow Zb$ in the PP $2\ell$ $0-1$J channel in the 0-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the case of 100$\%$ BR for $B \rightarrow Zb$ in the PP $2\ell$ $0-1$J channel in the 1-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the singlet model in the PP $2\ell$ $0-1$J channel in the 0-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the singlet model in the PP $2\ell$ $0-1$J channel in the 1-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the doublet model in the PP $2\ell$ $0-1$J channel in the 0-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the doublet model in the PP $2\ell$ $0-1$J channel in the 1-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the case of 100$\%$ BR for $T \rightarrow Zt$ in the PP $2\ell$ $0-1$J channel in the 0-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the case of 100$\%$ BR for $T \rightarrow Zt$ in the PP $2\ell$ $0-1$J channel in the 1-large-$R$ jet SR. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the singlet model in the PP $2\ell$ $\geq 2$J channel. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the singlet model in the PP $2\ell$ $\geq 2$J channel. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the doublet model in the PP $2\ell$ $\geq 2$J channel. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the doublet model in the PP $2\ell$ $\geq 2$J channel. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the case of 100$\%$ BR for $T \rightarrow Zt$ in the PP $2\ell$ $\geq 2$J channel. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the case of 100$\%$ BR for $B \rightarrow Zb$ in the PP $2\ell$ $\geq 2$J channel. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the singlet model in the PP $\geq 3\ell$ channel. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the singlet model in the PP $\geq 3\ell$ channel. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the doublet model in the PP $\geq 3\ell$ channel. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the doublet model in the PP $\geq 3\ell$ channel. Only statistical uncertainties are shown.

Signal cutflows for the $T\bar T$ process in the case of 100$\%$ BR for $T \rightarrow Zt$ in the PP $\geq 3\ell$ channel. Only statistical uncertainties are shown.

Signal cutflows for the $B\bar B$ process in the case of 100$\%$ BR for $B \rightarrow Zb$ in the PP $\geq 3\ell$ channel. Only statistical uncertainties are shown.

Signal cutflows for single-$T$-quark production with $\kappa_T = 0.5$ in the SP $2\ell$ channel. Only statistical uncertainties are shown.

Signal cutflows for single-$T$-quark production with $\kappa_T = 0.5$ in the SP $\geq 3\ell$ channel. Only statistical uncertainties are shown.

Signal efficiencies in $\%$ in the PP $2\ell$ $0-1$J channel in the 0-large-$R$ jet SR. Uncertainties are statistical only.

Signal efficiencies in $\%$ in the PP $2\ell$ $0-1$J channel in the 1-large-$R$ jet SR. Uncertainties are statistical only.

Signal efficiencies in $\%$ in the PP $2\ell$ $\geq 2$J channel. Uncertainties are statistical only.

Signal efficiencies in $\%$ in the PP $\geq 3\ell$ channel. Uncertainties are statistical only.

Signal efficiencies in $\%$ in the SP $2\ell$ channel for $\kappa_T = 0.5$. Uncertainties are statistical only.

Signal efficiencies in $\%$ in the SP $\geq 3\ell$ channel for $\kappa_T = 0.5$. Uncertainties are statistical only.

The systematic uncertainty on the unfolded distribution as a function of minimax-mbl, broken down by components.

The covariance matrix for the unfolded, normalized data in bins of minimax-mbl shown in Figure 2. Redundant entries are suppressed.

Leading lepton transverse momentum distribution for events passing the signal selection in bins of minimax-mbl. The data here are not unfolded.

Sub-leading lepton transverse momentum distribution for events passing the signal selection in bins of minimax-mbl. The data here are not unfolded.

Leading lepton pseudorapidity distribution for events passing the signal selection in bins of minimax-mbl. The data here are not unfolded.

Sub-leading lepton pseudorapidity distribution for events passing the signal selection in bins of minimax-mbl. The data here are not unfolded.

Leading b-jet transverse momentum distribution for events passing the signal selection in bins of minimax-mbl. The data here are not unfolded.

Sub-leading b-jet transverse momentum distribution for events passing the signal selection in bins of minimax-mbl. The data here are not unfolded.

Leading b-jet pseudorapidity distribution for events passing the signal selection in bins of minimax-mbl. The data here are not unfolded.

Sub-leading b-jet pseudorapidity distribution for events passing the signal selection in bins of minimax-mbl. The data here are not unfolded.

The detector-level minimax-m(bl) distribution for events entering the Z+jets control region.

The detector-level minimax-m(bl) distribution for events entering the same-sign lepton control region.

Exclusion limits on higgsino pair production divided by the theory cross-section.The results of the low-mass analysis are used below m(higgsino) = 300 GeV, while those of the high-mass analysis are used above. The figure shows the observed and expected 95% upper limits on the higgsino pair production cross-section as a function of m(higgsino).

Observed and expected 95% limits in the m(higgsino) vs BR(higgsino to higgs+gravitino) plane. The regions above the lines are excluded by the analyses.

The observed and expected 95% upper limits on the total pair production cross section for degenerate higgsinos as a function of m(higgsino) for the high-mass search. Only the high-mass analysis results are used in this figure.

The observed and expected 95% upper limits on the total pair production cross section for degenerate higgsinos as a function of m(higgsino) for the high-mass search, divided by the theory cross section. Only the high-mass analysis results are used in this figure.

The observed and expected 95% upper limits on the total pair production cross section for degenerate higgsinos as a function of m(higgsino) for the low-mass search. Only the low-mass analysis results are used in this figure.

The observed and expected 95% upper limits on the total pair production cross section for degenerate higgsinos as a function of m(higgsino) for the low-mass search, divided by the theory cross section. Only the low-mass analysis results are used in this figure.

Particle-level acceptance for the low-mass discovery signal regions low-SR-MET0-meff440 and low-SR-MET150-meff440, shown as a function of higgsino mass. The acceptance is defined as the fraction of signal events passing the particle-level event selection that emulates the detector-level selection. The acceptance calculation considers only those signal events where both higgsinos decay to Higgs bosons that subsequently both decay to b-quarks.

The experimental efficiency of the low-mass analysis, for the two discovery signal regions low-SR-MET0-meff440 and low-SR-MET150-meff440, as a function of higgsino mass. The experimental efficiency is defined as the number of events passing the detector-level event selection divided by the number of events passing the event selection for a perfect detector. The denominator is obtained by implementing particle-level event selection that emulate the detector-level selection. Such particle-level selection is not applied on the numerator.

Particle-level acceptance for the high-mass discovery signal regions SR-4b-meff1-A-disc and SR-3b-meff3-A, shown as a function of higgsino mass. The acceptance is defined as the fraction of signal events passing the particle-level event selection that emulates the detector-level selection. The acceptance calculation considers only those signal events where both higgsinos decay to Higgs bosons that subsequently both decay to b-quarks.

The experimental efficiency of the high-mass analysis, for the two discovery signal regions SR-4b-meff1-A-disc and SR-3b-meff3-A, as a function of higgsino mass. The experimental efficiency is defined as the number of events passing the detector-level event selection divided by the number of events passing the event selection for a perfect detector. The denominator is obtained by implementing particle-level event selection that emulate the detector-level selection. Such particle-level selection is not applied on the numerator.

Example cutflow for SR-3b-meff3-A.

Example cutflow for SR-4b-meff1-A-disc.

Cutflow for low-mass analysis for each signal mass point.

Observed ST distribution in the tau-tau signal region, compared to the expected SM background contributions. The distribution labeled electroweak contains the contributions from W+jets, Z+jets, and diboson processes. The signal distributions for single-leptoquark (LQ) production with mass 700 GeV are overlaid to illustrate the sensitivity. For the signal normalization, lambda = 1 and beta = 1 are assumed. The background uncertainty bands represent the sum in quadrature of statistical and systematic uncertainties obtained from the fit. The lower panels show the ratio between the observed and expected events in each bin. In all plots, the horizontal and vertical error bars on the data points represent the bin widths and the Poisson uncertainties, respectively.

Observed ST distribution in the e-mu control region, compared to the expected SM background contributions. The distribution labeled electroweak contains the contributions from W+jets, Z+jets, and diboson processes. The signal distributions for single-leptoquark (LQ) production with mass 700 GeV are overlaid to illustrate the sensitivity. For the signal normalization, lambda = 1 and beta = 1 are assumed. The background uncertainty bands represent the sum in quadrature of statistical and systematic uncertainties obtained from the fit. The lower panels show the ratio between the observed and expected events in each bin. In all plots, the horizontal and vertical error bars on the data points represent the bin widths and the Poisson uncertainties, respectively.

The covariance matrix of the bin contents of the background fit. em stands for e-mu channel, et for e-tau channel, mt for mu-tau channel, and tt for tau-tau channel. The numbers indicate the bin number in each final state.

Observed (black solid) and expected (black dotted) limits at 95% confidence level on the product of cross section (sigma) and branching fraction beta, obtained from the combination of the e-tau, mu-tau, and tau-tau signal regions, as well as from the emu control region, as a function of the leptoquark (LQ) mass. The green and yellow bands represent the one and two standard deviation uncertainties in the expected limits. The theory prediction is indicated by the blue solid line, together with systematic uncertainties due to the choice of PDF and renormalization and factorization scales, indicated by the blue band. For generator-level cuts, please refer to the resource.

Expected and observed exclusion limits at 95% confidence level on the Yukawa coupling (lambda) at the leptoquark(LQ)-lepton-quark vertex, as a function of the LQ mass. A unit branching fraction beta of the LQ to a tau lepton and a bottom quark is assumed. The orange vertical line indicates the limit obtained from a search for pair-produced LQs decaying to lepton-tau-bb. The left-hand side of the dotted (solid) line shows the expected (observed) exclusion region for the present analysis. The gray band shows +/- 1 standard deviations of the expected exclusion limit. The region with diagonal blue shading shows the parameter space preferred by one of the models proposed to explain anomalies observed in B physics.

Distributions of kinematic variables in the signal regions for the $2\ell$ channels after applying all selection requirements. The histograms show the post-fit background predictions. The last bin includes the overflow. The distribution for $R_{\textrm{ISR}}$ in SR$2\ell$_ISR is plotted. The expected distribution for a benchmark signal model, normalized to the NLO+NLL cross-section times integrated luminosity, is also shown for comparison.

Distributions of kinematic variables in the signal regions for the $3\ell$ channels after applying all selection requirements. The histograms show the post-fit background predictions. The last bin includes the overflow. The distribution for $H_{3,1}^{\textrm{PP}}$ in SR$3\ell$_Low is plotted. The expected distribution for a benchmark signal model, normalized to the NLO+NLL cross-section times integrated luminosity, is also shown for comparison.

Distributions of kinematic variables in the signal regions for the $3\ell$ channels after applying all selection requirements. The histograms show the post-fit background predictions. The last bin includes the overflow. The distribution of $p_{\textrm{T}}^{\ell_{1}}$ in SR$3\ell$_Low is plotted. The expected distribution for a benchmark signal model, normalized to the NLO+NLL cross-section times integrated luminosity, is also shown for comparison.

Distributions of kinematic variables in the signal regions for the $3\ell$ channels after applying all selection requirements. The histograms show the post-fit background predictions. The last bin includes the overflow. The distribution of $p_{\mathrm{T\ ISR}}^{~\textrm{CM}}$ in SR$3\ell$_ISR is plotted. The expected distribution for a benchmark signal model, normalized to the NLO+NLL cross-section times integrated luminosity, is also shown for comparison.

Distributions of kinematic variables in the signal regions for the $3\ell$ channels after applying all selection requirements. The histograms show the post-fit background predictions. The last bin includes the overflow. The distribution of $R_{\textrm{ISR}}$ in SR$3\ell$_ISR is plotted. The expected distribution for a benchmark signal model, normalized to the NLO+NLL cross-section times integrated luminosity, is also shown for comparison.

Expected exclusion limits at 95% CL on the masses of C1/N2 and N1 from the analysis of 36.1fb$^{-1}$ of 13 TeV $pp$ collision data obtained from the 2$\ell$ search, assuming 100\% branching ratio of the sparticles to decay to SM $W/Z$ bosons and LSP.

Observed exclusion limits at 95% CL on the masses of C1/N2 and N1 from the analysis of 36.1fb$^{-1}$ of 13 TeV $pp$ collision data obtained from the 2$\ell$ search, assuming 100% branching ratio of the sparticles to decay to SM $W/Z$ bosons and LSP.

Expected exclusion limits at 95% CL on the masses of C1/N2 and N1 from the analysis of 36.1fb$^{-1}$ of 13 TeV $pp$ collision data obtained from the the $3\ell$ search, assuming 100% branching ratio of the sparticles to decay to SM $W/Z$ bosons and LSP.

Observed exclusion limits at 95% CL on the masses of C1/N2 and N1 from the analysis of 36.1fb$^{-1}$ of 13 TeV $pp$ collision data obtained from the the $3\ell$ search, assuming 100% branching ratio of the sparticles to decay to SM $W/Z$ bosons and LSP.

Expected exclusion limits at 95% CL on the masses of C1/N2 and N1 from the analysis of 36.1fb$^{-1}$ of 13 TeV $pp$ collision data obtained from the statistical combination of the $2\ell$ and 3$\ell$ search channels, assuming 100\% branching ratio of the sparticles to decay to SM $W/Z$ bosons and LSP.

Observed exclusion limits at 95% CL on the masses of C1/N2 and N1 from the analysis of 36.1fb$^{-1}$ of 13 TeV $pp$ collision data obtained from the statistical combination of the $2\ell$ and 3$\ell$ search channels, assuming 100% branching ratio of the sparticles to decay to SM $W/Z$ bosons and LSP.

Signal region acceptance for chargino-neutralino production in SR2L_High. The acceptance is with respect to all possible $W/Z$ decays.

Signal region efficiency for chargino-neutralino production in SR2L_High.

Signal region acceptanceXefficiency for chargino-neutralino production in SR2L_High.

Signal region acceptance for chargino-neutralino production in SR2L_Int. The acceptance is with respect to all possible $W/Z$ decays.

Signal region efficiency for chargino-neutralino production in SR2L_Int.

Signal region acceptanceXefficiency for chargino-neutralino production in SR2L_Int.

Signal region acceptance for chargino-neutralino production in SR2L_Low. The acceptance is with respect to all possible $W/Z$ decays.

Signal region efficiency for chargino-neutralino production in SR2L_Low.

Signal region acceptanceXefficiency for chargino-neutralino production in SR2L_Low.

Signal region acceptance for chargino-neutralino production in SR2L_ISR. The acceptance is with respect to all possible $W/Z$ decays.

Signal region efficiency for chargino-neutralino production in SR2L_ISR.

Signal region acceptanceXefficiency for chargino-neutralino production in SR2L_ISR.

Signal region acceptance for chargino-neutralino production in SR3L_High. The acceptance is with respect to all possible $W/Z$ decays.

Signal region efficiency for chargino-neutralino production in SR3L_High.

Signal region acceptanceXefficiency for chargino-neutralino production in SR3L_High.

Signal region acceptance for chargino-neutralino production in SR3L_Int. The acceptance is with respect to all possible $W/Z$ decays.

Signal region efficiency for chargino-neutralino production in SR3L_Int.

Signal region acceptanceXefficiency for chargino-neutralino production in SR3L_Int.

Signal region acceptance for chargino-neutralino production in SR3L_Low. The acceptance is with respect to all possible $W/Z$ decays.

Signal region efficiency for chargino-neutralino production in SR3L_Low.

Signal region acceptanceXefficiency for chargino-neutralino production in SR3L_Low.

Signal region acceptance for chargino-neutralino production in SR3L_ISR. The acceptance is with respect to all possible $W/Z$ decays.

Signal region efficiency for chargino-neutralino production in SR3L_ISR.

Signal region acceptanceXefficiency for chargino-neutralino production in SR3L_ISR.

Observed 95% CL upper limit on the production cross-section of C1/N2 from the analysis of 36.1fb$^{-1}$ of 13 TeV $pp$ collision data obtained from the 2$\ell$ search, assuming 100% branching ratio of the sparticles to decay to SM $W/Z$ bosons and N1.

Observed 95% CL upper limit on the production cross-section of C1/N2 from the analysis of 36.1fb$^{-1}$ of 13 TeV $pp$ collision data obtained from the 3$\ell$ search, assuming 100% branching ratio of the sparticles to decay to SM $W/Z$ bosons and N1.

Observed 95% CL upper limit on the production cross-section of C1/N2 from the analysis of 36.1fb$^{-1}$ of 13 TeV $pp$ collision data obtained from the statistical combination of 2 and 3$\ell$ searches, assuming 100% branching ratio of the sparticles to decay to SM $W/Z$ bosons and N1.

Signal cutflow for SR2L_High and m[C1/N2,N1] = [500,0] GeV. 5000 events were generated for this point. The unweighted number of events correspond to the number of selected MC events in the produced sample, while the weighted yield is normalized to the luminosity of the data.

Signal cutflow for SR2L_Int and m[C1/N2,N1] = [400,200] GeV. 10000 events were generated for this point. The unweighted number of events correspond to the number of selected MC events in the produced sample, while the weighted yield is normalized to the luminosity of the data.

Signal cutflow for SR2L_Low and m[C1/N2,N1] = [200,100] GeV. 20000 events were generated for this point. The unweighted number of events correspond to the number of selected MC events in the produced sample, while the weighted yield is normalized to the luminosity of the data.

Signal cutflow for SR2L_ISR and m[C1/N2,N1] = [200,100] GeV. 20000 events were generated for this point. The unweighted number of events correspond to the number of selected MC events in the produced sample, while the weighted yield is normalized to the luminosity of the data.

Signal cutflow for SR3L_High and m[C1/N2,N1] = [500,0] GeV. 5000 events were generated for this point. The unweighted number of events correspond to the number of selected MC events in the produced sample, while the weighted yield is normalized to the luminosity of the data.

Signal cutflow for SR3L_Int and m[C1/N2,N1] = [400,200] GeV. 10000 events were generated for this point. The unweighted number of events correspond to the number of selected MC events in the produced sample, while the weighted yield is normalized to the luminosity of the data.

Signal cutflow for SR3L_Low and m[C1/N2,N1] = [200,100] GeV. 20000 events were generated for this point. The unweighted number of events correspond to the number of selected MC events in the produced sample, while the weighted yield is normalized to the luminosity of the data.

Signal cutflow for SR3L_ISR and m[C1/N2,N1] = [200,100] GeV. 20000 events were generated for this point. The unweighted number of events correspond to the number of selected MC events in the produced sample, while the weighted yield is normalized to the luminosity of the data.

The centrality dependence of e+e− invariant mass spectra within the STAR acceptance from Au+Au collisions and U+U collisions for pair pT < 0.15 GeV/c. The vertical bars on data points depict the statistical uncertainties, while the systematic uncertainties are shown as gray boxes. The hadronic cocktail yields from U+U collisions are ∼5%–12% higher than those from Au+Au collisions in given centrality bins; thus only cocktails for Au+Au collisions are shown here as solid lines, with shaded bands representing the systematic uncertainties for clarity.

The centrality dependence of e+e− invariant mass spectra within the STAR acceptance from Au+Au collisions and U+U collisions for pair pT < 0.15 GeV/c. The vertical bars on data points depict the statistical uncertainties, while the systematic uncertainties are shown as gray boxes. The hadronic cocktail yields from U+U collisions are ∼5%–12% higher than those from Au+Au collisions in given centrality bins; thus only cocktails for Au+Au collisions are shown here as solid lines, with shaded bands representing the systematic uncertainties for clarity.

The centrality dependence of e+e− invariant mass spectra within the STAR acceptance from Au+Au collisions and U+U collisions for pair pT < 0.15 GeV/c. The vertical bars on data points depict the statistical uncertainties, while the systematic uncertainties are shown as gray boxes. The hadronic cocktail yields from U+U collisions are ∼5%–12% higher than those from Au+Au collisions in given centrality bins; thus only cocktails for Au+Au collisions are shown here as solid lines, with shaded bands representing the systematic uncertainties for clarity.

The corresponding ratios of data over cocktail for Figure 1.

The corresponding ratios of data over cocktail for Figure 1.

The corresponding ratios of data over cocktail for Figure 1.

The corresponding ratios of data over cocktail for Figure 1.

The corresponding ratios of data over cocktail for Figure 1.

The corresponding ratios of data over cocktail for Figure 1.

The e+e− pair pT distributions within the STAR acceptance for different mass regions in 60%–80% Au+Au and U+U collisions compared to cocktails. The systematic uncertainties of the data are shown as gray boxes. The gray bands depict the systematic uncertainties of the cocktails.

The e+e− pair pT distributions within the STAR acceptance for different mass regions in 60%–80% Au+Au and U+U collisions compared to cocktails. The systematic uncertainties of the data are shown as gray boxes. The gray bands depict the systematic uncertainties of the cocktails.

The e+e− pair pT distributions within the STAR acceptance for different mass regions in 60%–80% Au+Au and U+U collisions compared to cocktails. The systematic uncertainties of the data are shown as gray boxes. The gray bands depict the systematic uncertainties of the cocktails.

The e+e− pair pT distributions within the STAR acceptance for different mass regions in 60%–80% Au+Au and U+U collisions compared to cocktails. The systematic uncertainties of the data are shown as gray boxes. The gray bands depict the systematic uncertainties of the cocktails.

The e+e− pair pT distributions within the STAR acceptance for different mass regions in 60%–80% Au+Au and U+U collisions compared to cocktails. The systematic uncertainties of the data are shown as gray boxes. The gray bands depict the systematic uncertainties of the cocktails.

The e+e− pair pT distributions within the STAR acceptance for different mass regions in 60%–80% Au+Au and U+U collisions compared to cocktails. The systematic uncertainties of the data are shown as gray boxes. The gray bands depict the systematic uncertainties of the cocktails.

The low-pT (pT < 0.15 GeV/c) e+e− excess mass spectra (data − cocktail) within the STAR acceptance in 60%–80% for Au + Au and U + U collisions, compared with a broadened rho meson model calculation [8]. The contributions of rho and phi from the photonuclear process are shown, as are the contributions of photon-photon process from two models [33,34]. The model calculations are for Au + Au collisions in the corresponding centrality bins.

The low-pT (pT < 0.15 GeV/c) e+e− excess mass spectra (data − cocktail) within the STAR acceptance in 60%–80% for Au + Au and U + U collisions, compared with a broadened rho meson model calculation [8]. The contributions of rho and phi from the photonuclear process are shown, as are the contributions of photon-photon process from two models [33,34]. The model calculations are for Au + Au collisions in the corresponding centrality bins.

The low-pT (pT < 0.15 GeV/c) e+e− excess mass spectra (data − cocktail) within the STAR acceptance in 40%–60% for Au + Au and U + U collisions, compared with a broadened rho meson model calculation [8]. The contributions of rho and phi from the photonuclear process are shown, as are the contributions of photon-photon process from two models [33,34]. The model calculations are for Au + Au collisions in the corresponding centrality bins.

The low-pT (pT < 0.15 GeV/c) e+e− excess mass spectra (data − cocktail) within the STAR acceptance in 40%–60% for Au + Au and U + U collisions, compared with a broadened rho meson model calculation [8]. The contributions of rho and phi from the photonuclear process are shown, as are the contributions of photon-photon process from two models [33,34]. The model calculations are for Au + Au collisions in the corresponding centrality bins.

The centrality dependence of integrated excess yields in the mass regions of 0.4–0.76, 0.76–1.2, and 1.2–2.6 GeV=c2 in Au + Au and U + U collisions. The centrality dependence of hadronic cocktail yields in the mass region of 0.76–1.2 GeV/c2 in both collisions is also shown for comparison. The systematic uncertainties are shown as gray boxes.

The centrality dependence of integrated excess yields in the mass regions of 0.4–0.76, 0.76–1.2, and 1.2–2.6 GeV=c2 in Au + Au and U + U collisions. The centrality dependence of hadronic cocktail yields in the mass region of 0.76–1.2 GeV/c2 in both collisions is also shown for comparison. The systematic uncertainties are shown as gray boxes.

The centrality dependence of integrated excess yields in the mass regions of 0.4–0.76, 0.76–1.2, and 1.2–2.6 GeV=c2 in Au + Au and U + U collisions. The centrality dependence of hadronic cocktail yields in the mass region of 0.76–1.2 GeV/c2 in both collisions is also shown for comparison. The systematic uncertainties are shown as gray boxes.

The centrality dependence of integrated excess yields in the mass regions of 0.4–0.76, 0.76–1.2, and 1.2–2.6 GeV=c2 in Au + Au and U + U collisions. The centrality dependence of hadronic cocktail yields in the mass region of 0.76–1.2 GeV/c2 in both collisions is also shown for comparison. The systematic uncertainties are shown as gray boxes.

The centrality dependence of integrated excess yields in the mass regions of 0.4–0.76, 0.76–1.2, and 1.2–2.6 GeV=c2 in Au + Au and U + U collisions. The centrality dependence of hadronic cocktail yields in the mass region of 0.76–1.2 GeV/c2 in both collisions is also shown for comparison. The systematic uncertainties are shown as gray boxes.

The centrality dependence of integrated excess yields in the mass regions of 0.4–0.76, 0.76–1.2, and 1.2–2.6 GeV/c2 in Au + Au and U + U collisions. The centrality dependence of hadronic cocktail yields in the mass region of 0.76–1.2 GeV/c2 in both collisions is also shown for comparison. The systematic uncertainties are shown as gray boxes.

The p2T distributions of excess yields within the STAR acceptance in the mass regions of 0.4–0.76 GeV/c2 in 60%–80% Au + Au and U + U collisions. The systematic uncertainties are shown as gray boxes. The solid and dotted lines are exponential fits to the data in Au + Au and U + U collisions, respectively. The dot-dashed and dot-dot- dashed lines represent the p2T distributions for the photon-photon process from two models [33,34] within the STAR acceptance in 60%–80% Au + Au collisions. The dashed lines illustrate the corresponding p2T distributions for e+e− pairs from the model [33] traversing 1 fm in a constant magnetic field of 10^14 T perpendicular to the beam line.

The p2T distributions of excess yields within the STAR acceptance in the mass regions of 0.4–0.76 GeV/c2 in 60%–80% Au + Au and U + U collisions. The systematic uncertainties are shown as gray boxes. The solid and dotted lines are exponential fits to the data in Au + Au and U + U collisions, respectively. The dot-dashed and dot-dot- dashed lines represent the p2T distributions for the photon-photon process from two models [33,34] within the STAR acceptance in 60%–80% Au + Au collisions. The dashed lines illustrate the corresponding p2T distributions for e+e− pairs from the model [33] traversing 1 fm in a constant magnetic field of 10^14 T perpendicular to the beam line.

The p2T distributions of excess yields within the STAR acceptance in the mass regions of 0.76–1.2 GeV/c2 in 60%–80% Au + Au and U + U collisions. The systematic uncertainties are shown as gray boxes. The solid and dotted lines are exponential fits to the data in Au + Au and U + U collisions, respectively. The dot-dashed and dot-dot- dashed lines represent the p2T distributions for the photon-photon process from two models [33,34] within the STAR acceptance in 60%–80% Au + Au collisions. The dashed lines illustrate the corresponding p2T distributions for e+e− pairs from the model [33] traversing 1 fm in a constant magnetic field of 10^14 T perpendicular to the beam line.

The p2T distributions of excess yields within the STAR acceptance in the mass regions of 0.76–1.2 GeV/c2 in 60%–80% Au + Au and U + U collisions. The systematic uncertainties are shown as gray boxes. The solid and dotted lines are exponential fits to the data in Au + Au and U + U collisions, respectively. The dot-dashed and dot-dot- dashed lines represent the p2T distributions for the photon-photon process from two models [33,34] within the STAR acceptance in 60%–80% Au + Au collisions. The dashed lines illustrate the corresponding p2T distributions for e+e− pairs from the model [33] traversing 1 fm in a constant magnetic field of 10^14 T perpendicular to the beam line.

The p2T distributions of excess yields within the STAR acceptance in the mass regions of 1.2–2.6 GeV/c2 in 60%–80% Au + Au and U + U collisions. The systematic uncertainties are shown as gray boxes. The solid and dotted lines are exponential fits to the data in Au + Au and U + U collisions, respectively. The dot-dashed and dot-dot- dashed lines represent the p2T distributions for the photon-photon process from two models [33,34] within the STAR acceptance in 60%–80% Au + Au collisions. The dashed lines illustrate the corresponding p2T distributions for e+e− pairs from the model [33] traversing 1 fm in a constant magnetic field of 10^14 T perpendicular to the beam line.

The p2T distributions of excess yields within the STAR acceptance in the mass regions of 1.2–2.6 GeV/c2 in 60%–80% Au + Au and U + U collisions. The systematic uncertainties are shown as gray boxes. The solid and dotted lines are exponential fits to the data in Au + Au and U + U collisions, respectively. The dot-dashed and dot-dot- dashed lines represent the p2T distributions for the photon-photon process from two models [33,34] within the STAR acceptance in 60%–80% Au + Au collisions. The dashed lines illustrate the corresponding p2T distributions for e+e− pairs from the model [33] traversing 1 fm in a constant magnetic field of 10^14 T perpendicular to the beam line.

The corresponding sqrt(p2T) of excess yields. The vertical bars on data points are the combined statistical and systematic uncertainties.

The corresponding sqrt(p2T) of excess yields. The vertical bars on data points are the combined statistical and systematic uncertainties.