CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-B2G-24-016 ; CERN-EP-2026-248
Search for pair-produced vector-like top quarks decaying into Lorentz-boosted top quarks and scalar bosons in proton-proton collisions at $ \sqrt{s}= $ 13 TeV
Submitted to Physical Review D
Abstract: A search for pair-produced vector-like top quarks, $ T $, decaying into a new scalar boson $ \phi $ and a standard model (SM) top quark is presented. The search targets events containing a photon pair originating from the decay of one of the scalar bosons, accompanied by all-hadronic decays of the SM top quark pair. The analysis is optimized for boosted top quark signatures in which the decay products of the top quarks are reconstructed as large-radius jets. The search is based on proton-proton collision data at $ \sqrt{s}= $ 13 TeV, collected by the CMS experiment at the CERN LHC during 2016--2018, corresponding to an integrated luminosity of 138 fb$ ^{-1} $. The search is performed for an array of T and $ \phi $ masses. No significant deviation from the predicted background is observed. The largest local excess has a significance of 2.7 standard deviations, reduced to 1.5 standard deviations after accounting for the look-elsewhere effect. Upper limits at 95% confidence level are set on the product of $ T \bar{T}$ production cross section and the branching fraction.
Figures & Tables Summary References CMS Publications
Figures

png pdf
Figure 1:
A representative leading order Feynman diagram for the pair production of vector-like top quarks $ T $ with subsequent $ T \to \phi\mathrm{t} $ decays.

png pdf
Figure 2:
Distribution of reconstructed T candidates in the $ m^\text{reco}_{T } $--$ m^\text{reco}_{\phi} $ plane for data (upper) and simulated $ \mathrm{t}\overline{\mathrm{t}} $ events (lower), shown for the SR (left) and CR (right). The simulated $ \mathrm{t}\overline{\mathrm{t}} $ events are normalized to sum to the total data events in the SR and the CR, respectively.

png pdf
Figure 2-a:
Distribution of reconstructed T candidates in the $ m^\text{reco}_{T } $--$ m^\text{reco}_{\phi} $ plane for data (upper) and simulated $ \mathrm{t}\overline{\mathrm{t}} $ events (lower), shown for the SR (left) and CR (right). The simulated $ \mathrm{t}\overline{\mathrm{t}} $ events are normalized to sum to the total data events in the SR and the CR, respectively.

png pdf
Figure 2-b:
Distribution of reconstructed T candidates in the $ m^\text{reco}_{T } $--$ m^\text{reco}_{\phi} $ plane for data (upper) and simulated $ \mathrm{t}\overline{\mathrm{t}} $ events (lower), shown for the SR (left) and CR (right). The simulated $ \mathrm{t}\overline{\mathrm{t}} $ events are normalized to sum to the total data events in the SR and the CR, respectively.

png pdf
Figure 2-c:
Distribution of reconstructed T candidates in the $ m^\text{reco}_{T } $--$ m^\text{reco}_{\phi} $ plane for data (upper) and simulated $ \mathrm{t}\overline{\mathrm{t}} $ events (lower), shown for the SR (left) and CR (right). The simulated $ \mathrm{t}\overline{\mathrm{t}} $ events are normalized to sum to the total data events in the SR and the CR, respectively.

png pdf
Figure 2-d:
Distribution of reconstructed T candidates in the $ m^\text{reco}_{T } $--$ m^\text{reco}_{\phi} $ plane for data (upper) and simulated $ \mathrm{t}\overline{\mathrm{t}} $ events (lower), shown for the SR (left) and CR (right). The simulated $ \mathrm{t}\overline{\mathrm{t}} $ events are normalized to sum to the total data events in the SR and the CR, respectively.

