Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the LL pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the LL pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the LL pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the TT pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of corrected HH mass (mHH) mHH, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

Slices of 2D distributions of observed events and the post-fit templates in the semi-resolved pass region, projected onto the plane of leading jet mass mJ1, including expected radion signal at 1.5 TeV.

The observed (solid black line) and expected (dashed black line) upper limits at 95% CL on $\sigma$(pp->X)B(X->HH->4b) for the narrow spin-0 radion model. The green (yellow) bands represent one (two) standard deviations from the expected limit. The predicted theoretical cross sections are also shown.

The observed (solid black line) and expected (dashed black line) upper limits at 95% CL on $\sigma$(pp->X)B(X->HH->4b) for the narrow width spin-2 bulk graviton model. The green (yellow) bands represent one (two) standard deviations from the expected limit. The predicted theoretical cross sections are also shown.

The observed (solid black line) and expected (dashed black line) upper limits at 95% CL on $\sigma$(pp->X)B(X->HH->4b) for the 10%-width spin-0 radion model. The green (yellow) bands represent one (two) standard deviations from the expected limit. The predicted theoretical cross sections are also shown.

The observed (solid black line) and expected (dashed black line) upper limits at 95% CL on $\sigma$(pp->X)B(X->HH->4b) for the 10%-width spin-2 bulk graviton model. The green (yellow) bands represent one (two) standard deviations from the expected limit. The predicted theoretical cross sections are also shown.

Upper 95% CL limits on HH signal cross section scanned over the $\kappa_{\lambda}$ parameter while fixing the $\kappa_{2V}$ and $\kappa_{V}$ to their SM-predicted values.

Upper 95% CL limits on VHH signal cross section scanned over the $\kappa_{2V}$ parameter while fixing the $\kappa_{\lambda}$ and $\kappa_{V}$ to their SM-predicted values.

Upper 95% CL limits on HH signal cross section scanned over the $\kappa_{2V}$ parameter while fixing the $\kappa_{\lambda}$ and $\kappa_{V}$ to their SM-predicted values.

Upper 95% CL limits on VHH signal cross section scanned over the $\kappa_{V}$ parameter while fixing the $\kappa_{2V}$ and $\kappa_{\lambda}$ to their SM-predicted values.

Upper 95% CL limits on HH signal cross section scanned over the $\kappa_{V}$ parameter while fixing the $\kappa_{2V}$ and $\kappa_{\lambda}$ to their SM-predicted values.

Differential cross sections normalized to the fiducial cross section as a function of the $p_T$ of the second-highest-$p_T$ jet

Differential cross sections normalized to the fiducial cross section as a function of the $|\eta|$ of the highest-$p_T$ jet

Differential cross sections normalized to the fiducial cross section as a function of the $|\eta|$ of the second-highest-$p_T$ jet

Differential cross sections normalized to the fiducial cross section as a function of the pseudorapidity difference between the two highest-$p_T$ jets in the event

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the two highest-$p_T$ jets in the event

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system, in the full $m_{4l}$ range

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with 0 jet, in the on-shell ZZ region

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with 1 jet, in the on-shell ZZ region

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with 2 jets, in the on-shell ZZ region

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with at least 3 jets, in the on-shell ZZ region

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with 0 jet, in the full $m_{4l}$ range

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with 1 jet, in the full $m_{4l}$ range

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with 2 jets, in the full $m_{4l}$ range

Differential cross sections normalized to the fiducial cross section as a function of the invariant mass of the four-lepton system in events with at least 3 jets, in the full $m_{4l}$ range

Observed and expected 95% CL upper limits on the nonresonant HH production via vector boson fusion cross section obtained for both single-lepton and dilepton channels, and from their combination

Observed and expected 95% CL upper limits on the nonresonant HH production cross section for two different benchmark scenarios “JHEP04(2016)01” and “JHEP03(2020)91”

Observed and expected 95% CL upper limits on the nonresonant HH production cross section as a function of the Higgs boson self-coupling strength modifier κλ All Higgs boson couplings other than λ are assumed to have the values predicted by the SM

Observed and expected 95% CL upper limits on the nonresonant HH production via VBF cross section as a function of the effective coupling κ2V. The ggF contribution in this case is set to the SM expectation. All other Higgs boson couplings are assumed to have the values predicted by the SM.

Observed and expected 95% CL upper limits on the nonresonant HH production cross section as a function of the effective coupling c2. All other Higgs boson couplings are assumed to have the values predicted in the SM

