PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of springeropenLink to Publisher's site
The European Physical Journal. C, Particles and Fields
 
Eur Phys J C Part Fields. 2017; 77(7): 475.
Published online 2017 July 17. doi:  10.1140/epjc/s10052-017-4968-5
PMCID: PMC5587080

Higgs physics at the CLIC electron–positron linear collider

H. Abramowicz,1 A. Abusleme,2 K. Afanaciev,3 N. Alipour Tehrani,4 C. Balázs,5 Y. Benhammou,1 M. Benoit,6 B. Bilki,7 J.-J. Blaising,8 M. J. Boland,9 M. Boronat,10 O. Borysov,1 I. Božović-Jelisavčić,11 M. Buckland,12 S. Bugiel,13 P. N. Burrows,14 T. K. Charles,5 W. Daniluk,15 D. Dannheim,4 R. Dasgupta,13 M. Demarteau,7 M. A. Díaz Gutierrez,2 G. Eigen,16 K. Elsener,4 U. Felzmann,9 M. Firlej,13 E. Firu,17 T. Fiutowski,13 J. Fuster,10 M. Gabriel,18 F. Gaede,4,19 I. García,10 V. Ghenescu,17 J. Goldstein,20 S. Green,21 C. Grefe,4 M. Hauschild,4 C. Hawkes,23 D. Hynds,4 M. Idzik,13 G. Kačarević,11 J. Kalinowski,24 S. Kananov,1 W. Klempt,4 M. Kopec,13 M. Krawczyk,24 B. Krupa,15 M. Kucharczyk,15 S. Kulis,4 T. Laštovička,25 T. Lesiak,15 A. Levy,1 I. Levy,1 L. Linssen,4 S. Lukić,11 A. A. Maier,4 V. Makarenko,3 J. S. Marshall,21 V. J. Martin,22 K. Mei,21 G. Milutinović-Dumbelović,11 J. Moroń,13 A. Moszczyński,15 D. Moya,26 R. M. Münker,4 A. Münnich,4 A. T. Neagu,17 N. Nikiforou,4 K. Nikolopoulos,23 A. Nürnberg,4 M. Pandurović,11 B. Pawlik,15 E. Perez Codina,4 I. Peric,27 M. Petric,4 F. Pitters,4 S. G. Poss,4 T. Preda,17 D. Protopopescu,28 R. Rassool,9 S. Redford,4 J. Repond,7 A. Robson,28 P. Roloff,corresponding author4 E. Ros,10 O. Rosenblat,1 A. Ruiz-Jimeno,26 A. Sailer,4 D. Schlatter,4 D. Schulte,4 N. Shumeiko,3 E. Sicking,4 F. Simon,18 R. Simoniello,4 P. Sopicki,15 S. Stapnes,4 R. Ström,4 J. Strube,4 K. P. Świentek,13 M. Szalay,18 M. Tesař,18 M. A. Thomson,21 J. Trenado,29 U. I. Uggerhøj,30 N. van der Kolk,18 E. van der Kraaij,16 M. Vicente Barreto Pinto,6 I. Vila,26 M. Vogel Gonzalez,2 M. Vos,10 J. Vossebeld,12 M. Watson,23 N. Watson,23 M. A. Weber,4 H. Weerts,7 J. D. Wells,31 L. Weuste,18 A. Winter,23 T. Wojtoń,15 L. Xia,7 B. Xu,21 A. F. Żarnecki,24 L. Zawiejski,15 and I.-S. Zgura17

Abstract

The Compact Linear Collider (CLIC) is an option for a future e+e- collider operating at centre-of-mass energies up to 3 TeV, providing sensitivity to a wide range of new physics phenomena and precision physics measurements at the energy frontier. This paper is the first comprehensive presentation of the Higgs physics reach of CLIC operating at three energy stages: s=350GeV, 1.4 and 3 TeV. The initial stage of operation allows the study of Higgs boson production in Higgsstrahlung (e+e- → ZH) and WW-fusion (e+e-Hνeν¯e), resulting in precise measurements of the production cross sections, the Higgs total decay width ΓH, and model-independent determinations of the Higgs couplings. Operation at s>1TeV provides high-statistics samples of Higgs bosons produced through WW-fusion, enabling tight constraints on the Higgs boson couplings. Studies of the rarer processes e+e-tt¯H and e+e-HHνeν¯e allow measurements of the top Yukawa coupling and the Higgs boson self-coupling. This paper presents detailed studies of the precision achievable with Higgs measurements at CLIC and describes the interpretation of these measurements in a global fit.

Introduction

The discovery of a Higgs boson [1, 2] at the Large Hadron Collider (LHC) provided confirmation of the electroweak symmetry breaking mechanism [38] of the Standard Model (SM). However, it is not yet known if the observed Higgs boson is the fundamental scalar of the SM or is either a more complex object or part of an extended Higgs sector. Precise studies of the properties of the Higgs boson at the LHC and future colliders are essential to understand its true nature.

The Compact Linear Collider (CLIC) is a mature option for a future multi-TeV high-luminosity linear e+e- collider that is currently under development at CERN. It is based on a novel two-beam acceleration technique providing accelerating gradients of 100 MV/m. Recent implementation studies for CLIC have converged towards a staged approach. In this scheme, CLIC provides high-luminosity e+e- collisions at centre-of-mass energies from a few 100 GeV up to 3 TeV. The ability of CLIC to collide e+e- up to multi-TeV energy scales is unique. For the current study, the nominal centre-of-mass energy of the first energy stage is s=350GeV. At this centre-of-mass energy, the Higgsstrahlung and WW-fusion processes have significant cross sections, providing access to precise measurement of the absolute values of the Higgs boson couplings to both fermions and bosons. Another advantage of operating CLIC at s350GeV is that it enables a programme of precision top quark physics, including a scan of the tt¯ cross section close to the production threshold. In practice, the centre-of-mass energy of the second stage of CLIC operation will be motivated by both the machine design and results from the LHC. In this paper, it is assumed that the second CLIC energy stage has s=1.4TeV and that the ultimate CLIC centre-of-mass energy is 3 TeV. In addition to direct and indirect searches for Beyond the Standard Model (BSM) phenomena, these higher energy stages of operation provide a rich potential for Higgs physics beyond that accessible at lower energies, such as the direct measurement of the top Yukawa coupling and a direct probe of the Higgs potential through the measurement of the Higgs self-coupling. Furthermore, rare Higgs boson decays become accessible due to the higher integrated luminosities at higher energies and the increasing cross section for Higgs production in WW-fusion. The proposed staged approach spans around twenty years of running.

The following sections describe the experimental conditions at CLIC, an overview of Higgs production at CLIC, and the Monte Carlo samples, detector simulation, and event reconstruction used for the subsequent studies. Thereafter, Higgs production at s=350GeV, Higgs production in WW-fusion at s>1TeV, Higgs production in ZZ-fusion, the measurement of the top Yukawa coupling, double Higgs production, and measurements of the Higgs boson mass are presented. The paper concludes with a discussion of the measurement precisions on the Higgs couplings obtained in a combined fit to the expected CLIC results, and the systematic uncertainties associated with the measurements.

The detailed study of the CLIC potential for Higgs physics presented here supersedes earlier preliminary estimates [9]. The work is carried out by the CLIC Detector and Physics (CLICdp) collaboration.

Experimental environment at CLIC

The experimental environment at CLIC is characterised by challenging conditions imposed by the CLIC accelerator technology, by detector concepts optimised for the precise reconstruction of complex final states in the multi-TeV energy range, and by the operation in several energy stages to maximise the physics potential.

Accelerator and beam conditions

The CLIC accelerator design is based on a two-beam acceleration scheme. It uses a high-intensity drive beam to efficiently generate radio frequency (RF) power at 12 GHz. The RF power is used to accelerate the main particle beam that runs parallel to the drive beam. CLIC uses normal-conducting accelerator structures, operated at room temperature. These structures permit high accelerating gradients, while the short pulse duration discussed below limits ohmic losses to tolerable levels. The initial drive beams and the main electron/positron beams are generated in the central complex and are then injected at the ends of the two linac arms. The feasibility of the CLIC accelerator has been demonstrated through prototyping, simulations and large-scale tests, as described in the Conceptual Design Report [10]. In particular, the two-beam acceleration at gradients exceeding 100 MV/m has been demonstrated in the CLIC test facility, CTF3. High luminosities are achievable by very small beam emittances, which are generated in the injector complex and maintained during transport to the interaction point.

CLIC will be operated with a bunch train repetition rate of 50 Hz. Each bunch train consists of 312 individual bunches, with 0.5 ns between bunch crossings at the interaction point. The average number of hard e+e- interactions in a single bunch train crossing is much less than one. However, for CLIC operation at s>1TeV, the highly-focussed intense beams lead to significant beamstrahlung (radiation of photons from electrons/positrons in the electric field of the other beam). Beamstrahlung results in high rates of incoherent electron–positron pairs and low-Q2 t-channel multi-peripheral γγ → hadron events, where Q2 is the negative of the four-momentum squared of the virtual space-like photon. In addition, the energy loss through beamstrahlung generates a long lower-energy tail to the luminosity spectrum that extends well below the nominal centre-of-mass energy, as shown in Fig. Fig.1.1. Both the CLIC detector design and the event reconstruction techniques employed are optimised to mitigate the influence of these backgrounds, which are most severe at the higher CLIC energies; this is discussed further in Sect. 4.2.

Fig. 1
The luminosity spectrum for CLIC operating at s=3TeV, where s is the effective centre-of-mass energy after beamstrahlung and initial state radiation [11]

The baseline machine design allows for up to ±80% longitudinal electron spin-polarisation by using GaAs-type cathodes [10]; and provisions have been made to allow positron polarisation as an upgrade option. Most studies presented in this paper are performed for zero beam polarisation and are subsequently scaled to account for the increased cross sections with left-handed polarisation for the electron beam.

Detectors at CLIC

The detector concepts used for the CLIC physics studies, described here and elsewhere [11], are based on the SiD [12, 13] and ILD [13, 14] detector concepts for the International Linear Collider (ILC). They were initially adapted for the CLIC 3 TeV operation, which constitutes the most challenging environment for the detectors in view of the high beam-induced background levels. For most sub-detector systems, the 3 TeV detector design is suitable at all energy stages, the only exception being the inner tracking detectors and the vertex detector, where the lower backgrounds at s=350GeV enable detectors to be deployed with a smaller inner radius.

The key performance parameters of the CLIC detector concepts with respect to the Higgs programme are:

  • excellent track-momentum resolution of σpT/pT22·10-5 GeV-1, required for a precise reconstruction of leptonic Z decays in ZH events;
  • precise impact parameter resolution, defined by a ≲ 5 μm and b ≲ 15 μm GeV in σd02=a2+b2/(p2sin3θ) to provide accurate vertex reconstruction, enabling flavour tagging with clean b-, c- and light-quark jet separation;
  • jet-energy resolution σE/E ≲ 3.5% for light-quark jet energies in the range 100 GeV to 1 TeV, required for the reconstruction of hadronic Z decays in ZH events and the separation of Zqq¯ and Hqq¯ based on the reconstructed di-jet invariant mass;
  • detector coverage for electrons extending to very low angles with respect to the beam axis, to maximise background rejection for WW-fusion events.

The main design driver for the CLIC (and ILC) detector concepts is the required jet-energy resolution. As a result, the CLIC detector concepts [11], CLIC_SiD and CLIC_ILD, are based on fine-grained electromagnetic and hadronic calorimeters (ECAL and HCAL), optimised for particle-flow reconstruction techniques. In the particle-flow approach, the aim is to reconstruct the individual final-state particles within a jet using information from the tracking detectors combined with that from the highly granular calorimeters [1518]. In addition, particle-flow event reconstruction provides a powerful tool for the rejection of beam-induced backgrounds [11]. The CLIC detector concepts employ strong central solenoid magnets, located outside the HCAL, providing an axial magnetic field of 5 T in CLIC_SiD and 4 T in CLIC_ILD. The CLIC_SiD concept employs central silicon-strip tracking detectors, whereas CLIC_ILD assumes a large central gaseous Time Projection Chamber. In both concepts, the central tracking system is augmented with silicon-based inner tracking detectors. The two detector concepts are shown schematically in Fig. Fig.22 and are described in detail in [11].

Fig. 2
Longitudinal cross section of the top right quadrant of the CLIC_ILD (a) and CLIC_SiD (b) detector concepts

Assumed staged running scenario

The studies presented in this paper are based on a scenario in which CLIC runs at three energy stages. The first stage is at s=350GeV, around the top-pair production threshold. The second stage is at s=1.4TeV; this energy is chosen because it can be reached with a single CLIC drive-beam complex. The third stage is at s=3TeV; the ultimate energy of CLIC. At each stage, four to five years of running with a fully commissioned accelerator is foreseen, providing integrated luminosities of 500 fb-1, 1.5 and 2 ab-1 at 350 GeV, 1.4 and 3 TeV, respectively.1 Cross sections and integrated luminosities for the three stages are summarised in Table Table11.

Table 1
Leading-order, unpolarised cross sections for Higgsstrahlung, WW-fusion, and ZZ-fusion processes for mH = 126 GeV at the three centre-of-mass energies discussed in this paper. s is the effective centre-of-mass energy of the e+e- collision. The ...

Overview of Higgs production at CLIC

A high-energy e+e- collider such as CLIC provides an experimental environment that allows the study of Higgs boson properties with high precision. The evolution of the leading-order e+e- Higgs production cross sections with centre-of-mass energy, as computed using the Whizard 1.95 [20] program, is shown in Fig. Fig.33 for a Higgs boson mass of 126 GeV [21].

Fig. 3
Cross section as a function of centre-of-mass energy for the main Higgs production processes at an e+e- collider for a Higgs mass of mH = 126 GeV. The values shown correspond to unpolarised beams and do not include the effect of beamstrahlung

The Feynman diagrams for the three highest cross section Higgs production processes at CLIC are shown in Fig. Fig.4.4. At s350GeV, the Higgsstrahlung process (e+e- → ZH) has the largest cross section, but the WW-fusion process (e+e-Hνeν¯e) is also significant. The combined study of these two processes probes the Higgs boson properties (width and branching ratios) in a model-independent manner. In the higher energy stages of CLIC operation (s=1.4TeV and 3 TeV), Higgs production is dominated by the WW-fusion process, with the ZZ-fusion process (e+e- → He+e-) also becoming significant. Here the increased WW-fusion cross section, combined with the high luminosity of CLIC, results in large data samples, allowing precise 𝒪(1%) measurements of the couplings of the Higgs boson to both fermions and gauge bosons. In addition to the main Higgs production channels, rarer processes such as e+e-tt¯H and e+e-HHνeν¯e, provide access to the top Yukawa coupling and the Higgs trilinear self-coupling. Feynman diagrams for these processes are shown in Fig. Fig.5.5. In all cases, the Higgs production cross sections can be increased with polarised electron (and positron) beams as discussed in Sect. 3.2.

Fig. 4
Leading-order Feynman diagrams of the highest cross section Higgs production processes at CLIC; Higgsstrahlung (a), WW-fusion (b) and ZZ-fusion (c)
Fig. 5
Feynman diagrams of the leading-order processes at CLIC involving (a) the top Yukawa coupling gHtt, and (b) the Higgs boson trilinear self-coupling λ

Table Table11 lists the expected numbers of ZH, Hνeν¯e and He+e- events for the three main CLIC centre-of-mass energy stages. These numbers account for the effect of beamstrahlung and initial state radiation (ISR), which result in a tail in the distribution of the effective centre-of-mass energy s. The impact of beamstrahlung on the expected numbers of events is mostly small. For example, it results in an approximately 10% reduction in the numbers of Hνeν¯e events at s>1TeV (compared to the beam spectrum with ISR alone), because the cross section rises relatively slowly with s. The reduction of the effective centre-of-mass energies due to ISR and beamstrahlung increases the ZH cross section at s=1.4 and 3 TeV.

The polar angle distributions for single Higgs production obtained using Whizard 1.95 [20] for the CLIC centre-of-mass energies are shown in Fig. Fig.6.6. Most Higgs bosons produced at s=350GeV can be reconstructed in the central parts of the detectors while Higgs bosons produced in the WW-fusion process and their decay products tend towards the beam axis with increasing energy. Hence good detector capabilities in the forward regions are crucial at s=1.4 and 3 TeV.

Fig. 6
Polar angle distributions for single Higgs events at s=350GeV, 1.4 and 3 TeV, including the effects of the CLIC beamstrahlung spectrum and ISR. The distributions are normalised to unity

A SM Higgs boson with mass of mH = 126 GeV has a wide range of decay modes, as listed in Table Table2,2, providing the possibility to test the SM predictions for the couplings of the Higgs to both gauge bosons and to fermions [22]. All the modes listed in Table Table22 are accessible at CLIC.

Table 2
The investigated SM Higgs decay modes and their branching ratios as well as the total Higgs width for mH = 126 GeV [22]

Motivation for s=350 GeV CLIC operation

The choice of the CLIC energy stages is motivated by the desire to pursue a programme of precision Higgs physics and to operate the machine above 1 TeV at the earliest possible time; no CLIC operation is foreseen below the top-pair production threshold.

