Dilepton mass spectra for the high mass $W\phi($ee$)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $W\phi($ee$)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $W\phi($ee$)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $Z\phi($ee$)$ SR event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $Z\phi($ee$)$ SR event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $t\bar{t}\phi($ee$)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $t\bar{t}\phi($ee$)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $t\bar{t}\phi($ee$)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $t\bar{t}\phi($ee$)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $t\bar{t}\phi($ee$)$ SR3 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $t\bar{t}\phi($ee$)$ SR3 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $W\phi(\mu\mu)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $W\phi(\mu\mu)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $W\phi(\mu\mu)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $W\phi(\mu\mu)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $Z\phi(\mu\mu)$ SR event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $Z\phi(\mu\mu)$ SR event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $t\bar{t}\phi(\mu\mu)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $t\bar{t}\phi(\mu\mu)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $t\bar{t}\phi(\mu\mu)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $t\bar{t}\phi(\mu\mu)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the low mass $t\bar{t}\phi(\mu\mu)$ SR3 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the high mass $t\bar{t}\phi(\mu\mu)$ SR3 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $W\phi(\tau\tau)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $Z\phi(\tau\tau)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $W\phi(\tau\tau)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $Z\phi(\tau\tau)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $W\phi(\tau\tau)$ SR3 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $Z\phi(\tau\tau)$ SR3 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $t\bar{t}\phi(\tau\tau)$ SR1 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $t\bar{t}\phi(\tau\tau)$ SR2 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $t\bar{t}\phi(\tau\tau)$ SR3 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $t\bar{t}\phi(\tau\tau)$ SR4 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $t\bar{t}\phi(\tau\tau)$ SR5 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $t\bar{t}\phi(\tau\tau)$ SR6 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

Dilepton mass spectra for the $t\bar{t}\phi(\tau\tau)$ SR7 event selections for the combined 2016–2018 data set. The lower panel shows the ratio of observed events to the total expected SM background prediction (Obs/Exp), and the gray band represents the sum of statistical and systematic uncertainties in the background prediction. The rightmost bin contains the overflow events in each distribution. The expected background distributions and the uncertainties are shown after the data is fit under the background-only hypothesis. For illustration, two example signal hypotheses for the production and decay of a scalar and a pseudoscalar $\phi$ boson are shown, and their masses (in units of GeV) are indicated in the legend. The signals are normalized to the product of the cross section and branching fraction of 10 fb.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $W\phi$ signal with scalar couplings in the ee decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $W\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $W\phi$ signal with pseudoscalar couplings in the ee decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $W\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $W\phi$ signal with scalar couplings in the $\mu\mu$ decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $W\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $W\phi$ signal with pseudoscalar couplings in the $\mu\mu$ decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $W\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $W\phi$ signal with scalar couplings in the $\tau\tau$ decay scenario. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $W\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $W\phi$ signal with pseudoscalar couplings in the $\tau\tau$ decay scenario. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $W\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $Z\phi$ signal with scalar couplings in the ee decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $Z\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $Z\phi$ signal with pseudoscalar couplings in the ee decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $Z\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $Z\phi$ signal with scalar couplings in the $\mu\mu$ decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $Z\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $Z\phi$ signal with pseudoscalar couplings in the $\mu\mu$ decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $Z\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $Z\phi$ signal with scalar couplings in the $\tau\tau$ decay scenario. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $Z\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $Z\phi$ signal with pseudoscalar couplings in the $\tau\tau$ decay scenario. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $Z\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $W\phi$ signal with H-like production in the ee decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $W\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $Z\phi$ signal with H-like production in the ee decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $Z\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $W\phi$ signal with H-like production in the $\mu\mu$ decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $W\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $Z\phi$ signal with H-like production in the $\mu\mu$ decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $Z\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $W\phi$ signal with H-like production in the $\tau\tau$ decay scenario. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $W\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $Z\phi$ signal with H-like production in the $\tau\tau$ decay scenario. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $Z\phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal with scalar couplings in the ee decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal with pseudoscalar couplings in the ee decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal with scalar couplings in the $\mu\mu$ decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal with pseudoscalar couplings in the $\mu\mu$ decay scenario. The vertical gray band indicates the mass region not considered in the analysis. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal with scalar couplings in the $\tau\tau$ decay scenario. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal.

