Post-fit distributon of BDT discriminants for $\rho_{tc}=1.0$ with $m_A$ = 350 GeV interfered with H.($m_A - m_H$ = 50 GeV)

Observed and expected 95% CL upper limit on the signal strength as a function of $m_A$ for the g2HDM signal model using different coupling assumptions: $\rho_{tu}$ = 0.1, 0.4, 1.0.

Observed and expected 95% CL upper limit on the signal strength as a function of $m_A$ for the g2HDM signal model using different coupling assumptions: $\rho_{tc}$ = 0.1, 0.4, 1.0.

Observed and expected 95% CL upper limit on the signal strength as a function of $m_A$ for the g2HDM signal model using different coupling assumptions: $\rho_{tu}$ = 0.1, 0.4, 1.0 with A–H interference assuming $m_A$ − $m_H$ = 50 GeV.

Observed and expected 95% CL upper limit on the signal strength as a function of $m_A$ for the g2HDM signal model using different coupling assumptions: $\rho_{tc}$ = 0.1, 0.4, 1.0 with A–H interference assuming $m_A$ − $m_H$ = 50 GeV.

Observed 95% CL upper limit on the signal strength as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model.

Expected exclusion as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model

Observed exclusion as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model

Observed 95% CL upper limit on the signal strength as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model.

Expected exclusion as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model

Observed exclusion as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model

Observed 95% CL upper limit on the signal strength as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV.

Expected exclusion as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV

Observed exclusion as a function of $m_A$ and $\rho_{tu}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV

Observed 95% CL upper limit on the signal strength as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV.

Expected exclusion as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV

Observed exclusion as a function of $m_A$ and $\rho_{tc}$ for the g2HDM signal model with A–H interference assuming $m_A$ − $m_H$ = 50 GeV

Product of detector efficiency and analysis acceptance for signal samples with various $m_a$ values for the muon channel.

Post-fit distributions of $M_{\ell v jj}$ in the $R2B_A$ subregion for electrons. All process yields and nuisance parameters are set to the values obtained from the background plus signal fit. The signal considered for the fit corresponds to the purely right-handed production of a W' with $m_{W'}$ of 3.6 TeV and a relative width of 1$\%$ of the $m_{W'}$, and is represented by the solid red line. The lower panels show the data minus the expected number of events, normalized to the statistical uncertainty of the data. The orange band represents the systematic uncertainties, also normalized to the statistical uncertainty of the data.

Post-fit distributions of $M_{\ell v jj}$ in the $RW'_A$ subregion for muons. All process yields and nuisance parameters are set to the values obtained from the background plus signal fit. The signal considered for the fit corresponds to the purely right-handed production of a W' with $m_{W'}$ of 3.6 TeV and a relative width of 1$\%$ of the $m_{W'}$, and is represented by the solid red line. The lower panels show the data minus the expected number of events, normalized to the statistical uncertainty of the data. The orange band represents the systematic uncertainties, also normalized to the statistical uncertainty of the data.

Post-fit distributions of $M_{\ell v jj}$ in the $RW'_A$ subregion for electrons. All process yields and nuisance parameters are set to the values obtained from the background plus signal fit. The signal considered for the fit corresponds to the purely right-handed production of a W' with $m_{W'}$ of 3.6 TeV and a relative width of 1$\%$ of the $m_{W'}$, and is represented by the solid red line. The lower panels show the data minus the expected number of events, normalized to the statistical uncertainty of the data. The orange band represents the systematic uncertainties, also normalized to the statistical uncertainty of the data.

Post-fit distributions of $M_{\ell v jj}$ in the $RT_A$ subregion for muons. All process yields and nuisance parameters are set to the values obtained from the background plus signal fit. The signal considered for the fit corresponds to the purely right-handed production of a W' with $m_{W'}$ of 3.6 TeV and a relative width of 1$\%$ of the $m_{W'}$, and is represented by the solid red line. The lower panels show the data minus the expected number of events, normalized to the statistical uncertainty of the data. The orange band represents the systematic uncertainties, also normalized to the statistical uncertainty of the data.

Post-fit distributions of $M_{\ell v jj}$ in the $RT_A$ subregion for electrons. All process yields and nuisance parameters are set to the values obtained from the background plus signal fit. The signal considered for the fit corresponds to the purely right-handed production of a W' with $m_{W'}$ of 3.6 TeV and a relative width of 1$\%$ of the $m_{W'}$, and is represented by the solid red line. The lower panels show the data minus the expected number of events, normalized to the statistical uncertainty of the data. The orange band represents the systematic uncertainties, also normalized to the statistical uncertainty of the data.

Observed and expected 95% CL upper limits on the product of the production cross section for a left-handed production of a tb quark pair in the $s$-channel, mediated by a W or W' bosons and including interference terms, given as functions of $m_{W'}$ for a relative width of 1%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a left-handed production of a tb quark pair in the $s$-channel, mediated by a W or W' bosons and including interference terms, given as functions of $m_{W'}$ for a relative width of 10%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a left-handed production of a tb quark pair in the $s$-channel, mediated by a W or W' bosons and including interference terms, given as functions of $m_{W'}$ for a relative width of 20%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a left-handed production of a tb quark pair in the $s$-channel, mediated by a W or W' bosons and including interference terms, given as functions of $m_{W'}$ for a relative width of 30%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of $m_{W'}$ for a relative width of 1%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of $m_{W'}$ for a relative width of 10%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of $m_{W'}$ for a relative width of 20%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of $m_{W'}$ for a relative width of 30%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid red curves show the theoretical expectation at LO.

