Physics potential for the measurement of σ ( H νν̄ ) × BR ( H → μ + μ − ) at the 1 . 4 TeV CLIC collider ?

The future Compact Linear Collider (CLIC) offers a possibility for a rich precision physics programme, in particular in the Higgs sector through the energy staging. This is the first paper addressing the measurement of the Standard Model Higgs boson decay into two muons at 1.4 TeV CLIC. With respect to similar studies at future linear colliders, this paper includes several novel contributions to the statistical uncertainty of the measurement. The later includes the Equivalent Photon Approximation and realistic forward electron tagging based on energy deposition maps in the forward calorimeters, as well as several processes with the Beamstrahlung photons that results in irreducible contribution to the signal. In addition, coincidence of the Bhabha scattering with the signal and background processes is considered, altering the signal selection efficiency. The study is performed using a fully simulated CLIC_ILD detector model. It is shown that the branching ratio for the Higgs decay into a pair of muons BR(H→ μ+μ−) times the Higgs production crosssection in WW -fusion σ(Hνν̄) can be measured with 38% statistical accuracy at √ s = 1.4 TeV, assuming an integrated luminosity of 1.5 ab−1 with unpolarised beams. If 80% electron beam polarisation is considered, the statistical uncertainty of the measurement is reduced to 25%. Systematic uncertainties are negligible in comparison to the statistical uncertainty.


Introduction
Measurements of Higgs branching ratios, and consequently Higgs couplings, provide a strong test of the Standard Model (SM) and possible physics beyond. Models that could possibly extend the SM Higgs sector (Two Higgs Doublet model, This work was carried out in the framework of the CLICdp collaboration a e-mail: gordanamd@vinca.rs Little Higgs models or Compositeness models) will require Higgs couplings to electroweak bosons and Higgs-fermion Yukawa couplings (coupling-mass linearity) to deviate from the SM predictions [1,2]. The Compact Linear Collider (CLIC) represents an excellent environment to study properties of the Higgs boson, including its couplings, with a very high precision [3,4]. Measurements of rare H → µ + µ − decays are particularly challenging because of the very low branching ratio of 2×10 −4 predicted in the SM [5] for a Higgs mass of 126 GeV. Current results indicate that the LHC was not able to access Higgs coupling to muons (g Hµ µ ), based on the runs at 7 TeV and 8 TeV centre-of-mass (CM) energies [6]. Projections for the HL-LHC, assuming 300 fb −1 and 3 ab −1 of data, predict uncertainties of 23% and 8% respectively for the g Hµ µ coupling [7]. In order to provide the best physics reach in the shortest time and for an optimal cost, the operation of the CLIC accelerator is foreseen in energy stages of 350 GeV, 1.4 TeV and 3 TeV [8]. At 1.4 TeV and 3 TeV, sufficiently large Higgs boson samples can be produced to allow studies of rare Higgs decays . A sample of 3.7 ×10 5 Higgs bosons can be produced at 1.4 TeV CM energy, for an integrated luminosity of 1.5 ab −1 with unpolarised beams. With the expected instantaneous luminosity of 3.2×10 34 cm −2 s −1 this can be achieved in approximately five years of detector operation, with 200 running days per year and an effective up-time of 50% [9]. The signal sample size will be doubled at 3 TeV CM energy due to rising cross-section for WWfusion [10].
A similar study has been performed at 3 TeV CM energy [10]. Compared to the study at 3 TeV, several challenges for the measurement of H → µ + µ − at CLIC are discussed for the first time in this paper. First, background processes with photons in the initial state simulated using both the expected Beamstrahlung spectrum at CLIC and the Equivalent Photon Approximation (EPA) [11,12], were considered. Forward 2 electron tagging (Section 5) leads to a rejection of 48% and 42% of the e − e + → e − e + µ + µ − and e ± γ → e ± µ + µ − background events, respectively. The impact of Bhabha scattering events on the rejection of events with forward electrons is investigated.
The paper is organised as follows. In Section 2 the simulation tools used for the analysis are listed and in Section 3 the CLIC_ILD detector model is briefly described. Signal and background processes and event samples are discussed in Section 4. Tagging of background high-energy electrons is described in Section 5. Event preselection and the final selection based on a Multivariate Analysis (MVA) approach are described in Section 6. The di-muon invariant mass fit and the extraction of the statistical uncertainty of the measurement are described in Section 7. In Section 8 the impact of electron polarisation on the statistical uncertainty of the σ (Hνν) × BR(H → µ + µ − ) measurement is described. Systematic uncertainties are discussed in Section 9, followed by the conclusions in Section 10.