The 95% confidence level upper limits on the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal with pseudoscalar couplings in the $\tau\tau$ decay scenario. The red line is the theoretical prediction for the product of the production cross section and branching fraction of the $t\bar{t} \phi$ signal.

The 95% confidence level upper limits on $g^2_{tS}$ for the dilaton-like $t\bar{t} \phi$ signal model. Masses of the $\phi$ boson above 300 GeV are not probed for the dilaton-like signal model as the $\phi$ branching fraction into top quark-antiquark pairs becomes nonnegligible.

The 95% confidence level upper limits on $g^2_{tPS}$ for the axion-like $t\bar{t} \phi$ signal model. Masses of the $\phi$ boson above 300 GeV are not probed for the axion-like signal model as the $\phi$ branching fraction into top quark-antiquark pairs becomes nonnegligible.

The 95% confidence level upper limits on the product of $sin^2 \theta$ and branching fraction for the H-like production of X$\phi \rightarrow$ ee. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level upper limits on the product of $sin^2 \theta$ and branching fraction for the H-like production of X$\phi \rightarrow \mu\mu$. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level upper limits on $sin^2 \theta$ for the H-like production and decay of X$\phi$ signal model.

Cross section in units of pb for the W$\phi$, Z$\phi$, and $t\bar{t}\phi$ signals as a function of the $\phi$ boson mass in GeV. All cross sections are inclusive of all W, Z, $t\bar{t}$ and $\phi$ decay modes.

The 95% confidence level expected and observed upper limits on the product of $g^{2}_{tS}$ and $\bf{\it{B}}(\phi \rightarrow $ee$)$ of the $t\bar{t} \phi$ signal with scalar couplings, where $g_{tS}$ denotes the coupling of the $\phi$ boson to the top quark and $\bf{\it{B}}(\phi \rightarrow $ee$)$ is the branching fraction of the $\phi$ boson into an electron pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $g^{2}_{tS}$ and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ of the $t\bar{t} \phi$ signal with scalar couplings, where $g_{tS}$ denotes the coupling of the $\phi$ boson to the top quark and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ is the branching fraction of the $\phi$ boson into a muon pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $g^{2}_{tS}$ and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ of the $t\bar{t} \phi$ signal with scalar couplings, where $g_{tS}$ denotes the coupling of the $\phi$ boson to the top quark and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ is the branching fraction of the $\phi$ boson into a tau pair.

The 95% confidence level expected and observed upper limits on the product of $g^{2}_{tPS}$ and $\bf{\it{B}}(\phi \rightarrow $ee$)$ of the $t\bar{t} \phi$ signal with pseudoscalar couplings, where $g_{tPS}$ denotes the coupling of the $\phi$ boson to the top quark and $\bf{\it{B}}(\phi \rightarrow $ee$)$ is the branching fraction of the $\phi$ boson into an electron pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $g^{2}_{tPS}$ and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ of the $t\bar{t} \phi$ signal with pseudoscalar couplings, where $g_{tPS}$ denotes the coupling of the $\phi$ boson to the top quark and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ is the branching fraction of the $\phi$ boson into a muon pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $g^{2}_{tPS}$ and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ of the $t\bar{t} \phi$ signal with pseudoscalar couplings, where $g_{tPS}$ denotes the coupling of the $\phi$ boson to the top quark and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ is the branching fraction of the $\phi$ boson into a tau pair.

