Transverse dilepton-$p_{T}^{miss}$ mass pre-fit distribution for selected events in SR2 of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton mass pre-fit distribution for selected events in SR3 of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton $p_T$ pre-fit distribution for selected events in SR3 of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dijet $p_T$ pre-fit distribution for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Lepton-dijet $p_T$ pre-fit distribution for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Maximum $p_T$ of the visible particles divided by the missing $p_T$, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between jets, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and missing $p_T$, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and the dijet system, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta R$ distance between the leptons, for selected events in the non-prompt validation region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the di-lepton system and the missing $p_T$, for selected events in the non-prompt validation region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and the di-jet system, for selected events in the non-prompt validation region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton-dijet system and the missing $p_T$, for selected events in the non-prompt validation region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Leading lepton $p_T$ pre-fit distribution for selected events in the top control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Trailing lepton $p_T$ pre-fit distribution for selected events in the top control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton mass pre-fit distribution for selected events in the Drell-Yan control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Missing $p_T$ pre-fit distribution for selected events in the Drell-Yan control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Missing $p_T$ pre-fit distribution for selected events in the $WW$ control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton $p_T$ pre-fit distribution for selected events in the $WW$ control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Minimum $p_T$ of the visible particles divided by the missing $p_T$, for selected events in the top control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between jets, for selected events in the top control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and missing $p_T$, for selected events in the $W$+jets control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and the dijet system, for selected events in the $W$+jets control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Unrolled $m_{ll}-m_{T}^{l min, p_{T}^{miss}}$ post-fit distributions in the di-leptonic channel in SR1 for the full data set. The histogram bins are spaced uniformly. Each group of five bins (from left to right) corresponds to $m_{T}^{l min, p_{T}^{miss}}$ distribution in a $m_{ll}$ region, placed in ascending order. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$.

Unrolled $m_{ll}-m_{T}^{l min, p_{T}^{miss}}$ post-fit distributions in the di-leptonic channel in SR2 for the full data set. The histogram bins are spaced uniformly. Each group of five bins (from left to right) corresponds to $m_{T}^{l min, p_{T}^{miss}}$ distribution in a $m_{ll}$ region, placed in ascending order. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$.

Unrolled $m_{ll}-m_{T}^{l min, p_{T}^{miss}}$ post-fit distributions in the di-leptonic channel in SR3 for the full data set. The histogram bins are spaced uniformly. Each group of five bins (from left to right) corresponds to $m_{T}^{l min, p_{T}^{miss}}$ distribution in a $m_{ll}$ region, placed in ascending order. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$.

Post-fit BDT distributions in the semi-leptonic channel for the full data set in the top control region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Post-fit BDT distributions in the semi-leptonic channel for the full data set in the $W$+jets control region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Post-fit BDT distributions in the semi-leptonic channel in the 2016 signal region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Post-fit BDT distributions in the semi-leptonic channel in the 2017+2018 signal region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

The number of events in data and simulations after the fit to data in the Top$\gamma$ CR, SS$WW\gamma$ CR, and signal region.

The distributions of the invariant mass of the dilepton-photon system after the signal region selection before the fit to data.The shaded band represents the Monte Carlo statistical unc

Upper limits on the cross section and derived limits in terms of Yukawa coupling at 95% CL for Hγ production initiated by light quarks.

The distributions of the invariant mass of the transverse mass of the WW system after the signal region selection before the fit to data.The shaded band represents the Monte Carlo statistical unc

The distributions of the invariant mass of the dilepton-photon system after the signal region selection before the fit to data.The shaded band represents the Monte Carlo statistical unc

The unrolled two-dimensional mTWW-mll$\gamma$ distributions with category 0 jet after the fit to data. The data are compared to the sum of the signal and expected background. The black points with error bars represent the data and their statistical uncertainties, whereas the hatched bands represent the total uncertainties in the predictions.

The distributions of the invariant mass of the transverse mass of the WW system after the signal region selection before the fit to data.The shaded band represents the Monte Carlo statistical unc

The unrolled two-dimensional mTWW-mll$\gamma$ distributions with category 1more jets after the fit to data. The data are compared to the sum of the signal and expected background. The black points with error bars represent the data and their statistical uncertainties, whereas the hatched bands represent the total uncertainties in the predictions.

The unrolled two-dimensional mTWW-mll$\gamma$ distributions with category 0 jet after the fit to data. The data are compared to the sum of the signal and expected background. The black points with error bars represent the data and their statistical uncertainties, whereas the hatched bands represent the total uncertainties in the predictions.

