Global analysis of the Higgs candidate with mass ∼ 125 GeV

: We analyze the properties of the Higgs candidate with mass ∼ 125 GeV discovered by the CMS and ATLAS Collaborations, constraining the possible deviations of its couplings from those of a Standard Model Higgs boson. The CMS, ATLAS and Tevatron data are compatible with Standard Model couplings to massive gauge bosons and fermions, and disfavour several types of composite Higgs models unless their couplings resemble those in the Standard Model. We show that the couplings of the Higgs candidate are consistent with a linear dependence on particle masses, scaled by the electroweak scale v ∼ 246 GeV, the power law and the mass scale both having uncertainties ∼ 20%.


JHEP09(2012)123
to photon and gluon pairs and assuming that the latter receive no significant contributions from particles beyond the Standard Model.
As reviewed below, one may parametrize the possible coupling deviations by coefficients a V and c f for vector bosons and fermions, respectively [51][52][53][54]. One possibility is that these coefficients are universal, i.e., a W = a Z ≡ a and c t = c b = c τ = c c = · · · ≡ c, with the Standard Model corresponding to a = c = 1. There has been some speculation that custodial symmetry might be broken with a W = a Z [55], and that couplings to some fermion species might be enhanced or suppressed. The present data are insufficient to probe these possibilities very precisely, and the overall quality of our global fit, presented below, indicates no need currently to adopt such hypotheses.
As already mentioned, a key prediction for the Standard Model Higgs boson is that its couplings to other particles are proportional to their masses. We probe this issue here by considering purely phenomenological parametrizations of the h couplings of the forms a V = v(M 2 V /M (1+2 ) ) and c f = v(m f /M 1+ ), where for a Standard Model Higgs boson = 0 and M = v = 246 GeV, the canonical Higgs vacuum expectation value (vev), corresponding to a = c = 1.
Figures 1, 2, 3 and 4 display our main results. They quantify the extent to which the h particle walks like a Higgs and quacks like a Higgs. Figure 1 shows the result in the (a, c) plane of our global fit to data on the h couplings from the Tevatron experiments and from the combined 7 and 8-TeV event samples of ATLAS and CMS. We see reasonable consistency with the Standard Model prediction: the overall best-fit region has c > 0 and, whilst the best fit has a > 1 and c < 1 (see also the marginalized one-dimensional likelihoods of our fit result projected on the a and c axes shown in figure 2), the Standard Model prediction lies within the 68% CL region. As we discuss in more detail below, figures 1 and 2 impose important constraints on composite Higgs models, disfavouring several such models unless their predictions resemble those of the Standard Model. Figure 3 displays the result of our global fit in the ( , M ) plane, where we see excellent consistency with the Higgs hypothesis: M = v, = 0. This is also seen in figure 4, which displays the marginalized one-dimensional likelihood projections of our fit result on the M and axes. The couplings of the h particle are clearly inconsistent with any massindependent scenario, which would require = −1. Figure 5 provides another way of understanding this observation. The points with vertical error bars represent the values of the h couplings to different particles found in our global fit to the parameters ( , M ). The diagonal dashed line is our best fit to ( , M ) and the dotted lines are given by the ±1σ ranges in these parameters, as given in the upper legend of the plot. The solid red line in figure 5 represents the Standard Model prediction (2.6), which is compatible within errors with the measurements, as already discussed.
In subsequent sections we describe how these results were obtained, and present more details of our analysis. In section 2 we review the phenomenological frameworks we employ, and in section 3 we describe our calculational procedure, which follows closely that in [56]. In section 4 we describe the data set we use, focusing in particular on the recent update from the Tevatron experiments as well as the recent preliminary results from ∼ 5/fb of JHEP09(2012)123 Figure 1. The constraints on the couplings (a, c) of the Higgs candidate h with mass ∼ 125 GeV obtained from our global analysis of the available CMS, ATLAS, CDF and D0 data. The Standard Model is represented by a black star, and the yellow lines represent various composite Higgs models described in the text, which are disfavoured if they deviate strongly from the Standard Model. In these and subsequent analogous plots, the most likely regions have the lightest shading, the dotted lines are 68% CL contours, the dashed lines are 95% CL contours, and the solid lines are 99% CL contours. 8-TeV LHC data in each of ATLAS and CMS. In section 5 we present in more detail our results in the (a, c) plane, discussing their implications for pseudo-dilatons [57][58][59][60][61][62][63][64][65][66][67][68][69][70][71][72] and other composite Higgs scenarios [73,74], 1 as well as fermiophobic (for a recent analysis see [78]) and gaugephobic [79] models. In section 6 we discuss in more detail our results in the ( , M ) plane, and in section 7 we present our conclusions and discuss the prospects that future data may soon clarify further the nature of the ∼ 125 GeV Higgs candidate h.