The 95% confidence level expected and observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow $ee$)$ of the $t\bar{t} \phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow $ee$)$ is the branching fraction of the $\phi$ boson into an electron pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ of the $t\bar{t} \phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ is the branching fraction of the $\phi$ boson into a muon pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ of the $t\bar{t} \phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ is the branching fraction of the $\phi$ boson into a tau pair.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{S}$ and $\bf{\it{B}}(\phi \rightarrow $ee$)$ of the $W\phi$ signal with scalar couplings, where $\Lambda_{S}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $ee$)$ is the branching fraction of the $\phi$ boson into an electron pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{S}$ and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ of the $W\phi$ signal with scalar couplings, where $\Lambda_{S}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ is the branching fraction of the $\phi$ boson into a muon pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{S}$ and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ of the $W\phi$ signal with scalar couplings, where $\Lambda_{S}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ is the branching fraction of the $\phi$ boson into a tau pair.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{PS}$ and $\bf{\it{B}}(\phi \rightarrow $ee$)$ of the $W\phi$ signal with pseudoscalar couplings, where $\Lambda_{PS}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $ee$)$ is the branching fraction of the $\phi$ boson into an electron pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{PS}$ and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ of the $W\phi$ signal with pseudoscalar couplings, where $\Lambda_{PS}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ is the branching fraction of the $\phi$ boson into a muon pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{PS}$ and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ of the $W\phi$ signal with pseudoscalar couplings, where $\Lambda_{PS}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ is the branching fraction of the $\phi$ boson into a tau pair.

The 95% confidence level expected and observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow $ee$)$ of the $W\phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow $ee$)$ is the branching fraction of the $\phi$ boson into an electron pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ of the $W\phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ is the branching fraction of the $\phi$ boson into a muon pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ of the $W\phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ is the branching fraction of the $\phi$ boson into a tau pair.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{S}$ and $\bf{\it{B}}(\phi \rightarrow $ee$)$ of the $Z\phi$ signal with scalar couplings, where $\Lambda_{S}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $ee$)$ is the branching fraction of the $\phi$ boson into an electron pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{S}$ and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ of the $Z\phi$ signal with scalar couplings, where $\Lambda_{S}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ is the branching fraction of the $\phi$ boson into a muon pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{S}$ and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ of the $Z\phi$ signal with scalar couplings, where $\Lambda_{S}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ is the branching fraction of the $\phi$ boson into a tau pair.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{PS}$ and $\bf{\it{B}}(\phi \rightarrow $ee$)$ of the $Z\phi$ signal with pseudoscalar couplings, where $\Lambda_{PS}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $ee$)$ is the branching fraction of the $\phi$ boson into an electron pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{PS}$ and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ of the $Z\phi$ signal with pseudoscalar couplings, where $\Lambda_{PS}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ is the branching fraction of the $\phi$ boson into a muon pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $\Lambda^{-2}_{PS}$ and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ of the $Z\phi$ signal with pseudoscalar couplings, where $\Lambda_{PS}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ is the branching fraction of the $\phi$ boson into a tau pair.

The 95% confidence level expected and observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow $ee$)$ of the $Z\phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow $ee$)$ is the branching fraction of the $\phi$ boson into an electron pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ of the $Z\phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow $$\mu\mu$$)$ is the branching fraction of the $\phi$ boson into a muon pair. The vertical gray band indicates the mass region not considered in the analysis.

The 95% confidence level expected and observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ of the $Z\phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow $$\tau\tau$$)$ is the branching fraction of the $\phi$ boson into a tau pair.

