Normalized $p_{\mathrm{T}}$-spectrum ratio as a function as a function of centrality in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $E_{\mathrm{T}} \in$ $0.5 \leq |\eta|\leq 0.8$.

Normalized $p_{\mathrm{T}}$-spectrum ratio as a function as a function of centrality in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{tracklets}} \in$ $0.5 \leq |\eta|\leq 0.8$.

Normalized $p_{\mathrm{T}}$-spectrum ratio as a function as a function of centrality in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}} \in$ $-3.7<\eta<-1.7$ and $2.8 < \eta <5.1$.

Event fraction distribution as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}} \in -3.7<\eta<-1.7$ and $2.8 < \eta <5.1$.

Event fraction distribution as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}} \in |\eta|\leq 0.8$.

Event fraction distribution as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $E_{\mathrm{T}} \in |\eta|\leq 0.8$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $|\eta| \leq 0.8$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $0.5 \leq |\eta| \leq 0.8$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $E_{\mathrm{T}}$ $\in$ $|\eta| \leq 0.8$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $E_{\mathrm{T}}$ $\in$ $0.5 \leq |\eta| \leq 0.8$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{tracklets}}$ $\in$ $|\eta| \leq 0.8$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{tracklets}}$ $\in$ $0.5 \leq |\eta| \leq 0.8$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{tracklets}}$ $\in$ $0.3 \leq |\eta| \leq 0.6$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{tracklets}}$ $\in$ $0.7 \leq |\eta| \leq 1$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $-3.7<\eta<-1.7$ and $2.8 < \eta <5.1$.

Extracted $c_{\mathrm{s}}^{2}$ as a function of the minimum pseudorapidity gap between the centrality estimation and transverse-momentum spectra measurement in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$.

$k_{2}/k_{2}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $|\eta| \leq 0.8$.

$k_{2}/k_{2}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $E_{\mathrm{T}}$ $\in$ $|\eta| \leq 0.8$.

$k_{2}/k_{2}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $-3.7<\eta<-1.7$ and $2.8 < \eta <5.1$.

$k_{3}/k_{3}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $|\eta| \leq 0.8$.

$k_{3}/k_{3}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $E_{\mathrm{T}}$ $\in$ $|\eta| \leq 0.8$.

$k_{3}/k_{3}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $-3.7<\eta<-1.7$ and $2.8 < \eta <5.1$.

$\gamma_{\langle [p_{\mathrm{T}}] \rangle}/\gamma_{\langle [p_{\mathrm{T}}] \rangle}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $|\eta| \leq 0.8$.

$\gamma_{\langle [p_{\mathrm{T}}] \rangle}/\gamma_{\langle [p_{\mathrm{T}}] \rangle}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $E_{\mathrm{T}}$ $\in$ $|\eta| \leq 0.8$.

$\gamma_{\langle [p_{\mathrm{T}}] \rangle}/\gamma_{\langle [p_{\mathrm{T}}] \rangle}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $-3.7<\eta<-1.7$ and $2.8 < \eta <5.1$.

$\Gamma_{\langle [p_{\mathrm{T}}] \rangle}/\Gamma_{\langle [p_{\mathrm{T}}] \rangle}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $|\eta| \leq 0.8$.

$\Gamma_{\langle [p_{\mathrm{T}}] \rangle}/\Gamma_{\langle [p_{\mathrm{T}}] angle}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $E_{\mathrm{T}}$ $\in$ $|\eta| \leq 0.8$.

$\Gamma_{\langle [p_{\mathrm{T}}] \rangle}/\Gamma_{\langle [p_{\mathrm{T}}] \rangle}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $-3.7<\eta<-1.7$ and $2.8 < \eta <5.1$.

$\kappa_{\langle [p_{\mathrm{T}}] \rangle}/\kappa_{\langle [p_{\mathrm{T}}] \rangle}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $|\eta| \leq 0.8$.

$\kappa_{\langle [p_{\mathrm{T}}] \rangle}/\kappa_{\langle [p_{\mathrm{T}}] \rangle}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $E_{\mathrm{T}}$ $\in$ $|\eta| \leq 0.8$.

$\kappa_{\langle [p_{\mathrm{T}}] \rangle}/\kappa_{\langle [p_{\mathrm{T}}] \rangle}^{0-5\%}$ as a function of $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $-3.7<\eta<-1.7$ and $2.8 < \eta <5.1$.