Phenomenological framework
We use the following nonlinear low-energy effective Lagrangian for the electroweak symmetry-breaking sector [51][52][53][54][57][58][59][60][61][62][63], see also [80][81][82][83][84][85]: where U is a unitary 2 × 2 matrix parametrizing the three Nambu-Goldstone fields that JHEP09(2012)123 are 'eaten' by the W ± and Z 0 , giving them masses, v ∼ 246 GeV is the conventional electroweak symmetry-breaking scale, h is a field describing the Higgs candidate with mass ∼ 125 GeV discovered by ATLAS and CMS, and λ f is the Yukawa coupling of the fermion flavour f in the Standard Model. The coefficients a and c f parametrize the deviations of the h couplings to massive vector bosons and to fermions, respectively, from those of the Higgs boson in the Standard Model. In writing (2.1), we have implicitly assumed a custodial symmetry: a W = a Z = a, an assumption whose plausibility can be judged from the overall quality of our fit. Also relevant for the phenomenology of the Higgs candidate h are its dimension-5 loop-induced couplings to gg and γγ [86][87][88][89]: We assume here that, as in the Standard Model, only the top quark makes a significant contribution to the coefficient b s , and only the top quark and the W ± contribute significantly to b em (with opposite signs in the Standard Model [86]). We recall that, in a scenario in which h is associated with a pseudo-dilaton field χ with vev V , one has One may also consider scenarios in which h is a pseudo-Goldstone boson of some higherorder chiral symmetry that is broken down to the SU(2) × SU(2) of the Standard Model
In addition to these theoretically-motivated models, we also consider the completely phenomenological possibility that h couples to other particles proportionally to some power of their masses. Thus, we generalize the Standard Model couplings to the following forms of couplings with anomalous scaling laws: The Standard Model is recovered in the double limit → 0, M → v, whereas the pseudodilaton/MCHM4 scenario would correspond to = 0 and M = V = v, in general. In terms of the parameterization (2.1), the parametrizations (2.7) correspond to After presenting our global fits to the parameters a, c, we shall explore the extent to which the data already indicate that h couples to other particles proportionally to masses, i.e., with = 0, and a normalization M similar to v = 246 GeV.

Calculational procedure
Assuming that the Higgs candidate h has no non-standard production or decay modes, its production cross-sections and decay widths are related to those of the Standard Model Higgs boson by simple factors of a and c. Assuming that gluon-gluon fusion (gg) and vectorboson fusion (VBF) dominate over the other processes, as at the LHC, one may combine

JHEP09(2012)123
Production sensitive to Decay sensitive to their respective production rescaling factors R ≡ σ/σ SM and cut efficiencies ξ gg,VBF to obtain a total production rescaling factor In the case of the Tevatron, where associated production (AP) is more important that VBF, one may use (3.1) with the replacement VBF → AP throughout. For the CMS diphoton subchannels, the collaboration provides a full breakdown of the percentage contribution from all production mechanisms which can be used directly instead of the ξ i F i factors above.
Similarly, relative to the Standard Model predictions, the decay widths R ≡ Γ/Γ SM to massive vector bosons, fermions and photons are given, respectively, by where the loop factors F t,w were given, e.g., in [56]. The principal dependences of the different Higgs-like signals on the rescaling factors (a, c) are summarized in table 1, which is adapted from [56]. It is important to emphasize that, since production mechanisms are in general also sensitive to both a and c, as well as decay branching ratios, their dependences also provide important constraints on model parameters. The signal strength modification factor µ i ≡ n i s /(n i s ) SM in any given channel i is the product of the production and decay rescalings: In the absence of more detailed experimental information, we follow [98] as in [56], assuming that in each JHEP09(2012)123 channel the underlying likelihood p(n obs |µn SM s + n b ) obeys a Poisson distribution, and use the approximation σ obs σ exp = µ 95% exp /1.96 for the standard deviation to solve for the central valueμ in the equation: The posterior probability density function is then given by with π(µ) generally assumed a priori to be flat within the range of interest and zero outside.