The 95% confidence level observed upper limits on the product of $\sigma$($t \bar{t} \phi$) and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $t \bar{t} \phi$ signal with scalar couplings, where $\sigma$ denotes the production cross section and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The red dash-dotted line is the theoretical prediction for $\sigma\bf{\it{B}}$ of the $t \bar{t} \phi$ signal. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\sigma$($t \bar{t} \phi$) and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $t \bar{t} \phi$ signal with pseudoscalar couplings, where $\sigma$ denotes the production cross section and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The red dash-dotted line is the theoretical prediction for $\sigma\bf{\it{B}}$ of the $t \bar{t} \phi$ signal. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\sigma$($W\phi$) and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $W\phi$ signal with scalar couplings, where $\sigma$ denotes the production cross section and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The red dash-dotted line is the theoretical prediction for $\sigma\bf{\it{B}}$ of the $W\phi$ signal. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\sigma$($W\phi$) and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $W\phi$ signal with pseudoscalar couplings, where $\sigma$ denotes the production cross section and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The red dash-dotted line is the theoretical prediction for $\sigma\bf{\it{B}}$ of the $W\phi$ signal. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\sigma$($W\phi$) and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $W\phi$ signal with H-like production, where $\sigma$ denotes the production cross section and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The red dash-dotted line is the theoretical prediction for $\sigma\bf{\it{B}}$ of the $W\phi$ signal. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\sigma$($Z\phi$) and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $Z\phi$ signal with scalar couplings, where $\sigma$ denotes the production cross section and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The red dash-dotted line is the theoretical prediction for $\sigma\bf{\it{B}}$ of the $Z\phi$ signal. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\sigma$($Z\phi$) and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $Z\phi$ signal with pseudoscalar couplings, where $\sigma$ denotes the production cross section and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The red dash-dotted line is the theoretical prediction for $\sigma\bf{\it{B}}$ of the $Z\phi$ signal. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\sigma$($Z\phi$) and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $Z\phi$ signal with H-like production, where $\sigma$ denotes the production cross section and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The red dash-dotted line is the theoretical prediction for $\sigma\bf{\it{B}}$ of the $Z\phi$ signal. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $g^{2}_{tS}$ and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $t \bar{t} \phi$ signal with scalar couplings, where $g_{tS}$ denotes the coupling of the $\phi$ boson to the top quark and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $g^{2}_{tPS}$ and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $t \bar{t} \phi$ signal with pseudoscalar couplings, where $g_{tPS}$ denotes the coupling of the $\phi$ boson to the top quark and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $t \bar{t} \phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\Lambda^{-2}_{S}$ and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $W\phi$ signal with scalar couplings, where $\Lambda_{S}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\Lambda^{-2}_{PS}$ and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $W\phi$ signal with pseudoscalar couplings, where $\Lambda_{PS}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $W\phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\Lambda^{-2}_{S}$ and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $Z\phi$ signal with scalar couplings, where $\Lambda_{S}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $\Lambda^{-2}_{PS}$ and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $Z\phi$ signal with pseudoscalar couplings, where $\Lambda_{PS}$ denotes the mass scale of the effective interaction and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The 95% confidence level observed upper limits on the product of $sin^2 \theta$ and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ for the $Z\phi$ signal with H-like production, where $\theta$ denotes the mixing angle of the Higgs boson with the $\phi$ boson and $\bf{\it{B}}(\phi \rightarrow \ell \ell)$ is the branching fraction of the $\phi$ boson into a lepton pair of given flavor. Exclusions on the dielectron, dimuon, and ditau decay scenarios of the $\phi$ boson are shown with the green, blue, and orange solid lines, respectively. The vertical gray band indicates the mass region not considered in the analysis in the dielectron and dimuon decay scenarios of the $\phi$ boson.