Simulation and analysis tools
Higgs production through WW -fusion is simulated in Whizard 1.95 [13,14] including a realistic CLIC beam spectrum and initial state radiation. The generator Pythia 6.4 [15] is used to simulate the Higgs decay into two muons. The CLIC luminosity spectrum and beam-induced processes are obtained by GuineaPig 1.4.4 [16]. Background events are also generated with Whizard using Pythia 6.4 to simulate the hadronization and fragmentation processes. Simulation of tau decays is done by Tauola [17]. The CLIC_ILD detector simulation is performed using Mokka [18] based on Geant4 [19]. Before digitisation of the detector signals, pile-up from γγ → hadrons interactions is overlaid on the physics events. The particle flow algorithm, PandoraPFA [20,21] is employed in reconstruction of the final-state particles within the Marlin reconstruction framework [22]. The TMVA package [23] is used to separate signal from background by MVA of signal and background kinematic properties.

The CLIC_ILD detector model
The ILD detector concept [24] is modified for CLIC according to the specific experimental conditions at higher energies [3]. The subsystems of particular relevance for the presented analysis are discussed here. A complete description of the CLIC_ILD detector can be found in [25] .
The main tracking device of CLIC_ILD is the Time Projection Chamber (TPC) providing a point resolution in the rφ plane better than 100 µm, for charged particles in the detector angular acceptance [3]. Additional silicon trackers cover polar angles down to 7°. They have a single point resolution of 7 µm, and together with the TPC improve the tracking accuracy in the rφ plane. In order to provide precision tracking and vertexing closer to the beam-pipe, a Vertex Detector capable of an impact parameter resolution of 3 µm [26] is foreseen. Calorimetry at CLIC is based on fine-grained sandwich calorimeters optimized for particleflow analysis (PFA). PFA is based on reconstruction of fourvectors of visible particles, combining the information from precise tracking with highly granular calorimetry. The detector comprises a central solenoid magnet, with a field of 4 T. High muon reconstruction efficiency of 99%, for muons above 7.5 GeV, is achieved by combining information from the central tracker (TPC plus silicon tracker) with information provided by the iron yoke instrumented with the 9 layers of Resistive Plate Chamber detectors.
In principle, hadrons produced in the interaction of the beam-induced photons affect the TPC occupancy and consequently the muon reconstruction efficiency. However, in the studied sample of muons from H → µ + µ − decays, muon reconstruction efficiency is above 99% in the barrel region, in the presence of γγ → hadrons.
The average muon transverse momentum resolution for the signal sample is ∆ (1/p T ) = 3.3×10 −5 GeV −1 in the barrel region. The impact of transverse momentum resolution on the statistical uncertainty of σ (Hνν) × BR(H → µ + µ − ) measurement is discussed in Section 9.
In the very forward region of the CLIC_ILD detector, below θ =8°, no tracking information is available. The region between 0.6°and 6.3°is instrumented with the two silicon-tungsten sampling calorimeters, LumiCal and Beam-Cal [27], for the luminosity measurement, beam-parameter control, as well as for the tagging of high-energy electrons escaping the main detector at low angles. Together with the very forward segments of the electromagnetic calorimeter covering the polar-angle region between 6.3°and 8°, it is possible to suppress the four-fermion SM background with the characteristic low-angle electron signature. The simulation of the very-forward electron tagging is described in Section 5.