png pdf
Figure 3:
Comparison of data and postfit background predictions from the background-only fit in the VR. The VR-CR (left) and VR-SR (right) regions are shown. The distributions are projected onto $ m^\text{reco}_{T } $ in three intervals of $ m^\text{reco}_{\phi} $: 25 $ < m^\text{reco}_{\phi} < $ 75 GeV (upper), 75 $ < m^\text{reco}_{\phi} < $ 325 GeV (middle), and 325 $ < m^\text{reco}_{\phi} < $ 825 GeV (lower). The $ \mathrm{t}\overline{\mathrm{t}} $ component scaled by the transfer function is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 3-a:
Comparison of data and postfit background predictions from the background-only fit in the VR. The VR-CR (left) and VR-SR (right) regions are shown. The distributions are projected onto $ m^\text{reco}_{T } $ in three intervals of $ m^\text{reco}_{\phi} $: 25 $ < m^\text{reco}_{\phi} < $ 75 GeV (upper), 75 $ < m^\text{reco}_{\phi} < $ 325 GeV (middle), and 325 $ < m^\text{reco}_{\phi} < $ 825 GeV (lower). The $ \mathrm{t}\overline{\mathrm{t}} $ component scaled by the transfer function is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 3-b:
Comparison of data and postfit background predictions from the background-only fit in the VR. The VR-CR (left) and VR-SR (right) regions are shown. The distributions are projected onto $ m^\text{reco}_{T } $ in three intervals of $ m^\text{reco}_{\phi} $: 25 $ < m^\text{reco}_{\phi} < $ 75 GeV (upper), 75 $ < m^\text{reco}_{\phi} < $ 325 GeV (middle), and 325 $ < m^\text{reco}_{\phi} < $ 825 GeV (lower). The $ \mathrm{t}\overline{\mathrm{t}} $ component scaled by the transfer function is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 3-c:
Comparison of data and postfit background predictions from the background-only fit in the VR. The VR-CR (left) and VR-SR (right) regions are shown. The distributions are projected onto $ m^\text{reco}_{T } $ in three intervals of $ m^\text{reco}_{\phi} $: 25 $ < m^\text{reco}_{\phi} < $ 75 GeV (upper), 75 $ < m^\text{reco}_{\phi} < $ 325 GeV (middle), and 325 $ < m^\text{reco}_{\phi} < $ 825 GeV (lower). The $ \mathrm{t}\overline{\mathrm{t}} $ component scaled by the transfer function is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 3-d:
Comparison of data and postfit background predictions from the background-only fit in the VR. The VR-CR (left) and VR-SR (right) regions are shown. The distributions are projected onto $ m^\text{reco}_{T } $ in three intervals of $ m^\text{reco}_{\phi} $: 25 $ < m^\text{reco}_{\phi} < $ 75 GeV (upper), 75 $ < m^\text{reco}_{\phi} < $ 325 GeV (middle), and 325 $ < m^\text{reco}_{\phi} < $ 825 GeV (lower). The $ \mathrm{t}\overline{\mathrm{t}} $ component scaled by the transfer function is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 3-e:
Comparison of data and postfit background predictions from the background-only fit in the VR. The VR-CR (left) and VR-SR (right) regions are shown. The distributions are projected onto $ m^\text{reco}_{T } $ in three intervals of $ m^\text{reco}_{\phi} $: 25 $ < m^\text{reco}_{\phi} < $ 75 GeV (upper), 75 $ < m^\text{reco}_{\phi} < $ 325 GeV (middle), and 325 $ < m^\text{reco}_{\phi} < $ 825 GeV (lower). The $ \mathrm{t}\overline{\mathrm{t}} $ component scaled by the transfer function is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 3-f:
Comparison of data and postfit background predictions from the background-only fit in the VR. The VR-CR (left) and VR-SR (right) regions are shown. The distributions are projected onto $ m^\text{reco}_{T } $ in three intervals of $ m^\text{reco}_{\phi} $: 25 $ < m^\text{reco}_{\phi} < $ 75 GeV (upper), 75 $ < m^\text{reco}_{\phi} < $ 325 GeV (middle), and 325 $ < m^\text{reco}_{\phi} < $ 825 GeV (lower). The $ \mathrm{t}\overline{\mathrm{t}} $ component scaled by the transfer function is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 4:
Comparison of data and postfit background predictions from the background-only fit in the VR. The VR-CR (left) and VR-SR (right) regions are shown. The distributions are projected onto $ m^\text{reco}_{\phi} $ in three intervals of $ m^\text{reco}_{T } $: 625 $ < m^\text{reco}_{T } < $ 825 GeV (upper), 825 $ < m^\text{reco}_{T } < $ 1125 GeV (middle), and 1125 $ < m^\text{reco}_{T } < $ 1525 GeV (lower). The $ \mathrm{t}\overline{\mathrm{t}} $ component scaled by the transfer function is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 4-a:
Comparison of data and postfit background predictions from the background-only fit in the VR. The VR-CR (left) and VR-SR (right) regions are shown. The distributions are projected onto $ m^\text{reco}_{\phi} $ in three intervals of $ m^\text{reco}_{T } $: 625 $ < m^\text{reco}_{T } < $ 825 GeV (upper), 825 $ < m^\text{reco}_{T } < $ 1125 GeV (middle), and 1125 $ < m^\text{reco}_{T } < $ 1525 GeV (lower). The $ \mathrm{t}\overline{\mathrm{t}} $ component scaled by the transfer function is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 4-b:
Comparison of data and postfit background predictions from the background-only fit in the VR. The VR-CR (left) and VR-SR (right) regions are shown. The distributions are projected onto $ m^\text{reco}_{\phi} $ in three intervals of $ m^\text{reco}_{T } $: 625 $ < m^\text{reco}_{T } < $ 825 GeV (upper), 825 $ < m^\text{reco}_{T } < $ 1125 GeV (middle), and 1125 $ < m^\text{reco}_{T } < $ 1525 GeV (lower). The $ \mathrm{t}\overline{\mathrm{t}} $ component scaled by the transfer function is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 4-c:
Comparison of data and postfit background predictions from the background-only fit in the VR. The VR-CR (left) and VR-SR (right) regions are shown. The distributions are projected onto $ m^\text{reco}_{\phi} $ in three intervals of $ m^\text{reco}_{T } $: 625 $ < m^\text{reco}_{T } < $ 825 GeV (upper), 825 $ < m^\text{reco}_{T } < $ 1125 GeV (middle), and 1125 $ < m^\text{reco}_{T } < $ 1525 GeV (lower). The $ \mathrm{t}\overline{\mathrm{t}} $ component scaled by the transfer function is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 4-d:
Comparison of data and postfit background predictions from the background-only fit in the VR. The VR-CR (left) and VR-SR (right) regions are shown. The distributions are projected onto $ m^\text{reco}_{\phi} $ in three intervals of $ m^\text{reco}_{T } $: 625 $ < m^\text{reco}_{T } < $ 825 GeV (upper), 825 $ < m^\text{reco}_{T } < $ 1125 GeV (middle), and 1125 $ < m^\text{reco}_{T } < $ 1525 GeV (lower). The $ \mathrm{t}\overline{\mathrm{t}} $ component scaled by the transfer function is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 4-e:
Comparison of data and postfit background predictions from the background-only fit in the VR. The VR-CR (left) and VR-SR (right) regions are shown. The distributions are projected onto $ m^\text{reco}_{\phi} $ in three intervals of $ m^\text{reco}_{T } $: 625 $ < m^\text{reco}_{T } < $ 825 GeV (upper), 825 $ < m^\text{reco}_{T } < $ 1125 GeV (middle), and 1125 $ < m^\text{reco}_{T } < $ 1525 GeV (lower). The $ \mathrm{t}\overline{\mathrm{t}} $ component scaled by the transfer function is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 4-f:
Comparison of data and postfit background predictions from the background-only fit in the VR. The VR-CR (left) and VR-SR (right) regions are shown. The distributions are projected onto $ m^\text{reco}_{\phi} $ in three intervals of $ m^\text{reco}_{T } $: 625 $ < m^\text{reco}_{T } < $ 825 GeV (upper), 825 $ < m^\text{reco}_{T } < $ 1125 GeV (middle), and 1125 $ < m^\text{reco}_{T } < $ 1525 GeV (lower). The $ \mathrm{t}\overline{\mathrm{t}} $ component scaled by the transfer function is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 5:
Comparison of the data and postfit background predictions projected onto $ m^\text{reco}_{T } $ for the signal hypothesis $ m_{T }= $ 900 GeV and $ m_{\phi}= $ 150 GeV. The CR (left) and SR (right) regions are shown. The distributions are shown in three intervals of $ m^\text{reco}_{\phi} $: 25 $ < m^\text{reco}_{\phi} < $ 75 GeV (upper), 75 $ < m^\text{reco}_{\phi} < $ 325 GeV (middle), and 325 $ < m^\text{reco}_{\phi} < $ 825 GeV (lower). The transferred $ \mathrm{t}\overline{\mathrm{t}} $ contribution is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The signal distribution is overlaid as a blue line. The shaded band represents the uncertainty in the total background predictions at 95% confidence level. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 5-a:
Comparison of the data and postfit background predictions projected onto $ m^\text{reco}_{T } $ for the signal hypothesis $ m_{T }= $ 900 GeV and $ m_{\phi}= $ 150 GeV. The CR (left) and SR (right) regions are shown. The distributions are shown in three intervals of $ m^\text{reco}_{\phi} $: 25 $ < m^\text{reco}_{\phi} < $ 75 GeV (upper), 75 $ < m^\text{reco}_{\phi} < $ 325 GeV (middle), and 325 $ < m^\text{reco}_{\phi} < $ 825 GeV (lower). The transferred $ \mathrm{t}\overline{\mathrm{t}} $ contribution is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The signal distribution is overlaid as a blue line. The shaded band represents the uncertainty in the total background predictions at 95% confidence level. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 5-b:
Comparison of the data and postfit background predictions projected onto $ m^\text{reco}_{T } $ for the signal hypothesis $ m_{T }= $ 900 GeV and $ m_{\phi}= $ 150 GeV. The CR (left) and SR (right) regions are shown. The distributions are shown in three intervals of $ m^\text{reco}_{\phi} $: 25 $ < m^\text{reco}_{\phi} < $ 75 GeV (upper), 75 $ < m^\text{reco}_{\phi} < $ 325 GeV (middle), and 325 $ < m^\text{reco}_{\phi} < $ 825 GeV (lower). The transferred $ \mathrm{t}\overline{\mathrm{t}} $ contribution is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The signal distribution is overlaid as a blue line. The shaded band represents the uncertainty in the total background predictions at 95% confidence level. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 5-c:
Comparison of the data and postfit background predictions projected onto $ m^\text{reco}_{T } $ for the signal hypothesis $ m_{T }= $ 900 GeV and $ m_{\phi}= $ 150 GeV. The CR (left) and SR (right) regions are shown. The distributions are shown in three intervals of $ m^\text{reco}_{\phi} $: 25 $ < m^\text{reco}_{\phi} < $ 75 GeV (upper), 75 $ < m^\text{reco}_{\phi} < $ 325 GeV (middle), and 325 $ < m^\text{reco}_{\phi} < $ 825 GeV (lower). The transferred $ \mathrm{t}\overline{\mathrm{t}} $ contribution is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The signal distribution is overlaid as a blue line. The shaded band represents the uncertainty in the total background predictions at 95% confidence level. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 5-d:
Comparison of the data and postfit background predictions projected onto $ m^\text{reco}_{T } $ for the signal hypothesis $ m_{T }= $ 900 GeV and $ m_{\phi}= $ 150 GeV. The CR (left) and SR (right) regions are shown. The distributions are shown in three intervals of $ m^\text{reco}_{\phi} $: 25 $ < m^\text{reco}_{\phi} < $ 75 GeV (upper), 75 $ < m^\text{reco}_{\phi} < $ 325 GeV (middle), and 325 $ < m^\text{reco}_{\phi} < $ 825 GeV (lower). The transferred $ \mathrm{t}\overline{\mathrm{t}} $ contribution is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The signal distribution is overlaid as a blue line. The shaded band represents the uncertainty in the total background predictions at 95% confidence level. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 5-e:
Comparison of the data and postfit background predictions projected onto $ m^\text{reco}_{T } $ for the signal hypothesis $ m_{T }= $ 900 GeV and $ m_{\phi}= $ 150 GeV. The CR (left) and SR (right) regions are shown. The distributions are shown in three intervals of $ m^\text{reco}_{\phi} $: 25 $ < m^\text{reco}_{\phi} < $ 75 GeV (upper), 75 $ < m^\text{reco}_{\phi} < $ 325 GeV (middle), and 325 $ < m^\text{reco}_{\phi} < $ 825 GeV (lower). The transferred $ \mathrm{t}\overline{\mathrm{t}} $ contribution is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The signal distribution is overlaid as a blue line. The shaded band represents the uncertainty in the total background predictions at 95% confidence level. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 5-f:
Comparison of the data and postfit background predictions projected onto $ m^\text{reco}_{T } $ for the signal hypothesis $ m_{T }= $ 900 GeV and $ m_{\phi}= $ 150 GeV. The CR (left) and SR (right) regions are shown. The distributions are shown in three intervals of $ m^\text{reco}_{\phi} $: 25 $ < m^\text{reco}_{\phi} < $ 75 GeV (upper), 75 $ < m^\text{reco}_{\phi} < $ 325 GeV (middle), and 325 $ < m^\text{reco}_{\phi} < $ 825 GeV (lower). The transferred $ \mathrm{t}\overline{\mathrm{t}} $ contribution is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The signal distribution is overlaid as a blue line. The shaded band represents the uncertainty in the total background predictions at 95% confidence level. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 6:
Comparison of the data and postfit background predictions projected onto $ m^\text{reco}_{\phi} $ for the signal hypothesis $ m_{T }= $ 900 GeV and $ m_{\phi}= $ 150 GeV. The CR (left) and SR (right) regions are shown. The distributions are shown in three intervals of $ m^\text{reco}_{T } $: 625 $ < m^\text{reco}_{T } < $ 825 GeV (upper), 825 $ < m^\text{reco}_{T } < $ 1125 GeV (middle), and 1125 $ < m^\text{reco}_{T } < $ 1525 GeV (lower). The transferred $ \mathrm{t}\overline{\mathrm{t}} $ contribution is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The signal distribution is overlaid as a blue line. The shaded band represents the uncertainty in the total background predictions at 95% confidence level. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 6-a:
Comparison of the data and postfit background predictions projected onto $ m^\text{reco}_{\phi} $ for the signal hypothesis $ m_{T }= $ 900 GeV and $ m_{\phi}= $ 150 GeV. The CR (left) and SR (right) regions are shown. The distributions are shown in three intervals of $ m^\text{reco}_{T } $: 625 $ < m^\text{reco}_{T } < $ 825 GeV (upper), 825 $ < m^\text{reco}_{T } < $ 1125 GeV (middle), and 1125 $ < m^\text{reco}_{T } < $ 1525 GeV (lower). The transferred $ \mathrm{t}\overline{\mathrm{t}} $ contribution is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The signal distribution is overlaid as a blue line. The shaded band represents the uncertainty in the total background predictions at 95% confidence level. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 6-b:
Comparison of the data and postfit background predictions projected onto $ m^\text{reco}_{\phi} $ for the signal hypothesis $ m_{T }= $ 900 GeV and $ m_{\phi}= $ 150 GeV. The CR (left) and SR (right) regions are shown. The distributions are shown in three intervals of $ m^\text{reco}_{T } $: 625 $ < m^\text{reco}_{T } < $ 825 GeV (upper), 825 $ < m^\text{reco}_{T } < $ 1125 GeV (middle), and 1125 $ < m^\text{reco}_{T } < $ 1525 GeV (lower). The transferred $ \mathrm{t}\overline{\mathrm{t}} $ contribution is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The signal distribution is overlaid as a blue line. The shaded band represents the uncertainty in the total background predictions at 95% confidence level. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 6-c:
Comparison of the data and postfit background predictions projected onto $ m^\text{reco}_{\phi} $ for the signal hypothesis $ m_{T }= $ 900 GeV and $ m_{\phi}= $ 150 GeV. The CR (left) and SR (right) regions are shown. The distributions are shown in three intervals of $ m^\text{reco}_{T } $: 625 $ < m^\text{reco}_{T } < $ 825 GeV (upper), 825 $ < m^\text{reco}_{T } < $ 1125 GeV (middle), and 1125 $ < m^\text{reco}_{T } < $ 1525 GeV (lower). The transferred $ \mathrm{t}\overline{\mathrm{t}} $ contribution is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The signal distribution is overlaid as a blue line. The shaded band represents the uncertainty in the total background predictions at 95% confidence level. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 6-d:
Comparison of the data and postfit background predictions projected onto $ m^\text{reco}_{\phi} $ for the signal hypothesis $ m_{T }= $ 900 GeV and $ m_{\phi}= $ 150 GeV. The CR (left) and SR (right) regions are shown. The distributions are shown in three intervals of $ m^\text{reco}_{T } $: 625 $ < m^\text{reco}_{T } < $ 825 GeV (upper), 825 $ < m^\text{reco}_{T } < $ 1125 GeV (middle), and 1125 $ < m^\text{reco}_{T } < $ 1525 GeV (lower). The transferred $ \mathrm{t}\overline{\mathrm{t}} $ contribution is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The signal distribution is overlaid as a blue line. The shaded band represents the uncertainty in the total background predictions at 95% confidence level. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 6-e:
Comparison of the data and postfit background predictions projected onto $ m^\text{reco}_{\phi} $ for the signal hypothesis $ m_{T }= $ 900 GeV and $ m_{\phi}= $ 150 GeV. The CR (left) and SR (right) regions are shown. The distributions are shown in three intervals of $ m^\text{reco}_{T } $: 625 $ < m^\text{reco}_{T } < $ 825 GeV (upper), 825 $ < m^\text{reco}_{T } < $ 1125 GeV (middle), and 1125 $ < m^\text{reco}_{T } < $ 1525 GeV (lower). The transferred $ \mathrm{t}\overline{\mathrm{t}} $ contribution is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The signal distribution is overlaid as a blue line. The shaded band represents the uncertainty in the total background predictions at 95% confidence level. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 6-f:
Comparison of the data and postfit background predictions projected onto $ m^\text{reco}_{\phi} $ for the signal hypothesis $ m_{T }= $ 900 GeV and $ m_{\phi}= $ 150 GeV. The CR (left) and SR (right) regions are shown. The distributions are shown in three intervals of $ m^\text{reco}_{T } $: 625 $ < m^\text{reco}_{T } < $ 825 GeV (upper), 825 $ < m^\text{reco}_{T } < $ 1125 GeV (middle), and 1125 $ < m^\text{reco}_{T } < $ 1525 GeV (lower). The transferred $ \mathrm{t}\overline{\mathrm{t}} $ contribution is shown in red and the $ \mathrm{W}/\mathrm{Z} $+jets contribution in green. The signal distribution is overlaid as a blue line. The shaded band represents the uncertainty in the total background predictions at 95% confidence level. The lower panels show the difference between the data and the background normalized by the uncertainty.