Experimental data set
We use the latest available information from ∼ 5/fb of LHC data obtained at each of 7 and 8 TeV in the centre of mass presented at ICHEP 2012 [15][16][17][18], and ∼ 10/fb of Tevatron data analyzed in [29,30]. In addition to [15,16], the CMS Collaboration provides additional information on its 7 and 8 TeV fits separately in public analysis notes, see below.
2. The searches in the h →bb VH channel are assumed to be dominated by associated production, with the Tevatron data updated from [56] to the latest results in [29,30], the 7 TeV Moriond results are used for ATLAS [5], and the CMS 7 and 8 TeV fits are obtained from [19]. In addition the CMS 7 TeVttH channel is included from [24].
3. The diphoton likelihoods in the ATLAS searches at 7 TeV were obtained from [7] as explained in [56]. In [20] CMS provides central values and one sigma error bars for both 7 and 8 TeV searches in four inclusive sub-channels dominated by gluon fusion and one or two di-jet categories. We treat the 8 TeV ATLAS results inclusively since the sub-channel best fit values are only provided for m h = 126.5 [26]. 2 The Tevatron search from [29,30] are also included.
4. The Tevatron results for h → W + W − are updated from [29,30]. The corresponding ATLAS results for 7 TeV [8] are supplemented by the 8 TeV data made public recently in [27]. CMS provide fits in the 0,1 and 2-jet categories for both 7 and 8 TeV [19,22].
5. The ATLAS τ τ searches at 7 TeV are treated as inclusive [10]. For CMS we use the best fits provided for the 0,1-jet and VBF channels at 7 and 8 TeV in [19,23]. At 7 TeV there is also an additional CMS search in the VH channel.

JHEP09(2012)123
As mentioned in the previous section, we use in our fit the CMS information on the percentage contribution from each production mechanism for all the diphoton sub-categories at 7 and 8 TeV separately. We treated the τ τ VBF categories assuming ∼30% contamination from gluon fusion in the production mechanism. As mentioned in [56], we expect that our global analysis is only accurate to ∼20% due to the limited experimental information available so far [99].

Tevatron data
We consider first the fit to the recent Tevatron data in terms of (a, c) that is shown in figure 6. We recall that the Tevatron experiments CDF and D0 provide information on the associated production (AP) of h followed by its decay intobb (upper left panel of figure 6), as well as inclusive measurements of h → W W * decay (upper right panel) and now also h → γγ decay (lower left panel). The central value of the h →bb signal is somewhat stronger than expected in the Standard Model, disfavouring fermiophobic models and corresponding to the possibility that either a and/or c > 1. However, the Tevatron h → W W * signal is relatively weak, disfavouring large a. The h → γγ signal is relatively strong, but very uncertain. In combination (lower right panel of figure 6), the Tevatron data are compatible with the Standard Model, while favouring slightly a, c > 1.

CMS data
We now turn to the analysis of the ICHEP 2012 CMS data shown in figure 7. The h →bb search (top left panel) was based on AP andtt Higgsstrahlung (HS) event selections, the former being more sensitive. The overall signal strength is somewhat below that expected in the Standard Model, slightly favouring a, c < 1, but very compatible with the Standard Model. The h → τ + τ − search (top right panel) was based on a combination of event selections favouring gluon-gluon fusion, VBF and AP production mechanisms. Once again the overall signal is weaker than expected in the Standard Model, but not very significantly. The h → ZZ * signal (middle left panel) has the strength expected in the Standard Model, disfavouring a 1. The h → W W * search shown in the middle right panel of figure 7 was based on a combination of event selections favouring gluon-gluon fusion and VBF production mechanisms, and the deficit compared to the Standard Model is not very significant. Finally, the h → γγ event selection includes samples with and without enhanced VBF contributions. As discussed in [56], since the h → γγ decay amplitude contains both t and W loops, which interfere, it provides unique discrimination between the cases a > and < 0, as seen in the bottom left panel of figure 7. Although the γγ signal strength is somewhat stronger than in the Standard Model, particularly in the VBF-enhanced sample, the discrepancy is not highly significant. Turning now to the overall combination of CMS data shown in the bottom right panel of figure 7, we see that the overall best fit is in a region with c < 1, driven by the γγ channel. The favoured region with c > 0 is compatible JHEP09(2012)123 with the Standard Model, with a ∼ 1 though c < 1 is somewhat favoured. 3 Comparing with the analogous panel in figure 1 of [56], we see that the accuracy in the determination of the h couplings has improved significantly.