The product of acceptance and efficiency, $A\varepsilon$, for a scalar $\phi$ boson in the $t\bar{t} \phi$ signal (with inclusive $t\bar{t}$ decay) in each signal region in the dielectron decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a scalar $\phi$ boson in the $t\bar{t} \phi$ signal (with inclusive $t\bar{t}$ decay) in each signal region in the dimuon decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a scalar $\phi$ boson in the $t\bar{t} \phi$ signal (with inclusive $t\bar{t}$ decay) in each signal region in the ditau decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a pseudoscalar $\phi$ boson in the $t\bar{t} \phi$ signal (with inclusive $t\bar{t}$ decay) in each signal region in the dielectron decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a pseudoscalar $\phi$ boson in the $t\bar{t} \phi$ signal (with inclusive $t\bar{t}$ decay) in each signal region in the dimuon decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a pseudoscalar $\phi$ boson in the $t\bar{t} \phi$ signal (with inclusive $t\bar{t}$ decay) in each signal region in the ditau decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a scalar $\phi$ boson in the $W\phi$ signal (with leptonic $W$ decay) in each signal region in the dielectron decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a scalar $\phi$ boson in the $W\phi$ signal (with leptonic $W$ decay) in each signal region in the dimuon decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a scalar $\phi$ boson in the $W\phi$ signal (with leptonic $W$ decay) in each signal region in the ditau decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a pseudoscalar $\phi$ boson in the $W\phi$ signal (with leptonic $W$ decay) in each signal region in the dielectron decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a pseudoscalar $\phi$ boson in the $W\phi$ signal (with leptonic $W$ decay) in each signal region in the dimuon decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a pseudoscalar $\phi$ boson in the $W\phi$ signal (with leptonic $W$ decay) in each signal region in the ditau decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for an H-like $\phi$ boson in the $W\phi$ signal (with leptonic $W$ decay) in each signal region in the dielectron decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for an H-like $\phi$ boson in the $W\phi$ signal (with leptonic $W$ decay) in each signal region in the dimuon decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for an H-like $\phi$ boson in the $W\phi$ signal (with leptonic $W$ decay) in each signal region in the ditau decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a scalar $\phi$ boson in the $Z\phi$ signal (with leptonic $Z$ decay) in each signal region in the dielectron decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a scalar $\phi$ boson in the $Z\phi$ signal (with leptonic $Z$ decay) in each signal region in the dimuon decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a scalar $\phi$ boson in the $Z\phi$ signal (with leptonic $Z$ decay) in each signal region in the ditau decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a pseudoscalar $\phi$ boson in the $Z\phi$ signal (with leptonic $Z$ decay) in each signal region in the dielectron decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a pseudoscalar $\phi$ boson in the $Z\phi$ signal (with leptonic $Z$ decay) in each signal region in the dimuon decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for a pseudoscalar $\phi$ boson in the $Z\phi$ signal (with leptonic $Z$ decay) in each signal region in the ditau decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for an H-like $\phi$ boson in the $Z\phi$ signal (with leptonic $Z$ decay) in each signal region in the dielectron decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for an H-like $\phi$ boson in the $Z\phi$ signal (with leptonic $Z$ decay) in each signal region in the dimuon decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

The product of acceptance and efficiency, $A\varepsilon$, for an H-like $\phi$ boson in the $Z\phi$ signal (with leptonic $Z$ decay) in each signal region in the ditau decay scenario. Each value is computed as the ratio of the number of simulated signal events passing all selection criteria to the total number of simulated signal events, and includes the data-to-simulation correction factors described in the paper.

Selected signal shapes of the $W\phi$(ee) signal for illustration purposes. All shape parametrizations for all coupling scenarios of the $X\phi$(ee) signal are provided in the SignalShapes_XPhiToEleEle.root file, and a README file with instructions is provided under Additional Resources.

Selected signal shapes of the $W\phi$$(\mu\mu)$ signal for illustration purposes. All shape parametrizations for all coupling scenarios of the $X\phi$$(\mu\mu)$ signal are provided in the SignalShapes_XPhiToMuMu.root file, and a README file with instructions is provided under Additional Resources.

Selected signal shapes of the $W\phi$$(\tau\tau)$ signal for illustration purposes. All shape parametrizations for all coupling scenarios of the $X\phi$$(\tau\tau)$ signal are provided in the SignalShapes_XPhiToTauTau.root file, and a README file with instructions is provided under Additional Resources.

Distributions of $\Delta\phi$ as obtained from simulation, requiring various $\textrm{t}$ tag multiplicities for the signal in two representative combinations of (gluino, neutralino) masses with large (2.2, 0.1)$\mathrm{TeV}$ and small (1.8, 1.3)$\mathrm{TeV}$ mass difference.

Results of fits to the $n_{\textrm{b}}$ multiplicity for control regions for the muon channel and with the requirements $3\leq n_{\textrm{jet}}\leq4$, $250<L_T<350~\mathrm{GeV}$, $500<H_T<750~\mathrm{GeV}$, $n_{\textrm{W}}\geq1$, $\Delta \phi<1$. The shaded area shows the fit uncertainty of the total background.

Results of fits to the $n_{\textrm{b}}$ multiplicity for control regions for the muon channel and with the requirements $3\leq n_{\textrm{jet}}\leq4$, $350<L_T<450~\mathrm{GeV}$, $H_T>1000~\mathrm{GeV}$, $n_{\textrm{W}}\geq0$, $\Delta \phi<1$. The shaded area shows the fit uncertainty of the total background.

Jet multiplicity distribution after the single-lepton baseline selection excluding the SRs for the multi-b analysis. The simulation is normalized to data with the SF mentioned in the plot.

