Full simulation study of the top Yukawa coupling at the ILC at \(\sqrt{s}=1\) TeV

The discovery of a standard model (SM)-like Higgs boson, announced on July 4, 2012 by the ATLAS and CMS CFor SiD a superconductiollaborations [1, 2], was celebrated as a major milestone in particle physics. In the SM, the coupling strength of the Higgs boson to a fermion is given by \(y_{f}=\sqrt{2} m_{f}/v\), where \(m_f\) is the fermion mass and \(v\approx 246~\mathrm{GeV}\) is the vacuum expectation value. Since the top quark is the heaviest known elementary particle, the measurement of the top Yukawa coupling, \(y_{\mathrm {t}}\), serves as the high endpoint to test this prediction. A sizable deviation in \(y_{\mathrm {t}}\) from the SM prediction is expected in various new physics scenarios, which motivates a precise measurement of \(y_{\mathrm {t}}\). For example, in composite Higgs models, where the Higgs boson is a pseudo-Nambu–Goldstone boson, \(y_{\mathrm {t}}\) could deviate up to tens of %, even in the scenario that no new particles are discovered in LHC Run 2 data [3].

A recent study of the prospects of measuring \(y_{\mathrm {t}}\) at the LHC [4] estimates that a precision of 14–15 % (7–10 %) is achievable with an integrated luminosity of \(0.3 ~\mathrm {ab}^{-1} \) (\(3 ~\mathrm {ab}^{-1} \)), including theoretical and systematic uncertainties. For \(\mathrm{e}^{+}\mathrm{e}^{-} \) collisions, detailed simulation studies have been carried out using the \(\mathrm{t}\bar{\mathrm{t}}\mathrm{H}\) process at various center-of-mass energies. At \(\sqrt{s}=500\) GeV [5–7], where the \(\mathrm{e}^{+}\mathrm{e}^{-} \rightarrow \mathrm{t}\bar{\mathrm{t}}\mathrm{H}\) cross section is sharply rising, the statistical precision is estimated to be about 10 % for an integrated luminosity of \(1 ~\mathrm {ab}^{-1} \). At \(\sqrt{s}=800\) GeV [8, 9], it is estimated that \(y_{\mathrm {t}}\) can be measured to a precision of 5–6 % for an integrated luminosity of \(1~ \mathrm {ab}^{-1} \), including the systematic uncertainties due to the background normalization.

The International Linear Collider (ILC) [10] is a proposed \(\mathrm{e}^{+}\mathrm{e}^{-} \) collider with a maximum center-of-mass energy \(\sqrt{s} = 1~\mathrm{TeV}\). It has a broad physics potential that is complementary to the LHC and precision measurements of the Higgs couplings are an integral part of the physics program at this machine. We present studies of the measurement of the top Yukawa coupling in direct observation at the 1 TeV stage of the ILC. The studies are carried out in ILD and SiD [11], the two detector concepts for the ILC. They are performed with detailed detector simulations taking into account the main beam-induced backgrounds at the collider as well as the dominant background from other physics processes. Two final states are considered – events where both \({\mathrm{W}}\) bosons from the top quarks decay hadronically, and events where exactly one of the two \({\mathrm{W}}\) bosons decays leptonically.

The studies performed for the two detector concepts have large overlaps, and we highlight significant differences between the two analyses wherever applicable. This document is organized as follows: Sect. 2 gives an overview over the signal sample and the considered physics background. Section 3 gives brief overviews over the two ILC detector models. The tools for the generation of physics processes and the detector simulation and reconstruction are listed in Sect. 4. The two dominant sources of machine-induced background in the detectors are introduced in Sect. 5. The techniques to reduce these backgrounds and reconstruct the top quarks and Higgs bosons are described in Sect. 6. Details of the event selection are given in Sect. 7 and the results are presented in Sect. 8. The dominant sources of systematic uncertainty are given in Sect. 9 and the two analyses are summarized in Sect. 10.

Figure 1 illustrates the lowest order Feynman diagrams for the process \({{\mathrm{e}}}^{+}{{\mathrm{e}}}^{-} \rightarrow \mathrm{t}\bar{\mathrm{t}}\mathrm{H}\).

The lowest order Feynman diagrams for the process \({{\mathrm{e}}}^{+}{{\mathrm{e}}}^{-} \rightarrow \mathrm{t}\bar{\mathrm{t}}\mathrm{H}\). In a the Higgs boson is radiated from a top quark and b is the background Higgs-strahlung process where the Higgs boson is radiated from the Z boson