From the Higgs physics perspective, operation at energies much below 1 TeV is motivated by the direct and model-independent measurement of the coupling of the Higgs boson to the Z, which can be obtained from the recoil mass distribution in ZH → e+e-H, ZH → μ+μ-H and ZHqq¯H production (see Sects. 5.1.1, 5.1.2). These measurements play a central role in the determination of the Higgs couplings at an e+e- collider.

However, from a Higgs physics perspective, there is no advantage to running CLIC at around s=250GeV where the ZH production cross section is larger, compared to running at s=350GeV. Firstly, the reduction in cross section at s=350GeV is compensated, in part, by the increased instantaneous luminosity achievable at a higher centre-of-mass energy. The instantaneous luminosity scales approximately linearly with the centre-of-mass energy, ℒ ∝ γe, where γe is the Lorentz factor for the beam electrons/positrons. For this reason, the precision on the coupling gHZZ at 350 GeV is comparable to that achievable at 250 GeV for the same period of operation. Secondly, the additional boost of the Z and H at s=350GeV provides greater separation between the final-state jets from Z and H decays. Consequently, the measurements of σ(ZH) × BR(H → X) are more precise at s=350GeV. Thirdly, and most importantly, operation of CLIC at s350GeV provides access to the e+e-Hνeν¯e fusion process; this improves the precision with which the total decay width ΓH can be determined at CLIC. For the above reasons, the preferred option for the first stage of CLIC operation is s350GeV.

Another advantage of s350GeV is that detailed studies of the top-pair production process can be performed in the initial stage of CLIC operation. Finally, the Higgs boson mass can be measured at s=350GeV with similar precision as at s=250GeV.

Impact of beam polarisation

The majority of CLIC Higgs physics studies presented in this paper are performed assuming unpolarised e+ and e- beams. However, in the baseline CLIC design, the electron beam can be polarised up to ±80%. There is also the possibility of positron polarisation at a lower level, although positron polarisation is not part of the baseline CLIC design. For an electron polarisation of P- and positron polarisation of P+, the relative fractions of collisions in the different helicity states are:

eR-eR+:14(1+P-)(1+P+),eR-eL+:14(1+P-)(1-P+),eL-eR+:14(1-P-)(1+P+),eL-eL+:14(1-P-)(1-P+).

By selecting different beam polarisations it is possible to enhance/suppress different physical processes. The chiral nature of the weak coupling to fermions results in significant possible enhancements in WW-fusion Higgs production, as indicated in Table Table3.3. The potential gains for the s-channel Higgsstrahlung process, e+e- → ZH, are less significant, and the dependence of the e+e- → He+e- cross section on beam polarisation is even smaller. In practice, the balance between operation with different beam polarisations will depend on the CLIC physics programme taken as a whole, including the searches for and potential measurements of BSM particle production.

Table 3
The dependence of the event rates for the s-channel e+e- → ZH process and the pure t-channel e+e-Hνeν¯e and e+e- → He+e- processes for several example beam polarisations. The scale factors assume an effective weak mixing ...

Overview of Higgs measurements at s=350 GeV

The Higgsstrahlung process, e+e- → ZH, provides an opportunity to study the couplings of the Higgs boson in an essentially model-independent manner. Such a model-independent measurement is unique to a lepton collider. Higgsstrahlung events can be selected based solely on the measurement of the four-momentum of the Z boson through its decay products, while the invariant mass of the system recoiling against the Z boson peaks at mH. The most distinct event topologies occur for Z → e+e- and Z → μ+μ- decays, which can be identified by requiring that the di-lepton invariant mass is consistent with mZ (see Sect. 5.1.1). SM background cross sections are relatively low. A slightly less clean, but more precise, measurement is obtained from the recoil mass analysis for Zqq¯ decays (see Sect. 5.1.2).

Recoil-mass studies provide an absolute measurement of the total ZH production cross section and a model-independent measurement of the coupling of the Higgs to the Z boson, gHZZ. The combination of the leptonic and hadronic decay channels allows gHZZ to be determined with a precision of 0.8%. In addition, the recoil mass from Zqq¯ decays provides a direct search for possible Higgs decays to invisible final states, and can be used to constrain the invisible decay width of the Higgs, Γinvis.

By identifying the individual final states for different Higgs decay modes, precise measurements of the Higgs boson branching fractions can be made. Because of the high flavour tagging efficiencies [11] achievable at CLIC, the Hbb¯ and Hcc¯ decays can be cleanly separated. Neglecting the Higgs decays into light quarks, the branching ratio of H → gg can also be inferred and H → τ+τ- decays can be identified.

Although the cross section is lower, the t-channel WW-fusion process e+e-Hνeν¯e is an important part of the CLIC Higgs physics programme at s350GeV. Because the visible final state consists of the Higgs boson decay products alone, the direct reconstruction of the invariant mass of the Higgs boson or its decay products plays a central role in the event selection. The combination of Higgs production and decay data from Higgsstrahlung and WW-fusion processes provides a model-independent extraction of Higgs couplings.

Extraction of Higgs couplings

At the LHC, only the ratios of the Higgs boson couplings can be inferred from the data in a model-independent way.

In contrast, at an electron–positron collider such as CLIC, absolute measurements of the couplings to the Higgs boson can be determined using the total e+e- → ZH cross section determined from recoil mass analyses. This allows the coupling of the Higgs boson to the Z to be determined with a precision of better than 1% in an essentially model-independent manner. Once the coupling to the Z is known, the Higgs coupling to the W can be determined from, for example, the ratios of Higgsstrahlung to WW-fusion cross sections:

σ(e+e-ZH)×BR(Hbb¯)σ(e+e-νeν¯eH)×BR(Hbb¯)gHZZgHWW2.

Knowledge of the Higgs total decay width, extracted from the data, allows absolute measurements of the other Higgs couplings.

For a Higgs boson mass of around 126 GeV, the total Higgs decay width in the SM (ΓH) is less than 5 MeV and cannot be measured directly at an e+e- linear collider. However, as the absolute couplings of the Higgs boson to the Z and W can be determined, the total decay width of the Higgs boson can be determined from H → WW or H → ZZ decays. For example, the measurement of the Higgs decay to WW in the WW-fusion process determines:

σ(Hνeν¯e)×BR(HWW)gHWW4ΓH,

and thus the total width can be determined utilising the model-independent measurement of gHWW. In practice, a fit (see Sect. 12) is performed to all of the experimental measurements involving the Higgs boson couplings.

Overview of Higgs measurements at s>1 TeV

For CLIC operation above 1 TeV, the large number of Higgs bosons produced in the WW-fusion process allow relative couplings of the Higgs boson to the W and Z bosons to be determined at the 𝒪(1%) level. These measurements provide a strong test of the SM prediction for:

gHWW/gHZZ = cos2θW

where θW is the weak-mixing angle. Furthermore, the exclusive Higgs decay modes can be studied with significantly higher precision than at s=350GeV. For example, CLIC operating at 3 TeV yields a statistical precision of 2% on the ratio gHcc/gHbb, providing a direct comparison of the SM coupling predictions for up-type and down-type quarks. In the context of the model-independent measurements of the Higgs branching ratios, the measurement of σ(Hνeν¯e)×BR(HWW) is particularly important. For CLIC operation at s1.4TeV, the large number of events allows this cross section to be determined with a precision of 1% (see Sect. 6.3). When combined with the measurements at s350GeV, this places a strong constraint on ΓH.

Although the WW-fusion process has the largest cross section for Higgs production above 1 TeV, other processes are also important. For example, measurements of the ZZ-fusion process provide further constraints on the gHZZ coupling. Moreover, CLIC operation at s=1.4TeV enables a determination of the top Yukawa coupling from the process e+e-tt¯HbW+b¯W-H with a precision of 4.2% (see Sect. 8). Finally, the self-coupling of the Higgs boson at the HHH vertex is measurable in 1.4 and 3 TeV operation.

In the SM, the Higgs boson originates from a doublet of complex scalar fields ϕ described by the potential:

V(ϕ) = μ2ϕϕλ(ϕϕ)2

where μ and λ are the parameters of the Higgs potential, with μ2 < 0 and λ > 0. The measurement of the strength of the Higgs self-coupling provides direct access to the coupling λ assumed in the Higgs mechanism. For mH of around 126 GeV, the measurement of the Higgs boson self-coupling at the LHC will be extremely challenging, even with 3000 fb-1 of data (see for example [24]). At a linear collider, the trilinear Higgs self-coupling can be measured through the e+e- → ZHH and e+e-HHνeν¯e processes. The e+e- → ZHH process at s=500GeV has been studied in the context of the ILC, where the results show that a very large integrated luminosity is required [25]. However for s1TeV, the sensitivity for the process e+e-HHνeν¯e increases with increasing centre-of-mass energy and the measurement of the Higgs boson self-coupling (see Sect. 9) forms a central part of the CLIC Higgs physics programme. Ultimately a precision of approximately 20% on λ can be achieved.

Event generation, detector simulation and reconstruction

The results presented in this paper are based on detailed Monte Carlo (MC) simulation studies including the generation of a complete set of relevant SM background processes, Geant4  [26, 27] based simulations of the CLIC detector concepts, and a full reconstruction of the simulated events.

Event generation

Because of the presence of beamstrahlung photons in the colliding electron and positron beams, it is necessary to generate MC event samples for e+e-, e+γ, γe-, and γγ interactions. The main physics backgrounds, with up to six particles in the final state, are generated using the Whizard 1.95 [20] program. In all cases the expected energy spectra for the CLIC beams, including the effects from beamstrahlung and the intrinsic machine energy spread, are used for the initial-state electrons, positrons and beamstrahlung photons. In addition, low-Q2 processes with quasi-real photons are described using the Weizsäcker-Williams approximation as implemented in Whizard. The process of fragmentation and hadronisation is simulated using Pythia 6.4 [28] with a parameter set tuned to OPAL e+e- data recorded at LEP [29] (see [11] for details). The decays of τ leptons are simulated using Tauola  [30]. The mass of the Higgs boson is taken to be 126 GeV 2 and the decays of the Higgs boson are simulated using Pythia with the branching fractions listed in [22]. The events from the different Higgs production channels are simulated separately. The background samples do not include Higgs processes. MC samples for the measurement of the top Yukawa coupling measurement (see Sect. 8) with eight final-state fermions are obtained using the PhysSim  [31] package; again Pythia is used for fragmentation, hadronisation and the Higgs boson decays.

Simulation and reconstruction

The Geant4 detector simulation toolkits Mokka  [32] and Slic  [33] are used to simulate the detector response to the generated events in the CLIC_ILD and CLIC_SiD concepts, respectively. The QGSP_BERT physics list is used to model the hadronic interactions of particles in the detectors. The digitisation, namely the translation of the raw simulated energy deposits into detector signals, and the event reconstruction are performed using the Marlin  [34] and org.lcsim [35] software packages. Particle flow reconstruction is performed using PandoraPFA  [15, 16, 36].

Vertex reconstruction and heavy flavour tagging are performed using the LcfiPlus program [37]. This consists of a topological vertex finder that reconstructs secondary interactions, and a multivariate classifier that combines several jet-related variables such as track impact parameter significance, decay length, number of tracks in vertices, and vertex masses, to tag bottom, charm, and light-quark jets. The detailed training of the multivariate classifiers for the flavour tagging is performed separately for each centre-of-mass energy and each final state of interest.

Because of the 0.5 ns bunch spacing in the CLIC beams, the pile-up of beam-induced backgrounds can affect the event reconstruction and needs to be accounted for. Realistic levels of pile-up from the most important beam-induced background, the γγ → hadrons process, are included in all the simulated event samples to ensure that the impact on the event reconstruction is correctly modelled. The γγ → hadrons events are simulated separately and a randomly chosen subset, corresponding to 60 bunch crossings, is superimposed on the physics event before the digitisation step [38]. 60 bunch crossings is equivalent to 30 ns, which is much longer than the assumed offline event reconstruction window of 10 ns around the hard physics event, so this is a good approximation [11]. For the s=350GeV samples, where the background rates are lower, 300 bunch crossings are overlaid on the physics event. The impact of the background is small at s=350GeV, and is most significant at s=3TeV, where approximately 1.2 TeV of energy is deposited in the calorimeters in a time window of 10 ns. A dedicated reconstruction algorithm identifies and removes approximately 90% of these out-of-time background particles using criteria based on the reconstructed transverse momentum pT of the particles and the calorimeter cluster time. A more detailed description can be found in [11].

Jet finding is performed on the objects reconstructed by particle flow, using the FastJet  [39] package. Because of the presence of pile-up from γγ → hadrons, it was found that the Durham [40] algorithm employed at LEP is not optimal for CLIC studies. Instead, the hadron-collider inspired kt algorithm [41, 42], with the distance parameter R based on Δη and Δϕ, is found to give better performance since it increases distances in the forward region, thus reducing the clustering of the (predominantly low transverse momentum) background particles together with those from the hard e+e- interaction. Instead, particles that are found by the kt algorithm to be closer to the beam axis than to any other particles, and that are thus likely to have originated from beam-beam backgrounds, are removed from the event. As a result of using the R-based kt algorithm, the impact of the pile-up from γγ → hadrons is largely mitigated, even without the timing and momentum cuts described above. Further details are given in [11]. The choice of R is optimised separately for different analyses. In many of the following studies, events are forced into a particular N-jet topology. The variable yij is the smallest kt distance when combining j jets to i = (j - 1) jets. These resolution parameters are widely used in a number of event selections, allowing events to be categorised into topologically different final states. In several studies it is found to be advantageous first to apply the kt algorithm to reduce the beam-beam backgrounds, and then to use only the remaining objects as input to the Durham algorithm.

To recover the effect of bremsstrahlung photons radiated from reconstructed leptons, all photons in a cone around the flight direction of a lepton candidate are added to its four-momentum. The impact of the bremsstrahlung recovery on the reconstruction of the Z → e+e- decays is illustrated in Fig. Fig.7.7. The bremsstrahlung effect leads to a tail at lower values in the Z candidate invariant mass distribution. This loss can be recovered by the procedure described above. It is also visible that a too large opening angle of the recovery cone leads to a tail at higher masses; typically, an opening angle of 3 is chosen.

Fig. 7
Reconstructed invariant mass of Z → e+e- candidates in e+e- → ZH → ZWW events at s=350GeV. Bremsstrahlung photons in cones of different opening angles around the electron direction are recovered as described in the text. All distributions are normalised ...

The event simulation and reconstruction of the large data samples used in this study was performed using the iLCDirac [43, 44] grid production tools.

Higgs production at s=350 GeV

The study of the Higgsstrahlung process is central to the precision Higgs physics programme at any future high-energy electron–positron collider [45]. This section presents studies of e+e- → ZH at s=350GeV with a focus on model-independent measurements of ZH production from the kinematic properties of the Z decay products. Complementary information obtained from Higgs production through WW-fusion at s=350GeV is also presented. All analyses at s=350GeV described in this paper use the CLIC_ILD detector model.

Recoil mass measurements of e+e- → ZH

In the process e+e- → ZH, it is possible to identify efficiently Z → e+e- and Z → μ+μ- decays with a selection efficiency that is essentially independent of the H decay mode. The four-momentum of the (Higgs boson) system recoiling against the Z can be obtained from Erec=s-EZ and prec =  - pZ, and the recoil mass, mrec, peaks sharply around mH. The recoil mass analysis for leptonic decays of the Z is described in Sect. 5.1.1. While these measurements provide a clean model-independent probe of ZH production, they are limited by the relatively small leptonic branching ratios of the Z. Studies of ZH production with Zqq¯ are inherently less clean, but are statistically more powerful. Despite the challenges related to the reconstruction of hadronic Z decays in the presence of various Higgs decay modes, a precise and nearly model-independent probe of ZH production can be obtained by analysing the recoil mass in hadronic Z decays, as detailed in Sect. 5.1.2. When all these measurements are taken together, a model-independent measurement of the gHZZ coupling constant with a precision of  < 1% can be inferred [45].

Leptonic decays: Z → e+e- and Z → μ+μ-

The signature for e+e- → ZH production with Z → e+e- or Z → μ+μ- is a pair of oppositely charged high-pT leptons, with an invariant mass consistent with that of the Z boson, mll ≈ mZ, and a recoil mass, calculated from the four-momenta of the leptons alone, consistent with the Higgs mass, mrec ≈ mH [46]. Backgrounds from two-fermion final states e+e- → l+l- (l = e, μ, τ) are trivial to remove. The dominant backgrounds are from four-fermion processes with final states consisting of a pair of oppositely-charged leptons and any other possible fermion pair. For both the μ+μ-X and e+e-X channels, the total four-fermion background cross section is approximately one thousand times greater than the signal cross section.

The event selection employs preselection cuts and a multivariate analysis. The preselection requires at least one negatively and one positively charged lepton of the lepton flavour of interest (muons or electrons) with an invariant mass loosely consistent with the mass of the Z boson, 40 GeV < mll < 126 GeV. For signal events, the lepton identification efficiencies are 99% for muons and 90% for electrons. Backgrounds from two-fermion processes are essentially eliminated by requiring that the di-lepton system has pT > 60 GeV. Four-fermion backgrounds are suppressed by requiring 95 GeV < mrec < 290 GeV. The lower bound suppresses e+e- → ZZ production. The upper bound is significantly greater than the Higgs boson mass, to allow for the possibility of ZH production with ISR or significant beamstrahlung, which, in the recoil mass analysis, results in a tail to the recoil mass distribution, as it is the mass of the system that is estimated.

Events passing the preselection cuts are categorised using a multivariate analysis of seven discriminating variables: the transverse momentum (pT) and invariant mass (mll) of the candidate Z; the cosine of the polar angle (|cosθ|) of the candidate Z; the acollinearity and acoplanarity of the leptons; the imbalance between the transverse momenta of the two selected leptons (pT1pT2); and the transverse momentum of the highest energy photon in the event. The event selection employs a Boosted Decision Tree (BDT) as implemented in Tmva  [47]. The resulting selection efficiencies are summarised in Table Table4.4. For both final states, the number of selected background events is less than twice the number of selected signal events. The impact of the background is reduced using a fit to the recoil mass distribution.

Table 4
Preselection and selection efficiencies for the ZH signal and most important background processes in the leptonic recoil mass analysis. The numbers of events correspond to 500 fb-1 at s=350GeV

A fit to the recoil mass distribution of the selected events (in both the Z → e+e- and Z → μ+μ- channels) is used to extract measurements of the ZH production cross section and the Higgs boson mass. The shape of the background contribution is parameterised using a fourth order polynomial and the shape of the signal distribution is modelled using Simplified Kernel Estimation [4850] that provides a description of the ZH recoil mass distribution in which the Higgs mass can subsequently be varied. The accuracy with which the Higgs mass and the number of signal events (and hence the ZH production cross section) can be measured, is determined using 1000 simulated test data samples. Each test sample was created by adding the high statistics selected signal sample (scaled to the correct normalisation) to the smooth fourth-order polynomial background,then applying Poisson fluctuations to individual bins to create a representative 500 fb-1 data sample. Each of the 1000 simulated data samples created in this way is fitted allowing the Higgs mass, the signal normalisation and the background normalisation to vary. Figure Figure8a8a displays the results of fitting a typical test sample for the μ+μ-X channel, while Fig. Fig.8b8b displays the results for the e+e-X channel. In the e+e-X channel fits are performed with, and without, applying an algorithm to recover bremsstrahlung photons. The resulting measurement precisions for the ZH cross section and the Higgs boson mass are summarised in Table Table5.5. In the e+e-X channel, the bremsstrahlung recovery leads to a moderate improvement on the expected precision for the cross section measurement and a similar degradation in the expected precision for the mass determination, because it significantly increases the number of events in the peak of the recoil mass distribution, but also increases the width of this peak. For an integrated luminosity of 500 fb-1 at s=350GeV, the combined precision on the Higgs boson mass is:

Δ(mH) = 110 MeV, 

and the combined precision on the ZH cross section is:

Δσ(ZH)σ(ZH)=3.8%.

The expected precision with (without) bremsstrahlung recovery in the e+e-X channel was used in the combination for the cross section (mass).

Fig. 8
Reconstructed recoil mass distributions of e+e- → ZH events at s=350GeV, where ZH → μ+μ-X (a) and ZH → e+e-X with bremsstrahlung recovery (b). All distributions are normalised to an integrated luminosity of 500 fb-1
Table 5
Summary of measurement precisions from the leptonic recoil mass analyses in the μ+μ-X and e+e-X channels for an integrated luminosity of 500 fb-1 at 350 GeV

Hadronic decays: Zqq¯

In the process e+e- → ZH, it is possible to cleanly identify Z → e+e- and Z → μ+μ- decays regardless of the decay mode of the Higgs boson and, consequently, the selection efficiency is almost independent of the Higgs decay mode. In contrast, for Zqq¯ decays, the selection efficiency shows a stronger dependence on the Higgs decay mode [45]. For example, e+e-(Zqq¯)(Hbb¯) events consist of four jets and the reconstruction of the Z boson is complicated by ambiguities in associations of particles with jets and the three-fold ambiguity in associating four jets with the hadronic decays of the Z and H. For this reason, it is much more difficult to construct a selection based only on the reconstructed Zqq¯ decay that has a selection efficiency independent of the Higgs decay mode. The strategy adopted is to first reject events consistent with a number of clear background topologies using the information from the whole event; and then to identify e+e-(Zqq¯)H events solely based on the properties from the candidate Zqq¯ decay.

The (Zqq¯)H event selection proceeds in three separate stages. In the first stage, to allow for possible BSM invisible Higgs decay modes, events are divided into candidate visible Higgs decays and candidate invisible Higgs decays, in both cases produced along with a Zqq¯. Events are categorised as potential visible Higgs decays if they are not compatible with a clear two-jet topology:

  • log10(y23) >  - 2.0    or    log10(y34) >  - 3.0.

All other events are considered as candidates for an invisible Higgs decay analysis, based on that described in Sect. 5.1.3, although with looser requirements to make the overall analysis more inclusive.

Preselection cuts then reduce the backgrounds from large cross section processes such as e+e-qq¯ and e+e-qq¯qq¯. The preselection variables are formed by forcing each event into three, four and five jets. In each case, the best candidate for being a hadronically decaying Z boson is chosen as the jet pair giving the di-jet invariant mass (mqq¯) closest to mZ, considering only jets with more than three charged particles. The invariant mass of the system recoiling against the Z boson candidate, mrec, is calculated assuming Erec=s-Eqq¯ and prec=-pqq¯. In addition, the invariant mass of all the visible particles not originating from the candidate Zqq¯ decay, mvis, is calculated. It is important to note that mvis is only used to reject specific background topologies in the preselection and is not used in the main selection as it depends strongly on the type of Higgs decay. The preselection cuts are:

  • 70GeV<mqq¯<110GeV and 80 GeV < mrec < 200 GeV;
  • the background from e+e-qq¯ is suppressed by removing events with overall pT < 20 GeV and either |cosθmis| > 0.90 or log10(y34) >  - 2.5, where θmis is the polar angle of the missing momentum vector;
  • events with little missing transverse momentum (pT < 20 GeV) are forced into four jets and are rejected if the reconstructed di-jet invariant masses (and particle types) are consistent with the expectations for e+e-qq¯ll, e+e-ZZqq¯qq¯, e+e-WWqq¯qq¯.

The final step in the event selection is a multivariate analysis. In order not to bias the event selection efficiencies for different Higgs decay modes, only variables related to the candidate Zqq¯ decay are used in the selection. Forcing the event into four jets is the right approach for (Zqq¯)H events where the Higgs decays to two-body final states, but not necessarily for final states such as HWWqq¯qq¯, where there is the chance that one of the jets from the WW decay will be merged with one of the jets from the Zqq¯, potentially biasing the selection against H → WW decays. To mitigate this effect, the Z candidate for the event selection can either be formed from the four-jet topology as described above, or can be formed from a jet pair after forcing the event into a five-jet topology. The latter case is only used when log10(y45) >  - 3.5 and the five-jet reconstruction gives better Z and H candidates than the four-jet reconstruction. Attempting to reconstruct events in the six-jet topology is not found to improve the overall analyses. Having chosen the best Z candidate in the event (from either the four-jet or five-jet reconstruction), it is used to form variables for the multivariate selection; information about the remainder of the event is not used.

A relative likelihood selection is used to classify all events passing the preselection cuts. Two event categories are considered: the e+e-ZHqq¯H signal and all non-Higgs background processes. The relative likelihood for an event being signal is defined as:

L=LsignalLsignal+Lback,

where the individual absolute likelihood L for each event type is estimated from normalised probability distributions, Pi(xi), of the discriminating variables xi for that event type:

L=σpresel×iNPi(xi),

where σpresel is the cross section after the preselection cuts. The discriminating variables used, all of which are based on the candidate Zqq¯ decay, are: the 2D distribution of mqq¯ and mrec; the polar angle of the Z candidate, |cosθZ|; and the modulus of angle of jets from the Z decay relative to its direction after boosting into its rest frame, |cosθq|. The clearest separation between signal and background is obtained from mqq¯ and the recoil mass mrec, as shown in Fig. Fig.99 for events passing the preselection. The signal is clearly peaked at mqq¯mZ and mrec ≈ mH. The use of 2D mass distributions accounts for the most significant correlations between the likelihood variables.

Fig. 9
Reconstructed di-jet invariant mass versus reconstructed recoil mass distributions for ZHqq¯X candidate events at s=350GeV, showing ZH signal events (a) and all background processes (b). In both cases the plots show all events passing ...

In this high-statistics limit, the fractional error on the number of signal events (where the Higgs decays to visible final states), svis, given a background b is:

Δsvissvis=svis+bsvis,

and this is minimised with the selection requirement ℒ > 0.65. The selection efficiencies and expected numbers of events for the signal dominated region, ℒ > 0.65, are listed in Table Table6,6, corresponding to a fractional error on the number of signal events of 1.9%. By fitting the shape of the likelihood distribution to signal and background contributions, this uncertainty is reduced to:

Δsvissvis=1.7%.

This is an example of a measurement for which it will be particularly important to tune the background modelling using high-statistics processes.

Table 6
Summary of the (Zqq¯)(Hvis.) event selection at s=350GeV, giving the raw cross sections, preselection efficiency, overall selection efficiency for a likelihood cut of ℒ > 0.65 and the expected numbers of events passing the event ...

Invisible Higgs decays

The above recoil mass analysis of leptonic decays of the Z boson in e+e- → ZH events provides a measurement of the Higgsstrahlung cross section, independent of the Higgs boson decay model. The recoil mass technique can also be used to search for BSM decay modes of the Higgs boson into long-lived neutral “invisible” final states [45]. At an e+e- collider, a search for invisible Higgs decays is possible by identification of e+e- → ZH events with a visible Zqq¯ decay and missing energy. Such events would typically produce a clear two-jet topology with invariant mass consistent with mZ, significant missing energy and a recoil mass corresponding to the Higgs mass. Higgsstrahlung events with leptonic Z decays, which have a much smaller branching ratio, are not included in the current analysis.

To identify candidate invisible Higgs decays, a loose preselection is imposed requiring: (i) a clear two-jet topology, defined by log10(y23) <  - 2.0 and log10(y34) <  - 3.0, using the minimal kt distances discussed in Sect. 4.2; (ii) a di-jet invariant mass consistent with mZ, 84GeV<mqq¯<104GeV; and (iii) the reconstructed momentum of the candidate Z boson pointing away from the beam direction, |cosθZ| < 0.7. After the preselection, a BDT multivariate analysis technique is applied using the Tmva package [47] to further separate the invisible Higgs signal from the SM background. In addition to mqq¯, |cosθZ| and log10(y23), four other discriminating variables are employed: mrec, the recoil mass of the invisible system recoiling against the observed Z boson; |cosθq|, the decay angle of one of the quarks in the Z rest frame, relative to the direction of flight of the Z boson; pT, the magnitude of the transverse momentum of the Z boson; and Evis, the visible energy in the event. As an example, Fig. Fig.1010 shows the recoil mass distribution for the simulated invisible Higgs decays and the total SM background. The reconstructed recoil mass for events with invisible Higgs decays peaks near mH. The cut applied on the BDT output is chosen to minimise the statistical uncertainty on the cross section for invisible Higgs decays.

Fig. 10
Reconstructed recoil mass distributions of e+e- → ZH events at s=350GeV, showing the H → invis. signal, assuming BR(H → invis.) = 100%, and SM backgrounds as stacked histograms. The distributions are normalised to an integrated luminosity of 500 fb-1 ...

In the case where the branching ratio to BSM invisible final states is zero (or very small), the uncertainty on the invisible branching ratio is determined by the statistical fluctuations on the background after the event selection:

ΔBR(Hinvis.)=bs(100%),

where b is the expected number of selected SM background events and s(100%) is the expected number of selected Higgsstrahlung events assuming all Higgs bosons decay invisibly, i.e. BR(H → invis.) = 100%. Table Table77 summarises the invisible Higgs decay event selection; the dominant background processes arise from the final states qq¯lν and qq¯νν¯. The resulting one sigma uncertainty on BR(H → invis.) is 0.57% (in the case where the invisible Higgs branching ratio is small) and the corresponding 90% C.L. upper limit (500 fb-1 at s =350 GeV) on the invisible Higgs branching ratio in the modified frequentist approach [51] is:

BR(H → invis.) < 0.97% at 90% C.L.

It should be noted that the SM Higgs decay chain HZZνν¯νν¯ has a combined branching ratio of 0.1% and is not measurable.

Table 7
Summary of the invisible Higgs decay event selection at s=350GeV, giving the raw cross sections, preselection efficiency, selection efficiency for a BDT cut of BDT > 0.088, and the expected numbers of events passing the event selection for an integrated ...

Model-independent ZH cross section

By combining the two analyses for ZH production where Zqq¯ and the Higgs decays either to invisible final states (see Sect. 5.1.3) or to visible final states (see Sect. 5.1.2), it is possible to determine the absolute cross section for e+e-  → ZH in an essentially model-independent manner:

σ(ZH)=σvis+σinvisBR(Zqq¯).

Here a slightly modified version of the invisible Higgs analysis is employed. With the exception of the cuts on y23 and y34, the invisible Higgs analysis employs the same preselection as for the visible Higgs analysis and a likelihood multivariate discriminant is used.

Since the fractional uncertainties on the total cross section from the visible and invisible cross sections are 1.7 and 0.6% respectively, the fractional uncertainty on the total cross section will be (at most) the quadrature sum of the two fractional uncertainties, namely 1.8%. This measurement is only truly model-independent if the overall selection efficiencies are independent of the Higgs decay mode. For all final state topologies, the combined (visible + invisible) selection efficiency lies is the range 19–26% regardless of the Higgs decay mode, covering a very wide range of event topologies. To assess the level of model independence, the Higgs decay modes in the MC samples are modified and the total (visible + invisible) cross section is extracted assuming the SM Higgs branching ratio. Table Table88 shows the resulting biases in the extracted total cross section for the case when a BR(H → X) → BR(H → X) + 0.05. Even for these very large modifications of the Higgs branching ratios over a wide range of final-state topologies – including the extreme cases highlighted at the bottom of Table Table88 such as HWWqq¯qq¯, which has six jets in the final state, and H → WW → τντν, which has a lot of missing energy – the resulting biases in the extracted total ZH cross section are less than 1% (compared to the 1.8% statistical uncertainty). However, such large deviations would have significant observable effects on exclusive Higgs branching ratio analyses (at both LHC and CLIC) and it is concluded that the analysis gives an effectively model-independent measurement of the (Zqq¯)H cross section.

Table 8
Biases in the extracted H(Zqq¯) cross section if the Higgs branching ratio to a specific final state is increased by 5%, i.e. BR(H → X) → BR(H → X) + 0.05

Combining the model-independent measurements of the ZH cross section from Z → l+l- and Zqq¯ gives an absolute measurement of the ZH cross section with a precision of:

Δσ(ZH)σ(ZH)=1.65%,

and, consequently, the absolute coupling of the H boson to the Z boson is determined to:

ΔgHZZgHZZ=0.8%.

The hadronic recoil mass analysis was repeated for collision energies of s=250GeV and s=420GeV [45]. Compared with s=350GeV, the sensitivity is significantly worse in both cases.

Exclusive Higgs branching ratio measurements at s=350GeV

The previous section described inclusive measurements of the e+e- → ZH production cross section, which provide a model-independent determination of the coupling at the HZZ vertex. In contrast, measurements of Higgs production and decay to exclusive final states provide a determination of the product σ(ZH) × BR(H → X), where X is a particular final state. This section focuses on the exclusive measurements of the Higgs decay branching ratios at s=350GeV. Higgs boson decays to bb¯, cc¯ and gg are studied in Sect. 5.2.1. The measurement of H → τ+τ- decays is described in Sect. 5.2.2, and the H → WW decay mode is described in Sect. 5.2.3.

Hbb¯,cc¯ and gg

As can be seen from Table Table1,1, at s=350GeV the cross section for e+e- → ZH (Higgsstrahlung) is approximately four times greater than the e+e-Hνeν¯e (mostly WW-fusion) cross section for unpolarised beams (or approximately a factor 2.5 with - 80% electron beam polarisation). For Higgsstrahlung, the signature of Hbb¯,cc¯,gg events depends on the Z decay mode.

To maximise the statistical power of the Hbb¯,cc¯,gg branching ratio measurements, two topologies are considered: four jets, and two jets plus missing momentum (from the unobserved neutrinos). The impact of Higgsstrahlung events with leptonic Z decays is found to be negligible. The jets plus missing momentum final state contains approximately equal contributions from Higgsstrahlung and WW-fusion events, although the event kinematics are very different. All events are initially reconstructed assuming both topologies; at a later stage of the event selection, events are assigned to either Hqq¯, Hνν¯, or background. To minimize the impact of ISR on the jet reconstruction, photons with a reconstructed energy higher than 15 GeV are removed from the events first.

The hadronic final states are reconstructed using the Durham algorithm. For the four-jet topology, the most probable Z and Higgs boson candidates are selected by choosing the jet combination that minimises:

χ2=(mij-mH)2/σH2+(mkl-mZ)2/σZ2,