ATLAS data
We now present a similar analysis of the available ATLAS data, which yields the results shown in figure 8. As in the previous figure, the top left panel displays the constraint in the (a, c) plane provided by the h →bb search, which in the ATLAS case is based on ∼ 5/fb of data at 7 TeV, as is the h → τ + τ − constraint shown in the top right panel of figure 8. These panels are the same as the corresponding panels in figure 5 of [56]. The middle left panel of figure 8 displays the h → ZZ * constraint including also ∼ 5/fb of data at 8 TeV: we see that the central value of the signal strength lies somewhat above the value expected JHEP09(2012)123

JHEP09(2012)123
in the Standard Model, corresponding to a > 1. The central value of the h → W W * signal shown in the middle right panel of figure 8 (which is based on ∼ 5/fb of data each at 7 TeV and 8 TeV) has a > 1, but is consistent with a = 1 at the 68% CL. Finally, the combined ATLAS 7-and 8-TeV search for h → γγ shown in the bottom right panel of figure 8 yields a central value of the strength lying somewhat above the Standard Model value, which is reflected in the preferred region of the (a, c) plane. The overall combination of the available ATLAS constraints, shown in the bottom right panel, indicates a general preference for a > 0, with values of a > 1, c < 1 being favoured.

Combined results in the (a, c) plane and implications for models
Looking at the bottom right panels of figures 6, 7 and 8, the general features of figure 1 can now be understood. Since they do not have high sensitivity to h → γγ, the Tevatron data are unable to discriminate between the upper and lower halves of the (a, c) plane. The CMS data do have high sensitivity to h → γγ, leading to some asymmetry between the upper and lower halves of the (a, c) plane, with a preference for c < 0. On the other hand, the ATLAS h → γγ and W W * data provide some preference for c > 0. Generally speaking, the Tevatron data prefer a, c > 1, whereas the CMS and ATLAS data prefer a > 1 and c < 1. The overall result, shown in figure 2, is that the global combination prefers a > 1 (left panel) and c > 0 (right panel), though not very significantly, and the favoured region with c > 0 has a slightly > 1 (left panel) and c slightly < 1 (right panel).
The absence of strongbb and τ + τ − signals at the LHC favours speculation that c 1, but in models with a universal coefficient c for all fermions, this is not the whole story. The fact the total h cross section is compatible with the Standard Model indicates that the htt couplings should be close to the Standard Model value, corresponding to c ∼ 1. Moreover, the Tevatron experiments report evidence for a strongbb decay signal, and as this is the dominant decay mode in the Standard Model the whole pattern of h decays would be very different if c 1. The right panel of figure 2 is the net result: no significant discrepancy with the Standard Model if c is universal.
Likewise, the absence of a strong h → W W * signal in the Tevatron, CMS and the ATLAS 7-TeV data might have led one to speculate that a < 1, or even that a W = a Z with custodial symmetry broken [55]. However, the fact that the h → γγ signals reported by both ATLAS and CMS are on the high side suggests that the hγγ loop amplitude in (3.2) must receive an important contribution from the W ± loops, which should be dominant. This and the ATLAS 8-TeV data suggest that a W cannot be very small, and favours a ≡ a W = a Z not ∼ 1 in our fit.