Full size image

The diagram for the Higgs-strahlung process \({{\mathrm{e}}}^{+}{{\mathrm{e}}}^{-} \rightarrow {{\mathrm{Z}}}^{*}{\mathrm{H}}{}{}\) with \({{\mathrm{Z}}}^{*}\rightarrow {\mathrm{t}}\bar{{\mathrm{t}}} \), which does not depend on \(y_{\mathrm {t}}\), has a small yet non-negligible contribution to the total cross section. The size of this effect is studied by evaluating the behavior of the \({{\mathrm{e}}}^{+}{{\mathrm{e}}}^{-} \rightarrow \mathrm{t}\bar{\mathrm{t}}\mathrm{H}\) cross section when changing \(y_{\mathrm {t}}\) from the SM value, using the linear approximation \({{\varDelta } y_{\mathrm {t}}}/{y_{\mathrm {t}}}=\kappa \cdot {{\varDelta } \sigma }/{\sigma }\). In the absence of the Higgs-strahlung diagram, we would find \(\kappa =0.5\). Instead, we find \(\kappa =0.52\), indicating a non-negligible contribution from the Higgs-strahlung diagram to the total cross section at \(\sqrt{s}=1~\mathrm{TeV}\). This factor is used in the extraction of the top Yukawa coupling precision. The correction will be known with good precision, because the Higgs coupling to the \({\mathrm{Z}}\) boson can be extracted from measurements of \(\mathrm{e}^{+}\mathrm{e}^{-} \rightarrow {\mathrm{Z}}{\mathrm{H}}{}{} \) events at \(\sqrt{s} = 250~\mathrm{GeV}\) with a statistical uncertainty of about 1.5 % [12].

For this study the semi-leptonic and hadronic decays of the \({\mathrm{t}}\bar{{\mathrm{t}}} \) system were studied with the Higgs decaying via the dominant decay mode into a \({\mathrm{b}}\bar{{\mathrm{b}}} \) pair. For the fully hadronic decay channel this leads to a signature of eight hadronic jets, four of which are \({\mathrm{b}}\) jets. In the semi-leptonic mode the final signal in the detector consists of six hadronic jets, four of which are \({\mathrm{b}}\) jets, an isolated lepton, and missing energy and momentum from a neutrino. For isolated leptons, only the prompt electrons and muons are reconstructed and considered as signal, neglecting the decays into \(\tau \) leptons. These two modes are reconstructed in independent samples and are combined statistically.

Irreducible backgrounds to these processes arise from the eight-fermion final states of \({\mathrm{t}}\bar{{\mathrm{t}}}{\mathrm{Z}} \) where the Z decays into a \({\mathrm{b}}\bar{{\mathrm{b}}} \) pair and \({\mathrm{t}}\bar{{\mathrm{t}}} {\mathrm{b}}\bar{{\mathrm{b}}} \) where the \({\mathrm{t}}\bar{{\mathrm{t}}} \) system radiates a hard gluon which forms a \({\mathrm{b}}\bar{{\mathrm{b}}} \) pair. A large background contribution also arises from \({\mathrm{t}}\bar{{\mathrm{t}}} \) due to the huge relative cross section compared to the signal. There is also a contribution from the other decay modes of the \(\mathrm{t}\bar{\mathrm{t}}\mathrm{H}\) system such as the Higgs decaying to final states other than a \({\mathrm{b}}\bar{{\mathrm{b}}} \) pair and the fully leptonic decays of the top quarks.

An overview of the cross sections (times branching ratio for the specifically listed final states) for the signal final states as well as for the considered backgrounds is shown in Table 1. For the measurement using the final state with six jets, all other \(\mathrm{t}\bar{\mathrm{t}}\mathrm{H}\) events, i.e., all events where both top quarks decay either leptonically or hadronically, or events where the Higgs boson does not decay into \({{\mathrm{b}}\bar{{\mathrm{b}}}} \), are treated as background. For the eight-jets final state events where at least one top quark decays leptonically or where the Higgs boson decays into final states other than \({{\mathrm{b}}\bar{{\mathrm{b}}}} \) are considered as background. The non-\(\mathrm{t}\bar{\mathrm{t}}\mathrm{H}\) backgrounds are considered for both measurements.