The unrolled two-dimensional mTWW-mll$\gamma$ distributions with category 1more jets after the fit to data. The data are compared to the sum of the signal and expected background. The black points with error bars represent the data and their statistical uncertainties, whereas the hatched bands represent the total uncertainties in the predictions.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 350 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 400 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 450 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 500 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 550 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 600 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 650 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 700 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 750 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 800 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 850 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 900 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 950 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Expected and observed 95% upper limits on the product of the production cross section and the branching fraction for a spin-0 resonance $X\rightarrow HY\rightarrow \gamma\gamma bb$, as a function of Y mass, for X mass = 1000 GeV. The $\pm1$ and $\pm2$ $\sigma$ uncertainty bands are given in addition to the expected median value. Numerical values provided in this table correspond to Figure 7 of the publication.

Full experimental (statistical and systematic) correlations (in \%) of the normalized partial decay rates over the $\overline{B}^0\to D^{*+} e^- \bar\nu_e$ and $\overline{B}^0\to D^{*+} \mu^- \bar\nu_\mu$ decays.

Measured partial decay rates $\Delta\Gamma$ (in units of $10^{-15}$ GeV) using the matrix inversion unfolding method

Average of normalized decay rates that are determined using the matrix inversion unfolding method

Full experimental (statistical and systematic) correlations (in \%) of the partial decay rates determined with the matrix inversion unfolding method for the $\overline{B}^0\to D^{*+} e^- \bar\nu_e$ and $\overline{B}^0\to D^{*+} \mu^- \bar\nu_\mu$ decays.

Full experimental (statistical and systematic) correlations (in \%) of the normalized partial decay rates determined using the matrix inversion unfolding method

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

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

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

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

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

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

Comparison between the data and SM background predictions for the number of muons in the long track final state.

Comparison between the data and SM background predictions for the m_{DTk; stopping power} in the long track final state.

Comparison between the data and SM background predictions for the number of electrons in the long track final state.

Comparison between the data and SM background predictions for the number of b-tagged jets in the short track final state.

Comparison between the data and SM background predictions for the number of jets in the short track final state.

Comparison between the data and SM background predictions for the hard missing transverse momentum in the short track final state.

Comparison between the data and SM background predictions for the number of muons in the short track final state.

Comparison between the data and SM background predictions for the m_{DTk; stopping power} in the short track final state.

Comparison between the data and SM background predictions for the number of electrons in the short track final state.

Comparison between the data and post-fit SM background predictions for the 49 search regions in the global search region.

95% CL upper limit on the cross section in pb for the T6bt 10 cm model.

Limit in the mass plane of the T6bt 10 cm model

95% CL upper limit on the cross section in pb for the T6bt 200 cm model.

Limit in the mass plane of the T6bt 200 cm model

95% CL upper limit on the cross section in pb for the T6tb 10 cm model.

Limit in the mass plane of the T6tb 10 cm model

95% CL upper limit on the cross section in pb for the T6tb 200 cm model.

Limit in the mass plane of the T6tb 200 cm model

95% CL upper limit on the cross section in pb for the T5btbt 10 cm model.

Limit in the mass plane of the T5btbt 10 cm model

95% CL upper limit on the cross section in pb for the T5btbt 200 cm model.

Limit in the mass plane of the T5btbt 200 cm model

95% CL upper limit on the cross section in pb for the pure Wino model.

Limit in the mass plane of the pure Wino model

95% CL upper limit on the cross section in pb for the pure Higgsino model.

Limit in the mass plane of the pure Higgsino model

Cut flow table showing object-level efficiency at each stage of the track selection. The denominator of the efficiency is the number of chargino-matched tracks that have been reconstructed and satisfy the pT requirements, for which the chargino decays within the pixel detector volume; the numerator is the number of such tracks which also satisfy all requirements up to and including the associated criterion.

Cut flow table showing object-level efficiency at each stage of the track selection. The denominator of the efficiency is the number of chargino-matched tracks that have been reconstructed and satisfy the pT requirements, for which the chargino decays within the strip detector volume; the numerator is the number of such tracks which also satisfy all requirements up to and including the associated criterion.

Percent efficiencies for sum of all SR bins, each SR bin, SR bin groups:

Percent efficiencies for sum of all SR bins, each SR bin, SR bin groups:

Percent efficiencies for sum of all SR bins, each SR bin, SR bin groups:

Percent efficiencies for sum of all SR bins, each SR bin, SR bin groups:

Percent efficiencies for sum of all SR bins, each SR bin, SR bin groups:

Percent efficiencies for sum of all SR bins, each SR bin, SR bin groups: