Study of single top production at high energy electron positron colliders

The effect of single top production on the study of top quark pair production in future high energy electron–positron colliders is evaluated. The rate of the single top quark production process is sizeable throughout a large range of center-of-mass energies and the final state cannot easily be distinguished from the dominant pair production process. We discuss the impact on the top quark mass extraction from a scan through the pair production threshold and the determination of top quark form factors in the continuum. These results advocate for the exploration of the inclusive e+e-→W+bW-b¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^+e^-\rightarrow W^+bW^-\bar{b}$$\end{document} process, that includes both top quark pair and single top quark production.


Introduction
A high-luminosity, high-energy, linear e + e − collider yields excellent opportunities for precision tests of the Standard Model of particle physics. The combination of precisely calculable electroweak production and strict control of the initial state with the relatively benign experimental environment and state-of-the-art detector systems allow for a characterization of Standard Model and new physics processes with a precision that goes well beyond what can be achieved at hadron colliders.
Two projects of linear electron-positron colliders are being considered: the International Linear Collider (ILC [1,2]) and the Compact Linear Collider (CLIC [3]). The physics case for a linear e + e − machine has been made in great detail in Refs. [4][5][6][7][8][9][10][11]. The specific case of a multi-TeV e + e − collider is discussed in Refs. [12][13][14]. In both cases, the center-of-mass energy will exceed √ s = 350 GeV, the threshold for top quark pair production. Unlike other quarks, the top quark has never been produced in e + e − machines, and therefore a precise measurement of electroweak top quark pair production is missing. The study of top quark properties is therefore one a e-mail: marcel.vos@ific.uv.es of the most exciting prospects for a future linear collider [15]. Detailed full-simulation studies have been made of the prospects for a precise top quark mass measurement [16,17] and characterization of the tt Z and ttγ vertices [18].
Single top production, through the e + e − → W − tb, W +t b process depicted in the central panel of Fig. 1, is abundant at e + e − colliders that operate at √ s > 300 GeV. Note that for the t → bW decay, this process gives rise to the same W + bW −b final state as top pair production. Using MADGRAPH we find ten leading-order diagrams of this type, against two for top quark production. A third group of processes gives rise to the same final state: W W Z, W W h and W W γ production, with Z /γ / h → bb. At treelevel there are approximately 50 diagrams that produce two W -bosons and two b-quarks that do not involve top quarks.
Ultimately, all three processes yield the same set of six-fermion final states. In many studies of the linear collider prospects for top physics, single top quark production has been neglected. Notable exceptions are found in Refs. [19,20].
A fully consistent analysis of the inclusive e + e − → W + bW −b is currently impossible. Even if event generation for the ILC TDR routinely included the full 2 → 6 matrix element at leading order, higher-order corrections for this process are not available for the extraction of top quark properties and couplings. In the continuum, results for the cross-section of top quark pair production are available at NNLO. Cross-sections for the process W + bW −b are for the moment only available at LO, but are required at least at NLO (note that NLO calculations are already available for the LHC).
In this note we investigate the impact of single top events in the study of top quark pair production and propose a more inclusive experimental strategy that compares precise predictions for the e + e − → W + bW −b process to measurements. As top quark pair production and single top quark production give rise to the same six-fermion final state, the question arises as to how one can distinguish both sources. At a fundamental level the single top and top quark pair production processes are entangled by interference between the different diagrams. No algorithm can ever separate them fully. However, one could hope to use some of the marked features of the e + e − → tt process to make it stand out among the other processes that give rise to the W + bW −b final state. One could then hope to isolate samples that are enriched in top quark pairs or single top quark events.
In Fig. 2 we present the invariant mass of the W + b and W −b combinations at truth level using e + e − → tt → W + bW −b events generated with WHIZARD [21,22] at √ s = 500 GeV, including the effect of initial state radiation and the beam energy spread expected at the ILC. The majority of events is found around However, in a significant fraction of events one of the Wboson b-quark pairs has an invariant far from the top quark mass. The cross is populated by a mixture of tt events with an off-shell top quark and genuine single-top production through diagrams such as that in the central panel of Fig. 1 In Fig. 3 we present the reconstructed beam energy of the W b decaying leptonically versus the W b decaying hadronically. The energy of the W b pair tends to cluster at about half the center-of-mass energy, as expected for double-top events. We observe however a significant fraction of events along the diagonal E lep + E had = √ s. Those are mostly single top events. This figure suggests another potential cri-terium for partial separation of single and double-top events using the reconstructed beam energy.
These figures suggest an (ad-hoc) truth-level categorisation of events according to the number of on-shell top quark candidates, that is used to quantify the non-tt contribution to e + e − → W + bW −b production in the remainder of this paper. We consider the event to correspond to top quark pair production whenever the following relation is satisfied for both W b combinations: Events that meet this criterium only for one of the W b pairs are labeled as single top quark events, and the remaining events are considered as non-top events.
With this criterion we find that at 500 GeV typically 90 % of e + e − → W + bW −b events is formed by top quark pairs, 9 % corresponds to single top production and only 1 % stems from processes that do not involve top quarks. These fractions depend strongly on the center-of-mass energy of the collider (more in Sect. 5) and on the beam polarization (more in Sect. 3). Over a broad range of energies and operating scenarios the fraction of single-tops is far from negligible and may have a significant impact on the measurement of top quark properties and the search for signs of new physics in tt production at lepton colliders.

Experimental study at √ s = 500 GeV
In this section we concentrate on the impact of realistic analysis cuts on the composition of the e + e − → W + bW −b sample at √ s = 500 GeV. The analysis is based on the study of cross-sections and asymmetries at √ s = 500 GeV in Ref. [18]. Final states of the type lνq 1 q 2 b 1 b 2 are generated using the WHIZARD Monte Carlo program. These events are mostly W + bW −b events, where one W decays hadronically and the other leptonically. Some 100,000 MC events are processed in a detailed simulation of the ILD detector [23,24] response based on GEANT4 [25]. Events are reconstructed, to the level of particle flow objects and jets, using the ILC software package. In these samples the beams are 100 % polarized.
First, we discuss briefly the single top content as a function of beam polarization. The fraction of W + bW −b events that corresponds to single top quark production is quite sensitive to the polarization of the electron and positron beams. We evaluate the fractions of top quark pair, single top, and non-top events using the criterion of Eq. (2.1) on e + e − → elνq 1 q 2 b 1 b 2 samples generated with WHIZARD, including initial state radiation and a realistic ILC luminosity profile. For a fully left-handed electron beam and fully right-handed positron beam (e − L e + R ) the total cross-section is nearly three times larger. The top quark pair fraction remains 90.2 %. The single top and non-top fractions are 8.9 % and 0.9 %, respectively. For the opposite (e − R e + L ) configuration, the total cross-section is similar to the unpolarized result. The top quark pairs make up 94.3 % of the sample, with only 5.6 % of single top quarks and about 0.1 % for non-top production.
The selection in Ref. [18] rejects background and improves the fraction of well-reconstructed events using a cut on a χ 2 formed by the reconstructed top mass, beam energy and bquark energy in the top rest frame. Only the hadronically decaying top is used in the analysis. Table 1 summarizes the selection efficiencies at various stages of the analysis. The preselection includes the lepton reconstruction and identification, and the requirement of two b-tagged jets. The results in the second row are obtained after the cut on the hadronic top χ 2 < 15 employed in Ref. [18]. We observe that for both beam polarizations the efficiency for single top events is significantly smaller than for top quark pairs. This is expected, since the W b pair mass and the beam energy are required to be compatible with the values expected for top quark pair production in the χ 2 cut.
The inclusion of cuts on the leptonic top in the χ 2 does not improve the rejection. The third row in Table 1 corresponds to cuts on both the hadronic and leptonic top χ 2 , where the χ 2 cut is adjusted to obtain the same efficiency on the total sample. This can be understood as follows: once an event passes the beam energy constraint, momentum conservation (forced by the neutrino reconstruction) imposes a similar invariant mass on both the hadronic and leptonic sides of the event, whether the event is a single or a double top. In other words, the leptonic side does not provide any additional background rejection.  The selection described here is insufficient to reduce the single top contamination to a safe level. The fraction of single top events continues to be large (5 % in the present example). At this level it is expected to have a non-negligible effect on the cross-section measurement. An attempt to describe the remaining contamination using Monte Carlo events is likely to bring sizeable systematic uncertainties.
The forward-backward asymmetry is even more sensitive, since the asymmetry of single top production is very small compared to that of top quark pair production. For the e − L e + R beam polarization the asymmetry of the W b system in inclusive W + bW −b production process is 6 % smaller than the forward-backward asymmetry in top quark pair production, a relative difference of 14 %. We conclude, therefore, that cuts on the properties of the reconstructed top candidates are unlikely to reduce single top quark production to the level where it has a negligible impact on observables such as the cross-section and forwardbackward asymmetry.

Analysis of top mass at threshold
As pointed out a long time ago, the cross-section of the e + e − → tt at threshold could be used to measure in a precise and well defined way the top mass [26]. Detailed calculations are present in Ref. [27] and a complete study of the extraction of top quark properties was performed in Ref. [28].
Recent reanalyses include realistic beam energy spectra for the ILC and CLIC [16] and beam polarization [17]. In the following we use the former analysis as a reference. The selection of events includes a kinematic fit, but as we showed in the previous section, the leptonic side of the event brings no additional rejection of single top events, since the neutrino can always be adjusted to fake a second top. This analysis is very detailed, but single top events have not been included. This paper includes for the case of ILC the program TOP-PIK, a NNLO calculation of the cross-section including the 1S resonance, initial state radiation and the ILC luminosity  [29] determines the 4-loop correction for the conversion from 1S to theM S scheme; the related theory systematic could be as small as 10 MeV, thus prompting a thorough evaluation of the experimental systematic uncertainties).
NLO calculations for the W + bW −b process have recently been implemented in WHIZARD around the double-top production threshold including QCD bound state effects. 1 We use these to estimate the single top fraction in W bW b production. In Table 2 and Fig. 4 the fraction of single top events is given for several center-of-mass energies around the pair production thresholds. The results include the effect of initial state radiation, of a realistic beam energy spectrum for the ILC at this energy, and the effect of polarized beams. The energy range and top quark mass follow the choices in Ref. [16].
We note that the presence of single top events modifies significantly the cross-section. Moreover, the fraction of single top events varies rapidly across the threshold region. This does not mean, however, that the top mass measurement is affected. According to the method proposed in [16], a fit is performed to the cross-section measurement, using a predefined function obtained from the double-top calculation at NNLO. As long as the shape of the curve is not modified, the result is not affected. The W + bW −b and tt curves are displayed in Fig. 4. We observe that according to the NLO calculation of the W + bW −b cross-section, the tt crosssection may be obtained quite exactly by shifting down the W + bW −b cross-section by some 0.022 pb. After fixing this offset with the measurement at 344 GeV we fit the W + bW −b shape with the curve for tt production. We find the best-fit mass is shifted by 30 MeV with respect to the input mass. This bias is of a similar size as the statistical uncertainty. We emphasize that this estimate should be considered as an upper limit. A more sophisticated treatment is likely to reduce the uncertainty due to the modelling of the single top contamination further. Ultimately, this uncertainty is absent in a fit of the prediction for the inclusive e + e − → W + bW −b process to the data.
This result is also supported by the NNLL result in Ref. [30], that we reproduce in Fig. 5. The calculations are not directly comparable, as initial state radiation and beam energy spread are not included in Ref. [30]. The authors find that the impact of a cut on the W b mass is an approximately constant shift. The cross-section for tt events with W b massin  (1) non-res Fig. 6 The effect on the cross-section of non-resonant corrections, according to Beneke et al. [31], for e + e − annihilation to W + bW −b in the tt threshold region. The cross-section for unpolarized beams is given as a function of center-of-mass energy. With respect to the original publication, the curve is shifted such that the 1S peak coincides approximately with those of Figs. 4 and 5 a range −15 < m W b − m t < 15 is, as before, shifted down from the inclusive cross-section by 0.050 pb. Finally, Fig. 6 shows the effect of non-resonant effects on the W bW b crosssection. Again, the shift in the cross section is essentially constant across the threshold region.
The single top content of the samples depends strongly on the energy, as seen in Table 2. In analyses where the top quark is reconstructed and cuts are applied on the mass and/or energy of the candidates, the selection efficiency can be quite different for top quark pair and single top events (up to 50 % smaller for single top in the example above). This means that, after correcting for efficiency, the shape of the curve will be deformed if single top events are ignored. Note that this does not apply to Ref. [16], since in this analysis no selection cuts based on reconstructed top quark mass and energy are applied.
As the most precise calculations are only available for the e + e − → tt process, none of the studies of future collider prospects so far has taken into account the full W + bW −b process. As soon as the necessary tools have been developed we suggest to repeat the study of the top quark mass extraction of future lepton colliders using W + bW −b events and polarized beams.

Energy dependence
The composition of the W + bW −b sample is energy dependent as shown in Fig. 7, that presents the unpolarized LO cross-sections for e + e − → W + bW −b . The solid line corresponds to the full 2 → 4 process, including single top  Fig. 7 Cross-section for W + bW −b and tt events produced in unpolarized electron-positron annihilations as a function of the center-of-mass energy. The curves are obtained using WHIZARD at LO production and non-top production. The dashed line presents the results for the e + e − → tt → W + bW −b process. All results have been obtained using WHIZARD. The difference between the tt cross-section and the full 2 → 4 result is 10 % at 500 GeV and rises to 20 % at 1 TeV and nearly 50 % at 3 TeV.
The relatively small contribution of the e + e − → tt process at large center-of-mass energy is confirmed by an analysis of the number of on-shell top quarks (|m W b − m MC t | < 15 GeV) in W + bW −b events at 3 TeV. The fraction of events with two on-shell top quarks is 48 %, events with a single on-shell top quark make up 35 % of the sample and the remaining 17 % of events has no on-shell top quarks. We conclude that while the rate for the e + e − → tt process drops at very large center-of-mass energy, single-top and nontop production increase rapidly. The three processes become comparable in size for e + e − colliders operating in the multi-TeV regime.

Summary and conclusions
The prospects for precision top quark physics at future lepton colliders relies on a comparison of measured cross-sections to very precise predictions of top quark pair production, including NLO and sometimes NNLO calculations. However, as we have shown in this note, single top events, that are ignored in most prospect studies, may lead to significant effects. A precise experimental separation of single top and top quark pair production processes seems difficult. For this reason we advocate the analysis of W + bW −b production, that includes single top quark production, and leads to increased statistics and possibly smaller systematic errors. To achieve the ultimate possible precision at the ILC, observables for the W + bW −b final state must be calculated with a precision well below 1 %, at NNLO, both at production threshold and in the continuum.