png pdf
Figure 7:
Expected (upper) and observed (lower) upper limits on $ \sigma_{T T } \ \mathcal{B}(T \to \phi \mathrm{t})^2 \ \mathcal{B}_{\gamma\gamma} (2 - \mathcal{B}_{\gamma\gamma}) $ in the $ (m_{T }, m_{\phi}) $ plane where $ \mathcal{B}_{\gamma\gamma} \equiv \mathcal{B}(\phi \to \gamma\gamma) $ is the branching fraction for $ \phi $ decay to two photons. The uncertainty for the expected upper limits is shown in parentheses.

png pdf
Figure 7-a:
Expected (upper) and observed (lower) upper limits on $ \sigma_{T T } \ \mathcal{B}(T \to \phi \mathrm{t})^2 \ \mathcal{B}_{\gamma\gamma} (2 - \mathcal{B}_{\gamma\gamma}) $ in the $ (m_{T }, m_{\phi}) $ plane where $ \mathcal{B}_{\gamma\gamma} \equiv \mathcal{B}(\phi \to \gamma\gamma) $ is the branching fraction for $ \phi $ decay to two photons. The uncertainty for the expected upper limits is shown in parentheses.

png pdf
Figure 7-b:
Expected (upper) and observed (lower) upper limits on $ \sigma_{T T } \ \mathcal{B}(T \to \phi \mathrm{t})^2 \ \mathcal{B}_{\gamma\gamma} (2 - \mathcal{B}_{\gamma\gamma}) $ in the $ (m_{T }, m_{\phi}) $ plane where $ \mathcal{B}_{\gamma\gamma} \equiv \mathcal{B}(\phi \to \gamma\gamma) $ is the branching fraction for $ \phi $ decay to two photons. The uncertainty for the expected upper limits is shown in parentheses.