Observed 95% CL upper limit on the production cross section for a left-handed W' boson in the tb final state, as functions of $m_{W'}$ and relative width $\Gamma/m_{W'}$. Numbers in red, written diagonally, represent values of the excluded cross sections that are lower than the theoretical ones for the analyzed model.

Observed 95% CL upper limit on the production cross section for a right-handed W' boson in the tb final state, as functions of $m_{W'}$ and relative width $\Gamma/m_{W'}$. Numbers in red, written diagonally, represent values of the excluded cross sections that are lower than the theoretical ones for the analyzed model.

Observed 95% CL upper limit on the production cross section for a generalized left-right coupling of the W' boson to a t and a b quark for a mass of the W' boson of 2 TeV.

Observed 95% CL upper limit on the production cross section for a generalized left-right coupling of the W' boson to a t and a b quark for a mass of the W' boson of 2.8 TeV.

Observed 95% CL upper limit on the production cross section for a generalized left-right coupling of the W' boson to a t and a b quark for a mass of the W' boson of 3.6 TeV.

Observed 95% CL upper limit on the production cross section for a generalized left-right coupling of the W' boson to a t and a b quark for a mass of the W' boson of 4.4 TeV.

Observed 95% CL upper limit on the production cross section for a generalized left-right coupling of the W' boson to a t and a b quark for a mass of the W' boson of 5.2 TeV.

Observed 95% CL upper limit on the production cross section for a generalized left-right coupling of the W' boson to a t and a b quark for a mass of the W' boson of 6 TeV.

Expected 95% CL lower limit on $m_{W'}$ for a generalized left-right coupling of the W' boson to a t and a b quark.

Observed 95% CL lower limit on $m_{W'}$ for a generalized left-right coupling of the W' boson to a t and a b quark.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of the $m_{W'}$ for a relative width of 10%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid curves show the theoretical expectation at LO in the case that the couplings, and thus the partial widths, are varied together with the total width. In this interpretation the branching fraction of the W' --> tb decays is the same for each value of the width of the W' boson.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of the $m_{W'}$ for a relative width of 20%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid curves show the theoretical expectation at LO in the case that the couplings, and thus the partial widths, are varied together with the total width. In this interpretation the branching fraction of the W' --> tb decays is the same for each value of the width of the W' boson.

Observed and expected 95% CL upper limits on the product of the production cross section for a right-handed W' boson and the W' --> tb branching fraction, as functions of the $m_{W'}$ for a relative width of 30%. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The solid curves show the theoretical expectation at LO in the case that the couplings, and thus the partial widths, are varied together with the total width. In this interpretation the branching fraction of the W' --> tb decays is the same for each value of the width of the W' boson.

Observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance. The acceptance $\mathcal{A}$ is defined as $\mathcal{A} = N$(events with $m_{X}^{GEN} > 85\% m_{X}^{input}$) / $N$(events generated in the full phase space defined by the CMS default generator settings). Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Mass exclusion ranges of the benchmark signal scenarios are depicted with hatched areas inside the black contours.

Observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance. The acceptance $\mathcal{A}$ is defined as $\mathcal{A} = N$(events with $m_{X}^{GEN} > 85\% m_{X}^{input}$) / $N$(events generated in the full phase space defined by the CMS default generator settings). Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Mass exclusion ranges of the benchmark signal scenarios are depicted with hatched areas inside the black contours.

Efficiencies of the selection requirements on the benchmark signal processes: $Z_{R} \to ggg$ with nominal width ($\Gamma_{X}/m_{X}\sim 3\%$). A value of -1 means that the corresponding efficiency is not calculated for this year.

Efficiencies of the selection requirements on the benchmark signal processes: $Z_{R} \to ggg$ with narrow width ($\Gamma_{X}/m_{X}\sim 0.01\%$). A value of -1 means that the corresponding efficiency is not calculated for this year.

Efficiencies of the selection requirements on the benchmark signal processes: $G_{KK} \to {\varphi}(gg)g$, where $\varphi$ is the radion. A value of -1 means that the corresponding efficiency is not calculated for this year.

Efficiencies of the selection requirements on the benchmark signal processes: $q^{*} \to V(qq)q$, where $V$ is a beyond-the-SM vector boson. A value of -1 means that the corresponding efficiency is not calculated for this year.

Acceptance of the signal selection requirement $m_{X}^{\text{GEN}}/m_{X}^{\text{input}} > 85\%$ on the benchmark signal process $Z_{R} \to ggg$. The acceptance is defined as $\mathcal{A} = N$(events with $m_{X}^{\text{GEN}}/m_{X}^{\text{input}} > 85\%$)/$N$(events generated in the full phase space defined by the CMS default generator settings).

Acceptance of the signal selection requirement $m_{X}^{\text{GEN}}/m_{X}^{\text{input}} > 85\%$ on the benchmark signal process $G_{KK} \to {\varphi}(gg)g$.

Acceptance of the signal selection requirement $m_{X}^{\text{GEN}}/m_{X}^{\text{input}} > 85\%$ on the benchmark signal process $q^{*} \to V(qq)q$.

Observed local significance for a 3-body decay $ggg$ resonance, shown for resonances with nominal width (blue solid line) and narrow width (red dashed line).

Observed local significance for a cascade decay $ggg$ resonance.

Observed local significance for a cascade decay $qqq$ resonance.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.2$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.3$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.4$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.5$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.6$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.7$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(gg)g) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.8$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance $m_{Y} / m_{X} = 0.2$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.3$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.4$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.5$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.6$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.7$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal.

Expected and observed limits at 95% CL on $\sigma \mathcal{B} (X \to Y(qq)q) \mathcal{A}$ for a cascade decay trijet resonance with $m_{Y} / m_{X} = 0.8$. Only 2016 data are used to derive limits below 2.0 TeV because of higher trigger thresholds in 2017 and 2018. Theoretical predictions of the benchmark are depicted with the red curve. A value of -1 in the table means that the corresponding theoretical prediction is not calculated for this signal

Cut flow table of the selection requirments for the $Z_{R}$ model scenarios. Values shown are absolute efficiencies, e.g., values shown in the column of 'Efficiency ($\Delta_{R}^{max} < 3.0$)' represent the cumulative efficiencies achieved through all selection requirements applied up to and including the current selection criterion.

Cut flow table of the selection requirments for the $G_{KK}$ model scenarios. Values shown are absolute efficiencies, e.g., values shown in the column of 'Efficiency ($\Delta_{R}^{max} < 3.0$)' represent the cumulative efficiencies achieved through all selection requirements applied up to and including the current selection criterion.

Cut flow table of the selection requirments for the $q^{*}$ model scenarios. Values shown are absolute efficiencies, e.g., values shown in the column of 'Efficiency ($\Delta_{R}^{max} < 3.0$)' represent the cumulative efficiencies achieved through all selection requirements applied up to and including the current selection criterion.

Predicted production cross section of the benchmark signal process $pp \to Z_{R} \to ggg$.

Predicted production cross section of the benchmark signal process $pp \to G_{KK} \to \varphi(gg)g$.

Predicted production cross section of the benchmark signal process $pp \to q^{*} \to V(qq)q$.

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$.

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 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.

Observed upper limits at 90% CL on the square of the kinetic mixing coefficient $\epsilon$ in the minimal model of a dark photon

Observed upper limits at 90% CL on the mixing angle $\theta_H$ for the type IV 2HDM+S scenario with fixed $tan\beta=0.5$

The combined efficiency of the dimuon scouting trigger and the MVA muon selection, averaged between 2017 and 2018, weighted by the integrated luminosity of each year. The trigger and selection efficiencies are also included in the table.

Theory cross-section times branching fraction to muons times acceptance for the dark photon and 2HDM+S models.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{1})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}^{\mathrm{add.}}_{2})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{2})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{add.}})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\Delta\mathrm{R}_{\mathrm{b}\mathrm{b}}^{\mathrm{avg}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}_{3})|$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}_{3})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}_{3})$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}_{3})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}_{4})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}_{4})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{b}}_{\mathrm{T}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{b}}_{\mathrm{T}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $|\eta(\mathrm{b}^{\mathrm{extra}}_{1})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{1})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}^{\mathrm{extra}}_{2})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{2})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{extra}})|$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $H^{\mathrm{j}}_{\mathrm{T}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{j}}_{\mathrm{T}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $\mathrm{m}_{\mathrm{b}\mathrm{b}}^{\mathrm{max}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $H^{\mathrm{light}}_{\mathrm{T}}$ in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $H^{\mathrm{light}}_{\mathrm{T}}$ in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space.

Normalized differential cross section of $N_{\mathrm{jets}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Normalized differential cross section of $N_{\mathrm{jets}}$ in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space.

Normalized differential cross section of $N_{\mathrm{b}}$ in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space.

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{add.}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\Delta\mathrm{R}_{\mathrm{b}\mathrm{b}}^{\mathrm{avg}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}_{4})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}_{4})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{extra}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $\mathrm{m}_{\mathrm{b}\mathrm{b}}^{\mathrm{max}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Correlation of parameters of interest in fit of $N_{\mathrm{jets}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $N_{\mathrm{jets}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Correlation of parameters of interest in fit of $N_{\mathrm{b}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{add.}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{add.}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{add.}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{add.}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\Delta\mathrm{R}_{\mathrm{b}\mathrm{b}}^{\mathrm{avg}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}_{3})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}_{3})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}_{4})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}_{4})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{b}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $|\Delta\phi(\mathrm{lj}^{\mathrm{extra}}_{1},\mathrm{b}_{\mathrm{soft}})|$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{1})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}^{\mathrm{extra}}_{2})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}^{\mathrm{extra}}_{2})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $|\eta(\mathrm{b}\mathrm{b}^{\mathrm{extra}})|$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\Delta\mathrm{R}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\mathrm{m}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{b}\mathrm{b}^{\mathrm{extra}})$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $p_{\mathrm{T}}(\mathrm{lj}^{\mathrm{extra}}_{1})$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{j}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $\mathrm{m}_{\mathrm{b}\mathrm{b}}^{\mathrm{max}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 6$ jets: $\geq 3 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $H^{\mathrm{light}}_{\mathrm{T}}$ observable in $\geq 7$ jets: $\geq 4 \mathrm{b}$, $\geq 3$ light phase space

Covariances of all nuisance parameters and POIs in fit of $N_{\mathrm{jets}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $N_{\mathrm{jets}}$ observable in $\geq 6$ jets: $\geq 4 \mathrm{b}$ phase space

Covariances of all nuisance parameters and POIs in fit of $N_{\mathrm{b}}$ observable in $\geq 5$ jets: $\geq 3 \mathrm{b}$ phase space

Postfit distributions of $S_\mathrm{T}^\mathrm{MET}$ in the $\mathrm{e}\mu$ channel of the $\geq$1b category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width. The last bin includes the overflow.

Postfit distributions of $S_\mathrm{T}^\mathrm{MET}$ in the $\ell\tau_\mathrm{h}$ channel of the 0b category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width. The last bin includes the overflow.

Postfit distributions of $S_\mathrm{T}^\mathrm{MET}$ in the $\ell\tau_\mathrm{h}$ channel of the $\geq$1b category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width. The last bin includes the overflow.

Postfit distributions of $S_\mathrm{T}^\mathrm{MET}$ in the $\tau_\mathrm{h}\tau_\mathrm{h}$ channel of the 0b category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width. The last bin includes the overflow.

Postfit distributions of $S_\mathrm{T}^\mathrm{MET}$ in the $\tau_\mathrm{h}\tau_\mathrm{h}$ channel of the $\geq$1b category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width. The last bin includes the overflow.

Postfit distributions of $\chi$ in the $\mathrm{e}\mu$ channel of the 0j, $400 < m_\mathrm{vis} < 600\,\mathrm{GeV}$ category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width.

Postfit distributions of $\chi$ in the $\mathrm{e}\mu$ channel of the 0j, $m_\mathrm{vis} > 600\,\mathrm{GeV}$ category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width.

Postfit distributions of $\chi$ in the $\ell\tau_\mathrm{h}$ channel of the 0j, $400 < m_\mathrm{vis} < 600\,\mathrm{GeV}$ category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width.

Postfit distributions of $\chi$ in the $\ell\tau_\mathrm{h}$ channel of the 0j, $m_\mathrm{vis} > 600\,\mathrm{GeV}$ category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width.

Postfit distributions of $\chi$ in the $\tau_\mathrm{h}\tau_\mathrm{h}$ channel of the 0j, $400 < m_\mathrm{vis} < 600\,\mathrm{GeV}$ category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width.

Postfit distributions of $\chi$ in the $\tau_\mathrm{h}\tau_\mathrm{h}$ channel of the 0j, $m_\mathrm{vis} > 600\,\mathrm{GeV}$ category for the combined 2016-2018 data set after a simultaneous fit of the background and vector LQ signal to the data. The number of events in each bin are divided by the respective bin width.

Exclusion limits on the total cross section of a scalar LQ as a function of LQ mass under the assumption of exclusive LQ couplings to b quarks and $\tau$ leptons.

Exclusion limits on the total cross section of a vector LQ as a function of LQ mass under the assumption of exclusive LQ couplings to b quarks and $\tau$ leptons.

Exclusion limits on the mass of a scalar LQ as a function of the coupling strength $\lambda$ under the assumption of exclusive LQ couplings to b quarks and $\tau$ leptons.

Exclusion limits on the mass of a vector LQ as a function of the coupling strength $\lambda$ under the assumption of exclusive LQ couplings to b quarks and $\tau$ leptons.

Exclusion limits on the coupling strength $\lambda$ of a scalar LQ as a function of the LQ mass under the assumption of exclusive LQ couplings to b quarks and $\tau$ leptons.

Exclusion limits on the coupling strength $\lambda$ of a vector LQ as a function of the LQ mass under the assumption of exclusive LQ couplings to b quarks and $\tau$ leptons.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet no-btag category with the BDT requirements selecting the most signal-like events. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet no-btag category with the BDT requirements selecting the most signal-like events. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet no-btag category with the BDT requirements selecting the most signal-like events. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to b quarks and $\tau$ leptons, using $\lambda = 1.5$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to light-flavor quarks and $\tau$ leptons, using $\lambda = 1.5$.

Expected upper limit at 95% CL on the coupling strength $\lambda$ of a scalar LQ to b quarks and $\tau$ leptons.

Expected upper limit at 95% CL on the coupling strength $\lambda$ of a scalar LQ to light-flavor quarks and $\tau$ leptons.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet btag category with the intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet no-btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet no-btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet btag category with the lowest BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet no-btag category with low BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet no-btag category with low BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet no-btag category with low BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\tau_\mathrm{h}$+jet no-btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the e+jet btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Observed and expected distributions of the collinear mass in the $\mu$+jet no-btag category with intermediate BDT requirement. The signal hypothesis corresponds to a scalar leptoquark coupled to u quarks and $\tau$ leptons with a coupling strength $\lambda$ equal to 1.5.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to b quarks and $\tau$ leptons, using $\lambda = 1.0$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to b quarks and $\tau$ leptons, using $\lambda = 2.0$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to b quarks and $\tau$ leptons, using $\lambda = 3.0$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to light-flavor quarks and $\tau$ leptons, using $\lambda = 0.5$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to light-flavor quarks and $\tau$ leptons, using $\lambda = 1.0$.

Expected and observed upper limits at 95% CL on the scalar lepton-induced LQ production cross section times branching fraction for a LQ coupled to light-flavor quarks and $\tau$ leptons, using $\lambda = 2.0$.

Observed and expected distributions of the jet multiplicity in the e+jet btag control region.

Observed and expected distributions of the jet multiplicity in the $\mu$+jet btag control region.

Inclusive cross section predictions up to QCD NNLO accuracy using the NNPDF3.1 PDF set

Measured production cross sections ratio $\sigma(W^+ + \overline{c})$ / $\sigma(W^- + c)$

Predictions for the cross sections ratio $\sigma(W^+ + \overline{c})$ / $\sigma(W^- + c)$ at QCD NLO accuracy from MCFM using different PDF sets

Predictions for the cross sections ratio $\sigma(W^+ + \overline{c})$ / $\sigma(W^- + c)$ up to QCD NNLO accuracy using the NNPDF3.1 PDF set

Measured diferential cross sections $\sigma(W^- + c) + \sigma(W^+ + \overline{c})$ as a function of the absolute value of the pseudorapidity of the lepton from the W decay, unfolded to the particle level.

Measured diferential cross sections as a function of the transverse momentum of the lepton from the W decay, unfolded to the particle level.

Measured diferential cross sections $\sigma(W^- + c) + \sigma(W^+ + \overline{c})$ as a function of the absolute value of the pseudorapidity of the lepton from the W decay, unfolded to the parton level.

Measured diferential cross sections as a function of the transverse momentum of the lepton from the W decay, unfolded to the parton level.

Measured diferential cross sections ratio $R=\sigma(W^+ + \overline{c}) / \sigma(W^- + c)$ as a function of the absolute value of the pseudorapidity of the lepton from the W decay.

Measured diferential cross sections ratio $R=\sigma(W^+ + \overline{c}) / \sigma(W^- + c)$ as a function of the transverse momentum of the lepton from the W decay.

Two dimensional observed likelihood profile scan as a function of mu_ggH and mu_qqH

Observed $2 \Delta NLL$ of a 2D scan for $c^{-}_{\varphi Q}$ and $c_{\varphi t}$ with the other WCs profiled.

Observed $2 \Delta NLL$ of a 2D scan for $c_{tG}$ and $c_{t\varphi}$ with the other WCs profiled.

Observed $2 \Delta NLL$ of a 2D scan for $c^{1}_{Qt}$ and $c^{8}_{Qt}$ with the other WCs profiled.

Observed $2 \Delta NLL$ of a 2D scan for $c^{1}_{QQ}$ and $c^{8}_{Qt}$ with the other WCs profiled.

Observed $2 \Delta NLL$ of a 2D scan for $c^{1}_{Qt}$ and $c^{1}_{tt}$ with the other WCs profiled.

Observed $2 \Delta NLL$ of a 2D scan for $c^{1}_{QQ}$ and $c^{1}_{Qt}$ with the other WCs profiled.

The nuclear modification factor $\mathrm{R_{AA}}$ versus $p_{\mathrm{T}}$ for prompt $\Lambda^+_c$ production in centrality regions of 0-90, 0-10, 10-30, 30-50 and 50-90% in PbPb collisions. The global uncertainty includes the uncertainties for the luminosity of pp collisions, number of MB events in PbPb collisions, and tracking efficiency. The global uncertainty for $\mathrm{R_{AA}}$ is 16.5%.

The ratio of the production cross sections of prompt $\Lambda^+_c$ to prompt $\mathrm{D_0}$ versus $p_{\mathrm{T}}$ from pp collisions. The global normalization uncertainty is 6.6%.

The ratio of the $\mathrm{T_{AA}}$-scaled prompt $\Lambda^+_c$ yields to prompt $\mathrm{D_0}$ versus $p_{\mathrm{T}}$ from PbPb collisions for 0-90 and 0-10% centrality classes. The global normalization uncertainty is 7.3%.

Feynman diagrams of $\mathrm{Z'}\to\mu^{-}\mu^{+}$ with a $\mathrm{Z'}$ boson produced via $\mathrm{s}\overline{\mathrm{b}}\to\mathrm{Z'}$, with one $\mathrm{b}$ quark and one $\overline{\mathrm{s}}$ quark in the final state.

Distribution of $\min(m_{\mu\mathrm{b}})$ as obtained from simulation in events with $N_{\mathrm{b}} \geq 1$ passing all the other selection requirements. In this search, we require $\min(m_{\mu\mathrm{b}})>175~\mathrm{GeV}$. The stacked histogram displays the expected distribution from the simulation of the SM backgrounds, while the overlaid open histograms illustrate the size and shape of the $\mathrm{Z'}$ contribution from the LFU model described in Eq. (1), for several $\mathrm{Z'}$ mass hypotheses. For illustrative purposes, we choose couplings $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$. The contribution of background processes other than DY and $\mathrm{t}\overline{\mathrm{t}}$ is so small that it is only barely visible at the bottom of the stacked histogram. The hatched region indicates the statistical uncertainty arising from the limited size of the SM simulated samples. Histograms are normalized to unit area.

Distributions of $m_{\mu\mu}$ in the $N_{\mathrm{b}}=1$ event category. The stacked histogram displays the expected distribution from the SM background simulation. The overlaid open distributions illustrate the $\mathrm{Z'}$ contribution from the LFU model at $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$ for a variety of $\mathrm{Z'}$ mass hypotheses. The observed data are shown as black points with statistical error bars. The hatched region indicates the statistical uncertainty arising from the limited size of the SM simulated samples. The size of the bins increases as a function of $m_{\mu\mu}$. In extracting the results of the search, the background is estimated directly from data, so the SM background simulation is only illustrative in these distributions.

Distributions of $m_{\mu\mu}$ in the $N_{\mathrm{b}}\geq 2$ event category. The stacked histogram displays the expected distribution from the SM background simulation. The overlaid open distributions illustrate the $\mathrm{Z'}$ contribution from the LFU model at $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$ for a variety of $\mathrm{Z'}$ mass hypotheses. The observed data are shown as black points with statistical error bars. The hatched region indicates the statistical uncertainty arising from the limited size of the SM simulated samples. The size of the bins increases as a function of $m_{\mu\mu}$. In extracting the results of the search, the background is estimated directly from data, so the SM background simulation is only illustrative in these distributions.

Invariant mass $m_{\mu\mu}$ distributions in the $N_{\mathrm{b}} = 1$ category, shown together with the corresponding selected background functional forms used as input to the $discrete~profiling$ method when probing the $m_{\mathrm{Z'}}=500~\mathrm{GeV}$ hypothesis. The expected signal distribution for the LFU model described in Eq. (1), with couplings $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$, is overlaid. The displayed mass range corresponds to the fit window used for this $m_{\mathrm{Z'}}$ hypothesis, which is $\pm 10\, \sigma_{\mathrm{mass}}$ around the probed $m_{\mathrm{Z'}}$ value. While the likelihood fits are performed on unbinned data, here we present the data in binned histograms with binning chosen to reflect the size of $\sigma_{\mathrm{mass}}$.

Invariant mass $m_{\mu\mu}$ distributions in the $N_{\mathrm{b}} \geq 2$ category, shown together with the corresponding selected background functional forms used as input to the $discrete~profiling$ method when probing the $m_{\mathrm{Z'}}=500~\mathrm{GeV}$ hypothesis. The expected signal distribution for the LFU model described in Eq. (1), with couplings $|g_{\ell}| = |g_{\nu}| = |g_{\mathrm{b}}| = 0.03$ and $\delta_{\mathrm{bs}}=0$, is overlaid. The displayed mass range corresponds to the fit window used for this $m_{\mathrm{Z'}}$ hypothesis, which is $\pm 10\, \sigma_{\mathrm{mass}}$ around the probed $m_{\mathrm{Z'}}$ value. While the likelihood fits are performed on unbinned data, here we present the data in binned histograms with binning chosen to reflect the size of $\sigma_{\mathrm{mass}}$.

Exclusion limits at $95\%$ CL on the number of selected BSM events with $N_{\mathrm{b}} \geq 1$ as functions of $m_{\mathrm{Z'}}$ for $f_{2\mathrm{b}}=0$. The quantity $f_{2\mathrm{b}}$ is the fraction of BSM events passing the analysis selection that have at least two $\mathrm{b}$ quark jets. The solid black (dashed red) curve represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing $68~(95)\%$ of the distribution of limits expected under the background-only hypothesis.

Exclusion limits at $95\%$ CL on the number of selected BSM events with $N_{\mathrm{b}} \geq 1$ as functions of $m_{\mathrm{Z'}}$ for $f_{2\mathrm{b}}=0.25$. The quantity $f_{2\mathrm{b}}$ is the fraction of BSM events passing the analysis selection that have at least two $\mathrm{b}$ quark jets. The solid black (dashed red) curve represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing $68~(95)\%$ of the distribution of limits expected under the background-only hypothesis.

Exclusion limits at $95\%$ CL on the number of selected BSM events with $N_{\mathrm{b}} \geq 1$ as functions of $m_{\mathrm{Z'}}$ for $f_{2\mathrm{b}}=0.75$. The quantity $f_{2\mathrm{b}}$ is the fraction of BSM events passing the analysis selection that have at least two $\mathrm{b}$ quark jets. The solid black (dashed red) curve represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing $68~(95)\%$ of the distribution of limits expected under the background-only hypothesis.

Exclusion limits at $95\%$ CL on the number of selected BSM events with $N_{\mathrm{b}} \geq 1$ as functions of $m_{\mathrm{Z'}}$ for $f_{2\mathrm{b}}=1$. The quantity $f_{2\mathrm{b}}$ is the fraction of BSM events passing the analysis selection that have at least two $\mathrm{b}$ quark jets. The solid black (dashed red) curve represents the observed (median expected) exclusion. The inner green (outer yellow) band indicates the region containing $68~(95)\%$ of the distribution of limits expected under the background-only hypothesis.

Observed (solid) and median expected (dashed) exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution, marked by the dotted curves. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $\delta_{\mathrm{bs}}=0$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed (solid) and median expected (dashed) exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution, marked by the dotted curves. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=350~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=700~\mathrm{GeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The exclusion limits are given up to coupling values at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid. The enclosed regions are excluded.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{b}}|$--$|g_{\ell}|$ plane for the LFU model, with $|\delta_{\mathrm{bs}}|=0.25$ and $|g_{\nu}|=|g_{\ell}|$, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width is equal to half of the $\mu\mu$ invariant mass resolution. Beyond these coupling values, the narrow width approximation intrinsic to the search strategy is not considered valid.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$. The solid black (dashed red) curves represent the observed (median expected) exclusions. The dotted curves denote the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and dotted curves is (expected to be) excluded. The dotted curve for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$ lies beyond the displayed $|g_{\mathrm{Z'}}|$ range and is, therefore, not shown. The shaded blue area represents the region preferred from the global fit in JHEP 04 (2023) 033 at $95\%$ CL. The region above the green dash-dotted curve is incompatible at $95\%$ CL with the measurement of the mass difference between the mass eigenstates of the neutral $\mathrm{B}_\mathrm{s}$ mesons.

Observed exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$.

Expected exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=500~\mathrm{GeV}$.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The solid black (dashed red) curves represent the observed (median expected) exclusions. The dotted curves denote the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and dotted curves is (expected to be) excluded. The shaded blue area represents the region preferred from the global fit in JHEP 04 (2023) 033 at $95\%$ CL. The region above the green dash-dotted curve is incompatible at $95\%$ CL with the measurement of the mass difference between the mass eigenstates of the neutral $\mathrm{B}_\mathrm{s}$ mesons.

Observed exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$.

Expected exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The solid black (dashed red) curves represent the observed (median expected) exclusions. The dotted curves denote the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and dotted curves is (expected to be) excluded. The shaded blue area represents the region preferred from the global fit in JHEP 04 (2023) 033 at $95\%$ CL. The region above the green dash-dotted curve is incompatible at $95\%$ CL with the measurement of the mass difference between the mass eigenstates of the neutral $\mathrm{B}_\mathrm{s}$ mesons.

Observed exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$.

Expected exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=1.5~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The solid black (dashed red) curves represent the observed (median expected) exclusions. The dotted curves denote the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and dotted curves is (expected to be) excluded. The shaded blue area represents the region preferred from the global fit in JHEP 04 (2023) 033 at $95\%$ CL. The region above the green dash-dotted curve is incompatible at $95\%$ CL with the measurement of the mass difference between the mass eigenstates of the neutral $\mathrm{B}_\mathrm{s}$ mesons.

Observed exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$.

Expected exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$.

Exclusion limits at $95\%$ CL in the $|\theta_{23}|$--$|g_{\mathrm{Z'}}|$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for $m_{\mathrm{Z'}}=2~\mathrm{TeV}$. The coupling values are shown at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{Z'}}|$--$m_{\mathrm{Z'}}$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for a fixed value of $\theta_{23}=0$. The solid black (dashed red) curve represents the observed (median expected) exclusion. The dotted curve denotes the coupling values at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid. The region enclosed between the solid black (dashed red) and the dotted curves is (expected to be) excluded. The shaded blue area represents the $|g_{\mathrm{Z'}}|$ range preferred from a global fit in JHEP 04 (2023) 033 at $95\%$ CL.

Observed exclusion limits at $95\%$ CL in the $|g_{\mathrm{Z'}}|$--$m_{\mathrm{Z'}}$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for a fixed value of $\theta_{23}=0$.

Expected exclusion limits at $95\%$ CL in the $|g_{\mathrm{Z'}}|$--$m_{\mathrm{Z'}}$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for a fixed value of $\theta_{23}=0$.

Exclusion limits at $95\%$ CL in the $|g_{\mathrm{Z'}}|$--$m_{\mathrm{Z'}}$ plane for the $B_{3}\!-\!L_{2}$ model from JHEP 04 (2023) 033, for a fixed value of $\theta_{23}=0$. The coupling values are shown at which the $\mathrm{Z'}$ width equals one half of the $\mu\mu$ invariant mass resolution. For larger values of the couplings, the narrow width approximation intrinsic to the search strategy is not considered valid.

Reconstructed mass of the Z' candidate in CR1, which consists of flavor sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Pre-fit (post-fit) results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR1, which consists of flavor sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Pre-fit (post-fit) results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR1, which consists of flavor sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Pre-fit (post-fit) results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR2, which consists of the dilepton mass sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Post-fit dielectron (dimuon) channel results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR2, which consists of the dilepton mass sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Post-fit dielectron (dimuon) channel results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR2, which consists of the dilepton mass sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Post-fit dielectron (dimuon) channel results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR2, which consists of the dilepton mass sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Post-fit dielectron (dimuon) channel results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR2, which consists of the dilepton mass sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Post-fit dielectron (dimuon) channel results are shown on the left (right). The fitting procedure is described in Section 7.

Reconstructed mass of the Z' candidate in CR2, which consists of the dilepton mass sidebands of SR1 (upper), SR2 (middle), and SR3 (lower) regions. Post-fit dielectron (dimuon) channel results are shown on the left (right). The fitting procedure is described in Section 7.

Distributions of the reconstructed Z' candidate mass in SR1 (upper row), SR2 (middle row) and SR3 (lower row). Left (right) column corresponds to the dielectron (dimuon) channel. Signal samples with ($m_{Z'}$, $m_{N}$) = (4000 GeV, 200 GeV)$ and ($m_{Z'}$, $m_{N}$) = (4000 GeV, 1200 GeV)$ are plotted together with a cross section times branching fraction of 1 fb

Distributions of the reconstructed Z' candidate mass in SR1 (upper row), SR2 (middle row) and SR3 (lower row). Left (right) column corresponds to the dielectron (dimuon) channel. Signal samples with ($m_{Z'}$, $m_{N}$) = (4000 GeV, 200 GeV)$ and ($m_{Z'}$, $m_{N}$) = (4000 GeV, 1200 GeV)$ are plotted together with a cross section times branching fraction of 1 fb

Distributions of the reconstructed Z' candidate mass in SR1 (upper row), SR2 (middle row) and SR3 (lower row). Left (right) column corresponds to the dielectron (dimuon) channel. Signal samples with ($m_{Z'}$, $m_{N}$) = (4000 GeV, 200 GeV)$ and ($m_{Z'}$, $m_{N}$) = (4000 GeV, 1200 GeV)$ are plotted together with a cross section times branching fraction of 1 fb

Distributions of the reconstructed Z' candidate mass in SR1 (upper row), SR2 (middle row) and SR3 (lower row). Left (right) column corresponds to the dielectron (dimuon) channel. Signal samples with ($m_{Z'}$, $m_{N}$) = (4000 GeV, 200 GeV)$ and ($m_{Z'}$, $m_{N}$) = (4000 GeV, 1200 GeV)$ are plotted together with a cross section times branching fraction of 1 fb

Distributions of the reconstructed Z' candidate mass in SR1 (upper row), SR2 (middle row) and SR3 (lower row). Left (right) column corresponds to the dielectron (dimuon) channel. Signal samples with ($m_{Z'}$, $m_{N}$) = (4000 GeV, 200 GeV)$ and ($m_{Z'}$, $m_{N}$) = (4000 GeV, 1200 GeV)$ are plotted together with a cross section times branching fraction of 1 fb

Distributions of the reconstructed Z' candidate mass in SR1 (upper row), SR2 (middle row) and SR3 (lower row). Left (right) column corresponds to the dielectron (dimuon) channel. Signal samples with ($m_{Z'}$, $m_{N}$) = (4000 GeV, 200 GeV)$ and ($m_{Z'}$, $m_{N}$) = (4000 GeV, 1200 GeV)$ are plotted together with a cross section times branching fraction of 1 fb

The observed and expected 95\% CL upper limits on the product of \Zp boson cross section times branching fraction are shown for the case of ${m}_{N}$ = ${m}_{Z'}/4$ (upper row) and ${m}_{N}$ = 100 GeV (lower row) for dielectron (left column) and dimuon (right column) channels. The green and yellow bands indicate the 68\% and 95\% CL regions around the expected limit.

The observed and expected 95\% CL upper limits on the product of \Zp boson cross section times branching fraction are shown for the case of ${m}_{N}$ = ${m}_{Z'}/4$ (upper row) and ${m}_{N}$ = 100 GeV (lower row) for dielectron (left column) and dimuon (right column) channels. The green and yellow bands indicate the 68\% and 95\% CL regions around the expected limit.

The observed and expected 95\% CL upper limits on the product of \Zp boson cross section times branching fraction are shown for the case of ${m}_{N}$ = ${m}_{Z'}/4$ (upper row) and ${m}_{N}$ = 100 GeV (lower row) for dielectron (left column) and dimuon (right column) channels. The green and yellow bands indicate the 68\% and 95\% CL regions around the expected limit.

The observed and expected 95\% CL upper limits on the product of \Zp boson cross section times branching fraction are shown for the case of ${m}_{N}$ = ${m}_{Z'}/4$ (upper row) and ${m}_{N}$ = 100 GeV (lower row) for dielectron (left column) and dimuon (right column) channels. The green and yellow bands indicate the 68\% and 95\% CL regions around the expected limit.

Observed and expected exclusion regions at 95% CL in the 2D phase space of ${m}_{Z'}$ vs. ${m}_{N}$ for dielectron (left) and dimuon (right) channels. One standard deviation (s.d.) limits are also shown. Both opposite- and same-sign (OS and SS) dileptons from decays of heavy neutrino pairs are considered.

Observed and expected exclusion regions at 95% CL in the 2D phase space of ${m}_{Z'}$ vs. ${m}_{N}$ for dielectron (left) and dimuon (right) channels. One standard deviation (s.d.) limits are also shown. Both opposite- and same-sign (OS and SS) dileptons from decays of heavy neutrino pairs are considered.

The Levy scale parameter $R$ versus the transverse mass $m_{\mathrm{T}}$ in different centrality classes for negatively and positively charged hadron pairs.

The $1/R^2$ distribution vs. transverse mass $m_{\mathrm{T}}$ in different centrality classes for negatively and positively charged hadron pairs.

The linear-fit parameter $A$ versus $\langle N_{\text{part}}\rangle$ for positively and negatively charged hadron pairs. The bin range on $\langle N_{\text{part}}\rangle$ represents the 68% confidence interval associated to its uncertainty.

The linear-fit parameter $B$ versus $\langle N_{\text{part}}\rangle$ for positively and negatively charged hadron pairs. The bin range on $\langle N_{\text{part}}\rangle$ represents the 68% confidence interval associated to its uncertainty.

The Levy scale parameter $R$ versus $\langle N_{\text{part}}\rangle^{1/3}$ in different $m_{\mathrm{T}}$ classes for negatively and positively charged hadron pairs. The bin range on $\langle N_{\text{part}}\rangle^{1/3}$ represents the 68% confidence interval associated to its uncertainty.

The Levy stability index $\alpha$ versus the transverse mass $m_{\mathrm{T}}$ in different centrality classes for negatively and positively charged hadron pairs.

The average Levy stability index $\langle \alpha \rangle$ versus $\langle N_{\text{part}}\rangle$ for both positively and negatively charged hadron pairs. The bin range on $\langle N_{\text{part}}\rangle$ represents the 68% confidence interval associated to its uncertainty.

The correlation strength $\lambda$ versus the transverse mass $m_{\mathrm{T}}$ in different centrality classes for negatively and positively charged hadron pairs.

The correlation strength (re-scaled with the square of the pion fraction) $\lambda ^*$ versus the transverse mass $m_{\mathrm{T}}$ in different centrality classes for positively and negatively charged hadron pairs.

The $m_{e\mu}$ distribution of the observed data is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in red solid line, and the background component of the fit in red dotted line. Events and fit in each category are weighted by $S/(S+B)$. The one and two standard deviation uncertainty bands of the background component are shown in green and yellow. The lower panel shows the residuals after subtracting the background component of the fit from data.

Observed (expected) 95% confidence level upper limits on $\sigma(p p \to X(146) \to e \mu)$ for each individual analysis category (as shown in the left axis label) and for the combination of all analysis categories assuming the relative SM-like production cross sections of the ggH and VBF production modes.

The $m_{e\mu}$ distribution of the observed data in the category ggH cat 0 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category ggH cat 1 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category ggH cat 2 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category ggH cat 3 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category VBF cat 0 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

The $m_{e\mu}$ distribution of the observed data in the category VBF cat 1 is shown, with the $S+B$ fit at $m_{X}=146\mathrm{GeV}$ in solid red line, and the background component of the fit in dotted red line. The green and yellow bands represent the one and two standard deviations of the background component. The lower panel shows the residuals after subtracting the background component of the fit from data.

Best fit of $\sigma(p p \to X(146) \to e \mu)$ for each analysis category (black point) compared to the best fit for the combination of all analysis categories (blue line). The one standard deviation uncertainty of the per-category best fit is shown in red line and the one standard deviation uncertainty of the combined best fit is shown in green band.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.0.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.1.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.2.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.3.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.4.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.5.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.6.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.7.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.8.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 0.9.

The observed and expected 95% CL upper limits on $\sigma_{ggH,f_{VBF}}(p p \to X \to e \mu)$ and $\sigma_{VBF,f_{VBF}}(p p \to X \to e \mu)$ with $f_{VBF}$ = 1.0.

The $v_{3}$ values as a function of centrality for prompt and nonprompt J/$\psi$ mesons.The kinematic range is 6.5 $< p_{\text{T}} <$ 50 GeV/c and $|y| < 2.4$.

The $v_{2}$ values of prompt J/$\psi$ mesons as a function of $p_{\text{T}}$ in the 10–60% centrality range. The results for 3 $< p_{\text{T}} <$ 6.5 and 6.5 $< p_{\text{T}} <$ 50 GeV/c are studied in the rapidity range of 1.6 $< |y| <$ 2.4 and $|y| <$ 2.4, respectively.

The $v_{2}$ values as functions of $p_{\text{T}}$ for prompt J/$\psi$ and prompt $\psi$(2S) mesons in the 10–60% centrality range. The results for $p_{\text{T}}$ < 6.5 and $p_{\text{T}}$ > 6.5 GeV/c are studied in the rapidity range of 1.6 $< |y| <$ 2.4 and $|y| <$ 2.4, respectively.

The $v_{2}$ values as functions of centrality for prompt J/$\psi$ and prompt $\psi$(2S) mesons. The kinematic range is 6.5 $< p_{\text{T}} <$ 50 GeV/c and $|y| <$ 2.4.

The $v_{3}$ values as functions of $p_{\text{T}}$ for prompt J/$\psi$ and prompt $\psi$(2S) mesons in the 10–60% centrality range. The results for $p_{\text{T}}$ < 6.5 and $p_{\text{T}}$ > 6.5 GeV/c are studied in the rapidity range of 1.6 $< |y| <$ 2.4 and $|y| <$ 2.4, respectively.

The $v_{3}$ values as functions of centrality for prompt J/$\psi$ and prompt $\psi$(2S) mesons. The kinematic range is 6.5 $< p_{\text{T}} <$ 50 GeV/c and $|y| <$ 2.4.