Event samples
At √ s = 1.4 TeV the SM Higgs boson is predominantly produced via WW -fusion ( Figure 1). The effective cross-section for Higgs production in WW -fusion is 244 fb without beam polarization. The Higgs production cross-section above 1 TeV can be measured with a statistical precision better than 1% as shown in [1]. The e + e − → Hνν, H → µ + µ − signal statistics are expected to be small (of the order of a few tens of events) because of the small branching fraction for this particular decay.
We have simulated a sample of 24000 signal events, roughly corresponding to 300 times the number of events expected in 1.5 ab −1 of data. This is needed in order to provide an adequate description of the signal Probability Density Function (PDF) (Section 7.1). The signal and the dominant background processes are listed in Table 1. For each of the background processes, samples of 2 ab −1 are generated.
In addition to the processes listed in Table 1, we have considered s-channel e + e − → µ + µ − production, as well as several processes with tau pair in the final state e + e − → τ + τ − , e + e − → ν τντ τ + τ − , e + e − → e − e + τ + τ − . Tau decays become relevant if both taus decay into two muons which happens in ∼3% of cases [28]. However, the invariant mass of the di-muon system will not match the Higgs mass window considered in this analysis (see Section 6.1). The same holds for e + e − → µ + µ − production.
In addition, processes with quasi-real photons in the initial state are simulated using the EPA. In this analysis, such events are grouped together with the processes involving Beamstrahlung photons. In this way, processes with roughly similar kinematic characteristics are grouped together , as shown in Table 1. The notation e ± γ represents the sum of cross-sections for the processes with either e − γ or e + γ in the inital state.

At
√ s = 1.4 TeV, the Higgs boson is also produced via ZZ-fusion, with a cross-section of about 10% of the Higgs production cross-section in WW -fusion. However, on a test sample of 300 ZZ-fusion events followed by the Higgs decay to a pair of muons, not a single event passed the selection described in Sections 6.1 and 6.2. This implies an efficiency smaller than 1.2% (95% CL) for this channel equivalent to less than 0.1 events passing the final selection. Therefore, the Higgs production through ZZ-fusion is not considered relevant for this analysis.
Photons, dominantly emitted by Beamstrahlung, produce incoherent pairs deposited mainly in the low-angle calorimeters. On average, 1.3 two-photon interactions producing hadronic final states occur per bunch crossing [29] which may affect the muon reconstruction in the tracking detectors. These hadrons are included in the analysis by overlaying 60 bunch crossings in the simulation, before the digitisation and event reconstruction phase. These events, as well as other physics events, are passed trough the full detector simulation [30]. Table 1 List of considered processes with corresponding crosssections. Cross-section values marked by * are generated with the additional kinematic requirements: 100 GeV< m µ µ <150 GeV, and 8°<θ µ <172°, where m µ µ stands for di-muon invariant mass and θ µ is the polar angle of the reconstructed muon. The cross-sections for all processes with photons in the initial state include both Beamstrahlung and processes with EPA photons. Cross-sections for processes e ± γ → e ± µ + µ − and e ± γ → e ± ν µνµ µ + µ − represent the sum of cross-sections for the processes with both initial states e − γ and e + γ. Process

Tagging of EM showers in the very forward region
In the polar angle region below θ =8°, tracking information and hadronic calorimetry are not available. The four-fermion background e + e − → e + e − µ − µ + of multiperipheral type and similar processes like e ± γ → e ± µ − µ + can fake the missing energy signature of the signal if the final state electrons (spectators), emitted at the polar angles smaller than θ =8°, escape undetected. Electron detection in the very forward region involves the reconstruction of electromagnetic showers in the presence of intense beam-induced background depositing in the very forward calorimeters a large number of low-energy particles, mostly incoherent pairs from Beamstrahlung [31]. This deposition amounts to several hundred thousand of e + e − pairs per bunch crossing [32].
Furthermore, Bhabha events where one or both electrons are detected in the very forward calorimeters may occur in coincidence with either signal or background, even within the 10 ns time stamp. Tagging of such Bhabha electrons will result in the rejection of signal (or background). In order to prevent significant loss of signal statistics, the electron tagging was optimized to identify showers with energy higher than 200 GeV and a polar angle above 1.7°only 1 . Under these requirements, the loss of the number of signal events due to tagging of Bhabha electrons amounts to 7%. Out of these 7%, in slightly more than a half of events one electron is added to the final state and, in the remainder two 4 Bhabha electrons are added. Table 2 shows rejection rates for signal and background obtained by the forward electron tagging due to Bhabha pile-up.
In conclusion, very forward tagging of high-energy electrons serves to half the fraction of background with spectator electrons, with a moderate loss of signal of 7% in the presence of Bhabha coincidence. Table 2 Rejection rates for signal and background by the forward electron tagging.