Jet multiplicity distribution after the single-lepton baseline selection excluding the SRs for the the zero-b analysis (right). The simulation is normalized to data with the SF mentioned in the plot.

Jet multiplicity distribution in the dilepton CRs for the multi-b analysis. The simulation is normalized to data with the SF mentioned in the plot.

Jet multiplicity distribution in the dilepton CRs for the zero-b analysis. The simulation is normalized to data with the SF mentioned in the plot.

The double ratio of the single-lepton and dilepton ratio between data and simulation together with fit results and their uncertainties is shown for the multi-b analysis. The fits are performed for each data-taking year; 2018 is shown as an example.

The double ratio of the single-lepton and dilepton ratio between data and simulation together with fit results and their uncertainties is shown for the zero-b analysis. The fits are performed for each data-taking year; 2018 is shown as an example.

The prefit $L_{\mathrm{P}}$ distribution for selected electron candidates in the baseline QCD selection, with modified requirements of $n_{\text{jet}}\in[3,4]$ and $n_{\mathrm{b}}=0$.

The prefit $L_{\mathrm{P}}$ distribution for anti-selected electron candidates in the baseline QCD selection, with modified requirements of $n_{\text{jet}}\in[3,4]$ and $n_{\mathrm{b}}=0$.

Observed event yields in the MB SRs of the multi-b analysis compared to signal and background predictions. The relative fraction of the different SM EW background contributions determined in simulation is shown by the stacked, colored histograms, normalized so that their sum is equal to the background estimated using data control regions. The QCD background is predicted using the $L_p$ method. The signal is shown for two representative combinations of (gluino, neutralino) masses with large (2.2, 0.1) $\mathrm{TeV}$ and small (1.8, 1.3) $\mathrm{TeV}$ mass differences.

Observed event yields in the MB SRs of the zero-b analysis compared to signal and background predictions. The $\textrm{W}$+jets, $\textrm{t}\bar{\textrm{t}}$+jets, and QCD predictions are extracted from data control samples, while the other background contributions are estimated from simulation. The signal is shown for two representative combinations of (gluino, neutralino) masses with large (2.2, 0.1) $\mathrm{TeV}$ and small (1.8, 1.3) $\mathrm{TeV}$ mass differences.