Correlation between $\langle p_{\mathrm{T}} \rangle / \langle p_{\mathrm{T}} \rangle ^{\mathrm{0-5\%}}$ and $\langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle / \langle \mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta \rangle ^{\mathrm{0-5\%}}$ in $\mathrm{Pb}-\mathrm{Pb}$ collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02~\mathrm{TeV}$. Data points are shown for centrality estimator based on $N_{\mathrm{ch}}$ $\in$ $0.5 \leq |\eta| \leq 0.8$ with $0.45 \leq p_{\mathrm{T}} \leq 50~\mathrm{GeV}/c.$.

Cutflows and signal efficiencies for the Stealth SYY SUSY model in the $1\ell$ channel corresponding to two values of $m_{\tilde{t}}$.

Cutflows and signal efficiencies for the RPV SUSY model in the $2\ell$ channel corresponding to two values of $m_{\tilde{t}}$.

Cutflows and signal efficiencies for the Stealth SYY SUSY model in the $2\ell$ channel corresponding to two values of $m_{\tilde{t}}$.

Signal efficiencies for the RPV SUSY model in the $0\ell$ channel corresponding to the four ABCD regions

Signal efficiencies for the SYY SUSY model in the $0\ell$ channel corresponding to the four ABCD regions

Signal efficiencies for the RPV SUSY model in the $1\ell$ channel corresponding to the four ABCD regions

Signal efficiencies for the SYY SUSY model in the $1\ell$ channel corresponding to the four ABCD regions

Signal efficiencies for the RPV SUSY model in the $2\ell$ channel corresponding to the four ABCD regions

Signal efficiencies for the SYY SUSY model in the $2\ell$ channel corresponding to the four ABCD regions

Expected and observed 95% CL upper limit on the top squark pair production cross section as a function of the top squark mass for the RPV SUSY model in the $0\ell$ channel

Expected and observed 95% CL upper limit on the top squark pair production cross section as a function of the top squark mass for the RPV SUSY model in the $1\ell$ channel

Expected and observed 95% CL upper limit on the top squark pair production cross section as a function of the top squark mass for the RPV SUSY model in the $2\ell$ channel

Expected and observed 95% CL upper limit on the top squark pair production cross section as a function of the top squark mass for the RPV SUSY model in the Combo channel

The $N_{jets}$ distributions from the background-only fits to data in the $0\ell$. The signal distribution overlaid corresponds to the RPV model with $m_{\tilde{t}} = 400$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $1\ell$. The signal distribution overlaid corresponds to the RPV model with $m_{\tilde{t}} = 400$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $2\ell$. The signal distribution overlaid corresponds to the RPV model with $m_{\tilde{t}} = 400$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $0\ell$. The signal distribution overlaid corresponds to the RPV model with $m_{\tilde{t}} = 800$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $1\ell$. The signal distribution overlaid corresponds to the RPV model with $m_{\tilde{t}} = 800$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $2\ell$. The signal distribution overlaid corresponds to the RPV model with $m_{\tilde{t}} = 800$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

Expected and observed 95% CL upper limit on the top squark pair production cross section as a function of the top squark mass for the SYY SUSY model in the $0\ell$ channel

Expected and observed 95% CL upper limit on the top squark pair production cross section as a function of the top squark mass for the SYY SUSY model in the $1\ell$ channel

Expected and observed 95% CL upper limit on the top squark pair production cross section as a function of the top squark mass for the SYY SUSY model in the $2\ell$ channel

Expected and observed 95% CL upper limit on the top squark pair production cross section as a function of the top squark mass for the SYY SUSY model in the Combo channel

The $N_{jets}$ distributions from the background-only fits to data in the $0\ell$. The signal distribution overlaid corresponds to the Stealth SYY model with $m_{\tilde{t}} = 400$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $1\ell$. The signal distribution overlaid corresponds to the Stealth SYY model with $m_{\tilde{t}} = 400$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $2\ell$. The signal distribution overlaid corresponds to the Stealth SYY model with $m_{\tilde{t}} = 400$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $0\ell$. The signal distribution overlaid corresponds to the Stealth SYY model with $m_{\tilde{t}} = 800$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $1\ell$. The signal distribution overlaid corresponds to the Stealth SYY model with $m_{\tilde{t}} = 800$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The $N_{jets}$ distributions from the background-only fits to data in the $2\ell$. The signal distribution overlaid corresponds to the Stealth SYY model with $m_{\tilde{t}} = 800$ GeV with the signal strength parameter $r$ set to one. The four pannels in each row show the $N_{jets}$ distributions for the A, B, C, and D regions (left to right).

The distributions of the $\mathrm{D}^{*}-\mathrm{D}^0$ mass difference $\Delta m$ for the normalization channel in data.

The distributions of the dimuon invariant mass $m_{\mu\mu}$ for the signal channel in data with the requirement $0.145<\Delta m<0.146$ GeV.

The distributions of the $\mathrm{D}^{*}-\mathrm{D}^0$ mass difference $\Delta m$ for the signal channel in data with the requirement $1.84<m_{\mu\mu}<1.89$ GeV.

The post-fit event yields for the signal, the combinatorial background, the $\mathrm{D}^0\to\pi^+\pi^-$ background, and the $\mathrm{D}^0\to\pi^-\mu^+\nu$ background. The observed numbers of events are given in the Data column. The subrange is in one dimension with a full range in the other dimension.

Distributions of the total transverse mass $M_{T}^{tot}$ in the SRs, comparing observed data with the SM prediction in the $\mu\tau_{h}$ final states in 2018 (center right) after the simultaneous maximum likelihood fit. Representative signal distributions are shown for the 2HDM+a (dashed red curve) and baryonic Z' (dashed black curve) models. The data points are shown with their statistical uncertainties, and the last bin includes overflow. The ``Other MC'' background contribution includes events from ggh, VBF, Wh, Zh, and electroweak vector boson production. The uncertainty band accounts for all systematic and statistical sources of uncertainty, after the fit to the data.

Distributions of the total transverse mass $M_{T}^{tot}$ in the SRs, comparing observed data with the SM prediction in the $\tau_{h}\tau_{h}$ final states in 2017 (lower left) after the simultaneous maximum likelihood fit. Representative signal distributions are shown for the 2HDM+a (dashed red curve) and baryonic Z' (dashed black curve) models. The data points are shown with their statistical uncertainties, and the last bin includes overflow. The ``Other MC'' background contribution includes events from ggh, VBF, Wh, Zh, and electroweak vector boson production. The uncertainty band accounts for all systematic and statistical sources of uncertainty, after the fit to the data.

Distributions of the total transverse mass $M_{T}^{tot}$ in the SRs, comparing observed data with the SM prediction in the $\tau_{h}\tau_{h}$ final states in 2018 (lower right) after the simultaneous maximum likelihood fit. Representative signal distributions are shown for the 2HDM+a (dashed red curve) and baryonic Z' (dashed black curve) models. The data points are shown with their statistical uncertainties, and the last bin includes overflow. The ``Other MC'' background contribution includes events from ggh, VBF, Wh, Zh, and electroweak vector boson production. The uncertainty band accounts for all systematic and statistical sources of uncertainty, after the fit to the data.

The $95\%$ CL upper limits on the $\mu$ for 2HDM+a model. Observed and expected $95\%$ CL upper limits on the $\mu$ as functions of $m_{a}$ using $m_A$ = 600 GeV, $\tan\beta$ = 1, and $m_\chi$ = 10 GeV. The interpolation between the points is linear.

The $95\%$ CL upper limits on the $\mu$ for 2HDM+a model. Observed and expected $95\%$ CL upper limits on the $\mu$ as functions of $m_{A}$ using $m_{a}$ = 150 GeV, $\sin\theta$ = 0.35, and $\tan\beta$ = 1. The interpolation between the points is linear.

The $95\%$ CL upper limits on the $\mu$ for 2HDM+a model. Observed and expected $95\%$ CL upper limits on the $\mu$ as functions of $\sin\theta$ using $m_{a}$ = 200 GeV, $m_A$ = 600 GeV, $\tan\beta$ = 1, and $m_\chi$ = 10 GeV. The interpolation between the points is linear.

The $95\%$ CL upper limits on the $\mu$ for 2HDM+a model. Observed and expected $95\%$ CL upper limits on the $\mu$ as functions of $\tan\beta$ using $m_a$ = 150 GeV, $m_A$ = 600 GeV, $\sin\theta$ = 0.35, and $m_\chi$ = 10 GeV. The interpolation between the points is linear.

The $95\%$ CL upper limits on the $\mu$ in the $m_a-m_A$ plane for 2HDM+a model. The regions inside the red and black curves correspond to the observed and expected exclusions at $95\%$ CL, respectively.

The $95\%$ CL upper limits on the $\mu$ for baryonic Z' model. Observed and expected $95\%$ CL upper limits on the $\mu$ as a function of $m_{Z'}$, using $m_{\chi}$ = 1 GeV. The interpolation between the points is linear.

The $95\%$ CL upper limits on the $\mu$ for baryonic Z' model. $95\%$ CL upper limits on the $\mu$ in $m_{Z'}-m_{\chi}$ plane. The regions inside the red and black curve correspond to the observed and expected exclusions at $95\%$ CL, respectively.

The unfolded charged particle thrust distribution of inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to 1.

The unfolded charged particle broadening distribution of inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to 1.

The unfolded charged particle isotropy distribution of inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to 1.

The unfolded charged particle transverse thrust distribution of inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to 1.

The unfolded charged particle transverse spherocity distribution of inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to 1.

The correlations of the uncertainties in the normalised charged particle multiplicity distribution of inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4.

The correlations of the uncertainties in the normalised charged particle invariant mass distribution of inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4.

The correlations of the uncertainties in the normalised charged particle sphericity distribution of inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4.

The correlations of the uncertainties in the normalised charged particle thrust distribution of inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4.

The correlations of the uncertainties in the normalised charged particle broadening distribution of inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4.

The correlations of the uncertainties in the normalised charged particle isotropy distribution of inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4.

The correlations of the uncertainties in the normalised charged particle transverse thrust distribution of inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4.

The correlations of the uncertainties in the normalised charged particle transverse spherocity distribution of inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4.

The unfolded charged particle invariant mass distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle invariant mass distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle invariant mass distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle invariant mass distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle invariant mass distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle sphericity distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle sphericity distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle sphericity distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle sphericity distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle sphericity distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle thrust distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle thrust distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle thrust distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle thrust distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle thrust distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle broadening distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle broadening distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle broadening distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle broadening distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle broadening distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle isotropy distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle isotropy distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle isotropy distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle isotropy distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle isotropy distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle transverse thrust distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle transverse thrust distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle transverse thrust distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle transverse thrust distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle transverse thrust distribution in a slice of multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle transverse spherocity distribution in a slice from multiplicity of inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle transverse spherocity distribution in a slice from multiplicity of inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle transverse spherocity distribution in a slice from multiplicity of inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle transverse spherocity distribution in a slice from multiplicity of inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded charged particle transverse spherocity distribution in a slice from multiplicity of inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to the fraction of events in this panel among all the panels.

The unfolded joint normalised distribution of charged particle invariant mass and multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to 1.

The unfolded joint normalised distribution of charged particle sphericity and multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to 1.

The unfolded joint normalised distribution of charged particle thrust and multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to 1.

The unfolded joint normalised distribution of charged particle broadening and multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to 1.

The unfolded joint normalised distribution of charged particle isotropy and multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to 1.

The unfolded joint normalised distribution of charged particle transverse thrust and multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to 1.

The unfolded joint normalised distribution of charged particle transverse spherocity and multiplicity of inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4. The total area of the histogram is normalised to 1.

The correlations of the uncertainties in the joint normalised distribution of charged particle invariant mass and multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4.

The correlations of the uncertainties in the joint normalised distribution of charged particle sphericity and multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4.

The correlations of the uncertainties in the joint normalised distribution of charged particle thrust and multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4.

The correlations of the uncertainties in the joint normalised distribution of charged particle broadening and multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4.

The correlations of the uncertainties in the joint normalised distribution of charged particle isotropy and multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4.

The correlations of the uncertainties in the joint normalised distribution of charged particle transverse thrust and multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4.

The correlations of the uncertainties in the joint normalised distribution of charged particle transverse spherocity and multiplicity from inelastic proton-proton collisions with at least three charged particles with transverse momentum higher than 0.5 GeV and pseudorapidity between -2.4 and 2.4.

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{lead.\,lep.}}$.

Fiducial differential cross-sections as a function of $p_{\mathrm{T}}^{\mathrm{sub-lead.\,lep.}}$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.\,lep.}}$.

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T}}^{\mathrm{sub-lead.\,lep.}}$.

Fiducial differential cross-sections as a function of $p_{\mathrm{T},e\mu}$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T},e\mu}$.

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T},e\mu}$.

Fiducial differential cross-sections as a function of $|y_{e\mu}|$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $|y_{e\mu}|$.

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $|y_{e\mu}|$.

Fiducial differential cross-sections as a function of $m_{e\mu}$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $m_{e\mu}$.

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $m_{e\mu}$.

Fiducial differential cross-sections as a function of $|\Delta\phi_{e\mu}|$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $|\Delta\phi_{e\mu}|$.

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $|\Delta\phi_{e\mu}|$.

Fiducial differential cross-sections as a function of $\cos(\theta^{*})$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $\cos(\theta^{*})$.

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $\cos(\theta^{*})$.

Fiducial differential cross-sections as a function of $E_{\mathrm{T}}^{\mathrm{miss}}$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $E_{\mathrm{T}}^{\mathrm{miss}}$.

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $E_{\mathrm{T}}^{\mathrm{miss}}$.

Fiducial differential cross-sections as a function of $H_{\mathrm{T}}^{\mathrm{lep}+\mathrm{MET}}$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $H_{\mathrm{T}}^{\mathrm{lep}+\mathrm{MET}}$.

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $H_{\mathrm{T}}^{\mathrm{lep}+\mathrm{MET}}$.

Fiducial differential cross-sections as a function of $m_{\mathrm{T},e\mu}$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $m_{\mathrm{T},e\mu}$.

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $m_{\mathrm{T},e\mu}$.

Fiducial differential cross-sections as a function of the number of jets. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable Number of jets ($p_{\mathrm{T}} > $ 30 GeV).

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable Number of jets ($p_{\mathrm{T}} > $ 30 GeV).

Fiducial differential cross-sections as a function of $S_{\mathrm{T}}$. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The results are compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1, with the NNLO + PS predictions from Powheg MiNNLO + Pythia8 and Geneva + Pythia8, as well as Sherpa2.2.12 NLO + PS predictions. The last three predictions are combined with Sherpa 2.2.2 for the $gg$ initial state and Sherpa 2.2.12 for electroweak $WWjj$ production. These contributions are modelled at LO but a NLO QCD $k$-factor of 1.7 is applied for gluon induced production. Theoretical predictions are indicated as markers with vertical lines denoting PDF, scale and parton shower uncertainties. Markers are staggered for better visibility.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $S_{\mathrm{T}}$.

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $S_{\mathrm{T}}$.

Fiducial differential cross-sections for $p_{\mathrm{T},e\mu}$ in the region with 85 $<m_{e\mu}<$ 150, for the nNNLO MATRIX prediction and various four flavour PDF sets. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The measurement is compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1 using four different four-flavour NNLO PDF sets: NNPDF3.0, NNPDF3.1, CT18, and MSHT20.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T},e\mu}$ (85 $<m_{e\mu}<$ 150).

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T},e\mu}$ (85 $<m_{e\mu}<$ 150).

Fiducial differential cross-sections for $|y_{e\mu}|$ in the region with 85 $<m_{e\mu}<$ 150, for the nNNLO MATRIX prediction and various four flavour PDF sets. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The measurement is compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1 using four different four-flavour NNLO PDF sets: NNPDF3.0, NNPDF3.1, CT18, and MSHT20.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $|y_{e\mu}|$ (85 $<m_{e\mu}<$ 150).

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $|y_{e\mu}|$ (85 $<m_{e\mu}<$ 150).

Fiducial differential cross-sections for $p_{\mathrm{T},e\mu}$ in the region with 150 $<m_{e\mu}<$ 250, for the nNNLO MATRIX prediction and various four flavour PDF sets. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The measurement is compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1 using four different four-flavour NNLO PDF sets: NNPDF3.0, NNPDF3.1, CT18, and MSHT20.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T},e\mu}$ (150 $<m_{e\mu}<$ 250).

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T},e\mu}$ (150 $<m_{e\mu}<$ 250).

Fiducial differential cross-sections for $|y_{e\mu}|$ in the region with 150 $<m_{e\mu}<$ 250, for the nNNLO MATRIX prediction and various four flavour PDF sets. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The measurement is compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1 using four different four-flavour NNLO PDF sets: NNPDF3.0, NNPDF3.1, CT18, and MSHT20.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $|y_{e\mu}|$ (150 $<m_{e\mu}<$ 250).

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $|y_{e\mu}|$ (150 $<m_{e\mu}<$ 250).

Fiducial differential cross-sections for $p_{\mathrm{T},e\mu}$ in the region with $m_{e\mu}>$ 250, for the nNNLO MATRIX prediction and various four flavour PDF sets. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The right-hand-side axis indicates the integrated cross-section of the rightmost bin. The measurement is compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1 using four different four-flavour NNLO PDF sets: NNPDF3.0, NNPDF3.1, CT18, and MSHT20.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T},e\mu}$ ($m_{e\mu}>$ 250).

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $p_{\mathrm{T},e\mu}$ ($m_{e\mu}>$ 250).