SiD [11, chapter 2] and ILD concepts [11, chapter 3] are designed to be the two general-purpose detectors for the ILC, with a \(4 \pi \) coverage, employing highly granular calorimeters for particle flow calorimetry.

For SiD a superconducting solenoid with an inner radius of 2.6 m provides a central magnetic field of 5 T. The calorimeters are placed inside the coil and consist of a 30 layer tungsten–silicon electromagnetic calorimeter (ECAL) with \(13~\mathrm{mm}^2\) segmentation, followed by a hadronic calorimeter (HCAL) with steel absorber and instrumented with resistive plate chambers (RPC) – 40 layers in the barrel region and 45 layers in the endcaps. The read-out cell size in the HCAL is \(10\times 10 ~\mathrm{mm}^2\). The iron return yoke outside of the coil is instrumented with 11 RPC layers with \(30\times 30~\mathrm{mm}^2\) read-out cells for muon identification. The silicon-only tracking system consists of five layers of \(20\times 20~\upmu \mathrm{m}^2\) pixels followed by five strip layers with a pitch of \(25~\upmu \mathrm{m}\), a read-out pitch of \(50~\upmu \mathrm{m}\) and a length of 92 mm per module in the barrel region. The tracking system in the endcap consists of four stereo-strip disks with similar pitch and a stereo angle of \(12^\circ \), complemented by four pixelated disks in the vertex region with a pixel size of \(20\times 20~\upmu \mathrm{m}^2\) and three disks in the far-forward region at lower radii with a pixel size of \(50\,{\times }\,50~\upmu \mathrm{m}^2\). All sub-detectors have the capability of time-stamping at the level of individual bunches, 337 ns apart, \(\approx \)1300 to a train. This allows to separate hits originating from different bunch crossings. The whole detector will be read out in the 200 ms between bunch trains.

The ILD detector model is designed around a different optimization with a larger size. The ECAL and HCAL are placed inside a superconducting solenoid, which provide a magnetic field of 3.5 T. The silicon-tungsten ECAL has an inner radius of 1.8 m and a total thickness of 20 cm, with \(5\times 5~\mathrm{mm}^2\) transverse cell size and 30 layers of longitudinal segmentation. The steel-scintillator HCAL has an outer radius of 3.4 m with \(3\times 3 \, \mathrm{cm}^2\) transverse tiles and 48 layers longitudinal segmentation. ILD employs a hybrid tracking system consisting of a time projection chamber (TPC) which provides up to 224 points per track and silicon-strip sensors for improved track momentum resolution, which are placed in the barrel region both inside and outside the TPC and in the endcap region outside the TPC. The vertex detector consists of three double layers of silicon pixel sensors with radii ranging from 15 to 60 mm, providing a spatial resolution of \(2.8~\upmu \mathrm{m}\). An iron return yoke instrumented with a muon detector and a tail catcher is placed outside the yoke. In addition, silicon trackers and beam/luminosity calorimeters are installed in the forward region.

The \(\mathrm{t}\bar{\mathrm{t}}\mathrm{H}\), \({\mathrm{t}}\bar{{\mathrm{t}}}{\mathrm{Z}} \), and \({\mathrm{t}}\bar{{\mathrm{t}}} {\mathrm{b}}\bar{{\mathrm{b}}} \) samples were generated using the physsim [13] event generator. The sample referred to as \({\mathrm{t}}\bar{{\mathrm{t}}} \) in the following includes six-fermion final states consistent with the \({\mathrm{t}}\bar{{\mathrm{t}}} \) decays but not limited to the resonant \({\mathrm{t}}\bar{{\mathrm{t}}} \) production. The \({\mathrm{t}}\bar{{\mathrm{t}}} \) events were generated using the whizard 1.95 [14, 15] event generator. All samples were generated taking into account the expected beam energy spectrum at the \(\sqrt{s}=1~\mathrm{TeV}\) ILC, including initial state radiation [16] and beamstrahlung. The spectrum was sampled from a simulation of beam events [17]. The model for the hadronization in pythia 6.4 [18] uses a tune based on opal data [19, Appendix B.3].

Detailed detector simulations based on geant4 [20, 21] are performed. In the SiD analysis, the event reconstruction is performed in the org.lcsim [22] package. The ILD analysis uses the Marlin [23, 24] framework. Both analyses use the PandoraPFA [25] algorithm for calorimeter clustering and combined analysis of track and calorimeter information based on the particle flow approach. The LCFIPlus package [11, Section 2.2.2.3] is used for the tagging of heavy flavor jets. The assumed integrated luminosity of the analysis is \(1~\mathrm {ab}^{-1} \), which is split equally between the two polarization configurations \((+80~\%, -20~\%)\) and \((-80~\%, +20~\%)\) for the polarization of the electron and positron beams \(P({{\mathrm{e}}}^{-},{{\mathrm{e}}}^{+})\). Detector hits from Beam-induced backgrounds from processes described in Sect. 5 are treated correctly in the simulation of the detector readout and in the reconstruction.

The ILC operating at \(\sqrt{s}=1~\mathrm{TeV}\) has an instantaneous luminosity of \(4.2 \times 10^{34}~{\text {cm}^{-2}~\mathrm {s}^{-1}}\). During the collision, a number of processes occur in addition to the primary scattering event. The production of incoherent electron-positron pairs results in an average of \(4.5\times 10^5\) low-momentum particles per bunch crossing. We assume an average of 4.1 hadronic events from two-photon processes (\({{\gamma }{\gamma }\rightarrow \mathrm {hadrons}} \)) with a diphoton-invariant mass greater than 300 MeV. The distributions of the particles originating from these processes in the (polar angle, energy) plane are shown in Fig. 2. They do not affect the reconstruction significantly, but present a challenge to the sub-detector occupancies and pattern recognition. The SiD analysis includes both effects, while the ILD analysis includes only the \({{\gamma }{\gamma }\rightarrow \mathrm {hadrons}} \) processes. The SiD analysis shows that incoherent pairs are under control by incorporating a detector design capable of time stamping individual bunch crossings. The baseline technology of the ILD barrel vertex detector integrates over 18 bunch crossings in layer 2 up to 180 bunch crossings in layers 3–6. Preliminary results of studies that take advantage of recent progress in the track reconstruction show a relative reduction of the b-tag purity at 80 % efficiency by about 7 %, while the c-tag purity at 60 % efficiency suffers a relative loss of 13 % [26]. Reducing the impact using advanced pattern recognition techniques in a detector design with smaller pixels that integrate over a whole bunch train is the subject of ongoing efforts.

The distributions of energy versus polar angle for the dominant beam-induced backgrounds. The contributions from incoherent \({{\mathrm{e}}}^{+}{{\mathrm{e}}}^{-}\) production (a) and from \({{\gamma }{\gamma }\rightarrow \mathrm {hadrons}} \) processes (b) are shown separately. The number of entries in the histograms corresponds to the \(4.5\times 10^{5}\) (a) and 4.1 (b) events expected for one bunch crossing. The number of events in each bin in \((\theta , E)\) is shown in a color scale

Full size image

While the most energetic particles from incoherent pair production are primarily outside of the detector acceptance of both detectors, some low-\(p_\mathrm {T}\) particles lead to an occupancy of up to 0.06 \(\text {hits} / \mathrm {mm}^{2}\) per bunch crossing in the vertex detector and up to \(5\times 10^{-5}\) \(\text {hits}/\mathrm {mm}^{2}\) per bunch crossing in the main tracker for the SiD detector model. They do not, however, impact on the energy reconstruction. Particles from \({{\gamma }{\gamma }\rightarrow \mathrm {hadrons}} \) processes on the other hand can have sizable values of \(p_\mathrm {T}\) and reach the calorimeters, affecting the jet energy resolution. The beam-induced backgrounds do not degrade the tracking performance significantly [11].

The primary vertices of the beam-induced backgrounds are distributed with a Gaussian profile along the beam direction across the luminous region of 225 \(\upmu \mathrm {m}\), taking into account the bunch length along the beam direction.

6.1 Reconstruction of isolated leptons

Signal events with six jets contain one high-energy isolated lepton from the leptonic \({\mathrm{W}}\) boson decay. No isolated leptons are expected in signal events with two hadronic \({\mathrm{W}}\) decays. Hence the number of isolated leptons is an important observable in the signal selections for both final states.

The electron and muon identification criteria used in this study are based on the energy deposition in the ECAL and HCAL and the momentum measured by the tracker. Electron candidates are selected by requiring that almost all of the energy deposition is in the ECAL and that the total calorimetric energy deposition is consistent with the momentum measured by the tracker. For the muon candidates, most of the energy deposition is in the HCAL, while the calorimetric energy is required to be small compared to the corresponding momentum measured by the tracker. A selection on the impact parameter reduces non-prompt leptons.

The SiD analysis uses the IsolatedLeptonFinder processor implemented in MarlinReco [23] to identify leptons in regions with otherwise little calorimetric activity. The ILD analysis additionally exploits the transverse distance from the jet axis to identify leptons from leptonic \({\mathrm{W}}\) decays.

The electron and muon identification capabilities of the reconstruction in a multi-jet environment were tested in a sample of four jets, one lepton and missing energy. The efficiency is defined as the fraction of leptons with correctly identified flavor. The purity is defined as the ratio of the number of leptonic \({\mathrm{W}}\) decays of a given flavor to the number of reconstructed isolated leptons of that flavor. Leptons from heavy flavor meson decays and charged pions mis-identified as leptons are considered as background. An efficiency of 82 % (89 %) and purity of 95 % (97 %) for electrons (muons) is observed in ILD and 86 % (86 %) efficiency and 94 % (95 %) purity for electrons (muons) in SiD.

6.2 Suppression of beam-induced backgrounds

The particles from beam-induced backgrounds as described in Sect. 5 tend towards low transverse momenta and small angles with respect to the beam axis. Different approaches are used to suppress the impact of the beam-induced backgrounds. For the SiD analysis, only the reconstructed objects in the range \(20^{\circ } < \theta < 160^{\circ }\) are considered, because the \(\mathrm{t}\bar{\mathrm{t}}\mathrm{H}\) final state is produced via s-channel exchange and is not suppressed by this selection. In the ILD analysis, the longitudinally-invariant \(k_T\) jet algorithm [27, 28] with a value of 1.2 for the R parameter is employed to suppress the particles close to the beam axis. Only the particles grouped into the physics jets by the \(k_T\) algorithm are considered further in the analysis. Figure 3 shows how the impact of the beam-induced backgrounds on the reconstructed Higgs mass is mitigated by the removal procedure. A modified version of the Durham jet finding algorithm [29] then groups all particles in the event into a specified number of jets, without splitting decay products of secondary vertices across different jets.

Distribution of the reconstructed Higgs mass in the ILD detector for signal events (six jets + lepton mode), without beam-induced backgrounds using the Durham algorithm, and with background, comparing the Durham and \(k_T\) algorithms. By itself, the \(k_T\) algorithm performs better than the Durham algorithm in terms of mitigating the effects of beam-induced backgrounds

Full size image

6.3 Jet clustering and flavor identification

Depending on the signal definition for the semi-leptonic or hadronic final state, the Durham jet clustering algorithm is used in the exclusive mode to cluster the event into six or eight jets, respectively. In either case, the isolated leptons described in Sect. 6.1 are removed before the jet finding steps.

Heavy flavor identification is primarily used to remove the \({\mathrm{t}}\bar{{\mathrm{t}}} \) background. Both the six-jet and eight-jet final states contain four \({\mathrm{b}}\) jets. The \({\mathrm{t}}\bar{{\mathrm{t}}} \) events contain no more than two \({\mathrm{b}}\) jets from the top decays as do \(\sim \)80 % of \({\mathrm{t}}\bar{{\mathrm{t}}}{\mathrm{Z}} \). The flavor tagging classifier for the measurement of \(\mathrm{t}\bar{\mathrm{t}}\mathrm{H}\) production was trained on events with six quarks of the same flavor produced in electron-positron annihilation. For the training, 60000 cand \({\mathrm{b}}\) jets, and 180000 light flavor quark jets are used. These samples were chosen since the jets have similar kinematic properties as those in \(\mathrm{t}\bar{\mathrm{t}}\mathrm{H}\) signal events. For a b-jet identification efficiency of 50 %, the misidentification fraction is found to be about 0.12 % for c-jets and about 0.05 % for light quark (uds) jets for both detectors, evaluated using two-jet final states at \(\sqrt{s}=91~\mathrm{GeV}\) [11]. The misidentification fraction for c-jets (uds-jets) increases by 100 % (60 %) when incoherent electron-positron pairs and \({{\gamma }{\gamma }\rightarrow \mathrm {hadrons}} \) processes are included in the simulation as described in Sect. 5. For the six-quark final states at \(\sqrt{s}=1~\mathrm{TeV}\), the misidentification fractions typically degrade by a factor of two in addition due to the varying jet energies and the confusion in the jet clustering due to the increased number of jets.