Process
Rejection rate The event selection is done in two steps. First, a preselection is performed aiming to suppress background originating from beamstrahlung as well as the processes with spectator electrons described in Section 5. The final event selection uses a multivariate classifier based on boosted decision trees (BDT) to suppress remain background processes on the basis of their kinematic properties.

Preselection
In order to suppress the impact of the beam-induced background, only reconstructed particles with transverse momenta p T > 5 GeV are used in the analysis. Furthermore, the preselection of events was made by requiring a reconstruction of exactly two muons in the event, with an invariant mass of the di-muon system in the window centered around the Higgs mass 105-145 GeV. In addition, the absence of tagged electrons with energy above 200 GeV and polar angle above 1.7°is required in order to suppress background with spectator electons emitted in the very forward region of the detector.

MVA selection
As a second step in the event selection, MVA techniques are used based on the BDT classifier implemented in the TMVA package. From the signal sample, quater of all events are reserved for TMVA training, as well as 0.5 ab −1 of each background. The following observables were used for the classification of events , similar to the CLIC study at √ s = 3 TeV [10]: visible energy of the event excluding the energy of the di-muon system, E vis , transverse momentum of the di-muon system, p T (µ µ), scalar sum of the transverse momenta of the two selected muons, p T (µ 1 ) + p T (µ 2 ), boost of the di-muon system, β µ µ = p µ µ /E µ µ , polar angle of the di-muon system, θ µ µ , cosine of the helicity angle, cos θ * .
The process e + e − → νν µ + µ − , with the same final state as the signal, represents an irreducible background and can not be substantially suppressed before the invariant mass fit of the di-muon system. The process γγ → ν µνµ µ + µ − has a similar final state, but a different CM energy distribution in the initial state, since it involves Beamstrahlung or EPA photons rather than initial electrons. This leads to a different distribution of the boost of the di-muon system, allowing separation from the signal to some extent. The processes e + e − → νν µ + µ − and γγ → ν µνµ µ + µ − , have slightly different distributions of the helicity angle from the signal. All processes with one or two spectator electrons show significant differences from the signal, primarily in the distribution of the visible energy ( Figure 2). These processes are also effectively suppressed by the p T (µ 1 ) + p T (µ 2 ) observable. In addition, for the e + e − → e + e − µ + µ − process, the distribution of p T (µ µ) exhibits a peak at lower values than the signal (Figure 3). This peak corresponds to events in which the di-muon system recoils against electron spectators or outgoing photons that are emitted below the angular cut of the very forward EM-shower tagging. The above is illustrated in Figure 3 showing the p T distributions for representative background processes. Fig. 2 Distribution of the visible energy for the signal and e ± γ → e ± µ + µ − , e + e − → νν µ + µ − , e ± γ → e ± ν µνµ µ + µ − and e + e − → e + e − µ + µ − background. The distribution of the BDT classifier variable for the signal and the main background processes is shown in Figure 4 (a). The classifier cut position was selected to maximise the significance, defined as N s / √ N s + N b , where N s and N b are the number of selected signal and background events, respectively. A plot of significance as a function of the position of the BDT cut is shown in Figure 4 (b). The optimal cut position was found at BDT = 0.23. Distributions of the di-muon invariant mass before and after the MVA selection are shown in Figure 5. Figure 5 (a) includes all events that pass the preselection, while Figure 5 (b) shows all events passing the MVA selection. All samples are normalised to the integrated luminosity of 1.5 ab −1 . The signal preselection efficiency is 82%. The MVA selection efficiency for the signal is 32%, reflecting the fact that sensitive observables have limited power to discriminate between the signal and background. The overall signal efficiency including reconstruction, preselection, losses due to coincident tagging of Bhabha particles and the MVA is 24%, resulting in an expected number of 19 signal events.