where mij and mkl are the invariant masses of the jet pairs used to reconstruct the Higgs and Z boson candidates, respectively, and σH,Z are the estimated invariant mass resolutions for Higgs and Z boson candidates. In the case of the two jets plus missing energy final state, either from ZH with Zνν¯ or from Hνν¯, the event is clustered into two jets forming the H candidate.

To help veto backgrounds with leptonic final states, isolated electrons or muons with E > 10 GeV are identified with the additional requirement that there should be less than 20 GeV of energy from other particles within a cone with an opening angle of 20 around the lepton direction. All events are then classified by gradient boost decision trees employing reconstructed kinematic variables from each of the two event topology hypotheses described above. The variables used include jet energies, event shape variables (such as thrust and sphericity), the masses of H and Z candidates, their decay angles and transverse momenta, and the number of isolated leptons in the final state. The total number of variables is about 50, which is larger than in other studies presented in this paper, because each event is reconstructed assuming two different final state configurations and information from the H candidate decay can be included here, in contrast with the recoil mass analyses described in Sect. 5.1.

Two separate BDT classifiers are used, one for each signal final state (Hqq¯ and Hνν¯), irrespective of the nature of the hadronic Higgs decay mode. Two-fermion (qq¯) and four-fermion (qq¯νν¯, qq¯lν, qq¯ll and qq¯qq¯) final states and other Higgs decay modes are taken as background for both classifiers. In addition, the other signal mode is included in the background for a given classifier. The training is performed using a dedicated training sample, simultaneously training both classifiers. At this point, no flavour tagging information is used.

Each event is evaluated with both classifiers. An event is only accepted if exactly one of the signal classifiers is above a positive threshold and the other classifier is below a corresponding negative threshold. The event is then tagged as a candidate for the corresponding signal process. If none of the classifiers passes the selection threshold, the event is considered as background and is rejected from the analysis. The number of events for which both signal classifiers are above the positive threshold is negligible. Table Table99 summarises the classification of all events into the two signal categories, with event numbers based on an integrated luminosity of 500 fb-1.

Table 9
Summary of the expected numbers of events for the different Higgs and non-Higgs final states passing the hadronic Higgs decay signal selection for 500 fb-1 at s=350GeV (unpolarised beams). No preselection is applied in this analysis

The second stage of the analysis is to measure the contributions of the hadronic Higgs decays into the Hbb¯, Hcc¯ and H → gg exclusive final states, separated into the two production modes Higgsstrahlung and WW-fusion. This is achieved by a multi-dimensional template fit using flavour tagging information and, in the case of the Hνν¯ final state, the transverse momentum of the Higgs candidate.

The jets forming the Higgs candidate are classified with the LcfiPlus flavour tagging package. Each jet pair is assigned a bb¯ likelihood and a cc¯ likelihood:

bb¯likelihood=b1b2b1b2+(1-b1)(1-b2),

cc¯likelihood=c1c2c1c2+(1-c1)(1-c2),

where b1 and b2 (c1 and c2) are the b-tag (c-tag) values obtained for the two jets forming the Higgs candidate.

The resulting two-dimensional distributions of the bb¯ and cc¯ likelihoods in Hqq¯ events are shown in Fig. Fig.11,11, where separation between the different event categories can be seen. These distributions form the templates used to determine the contribution of the different signal categories for the Hqq¯ final states.

Fig. 11
bb¯ likelihood versus cc¯ likelihood distributions for e+e- → ZH events at s=350GeV, for (a) all events and for the different event classes: (b) H → bb¯, (c) H → cc¯, (d) H → gg, background from (e) other Higgs decays ...

Signal and background templates are also obtained for the Hνν¯ final state. As Hνν¯ has roughly equal contributions from the Higgsstrahlung and the WW-fusion process, separation into the two production processes is required, in addition to separation into the different signal and background final states. This is achieved by adding the transverse momentum of the Higgs candidate to the templates as a third dimension. This exploits the fact that the transverse momentum of the Higgs candidate is substantially different for Higgsstrahlung and WW-fusion events, as illustrated in Fig. Fig.1212 for events with a high bb¯ likelihood, which provides a high signal purity.

Fig. 12
Reconstructed Higgs candidate transverse momentum distributions for selected Hνν¯ events at s=350GeV, showing the contributions from Higgsstrahlung, WW-fusion and non-Higgs background. The distributions are normalised to an integrated ...

Contributions from events with Hbb¯, Hcc¯ and H → gg decays, separated by production mode, are extracted in a template fit maximizing the combined likelihood of the Hqq¯ and Hνν¯ templates. It is assumed that the contributions from other Higgs decay modes are determined from independent measurements and therefore these contributions are fixed in the fit.

The results of the above analysis are summarised in Table Table10,10, giving the statistical uncertainties of the various σ × BR measurements. Since the parameters in this analysis are determined in a combined extraction from overlapping distributions, the results are correlated. In particular the Higgsstrahlung and WW-fusion results for the same final states show sizeable anti-correlations, as large as - 38% for the cases of Hcc¯ and H → gg. These correlations are taken into account in the global fits described in Sect. 12.

Table 10
Summary of statistical uncertainties for events with a Hbb¯, Hcc¯ or H → gg decay, where the Higgs boson is produced by Higgsstrahlung or WW-fusion, at s=350GeV derived from the template fit as described in the text. ...

H → τ+τ-

Because of the neutrino(s) produced in τ decays, the signature for H → τ+τ- is less distinct than that for other decay modes. The invariant mass of the visible decay products of the τ+τ- system will be less than mH, and it is difficult to identify H → τ+τ- decays from the WW-fusion process or from Higgsstrahlung events where Zνν¯. For this reason, the product of σ(ZH) × BR(H → τ+τ-) is only determined for the case of hadronic Z decays at s=350GeV. In this analysis only hadronic τ decays are considered, so the experimental signature is two hadronic jets from Zqq¯ and two isolated low-multiplicity narrow jets from the two tau decays [52]. Candidate τ leptons are identified using the TauFinder algorithm [53], which is a seeded-cone based jet-clustering algorithm. The algorithm was optimised to distinguish the tau lepton decay products from hadronic gluon or quark jets. Tau cones are seeded from single tracks (pT > 5 GeV). The seeds are used to define narrow cones of 0.05 rad. The cones are required to contain either one or three charged particles (from one- and three-prong tau decays) and further rejection of background from hadronic jets is implemented using cuts on isolation-related variables. Tau cones which contain identified electrons or muons are rejected and only the hadronic one- and three-prong τ decays are retained. The τ identification efficiency for hadronic tau decays is found to be 73% and the fake rate to mistake a quark for a τ is 5%. The fake rate is relatively high, but is acceptable as the background from final states with quarks can be suppressed using global event properties.

Events with two identified hadronic tau candidates (with opposite net charge) are considered as H → τ+τ- decays. Further separation of the signal and background events is achieved using a BDT classifier based on the properties of the tau candidates and global event properties. Seventeen discriminating variables are used as BDT inputs, including the thrust and oblateness of the quark and tau systems, and masses, transverse momenta, and angles in the events. A full list is given in [52]. The resulting BDT distributions for the signal and the backgrounds are shown in Fig. Fig.13.13. Events passing a cut on the BDT output maximising the significance of the measurement are selected. The cross sections and numbers of selected events for the signal and the dominant background processes are listed in Table Table11.11. The contribution from background processes with photons in the initial state is negligible after the event selection. A template fit to the BDT output distributions leads to:

Δ[σ(ZH)×BR(Hτ+τ-)]σ(ZH)×BR(Hτ+τ-)=6.2%.

Fig. 13
BDT classifier distributions for H → τ+τ- events at s=350GeV, showing the signal and main backgrounds as stacked histograms. The distributions are normalised to an integrated luminosity of 500 fb-1
Table 11
Cross sections and numbers of preselected and selected events with BDT > 0.08 (see Fig. Fig.13)13) for e+e-ZH(Zqq¯,Hτ+τ-) signal events and the dominant backgrounds at s=350GeV assuming ...

H → WW

In case the Higgs boson decays to a pair of W bosons, only the fully hadronic channel, HWWqq¯qq¯, allows the reconstruction of the Higgs invariant mass. Two final states in e+e- → ZH events have been studied depending on the Z boson decay mode: Z → l+l-, where l is an electron or muon, and Zqq¯.

First, isolated electrons and muons from Z decays are identified. Photons in a cone with an opening angle of 3 around the lepton candidates are added to their four-momentum as described in Sect. 4.2.

If a leptonic Z candidate is found, four jets are reconstructed from all particles not originating from the Z decay. The jets are paired, with the pair that gives the mass closest to the W boson mass being taken as one W boson candidate, and the other pair taken as the W. The events are considered further if the invariant mass of the Z boson candidate is in the range between 70 and 110 GeV and at least 20 particles are reconstructed.

In events without a leptonic Z candidate, six jets are reconstructed. The jets are grouped into W, Z and Higgs boson candidates by minimising:

χ2=(mij-mW)2σW2+(mkl-mZ)2σZ2+(mijmn-mH)2σH2,

where mij is the invariant mass of the jet pair used to reconstruct the W candidate, mkl is the invariant mass of the jet pair used to reconstruct the Z candidate, mijmn is the invariant mass of the four jets used to reconstruct the Higgs candidate and σW,Z,H are the estimated invariant mass resolutions for W, Z and Higgs boson candidates. The preselection cuts for this final state are:

  • invariant mass of the Z candidate greater than 40 GeV;
  • at least 50 reconstructed particles;
  • event thrust of less than 0.95;
  • no jet with a b-tag probability of more than 0.95;
  • topology of the hadronic system consistent with six jets: log10(y12) >  - 2.0, log10(y23) >  - 2.6, log10(y34) >  - 3.0, log10(y45) >  - 3.5 and log10(y56) >  - 4.0.

For both final states, BDT classifiers are used to suppress the backgrounds further. The event selection for the signal processes and the most relevant background samples is summarised in Table Table12.12. The expected precisions for the measurement of the investigated processes are summarised in Table Table13.13. The best precision is achieved using the Zqq¯ decay due to its large branching ratio compared to leptonic decays. The selection of Z → e+e- events is more difficult compared to Z → μ+μ- events because the e+e-qq¯ll background sample contains more events with electron pairs than events with muon pairs. Hence the precision achieved using Z → μ+μ- decays is somewhat better compared to that obtained using Z → e+e- decays. The combined precision for an integrated luminosity of 500 fb-1 is:

Δ[σ(ZH)×BR(HWW)]σ(ZH)×BR(HWW)=5.1%,

which is dominated by the final state with hadronic Z boson decays.

Table 12
Preselection and selection efficiencies for the ZH signal and most important background processes of the H → WW analysis in all three considered Z decay channels. The numbers assume an integrated luminosity of 500 fb-1 at s=350GeV
Table 13
Statistical precisions for the listed processes at s=350GeV for an integrated luminosity of 500 fb-1

WW-fusion at s>1TeV

This section presents measurements of Higgs decays from the WW-fusion process at CLIC with centre-of-mass energies of 1.4 TeV and 3 TeV. The Higgs self-coupling measurement, which is also accessed in WW-fusion production, is discussed in Sect. 9. The cross section of the Higgs production via the vector boson fusion process e+e-Hνeν¯e scales with log(s) and becomes the dominating Higgs production process in e+e- collisions with s>500GeV. The respective cross sections for e+e-Hνeν¯e at s=1.4TeV and 3 TeV are approximately 244 fb and 415 fb, respectively, including the effects of the CLIC beamstrahlung spectrum and ISR. The relatively large cross sections at the higher energies allow the Higgs decay modes to be probed with high statistical precision and provide access to rarer Higgs decays, such as H → μ+μ-.

Since WW-fusion e+e-Hνeν¯e proceeds through the t-channel, the Higgs boson is typically boosted along the beam direction and the presence of neutrinos in the final state can result in significant missing pT. Because of the missing transverse and longitudinal momentum, the experimental signatures for Hνeν¯e production are relatively well separated from most SM backgrounds. At s=350GeV, the main SM background processes are two- and four-fermion production, e+e- → 2f and e+e- → 4f. At higher energies, backgrounds from γγ and γe± hard interactions become increasingly relevant for measurements of Higgs boson production in WW-fusion. Additionally, pile-up of relatively soft γγ → hadrons events with the primary interaction occurs. However, this background of relatively low-pT particles is largely mitigated through the timing cuts and jet finding strategy outlined in Sect. 4.

Hbb¯,cc¯,gg

The physics potential for the measurement of hadronic Higgs decays at the centre-of-mass energies of 1.4 and 3 TeV was studied using the CLIC_SiD detector model. The signatures for Hbb¯, Hcc¯ and H → gg decays in e+e-Hνeν¯e events are two jets and missing energy. Flavour tagging information from LcfiPlus is used to separate the investigated Higgs boson decay modes in the selected event sample. The invariant mass of the reconstructed di-jet system provides rejection against background processes, e.g. hadronic Z boson decays.

At both centre-of-mass energies, an invariant mass of the di-jet system in the range from 60 to 160 GeV and a distance between both jets in the ηϕ plane of less than 4 are required. The energy sum of the two jets must exceed 75 GeV and a missing momentum of at least 20 GeV is required. The efficiencies of these preselection cuts on the signal and dominant background samples are listed in Tables Tables1414 and and1515 for the centre-of-mass energies of 1.4 and 3 TeV, respectively.

Table 14
Preselection and selection efficiencies for the signal and most important background processes in the Hbb¯, Hcc¯ and H → gg analysis. The numbers of events correspond to 1.5 ab-1 at s=1.4TeV
Table 15
Preselection and selection efficiencies for the signal and most important background processes in the Hbb¯, Hcc¯ and H → gg analysis. The numbers of events correspond to 2 ab-1 at s=3TeV

The backgrounds are suppressed further using a single BDT at each energy. The samples of signal events used to train these classifiers consist of equal amounts of Hbb¯, Hcc¯, and H → gg events, while the different processes in the background sample were normalised according to their respective cross sections. No flavour tagging information is used in the event selection. This leads to classifiers with similar selection efficiencies for events with the different signal Higgs decays.

The fractions of signal events with Hbb¯, Hcc¯ and H → gg decays in the selected event samples are extracted from the two-dimensional distributions of the bb¯ versus cc¯ likelihood variables for the two reconstructed jets as defined in Sect. 5.2. The normalisations of the backgrounds from other Higgs decays and non-Higgs events are fixed and expected to be provided by other measurements. The results of these fits are shown in Tables Tables1616 and and1717 at 1.4 and 3 TeV, respectively.

Table 16
Statistical precisions for the listed processes from the fit described in the text at s=1.4TeV for an integrated luminosity of 1.5 ab-1
Table 17
Statistical precisions for the listed processes from the fit described in the text at s=3TeV for an integrated luminosity of 2 ab-1

The expected precisions obtained at 1.4 and 3 TeV are similar although the number of signal events is about twice as large at 3 TeV compared to 1.4 TeV. The main reasons for this are that the jet reconstruction and flavour tagging are more challenging at 3 TeV, since the jets from the Higgs decay tend more towards the beam axis, and the impact of the beam-induced backgrounds is larger compared to 1.4 TeV. In addition, the cross sections for the most important background processes rise with s (see Tables Tables1414 and and1515).

H → τ+τ-

The sensitivity for the measurement of σ(e+e-Hνeν¯e)×BR(Hτ+τ-) at CLIC has been studied using the CLIC_ILD detector model at centre-of-mass energies of 1.4 and 3 TeV  [54]. For a SM Higgs with a mass of 126 GeV, BR(H → τ+τ-) = 6.2%, resulting in an effective signal cross section of 15.0 fb at s=1.4TeV and 25.5 fb at s=3TeV.

The experimental signature is two relatively high-momentum narrow jets from the two tau decays and significant missing transverse and longitudinal momenta. A typical event display is shown in Fig. Fig.14.14. The analysis is restricted to hadronic τ decays, which are identified using the TauFinder algorithm, as described in Sect. 5.2.2. The TauFinder algorithm parameters were tuned using the H → τ+τ- signal events and e+e-qq¯νν¯ background events. The working point has a τ selection efficiency of 70% (60%) with a quark jet fake rate of 7% (9%) at s=1.4TeV (s=3TeV). All relevant SM backgrounds are taken into account, including γγ and γe± collisions. The most significant backgrounds are e+e-τ+τ-νν¯, e±γ → τ+τ-e± and γγτ+τ-νν¯. The latter two processes become increasingly important at higher s, due to the increasing number of beamstrahlung photons. Backgrounds from Higgs decays other than H → τ+τ- are expected to be negligible [55].

Fig. 14
Event display of a H → τ+τ- event at s=1.4TeV in the CLIC_ILD detector. A 1-prong tau decay is visible in the central part of the detector (blue). The other tau lepton decays to three charged particles and is reconstructed in the forward direction (red ...

The event preselection requires two identified τ leptons, both of which must be within the polar angle range 15 < θ(τ) < 165 and have pT(τ) > 25 GeV. To reject back-to-back or nearby tau leptons, the angle between the two tau candidates must satisfy 29 < Δθ(ττ) < 177. The visible invariant mass m(ττ) and the visible transverse mass mT(ττ) of the two tau candidates must satisfy 45 GeV < m(ττ) < 130 GeV and mT(ττ) < 20 GeV. Finally the event thrust must be less than 0.99.

Events passing the preselection are classified as either signal or SM background using a BDT classifier. The kinematic variables used in the classifier are m(ττ), mT(ττ), event shape variables (such as thrust and oblateness), the missing pT, the polar angle of the missing momentum vector |cosθmiss| and the total reconstructed energy excluding the Higgs candidate. The event selection for the signal and the most relevant background processes is summarised in Table Table1818 for s=1.4TeV and in Table Table1919 for s=3TeV. Rather than applying a simple cut, the full BDT shape information is used in a template fit. The resulting statistical uncertainties for 1.5 ab-1 at s=1.4TeV and 2.0 ab-1 at s=3TeV are:

Δ[σ(Hνeν¯e)×BR(Hτ+τ-)]σ(Hνeν¯e)×BR(Hτ+τ-)=4.2%at1.4TeV,Δ[σ(Hνeν¯e)×BR(Hτ+τ-)]σ(Hνeν¯e)×BR(Hτ+τ-)=4.4%at3TeV.

Table 18
Preselection and selection efficiencies for the signal and most important background processes in the H → τ+τ- analysis. The numbers of events correspond to 1.5 ab-1 at s=1.4TeV. The cross sections for the backgrounds include cuts on the kinematic ...
Table 19
Preselection and selection efficiencies for the signal and most important background processes in the H → τ+τ- analysis. The numbers of events correspond to 2 ab-1 at s=3TeV. The cross sections for the backgrounds include cuts on the kinematic ...

Similar to the observations described in Sect. 6.1, the expected precisions at 1.4 TeV and 3 TeV are similar. The identification of tau leptons is more challenging at 3 TeV where the impact of the beam-induced backgrounds is larger and the tau leptons from Higgs decays in signal events tend more towards the beam axis.

H → WW

The signature for H → WW decays in e+e-Hνeν¯e depends on the WW decay modes. As mH < 2mW, at least one of the W-bosons is off mass-shell. Studies for two different final states are described in the following. The presence of a charged lepton in the WWqq¯lν final state suppresses backgrounds from other Higgs decays. However, the invariant mass of the Higgs boson in H → WW decays can be reconstructed for fully-hadronic decays alone, WWqq¯qq¯.

WWqq¯qq¯

The experimental signature for Hνeν¯e production with HWWqq¯qq¯ is a four-jet final state with missing pT and a total invariant mass consistent with the Higgs mass, where one pair of jets has a mass consistent with mW.

The H → WW event selection has been studied at s=1.4TeV using the CLIC_ILD detector model. It proceeds in two separate stages: a set of preselection cuts designed to reduce the backgrounds from large cross section processes such as e+e-qq¯ and e+e-qq¯qq¯; followed by a likelihood-based multivariate event selection. The preselection variables are formed by forcing each event into four jets using the Durham jet finder. Of the three possible jet associations with candidate W bosons, (12)(34), (13)(24) or (14)(23), the one giving a di-jet invariant mass closest to mW is selected. The preselection requires that there is no high-energy electron or muon with E > 30 GeV. Further preselection cuts are made on the properties of the jets, the invariant masses of the off-shell and on-shell W boson candidates, the Higgs boson candidate, the total visible energy and the missing transverse momentum. In addition, in order to reject Hbb¯ decays, the event is forced into a two-jet topology and flavour tagging is applied to the two jets. Events where at least one jet has a b-tag probability above 0.95 are rejected as part of the preselection. The cross sections and preselection efficiencies for the signal and main background processes are listed in Table Table20.20. After the preselection, the main backgrounds are e+e-qq¯νν¯, γe±qq¯qq¯ν and other Higgs decay modes, predominantly Hbb¯ and H → gg, where QCD radiation in the parton shower can lead to a four-jet topology.

Table 20
Summary of the HWWqq¯qq¯ event selection at s=1.4TeV, giving the raw cross sections, preselection efficiency, selection efficiency for a likelihood cut of ℒ > 0.35, and the expected numbers of events passing ...

A relative likelihood selection is used to classify all events passing the preselection cuts. Five event categories including the signal are considered. The relative likelihood of an event being signal is estimated as:

L=L(HWWqq¯qq¯)L(HWWqq¯qq¯)+L1+L2+L3+L4,

where Li represents the likelihood for four background categories: Hbb¯, H → gg, e+e-qq¯νν¯ and γe±qq¯qq¯ν. The absolute likelihood L for each event type is formed from normalised probability distributions Pi(xi) of the N likelihood discriminating variables xi for that event type. For example, the distribution of the reconstructed Higgs mass for all events passing the preselection is shown in Fig. Fig.15;15; it can be seen that good separation between signal and background is achievable. The discriminating variables are: the 2D distribution of reconstructed invariant masses mH and mW, the 2D distribution of minimal kt distances y23, y34, and the 2D distribution of b-tag probabilities when the event is forced into two jets. The use of 2D distributions accounts for the most significant correlations between the likelihood variables. The selection efficiencies and expected numbers of events for the signal dominated region, ℒ > 0.35, are listed in Table Table2020.

Fig. 15
Reconstructed Higgs invariant mass distributions for preselected HWWqq¯qq¯ events at s=1.4TeV, showing the signal and main backgrounds as stacked histograms. The distributions are normalised to an integrated ...

The expected precision on BR(H → WW) is extracted from a fit to the likelihood distribution. Given the non-negligible backgrounds from other Higgs decays, it is necessary to simultaneously fit the different components. A χ2 fit to the expected distribution is performed by scaling independently five components: the H → WW signal, the Hbb¯, Hcc¯ and H → gg backgrounds, and all other backgrounds (dominated by qq¯νν¯ and qq¯qq¯ν). The constraints on the Hbb¯, Hcc¯ and H → gg branching ratios, as described in Sect. 6.1, are implemented by modifying the χ2 function to include penalty terms:

χ2χ2+(sbb¯-1)2σbb¯2+(scc¯-1)2σcc¯2+(sgg-1)2σgg2+(sZZ-1)2σZZ2+(b-1)2σb2.

Here, for example, sgg is the amount by which the H → gg complement is scaled in the fit and σgg is the expected statistical error on BR(H → gg) from the analysis of Sect. 6.1. The expected uncertainties on the contributions from Hbb¯ and Hcc¯ are taken from the same analysis. The background from H → ZZ is assumed here to be known to 1% from other measurements of gHZZ2 and gHZZ2/gHWW2. The systematic uncertainty in the non-H background, denoted by b, is taken to be 1%. This has a small effect on the resulting uncertainty on the H → WW branching ratio, which is:

Δ[σ(Hνeν¯e)×BR(HWW)]σ(Hνeν¯e)×BR(HWW)=1.5%.

WWqq¯lν

As a second channel, the HWWqq¯lν decay is investigated [56]. The study is performed at s=1.4TeV using the CLIC_ILD detector model.

As a first step, isolated electrons or muons from W boson decay are identified. An efficiency of 93% is achieved for the identification of electrons and muons in signal events including the geometrical acceptance of the detector. Two jets are reconstructed from the remaining particles, excluding the isolated electron or muon. Flavour tagging information is obtained from the LcfiPlus package.

The following preselection cuts are imposed:

  • energy of the W candidate less than 590 GeV;
  • mass of the W candidate less than 230 GeV;
  • energy of the H candidate less than 310 GeV;
  • total missing energy of the event in the range between 670 GeV and 1.4 TeV.

Nearly all signal events pass this preselection, while more than 30% of the critical e+e-qq¯lν background events are rejected. The background processes are suppressed further using a BDT classifier with 19 input variables including the number of isolated leptons. The event selection is summarised in Table Table21.21. The resulting statistical precision for 1.5 ab-1 is:

Δ[σ(Hνeν¯e)×BR(HWW)]σ(Hνeν¯e)×BR(HWW)=1.3%.

The combined precision for HWWqq¯qq¯ and HWWqq¯lν decays at s=1.4TeV for an integrated luminosity of 1.5 ab-1 is 1.0%.

Table 21
Preselection and selection efficiencies for the signal and most important background processes in the HWWqq¯lν analysis. Numbers of events correspond to 1.5 ab-1 at s=1.4TeV

H → ZZ

The decay H → ZZ in e+e-Hνeν¯e events is studied using Z()qq¯ and Z(∗) → l+l- decays at s=1.4TeV using the CLIC_ILD detector model. The experimental signature is two jets, a pair of oppositely charged leptons and missing pT. The total invariant mass of all visible final-state particles is equal to the Higgs mass, while either the quarks or the charged lepton pair have a mass consistent with mZ. Due to the large background from H → WW, the ZZqq¯qq¯ final state is not considered here. The ZZ → l+l-l+l- signature is not expected to be competitive at CLIC due to the small number of expected events and is not further considered.

The analysis is performed in several steps. First, isolated electrons and muons with an impact parameter of less than 0.02 mm are searched for. Hadronic τ lepton decays are identified using the TauFinder algorithm described in Sect. 5.2.2, with the requirement pT > 10 GeV for the seed track and pT > 4 GeV for all other tracks within a search cone of 0.15 radian. In signal events, 87% of the electron or muon pairs and 37% of the tau lepton pairs are found, including the effect of the geometrical acceptance of the detector in the forward direction.

In events with exactly two identified leptons of the same flavour and opposite charge, two jets are reconstructed from the remaining particles. No other preselection cuts are applied. Flavour tagging information is obtained from the LcfiPlus package.

A BDT classifier is used to suppress the background processes using 17 input variables, including:

  • the invariant masses of the H, Z and Zcandidates;
  • the topology of the hadronic system: - log10(y34), - log10(y23) and - log10(y12);
  • the b-tag and c-tag probabilities for both jets;
  • the visible energy and the missing transverse momentum of the event;
  • the number of particles in the event.

The event selection is summarised in Table Table22.22. Only backgrounds from other Higgs decays pass the event selection, while all other background processes are fully rejected. The invariant mass distribution of the Higgs candidates in events with two isolated leptons after the full selection chain, including the BDT classifier, is shown in Fig. Fig.16.16. The resulting statistical uncertainty is:

Δ[σ(Hνeν¯e)×BR(HZZ)]σ(Hνeν¯e)×BR(HZZ)=5.6%.

Table 22
Preselection and selection efficiencies for the signal and the relevant background processes in the H → ZZ analysis. The numbers of events correspond to 1.5 ab-1 at s=1.4TeV. All background processes other than Higgs production are completely ...
Fig. 16
Reconstructed Higgs invariant mass distributions of HZZqq¯l+l- events at s=1.4TeV, showing the signal and main backgrounds as stacked histograms a after preselection, and b after the full event selection including a ...

H → γγ

The measurement of the H → γγ decay played a central role in the discovery of the Higgs boson at the LHC [1, 2]. In the SM, this decay is induced via loops of heavy charged particles, with dominant contributions from W bosons and t quarks. For BSM scenarios, other heavy charged particles can appear in the loops, modifying the expected effective H → γγ branching ratio. The sensitivity for the measurement of BR(H → γγ) at CLIC has been studied using the CLIC_SiD detector model for s=1.4TeV and an integrated luminosity of 1.5 ab-1. The SM branching ratio for mH = 126 GeV is 0.23% which results in approximately 840 signal events. The experimental signature for e+e-Hνeν¯e; H → γγ is two high pT photons with invariant mass m(γγ) consistent with mH, and missing momentum from the νeν¯e system. All relevant SM background processes with one or two photons in the final state have been considered. In addition to the photons from the hard interaction, the MC samples include additional ISR and FSR photons.

The following preselection cuts are applied to restrict the analysis to relevant events. At least two reconstructed photons each with energy Eγ > 15 GeV and pT > 10 GeV are required. The two highest energy photons passing these requirements are used to form the H candidate and the preselection requires an invariant mass consistent with mH, 115 GeV < m(γγ) < 140 GeV. The highest energy photon in the event is required to have pT > 40 GeV. In addition, to remove contributions from FSR, both photons are required to be isolated with no reconstructed particle with pT > 5 GeV within a cone of radius 500 mrad centred on the photon. Furthermore, the remaining reconstructed energy after excluding the Higgs candidate has to be less than 250 GeV. The cross sections and efficiencies of the preselection cuts for the signal and the main backgrounds are listed in Table Table23.23. At this stage in the event selection the background dominates.

Table 23
Signal and relevant background processes used in the H → γγ analysis. Additional photons from ISR and FSR are present in each sample. The cross sections for the backgrounds include cuts applied at generator level that are slightly looser than the preselection ...

To illustrate the photon reconstruction capabilities of the CLIC_SiD detector concept, the invariant mass of Higgs candidates in signal events after the preselection is shown in Fig. Fig.17.17. A fit to the distribution using a Gaussian function indicates a mass resolution in the signal sample of σ = 3.3 GeV.

Fig. 17
Reconstructed di-photon invariant mass distribution of preselected signal H → γγ events at s=1.4TeV. The distribution is normalised to an integrated luminosity of 1.5 ab-1. The statistical uncertainties correspond to the size of the simulated event sample. ...

The signal and background events are classified using a BDT. The 13 variables used to distinguish the signal from the backgrounds include:

  • the invariant mass of the Higgs candidate;
  • kinematic properties of the Higgs candidate;
  • kinematic properties of the two photons;
  • the angle between the two photons and the helicity angle of the Higgs candidate;
  • the remaining reconstructed energy excluding the Higgs candidate.

For the optimal BDT cut, the total signal selection efficiency is 40%, corresponding to 337 selected signal events in 1.5 ab-1. The event selection for the signal and the main backgrounds is summarised in Table Table23,23, leading to a statistical uncertainty of:

Δ[σ(Hνeν¯e)×BR(Hγγ)]σ(Hνeν¯e)×BR(Hγγ)=15%.

H → Zγ

As is the case for H → γγ, at lowest order, the SM decay H → Zγ is induced by loops of heavy charged particles. Contributions from BSM particles would lead to deviations from the SM expectation for BR(H → Zγ). For mH = 126 GeV, the decay H → Zγ is expected to have a branching ratio of BR(H → Zγ) = 0.16%. The potential to measure σ(e+e-Hνeν¯e)×BR(HZγ) at CLIC has been studied at s=1.4TeV with the CLIC_SiD detector model, where 585 H → Zγ events are expected in 1.5 ab-1 of data [57]. For the purpose of the event selection, only Zqq¯ and Z → l+l- (with l = e, μ) are useful, giving small event samples of 409 qq¯γ, 21 e+e-γ and 21 μ+μ-γ events from H → Zγ in 1.5 ab-1 at s=1.4TeV. A typical event display is shown in Fig. Fig.1818.

Fig. 18
Event display of a HZγqq¯γ event at s=1.4TeV in the CLIC_SiD detector. Both jets are visible. The photon creates a cluster in the central part of the electromagnetic calorimeter (blue)

The visible final states of the signal channels qq¯γ or l+l-γ are also produced in several background processes, some of which have much larger cross sections than the signal. In addition to background with photons from the hard process, e+e-qq¯ or e+e- → l+l- events with a FSR or ISR photon can mimic the signal.

The H → Zγ event selection requires at least one identified high-pT photon and either two electrons, muons or quarks consistent with a Z decay. The photon with the highest energy in the event is identified. Events are considered as either e+e-γ, μ+μ-γ or qq¯γ candidates. In the case where an e+e- or μ+μ- pair is found, photons nearly collinear with the lepton trajectories (within 0.3) are combined with the leptons under the assumption that these photons originate from bremsstrahlung. If neither an e+e- nor a μ+μ- pair is found, all reconstructed particles except for the photon of highest energy are clustered into two jets using a jet radius of R = 1.2. In all cases, the selected Z decay candidate and the highest energy photon are combined to form the H candidate.

In order to reduce the number of background process events, two selection steps are performed. First, preselection cuts are applied: the Higgs candidate daughter photon and jets, electrons, or muons are only accepted if they have an energy of E > 20 GeV and pT > 15 GeV. In the qq¯γ channel, only jets with at least 5 particles are considered in order to suppress hadronic τ decays. In addition, the reconstructed Z and H masses in the event are required to be consistent with a H → Zγ decay. The second step in the event selection is three BDT selections (one for each signal final state). The input variables are the properties of the reconstructed H, Z, and γ such as mass, energy, momentum, and polar angle, event shapes such as sphericity and aplanarity, as well as missing energy distributions and particle multiplicity distributions.

For the optimal BDT cuts, expected statistical significances of 2.2, 0.54 and 0.78 (in units of standard deviations) are found for the qq¯γ, e+e-γ and μ+μ-γ channels respectively. The signal selection efficiencies and contributions from the most important backgrounds are summarised in Table Table24.24. When the results from all three channels are combined, the expected statistical precision at s=1.4TeV for an integrated luminosity of 1.5 ab-1 is:

Δ[σ(Hνeν¯e)×BR(HZγ)]σ(Hνeν¯e)×BR(HZγ)=42%.

With electron polarisation the statistical precision can be increased, for example with - 80% electron polarisation, Δ[σ(e+e-Hνeν¯e)×BR(HZγ)]31%. Further gains are expected at higher centre-of-mass energies, as the Higgs production cross section at s=3TeV is 70% higher than at 1.4 TeV.

Table 24
Preselection and selection efficiencies for H → Zγ events in all three considered Z decay channels. The cross sections for the backgrounds include kinematic cuts applied at generator level. All numbers assume an integrated luminosity of 1.5 ab-1 at 1.4TeV ...

H → μ+μ-

The measurement of the rare H → μ+μ- decay is challenging due to the very low SM branching ratio of 2 × 10-4. In e+e-Hνeν¯e production, the signature for H → μ+μ- decay is a μ+μ- pair with invariant mass consistent with mH and missing momentum. The efficient rejection of background relies on the excellent detector momentum resolution, which directly influences the width of the reconstructed di-muon invariant mass peak. Signal and background events have been simulated at s=1.4TeV and 3 TeV using the CLIC_ILD and CLIC_SiD detector models respectively [58, 59]. In contrast with other studies presented in this paper, an electron beam polarisation of - 80% is assumed owing to the very small branching ratio for the H → μ+μ- decay. The two analyses were performed independently. They follow the same strategy but differ in some of the observables that are used in the event selection.

The most important background processes include μ+μ-νν¯ in the final state, as shown in Table Table2525 for 1.4 TeV and in Table Table2626 for 3 TeV. A significant fraction of these events are also produced from interactions involving beamstrahlung photons. Another important background is e+e- → e+e-μ+μ-, where both electrons are usually emitted at very low polar angles and thus might not be detected. Tagging of these low angle electrons in the very forward calorimeters – LumiCal and BeamCal – is essential to keep this background under control.

Table 25
The signal and main backgrounds in the H → μ+μ- analysis at s=1.4TeV with the corresponding cross sections. The numbers of selected events assume an integrated luminosity of 1.5 ab-1 and - 80% polarisation of the electron beam. Other processes, including ...
Table 26
The signal and most important background processes in the H → μ+μ- analysis at s=3TeV with the corresponding cross sections. The numbers of selected events assume an integrated luminosity of 2 ab-1 and - 80% polarisation of the electron beam. All other ...

The event selection requires two reconstructed, oppositely charged muons with a di-muon invariant mass within the relevant mass region of 105 - 145 GeV. Events with one or more detected high-energy electrons (E > 200 GeV at 1.4 TeV, E > 250 GeV at 3 TeV) in the very forward calorimeters are vetoed. This introduces the possibility of vetoing signal events if they coincide with Bhabha scattering events. The e+e- → e+e- cross section is sufficiently high that the probability of such a coincidence within 20 bunch crossings (10 ns) is about 7% in both analyses. The cuts on the minimum energy and the minimum polar angle for vetoing forward electrons need to be chosen carefully. e+e- → e+e-μ+μ- and e±γ → e±μ+μ- events need to be rejected efficiently while a low probability for coincidence with Bhabha scattering events needs to be maintained.

The 3 TeV analysis includes some additional preselection cuts to remove phase space regions that do not include any signal events. These cuts reject events that contain a reconstructed non-muon object with an energy greater than 100 GeV; in addition, events containing electrons in the central region of the detector with an energy above 20 GeV are also rejected. The sum of the transverse momenta of the two muons, pT-) + pT+), is required to be above 50 GeV and the transverse momentum of the di-muon system should be above 25 GeV.

The final event selection uses a BDT classifier using various kinematic variables, excluding the invariant mass of the di-muon system. The 1.4 TeV analysis uses the visible energy of the event after removal of the di-muon system Evis, the transverse momentum of the di-muon system pT(μμ), the sum of the transverse momenta of the two muons pT-) + pT+), the polar angle of the di-muon system θμμ, the boost of the di-muon system βμμ, and the cosine of the helicity angle cosθ. The 3 TeV analysis uses the energy of the hardest non-muon object instead of the total visible energy and also includes the energy, transverse momentum, polar angle and azimuthal angle of both individual muons. This event selection reduces background from four-fermion processes by several orders of magnitude, while maintaining an overall signal selection efficiency of ϵ = 30.5% and ϵ = 26.3% at 1.4 and 3 TeV respectively.

The number of signal events is extracted from the reconstructed invariant mass distribution after the event selection, as shown in Fig. Fig.19.19. Using a large MC sample, the signal and background shapes are extracted. The signal is described by a Gaussian distribution with asymmetric exponential tails. The combined background is parameterised as the sum of an exponential and a constant function. To assess the expected statistical precision, a large number of trial samples are generated from the expected reconstructed mass distributions of signal and background and are then fitted to the signal and background components. For P(e-) =  - 80%, the expected relative uncertainty on the σ(e+e-Hνeν¯e)×BR(Hμ+μ-) is 27%, corresponding to a significance of 3.7, at 1.4 TeV, and 19%, corresponding to a significance of 5.2, at 3 TeV. The corresponding uncertainties for unpolarised beams are:

Δ[σ(Hνeν¯e)×BR(Hμ+μ-)]σ(Hνeν¯e)×BR(Hμ+μ-)=38%at1.4TeV,Δ[σ(Hνeν¯e)×BR(Hμ+μ-)]σ(Hνeν¯e)×BR(Hμ+μ-)=25%at3TeV.

Fig. 19
Reconstructed di-muon invariant mass distribution of selected H → μ+μ- events at s=3TeV. The simulated data are shown as dots while the solid line represents the fit function described in the text. The dotted line shows the background contribution of ...

ZZ-fusion

Higgs boson production through the t-channel fusion of two Z bosons, e+e- → He+e-, is analogous to the WW-fusion process but gives access to gHZZ and gHbb using a complementary technique. At s=1.4TeV, ZZ-fusion is the sub-leading Higgs production process, with a cross section of around 25 fb, which is 10% of that for the WW-fusion process. The potential for the measurement of the ZZ-fusion process has been investigated at s=1.4TeV using the CLIC_ILD detector.

The characteristic signature of the ZZ-fusion process is two scattered beam electrons reconstructed in the forward regions of the detector, plus the Higgs boson decay products. Here, the scattered beam electrons are required to be fully reconstructed, and the final state Hbb¯ is considered.

Events are clustered into a four-jet topology using a kt exclusive clustering algorithm with R = 1.0. For a well-reconstructed signal event, two of the resulting ‘jets’ are expected to be the reconstructed electrons, and the remaining two jets originate from the Higgs decay to bb¯. The event selection requires two oppositely-charged electron candidates, separated by |Δη| > 1, each with E > 100 GeV. This preselection preserves 27% of the e+e-He+e-bb¯e+e- signal (3.6 fb), with the lost events almost entirely due to the scattered electrons falling outside the detector acceptance, as shown in Fig. Fig.20.20. After the preselection, the SM background consists mainly of events that have two real electrons and a qq¯ pair, either from the continuum or from the decay of Z bosons. Although the preselection suppresses 98% of the e+e-qq¯e+e- background, the accepted cross section is 48 fb, which is thirteen times larger than that for the remaining signal. A further requirement that one of the two jets associated with the Higgs decay has a b-tag value  > 0.4 preserves 80% of the remaining signal and rejects 80% of the remaining background.

Fig. 20
Generated electron pseudorapidity (η=-lntanθ2) distributions for e+e- → He+e- events at s=1.4 and 3 TeV. The distributions are normalised to 1.5 and 2 ab-1 respectively. The vertical arrows show the detector acceptance

A relative likelihood classifier 1, which treats ZZ-fusion events with Hbb¯ as signal and H → WW and H → ZZ as background, is used to reduce contributions from other Higgs decays. Seven variables are used to construct the likelihood: the jet clustering variable y45; the invariant mass of the two jets associated with the Higgs decay; the visible mass of the event with the scattered beam electrons removed; the higher of the b-tag values of the two jets associated with the Higgs decay; the c-tag value corresponding to the same jet; and the b-c-separation returned by the tagger, for both Higgs decay jets. Requiring a high signal likelihood, 1 > 0.8, reduces the Hbb¯ signal to 3000 events but leaves only 90 events from other Higgs decays, while also reducing the non-Higgs backgrounds to 4700 events.

Finally, to separate the signal from all backgrounds, a further relative likelihood classifier 2 is constructed using four variables that provide separation power between signal and background: the opening between the reconstructed electrons ΔR; the recoil mass of the event determined from the momenta of the reconstructed electrons, mrec; the jet clustering variable y34; and the invariant mass of the two jets associated with the Higgs decay.

The resulting likelihood is shown in Fig. Fig.2121 and gives good separation between signal and background. The likelihood distribution is fitted by signal and background components (where the normalisation is allowed to vary), giving:

Δ[σ(He+e-)×BR(Hbb¯)]σ(He+e-)×BR(Hbb¯)=1.8%

for 1.5 ab-1 at s=1.4TeV.

Fig. 21
Likelihood distributions for Hbb¯ events in the ZZ-fusion analysis at s=1.4TeV, shown for the signal and main background. The distributions are normalised to an integrated luminosity of 1.5 ab-1

Top Yukawa coupling

At an e+e- collider the top Yukawa coupling, yt, can be determined from the production rate in the process where a Higgs boson is produced in association with a top quark pair, e+e-tt¯H. The top quarks decay almost exclusively by t → bW. The signal event topology thus depends on the nature of the W and Higgs boson decays. Here Hbb¯ decays have been studied for two tt¯H decay channels at s=1.4TeV using the CLIC_SiD detector model [60, 61]:

  • the fully-hadronic channel (where both W bosons decay hadronically), giving a tt¯H final state of eight jets, including four b jets;
  • the semi-leptonic channel (where one W boson decays leptonically), giving a tt¯H final state of six jets (four b jets), one lepton and one neutrino.

The two channels are distinguished by first searching for isolated leptons (muons and electrons with an energy of at least 15 GeV and tau candidates from TauFinder containing a track with pT > 10 GeV). If zero leptons are found, the event is classified as fully-hadronic. If one lepton is found, the event is classified as semi-leptonic. Events in which more than one lepton is found are not analysed further. The kt algorithm is used to cluster the particles of each event into a specific number of jets, and remove particles arising from beam-beam interactions that are closer to the beam axis than to a hard jet as described in Sect. 4.2. Events classified as fully-hadronic are clustered into eight jets. In semi-leptonic events, the lepton is removed and the remaining particles are clustered into six jets. A semi-leptonic event is shown in Fig. Fig.22.22. The particles not clustered into jets by the kt algorithm are removed from the event and the remaining particles are then re-clustered using the e+e- Durham algorithm in LcfiPlus, which performs flavour tagging for each jet, and prevents particles from displaced vertices being split between two or more jets. The jets are combined to form candidate primary particles in such a way so as to minimise a χ2 function expressing the consistency of the reconstructed di- and tri-jet invariant masses with the tt¯(Hbb¯) hypothesis. For example, in the case of the semi-leptonic channel, the jet assignment with the minimum of:

χ2=(mij-mW)2σW2+(mijk-mt)2σt2+(mlm-mH)2σH2,

gives the W, top and Higgs candidates, where mij is the invariant mass of the jet pair used to reconstruct the W candidate, mijk is the invariant mass of the three jets used to reconstruct the top quark candidate and mlm is the invariant mass of the jet pair used to reconstruct the Higgs candidate. The expected invariant mass resolutions σW,t,H were estimated from combinations of two or three reconstructed jets matched to W, top and Higgs particles on generator level.

Fig. 22
Event display of a tt¯Hbb¯bb¯qq¯τ-ν¯τ event at s=1.4TeV in the CLIC_SiD detector. The tau lepton decays hadronically

Having forced each event into one of the two signal-like topologies, multivariate BDT classifiers (one for fully-hadronic events and one for semi-leptonic events) are used to separate signal and background. The discriminating variables include: kinematic quantities such as the reconstructed Higgs mass, the visible energy in the jets and the missing pT; angular variables such as the angles between the Higgs decay products in the rest frame of the Higgs candidate with respect to its flight direction and the angle between the momenta of the top and Higgs candidates; event variables such as thrust, sphericity and the number of particles in the event; and flavour tag variables for the four most likely b-jets. As an example, the BDT response distributions for the fully-hadronic channel are shown in Fig. Fig.23.23. The selection is chosen to maximise the signal significance. The expected numbers of selected events for 1.5 ab-1 of s=1.4TeV data are listed in Table Table27.27. The contributions from other investigated background processes were found to be negligible. The tt¯H cross section can be measured with an accuracy of 12% in the semi-leptonic channel and 11% in the hadronic channel. The combined precision of the two channels is 8%.

Fig. 23
BDT classifier distributions for fully-hadronic tt¯H events at s=1.4TeV, shown for the tt¯H signal and main backgrounds. The distributions are normalised to an integrated luminosity of 1.5 ab-1. The vertical arrow shows the value of the ...
Table 27
Expected numbers of signal and background events in the fully-hadronic (HAD) and semi-leptonic (SL) channels for 1.5 ab-1 at s=1.4TeV. The columns show the total numbers of events before selection and the numbers of events passing the fully-hadronic and ...

To translate the measurement of the tt¯H cross section into a measurement of the top Yukawa coupling, a correction is applied to take into account the contribution from the Higgsstrahlung diagram, where the Higgs boson is radiated off the intermediate Z boson in e+e-tt¯ [62, 63]. To evaluate the small degradation in sensitivity, the WHIZARD program is used to calculate the cross section for the inclusive process e+e-tt¯H as a function of the value of the top Yukawa coupling. The factor required to translate the measured cross section uncertainty into a coupling uncertainty is determined from the slope of the cross section at the SM value of the top Yukawa coupling, and is found to be:

Δytyt=0.53Δσσ,

which is slightly larger than the factor of 0.50 expected without the Higgsstrahlung diagram. Thus, the expected precision on the top Yukawa coupling is:

Δytyt=4.2%,

for 1.5 ab-1 of data at s=1.4TeV without beam polarisation. This value is expected to improve to about 4.0% for the same amount of data collected using the P(e-) =  - 80% polarisation configuration [64]. Since the cross section for the tt¯H cross section falls with increasing s (see Fig. Fig.3),3), the precision with 2 ab-1 at 3 TeV is not expected to be better than the result presented here.

Double Higgs production

In e+e- collisions at high energy, double Higgs production, e+e-HHνeν¯e, can occur through the processes shown in Fig. Fig.24.24. Despite the small cross section (0.15 and 0.59 fb for CLIC operated at s=1.4 and 3 TeV, respectively), measurements of the double Higgs production rate can be used to extract the Higgs boson trilinear self-coupling parameter λ, that determines the shape of the fundamental Higgs potential. BSM physics scenarios can introduce deviations of λ from its SM value of up to tens of percent [65]. The physics potential for the measurement of this coupling has been studied using the CLIC_ILD detector model for 1.5 ab-1 of data at s=1.4TeV and for 2 ab-1 of data at s=3TeV. The process e+e- → HHe+e- has not been included as its cross section is about an order of magnitude smaller compared to e+e-HHνeν¯e.

Fig. 24
Feynman diagrams of leading-order processes that produce two Higgs bosons and missing energy at CLIC at s=1.4TeV and 3 TeV. The diagram (a) is sensitive to the trilinear Higgs self-coupling λ. The diagram (b) is sensitive to the quartic coupling gHHWW ...

Two signatures for e+e-HHνeν¯e production are considered in the following: HHbb¯bb¯ and HHbb¯WWbb¯qq¯qq¯. All events without isolated leptons are considered for the analysis. These events are clustered into four jets using the kt algorithm. Flavour tagging information is obtained from the LcfiPlus package. Events where the sum of the b-tag values of the four jets is smaller than 2.3 and the hadronic system fulfills the requirement - log10(y34) < 3.7(3.6) at 1.4 TeV (3 TeV) are considered as bb¯WW candidates, while all other events are considered as bb¯bb¯ candidates. The following steps of the analysis are performed separately for the two final states.

At 1.4 TeV, a cut on the sum of the four b-tag values of at least 1.5 is imposed for bb¯bb¯ candidate events. Those events with a sum of the four b-tag values less than 2.3 are required to have a sum of the jet energies of at least 150 GeV and a second highest jet transverse momentum of at least 25 GeV. A cut on the sum of the four b-tag values of at least 2.3 is imposed for all events at 3 TeV. The jets are grouped into two Higgs boson candidates by minimising |mijmkl|, where mij and mkl are the invariant masses of the jet pairs used to reconstruct the Higgs candidates. For events passing the preselection cuts, at both energies BDT classifiers with the same 10 input variables are used to suppress the backgrounds further.

For the bb¯WW final state, the events are re-clustered into six jets. These jets are then grouped into W and H candidates by minimising:

χ2=(mij-mH)2σHbb¯2+(mklmn-mH)2σHWW2+(mkl-mW)2σW2,

where mij and mklmn are the jet combinations used to reconstruct the Higgs candidates, mkl is the invariant mass of the jet pair used to reconstruct the W candidate and σHbb¯, σH→WW, σW are the estimated invariant mass resolutions for the reconstruction of Hbb¯, H → WW and W decays. Events with an invariant mass of the two H boson candidates above 150 GeV are considered further. At 3 TeV a highest b-tag value of at least 0.7 is required while at 1.4 TeV the second highest b-tag values has to be larger than 0.2 and the visible transverse momentum has to be larger than 30 GeV. After this preselection, BDT classifiers using 32 input variables are used to suppress the backgrounds further.

The event selections for both studies at 1.4 and 3 TeV are summarised in Tables Tables2828 and and29,29, respectively. Combining the expected precisions on the cross sections for both signatures leads to:

Δ[σ(HHνeν¯e)]σ(HHνeν¯e)=44%at1.4TeV,Δ[σ(HHνeν¯e)]σ(HHνeν¯e)=20%at3TeV.

Table 28
Preselection and selection efficiencies for the double Higgs signal and most important background processes in both considered decay channels at s=1.4TeV. The numbers of events correspond to 1.5 ab-1. Contributions from all other backgrounds are found ...
Table 29
Preselection and selection efficiencies for the double Higgs signal and most important background processes in both considered decay channels at s=3TeV. The numbers of events correspond to 2 ab-1. Contributions from all other backgrounds are found to ...

The double Higgs production cross section is sensitive to the trilinear Higgs self-coupling λ. Since diagrams not involving λ also contribute to the e+e-HHνeν¯e process, their effect must be taken into account. The relation between the relative uncertainty on the cross section and the relative uncertainty of the Higgs trilinear coupling can be approximated as:

Δλλκ·Δ[σ(HHνeν¯e)]σ(HHνeν¯e).

The value of κ can be determined from the WHIZARD generator by parameterising the e+e-HHνeν¯e cross section as a function of the input value for λ, as indicated in Fig. Fig.25.25. The fact that the slope is negative indicates that the main dependence on λ enters through interference with other SM diagrams. The value of κ is determined from the derivative of the cross section dependence as a function of λ, evaluated at its SM value, giving κ = 1.22 and κ = 1.47 at 1.4 TeV and 3 TeV, respectively. However, this method does not account for the possibility that the event selection might preferentially favour some diagrams over others, and hence change the analysis sensitivity to λ.

Fig. 25
Cross section for the e+e-HHνeν¯e process as a function of the ratio λ/λSM at s=1.4 and 3 TeV

In the case of zero beam polarisation, the combined cross sections for double Higgs production give:

Δλ/λ=54%ats=1.4TeV,Δλ/λ=29%ats=3TeV.

Because the process involving the trilinear Higgs coupling involves t-channel WW-fusion, it can be enhanced by operating with polarised beams. For the case of P(e-) =  - 80%, this yields:

Δλ/λ=40%ats=1.4TeV,Δλ/λ=22%ats=3TeV.

The statistical precision on λ improves to 26% for unpolarised beams and to 19% for P(e - ) =  - 80% when combining both energy stages. These results will be improved further using template fits to the BDT output distributions as the different diagrams contributing to double Higgs production lead to different event topologies.

Higgs mass

At a centre-of-mass energy of s=350GeV, the Higgs boson mass can be measured in the e+e- → ZH process. The Higgs boson mass can be extracted from the four-momentum recoiling against in Z boson using Z → e+e- or Z → μ+μ- events as described in Sect. 5. Due to the small branching ratios for leptonic Z boson decay channels and the impact of the CLIC beamstrahlung spectrum, the achievable precision is limited to 110 MeV.

In a different approach, the Higgs mass is reconstructed from the measured four-vectors of its decay products. The best precision is expected using Hbb¯ decays in e+e-Hνeν¯e events at high energy. For this purpose, the analysis described in Sect. 6.1 has been modified. After the preselection, a single BDT is used at each energy to select Hbb¯ decays. In contrast to the coupling measurement, the flavour tagging information is included in the BDT classifier.

The invariant mass distribution for selected events at s=1.4TeV is shown in Fig. Fig.26.26. The Higgs mass is extracted in the range 105GeV<mbb¯<145GeV where good purity of the signal channel is achieved. At the nominal Z boson mass, a second peak from e+e-Zνeν¯e;Zbb¯ events is visible. These events can be used to calibrate the jet energy scale for the precision measurement of the Higgs boson mass.

Fig. 26
Reconstructed di-jet invariant mass distribution of selected Hbb¯ events at s=1.4TeV, showing the signal and backgrounds as stacked histograms. The distributions are normalised to an integrated luminosity of 1.5 ab-1

A template fit using e+e-Hνeν¯e;Hbb¯ event samples generated using slightly shifted values for the Higgs mass parameter is performed. The Higgs mass and production cross section are extracted simultaneously. The following statistical precisions on the Higgs mass are achieved:

Δ(mH)=47MeVat1.4TeV,Δ(mH)=44MeVat3TeV.

A combination of both energy stages would lead to a precision of 32 MeV.

Systematic uncertainties

The complete Higgs physics potential of a CLIC collider implemented in three energy stages is described in this paper. The expected statistical uncertainties given in the previous sections do not include potential sources of systematic uncertainty. The obtained results therefore illustrate the level of precision desirable for the control of systematic effects. This is crucial input for the choice of detector technologies and the development of calibration procedures in the coming years.

A comprehensive study of systematic uncertainties requires more knowledge on the technical implementation of the detector. This is beyond the scope of this paper. At this stage, the impact of potentially relevant sources of systematic uncertainty is discussed. The measurements of σ(Hνeν¯e)×BR(Hbb¯) and the Higgs mass at s=3 TeV, described in Sects. 6.1 and 10, are used as examples. These measurements are the most challenging test cases for many systematic effects due to the very small expected statistical uncertainties of 0.3% and 44 MeV, respectively. In addition, the experimental conditions are most challenging at 3 TeV.

The impact of theoretical uncertainties on the Higgs branching fractions is discussed in Sect. 12 in the context of a combined fit.

  • Luminosity spectrum: A good knowledge of the luminosity spectrum is mandatory for precision Higgs physics at CLIC. The reconstruction of the CLIC luminosity spectrum from Bhabha scattering events is described in [66]. A model of the CLIC luminosity spectrum with 19 free parameters is assumed. The expected uncertainties of these parameters and their correlations are propagated to the measurement of σ(Hνeν¯e)×BR(Hbb¯) and lead to a systematic uncertainty of 0.15%. The luminosity spectrum affects the event rate more than the observed invariant mass of the two jets. Concerning the Higgs mass extraction, the luminosity spectrum is not expected to represent a dominant source of systematic uncertainty since the cross section is a free parameter in the template fit.
  • Total luminosity: The expected statistical precision of the σ(Hνeν¯e)×BR(Hbb¯) measurement indicates the desired precision for the knowledge of the total luminosity. It is expected that an accuracy of a few permille can be achieved using the luminometer envisaged for CLIC [67, 68].
  • Beam polarisation: The knowledge of the beam polarisation at the interaction point is most important for the measurement of WW-fusion events at high energy. The beam polarisation can be controlled to a level of 0.2% using single W, Z and γ events with missing energy [69]. The resulting systematic uncertainty on σ(Hνeν¯e)×BR(Hbb¯) is 0.1%. For the Higgs mass measurement, the effect of the estimated beam polarisation uncertainty is negligible.
  • Jet energy scale: The measurement of the Higgs boson mass using Hbb¯ decays requires a precise knowledge of the energy scale correction for b-jets. An uncertainty on the jet energy scale of 3.5 × 10-4 leads to a systematic uncertainty on the Higgs mass similar to the statistical error at 3 TeV. The same jet energy scale uncertainty would have negligible impact on σ(Hνeν¯e)×BR(Hbb¯). A suitable process for the calibration is e+e-Zνeν¯e;Zbb¯ which is kinematically similar to Higgs production in WW-fusion. σ(Zνeν¯e)×BR(Zbb¯)=276 fb leads to an expected number of events for calibration which is slightly larger than the signal event sample. To improve the precision further, additional high-statistics Z boson samples would be needed. Generator-level studies show that e±γZe±;Zbb¯ with a cross section about one order of magnitude larger compared to the signal process is a promising channel for this purpose.
  • Flavour tagging: Several of the precision measurements discussed in this paper rely on b-tagging information. The calibration of the flavour tagging at CLIC is a topic for future study. To illustrate the impact of a non-perfect understanding of the mistag rate for charm and light quark jets, an ad hoc variation of the b-tag distributions for jets in background events is performed. Even after the BDT selection, the background contains only very few b-jets in the σ(Hνeν¯e)×BR(Hbb¯) analysis. First, the b-tag distributions for both jets were decreased (increased) by 0.5% using event reweighting for values below (above) the median keeping the overall number of background events constant. The opposite variation is applied in a second step. These variations lead to a ±0.25% change of the result. As the flavour tagging efficiency mostly affects the event rate, it is not expected to be a dominant source of systematic uncertainty for the Higgs mass measurement.

In summary, it seems possible to control the systematic uncertainties discussed above with similar or better precision compared to the statistical uncertainty for the measurement of σ(Hνeν¯e)×BR(Hbb¯). An excellent understanding of the b-jet energy scale is necessary for a competitive Higgs mass measurement at CLIC.

Many of the analyses described in this paper, especially where harmonization is relevant, will require a careful tuning of the Monte Carlo models using other high-precision processes. Such an investigation is beyond the scope of this first study of Higgs physics at CLIC presented here.

Combined fits

The results discussed in the preceding sections are summarised in Tables Tables3030 and and31.31. From the σ and σ × BR measurements given in the tables the Higgs coupling parameters and total width are extracted by a global fit as described below. Here, a - 80% electron polarisation is assumed for the 1.4 and the 3 TeV stages. The increase in cross section is taken into account by multiplying the event rates with a factor of 1.8 for all WW-fusion measurements (see Table Table3),3), resulting in a reduction of the uncertainties by a factor of 1.8. This approach is conservative since it assumes that all backgrounds including those from s-channel processes, which do not receive the same enhancement by polarisation, scale with the same factor.

Table 30
Summary of the precisions obtainable for the Higgs observables in the first stage of CLIC for an integrated luminosity of 500 fb-1 at s=350GeV, assuming unpolarised beams. For the branching ratios, the measurement precision refers to the expected statistical ...
Table 31
Summary of the precisions obtainable for the Higgs observables in the higher-energy CLIC stages for integrated luminosities of 1.5 ab-1 at s=1.4TeV, and 2.0 ab-1 at s=3TeV. In both cases unpolarised beams have been assumed. For gHtt, the 3 TeV case has ...

A few of the observables listed in Table Table3131 were studied only at s=1.4TeV, but not at s=3TeV. In cases where those observables have a significant impact on the combined fits described in this section, the precisions obtained at s=1.4TeV were extrapolated to s=3TeV. The extrapolation is based on the number of signal events within the detector acceptance at 1.4 and 3 TeV. It is assumed that the background processes scale in the same way with s as the signal events. However, in fact the signal Higgs bosons are produced in vector boson fusion which increases with increasing s, while several backgrounds are dominated by s-channel diagrams which decrease with increasing s.

Since the physical observables (σ or σ × BR) typically depend on several coupling parameters and on the total width, these parameters are extracted with a combined fit of all measurements. To provide a first indication of the overall impact of the CLIC physics programme, simple fits considering only the statistical uncertainties of the measurements are performed. Two types of fits are used: A model-independent fit making minimal theoretical assumptions, and a model-dependent fit following the strategies used for the interpretation of LHC Higgs results.

Both fits are based on a χ2 minimisation using the Minuit package [70]. The measurements which serve as input to the fit, presented in detail in the preceding sections, are either a total cross section σ in the case of the measurement of e+e- → ZH via the recoil mass technique, or a cross section  ×  branching ratio σ × BR for specific Higgs production modes and decays. To obtain the expected sensitivity for CLIC it is assumed that for all measurements the value expected in the SM has been measured, so only the statistical uncertainties of each measurement are used in the χ2 calculation. In the absence of correlations, the contribution of a single measurement is given by

χi2=(Ci/CiSM-1)2ΔFi2,

where Ci is the fitted value of the relevant combination of relevant Higgs couplings (and total width) describing the particular measurement, CiSM is the SM expectation, and ΔFi is the statistical uncertainty of the measurement of the considered process. Since this simplified description does not allow the accurate treatment of correlations between measurements, nor the inclusion of correlated theory systematics in the model-dependent fit, the global χ2 of the fit is constructed from the covariance matrix of all measurements. It is given by

χ2ζTV-1ζ

where V is the covariance matrix and ζ is the vector of deviations of fitted values of the relevant combination of Higgs couplings and total width describing the particular measurement deviation from the SM expectation as introduced above, ζi=Ci/CiSM-1.

The Ci’s depend on the particular measurements and on the type of fit (model-independent or model-dependent), given in detail below. In the absence of systematic uncertainties, the diagonal elements of V are given by the statistical uncertainty of the measurement,

Vii=ΔFi2,

while the off-diagonal elements represent the correlations between measurements. In the fit, correlations are taken into account in cases where they are expected to be large. This applies to the measurements of σ × BR for Hbb¯,cc¯,gg in Higgsstrahlung and WW-fusion events at 350 GeV and in WW-fusion events only at 1.4 and 3 TeV, which are extracted in a combined fitting procedure at each energy. These measurements show correlation coefficients with absolute values as large as 0.32.

In signal channels with substantial contaminations from other Higgs decays, penalty terms were added to the χ2 to take into account the normalisation of the other channels. These additional uncertainties, which are also of a statistical nature, are derived from the statistical uncertainties of the respective Higgs final state analysis, taking the level of contamination into account. The channels where this results in non-negligible effects are the H → WW analyses at all energies, in particular in the all-hadronic decay modes, with corrections to the statistical uncertainties as large as 8% at 350 GeV.

Model-independent fit

The model-independent fit uses the zero-width approximation to describe the individual measurements in terms of Higgs couplings and the total width, ΓH. Here, the total cross section of e+e- → ZH depends on:

CZH=gHZZ2,

while for specific final states such as e+e- → ZH; Hbb¯ and e+e-Hνeν¯e; Hbb¯:

CZH,Hbb¯=gHZZ2gHbb2ΓH

and:

CHνeν¯e,Hbb¯=gHWW2gHbb2ΓH,

respectively.

The fit is performed with 11 free parameters: gHZZ, gHWW, gHbb, gHcc, gHττ, gHμμ, gHtt and ΓH, as well as the three effective couplings gHgg, gHγγ and gHZγ. The latter three parameters are treated in the same way as the physical Higgs couplings in the fit.

The fit is performed in three stages, taking the statistical uncertainties obtainable from CLIC at the three considered energy stages (350 GeV, 1.4, 3 TeV) successively into account. Each new stage also includes all measurements of the previous stages. Table Table3232 summarises the results. They are graphically illustrated in Fig. Fig.27.27. Since the model-independence of the analysis hinges on the absolute measurement of σ(ZH) at 350 GeV, which provides the coupling gHZZ, the precision of all other couplings is ultimately limited by this uncertainty.

Table 32
Results of the model-independent fit. Values marked “−” can not be measured with sufficient precision at the given energy. For gHtt, the 3 TeV case has not yet been studied, but is not expected to result in substantial improvement ...
Fig. 27
Illustration of the precision of the Higgs couplings of the three-stage CLIC programme determined in a model-independent fit without systematic or theoretical uncertainties. The dotted lines show the relative precisions of 1 and 5%

Model-dependent fit

For the model-dependent fit, it is assumed that the Higgs decay properties can be described by ten independent parameters κHZZ, κHWW, κHbb, κHcc, κHττ, κHμμ, κHtt, κHgg, κHγγ and κHZγ. These factors are defined by the ratio of the Higgs partial width divided by the partial width expected in the Standard Model as:

κi2=Γi/ΓiSM.

In this scenario, the total width is given by the sum of the ten partial widths considered, which is equivalent to assuming no non-Standard-Model Higgs decays such as decays into new invisible particles. The ratio of the total width to its SM value is thus given by:

ΓH,mdΓHSM=iκi2BRi,
1

where BRi is the SM branching fraction for the respective final state and the subscript “md” stands for “model-dependent”. To obtain these branching fractions, a fixed value for the Higgs mass has to be imposed. For the purpose of this study, 126 GeV is assumed. The branching ratios are taken from the LHC Higgs cross section working group [22]. To exclude effects from numerical rounding errors, the total sum of BR’s is normalised to unity.

With these definitions, the Ci’s in the χ2 take the following forms: for the total e+e- → ZH cross section:

CZH=κHZZ2;

while for specific final states such as e+e- → ZH; Hbb¯ and e+e-Hνeν¯e; Hbb¯:

CZH,Hbb¯=κHZZ2κHbb2ΓH,md/ΓHSM

and:

CHνeν¯e,Hbb¯=κHWW2κHbb2ΓH,md/ΓHSM,

respectively.

Since at the first energy stage of CLIC no significant measurements of the H → μ+μ-, H → γγ and H → Zγ decays are possible, the fit is reduced to six free parameters (the coupling to top is also not constrained, but this is without effect on the total width) by setting H → μ+μ-, H → γγ and H → Zγ to zero. These branching ratios are much smaller than the derived uncertainty on the total width.

Two versions of the model-dependent fit are performed, one ignoring theoretical uncertainties to illustrate the full potential of the constrained fit, and one taking the present theoretical uncertainties of the branching fractions into account [22]. To avoid systematic biases in the fit results, the uncertainties are symmetrised, preserving the overall size of the uncertainties. Theoretical uncertainties on the production are assumed to be substantially smaller than in the decay, and are ignored in the present study. Depending on the concrete Higgs decay, multiple measurements may enter in the fit, originating from different centre-of-mass energies, different production channels or different signal final states. To account for this, the theoretical uncertainties are treated as fully correlated for each given Higgs decay.

As in the model-independent case the fit is performed in three stages, taking the statistical errors of CLIC at the three considered energy stages (350 GeV, 1.4, 3 TeV) successively into account. Each new stage also includes all measurements of the previous stages. The total width is not a free parameter of the fit. Instead, its uncertainty, based on the assumption given in Eq. 1, is calculated from the fit results, taking the full correlation of all parameters into account. Table Table3333 summarises the results of the fit without taking theoretical uncertainties into account, and Fig. Fig.2828 illustrates the evolution of the precision over the full CLIC programme. Table Table3434 summarises the results of the model-dependent fit with theoretical uncertainties of the branching fractions.

Table 33
Results of the model-dependent fit without theoretical uncertainties. Values marked “−” can not be measured with sufficient precision at the given energy. For gHtt, the 3 TeV case has not yet been studied, but is not expected to ...
Fig. 28
Illustration of the precision of the Higgs couplings of the three-stage CLIC programme determined in a model-dependent fit without systematic or theoretical uncertainties. The dotted lines show the relative precisions of 0.5 and 2.5%
Table 34
Results of the model-dependent fit with the current theoretical uncertainties on the decay branching fractions. Values marked “−” can not be measured with sufficient precision at the given energy. For gHtt, the 3 TeV case has not ...

Discussion of fit results

The full Higgs physics programme of CLIC, interpreted with a combined fit of the couplings to fermions and gauge bosons as well as the total width, and combined with the measurement of the self-coupling, will provide a comprehensive picture of the properties of this recently discovered particle. Figure Figure2929 illustrates the expected uncertainties of the various couplings determined in the model-independent fit as well as the self-coupling as a function of the particle mass. Combined with the quasi model-independent measurement of the total width with a precision of 3.5%, this illustrates the power of the three-stage CLIC programme. Each of the stages contributes significantly to the total precision, with the first stage at 350 GeV providing the model-independent “anchor” of the coupling to the Z boson, as well as a first measurement of the total width and coupling measurements to most fermions and bosons. The higher-energy stages add direct measurements of the coupling to top quarks, to muons and photons as well as overall improvements of the branching ratio measurements and with that of the total width and all couplings except the one to the Z already measured in the first stage. They also provide a measurement of the self-coupling of the Higgs boson. In a model-dependent analysis, the improvement with increasing energy is even more significant than in the model-independent fit, since the overall limit of all couplings imposed by the model-independent measurement of the ZH recoil process is removed.

Fig. 29
Illustration of the precision of the model-independent Higgs couplings and of the self-coupling as a function of particle mass. The line shows the SM prediction that the Higgs coupling of each particle is proportional to its mass

Summary and conclusions

A detailed study of the Higgs physics reach of CLIC has been presented in the context of CLIC operating in three energy stages, s=350GeV, 1.4 and 3 TeV. The initial stage of operation, 500 fb-1 at s=350GeV, allows the study of Higgs production from both the e+e- → ZH and the WW-fusion process. These data yield precise model-independent measurements of the Higgs boson couplings, in particular Δ(gHZZ) = 0.8%, Δ(gHWW) = 1.4% and Δ(gHbb) = 3.0%. In addition, the branching ratio to invisible decay modes is constrained to Γinvis/ΓH < 0.01 at 90% C.L. and the total Higgs width is measured to Δ(ΓH) = 6.7%. Operation of CLIC at s>1TeV provides high-statistics samples of Higgs bosons produced through the WW-fusion process and give access to rarer processes such as e+e-tt¯H and e+e-HHνeν¯e. Studies of these rare processes provide measurements of the top Yukawa coupling to 4.2% and the Higgs boson self-coupling to about 20%. Furthermore, the full data sample leads to very strong constraints on the Higgs couplings to vector bosons and fermions. In a model-independent treatment, many of the accessible couplings are measured to better than 2%, and the model-dependent κ parameters are determined with a precision of between 0.1 and 1%.

Acknowledgements

This work benefited from services provided by the ILC Virtual Organisation, supported by the national resource providers of the EGI Federation. This research was done using resources provided by the Open Science Grid, which is supported by the National Science Foundation and the U.S. Department of Energy’s Office of Science. The authors would like to acknowledge the use of the Oxford Particle Physics Computing Cluster. This work was supported by the Comisión Nacional de Investigación Científica y Tecnológica (CONICYT), Chile; the Ministry of Education, Youth and Sports, Czech Republic, under Grant INGO II-LG 14033; the DFG cluster of excellence “Origin and Structure of the Universe”, Germany; the EC HIGGSTOOLS Project, under Contract PITN-GA-2012-316704; the European Union’s Horizon 2020 Research and Innovation programme under Grant Agreement No. 654168; the German - Israel Foundation (GIF); the Israel Science Foundation (ISF); the I-CORE programme of VATAT, ISF and the Israel Academy of Sciences, Israel; the Research Council of Norway; the Ministry of Education, Science and Technological Development of the Republic of Serbia through the National Project OI171012; the Polish Ministry of Science and Higher Education under Contract No. 3501/H2020/2016/2; the National Science Centre, Poland, HARMONIA Project, under Contracts 2013/10/M/ST2/00629 and UMO-2015/18/M/ST2/00518; the Romanian agencies UEFISCDI and ROSA; the Secretary of State of Research, Development and Innovation of Spain, under Project FPA2011-15330-E, FPA2015-71956-REDT; the Gates Foundation, United Kingdom; the UK Science and Technology Facilities Council (STFC), United Kingdom; and the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences and Office of High Energy Physics under Contract DE-AC02-06CH11357.

Footnotes

1As a result of this paper and other studies, a slightly different staging scenario for CLIC, with a first stage at s=380GeV to include precise measurements of top quark properties as a probe for BSM physics, and the next stage at 1.5 TeV, has recently been adopted and will be used for future studies [19].

2A Higgs boson of 125 GeV is used in the process e+e-tt¯H.

References

1. Aad G, et al. Phys. Lett. B. 2012;716:1. doi: 10.1016/j.physletb.2012.08.020. [Cross Ref]
2. Chatrchyan S, et al. Phys. Lett. B. 2012;716:30. doi: 10.1016/j.physletb.2012.08.021. [Cross Ref]
3. Englert F, Brout R. Phys. Rev. Lett. 1964;13:321. doi: 10.1103/PhysRevLett.13.321. [Cross Ref]
4. Higgs P. Phys. Lett. 1964;12:132. doi: 10.1016/0031-9163(64)91136-9. [Cross Ref]
5. Higgs P. Phys. Rev. Lett. 1964;13:508. doi: 10.1103/PhysRevLett.13.508. [Cross Ref]
6. Guralnik G, Hagen C, Kibble T. Phys. Rev. Lett. 1964;13:585. doi: 10.1103/PhysRevLett.13.585. [Cross Ref]
7. Higgs P. Phys. Rev. 1966;145:1156. doi: 10.1103/PhysRev.145.1156. [Cross Ref]
8. Kibble T. Phys. Rev. 1967;155:1554. doi: 10.1103/PhysRev.155.1554. [Cross Ref]
9. H. Abramowicz et al., Physics at the CLICe+e-e+e-Linear Collider – Input to the Snowmass process 2013 (2013). arxiv:1307.5288
10. M. Aicheler et al. (eds.), A Multi-TeV Linear Collider based on CLIC Technology: CLIC Conceptual Design Report (CERN) (2012). doi:10.5170/CERN-2012-007
11. L. Linssen et al. (eds.), Physics and Detectors at CLIC: CLIC Conceptual Design Report (CERN) (2012). doi:10.5170/CERN-2012-003
12. H. Aihara et al. (eds.), SiD Letter of Intent (2009). SLAC-R-989. arxiv:0911.0006
13. T. Behnke et al. (eds.), The International Linear Collider Technical Design Report - Volume 4: Detectors (2013). ILC-Report-2013-040. arxiv:1306.6329
14. T. Abe et al., The International Large Detector: Letter of Intent (2010). DESY 2009-87. arxiv:1006.3396
15. Thomson M. Nucl. Instrum. Methods A. 2009;611:25. doi: 10.1016/j.nima.2009.09.009. [Cross Ref]
16. Marshall J, Münnich A, Thomson M. Nucl. Instrum. Methods A. 2013;700:153. doi: 10.1016/j.nima.2012.10.038. [Cross Ref]
17. Buskulic D, et al. Nucl. Instrum. Methods A. 1995;360(3):481. doi: 10.1016/0168-9002(95)00138-7. [Cross Ref]
18. The CMS Collaboration, Particle-flow event reconstruction in CMS and performance for jets, taus, and missing ET (2009). https://cds.cern.ch/record/1194487/. CMS PAS PFT-09-001
19. P. Burrows et al. (eds.), Updated baseline for a staged Compact Linear Collider (CERN) (2016). doi:10.5170/CERN-2016-004
20. Kilian W, Ohl T, Reuter J. Eur. Phys. J. C. 2011;71:1742. doi: 10.1140/epjc/s10052-011-1742-y. [Cross Ref]
21. Olive KA, et al. Chin. Phys. C. 2014;38:090001. doi: 10.1088/1674-1137/38/9/090001. [Cross Ref]
22. S. Dittmaier et al., Handbook of LHC Higgs cross sections: 2. Differential distributions (2012). doi:10.5170/CERN-2012-002
23. Beringer J, et al. Phys. Rev. D. 2012;86:010001. doi: 10.1103/PhysRevD.86.010001. [Cross Ref]
24. S. Dawson et al., Working Group Report: Higgs Boson. FERMILAB-CONF-13-671-T. arxiv:1310.8361
25. H. Baer et al. (eds.), The International Linear Collider Technical Design Report - Volume 2: Physics (2013). ILC-Report-2013-040. arxiv:1306.6352
26. Agostinelli S, et al. Nucl. Instrum. Methods A. 2003;506:250. doi: 10.1016/S0168-9002(03)01368-8. [Cross Ref]
27. Allison J, et al. IEEE Trans. Nucl. Sci. 2006;53:270. doi: 10.1109/TNS.2006.869826. [Cross Ref]
28. Sjöstrand T, Mrenna S, Skands PZ. JHEP. 2006;0605:026. doi: 10.1088/1126-6708/2006/05/026. [Cross Ref]
29. Alexander G, et al. Z. Phys. C. 1996;69:543. doi: 10.1007/s002880050059. [Cross Ref]
30. Was Z. Nucl. Phys. Proc. Suppl. 2001;98:96. doi: 10.1016/S0920-5632(01)01200-2. [Cross Ref]
31. K. Fujii, Physics Study Libraries. http://www-jlc.kek.jp/subg/offl/physsim/. Accessed 25 Aug 2016
32. P. Mora de Freitas, H. Videau, Detector simulation with MOKKA / GEANT4: present and future (2002). LC-TOOL-2003-010
33. Graf N, McCormick J. AIP Conf. Proc. 2006;867:503. doi: 10.1063/1.2396991. [Cross Ref]
34. Gaede F. Nucl. Instrum. Methods A. 2006;559:177. doi: 10.1016/j.nima.2005.11.138. [Cross Ref]
35. Graf NA. J. Phys. Conf. Ser. 2011;331:032012. doi: 10.1088/1742-6596/331/3/032012. [Cross Ref]
36. Marshall JS, Thomson MA. Eur. Phys. J. C. 2015;75:439. doi: 10.1140/epjc/s10052-015-3659-3. [Cross Ref]
37. Suehara T, Tanabe T. Nucl. Instrum. Methods A. 2016;808:109. doi: 10.1016/j.nima.2015.11.054. [Cross Ref]
38. P. Schade, A. Lucaci-Timoce, Description of the signal and background event mixing as implemented in the Marlin processor OverlayTiming (2011). http://cds.cern.ch/record/1443537. LCD-Note-2011-006
39. Cacciari M, Salam GP, Soyez G. Eur. Phys. J. C. 2012;72:1896. doi: 10.1140/epjc/s10052-012-1896-2. [Cross Ref]
40. Catani S, et al. Phys. Lett. B. 1991;269:432. doi: 10.1016/0370-2693(91)90196-W. [Cross Ref]
41. Catani S, et al. Nucl. Phys. B. 1993;406:187. doi: 10.1016/0550-3213(93)90166-M. [Cross Ref]
42. Ellis SD, Soper DE. Phys. Rev. D. 1993;48:3160. doi: 10.1103/PhysRevD.48.3160. [PubMed] [Cross Ref]
43. Grefe C, et al. J. Phys. Conf. Ser. 2014;513:032077. doi: 10.1088/1742-6596/513/3/032077. [Cross Ref]
44. Tsaregorodtsev A, et al. J. Phys. Conf. Ser. 2008;119:062048. doi: 10.1088/1742-6596/119/6/062048. [Cross Ref]
45. Thomson M. Eur. Phys. J. C. 2016;76:72. doi: 10.1140/epjc/s10052-016-3911-5. [Cross Ref]
46. J. Marshall, Higgs Mass and Cross-Section Measurements at a 500 GeVCLIC Machine, Operating ats=350GeVand 500 GeV (2012). https://cds.cern.ch/record/1498592. LCD-Note-2012-015
47. Therhaag J. AIP Conf. Proc. 2009;1504:1013.
48. Cranmer KS. Comput. Phys. Commun. 2001;136:198. doi: 10.1016/S0010-4655(00)00243-5. [Cross Ref]
49. K.S. Cranmer, Kernel Estimation for Parameterizing Distributions of Discriminating Variables (1999). ALEPH Note 99-144
50. The OPAL Collaboration, Search for Neutral Higgs Bosons in e+e-e+e-Collisions ats=192-202GeV (2000). OPAL Physics Note PN426
51. A.L. Read, Modified frequentist analysis of search results (the CL_s method) (2000). doi:10.5170/CERN-2000-005.81. http://cds.cern.ch/record/451614/. CERN-OPEN-2000-205
52. A. Münnich, Measurement ofσ(e+e-Hνν¯)×BR(Hτtau)at CLIC @350GeV350GeV (2013). http://cds.cern.ch/record/1504764. LCD-Note-2012-019
53. A. Münnich, Taufinder: A Reconstruction Algorithm for tau Leptons at Linear Colliders (2010). http://cds.cern.ch/record/1443551/. LCD-Note-2010-009
54. A. Münnich, Measurement ofσ(e+e-Hνν¯)×BR(Hττ)at CLIC @1.4TeV1.4TeV (2013). http://cds.cern.ch/record/1558360/. LCD-Note-2012-010
55. Kawada S, et al. Eur. Phys. J. C. 2015;75:617. doi: 10.1140/epjc/s10052-015-3854-2. [Cross Ref]
56. A. Winter, N. Watson, Measurement of theH → WWHWWBranching Ratio at1.4TeV1.4TeVusing the semileptonic final state at CLIC (2016). http://cds.cern.ch/record/2157035/. CLICdp-Note-2016-003
57. C. Grefe, E. Sicking, Physics potential of theσe+e-Hνeν¯e×BRHZγmeasurement at a 1.4 TeV Compact Linear Collider (2014). http://cds.cern.ch/record/1749741. CLICdp-Note-2014-003
58. Grefe C, Lastovicka T, Strube J. Eur. Phys. J. C. 2013;73(2):2290. doi: 10.1140/epjc/s10052-013-2290-4. [Cross Ref]
59. Milutinović-Dumbelović G, et al. Eur. Phys. J. C. 2015;75:515. doi: 10.1140/epjc/s10052-015-3742-9. [Cross Ref]
60. S. Redford, P. Roloff, M. Vogel, Physics potential of the top Yukawa coupling measurement at a 1.4 TeV Compact Linear Collider using theCLICSiDCLICSiDdetector (2014). http://cds.cern.ch/record/1690648. CLICdp-Note-2014-001
61. S. Redford, P. Roloff, M. Vogel, Study of the effect of additional background channels on the top Yukawa coupling measurement at a 1.4 TeV CLIC (2015). http://cds.cern.ch/record/1982243/. CLICdp-Note-2015-001
62. Djouadi A, Kalinowski J, Zerwas PM. Z. Phys. C. 1992;54:255. doi: 10.1007/BF01566654. [Cross Ref]
63. Djouadi A, Kalinowski J, Zerwas PM. Mod. Phys. Lett. A. 1992;7:1765. doi: 10.1142/S0217732392001464. [Cross Ref]
64. Price T, et al. Eur. Phys. J. C. 2015;75:309. doi: 10.1140/epjc/s10052-015-3532-4. [Cross Ref]
65. Gupta RS, Rzehak H, Wells JD. Phys. Rev. D. 2013;88:055024. doi: 10.1103/PhysRevD.88.055024. [Cross Ref]
66. Poss S, Sailer A. Eur. Phys. J. C. 2014;74:2833. doi: 10.1140/epjc/s10052-014-2833-3. [Cross Ref]
67. Lukić S, et al. JINST. 2013;8:P05008. doi: 10.1088/1748-0221/8/05/P05008. [Cross Ref]
68. Božović-Jelisavčić I, et al. JINST. 2013;8:P08012. doi: 10.1088/1748-0221/8/08/P08012. [Cross Ref]
69. G. Wilson, Beam Polarization Measurement Using Single Bosons with Missing Energy (2012). International Workshop on Future Linear Colliders (LCWS12)
70. James F, Roos M. Comput. Phys. Commun. 1975;10(6):343. doi: 10.1016/0010-4655(75)90039-9. [Cross Ref]

Articles from Springer Open Choice are provided here courtesy of Springer