Cross section limits at 95% CL for the T1tttt (left) model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T1tttt model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the observed (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Cross section limits at 95% CL for the T5qqqqWW model, as functions of the gluino and LSP masses, assuming a branching fraction of 100%. The mass of the intermediate chargino is taken to be halfway between the gluino and the neutralino masses. The solid black (dashed red) lines correspond to the experved (expected) mass limits, with the thicker lines representing the central values and the thinner lines representing the $\pm1\sigma$ uncertainty bands related to the theoretical (experimental) uncertainties.

Observed number of events in the MB SR bins of the multi-b analysis. All bins are defined with $\Delta\phi > 0.75$.

Observed number of events in the MB SR bins of the zero-b analysis.

The simplified likelihood can be constructed using this covariance matrix of nuisance parameters in the signal region (MB SR). The bin contents have been aggregated accross the three years to simplify the likelihood.

Selection efficiency for the multi-b signal for the SRs.

The simplified likelihood can be constructed using this covariance matrix of nuisance parameters in the signal region (MB SR). The bins were aggregated across HT to stabilize the simplified likelihood.

Selection efficiency for the zero-b signal for the SRs. The bins correspond to the aggregated bins that were used for the covariance matrix.

95% CL exclusion limits on $\mathcal{B}(H\rightarrow aa \rightarrow b\bar{b}b\bar{b})$ for $m_a = 45$ GeV.

95% CL exclusion limits on $\mathcal{B}(H\rightarrow aa \rightarrow b\bar{b}b\bar{b})$ for $m_a = 55$ GeV.

The fraction of $a$ boson decays matched to reconstructed displaced vertices passing all vertex selections in signal MC.

The extrapolated signal selection efficiency as a function of $c\tau_{a}$.

Acceptances on TRUTH level of the DMt samples in the tW1L channel for all bins in the SR. The acceptance is defined as the number of weighted TRUTH events in the SR over the number of expected events without any selections. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Efficiencies of the DMt samples in the tW2L SR. The efficiency is defined as the number of weighted reconstructed events over the number of weighted TRUTH events in the SR. The maps include all samples in the $m_a - m_H$ plane with $tan\beta = 1$.

Acceptances on TRUTH level of the DMt samples in the tW2L SR. The acceptance is defined as the number of weighted TRUTH events in the SR over the number of expected events without any selections. The maps include all samples in the $m_a - m_H$ plane with $tan\beta = 1$.

Efficiencies of the DMt samples in the tW2L SR. The efficiency is defined as the number of weighted reconstructed events over the number of weighted TRUTH events in the SR. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Acceptances on TRUTH level of the DMt samples in the tW2L SR. The acceptance is defined as the number of weighted TRUTH events in the SR over the number of expected events without any selections. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Efficiencies of the DMt samples in the tj1L channel for all bins in the SR. The efficiency is defined as the number of weighted reconstructed events over the number of weighted TRUTH events in the SR. The map includes all used samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Acceptances on TRUTH level of the DMt samples in the tj1L channel for all bins in the SR. The acceptance is defined as the number of weighted TRUTH events in the SR over the number of expected events without any selections. The map includes all used samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW1L analysis considering only DMt signal.

Upper limits on excluded cross sections of the tW1L analysis considering only the DMt signal.

The expected exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming only $tW$+DM contributions, for the tW1L analysis channel.

The observed exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming only $tW$+DM contributions, for the tW1L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW1L analysis considering only DMt signal.

Upper limits on excluded cross sections of the tW1L analysis considering only the DMt signal.

The expected exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming only $tW$+DM contributions, for the tW1L analysis channel.

The observed exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming only $tW$+DM contributions, for the tW1L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW2L analysis considering only DMt signal.

Upper limits on excluded cross sections of the tW2L analysis considering only the DMt signal.

The expected exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming only $tW$+DM contributions, for the tW2L analysis channel.

The observed exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming only $tW$+DM contributions, for the tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW2L analysis considering only DMt signal.

Upper limits on excluded cross sections of the tW2L analysis considering only the DMt signal.

The expected exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming only $tW$+DM contributions, for the tW2L analysis channel.

The observed exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming only $tW$+DM contributions, for the tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the combined tW1L and tW2L analyses considering only the DMt signal.

Upper limits on excluded cross sections of the combined tW1L and tW2L analyses considering only the DMt signal.

The expected exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming only $tW$+DM contributions, for the statistical combination of the tW1L and tW2L analysis channel.

The observed exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming only $tW$+DM contributions, for the statistical combination of the tW1L and tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the combined tW1L and tW2L analyses considering only the DMt signal.

Upper limits on excluded cross sections of the combined tW1L and tW2L analyses considering only the DMt signal.

The expected exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming only $tW$+DM contributions, for the statistical combination of the tW1L and tW2L analysis channel.

The observed exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming only $tW$+DM contributions, for the statistical combination of the tW1L and tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW1L analysis considering the DMt$\bar{t}$+DMt signal.

The expected exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW1L analysis channel.

The observed exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW1L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW1L analysis considering the DMt$\bar{t}$+DMt signal.

The expected exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW1L analysis channel.

The observed exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW1L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW2L analysis considering the DMt$\bar{t}$+DMt signal.

The expected exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW2L analysis channel.

The observed exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tW2L analysis considering the DMt$\bar{t}$+DMt signal.

The expected exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW2L analysis channel.

The observed exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the combined tW1L and tW2L analyses considering the DMt$\bar{t}$+DMt signal.

The expected exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the statistical combination of the tW1L and tW2L analysis channel.

The observed exclusion contours as a function of $(m_a, m_{H^{\pm}})$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the statistical combination of the tW1L and tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the combined tW1L and tW2L analyses considering the DMt$\bar{t}$+DMt signal.

The expected exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the statistical combination of the tW1L and tW2L analysis channel.

The observed exclusion contours as a function of $(m_{H^{\pm}}, \tan\beta)$, assuming DM$t\bar{t}$ and DM$t$ contributions, for the statistical combination of the tW1L and tW2L analysis channel.

Upper limits on signal strength (excluded cross section over theoretical cross section) of the tj1L analysis considering only the DMt signal.

Upper limits on upper limits on excluded cross sections of the tj1L analysis considering only the DMt signal.

The expected and observed cross section exclusion limits as a function of $m_{H^{\pm}}$ in the tj1L analysis channel for signal models with $m_a = 250~GeV$, and $\tan\beta=0.3$. The $\sigma^{}_\mathrm{BSM}$ is the cross section of the $t$-channel DM production process.

The expected and observed cross section exclusion limits as a function of $m_{H^{\pm}}$ in the tj1L analysis channel for signal models with $m_a = 250~GeV$, and $\tan\beta=0.5$. The $\sigma^{}_\mathrm{BSM}$ is the cross section of the $t$-channel DM production process.

Cross sections of the DMt samples in the tW1L channel. The maps include all samples in the $m_a - m_H$ plane with $tan\beta = 1$.

Cross sections of the DMt samples in the tW1L channel. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Cross sections times branching ratio of the DMt samples in the tW2L channel. The maps include all samples in the $m_a - m_H$ plane with $tan\beta = 1$.

Cross sections times branching ratio of the DMt samples in the tW2L channel. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Cross sections of the DMt samples in the tj1L channel. The map includes all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

MC generator filter efficiencies of the DMt samples in the tW1L channel. The maps include all samples in the $m_a - m_H$ plane with $tan\beta = 1$.

MC generator filter efficiencies of the DMt samples in the tW1L channel. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

MC generator filter efficiencies of the DMt samples in the tW2L channel. The maps include all samples in the $m_a - m_H$ plane with $tan\beta = 1$.

MC generator filter efficiencies of the DMt samples in the tW2L channel. The maps include all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

MC generator filter efficiencies of the DMt samples in the tj1L channel. The map includes all samples in the $m_H - tan\beta$ plane with $m_a = 250~GeV$.

Background-only fit results for the tW1L and tW2L signal regions. The backgrounds which contribute only a small amount (rare processes such as triboson, Higgs boson production processes, $t\bar{t}t\bar{t}$, $t\bar{t}WW$ and non-prompt or misidentified leptons background) are grouped and labelled as ``Others´´. The quoted uncertainties on the fitted SM background include both the statistical and systematic uncertainties.

Background-only fit results for the tj1L signal regions. The backgrounds which contribute only a small amount ($Z$+jets, rare processes such as $tWZ$, triboson, Higgs boson production processes, ,$t\bar{t}t\bar{t}$, $t\bar{t}WW$) are grouped and labelled as ``Others´´. The quoted uncertainties on the fitted SM background include both the statistical and systematic uncertainties.

Cutflow of the weighted events with statistical uncertainties for two DMt samples in all bins of the tW1L channel. The PreSelection includes at least 1 lepton in the event, at least 1 $b$-jet with $p_{\mathrm{T}} > 50~GeV$, $m\mathrm{_{T}^{lep}} > 30~GeV$, $\Delta\phi\mathrm{_{4jets, MET}^{min}} > 0.5$ and $E\mathrm{_{T}^{miss}} > 200~GeV$.

Cutflow of the weighted events with statistical uncertainties for two DMt samples in the tW2L channel. The PreSelection includes at least 2 leptons in the event, at least 1 $b$-jet with $p_{\mathrm{T}} > 40~GeV$, $m_{ll} > 40~GeV$, $m\mathrm{_{T2}} > 40~GeV$, $\Delta\phi\mathrm{_{4jets, MET}^{min}} > 0.5$ and $E\mathrm{_{T}^{miss}} > 200~GeV$.

Cutflow of the weighted events with the statistical uncertainties (except for the first cuts) for two DMt samples in all bins off the tj1L channel. The PreSelection includes at least 1 lepton in the event, at least 1 $b$-jet with $p_{\mathrm{T}} > 50~GeV$, $m\mathrm{_{T}^{lep}} > 30~GeV$, $\Delta\phi\mathrm{_{4jets, MET}^{min}} > 0.5$ and $E\mathrm{_{T}^{miss}} > 200~GeV$.