Di-muon invariant mass fit
The quantity σ (Hνν) × BR(H → µ + µ − ) is determined from the equation: where L stands for the integrated luminosity and ε s is the total counting efficiency for the signal, including the reconstruction, preselection and MVA selection. In the experiment, the number of signal events N s will be determined by fitting the di-muon invariant mass distribution with a function f (m µ µ ): where f s,b are probability density functions (PDF) used to describe the signal and the sum of all background processes, and N s and N b are the respective numbers of signal and background events in the fitting mass window. In this analysis, an unbinned likelihood fit, with all parameters of f s,b (m µ µ ) fixed, is performed on simulated signal and background samples. N s and N b are left as free parameters determined from the fit. The way the signal and background PDFs are obtained is discussed in Section 7.1.
In order to estimate the statistical uncertainty of the signal count, 5000 toy Monte Carlo (MC) experiments are performed, where pseudo-data are obtained by randomly picking the signal m µ µ values from the fully simulated signal sample, while background m µ µ values are randomly generated from the total background PDF f b (m µ µ ). The size of the signal sample N s and sample sizes N b,i of individual backgrounds considered, are obtained from the Poisson distribution for the integrated luminosity of 1.5 ab −1 , taking into consideration corresponding cross-sections σ and the selection efficiencies ε ( N s = L · σ s · ε s , N b,i = L · σ i · ε i , where i is indexing the different background processes listed in Table 1).
For each toy MC experiment, the m µ µ distribution is fitted by the function f (m µ µ ) given in Eq.2, and the standard deviation of the resulting distribution of N s over all toy MC experiments is taken as the estimate of the statistical uncertainty of the σ (Hνν) × BR(H → µ + µ − ) measurement.
As will be discussed in Section 9.1, the di-muon invariant mass distribution is sensitive to the detector p T resolution, while the Higgs width Γ H can be considered negligible in comparison to the detector energy resolution.

Signal and background PDFs
Fully simulated samples of signal and background (Table 1) are fitted to extract the PDFs. The sizes of the samples vary from several tens of thousands of events for the signal, up to a few million of events for various background processes.
The signal PDF was defined as a linear combination of a Gaussian function with exponential tails, f exp and a Gaussian function with tails that asymptotically approach con- The parameters of Eq. 3 are determined by fitting the dimuon invariant mass distribution for the signal (Figure 6).
The total background PDF is defined as a linear combination of a constant and exponential term: The di-muon invariant mass fit of the total background is shown in Figure 7, together with the fit results for the free parameters in Eq.4. As the normalisation to the common integrated luminosity requires different normalisation coefficients for different processes, binned data were used 7 to combine the background processes in a straightforward manner and a binned χ 2 fit was performed. The χ 2 /N d f of the background fit was 62/61. The overall function f (m µ µ ) (Eq.2) is fitted to the pseudodata of each toy MC experiment using the unbinned likelihood fit. An example of a toy MC fit is given in Figure 8.
The standard deviation of the resulting signal count distribution in 5000 repeated toy MC experiments corresponds to the statistical uncertainty of the measurement and is 38%. (Figure 9). According to Eq.1 it translates into the statistical uncertainty of the σ (Hνν) × BR(H → µ + µ − ) measurement, having in mind that the total uncertainty of the integrated luminosity can be determined at the permille level [34].
The statistical uncertainty of the signal counting is dominated by contributions from the limited signal statistics and from a presence of irreducible backgrounds. To estimate the significance of the signal against the null-hypothesis, another set of 5000 toy MC experiments was performed with zero signal count, and f (m µ µ ) (Eq.2) was fitted to the pseudodata. The resulting N s distribution was centered on zero with a standard deviation of 5.4. Thus, in the case where the SM expected number of 19 signal events are found in an experiment, the corresponding signal significance would be 3.7 σ .
The Higgs coupling to muons, g Hµ µ , is optimally extracted in a global fit procedure taking into account all Higgs measurements at the 350 GeV, 1.4 TeV and 3 TeV stages. The global fit serves to extract Higgs couplings from all measurements, as well as the experimental Higgs width Γ H . and having access to Γ H and g HWW from other measurements, extraction of g Hµ µ is possible solely from the measurement presented here. An example of a minimal set of measurements giving a model-independent access to Γ H and g HWW is the following: the H → bb measurements at both 350 GeV and 1.4 TeV give access to the ratio g HWW g HZZ , the recoil mass measurement at 350 GeV CM energy gives access to g HZZ , and the H → W + W − measurement at 1.4 TeV gives access to the ratio g 4 HWW Γ H [4]. The contributions of these measurements towards the final ∆ g Hµ µ is negligible at the second significant digit. The dominant contribution to the g Hµ µ coupling uncertainty is the statistical uncertainty of the measurement presented here. Systematic uncertainties affect the total uncertainty of g Hµ µ determination only at the third significant digit, and thus can be neglected (Section 9). Under these assumptions, the relative uncertainty of g Hµ µ is approximated to be 19%.