Fiducial differential cross-sections for $|y_{e\mu}|$ in the region with $m_{e\mu}>$ 250, for the nNNLO MATRIX prediction and various four flavour PDF sets. The measured cross-section values are shown as points with error bars giving the statistical uncertainty and solid bands indicating the size of the total uncertainty. The measurement is compared to fixed-order nNNLO QCD + NLO EW predictions of Matrix 2.1 using four different four-flavour NNLO PDF sets: NNPDF3.0, NNPDF3.1, CT18, and MSHT20.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $|y_{e\mu}|$ ($m_{e\mu}>$ 250).

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $|y_{e\mu}|$ ($m_{e\mu}>$ 250).

Measured value of the charge asymmetry, $A_{C}$.

Charge asymmetry $A_C$ measured as a function of the absolute difference of the absolute lepton rapidities $||\eta_{l+}|-|\eta_{l-}||$. The measured values are shown as points, whose error bars represent the statistical uncertainty, while the solid bands indicate the size of the total uncertainty. In the bottom panel, the central value of the charge asymmetry divided by the total uncertainty of the $A_C$ measurement is shown. The measurement is compared with the NNLO+PS predictions from Powheg MiNNLO+Pythia8 with vertical lines denoting PDF uncertainties.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $A_{C}$ vs. $||\eta_{l+}|-|\eta_{l-}||$.

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $A_{C}$ vs. $||\eta_{l+}|-|\eta_{l-}||$.

Charge asymmetry $A_C$ measured as a function of the invariant mass of the dilepton system $m_{e\mu}$. The measured values are shown as points, whose error bars represent the statistical uncertainty, while the solid bands indicate the size of the total uncertainty. In the bottom panel, the central value of the charge asymmetry divided by the total uncertainty of the $A_C$ measurement is shown. The measurement is compared with the NNLO+PS predictions from Powheg MiNNLO+Pythia8 with vertical lines denoting PDF uncertainties.

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $A_{C}$ vs. $m_{e\mu}$.

Correlation matrix of the total uncertainties in the measured fiducial cross section for the observable $A_{C}$ vs. $m_{e\mu}$.

Measured value of the asymmetry in the CP-odd observable $O_{z}$, $A_{CP}$.

Corrected Yield R=0.5 pi0+jet 15-20 pp at sqrt{s_{NN}}=200 GeV

Corrected Yield R=0.2 gamma+jet 10-15 pp at sqrt{s_{NN}}=200 GeV

Corrected Yield R=0.5 gamma+jet 10-15 pp at sqrt{s_{NN}}=200 GeV

Corrected Yield R=0.5 pi0+jet 10-15 AuAu 0-15% at sqrt{s_{NN}}=200 GeV

Corrected Yield R=0.2 pi0+jet 10-15 AuAu 0-15% at sqrt{s_{NN}}=200 GeV

Corrected Yield R=0.5 pi0+jet 15-20 AuAu 0-15% at sqrt{s_{NN}}=200 GeV

Corrected Yield R=0.2 pi0+jet 15-20 AuAu 0-15% at sqrt{s_{NN}}=200 GeV

Corrected Yield R=0.2 gamma+jet 10-15 AuAu 0-15% at sqrt{s_{NN}}=200 GeV

Corrected Yield R=0.5 gamma+jet 10-15 AuAu 0-15% at sqrt{s_{NN}}=200 GeV

IAA(Δφ) R=0.5 gamma+jet 10-15 GeV/c

IAA(Δφ) R=0.2 gamma+jet 10-15 GeV/c

IAA(Δφ) R=0.5 pi0+jet 10-15 GeV/c

IAA(Δφ) R=0.2 pi0+jet 10-15 GeV/c

IAA(Δφ) R=0.5 pi0+jet 15-20 GeV/c

IAA(Δφ) R=0.2 pi0+jet 15-20 GeV/c

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the 1 b jet, 0 forward jet, $t\leq 0$ signal region in the single lepton channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the 1 b jet, 0 forward jet, $t>0$ signal region in the single lepton channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the 1 b jet, $\geq 1$ forward jet, $t\leq 0$ signal region in the single lepton channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the 1 b jet, $\geq 1$ forward jet, $t>0$ signal region in the single lepton channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the $\geq 2$ b jet, $t\leq 0$ signal region in the single lepton channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the $\geq 2$ b jet, $t>0$ signal region in the single lepton channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit NN output score distribution of the $\geq 2$ b jet, different flavor signal region in the dileptonic channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit NN output score distribution of the $\geq 2$ b jet, same flavor signal region in the dileptonic channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit NN output score distribution of the 1 b jet, different flavor signal region in the dileptonic channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The post-fit NN output score distribution of the 1 b jet, same flavor signal region in the dileptonic channel. A representative signal model distribution is shown for the scalar mediator interaction with $(m_{\chi},m_{\phi})=(1,100)$GeV and couplings set to unity. The grey dashed area in the upper panel represents the total uncertainty in all of the backgrounds and the chosen signal model, while in the lower panel it represents only the total uncertainty in the backgrounds.