png pdf
Figure 8:
Expected (dashed) and observed (solid) upper limits on $ \sigma_{T T } \ \mathcal{B}(T \to \phi \mathrm{t})^2 \ \mathcal{B}_{\gamma\gamma} (2 - \mathcal{B}_{\gamma\gamma}) $ as a function of $ m_{T } $ for different values of $ m_{\phi} $. The first row shows $ m_{\phi}= $ 100 and 150 GeV (left to right), the second row shows $ m_{\phi}= $ 200 and 250 GeV (left to right), and the third row shows $ m_{\phi}= $ 300 GeV. The product of the theoretical cross section and the $ \phi $ branching fraction is shown for $ \mathcal{B}_{\gamma\gamma} = $ 0.004 (magenta).

png pdf
Figure 8-a:
Expected (dashed) and observed (solid) upper limits on $ \sigma_{T T } \ \mathcal{B}(T \to \phi \mathrm{t})^2 \ \mathcal{B}_{\gamma\gamma} (2 - \mathcal{B}_{\gamma\gamma}) $ as a function of $ m_{T } $ for different values of $ m_{\phi} $. The first row shows $ m_{\phi}= $ 100 and 150 GeV (left to right), the second row shows $ m_{\phi}= $ 200 and 250 GeV (left to right), and the third row shows $ m_{\phi}= $ 300 GeV. The product of the theoretical cross section and the $ \phi $ branching fraction is shown for $ \mathcal{B}_{\gamma\gamma} = $ 0.004 (magenta).

