CMSPASB2G19001  
Search for production of single vectorlike quarks decaying to tH or tZ in the allhadronic final state in pp collisions at $ \sqrt{s} = $ 13 TeV  
CMS Collaboration  
1 September 2023  
Abstract: A search for electroweak production of a single vectorlike quark T in association with a bottom quark (b) and decaying fully hadronically is presented. This search uses protonproton collision data at $ \sqrt{s} = $ 13 TeV collected by the CMS experiment at the CERN LHC during 20162018, corresponding to an integrated luminosity of 138 fb$ ^{1} $. The T can have charge 2/3 and can decay to a top quark (t) and a Higgs or Z boson. Event kinematics and the presence of jets containing b hadrons are used to reconstruct the hadronic decays of the t and Higgs or Z boson. No discrepancy from the standard model prediction is observed in the data. For T masses from 600 to 1200 GeV, the 95% confidence level upper limits on the production cross section of a T produced in association with a b and decaying via a t and a Higgs or Z boson range from 1260 to 68 fb.  
Links:
CDS record (PDF) ;
CADI line (restricted) ;
These preliminary results are superseded in this paper, Submitted to PRD. The superseded preliminary plots can be found here. 
Figures & Tables  Summary  Additional Figures  References  CMS Publications 

Figures  
png pdf 
Figure 1:
Example Feynman diagram for electroweak production of a vectorlike T. 
png pdf 
Figure 2:
The fivejet invariant mass distribution in the 2M1L region after the highmass (green crosses) and lowmass (black circles) selections in 2018 dataset. The lowmass selection results in a mass distribution that is smoothly falling, unlike the highmass selection. The highmass selection is more efficient for signal T masses above 700 GeV. 
png pdf 
Figure 3:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2018 data for the lowmass (upper) and highmass (lower) selections. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel and other years, similar functions are derived. 
png pdf 
Figure 3a:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2018 data for the lowmass (upper) and highmass (lower) selections. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel and other years, similar functions are derived. 
png pdf 
Figure 3b:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2018 data for the lowmass (upper) and highmass (lower) selections. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel and other years, similar functions are derived. 
png pdf 
Figure 3c:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2018 data for the lowmass (upper) and highmass (lower) selections. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel and other years, similar functions are derived. 
png pdf 
Figure 3d:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2018 data for the lowmass (upper) and highmass (lower) selections. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel and other years, similar functions are derived. 
png pdf 
Figure 4:
The fivejet invariant mass distribution in the tH channel (black markers) after the highmass selection in the multijet 3T validation region (upper left), the $ \mathrm{t} \overline{\mathrm{t}} $ 2T1L validation region (upper right) and the 3M signal region (lower) for 2018 dataset. The histograms are the corresponding reweighted 2M1L distributions. The background distribution is normalized to the number of entries in the data. The shaded area corresponds to the statistical uncertainties in the 2M1L regions. A potential 900 GeV T signal (red crosshatched histogram) is added to the background histogram demonstrating a negligible contribution. Similar results are observed in the tZ channel, and for the other years, but with slightly larger statistical uncertainties. 
png pdf 
Figure 4a:
The fivejet invariant mass distribution in the tH channel (black markers) after the highmass selection in the multijet 3T validation region (upper left), the $ \mathrm{t} \overline{\mathrm{t}} $ 2T1L validation region (upper right) and the 3M signal region (lower) for 2018 dataset. The histograms are the corresponding reweighted 2M1L distributions. The background distribution is normalized to the number of entries in the data. The shaded area corresponds to the statistical uncertainties in the 2M1L regions. A potential 900 GeV T signal (red crosshatched histogram) is added to the background histogram demonstrating a negligible contribution. Similar results are observed in the tZ channel, and for the other years, but with slightly larger statistical uncertainties. 
png pdf 
Figure 4b:
The fivejet invariant mass distribution in the tH channel (black markers) after the highmass selection in the multijet 3T validation region (upper left), the $ \mathrm{t} \overline{\mathrm{t}} $ 2T1L validation region (upper right) and the 3M signal region (lower) for 2018 dataset. The histograms are the corresponding reweighted 2M1L distributions. The background distribution is normalized to the number of entries in the data. The shaded area corresponds to the statistical uncertainties in the 2M1L regions. A potential 900 GeV T signal (red crosshatched histogram) is added to the background histogram demonstrating a negligible contribution. Similar results are observed in the tZ channel, and for the other years, but with slightly larger statistical uncertainties. 
png pdf 
Figure 4c:
The fivejet invariant mass distribution in the tH channel (black markers) after the highmass selection in the multijet 3T validation region (upper left), the $ \mathrm{t} \overline{\mathrm{t}} $ 2T1L validation region (upper right) and the 3M signal region (lower) for 2018 dataset. The histograms are the corresponding reweighted 2M1L distributions. The background distribution is normalized to the number of entries in the data. The shaded area corresponds to the statistical uncertainties in the 2M1L regions. A potential 900 GeV T signal (red crosshatched histogram) is added to the background histogram demonstrating a negligible contribution. Similar results are observed in the tZ channel, and for the other years, but with slightly larger statistical uncertainties. 
png pdf 
Figure 5:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The fit is performed on the combined data from all 3 years in the alltH channel. 
png pdf 
Figure 5a:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The fit is performed on the combined data from all 3 years in the alltH channel. 
png pdf 
Figure 5b:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The fit is performed on the combined data from all 3 years in the alltH channel. 
png pdf 
Figure 5c:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The fit is performed on the combined data from all 3 years in the alltH channel. 
png pdf 
Figure 5d:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The fit is performed on the combined data from all 3 years in the alltH channel. 
png pdf 
Figure 5e:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The fit is performed on the combined data from all 3 years in the alltH channel. 
png pdf 
Figure 5f:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The fit is performed on the combined data from all 3 years in the alltH channel. 
png pdf 
Figure 6:
Observed pvalues when considering the tH channel for each year and their combination. The excess observed in 2016 data is not confirmed in the other years. 
png pdf 
Figure 7:
The observed and expected 95% CL limits on the cross section for associated production with a b for final states $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $ (upper left), $ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (upper right), their sum $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (lower left), and $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ including the leakage of tH events into the tZ channel (bottom right) for different assumed values of the T mass. The vertical dashed line represents the crossover point in sensitivity: for masses to the left, the lowmass selection is used to set limits, while for masses to the right, the highmass selection is used to set limits. The red lines indicate the theoretical cross section for the singlet model. 
png pdf 
Figure 7a:
The observed and expected 95% CL limits on the cross section for associated production with a b for final states $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $ (upper left), $ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (upper right), their sum $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (lower left), and $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ including the leakage of tH events into the tZ channel (bottom right) for different assumed values of the T mass. The vertical dashed line represents the crossover point in sensitivity: for masses to the left, the lowmass selection is used to set limits, while for masses to the right, the highmass selection is used to set limits. The red lines indicate the theoretical cross section for the singlet model. 
png pdf 
Figure 7b:
The observed and expected 95% CL limits on the cross section for associated production with a b for final states $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $ (upper left), $ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (upper right), their sum $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (lower left), and $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ including the leakage of tH events into the tZ channel (bottom right) for different assumed values of the T mass. The vertical dashed line represents the crossover point in sensitivity: for masses to the left, the lowmass selection is used to set limits, while for masses to the right, the highmass selection is used to set limits. The red lines indicate the theoretical cross section for the singlet model. 
png pdf 
Figure 7c:
The observed and expected 95% CL limits on the cross section for associated production with a b for final states $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $ (upper left), $ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (upper right), their sum $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (lower left), and $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ including the leakage of tH events into the tZ channel (bottom right) for different assumed values of the T mass. The vertical dashed line represents the crossover point in sensitivity: for masses to the left, the lowmass selection is used to set limits, while for masses to the right, the highmass selection is used to set limits. The red lines indicate the theoretical cross section for the singlet model. 
png pdf 
Figure 7d:
The observed and expected 95% CL limits on the cross section for associated production with a b for final states $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $ (upper left), $ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (upper right), their sum $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (lower left), and $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ including the leakage of tH events into the tZ channel (bottom right) for different assumed values of the T mass. The vertical dashed line represents the crossover point in sensitivity: for masses to the left, the lowmass selection is used to set limits, while for masses to the right, the highmass selection is used to set limits. The red lines indicate the theoretical cross section for the singlet model. 
png pdf 
Figure 8:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the lowmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Figure 8a:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the lowmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Figure 8b:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the lowmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Figure 8c:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the lowmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Figure 8d:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the lowmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Figure 9:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the highmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Figure 9a:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the highmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Figure 9b:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the highmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Figure 9c:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the highmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Figure 9d:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the highmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Figure 10:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 10a:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 10b:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 10c:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 10d:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 10e:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 10f:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 11:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 11a:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 11b:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 11c:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 11d:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 11e:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 11f:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 12:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 12a:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 12b:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 12c:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 12d:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 12e:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 12f:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 13:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 13a:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 13b:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 13c:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 13d:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 13e:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 13f:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Figure 14:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass selection (left) and the highmass selection (right). The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 14a:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass selection (left) and the highmass selection (right). The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 14b:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass selection (left) and the highmass selection (right). The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 14c:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass selection (left) and the highmass selection (right). The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 14d:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass selection (left) and the highmass selection (right). The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 14e:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass selection (left) and the highmass selection (right). The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 14f:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass selection (left) and the highmass selection (right). The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 15:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 15a:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 15b:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 15c:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 15d:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 15e:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 15f:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 16:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 16a:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 16b:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 16c:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 16d:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 16e:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 16f:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Figure 17:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The tZ channel is shown when the fit is performed on the combined data from all three years in the tZ and tH channels. 
png pdf 
Figure 17a:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The tZ channel is shown when the fit is performed on the combined data from all three years in the tZ and tH channels. 
png pdf 
Figure 17b:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The tZ channel is shown when the fit is performed on the combined data from all three years in the tZ and tH channels. 
png pdf 
Figure 17c:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The tZ channel is shown when the fit is performed on the combined data from all three years in the tZ and tH channels. 
png pdf 
Figure 17d:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The tZ channel is shown when the fit is performed on the combined data from all three years in the tZ and tH channels. 
png pdf 
Figure 17e:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The tZ channel is shown when the fit is performed on the combined data from all three years in the tZ and tH channels. 
png pdf 
Figure 17f:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The tZ channel is shown when the fit is performed on the combined data from all three years in the tZ and tH channels. 
png pdf 
Figure 18:
The observed and expected 95% CL limits on the cross section for associated production with a b for final states $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $ (upper left), $ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (upper right), their sum $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (lower left), and $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ including the leakage of tH events into the tZ channel (bottom right) for different assumed values of the T mass in logarithmic scale. The vertical dashed line represents the crossover point in sensitivity: for masses to the left, the lowmass selection is used to set limits, while for masses to the right, the highmass selection is used to set limits. The red lines indicate the theoretical cross section for the singlet model. 
png pdf 
Figure 18a:
The observed and expected 95% CL limits on the cross section for associated production with a b for final states $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $ (upper left), $ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (upper right), their sum $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (lower left), and $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ including the leakage of tH events into the tZ channel (bottom right) for different assumed values of the T mass in logarithmic scale. The vertical dashed line represents the crossover point in sensitivity: for masses to the left, the lowmass selection is used to set limits, while for masses to the right, the highmass selection is used to set limits. The red lines indicate the theoretical cross section for the singlet model. 
png pdf 
Figure 18b:
The observed and expected 95% CL limits on the cross section for associated production with a b for final states $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $ (upper left), $ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (upper right), their sum $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (lower left), and $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ including the leakage of tH events into the tZ channel (bottom right) for different assumed values of the T mass in logarithmic scale. The vertical dashed line represents the crossover point in sensitivity: for masses to the left, the lowmass selection is used to set limits, while for masses to the right, the highmass selection is used to set limits. The red lines indicate the theoretical cross section for the singlet model. 
png pdf 
Figure 18c:
The observed and expected 95% CL limits on the cross section for associated production with a b for final states $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $ (upper left), $ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (upper right), their sum $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (lower left), and $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ including the leakage of tH events into the tZ channel (bottom right) for different assumed values of the T mass in logarithmic scale. The vertical dashed line represents the crossover point in sensitivity: for masses to the left, the lowmass selection is used to set limits, while for masses to the right, the highmass selection is used to set limits. The red lines indicate the theoretical cross section for the singlet model. 
png pdf 
Figure 18d:
The observed and expected 95% CL limits on the cross section for associated production with a b for final states $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $ (upper left), $ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (upper right), their sum $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (lower left), and $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ including the leakage of tH events into the tZ channel (bottom right) for different assumed values of the T mass in logarithmic scale. The vertical dashed line represents the crossover point in sensitivity: for masses to the left, the lowmass selection is used to set limits, while for masses to the right, the highmass selection is used to set limits. The red lines indicate the theoretical cross section for the singlet model. 
Tables  
png pdf 
Table 1:
Mean and standard deviation values from a Gaussian fit of the Higgs, Z, and W boson and top quark mass distributions in the 700 GeV T quark sample, requiring the jet kinematic criteria described above and matching to generated particles. All quantities are in units of GeVns. The yeartoyear variations are within the jet energy scale uncertainties. 
png pdf 
Table 2:
Cumulative efficiencies of the highmass selection criteria for signal and various simulated backgrounds. The first and last lines indicate the expected number of events normalized to an integrated luminosity of 35.9 fb$ ^{1} $. Only statistical uncertainties are reported. The ``Other backgrounds'' column includes W/Z+jets, single t, $ {\mathrm{t}\overline{\mathrm{t}}} $H, and $ {\mathrm{t}\overline{\mathrm{t}}} $Z background processes. The $ \mathrm{t} \overline{\mathrm{t}} $H and $ \mathrm{t} \overline{\mathrm{t}} $Z processes do not form a resonance and have production rates roughly the same as the signal rate. 
png pdf 
Table 3:
Summary of the signal regions (highmass selection) and validation samples definitions. The changes with respect to the 3T signal region are highlighted in bold. The ``T', ``M'', and ``L'' labels represent tight, medium, and loose btagged jet requirements, respectively. 
Summary 
A search for a vectorlike top quark partner T in the single production mode was performed using protonproton collision events at $ \sqrt{s} = $ 13 TeV collected by the CMS experiment in 2016, 2017, and 2018. In this search, the T is assumed to only couple to standard model thirdgeneration quarks. We consider signatures containing a top quark and a Higgs (tH) or Z (tZ) boson decaying to a bottom quarkantiquark pair. The major background processes are top quarkantiquark pair and multijet production. The excess in the tH final state found in an earlier search [23] is not observed with a larger dataset and improved event selection. We set 95% confidence level upper limits on product of the cross section $ \sigma $ and branching fractions $ \mathcal{B}(\mathrm{T}\to\mathrm{t}\mathrm{H}) $ and $ \mathcal{B}(\mathrm{T}\to\mathrm{t}\mathrm{Z} $) ranging from 1260 to 68 fb for T masses from 600 to 1200 GeV in the $ \mathrm{p}\mathrm{p}\to\mathrm{T}\mathrm{b}\mathrm{q} $ production channel. The limits are stronger than those in the previous search by at least a factor of three. 
Additional Figures  
png pdf 
Additional Figure 1:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the lowmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Additional Figure 1a:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the lowmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Additional Figure 1b:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the lowmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Additional Figure 1c:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the lowmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Additional Figure 1d:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the lowmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Additional Figure 2:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the highmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Additional Figure 2a:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the highmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Additional Figure 2b:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the highmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Additional Figure 2c:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the highmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Additional Figure 2d:
Weights from b tagging efficiency ratios (open markers) as a function of the fivejet invariant mass in 2016 (upper) and 2017 (lower) data for the highmass selection. The left graph shows weights connecting the 2M1L and 3M regions, and the right graph shows weights connecting the 3M and 3T regions. The red line corresponds to the central value of the transfer function and the shaded area represents the 95% confidence level uncertainty band. For the lowmass (highmass) analysis only signals with mass below (above) 800 GeV are tested, so primarily the lower (upper) part of the distribution contributes to the final result. In the tZ channel, similar functions are derived. 
png pdf 
Additional Figure 3:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 3a:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 3b:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 3c:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 3d:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 3e:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 3f:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 4:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 4a:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 4b:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 4c:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 4d:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 4e:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 4f:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 5:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 5a:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 5b:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 5c:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 5d:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 5e:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 5f:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 6:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 6a:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 6b:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 6c:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 6d:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 6e:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 6f:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the (top) 2M1L, (middle) 3M, and (bottom) 3T regions for the lowmass (left) and highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. The dashed vertical line denotes the division between the lowmass and highmass selections. 
png pdf 
Additional Figure 7:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass selection (left) and the highmass selection (right). The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 7a:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass selection (left) and the highmass selection (right). The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 7b:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass selection (left) and the highmass selection (right). The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 7c:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass selection (left) and the highmass selection (right). The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 7d:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass selection (left) and the highmass selection (right). The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 7e:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass selection (left) and the highmass selection (right). The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 7f:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2016 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass selection (left) and the highmass selection (right). The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 8:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 8a:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 8b:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 8c:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 8d:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 8e:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 8f:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2017 dataset (black markers) in the 2M1L (lower), 3M (middle), and 3T (upper) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 9:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 9a:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 9b:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 9c:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 9d:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 9e:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 9f:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the 2018 dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the bestfit signal distribution from a 700 GeV T. The fit is performed on the combined data from all three years in the alltH channel. 
png pdf 
Additional Figure 10:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The tZ channel is shown when the fit is performed on the combined data from all three years in the tZ and tH channels. 
png pdf 
Additional Figure 10a:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The tZ channel is shown when the fit is performed on the combined data from all three years in the tZ and tH channels. 
png pdf 
Additional Figure 10b:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The tZ channel is shown when the fit is performed on the combined data from all three years in the tZ and tH channels. 
png pdf 
Additional Figure 10c:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The tZ channel is shown when the fit is performed on the combined data from all three years in the tZ and tH channels. 
png pdf 
Additional Figure 10d:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The tZ channel is shown when the fit is performed on the combined data from all three years in the tZ and tH channels. 
png pdf 
Additional Figure 10e:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The tZ channel is shown when the fit is performed on the combined data from all three years in the tZ and tH channels. 
png pdf 
Additional Figure 10f:
Fivejet invariant mass distributions after a backgroundonly fit (blue histogram) to the complete dataset (black markers) in the 2M1L (upper), 3M (middle), and 3T (lower) regions for lowmass (left) and the highmass (right) selections. The dashed blue band represents the uncertainty on the fitted background estimate, and red dashed line shows the expected signal distribution for a 700 GeV (lowmass selection) and a 900 GeV (highmass selection) T. The tZ channel is shown when the fit is performed on the combined data from all three years in the tZ and tH channels. 
png pdf 
Additional Figure 11:
The observed and expected 95% CL limits on the cross section for associated production with a b for final states $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $ (upper left), $ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (upper right), their sum $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (lower left), and $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ including the leakage of tH events into the tZ channel (bottom right) for different assumed values of the T mass in logarithmic scale. The vertical dashed line represents the crossover point in sensitivity: for masses to the left, the lowmass selection is used to set limits, while for masses to the right, the highmass selection is used to set limits. The red lines indicate the theoretical cross section for the singlet model. 
png pdf 
Additional Figure 11a:
The observed and expected 95% CL limits on the cross section for associated production with a b for final states $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $ (upper left), $ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (upper right), their sum $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (lower left), and $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ including the leakage of tH events into the tZ channel (bottom right) for different assumed values of the T mass in logarithmic scale. The vertical dashed line represents the crossover point in sensitivity: for masses to the left, the lowmass selection is used to set limits, while for masses to the right, the highmass selection is used to set limits. The red lines indicate the theoretical cross section for the singlet model. 
png pdf 
Additional Figure 11b:
The observed and expected 95% CL limits on the cross section for associated production with a b for final states $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $ (upper left), $ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (upper right), their sum $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (lower left), and $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ including the leakage of tH events into the tZ channel (bottom right) for different assumed values of the T mass in logarithmic scale. The vertical dashed line represents the crossover point in sensitivity: for masses to the left, the lowmass selection is used to set limits, while for masses to the right, the highmass selection is used to set limits. The red lines indicate the theoretical cross section for the singlet model. 
png pdf 
Additional Figure 11c:
The observed and expected 95% CL limits on the cross section for associated production with a b for final states $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $ (upper left), $ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (upper right), their sum $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (lower left), and $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ including the leakage of tH events into the tZ channel (bottom right) for different assumed values of the T mass in logarithmic scale. The vertical dashed line represents the crossover point in sensitivity: for masses to the left, the lowmass selection is used to set limits, while for masses to the right, the highmass selection is used to set limits. The red lines indicate the theoretical cross section for the singlet model. 
png pdf 
Additional Figure 11d:
The observed and expected 95% CL limits on the cross section for associated production with a b for final states $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $ (upper left), $ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (upper right), their sum $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ (lower left), and $ \mathrm{t}\mathrm{H}\mathrm{b}\mathrm{q} $+$ \mathrm{t}\mathrm{Z}\mathrm{b}\mathrm{q} $ including the leakage of tH events into the tZ channel (bottom right) for different assumed values of the T mass in logarithmic scale. The vertical dashed line represents the crossover point in sensitivity: for masses to the left, the lowmass selection is used to set limits, while for masses to the right, the highmass selection is used to set limits. The red lines indicate the theoretical cross section for the singlet model. 
References  
1  Y. Okada and L. Panizzi  LHC signatures of vectorlike quarks  Adv. High Energy Phys. 2013 (2013) 364936  1207.5607 
2  J. A. AguilarSaavedra, R. Benbrik, S. Heinemeyer, and M. PérezVictoria  Handbook of vectorlike quarks: Mixing and single production  PRD 88 (2013) 094010  1306.0572 
3  A. De Simone, O. Matsedonskyi, R. Rattazzi, and A. Wulzer  A first top partner hunter's guide  JHEP 04 (2013) 004  1211.5663 
4  M. Buchkremer, G. Cacciapaglia, A. Deandrea, and L. Panizzi  Model independent framework for searches of top partners  NPB 876 (2013) 376  1305.4172 
5  A. Djouadi and A. Lenz  Sealing the fate of a fourth generation of fermions  PLB 715 (2012) 310  1204.1252 
6  A. Lenz  Constraints on a fourth generation of fermions from Higgs Boson searches  Adv. High Energy Phys. 2013 (2013) 910275  
7  CMS Collaboration  Combined search for the quarks of a sequential fourth generation  PRD 86 (2012) 112003  CMSEXO11098 1209.1062 
8  N. Vignaroli  Early discovery of top partners and test of the Higgs nature  PRD 86 (2012) 075017  1207.0830 
9  J. A. AguilarSaavedra  Identifying top partners at LHC  JHEP 11 (2009) 030  0907.3155 
10  CMS Collaboration  Search for pair production of vectorlike quarks in the $ \mathrm{b}\mathrm{W}\overline{\mathrm{b}}\mathrm{W} $ channel from protonproton collisions at $ \sqrt{s} = $ 13 TeV  PLB 779 (2018) 82  1710.01539 
11  ATLAS Collaboration  Search for pair production of heavy vectorlike quarks decaying into hadronic final states in $ pp $ collisions at $ \sqrt{s} = $ 13 TeV with the ATLAS detector  PRD 98 (2018) 092005  1808.01771 
12  ATLAS Collaboration  Search for pairproduced vectorlike top and bottom partners in events with large missing transverse momentum in pp collisions with the ATLAS detector  EPJC 83 (2023) 719  2212.05263 
13  CMS Collaboration  Search for pair production of vectorlike quarks in leptonic final states in protonproton collisions at $ \sqrt{s} $ = 13 TeV  JHEP 07 (2023) 020  2209.07327 
14  ATLAS Collaboration  Search for pairproduction of vectorlike quarks in pp collision events at s=13 TeV with at least one leptonically decaying Z boson and a thirdgeneration quark with the ATLAS detector  PLB 843 (2023) 138019  2210.15413 
15  ATLAS Collaboration  Search for single production of vectorlike quarks decaying into $ Wb $ in $ pp $ collisions at $ \sqrt{s} = $ 13 TeV with the ATLAS detector  JHEP 05 (2019) 164  1812.07343 
16  ATLAS Collaboration  Search for large missing transverse momentum in association with one topquark in protonproton collisions at $ \sqrt{s} $ = 13 TeV with the ATLAS detector  JHEP 05 (2019) 041  1812.09743 
17  ATLAS Collaboration  Search for single production of a vectorlike $ T $ quark decaying into a Higgs boson and top quark with fully hadronic final states using the ATLAS detector  PRD 105 (2022) 092012  2201.07045 
18  CMS Collaboration  Search for single production of vectorlike quarks decaying into a b quark and a W boson in protonproton collisions at $ \sqrt{s} = $ 13 TeV  PLB 772 (2017) 634  1701.08328 
19  CMS Collaboration  Search for single production of a vectorlike T quark decaying to a Z boson and a top quark in protonproton collisions at $ \sqrt{s} = $ 13 TeV  PLB 781 (2018) 574  1708.01062 
20  CMS Collaboration  Search for single production of a vectorlike T quark decaying to a top quark and a Z boson in the final state with jets and missing transverse momentum at $ \sqrt{s} $ = 13 TeV  JHEP 05 (2022) 093  2201.02227 
21  CMS Collaboration  Search for single production of a heavy vectorlike T quark decaying to a Higgs boson and a top quark with a lepton and jets in the final state  PLB 771 (2017) 80  1612.00999 
22  CMS Collaboration  Search for electroweak production of a vectorlike quark decaying to a top quark and a Higgs boson using boosted topologies in fully hadronic final states  JHEP 04 (2017) 136  1612.05336 
23  CMS Collaboration  Search for electroweak production of a vectorlike T quark using fully hadronic final states  JHEP 01 (2020) 036  1909.04721 
24  CMS Collaboration  Description and performance of track and primaryvertex reconstruction with the CMS tracker  JINST 9 (2014) P10009  CMSTRK11001 1405.6569 
25  CMS Collaboration  Performance of photon reconstruction and identification with the CMS detector in protonproton collisions at $ \sqrt{s} = $ 8 TeV  JINST 10 (2015) P08010  CMSEGM14001 1502.02702 
26  CMS Collaboration  The CMS experiment at the CERN LHC  JINST 3 (2008) S08004  
27  CMS Collaboration  The CMS trigger system  JINST 12 (2017) P01020  CMSTRG12001 1609.02366 
28  CMS Collaboration  Technical proposal for the PhaseII upgrade of the Compact Muon Solenoid  CMS Technical Proposal CERNLHCC2015010, CMSTDR1502, 2015 CDS 