The model-independent 95\% CL limits on production cross section for new physics processes for the scalar interaction. The expected limit is shown by the black dashed line with the 68 and 95\% CL uncertainty bands in green and yellow, respectively, while the observed limit is shown by the solid black line. Theoretical LO cross section values for the DM model and their associated uncertainties are also presented (grey line).

The model-independent 95\% CL limits on production cross section for new physics processes for the pseudoscalar interaction. The expected limit is shown by the black dashed line with the 68 and 95\% CL uncertainty bands in green and yellow, respectively, while the observed limit is shown by the solid black line. Theoretical LO cross section values for the DM model and their associated uncertainties are also presented (grey line).

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the QCD control region in the all hadronic channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the tt(1l) control region in the all hadronic channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the W+jets control region in the all hadronic channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the Z(2l) control region in the all hadronic channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the W+jets control region in the single lepton channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit $p_{\mathrm{T}}^{\text{miss}}$ distribution of the tt(2l) control region in the single lepton channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit $(N_{\text{jet}},N_{\text{bjet}})$ distribution of the ttZ control region in the dileptonic channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

The post-fit NN distribution of the Drell-Yan control region in the dileptonic channel. The grey dashed area in the upper and lower panels represents the total uncertainty in all of the backgrounds.

Predicted signal yields for the 2016, 2017, and 2018 data-taking periods, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{\phi}) = (1,50)$ GeV scalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{\phi}) = (1,500)$ GeV scalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{a}) = (1,50)$ GeV pseudoscalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{a}) = (1,500)$ GeV pseudoscalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{\phi}) = (1,50)$ GeV scalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{\phi}) = (1,500)$ GeV scalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{a}) = (1,50)$ GeV pseudoscalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{a}) = (1,500)$ GeV pseudoscalar mediator signal hypothesis for the all hadronic channel. The requirements listed up to and including the $M_{T}^{b}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last three requirements show the separate yields after the cumulative requirements up to the $M_{T}^{b}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{\phi}) = (1,50)$ GeV scalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{\phi}) = (1,500)$ GeV scalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{a}) = (1,50)$ GeV pseudoscalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{a}) = (1,500)$ GeV pseudoscalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{\phi}) = (1,50)$ GeV scalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{\phi}) = (1,500)$ GeV scalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{a}) = (1,50)$ GeV pseudoscalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{a}) = (1,500)$ GeV pseudoscalar mediator signal hypothesis for the single lepton channel. The requirements listed up to and including the $M_{T2}^{W}$ cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last six requirements show the separate yields after the cumulative requirements up to the $M_{T2}^{W}$ cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{\phi}) = (1,50)$ GeV scalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{\phi}) = (1,500)$ GeV scalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{a}) = (1,50)$ GeV pseudoscalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the t+DM $(m_{\chi},m_{a}) = (1,500)$ GeV pseudoscalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{\phi}) = (1,50)$ GeV scalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{\phi}) = (1,500)$ GeV scalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{a}) = (1,50)$ GeV pseudoscalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Predicted signal yields for the 2016, 2017, and 2018 data-taking period, corresponding to $36.3\,\text{fb}^{-1}$, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, $41.5\,\text{fb}^{-1}$, and $59.7\,\text{fb}^{-1}$, respectively, after the application of each search requirement for the tt+DM $(m_{\chi},m_{a}) = (1,500)$ GeV pseudoscalar mediator signal hypothesis for the dileptonic channel. The requirements listed up to and including the Z window cut are cumulative such that only the events and objects passing a given requirement are considered in subsequent requirements. The last four requirements show the separate yields after the cumulative requirements up to the Z window cut.