png pdf
Figure 8-b:
Expected (dashed) and observed (solid) upper limits on $ \sigma_{T T } \ \mathcal{B}(T \to \phi \mathrm{t})^2 \ \mathcal{B}_{\gamma\gamma} (2 - \mathcal{B}_{\gamma\gamma}) $ as a function of $ m_{T } $ for different values of $ m_{\phi} $. The first row shows $ m_{\phi}= $ 100 and 150 GeV (left to right), the second row shows $ m_{\phi}= $ 200 and 250 GeV (left to right), and the third row shows $ m_{\phi}= $ 300 GeV. The product of the theoretical cross section and the $ \phi $ branching fraction is shown for $ \mathcal{B}_{\gamma\gamma} = $ 0.004 (magenta).

png pdf
Figure 8-c:
Expected (dashed) and observed (solid) upper limits on $ \sigma_{T T } \ \mathcal{B}(T \to \phi \mathrm{t})^2 \ \mathcal{B}_{\gamma\gamma} (2 - \mathcal{B}_{\gamma\gamma}) $ as a function of $ m_{T } $ for different values of $ m_{\phi} $. The first row shows $ m_{\phi}= $ 100 and 150 GeV (left to right), the second row shows $ m_{\phi}= $ 200 and 250 GeV (left to right), and the third row shows $ m_{\phi}= $ 300 GeV. The product of the theoretical cross section and the $ \phi $ branching fraction is shown for $ \mathcal{B}_{\gamma\gamma} = $ 0.004 (magenta).