29  CMS Collaboration  Particleflow reconstruction and global event description with the CMS detector  JINST 12 (2017) P10003  CMSPRF14001 1706.04965 
30  M. Cacciari, G. P. Salam, and G. Soyez  The anti$ k_{\mathrm{T}} $ jet clustering algorithm  JHEP 04 (2008) 063  0802.1189 
31  M. Cacciari, G. P. Salam, and G. Soyez  FastJet user manual  EPJC 72 (2012) 1896  1111.6097 
32  CMS Collaboration  Pileup mitigation at CMS in 13 TeV data  JINST 15 (2020) P09018  CMSJME18001 2003.00503 
33  CMS Collaboration  Jet energy scale and resolution in the CMS experiment in pp collisions at 8 TeV  JINST 12 (2017) P02014  CMSJME13004 1607.03663 
34  CMS Collaboration  Identification of heavyflavour jets with the CMS detector in pp collisions at 13 TeV  JINST 13 (2018) P05011  CMSBTV16002 1712.07158 
35  J. Alwall et al.  The automated computation of treelevel and nexttoleading order differential cross sections, and their matching to parton shower simulations  JHEP 07 (2014) 079  1405.0301 
36  T. Sjöstrand et al.  An Introduction to PYTHIA8.2  Comput. Phys. Commun. 191 (2015) 159  1410.3012 
37  CMS Collaboration  Event generator tunes obtained from underlying event and multiparton scattering measurements  EPJC 76 (2016) 155  CMSGEN14001 1512.00815 
38  CMS Collaboration  Extraction and validation of a new set of CMS PYTHIA8 tunes from underlyingevent measurements  EPJC 80 (2020) 4  CMSGEN17001 1903.12179 
39  NNPDF Collaboration  Parton distributions for the LHC Run II  JHEP 04 (2015) 040  1410.8849 
40  NNPDF Collaboration  Parton distributions from highprecision collider data  EPJC 77 (2017) 663  1706.00428 
41  GEANT4 Collaboration  GEANT 4a simulation toolkit  NIM A 506 (2003) 250  
42  CMS Collaboration  Precision luminosity measurement in protonproton collisions at $ \sqrt{s} = $ 13 TeV in 2015 and 2016 at CMS  EPJC 81 (2021) 800  CMSLUM17003 2104.01927 
43  CMS Collaboration  CMS luminosity measurements for the 2017 data taking period  Technical Report , CERN, 2017 CMSPASLUM17004 
CMSPASLUM17004 
44  CMS Collaboration  CMS luminosity measurements for the 2018 data taking period  Technical Report , CERN, 2018 CMSPASLUM18002 
CMSPASLUM18002 
45  P. Nason  A new method for combining NLO QCD with shower Monte Carlo algorithms  JHEP 11 (2004) 040  hepph/0409146 
46  S. Frixione, P. Nason, and C. Oleari  Matching NLO QCD computations with parton shower simulations: the POWHEG method  JHEP 11 (2007) 070  0709.2092 
47  S. Alioli, P. Nason, C. Oleari, and E. Re  A general framework for implementing NLO calculations in shower Monte Carlo programs: the POWHEG BOX  JHEP 06 (2010) 043  1002.2581 
48  S. Frixione, P. Nason, and G. Ridolfi  A positiveweight nexttoleadingorder Monte Carlo for heavy flavour hadroproduction  JHEP 09 (2007) 126  0707.3088 
49  CMS Tracker Group Collaboration  The CMS phase1 pixel detector upgrade  JINST 16 (2021) P02027  2012.14304 
50  F. J. Massey Jr.  The KolmogorovSmirnov test for goodness of fit  J. Am. Stat. Assoc. 46 (1951) 68  
51  J. Butterworth et al.  PDF4LHC recommendations for LHC Run II  JPG 43 (2016) 023001  1510.03865 
52  L. Demortier  P values and nuisance parameters  in Statistical issues for LHC physics. Proceedings, Workshop, PHYSTATLHC, Geneva, 2007 link 

53  T. Junk  Confidence level computation for combining searches with small statistics  NIM A 434 (1999) 435  hepex/9902006 
54  A. L. Read  Presentation of search results: The $ \text{CL}_\text{s} $ technique  JPG 28 (2002) 2693  
55  ATLAS and CMS Collaborations, and LHC Higgs Combination Group  Procedure for the LHC Higgs boson search combination in Summer 2011  Technical Report CMSNOTE2011005, ATLPHYSPUB201111, 2011  
56  G. Cowan, K. Cranmer, E. Gross, and O. Vitells  Asymptotic formulae for likelihoodbased tests of new physics  EPJC 71 (2011) 1554  1007.1727 
Compact Muon Solenoid LHC, CERN 