Figure 4a shows the distribution of the response from the flavor-tagging multivariate selection for the jet that has the third-highest tagging probability. In both analyses, the shape of the distribution of the flavor tagging response, rather than a simple cut, is used. The background channels, in particular \({\mathrm{t}}\bar{{\mathrm{t}}} \), are dominated by the peak at low values. The peak at higher values in the \({\mathrm{t}}\bar{{\mathrm{t}}}{\mathrm{Z}} \) channel is due to events with four genuine \({\mathrm{b}}\) jets.

Distributions of the event selection variables for the different signal and background processes in the ILD detector: a the third highest b-tag in the event; b the event thrust; c the jet resolution parameter \(Y_{65}\); d the number of reconstructed particles in the event. All histograms are normalized to unit area

Full size image

6.4 Reconstruction of W, top and Higgs candidates

To form \({\mathrm{W}}\), top and Higgs candidates from the reconstructed jets, the following function is minimized for the final state with eight jets:

where \(M_{12}\) and \(M_{45}\) are the invariant masses of the jet pairs used to reconstructed the \({\mathrm{W}}\) candidates, \(M_{123}\) and \(M_{456}\) are the invariant masses of the three jets used to reconstruct the top candidates and \(M_{78}\) is the invariant mass of the jet pair used to reconstruct the Higgs candidate. \(M_{{\mathrm{W}}}\), \(M_{{\mathrm{t}}}\) and \(M_{{\mathrm{H}}{}{}}\) are the nominal \({\mathrm{W}}\), top and Higgs masses. The resolutions \(\sigma _{{\mathrm{W}}}\), \(\sigma _{{\mathrm{t}}}\) and \(\sigma _{{\mathrm{H}}{}{}}\) were obtained from reconstructed jet combinations matched to \({\mathrm{W}}\), top and Higgs particles at generator level. The corresponding function minimized for the six-jets final state is given by:

In the ILD analysis, the \({\mathrm{b}}\) tagging information is also used to reduce the number of combinations by forming the Higgs candidate from two of the four jets with the highest value of the \({\mathrm{b}}\)-tagging classifier. The other jets in the event are used to form the hadronic top candidates.

Events were selected using boosted decision trees (BDTs) as implemented in TMVA [30]. The BDTs were trained separately for the eight- and six-jets final states. The following input variables were used:

• the four highest values of the \({\mathrm{b}}\)-tagging classifier. The third (see Fig. 4a) and fourth highest \({\mathrm{b}}\)-tag values are especially suited to reject \({\mathrm{t}}\bar{{\mathrm{t}}} \) and most of the \({\mathrm{t}}\bar{{\mathrm{t}}}{\mathrm{Z}} \) events, both of which contain only two \({\mathrm{b}}\) jets;

• the event thrust [31] (see Fig. 4b) defined as

  $$\begin{aligned} T = \max \frac{\sum _{i}|\hat{n}\cdot \mathbf {p_{i}}|}{\sum _{i}|\mathbf {p_{i}}|}, \end{aligned}$$

  (3)

  where \(p_{i}\) is the momentum of the jet. Since the top quarks in \({\mathrm{t}}\bar{{\mathrm{t}}} \) events are produced back-to-back, the thrust variable has larger values in \({\mathrm{t}}\bar{{\mathrm{t}}} \) events compared to \(\mathrm{t}\bar{\mathrm{t}}\mathrm{H}\), \({\mathrm{t}}\bar{{\mathrm{t}}} Z\) or \({\mathrm{t}}\bar{{\mathrm{t}}} {{\mathrm{b}}\bar{{\mathrm{b}}}} \) events;

• the jet resolution parameter from the Durham algorithm in the E recombination scheme \(Y_{ij}\), when combining i jets to \(j = (i - 1)\) jets. For the six-jets final state \(Y_{54}\) and \(Y_{65}\) (see Fig. 4c) are found to be effective, while \(Y_{76}\) and \(Y_{87}\) are used for the eight-jets final state. Isolated leptons are removed prior to the jet clustering;

• the number of identified isolated electrons and muons (ILD only);

• the missing transverse momentum, \(p_\mathrm {T}^{\text {miss}}\). Due to the leptonic \({\mathrm{W}}\) boson decay, finite values of \(p_\mathrm {T}^{\text {miss}}\) are reconstructed for six-jets signal events while \(p_\mathrm {T}^{\text {miss}}\) tends towards zero for eight-jets signal events;

• the visible energy of the event defined as the scalar sum of all jet energies;