Observed upper limits at 95% CL on B($\text{H} \rightarrow \text{a}_{1}\text{a}_{1}$) in percent, obtained from the combination of the $\mu\mu$bb and $\tau\tau$bb channels. The results are obtained as functions of $m_{\text{a}_{1}}$ for 2HDM+S Type I (independent of tan$\beta$), Type II (tan$\beta$=2.0), Type III (tan$\beta$=2.0), and Type IV (tan$\beta$=0.6), respectively.

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.

Signal yields for the charge scenarios considered in the analysis, as well as their associateds uncertainties.

Signal yields for two charge scenarios considered in the analysis, as well as their associated uncertainties.

Signal yields for the charge scenarios considered in the analysis, as well as their associateds uncertainties.

Signal yields for two charge scenarios considered in the analysis, as well as their associated uncertainties.

Signal yields for the charge scenarios considered in the analysis, as well as their associateds uncertainties.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the early 2016 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 6%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text. Comparing with respect to the binomial fit starting at $N_{\text{hits}}^{\text{low dE/dx}} = 2$, and not $N_{\text{hits}}^{\text{low dE/dx}} = 1$, is needed to account for the fact that early 2016 data is more strongly affected low dE/dx by instrumental effects that widen the N hits distribution.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the early 2016 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 6%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text. Comparing with respect to the binomial fit starting at $N_{\text{hits}}^{\text{low dE/dx}} = 2$, and not $N_{\text{hits}}^{\text{low dE/dx}} = 1$, is needed to account for the fact that early 2016 data is more strongly affected low dE/dx by instrumental effects that widen the N hits distribution.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the late 2016 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 78%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the late 2016 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 78%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the 2017 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 65%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the 2017 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 65%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the 2018 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 9%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Distribution of $N_{\text{hits}}^{\text{low dE/dx}}$ in the SR and the CR for the 2018 data set, as well as for an FCP signal at a mass of 100 GeV and different charge scenarios. The vertical bars and the shaded area correspond to the statistical uncertainty in the SR and the CR, respectively. The p-value of the fit is 9%. The two lower panels show the ratio of the number of tracks observed in the CR (upper) and SR (lower), and the fit function. The vertical bars correspond to the uncertainty from statistical sources, while the shaded area shows the systematic uncertainty in the fit due to the choice of the fitting function and the binomial fit range as explained in low dE/dx the main text.

Exclusion region (hatched) at 95% CL in the FCP charge-mass plane for the considered signal. The expected exclusion is shown with the associated 1 and 2 standard deviations $\sigma$ bands. Signal points at charges 0.9, 0.8, 2/3, 0.5, and 1/3 e are connected by straight lines to guide the eye. This is a conservative interpolation. Previous exclusions from CMS [Phys. Rev. D 87 (2013) 092008), JHEP 07 (2013) 122] as well as OPAL [Phys. Lett. B 572 (2003) 8] are given for comparison.

Exclusion region (hatched) at 95% CL in the FCP charge-mass plane for the considered signal. The expected exclusion is shown with the associated 1 and 2 standard deviations $\sigma$ bands. Signal points at charges 0.9, 0.8, 2/3, 0.5, and 1/3 e are connected by straight lines to guide the eye. This is a conservative interpolation. Previous exclusions from CMS [Phys. Rev. D 87 (2013) 092008), JHEP 07 (2013) 122] as well as OPAL [Phys. Lett. B 572 (2003) 8] are given for comparison.

Absolute differential ttbar production cross section measured as function of ttbar rapidity in bins of ttbar mass at the parton level in the full phase space.

The JEC, JER, and total systematic uncertainties in unfolded cross sections as functions of transverse momentum, for 1.5<|y|<2. The total systematic uncertainty includes also the luminosity, jet identification and trigger efficiency uncertainties.

The unfolded measured particle-level inclusive jet cross section as functions of jet pT (markers), for |y|<0.5, compared to the NLO perturbative QCD prediction (histogram), using the CT14NLO PDF set, with muR = muF = HT, and corrected for the NP effects. The experimental and theoretical systematic uncertainties are shown.

The unfolded measured particle-level inclusive jet cross section as functions of jet pT (markers), for 0.5<|y|<1, compared to the NLO perturbative QCD prediction (histogram), using the CT14NLO PDF set, with muR = muF = HT, and corrected for the NP effects. The experimental and theoretical systematic uncertainties are shown.

The unfolded measured particle-level inclusive jet cross section as functions of jet pT (markers), for 1<|y|<1.5, compared to the NLO perturbative QCD prediction (histogram), using the CT14NLO PDF set, with muR = muF = HT, and corrected for the NP effects. The experimental and theoretical systematic uncertainties are shown.