Impact of electron polarization
If 80% left-handed polarisation of the electron beam is assumed during the entire operation time at 1.4 TeV, the Higgs production cross-section through WW -fusion would be enhanced by a factor 1.8 [4]. The most important background contribution after the MVA selection, the e + e − → ν eνe µ + µ − process, is enhanced by the same factor because it is also mediated by W bosons which have only left-handed interactions. The process e ± γ → e ± µ + µ − is enhanced by a factor 1.32, while cross-sections for other background processes are not significantly changed w.r.t. the unpolarised case. The overall selection efficiency of the signal is 30%, because the classifier cut position is moved to a lower value which consequently leads to a higher signal efficiency. The final statistical uncertainty of the σ (Hνν) × BR(H → µ + µ − ) measurement is 25%. The corresponding uncertainty of g Hµ µ is 13%. A summary of the results of the presented analysis is given in Table 3. It is important to note that all kinematic variables are unaffected by the beam polarization.

Systematic uncertainties
From Eq.1 it is clear that uncertainties of the integrated luminosity and muon identification efficiency influence the uncertainty of the H → µ + µ − branching ratio measurement at the systematic level. It has been shown that at 3 TeV CLIC [35], where the impact of the beam-induced processes is the most severe, the luminosity above 75% of the nominal CM energy can be determined at the permille level, using lowangle Bhabha scattering. Below 75% of the nominal CM energy, the luminosity spectrum can be measured with a precision of a few percent using wide-angle Bhabha scattering [36]. About 17% of all Higgs production events occur at a CM energy below 75% of the nominal CM energy. Having in mind the intrinsic statistical limitations of the signal sample, this source of systematic uncertainty can be considered negligible.
On the detector side, an important systematic effect is the uncertainty on the transverse momentum resolution, because it directly influences the expected shape of the signal m µ µ distribution. The sensitivity of the signal count to the accuracy of the knowledge of the p T resolution σ p T has been studied by performing the analysis with an artificially introduced uncertainty of an exaggerated magnitude on the assumed p T resolution used to extract the signal PDF. Results of the relative shift in signal counts w.r.t. the relative shift of σ p T are shown in Figure 10. The relative bias in signal counting per one percent change of σ p T is 0.35%.
The uncertainty of the muon identification efficiency will directly influence the signal selection efficiency. In addition, the uncertainty of the muon polar angle resolution impacts the m µ µ reconstruction. Based on the results of the LEP experiments [37], it can be assumed that these detector related uncertainties are below a percent.
The systematic uncertainty of the signal count caused by the fit with f m µ µ defined in Eq.2, was found to be about 1% which is small compared to the statistical error.
Because of the forward EM shower tagging, about 7% of all events are rejected by coincident detection of Bhabha events. This fraction must be precisely calculated taking into account Bhabha event distributions, beam-beam effects, as well as the dependence of the tagging efficiency on energy and angle of the incident electrons and photons. This is work in progress [38][39][40], but the uncertainty of this effect is also expected to be negligible compared to the statistical uncertainty of the measurement. To estimate the benefit of a better p T resolution, the analysis was repeated by substituting the muon four-momenta reconstructed in the full simulation of the signal by the fourmomenta obtained by a parametrisation of the momentum resolution for several different values of the detector resolution. Figure 11 displays the approximate dependence of the statistical uncertainty of the measurement on the average transverse momentum resolution in the whole detector. Even a large improvement of the muon momentum resolution would result in only a moderate improvement of the statistical uncertainty of the measured product of the Higgs production cross-section and the branching ratio for the H → µ + µ − decay.

Conclusions
It has been shown that the measurement of the cross-section times the branching ratio for the SM Higgs decay into two muons can be performed with a relative statistical uncertainty of 38% at 1.4 TeV CLIC, assuming 1.5 ab −1 integrated luminosity with unpolarised beams. The result is dominated by the limited signal statistics and the irreducible background. The systematic uncertainties are negligible in comparison to the statistical one. This translates into a relative uncertainty of the coupling of Higgs to muons g Hµ µ of aproximately 19%. If the same integrated luminosity is collected with 80% left-handed polarisation for the electrons, the relative statistical uncertainty improves to 25% and 13% for σ (Hνν) × BR(H → µ + µ − ) and g Hµ µ , respectively.