• the masses \(M_{12}\), \(M_{123}\) and \(M_{45}\) as defined in Sect. 6.4.

For the eight-jets final state additionally the two variables \(M_{456}\) and \(M_{78}\) as defined in Sect. 6.4 are included.

The ILD analysis includes the helicity angle of the Higgs candidate as defined by the angle between the two \({\mathrm{b}}\) jet momenta in the dijet rest frame.

To select events, cuts on the BDT response are applied. The cuts were optimized by maximizing the signal significance given by: \(\frac{S}{\sqrt{S+B}}\), where S is the number of signal events and B is the number of background events. As an example, the reconstructed top and Higgs masses in six-jets events after the cut on the BDT output are shown in Fig. 5. The selection efficiencies (purities) for signal events are 33.1 % (27.7 %) and 56.0 % (25.2 %) for the six- and eight-jets analyses in ILD, respectively, and 30.5 % (28.9 %) and 45.9 % (26.7 %) in SiD. In Table 2 the expected yields are shown separately for all investigated final states.

Reconstructed top (a) and Higgs (b) masses for selected six-jets events in the SiD detector. All histograms were normalized to an integrated luminosity of 1 \(\mathrm {ab}^{-1} \). The distribution for \({\mathrm{t}}\bar{{\mathrm{t}}} \) was scaled by a factor of 0.5

Full size image

The cross section can be directly obtained from the number of background-subtracted signal events after the selection. The uncertainty of the cross section measurement was estimated using the number of selected signal and background events. Assuming an integrated luminosity of 1 \(\mathrm {ab}^{-1} \) split equally between the \(P({{\mathrm{e}}}^{-},{{\mathrm{e}}}^{+}) = (-80~\%,+20~\%)\) and \(P({{\mathrm{e}}}^{-},{{\mathrm{e}}}^{+}) = (+80~\%,-20~\%)\) beam polarization configurations, the cross section can be measured with a statistical precision of 10–11 % using the eight-jets final state and with a statistical precision of \(\approx \)13 % for the six-jets final state.

The uncertainties of the measured cross sections translate to precisions on the top Yukawa coupling of 5–6 and \(\approx \)7 % from the eight- and six-jets final states, respectively. If both measurements are combined, the top Yukawa coupling can be extracted with a statistical precision of better than 4.5 %.

For 1 \(\mathrm {ab}^{-1} \) of data with only \(P({{\mathrm{e}}}^{-},{{\mathrm{e}}}^{+}) = (-80~\%,+20~\%)\) polarization, this number improves to \(4~\%\). The precision for the six-jets final state could be improved further if \(\tau \) leptons were included in the reconstruction. Additional improvement is also foreseen by employing kinematic fitting. The achieved precision of both analyses indicates that the reconstruction of the investigated final states is not significantly affected by the differences in subdetector performance [11] between the two investigated concepts. This is consistent with the findings of the study of the Higgs self-coupling in the six-jets final state of the \({\mathrm{Z}}{\mathrm{H}}{}{}{\mathrm{H}}{}{}\) channel [32, Chapter 2.5.2], where the confusion in the jet clustering was a dominant contribution to the mass resolution.

Given the low cross section and relatively clean environment at a \(\sqrt{s}=1~\mathrm{TeV}\) ILC, it is expected that the statistical uncertainty of the measurement of the top Yukawa coupling in direct observation dominates over systematic uncertainty. In the following we estimate the contributions from the main sources of systematic errors to this measurement.

The number of background events in the final selection is comparable to the number of signal events, making the estimation of normalization and shapes of the background an important source of systematic uncertainty. The total cross sections calculable from theory for the \({\mathrm{t}}\bar{{\mathrm{t}}}{\mathrm{Z}} \) and \({\mathrm{t}}\bar{{\mathrm{t}}} \) processes are expected to improve in the coming years. For the \({\mathrm{t}}\bar{{\mathrm{t}}}{\mathrm{Z}} \) process, the QCD and electroweak corrections are known at the 1-loop level [33]. For the \({\mathrm{t}}\bar{{\mathrm{t}}} \) process, the electroweak corrections are known to the 1-loop level [34], while the QCD corrections are known at the 3-loop level [35–39]. QCD contributions to the \({\mathrm{t}}\bar{{\mathrm{t}}} {\mathrm{b}}\bar{{\mathrm{b}}} \) cross section make this value more challenging to compute precisely; in principle the measurement of the gluon splitting rate at relevant energies will provide a handle to estimate its size. A crucial aspect in the estimation of the efficiencies is the accurate modeling of the event selection variables. Here we illustrate how one might arrive at control samples for different background sources in order to estimate the efficiency of each component accurately.