The unfolded measured particle-level inclusive jet cross section as functions of jet pT (markers), for 1.5<|y|<2, compared to the NLO perturbative QCD prediction (histogram), using the CT14NLO PDF set, with muR = muF = HT, and corrected for the NP effects. The experimental and theoretical systematic uncertainties are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = pT, for |y|<0.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = pT, for 0.5<|y|<1. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = pT, for 1<|y|<1.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = pT, for 1.5<|y|<2. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = HT, for |y|<0.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = HT, for 0.5<|y|<1. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = HT, for 1<|y|<1.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NLO theoretical predictions, using the CT14NLO PDF set, with mu = HT, for 1.5<|y|<2. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the CT14NNLO PDF set, with mu = HT, for |y|<0.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the CT14NNLO PDF set, with mu = HT, for 0.5<|y|<1. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the CT14NNLO PDF set, with mu = HT, for 1<|y|<1.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the CT14NNLO PDF set, with mu = HT, for 1.5<|y|<2. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the NNPDF31NNLO PDF set, with mu = HT, for |y|<0.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the NNPDF31NNLO PDF set, with mu = HT, for 0.5<|y|<1. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the NNPDF31NNLO PDF set, with mu = HT, for 1<|y|<1.5. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

Ratios (points) of the unfolded measured cross sections to the NNLO theoretical predictions, using the NNPDF31NNLO PDF set, with mu = HT, for 1.5<|y|<2. The vertical error bars show the statistical experimental uncertainty. The systematic experimental uncertainty, the total theoretical uncertainty, and the individual sources of theoretical uncertainty are shown.

The effect of aS(MZ) variation, for |y|<0.5. The NNLO theoretical cross section predictions using the NNPDF31NNLO PDF with m = HT, calculated for different choices of aS (0.108, 0.110, 0.112, 0.114, 0.116, 0.117, 0.118, 0.119, 0.120, 0.122, and 0.124), are divided by the benchmark NNLO prediction for aS = 0.118 and the same choice of PDF set, muR, and muF. Also shown is the experimental unfolded measurement divided by the same benchmark prediction. The width of the unity line corresponds to the statistical uncertainty from the MC integration for the determination of the NNLO prediction. The error bars on the unfolded data correspond to the total experimental statistical and systematic uncertainty added in quadrature.

The effect of aS(MZ) variation, for 0.5<|y|<1. The NNLO theoretical cross section predictions using the NNPDF31NNLO PDF with m = HT, calculated for different choices of aS (0.108, 0.110, 0.112, 0.114, 0.116, 0.117, 0.118, 0.119, 0.120, 0.122, and 0.124), are divided by the benchmark NNLO prediction for aS = 0.118 and the same choice of PDF set, muR, and muF. Also shown is the experimental unfolded measurement divided by the same benchmark prediction. The width of the unity line corresponds to the statistical uncertainty from the MC integration for the determination of the NNLO prediction. The error bars on the unfolded data correspond to the total experimental statistical and systematic uncertainty added in quadrature.

The effect of aS(MZ) variation, for 1<|y|<1.5. The NNLO theoretical cross section predictions using the NNPDF31NNLO PDF with m = HT, calculated for different choices of aS (0.108, 0.110, 0.112, 0.114, 0.116, 0.117, 0.118, 0.119, 0.120, 0.122, and 0.124), are divided by the benchmark NNLO prediction for aS = 0.118 and the same choice of PDF set, muR, and muF. Also shown is the experimental unfolded measurement divided by the same benchmark prediction. The width of the unity line corresponds to the statistical uncertainty from the MC integration for the determination of the NNLO prediction. The error bars on the unfolded data correspond to the total experimental statistical and systematic uncertainty added in quadrature.

The effect of aS(MZ) variation, for 1.5<|y|<2. The NNLO theoretical cross section predictions using the NNPDF31NNLO PDF with m = HT, calculated for different choices of aS (0.108, 0.110, 0.112, 0.114, 0.116, 0.117, 0.118, 0.119, 0.120, 0.122, and 0.124), are divided by the benchmark NNLO prediction for aS = 0.118 and the same choice of PDF set, muR, and muF. Also shown is the experimental unfolded measurement divided by the same benchmark prediction. The width of the unity line corresponds to the statistical uncertainty from the MC integration for the determination of the NNLO prediction. The error bars on the unfolded data correspond to the total experimental statistical and systematic uncertainty added in quadrature.

Unfolded correlation matrix for |y|<0.5.

Unfolded correlation matrix for 0.5<|y|<1.

Unfolded correlation matrix for 1<|y|<1.5.

Unfolded correlation matrix for 1.5<|y|<2.