png pdf
Figure 8-d:
Expected (dashed) and observed (solid) upper limits on $ \sigma_{T T } \ \mathcal{B}(T \to \phi \mathrm{t})^2 \ \mathcal{B}_{\gamma\gamma} (2 - \mathcal{B}_{\gamma\gamma}) $ as a function of $ m_{T } $ for different values of $ m_{\phi} $. The first row shows $ m_{\phi}= $ 100 and 150 GeV (left to right), the second row shows $ m_{\phi}= $ 200 and 250 GeV (left to right), and the third row shows $ m_{\phi}= $ 300 GeV. The product of the theoretical cross section and the $ \phi $ branching fraction is shown for $ \mathcal{B}_{\gamma\gamma} = $ 0.004 (magenta).

png pdf
Figure 8-e:
Expected (dashed) and observed (solid) upper limits on $ \sigma_{T T } \ \mathcal{B}(T \to \phi \mathrm{t})^2 \ \mathcal{B}_{\gamma\gamma} (2 - \mathcal{B}_{\gamma\gamma}) $ as a function of $ m_{T } $ for different values of $ m_{\phi} $. The first row shows $ m_{\phi}= $ 100 and 150 GeV (left to right), the second row shows $ m_{\phi}= $ 200 and 250 GeV (left to right), and the third row shows $ m_{\phi}= $ 300 GeV. The product of the theoretical cross section and the $ \phi $ branching fraction is shown for $ \mathcal{B}_{\gamma\gamma} = $ 0.004 (magenta).
Tables

png pdf
Table 1:
The baseline, signal region, and control region event selection.