Cross sections calculated to leading order (LO) accuracy for the tt+DM scalar mediator signal process considering all possible decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the tt+DM pseudoscalar mediator signal process considering all possible decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the tt+DM scalar mediator signal process considering dileptonic decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the tt+DM pseudoscalar mediator signal process considering dileptonic decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the t+DM t-channel scalar mediator signal process considering all possible decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the t+DM tW-channel scalar mediator signal process considering all possible decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the t+DM t-channel pseudoscalar mediator signal process considering all possible decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the t+DM tW-channel pseudoscalar mediator signal process considering all possible decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the t+DM tW-channel scalar mediator signal process considering dileptonic decays from the visible particles in the final state.

Cross sections calculated to leading order (LO) accuracy for the t+DM tW-channel pseudoscalar mediator signal process considering dileptonic decays from the visible particles in the final state.

Global polarization of $\Lambda$ and $\bar\Lambda$ and their difference in 20-50\% centrality in Ru+Ru collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Global polarization of $\Lambda$ and $\bar\Lambda$ and their difference in 20-50\% centrality in Zr+Zr collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Global polarization of $\Lambda$ and $\bar\Lambda$ and their difference in 20-50\% centrality in Ru+Ru&Zr+Zr collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Global polarization of $\Lambda$+$\bar\Lambda$ as a function of centrality in the combined Ru+Ru and Zr+Zr collisions at $\sqrt{s_{NN}}$ = 200 GeV (left panel). The data in Au+Au collisions is from Phys. Rev. C 98 (2018) 014910.

Global polarization of $\Lambda$+$\bar\Lambda$ in 20-50\% centrality centrality in the combined Ru+Ru and Zr+Zr collisions at $\sqrt{s_{NN}}$ = 200 GeV (right panel). The data in Au+Au collisions is from Phys. Rev. C 98 (2018) 014910.

Global polarization of $\Lambda$+$\bar\Lambda$ as a function of the average number of participants $\langle N_{\rm part} \rangle$ in the combined Ru+Ru and Zr+Zr collisions at $\sqrt{s_{NN}}$ = 200 GeV. The data in Au+Au collisions is from Phys. Rev. C 98 (2018) 014910.

The transverse momentum dependence of $\Lambda$+$\bar\Lambda$ $P_{\rm H}$ for 20\%-60\% centrality in the combined Ru+Ru and Zr+Zr collisions at $\sqrt{s_{NN}}$ = 200 GeV, together with the results for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The data in Au+Au collisions is from Phys. Rev. C 98 (2018) 014910.

The pseudorapidity dependence of $\Lambda$+$\bar\Lambda$ $P_{\rm H}$ for 20\%-60\% centrality in the combined Ru+Ru and Zr+Zr collisions at $\sqrt{s_{NN}}$ = 200 GeV, together with the results for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The data in Au+Au collisions is from Phys. Rev. C 98 (2018) 014910.

The polarization coefficient $P_{y,c0}$ of $\Lambda\!+\!\bar\Lambda$ from method-1 and method-2 as a function of centrality (left panel) in isobar Ru+Ru and Zr+Zr collisions at $\sqrt{s_{NN}}$ = 200 GeV. In this figure, $P_{H}$ is same data as in Figure 4.

The polarization coefficient $P_{y,c0}$ of $\Lambda\!+\!\bar\Lambda$ from method-1 and method-2 for 20-50\% centrality (right panel) in isobar Ru+Ru and Zr+Zr collisions at $\sqrt{s_{NN}}$ = 200 GeV. In this figure, $P_{H}$ is same data as in Figure 4.

The polarization coefficient $P_{y,c2}$ of $\Lambda\!+\!\bar\Lambda$ from method-1 and method-2 as a function of centrality (left panel) in isobar Ru+Ru and Zr+Zr collisions at $\sqrt{s_{NN}}$ = 200 GeV.

The polarization coefficient $P_{y,c2}$ of $\Lambda\!+\!\bar\Lambda$ from method-1 and method-2 for 20-50\% centrality (right panel) in isobar Ru+Ru and Zr+Zr collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Global polarization of $\Xi^{-}$+$\bar\Xi^{+}$ from polarization transfer method and direct measurement method as a function of centrality (left panel) in Ru+Ru and Zr+Zr collisions at $\sqrt{s_{NN}}$ = 200 GeV. In this figure, $P_{H}$ of inclusive $\Lambda+\bar\Lambda$ is same data as in Figure 4.

Global polarization of $\Xi^{-}$+$\bar\Xi^{+}$ from polarization transfer method and direct measurement method for 20-50\% centrality (right panel) in Ru+Ru and Zr+Zr collisions at $\sqrt{s_{NN}}$ = 200 GeV . In this figure, $P_{H}$ of inclusive $\Lambda+\bar\Lambda$ is same data as in Figure 4.