JHEP09(2012)123
are strongly favoured. These correspond to the cases V ∼ v in pseudo-dilaton models and ξ ∼ 0 in the MCHM4 model. However, in pseudo-dilaton models, in particular, there may be additional heavy particles contributing to the loop coefficients b s , b em in (2.2) [64][65][66][67][68][69][70][71][72], so this observation is model-dependent and further analysis is needed [100,101]. We also see in figure 1 that the ξ → 0 limit of the MCHM5 model is preferred, while 'antidilaton' models are slightly disfavoured compared to pseudo-dilaton models, and would prefer a = −c < 1: see also the right panel of figure 2. Finally, we observe that the fermiophobic and gaugephobic models are strongly disfavoured.
7 Combined results in the ( , M ) plane The h particle is clearly very different from any other known 'fundamental' particle. The fact that it decays into γγ implies that it cannot have spin one, and hence cannot be a gauge boson. It may well have spin zero: this remains to be demonstrated, though this hypothesis has been used in the W W * and ZZ * event selections. If it has spin 2, that would make it an even more remarkable discovery. If it does have spin zero, there is no reason why its couplings to different fermion generations (for example) should be universal, and the Higgs hypothesis suggests that its couplings to other particles should be proportional to their masses.
We now discuss the light on the nature of the h particle that is cast by figures 3, 4 and 5. In particular, figure 3 suggests that the data are heading straight towards the Higgs 'bull's eye' at the cross-hairs where M ∼ v = 246 GeV and ∼ 0, corresponding to couplings scaling with masses. As we see in the right panel of figure  At first sight, one might be surprised that it is already possible to obtain such a tight constraint on . The essential reason is that, because it is so much heavier than all the other fermions, the coupling to the top quark provides a long lever arm, and similarly for the W ± and Z because they are also much heavier than the other fermions. As already commented in section 1, figure 5 provides another way of visualizing this observation. The diagonal line in figure 5 that represents the mass dependence of the Higgs couplings expected in the Standard Model (2.6) is completely compatible within errors with the measurements. We consider figures 3, 4 and 5 to be the most remarkable results of our analysis.

Overview and prospects
A new particle has been discovered: how closely does it resemble the Higgs boson of the Standard Model? In this paper we have presented a global analysis of the data from the Tevatron experiments [29,30] CDF, D0, ATLAS [17,18] and CMS [15,16] made available before the ICHEP 2012 conference, making two types of fit. One is in terms of universal coefficients (a, c) that parametrize the deviations of the h couplings to fermions and bosons JHEP09(2012)123 in a way well adapted to constraining composite Higgs models such as the pseudo-dilaton, MCHM4 and MCHM5 models, as well as the 'anti-dilaton', fermiophobic and gaugephobic scenarios. As seen in figures 1 and 2, the only models favoured in this fit are pseudodilaton and MCHM4 models with parameters close to the Standard Model case. We have also made a fit with h couplings to fermions and bosons scaling as some power 1 + of the particle masses, with a normalization scale M = v in general. As seen in figures 3, 4 and 5, this fit favours ∼ 0 and M ∼ v, as expected for the Standard Model Higgs boson. As seen in figure 9, the overall quality of a fit to the Standard Model Higgs boson is good, and does not depend strongly on m h . This plot was made by calculating the h production cross-sections and decay branching rates assuming Standard Model couplings, i.e., a = c = 1, = 0, M = v, while leaving m h as a free parameter. Note that information on the shapes of the h signal in, e.g., the high-resolution γγ and ZZ * channels as functions of m h was not used in this exercise. The value of the global χ 2 function at the minimum, namely 34.1, is comparable to the number of degrees of freedom in the fit. Better understanding of the correlations in the data are needed, but the overall quality of the Standard Model Higgs fit is clearly good. We also see in figure 9 that the quality of this fit does not vary significantly over the range [124, 127] GeV, which brackets the central values of m h found in the high-resolution γγ and ZZ * channels by CMS and ATLAS. Within this range, our other results are insensitive to the value m h = 125 GeV assumed in our global fits.
We anticipate that the LHC experiments will be able to constrain the h couplings significantly further in the coming months, with improved analyses of the channels already studied, analyses of more channels using the data accumulated so far, and the prospect JHEP09(2012)123 of more data on the way. We expect that these improvements will enable the ranges of parameters in simple two-parameter fits such as those presented here to be reduced to the 10% level [102][103][104][105][106]. The upcoming data should also make possible more detailed fits incorporating more parameters, leading eventually to individual determinations of the h couplings to different bosons and fermions. In this way, we shall see whether the indication of h couplings depending linearly on other particle masses seen in figures 3, 4 and 5 will be confirmed.
So far, the h particle does indeed walk and quack very much like a Higgs boson. [5] ATLAS collaboration, Search for the standard model Higgs boson produced in association with a vector boson and decaying to a b-quark pair using up to 4.7 fb −1 pp collision data at √ s = 7 TeV with the ATLAS detector at the LHC, ATLAS-CONF-2012-015 (2012).