png pdf
Table 2:
Largest postfit relative impacts of the dominant systematic uncertainties on the fitted signal-strength parameter $ r $ for mass point ($ m_{T }= $ 900 GeV,$ m_{\phi}= $ 150 GeV).
Summary
A search for the pair production of vector-like top quarks $ T $ in a scenario involving an additional scalar boson has been presented. The search uses proton-proton collision data collected by the CMS experiment at the LHC during 2016--2018, at center-of-mass energy 13 TeV, corresponding to an integrated luminosity of 138 fb$ ^{-1} $. The analysis considers $ T \bar{T}$ production, with each $ T $ decaying into a standard model top quark and a scalar boson $ \phi $. The top quarks are required to decay hadronically, and one of the scalar bosons is required to decay to a pair of photons. The top quarks are reconstructed as large-radius jets and identified using a ParticleNet discriminant, while a multivariate discriminant is used for photon identification. The background in the signal region is predicted from the distribution of data in a $ \mathrm{t}\overline{\mathrm{t}} $ control region via a transfer function. No significant deviation from the background prediction is observed. The largest local excess has a significance of 2.7 standard deviations, reduced to 1.5 standard deviations after accounting for the look-elsewhere effect. This is the first exclusive search for pair-produced $ T $ quarks in the boosted regime.
References
1 ATLAS Collaboration Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS detector at the LHC PLB 716 (2012) 1 1207.7214
2 CMS Collaboration Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC PLB 716 (2012) 30 CMS-HIG-12-028
1207.7235
3 CMS Collaboration A portrait of the Higgs boson by the CMS experiment ten years after the discovery Nature 607 (2022) 60 CMS-HIG-22-001
2207.00043
4 R. Barbieri and G. F. Giudice Upper bounds on supersymmetric particle masses NPB 306 (1988) 63
5 M. E. Peskin What is the hierarchy problem? NPB 1018 (2025) 116971
6 B. Guo, Y.-X. Liu, K. Yang, and S.-W. Wei Hierarchy problem and new warped extra dimension PRD 98 (2018) 085022 1805.03976
7 G. F. Giudice, C. Grojean, A. Pomarol, and R. Rattazzi The strongly-interacting light Higgs JHEP 06 (2007) 045 hep-ph/0703164
8 P. Lodone Vector-like quarks in a 'composite' Higgs model JHEP 12 (2008) 029 0806.1472
9 R. Benbrik et al. Signatures of vector-like top partners decaying into new neutral scalar or pseudoscalar bosons JHEP 05 (2020) 028 1907.05929
10 CMS Collaboration Search for pair production of vector-like quarks in leptonic final states in proton-proton collisions at $ \sqrt{s} = $ 13 TeV JHEP 07 (2023) 020 2209.07327
11 R. Benbrik et al. Vector-like quarks at the LHC: A unified perspective from ATLAS and CMS exclusion limits JHEP 03 (2025) 020 2412.01761
12 ATLAS Collaboration Search for pair-production of vector-like $ T $ quarks decaying into a top quark and a spin-0 particle in the diphoton final state in proton-proton collisions at $ \sqrt{s}= $ 13 TeV with the ATLAS detector 2607.28381
13 D. Stremmer and M. Worek Associated production of a top-quark pair with two isolated photons at the LHC through NLO in QCD 2306.16968
14 ATLAS Collaboration Exploration at the high-energy frontier: ATLAS Run 2 searches investigating the exotic jungle beyond the Standard Model Phys. Rept. 1116 (2025) 301 2403.09292
15 CMS Collaboration Review of searches for vector-like quarks, vector-like leptons, and heavy neutral leptons in proton-proton collisions at \ensuremath\sqrts = 13 TeV at the CMS experiment Phys. Rept. 1115 (2025) 570 CMS-EXO-23-006
2405.17605
16 H. Qu and L. Gouskos ParticleNet: Jet tagging via particle clouds PRD 101 (2020) 056019 1902.08570
17 CMS Collaboration Electron and photon reconstruction and identification with the CMS experiment at the CERN LHC JINST 16 (2021) P05014 CMS-EGM-17-001
2012.06888
18 CMS Collaboration HEPData record for this analysis link
19 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004
20 CMS Collaboration Development of the CMS detector for the CERN LHC Run 3 JINST 19 (2024) P05064
21 CMS Collaboration Performance of the CMS Level-1 trigger in proton-proton collisions at $ \sqrt{s} = $ 13 TeV JINST 15 (2020) P10017 CMS-TRG-17-001
2006.10165
22 CMS Collaboration The CMS trigger system JINST 12 (2017) P01020 CMS-TRG-12-001
1609.02366
23 CMS Collaboration Performance of the CMS high-level trigger during LHC Run 2 JINST 19 (2024) P11021 CMS-TRG-19-001
2410.17038
24 CMS Collaboration Performance of the CMS muon detector and muon reconstruction with proton-proton collisions at $ \sqrt{s}= $ 13 TeV JINST 13 (2018) P06015 CMS-MUO-16-001
1804.04528
25 CMS Collaboration Description and performance of track and primary-vertex reconstruction with the CMS tracker JINST 9 (2014) P10009 CMS-TRK-11-001
1405.6569
26 CMS Collaboration Search for resonant and nonresonant new phenomena in high-mass dilepton final states at $ \sqrt{s}= $ 13 TeV JHEP 07 (2021) 208 CMS-EXO-19-019
2103.02708
27 CMS Collaboration Particle-flow reconstruction and global event description with the CMS detector JINST 12 (2017) P10003 CMS-PRF-14-001
1706.04965
28 A. Hoecker et al. TMVA: Toolkit for multivariate data analysis PoS ACAT 04 (2007) 0 physics/0703039
29 CMS Collaboration Performance of photon reconstruction and identification with the CMS detector in proton-proton collisions at $ \sqrt{s} = $ 8 TeV JINST 10 (2015) P08010 CMS-EGM-14-001
1502.02702
30 CMS Collaboration A measurement of the Higgs boson mass in the diphoton decay channel PLB 805 (2020) 135425 CMS-HIG-19-004
2002.06398
31 M. Cacciari, G. P. Salam, and G. Soyez The anti-$ k_{\mathrm{T}} $ jet clustering algorithm JHEP 04 (2008) 063 0802.1189
32 M. Cacciari, G. P. Salam, and G. Soyez FastJet user manual EPJC 72 (2012) 1896 1111.6097
33 CMS Collaboration Pileup mitigation at CMS in 13 TeV data JINST 15 (2020) P09018 CMS-JME-18-001
2003.00503
34 CMS Collaboration Jet energy scale and resolution in the CMS experiment in pp collisions at 8 TeV JINST 12 (2017) P02014 CMS-JME-13-004
1607.03663
35 M. Dasgupta, A. Fregoso, S. Marzani, and G. P. Salam Towards an understanding of jet substructure JHEP 09 (2013) 029 1307.0007
36 J. M. Butterworth, A. R. Davison, M. Rubin, and G. P. Salam Jet substructure as a new Higgs search channel at the LHC PRL 100 (2008) 242001 0802.2470
37 A. J. Larkoski, S. Marzani, G. Soyez, and J. Thaler Soft drop JHEP 05 (2014) 146 1402.2657
38 J. Alwall et al. The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations JHEP 07 (2014) 079 1405.0301
39 A. Banerjee et al. Phenomenological aspects of composite Higgs scenarios: exotic scalars and vector-like quarks Snowmass 2021 2203.07270
40 S. Frixione, P. Nason, and C. Oleari Matching NLO QCD computations with parton shower simulations: the POWHEG method JHEP 11 (2007) 070 0709.2092
41 N. Collaboration Parton distributions from high-precision collider data EPJC 77 (2017) 663 1706.00428
42 T. Sjostrand et al. An introduction to PYTHIA 8.2 Comput. Phys. Commun. 191 (2015) 159 1410.3012
43 CMS Collaboration Extraction and validation of a new set of CMS PYTHIA8 tunes from underlying-event measurements EPJC 80 (2020) 4 CMS-GEN-17-001
1903.12179
44 GEANT4 Collaboration GEANT 4---a simulation toolkit NIM A 506 (2003) 250
45 CMS Collaboration The CMS statistical analysis and combination tool: COMBINE Comput. Softw. Big Sci. 8 (2024) 19 CMS-CAT-23-001
2404.06614
46 W. Verkerke and D. Kirkby The RooFit toolkit for data modeling in 13th Int. Conf. on Computing in High Energy and Nuclear Physics, La Jolla, CA, March 24--28
Proc. 1 (2003) 3
physics/0306116
47 L. Moneta et al. The RooStats project in 13th Int. Workshop on Advanced Computing and Analysis Techniques in Physics Research, Jaipur, India, February 22--27
Proc. 1 (2010) 3
1009.1003
48 CMS Collaboration Precision luminosity measurement in proton-proton collisions at $ \sqrt{s} = $ 13 TeV in 2015 and 2016 at CMS EPJC 81 (2021) 800 CMS-LUM-17-003
2104.01927
49 CMS Collaboration Precision luminosity measurement in proton-proton collisions at a center-of-mass energy of 13 TeV with the CMS detector at the Large Hadron Collider Submitted to Phys.Rev.X, 2026 CMS-LUM-20-001
2606.26832
50 CMS Collaboration Measurement of the inelastic proton-proton cross section at $ \sqrt{s}= $ 13 TeV JHEP 07 (2018) 161 CMS-FSQ-15-005
1802.02613
51 J. Butterworth et al. PDF4LHC recommendations for LHC Run II JPG 43 (2016) 023001 1510.03865
52 C. Clopper and E. Pearson The use of confidence or fiducial limits illustrated in the case of the binomial Biometrika 26 (1934) 404
53 A. L. Read Presentation of search results: The CL$ _{\text{s}} $ technique JPG 28 (2002) 2693
54 T. Junk Confidence level computation for combining searches with small statistics NIM A 434 (1999) 435 hep-ex/9902006
55 G. Cowan, K. Cranmer, E. Gross, and O. Vitells Asymptotic formulae for likelihood-based tests of new physics Erratum:
EPJC 71 (2011) 1554
1007.1727
56 Particle Data Group Collaboration Review of particle physics PRD 110 (2024) 030001
57 M. Czakon and A. Mitov Top++: A program for the calculation of the top-pair cross-section at hadron colliders Comput. Phys. Commun. 185 (2014) 2930 1112.5675
58 M. Czakon, P. Fiedler, and A. Mitov Total top-quark pair-production cross section at hadron colliders PRL 110 (2013) 252004 1303.6254
Compact Muon Solenoid
LHC, CERN