The \({\mathrm{t}}\bar{{\mathrm{t}}}{\mathrm{Z}} \) final state can be reconstructed in a similar fashion to the \(\mathrm{t}\bar{\mathrm{t}}\mathrm{H}\) final state. For hadronically decaying \({\mathrm{Z}}\), the number of jets in the final state will be the same as in the \(\mathrm{t}\bar{\mathrm{t}}\mathrm{H}\) analysis. For our nominal integrated luminosities of \(0.5 ~\mathrm {ab}^{-1} \) for each of the two polarization states, nearly 1300 signal events \(\mathrm{t}\bar{\mathrm{t}}\mathrm{H}(\rightarrow {\mathrm{b}}\bar{{\mathrm{b}}})\) are expected and 800 events are expected for \({\mathrm{t}}\bar{{\mathrm{t}}}{\mathrm{Z}} (\rightarrow {\mathrm{b}}\bar{{\mathrm{b}}})\), taking into account the \({\mathrm{Z}}\rightarrow {\mathrm{b}}\bar{{\mathrm{b}}} \) branching ratio. Other hadronic decays of the \({\mathrm{Z}}\) boson will have large \({\mathrm{t}}\bar{{\mathrm{t}}} \) background due to the absence of the two \({\mathrm{b}}\) tags. Including leptonic decays of the \({\mathrm{Z}}\) boson will help increase the sensitivity to this channel. Overall, one can expect that the statistical uncertainty for \({\mathrm{t}}\bar{{\mathrm{t}}}{\mathrm{Z}} \) will be similar to that of \(\mathrm{t}\bar{\mathrm{t}}\mathrm{H}\), i.e. at the few percent level.

The large cross section of \({\mathrm{t}}\bar{{\mathrm{t}}} \) events will allow for detailed systematic studies. While only a certain class of these events may enter the final selection, we estimate that the systematic uncertainty to the measurement of the top Yukawa coupling is comparable to that of \({\mathrm{t}}\bar{{\mathrm{t}}}{\mathrm{Z}} \).

Other sources of systematic uncertainty such as the luminosity measurement, jet energy scale, and flavor tagging are typically at the 1 % level or better for \({{\mathrm{e}}}^{+}{{\mathrm{e}}}^{-} \) colliders [40]. The uncertainty on BR(\({\mathrm{H}}{}{} \rightarrow {{\mathrm{b}}\bar{{\mathrm{b}}}} \)) is not taken into account in our calculation of the top Yukawa coupling from the \(\mathrm{t}\bar{\mathrm{t}}\mathrm{H}\) production cross section. It is expected that this quantity can be measured with a precision of better than \(1~\%\) using \(\mathrm{e}^{+}\mathrm{e}^{-} \rightarrow {\nu }\bar{\nu }{\mathrm{H}}{}{} \) events [41, 42].

The physics potential for a measurement of the top Yukawa coupling at 1 TeV at the ILC is investigated. The study is based on detailed detector simulations using both the SiD and ILD detector concepts. Beam-induced backgrounds are considered in the analysis. The combination of results obtained for two different final states leads to a statistical uncertainty on the top Yukawa coupling of better than 4.5 % for an integrated luminosity of 0.5 \(\mathrm {ab}^{-1} \) with the \(P({{\mathrm{e}}}^{-},{{\mathrm{e}}}^{+}) = (-80~\%,+20~\%)\) beam polarization configuration and 0.5 \(\mathrm {ab}^{-1} \) with \(P({{\mathrm{e}}}^{-},{{\mathrm{e}}}^{+}) = (+80~\%,-20~\%)\) polarization. If 1 \(\mathrm {ab}^{-1} \) of data were recorded with only the \(P({{\mathrm{e}}}^{-},{{\mathrm{e}}}^{+}) = (-80~\%,+20~\%)\) beam polarization configuration, the expected precision would improve to \(4~\%\).

The results from the studies presented in this paper demonstrate the robustness of the physics reconstruction of high jet multiplicity final states at \(\sqrt{s}=1\) TeV under realistic simulation conditions. The expected precisions for measurements of the top Yukawa coupling were found to be very similar for two different detector concepts.