CERN-PH-EP-2012-218 Submitted to: Physics Letters B
The ATLAS Collaboration
This paper is dedicated to the memory of our ATLAS colleagues who did not live to see the full impact and significance of their contributions to the experiment.
Abstract
A search for the Standard Model Higgs boson in proton-proton collisions with the ATLAS detector at the LHC is presented. The datasets used correspond to integrated luminosities of approximately
√√
4.8 fb−1 collected at s = 7 TeV in 2011 and 5.8 fb−1 at s = 8 TeV in 2012. Individual searches in the channels H→ ZZ(∗)→ 4e, H→ γγ and H→ WW(∗)→ eνµν in the 8 TeV data are combined with previously published results of searches for H→ ZZ(∗), WW(∗), bb¯and τ+τ− in the 7 TeV data and results from improved analyses of the H→ ZZ(∗)→ 4e and H→ γγ channels in the 7 TeV data. Clear evidence for the production of a neutral boson with a measured mass of 126.0 ± 0.4 (stat) ± 0.4 (sys) GeV is presented. This observation, which has a significance of 5.9 standard deviations, corresponding to a background fluctuation probability of 1.7 ×10−9, is compatible with the production and decay of the Standard Model Higgs boson.
The ATLAS Collaboration
This paper is dedicated to the memory of our ATLAS colleagues who did not live to see the full impact and significance of their contributions to the experiment.
Abstract
A search for the Standard Model Higgs boson in proton-proton collisions with the ATLAS detector at the LHC is presented. The datasets used correspond to integrated luminosities of approximately 4.8 fb−1 collected at √ s = 7 TeV in 2011 and 5.8 fb−1 at √ s = 8 TeV in 2012. Individual searches in the channels H→ ZZ(∗)→ 4e, H→ γγ and H→ WW(∗)→ eνµν in the 8 TeV data are combined with previously published results of searches for H→ ZZ(∗), WW(∗), bb¯and τ+τ− in the 7 TeV data and results from improved analyses of the H→ ZZ(∗)→ 4e and H→ γγ channels in the 7 TeV data. Clear evidence for the production of a neutral boson with a measured mass of
126.0 ± 0.4 (stat) ± 0.4 (sys) GeV is presented. This observation, which has a significance of 5.9 standard deviations, corresponding to a background fluctuation probability of 1.7 × 10−9, is compatible with the production and decay of the Standard Model Higgs boson.
1. Introduction
The Standard Model (SM) of particle physics [1–4] has been tested by many experiments over the last four decades and has been shown to successfully describe high energy particle interactions. However, the mechanism that breaks electroweak symmetry in the SM has not been verified experimentally. This mechanism [5– 10], which gives mass to massive elementary particles, implies the existence of a scalar particle, the SM Higgs boson. The search for the Higgs boson, the only elementary particle in the SM that has not yet been observed, is one of the highlights of the Large Hadron Col lider [11] (LHC) physics programme.
Indirect limits on the SM Higgs boson mass of mH < 158 GeV at 95% confidence level (CL) have been set using global fits to precision electroweak results [12]. Direct searches at LEP [13], the Tevatron [14–16] and the LHC [17, 18] have previously excluded, at 95% CL, a SM Higgs boson with mass below 600 GeV, apart from some mass regions between 116 GeV and 127 GeV.
Both the ATLAS and CMS Collaborations reported excesses of events in their 2011 datasets of proton
√
proton (pp) collisions at centre-of-mass energy s = 7 TeV at the LHC, which were compatible with SM Higgs boson production and decay in the mass region 124–126 GeV, with significances of 2.9 and 3.1 standard deviations (σ ), respectively [17, 18]. The CDF and DØ experiments at the Tevatron have also recently reported a broad excess in the mass region 120–135 GeV; using the existing LHC constraints, the observed local significances for mH = 125 GeV are 2.7 σ for CDF [14], 1.1 σ for DØ [15] and 2.8 σ for their combination [16].
The previous ATLAS searches in 4.6–4.8 fb−1 of data
√
at s = 7 TeV are combined here with new searches for
H→ ZZ(∗)→ 4e, H→ γγ and H→ WW(∗)→√eνµν in the 5.8–5.9 fb−1 of pp collision data taken at s = 8 TeV between April and June 2012. In the H→ WW(∗)→ eνµν channel, the kinematic region in which a SM Higgs boson with a mass between 110 GeV and 140 GeV is searched for was kept blinded until satisfactory agreement was found between the observed and predicted number of events in control samples dominated by the principal backgrounds. The analyses of the H→ ZZ(∗)→ 4e and H→ γγ channels were re-optimised with simulation and frozen before looking at the 8 TeV data.
The data were recorded with instantaneous luminosi
−2
ties up to 6.8 × 1033 cms−1; they are therefore affected by multiple pp collisions occurring in the same or neighbouring bunch crossings (pile-up). In the 7 TeV
1The symbol e stands for electron or muon.
data, the average number of interactions per bunch crossing was approximately 10; the average increased to approximately 20 in the 8 TeV data. The reconstruction, identification and isolation criteria used for electrons and photons in the 8 TeV data are improved, making the H→ ZZ(∗)→ 4e and H→ γγ analyses more robust against the increased pile-up. In the H→ WW(∗)→ eνeν channel, the increased pile-up deteriorates the event missing transverse momentum, Emiss, resolution, which
T
results in significantly larger Drell-Yan background in the same-flavour final states. Since the eµ final state provides most of the sensitivity of the search, only this final state has been used in the analysis of the 8 TeV data.
This letter is organised as follows. The ATLAS de tector is briefly described in Section 2. The simulation samples and the signal predictions are presented in Section 3. The analyses of the H→ ZZ(∗)→ 4e, H→ γγ and H→ WW(∗)→ eνµν channels are described in Sec tions 4–6, respectively. The statistical procedure used to analyse the results is summarised in Section 7. The systematic uncertainties which are correlated between datasets and search channels are described in Section 8. The results of the combination of all channels are re ported in Section 9, while Section 10 provides the con clusions.
2. The ATLAS detector
The ATLAS detector [19–21] is a multipurpose parti cle physics apparatus with forward-backward symmetric cylindrical geometry. The inner tracking detector (ID) consists of a silicon pixel detector, a silicon microstrip detector (SCT), and a straw-tube transition radiation tracker (TRT). The ID is surrounded by a thin superconducting solenoid which provides a 2 T magnetic field, and by high-granularity liquid-argon (LAr) sampling electromagnetic calorimetry. The electromagnetic calorimeter is divided into a central barrel (pseudorapidity 2 |η| < 1.475) and end-cap regions on either end of the detector (1.375 < |η| < 2.5 for the outer wheel and 2.5 < |η| < 3.2 for the inner wheel). In the region matched to the ID (|η| < 2.5), it is radially segmented
2ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector, and the z-axis along the beam line. The x-axis points from the IP to the centre of the LHC ring, and the y-axis points upwards. Cylindrical coordinates (r,φ) are used in the transverse plane, φ being the azimuthal angle around the beam line. Observables labeled “transverse” are projected into the x − y plane. The pseudorapidity is defined in terms of the polar angle θ as η = − ln tan(θ/2).
into three layers. The first layer has a fine segmentation in η to facilitate e/γ separation from π0 and to improve the resolution of the shower position and direction measurements. In the region |η| < 1.8, the electromagnetic calorimeter is preceded by a presampler detector to correct for upstream energy losses. An iron-scintillator/tile calorimeter gives hadronic coverage in the central rapidity range (|η| < 1.7), while a LAr hadronic end-cap calorimeter provides coverage over 1.5 < |η| < 3.2. The forward regions (3.2 < |η| < 4.9) are instrumented with LAr calorimeters for both electromagnetic and hadronic measurements. The muon spectrometer (MS) surrounds the calorimeters and consists of three large air-core superconducting magnets providing a toroidal field, each with eight coils, a system of precision tracking chambers, and fast detectors for triggering. The combination of all these systems provides charged particle measurements together with efficient and precise lepton and photon measurements in the pseudorapidity range |η| < 2.5. Jets and Emiss are reconstructed using en-
T
ergy deposits over the full coverage of the calorimeters, |η| < 4.9.
3. Signal and background simulation samples
The SM Higgs boson production processes considered in this analysis are the dominant gluon fusion (gg → H, denoted ggF), vector-boson fusion
'
(qq→ qq'H, denoted VBF) and Higgs-strahlung
'
(qq→ WH, ZH, denoted WH/ZH). The small contribution from the associated production with a tt pair (qq¯/gg → t¯tH) is taken into account only
tH, denoted t¯in the H→ γγ analysis.
For the ggF process, the signal cross section is computed at up to next-to-next-to-leading order (NNLO) in QCD [22–28]. Next-to-leading order (NLO) elec troweak (EW) corrections are applied [29, 30], as well as QCD soft-gluon re-summations at up to next-to-next to-leading logarithm (NNLL) [31]. These calculations, which are described in Refs. [32–35], assume factori sation between QCD and EW corrections. The transverse momentum, pT, spectrum of the Higgs boson in the ggF process follows the HqT calculation [36], which includes QCD corrections at NLO and QCD soft-gluon re-summations up to NNLL; the effects of finite quark masses are also taken into account [37].
For the VBF process, full QCD and EW correc tions up to NLO [38–41] and approximate NNLO QCD corrections [42] are used to calculate the cross sec tion. Cross sections of the associated WH/ZH processes (VH) are calculated including QCD corrections up to NNLO [43–45] and EW corrections up to NLO [46].
The cross sections for the t¯
tH process are estimated up to NLO QCD [47–51]. The total cross sections for SM Higgs boson production at the LHC with mH = 125 GeV are predicted to
√√
be 17.5 pb for s = 7 TeV and 22.3 pb for s = 8 TeV [52, 53].
The branching ratios of the SM Higgs boson as a function of mH, as well as their uncertainties, are calcu lated using the HDECAY [54] and PROPHECY4F [55, 56 ] programs and are taken from Refs. [52, 53]. The interference in the H→ ZZ(∗)→ 4e final states with iden tical leptons is taken into account [53, 55, 56].
Table 1: Event generators used to model the signal and background processes. “PYTHIA” indicates that PYTHIA6 and PYTHIA8 are
√√
used for simulations of s = 7 TeV and s = 8 TeV data, respectively.
Process Generator
ggF, VBF POWHEG [57, 58] +PYTHIA WH, ZH, t¯PYTHIA
tH
W+jets, Z/γ∗ +jets ALPGEN [59] +HERWIG
tt, tW, tb MC@NLO [60] +HERWIG
tqb AcerMC [61] +PYTHIA
qq¯→ WW MC@NLO+HERWIG
gg → WW gg2WW [62] +HERWIG
qq¯→ ZZ POWHEG [63] +PYTHIA
gg → ZZ gg2ZZ [64] +HERWIG
WZ MadGraph+PYTHIA, HERWIG
Wγ+jets ALPGEN+HERWIG
Wγ∗ [65] MadGraph+PYTHIA
qq¯/gg → γγ SHERPA
The event generators used to model signal and background processes in samples of Monte Carlo (MC) sim ulated events are listed in Table 1. The normalisations of the generated samples are obtained from the state of the art calculations described above. Several different programs are used to generate the hard-scattering processes. To generate parton showers and their hadroni sation, and to simulate the underlying event [66–68], PYTHIA6 [69] (for 7 TeV samples and 8 TeV sam ples produced with MadGraph [70, 71] or AcerMC) or PYTHIA8 [72] (for other 8 TeV samples) are used. Al ternatively, HERWIG [73] or SHERPA [74] are used to generate and hadronise parton showers, with the HERWIG underlying event simulation performed using JIMMY [75]. When PYTHIA6 or HERWIG are used, TAUOLA [76] and PHOTOS [77] are employed to de scribe tau lepton decays and additional photon radiation from charged leptons, respectively.
The following parton distribution function (PDF) sets are used: CT10 [78] for the POWHEG, MC@NLO, SHERPA, gg2WW and gg2ZZ samples; CTEQ6L1 [79] for the ALPGEN, MadGraph and HERWIG samples; and MRSTMCal [80] for the PYTHIA6, PYTHIA8 and AcerMC samples.
Acceptances and efficiencies are obtained mostly from full simulations of the ATLAS detector [81] us ing Geant 4 [82]. These simulations include a realistic modelling of the pile-up conditions observed in the data. Corrections obtained from measurements in data are applied to account for small differences between data and simulation (e.g. large samples of W, Z and J/ψ decays are used to compare lepton reconstruction and identification efficiencies).
4. H → ZZ(∗) → 4t channel
The search for the SM Higgs boson through the decay H → ZZ(∗) → 4e, where e = e or µ, provides good sensitivity over a wide mass range (110600 GeV), largely due to the excellent momentum resolution of the ATLAS detector. This analysis searches for Higgs boson candidates by selecting two pairs of isolated leptons, each of which is comprised of two leptons with the same flavour and opposite charge. The expected cross section times branching ratio for the process H → ZZ(∗) → 4e with mH = 125 GeV is 2.2 fb for
√√
s = 7 TeV and 2.8 fb for s = 8 TeV.
The largest background comes from continuum (Z(∗)/γ∗)(Z(∗)/γ∗) production, referred to hereafter as ZZ(∗). For low masses there are also important background contributions from Z + jets and tt¯production, where charged lepton candidates arise either from decays of hadrons with b-or c-quark content or from misidentification of jets.
The 7 TeV data have been re-analysed and combined with the 8 TeV data. The analysis is improved in several aspects with respect to Ref. [83] to enhance the sensitiv ity to a low-mass Higgs boson. In particular, the kinematic selections are revised, and the 8 TeV data analysis benefits from improvements in the electron reconstruction and identification. The expected signal significances for a Higgs boson with mH = 125 GeV are
1.6 σ for the 7 TeV data (to be compared with 1.25 σ in Ref. [83]) and 2.1 σ for the 8 TeV data.
The data are selected using single-lepton or dilepton triggers. For the single-muon trigger, the pT threshold is 18 GeV for the 7 TeV data and 24 GeV for the 8 TeV data, while for the single-electron trigger the transverse energy, ET, threshold varies from 20 GeV to 22 GeV for the 7 TeV data and is 24 GeV for the 8 TeV data. For the dielectron triggers, the thresholds are 12 GeV for both electrons. For the dimuon triggers, the thresholds for the 7 TeV data are 10 GeV for each muon, while for the 8 TeV data the thresholds are 13 GeV. An additional asymmetric dimuon trigger is used in the 8 TeV data with thresholds 18 GeV and 8 GeV for the leading and sub-leading muon, respectively.
Muon candidates are formed by matching reconstructed ID tracks with either a complete track or a track-segment reconstructed in the MS [84]. The muon acceptance is extended with respect to Ref. [83] us ing tracks reconstructed in the forward region of the MS (2.5 < |η| < 2.7), which is outside the ID coverage. If both an ID and a complete MS track are present, the two independent momentum measurements are combined; otherwise the information of the ID or the MS is used alone. Electron candidates must have a well-reconstructed ID track pointing to an electromagnetic calorimeter cluster and the cluster should satisfy a set of identification criteria [85] that require the lon gitudinal and transverse shower profiles to be consistent with those expected for electromagnetic showers. Tracks associated with electromagnetic clusters are fit ted using a Gaussian-Sum Filter [86], which allows for bremsstrahlung energy losses to be taken into account.
Each electron (muon) must satisfy pT > 7 GeV (pT > 6 GeV) and be measured in the pseudorapidity range |η| < 2.47 (|η| < 2.7). All possible quadruplet combinations with same-flavour opposite-charge lepton pairs are then formed. The most energetic lepton in the quadruplet must satisfy pT > 20 GeV, and the second (third) lepton in pT order must satisfy pT > 15 GeV (pT > 10 GeV). At least one of the leptons must satisfy the single-lepton trigger or one pair must satisfy the dilepton trigger requirements. The leptons are required to be separated from each other by ΔR =
r
(Δη)2 + (Δφ)2 > 0.1 if they are of the same flavour and by ΔR > 0.2 otherwise. The longitudinal impact parameters of the leptons along the beam axis are required to be within 10 mm of the reconstructed primary vertex. The primary vertex used for the event is defined as the
reconstructed vertex with the highestp2 of associated
T
tracks and is required to have at least three tracks with pT > 0.4 GeV. To reject cosmic rays, muon tracks are required to have a transverse impact parameter, defined as the distance of closest approach to the primary vertex in the transverse plane, of less than 1 mm.
The same-flavour and opposite-charge lepton pair with an invariant mass closest to the Z boson mass (mZ ) in the quadruplet is referred to as the leading lepton pair. Its invariant mass, denoted by m12, is required to be between 50 GeV and 106 GeV. The remaining same-flavour, opposite-charge lepton pair is the sub-leading lepton pair. Its invariant mass, m34, is required to be in the range mmin < m34 < 115 GeV, where the value of mmin depends on the reconstructed four-lepton invariant mass, m4e. The value of mmin varies monotonically from 17.5 GeV at m4e = 120 GeV to 50 GeV for m4e> 190 GeV [87]. All possible lepton pairs in the quadruplet that have the same flavour and opposite charge must satisfy mee > 5 GeV in order to reject backgrounds involving the production and decay of J/ψ mesons. If two or more quadruplets satisfy the above selection, the one with the highest value of m34 is selected. Four different analysis sub-channels, 4e,2e2µ, 2µ2e and 4µ, arranged by the flavour of the leading lepton pair, are defined.
Non-prompt leptons from heavy flavour decays, electrons from photon conversions and jets mis-identified as electrons have broader transverse impact parameter distributions than prompt leptons from Z boson decays and/or are non-isolated. Thus, the Z+jets and tt¯background contributions are reduced by applying a cut on the transverse impact parameter significance, defined as the transverse impact parameter divided by its uncertainty, d0/σd0. This is required to be less than 3.5 (6.5) for muons (electrons). The electron impact parameter is affected by bremsstrahlung and thus has a broader distribution.
In addition, leptons must satisfy isolation requirements based on tracking and calorimetric information. The normalised track isolation discriminant is defined as the sum of the transverse momenta of tracks inside a cone of size ΔR = 0.2 around the lepton direction, excluding the lepton track, divided by the lepton pT. The tracks considered in the sum are those compatible with the lepton vertex and have pT > 0.4 GeV (pT > 1 GeV) in the case of electron (muon) candidates. Each lepton is required to have a normalised track isolation smaller than 0.15. The normalised calorimetric isolation for electrons is computed as the sum of the ET of positive- energy topological clusters [88] with a reconstructed barycentre falling within a cone of size ΔR = 0.2 around the candidate electron cluster, divided by the electron ET. The summed energy of the cells assigned to the electron cluster is excluded, while a correction is applied to account for the electron energy deposited outside the cluster. The algorithm for topological clustering suppresses noise by keeping cells with a significant energy deposit and their neighbours. The ambient energy deposition in the event from pile-up and the underlying event is accounted for using a calculation of the median transverse energy density from low-pT jets [89, 90]. The normalised calorimetric isolation for electrons is required to be less than 0.20. The normalised calorimetric isolation discriminant for muons is defined by the ratio to the pT of the muon of the ET sum of the calorimeter cells inside a cone of size ΔR = 0.2 around the muon direction minus the energy deposited by the muon. Muons are required to have a normalised calorimetric isolation less than 0.30 (0.15 for muons without an associated ID track). For both the track-and calorimeter-based isolation, any contributions arising from other leptons of the quadruplet are subtracted.
The combined signal reconstruction and selection efficiencies for a SM Higgs with mH = 125 GeV for the 7 TeV (8 TeV) data are 37% (36%) for the 4µ channel, 20% (22%) for the 2e2µ/2µ2e channels and 15% (20%) for the 4e channel.
The 4e invariant mass resolution is improved by applying a Z-mass constrained kinematic fit to the leading lepton pair for m4e< 190 GeV and to both lepton pairs for higher masses. The expected width of the reconstructed mass distribution is dominated by the experimental resolution for mH < 350 GeV, and by the natural width of the Higgs boson for higher masses (30 GeV at mH = 400 GeV). The typical mass resolutions for mH = 125 GeV are 1.7 GeV, 1.7 GeV/2.2 GeV and
2.3 GeV for the 4µ,2e2µ/2µ2e and 4e sub-channels, respectively.
The expected background yield and composition are estimated using the MC simulation normalised to the theoretical cross section for ZZ(∗) production and by methods using control regions from data for the Z + jets and tt¯processes. Since the background composition depends on the flavour of the sub-leading lepton pair, different approaches are taken for the ee+µµ and the ee+ee final states. The transfer factors needed to extrapolate the background yields from the control regions defined below to the signal region are obtained from the MC simulation. The MC description of the selection efficiencies for the different background components has been verified with data.
The reducible ee + µµ background is dominated by tt¯and Z + jets (mostly Zbb¯) events. A control region is defined by removing the isolation requirement on the leptons in the sub-leading pair, and by requiring that at least one of the sub-leading muons fails the transverse impact parameter significance selection. These modifications remove ZZ(∗) contributions, and allow both the tt¯and Z + jets backgrounds to be estimated simultaneously using a fit to the m12 distribution. The tt¯background contribution is cross-checked by selecting a control sample of events with an opposite charge eµ pair with an invariant mass between 50 GeV and 106 GeV, accompanied by an opposite-charge muon pair. Events with a Z candidate decaying to a pair of electrons or muons in the aforementioned mass range are excluded. Isolation and transverse impact parameter significance requirements are applied only to the leptons of the eµ pair.
In order to estimate the reducible ee + ee background, a control region is formed by relaxing the selection criteria for the electrons of the sub-leading pair. The different sources of electron background are then separated into categories consisting of non-prompt leptons from heavy flavour decays, electrons from photon conversions and jets mis-identified as electrons, using appro priate discriminating variables [91]. This method allows the sum of the Z + jets and tt¯background contributions to be estimated. As a cross-check, the same method is also applied to a similar control region containing same-charge sub-leading electron pairs. An additional crosscheck of the ee + ee background estimation is performed by using a control region with same-charge sub-leading electron pairs, where the three highest pT leptons satisfy all the analysis criteria whereas the selection cuts are relaxed for the remaining electrons. All the cross-checks yield consistent results.
Table 2: Summary of the estimated numbers of Z + jets and tt¯back
√√
ground events, for the s = 7 TeV and s = 8 TeV data in the entire phase-space of the analysis after the kinematic selections described in the text. The backgrounds are combined for the 2µ2e and 4e channels, as discussed in the text. The first uncertainty is statistical, while the second is systematic.
Background Estimated numbers of events
√√
s = 7 TeV s = 8 TeV
4µ
Z+jets 0.3± 0.1 ±0.1 0.5± 0.1 ±0.2 tt¯0.02±0.02±0.01 0.04±0.02±0.02
2e2µ
Z+jets 0.2± 0.1 ±0.1 0.4± 0.1 ±0.1 tt¯0.02±0.01±0.01 0.04±0.01±0.01
2µ2e
Z+jets, t¯t | 2.6± 0.4 ±0.4 | 4.9± 0.8 ±0.7 | |
4e | |||
Z+jets, t¯t | 3.1± 0.6 ±0.5 | 3.9± 0.7 ±0.8 |
The data-driven background estimates are sum marised in Table 2. The distribution of m34, for events selected by the analysis except that the isolation and transverse impact parameter requirements for the sub leading lepton pair are removed, is presented in Fig. 1.
Figure 1: Invariant mass distribution of the sub-leading lepton pair (m34) for a sample defined by the presence of a Z boson candidate and an additional same-flavour electron or muon pair, for the combination
√√
of s = 7 TeV and s = 8 TeV data in the entire phase-space of the analysis after the kinematic selections described in the text. Isolation and transverse impact parameter significance requirements are applied to the leading lepton pair only. The MC is normalised to the data-driven background estimations. The relativelly small contribution of a SM Higgs with mH = 125 GeV in this sample is also shown.
The uncertainties on the integrated luminosities are determined to be 1.8% for the 7 TeV data and 3.6% for the 8 TeV data using the techniques described in Ref. [92].
The uncertainties on the lepton reconstruction and identification efficiencies and on the momentum scale and resolution are determined using samples of W, Z and J/ψ decays [84, 85]. The relative uncertainty on the signal acceptance due to the uncertainty on the muon reconstruction and identification efficiency is ±0.7% (±0.5%/±0.5%) for the 4µ (2e2µ/2µ2e) channel for m4e = 600 GeV and increases to ±0.9% (±0.8%/±0.5%) for m4e = 115 GeV. Similarly, the relative uncertainty on the signal acceptance due to the uncertainty on the electron reconstruction and identification efficiency is ±2.6% (±1.7%/±1.8%) for the 4e (2e2µ/2µ2e) channel for m4e = 600 GeV and reaches ±8.0% (±2.3%/±7.6%) for m4e = 115 GeV. The uncertainty on the electron energy scale results in an uncertainty of ±0.7% (±0.5%/±0.2%) on the mass scale of the m4e distribution for the 4e (2e2µ/2µ2e) channel. The impact of the uncertainties on the electron energy resolution and on the muon momentum resolution and scale are found to be negligible.
The theoretical uncertainties associated with the sig nal are described in detail in Section 8. For the SM ZZ(∗) background, which is estimated from MC simulation, the uncertainty on the total yield due to the QCD scale uncertainty is ±5%, while the effect of the PDF and αs uncertainties is ±4% (±8%) for processes initi ated by quarks (gluons) [53]. In addition, the dependence of these uncertainties on the four-lepton invariant mass spectrum has been taken into account as discussed in Ref. [53]. Though a small excess of events is ob served for m4l > 180 GeV, the measured ZZ(∗) → 4e cross section [93] is consistent with the SM theoreti cal prediction. The impact of not using the theoretical constraints on the ZZ(∗) yield on the search for a Higgs boson with mH < 2mZ has been studied in Ref. [87] and has been found to be negligible . The impact of the interference between a Higgs signal and the non-resonant gg → ZZ(∗) background is small and becomes negligible for mH < 2mZ [94].
Figure 2: The distribution of the four-lepton invariant mass, m4e , for the selected candidates, compared to the background expectation in
√
the 80–250 GeV mass range, for the combination of the s = 7 TeV
√
and s = 8 TeV data. The signal expectation for a SM Higgs with mH = 125 GeV is also shown.
4.4. Results The expected distributions of m4e for the background and for a Higgs boson signal with mH = 125 GeV are compared to the data in Fig. 2. The numbers of ob
served and expected events in a window of ±5 GeV around mH = 125 GeV are presented for the combined 7 TeV and 8 TeV data in Table 3. The distribution of the m34 versus m12 invariant mass is shown in Fig. 3. The statistical interpretation of the excess of events near m4e = 125 GeV in Fig. 2 is presented in Section 9.
Table 3: The numbers of expected signal (mH = 125 GeV) and background events, together with the numbers of observed events in the data, in a window of size ±5 GeV around 125 GeV, for the combined
√√
s = 7 TeV and s = 8 TeV data.
Signal ZZ(∗) Z + jets, tt¯Observed
4µ | 2.09±0.30 | 1.12±0.05 | 0.13±0.04 | 6 |
---|---|---|---|---|
2e2µ/2µ2e | 2.29± 0.33 | 0.80±0.05 | 1.27±0.19 | 5 |
4e | 0.90±0.14 | 0.44±0.04 | 1.09±0.20 | 2 |
Figure 3: Distribution of the m34 versus the m12 invariant mass, before the application of the Z-mass constrained kinematic fit, for the selected candidates in the m4e range 120–130 GeV. The expected distributions for a SM Higgs with mH = 125 GeV (the sizes of the boxes indicate the relative density) and for the total background (the intensity of the shading indicates the relative density) are also shown.
5. H→ γγ channel
The search for the SM Higgs boson through the decay H→ γγ is performed in the mass range between 110 GeV and 150 GeV. The dominant background is SM diphoton production (γγ); contributions also come from γ+jet and jet+jet production with one or two jets mis-identified as photons (γ j and jj) and from the Drell-Yan process. The 7 TeV data have been re-analysed and the results combined with those from the 8 TeV data. Among other changes to the analysis, a new category of events with two jets is introduced, which enhances the sensitivity to the VBF process. Higgs boson events produced by the VBF process have two forward jets, originating from the two scattered quarks, and tend to be devoid of jets in the central region. Overall, the sensitivity of the analysis has been improved by about 20% with respect to that described in Ref. [95].
The data used in this channel are selected using a diphoton trigger [96], which requires two clusters formed from energy depositions in the electromagnetic calorimeter. An ET threshold of 20 GeV is applied to each cluster for the 7 TeV data, while for the 8 TeV data the thresholds are increased to 35 GeV on the leading (the highest ET) cluster and to 25 GeV on the sub-leading (the next-highest ET) cluster. In addition, loose criteria are applied to the shapes of the clusters to match the expectations for electromagnetic showers initiated by photons. The efficiency of the trigger is greater than 99% for events passing the final event selection.
Events are required to contain at least one reconstructed vertex with at least two associated tracks with pT > 0.4 GeV, as well as two photon candidates. Photon candidates are reconstructed in the fiducial region |η| < 2.37, excluding the calorimeter barrel/end-cap transition region 1.37 ≤|η| < 1.52. Photons that convert to electron-positron pairs in the ID material can have one or two reconstructed tracks matched to the clusters in the calorimeter. The photon reconstruction efficiency is about 97% for ET > 30 GeV.
In order to account for energy losses upstream of the calorimeter and energy leakage outside of the cluster, MC simulation results are used to calibrate the energies of the photon candidates; there are separate calibrations for unconverted and converted candidates. The calibration is refined by applying η-dependent correction factors, which are of the order of ±1%, determined from measured Z→ e+e− events. The leading (sub-leading) photon candidate is required to have ET > 40 GeV (30 GeV).
Photon candidates are required to pass identification criteria based on shower shapes in the electromagnetic calorimeter and on energy leakage into the hadronic calorimeter [97]. For the 7 TeV data, this information is combined in a neural network, tuned to achieve a similar jet rejection as the cut-based selection described in Ref. [95], but with higher photon e fficiency. For the 8 TeV data, cut-based criteria are used to ensure reliable photon performance for recently-recorded data. This cut-based selection has been tuned to be robust against pile-up by relaxing requirements on shower shape criteria more susceptible to pile-up, and tightening others.
The photon identification efficiencies, averaged over η, range from 85% to above 95% for the ET range under consideration.
To further suppress the jet background, an isolation requirement is applied. The isolation transverse energy is defined as the sum of the transverse energy of positive-energy topological clusters, as described in Section 4, within a cone of size ΔR = 0.4 around the photon candidate, excluding the region within 0.125 ×
0.175 in Δη×Δφ around the photon barycentre. The distributions of the isolation transverse energy in data and simulation have been found to be in good agreement using electrons from Z→ e+e− events and photons from
+
Z → ee−γ events. Remaining small differences are taken into account as a systematic uncertainty. Photon candidates are required to have an isolation transverse energy of less than 4 GeV.
The invariant mass of the two photons is evaluated using the photon energies measured in the calorimeter, the azimuthal angle φ between the photons as determined from the positions of the photons in the calorimeter, and the values of η calculated from the position of the identified primary vertex and the impact points of the photons in the calorimeter.
The primary vertex of the hard interaction is identified by combining the following information in a global likelihood: the directions of flight of the photons as determined using the longitudinal segmentation of the electromagnetic calorimeter (calorimeter pointing), the
2
parameters of the beam spot, and the pT of the tracks associated with each reconstructed vertex. In addition, for the 7 TeV data analysis, the reconstructed conversion vertex is used in the likelihood for converted photons with tracks containing hits in the silicon layers of the ID. The calorimeter pointing is sufficient to ensure that the contribution of the opening angle between the photons to the mass resolution is negligible. Using the calorimeter pointing alone, the resolution of the vertex z coordinate is ∼ 15 mm, improving to ∼ 6 mm for events with two reconstructed converted photons. The tracking information from the ID improves the identification of the vertex of the hard interaction, which is needed for the jet selection in the 2-jet category.
With the selection described in Section 5.1, in the diphoton invariant mass range between 100 GeV and 160 GeV, 23788 and 35251 diphoton candidates are observed in the 7 TeV and 8 TeV data samples, respectively.
Data-driven techniques [98] are used to estimate the numbers of γγ, γ j and jj events in the selected sample. The contribution from the Drell-Yan background is determined from a sample of Z→ e+e− decays in data where either one or both electrons pass the photon selection. The measured composition of the selected sample is approximately 74%, 22%, 3% and 1% for the γγ, γ j, jj and Drell-Yan processes, respectively, demonstrating the dominance of the irreducible diphoton production. This decomposition is not directly used in the signal search; however, it is used to study the parameterisation of the background modelling.
To increase the sensitivity to a Higgs boson signal, the events are separated into ten mutually exclusive categories having different mass resolutions and signal-tobackground ratios. An exclusive category of events containing two jets improves the sensitivity to VBF. The other nine categories are defined by the presence or not of converted photons, η of the selected photons, and pTt, the component 3 of the diphoton pT that is orthogonal to the axis defined by the difference between the two pho ton momenta [99, 100].
Jets are reconstructed [101] using the anti- kt algo rithm [102] with radius parameter R = 0.4. At least two jets with |η| < 4.5 and pT > 25 GeV are required in the 2-jet selection. In the analysis of the 8 TeV data, the pT threshold is raised to 30 GeV for jets with
2.5 < |η| < 4.5. For jets in the ID acceptance (|η| < 2.5), the fraction of the sum of the pT of tracks, associated with the jet and matched to the selected primary vertex, with respect to the sum of the pT of tracks associated with the jet (jet vertex fraction, JVF) is required to be at least 0.75. This requirement on the JVF reduces the number of jets from proton-proton interactions not associated with the primary vertex. Motivated by the VBF topology, three additional cuts are applied in the 2-jet selection: the difference of the pseudorapidity between the leading and sub-leading jets (tag jets) is required to be larger than 2.8, the invariant mass of the tag jets has to be larger than 400 GeV, and the azimuthal angle difference between the diphoton system and the system of the tag jets has to be larger than 2.6. About 70% of the signal events in the 2-jet category come from the VBF process.
The other nine categories are defined as follows: events with two unconverted photons are separated into unconverted central (|η| < 0.75 for both candidates) and unconverted rest (all other events), events with at least
γ2 γ2 γ2 γ1
3 pTt = (p γ1 + p ) × (p γ1 − p ) / p γ1 − p , where p and p γ2
TTTTTTT T
are the transverse momenta of the two photons.
Table 4: Number of events in the data (ND) and expected number of signal events (NS) for mH = 126.5 GeV from the H→ γγ analysis, for each category in the mass range 100−160 GeV. The mass resolution FWHM (see text) is also given for the 8 TeVdata. The Higgs boson production cross section multiplied by the branching ratio into two photons (σ×B(H → γγ)) is listed for mH = 126.5 GeV. The statistical uncertainties on NS and FWHM are less than 1 %.
√
s | 7 TeV | 8 TeV | ||
σ × B(H → γγ) [fb] | 39 | 50 | FWHM [GeV] | |
Category | ND NS | ND | NS | |
Unconv. central, low pTt Unconv. central, high pTt Unconv. rest, low pTt Unconv. rest, high pTt Conv. central, low pTt Conv. central, high pTt Conv. rest, low pTt Conv. rest, high pTt Conv. transition 2-jet | 2054 10.5 97 1.5 7129 21.6 444 2.8 1493 6.7 77 1.0 8313 21.1 501 2.7 3591 9.5 89 2.2 | 2945 173 12136 785 2015 113 11099 706 5140 139 | 14.2 2.5 30.9 5.2 8.9 1.6 26.9 4.5 12.8 3.0 | 3.4 3.2 3.7 3.6 3.9 3.5 4.5 3.9 6.1 3.7 |
All categories (inclusive) | 23788 79.6 | 35251 | 110.5 | 3.9 |
one converted photon are separated into converted central (|η| < 0.75 for both candidates), converted transition (at least one photon with 1.3 < |η| < 1.75) and converted rest (all other events). Except for the converted transition category, each category is further divided by a cut at pTt= 60 GeV into two categories, low pTt and high pTt. MC studies show that signal events, particularly those produced via VBF or associated production (WH/ZH and t¯
tH), have on average larger pTt than background events. The number of data events in each category, as well as the sum of all the categories, which is denoted inclusive , are given in Table 4.
The description of the Higgs boson signal is obtained from MC, as described in Section 3. The cross sections multiplied by the branching ratio into two photons are given in Table 4 for mH = 126.5 GeV. The number of signal events produced via the ggF process is rescaled to take into account the expected destructive interference between the gg → γγ continuum background and ggF [103], leading to a reduction of the production rate by 2−5% depending on mH and the event category. For both the 7 TeV and 8 TeV MC samples, the fractions of ggF, VBF, WH, ZH and t¯
tH production are approximately 88%, 7%, 3%, 2% and 0.5%, respectively, for mH = 126.5 GeV.
In the simulation, the shower shape distributions are shifted slightly to improve the agreement with the data [97], and the photon energy resolution is broad ened (by approximately 1% in the barrel calorimeter and 1.2−2.1% in the end-cap regions) to account for small differences observed between Z→ e+e− data and MC events. The signal yields expected for the 7 TeV and 8 TeV data samples are given in Table 4. The over all selection efficiency is about 40%.
The shape of the invariant mass of the signal in each category is modelled by the sum of a Crystal Ball func tion [104], describing the core of the distribution with a width σCB, and a Gaussian contribution describing the tails (amounting to <10%) of the mass distribution. The expected full-width-at-half-maximum (FWHM) is
3.9 GeV and σCB is 1.6 GeV for the inclusive sample. The resolution varies with event category (see Table 4); the FWHM is typically a factor 2.3 larger than σCB.
The background in each category is estimated from data by fitting the diphoton mass spectrum in the mass range 100−160 GeV with a selected model with free parameters of shape and normalisation. Different models are chosen for the different categories to achieve a good compromise between limiting the size of a potential bias while retaining good statistical power. A fourth-order Bernstein polynomial function [105] is used for the unconverted rest (low pTt), converted rest (low pTt) and inclusive categories, an exponential function of a second-order polynomial for the unconverted central (low pTt), converted central (low pTt) and converted transition categories, and an exponential function for all others.
Studies to determine the potential bias have been performed using large samples of simulated background events complemented by data-driven estimates. The background shapes in the simulation have been crosschecked using data from control regions. The potential bias for a given model is estimated, separately for each category, by performing a maximum likelihood fit to large samples of simulated background events in the mass range 100−160 GeV, of the sum of a signal plus the given background model. The signal shape is taken to follow the expectation for a SM Higgs boson; the signal yield is a free parameter of the fit. The potential bias is defined by the largest absolute signal yield obtained from the likelihood fit to the simulated background samples for hypothesized Higgs boson masses in the range 110−150 GeV. A pre-selection of background parameterisations is made by requiring that the potential bias, as defined above, is less than 20% of the statistical uncertainty on the fitted signal yield. The pre-selected parameterisation in each category with the best expected sensitivity for mH = 125 GeV is selected as the background model.
The largest absolute signal yield as defined above is taken as the systematic uncertainty on the background model. It amounts to ±(0.2−4.6) and ±(0.3−6.8) events, depending on the category for the 7 TeV and 8 TeV data samples, respectively. In the final fit to the data (see Section 5.7) a signal-like term is included in the likeli hood function for each category. This term incorporates the estimated potential bias, thus providing a conservative estimate of the uncertainty due to the background modeling.
Hereafter, in cases where two uncertainties are quoted, they refer to the 7 TeV and 8 TeV data, respectively. The dominant experimental uncertainty on the signal yield (±8%, ±11%) comes from the photon reconstruction and identification efficiency, which is estimated with data using electrons from Z decays and
+
photons from Z → ee−γ events. Pile-up modelling also affects the expected yields and contributes to the uncertainty (±4%). Further uncertainties on the signal yield are related to the trigger (±1%), photon isolation (±0.4%, ±0.5%) and luminosity (±1.8%, ±3.6%). Uncertainties due to the modelling of the underlying event are ±6% for VBF and ±30% for other production processes in the 2-jet category. Uncertainties on the predicted cross sections and branching ratio are sum marised in Section 8.
The uncertainty on the expected fractions of signal events in each category is described in the following. The uncertainty on the knowledge of the material in front of the calorimeter is used to derive the amount of possible event migration between the converted and unconverted categories (±4%). The uncertainty from pileup on the population of the converted and unconverted categories is ±2%. The uncertainty from the jet energy scale (JES) amounts to up to ±19% for the 2-jet category, and up to ±4% for the other categories. Uncertainties from the JVF modelling are ±12% (for the 8 TeV data) for the 2-jet category, estimated from Z+2-jets events by comparing data and MC. Different PDFs and scale variations in the HqT calculations are used to derive possible event migration among categories (±9%) due to the modelling of the Higgs boson kinematics.
The total uncertainty on the mass resolution is ±14%. The dominant contribution (±12%) comes from the uncertainty on the energy resolution of the calorimeter, which is determined from Z→ e+e− events. Smaller contributions come from the imperfect knowledge of the material in front of the calorimeter, which affects the extrapolation of the calibration from electrons to photons (±6%), and from pile-up (±4%).
Figure 4: The distributions of the invariant mass of diphoton candidates after all selections for the combined 7 TeV and 8 TeV data sample. The inclusive sample is shown in a) and a weighted version of the same sample in c); the weights are explained in the text. The result of a fit to the data of the sum of a signal component fixed to mH = 126.5 GeV and a background component described by a fourth-order Bernstein polynomial is superimposed. The residuals of the data and weighted data with respect to the respective fitted background component are displayed in b) and d).
The distributions of the invariant mass, mγγ, of the diphoton events, summed over all categories, are shown in Fig. 4(a) and (b). The result of a fit including a signal component fixed to mH = 126.5 GeV and a background component described by a fourth-order Bernstein polynomial is superimposed.
The statistical analysis of the data employs an unbinned likelihood function constructed from those of the ten categories of the 7 TeV and 8 TeV data samples. To demonstrate the sensitivity of this likelihood analy sis, Fig. 4(c) and (d) also show the mass spectrum ob tained after weighting events with category-dependent factors reflecting the signal-to-background ratios. The weight wi for events in category i ∈ [1, 10] for the 7 TeV and 8 TeV data samples is defined to be ln (1 + Si/Bi), where Si is 90% of the expected signal for mH = 126.5 GeV, and Bi is the integral, in a window containing Si, of a background-only fit to the data. The values Si/Bi have only a mild dependence on mH .
The statistical interpretation of the excess of events near mγγ = 126. 5 GeV in Fig. 4 is presented in Sec tion 9.
6. H→ WW(∗)→ eνµν channel
The signature for this channel is two opposite-charge leptons with large transverse momentum and a large momentum imbalance in the event due to the escaping neutrinos. The dominant backgrounds are non-resonant WW, tt¯, and Wt production, all of which have real W pairs in the final state. Other important backgrounds include Drell-Yan events (pp→ Z/γ(∗)→ ee) with Emiss
T
that may arise from mismeasurement, W+jets events in which a jet produces an object reconstructed as the second electron or muon, and Wγ events in which the photon undergoes a conversion. Boson pair production (Wγ∗/WZ(∗) and ZZ(∗)) can also produce opposite-charge lepton pairs with additional leptons that are not detected.
The analysis of the 8 TeV data presented here is focused on the mass range 110 < mH < 200 GeV. It follows the procedure used for the 7 TeV data, described in Ref. [106], except that more stringent criteria are ap plied to reduce the W+jets background and some selections have been modified to mitigate the impact of the higher instantaneous luminosity at the LHC in 2012. In particular, the higher luminosity results in a larger Drell-Yan background to the same-flavour final states, due to the deterioration of the missing transverse momentum resolution. For this reason, and the fact that the eµ final state provides more than 85% of the sensitivity of the search, the same-flavour final states have not been used in the analysis described here.
For the 8 TeV H→ WW(∗)→ eνµν search, the data are selected using inclusive single-muon and single-electron triggers. Both triggers require an isolated lepton with pT > 24 GeV. Quality criteria are applied to suppress non-collision backgrounds such as cosmic-ray muons, beam-related backgrounds, and noise in the calorimeters. The primary vertex selection fol lows that described in Section 4. Candidates for the H→ WW(∗)→ eνµν search are pre-selected by requiring exactly two opposite-charge leptons of different flavours, with pT thresholds of 25 GeV for the leading lepton and 15 GeV for the sub-leading lepton. Events are classified into two exclusive lepton channels depending on the flavour of the leading lepton, where eµ (µe) refers to events with a leading electron (muon). The dilepton invariant mass is required to be greater than 10 GeV.
The lepton selection and isolation have more stringent requirements than those used for the H → ZZ(∗) → 4e analysis (see Section 4), to reduce the larger back ground from non-prompt leptons in the eνeν final state. Electron candidates are selected using a combination of tracking and calorimetric information [85]; the criteria are optimised for background rejection, at the expense of some reduced efficiency. Muon candidates are re stricted to those with matching MS and ID tracks [84], and therefore are reconstructed over |η| < 2.5. The isolation criteria require the scalar sums of the pT of charged particles and of calorimeter topological clusters within ΔR = 0.3 of the lepton direction (excluding the lepton itself) each to be less than 0.12-0.20 times the lepton pT. The exact value differs between the criteria for tracks and calorimeter clusters, for both electrons and muons, and depends on the lepton pT. Jet selec tions follow those described in Section 5.3, except that the JVF is required to be greater than 0.5.
Since two neutrinos are present in the signal final
Emiss
state, events are required to have large Emiss. is
TT
the negative vector sum of the transverse momenta of the reconstructed objects, including muons, electrons, photons, jets, and clusters of calorimeter cells not associated with these objects. The quantity Emiss used
T,rel
in this analysis is required to be greater than 25 GeV
Emiss
and is defined as: Emiss = sin Δφmin, where Δφmin
T,rel T is min(Δφ, π ), and Emiss is the magnitude of the vec
2T tor Emiss
. Here, Δφ is the angle between Emiss and the
TT transverse momentum of the nearest lepton or jet with
, Emiss
pT > 25 GeV. Compared to Emiss has increased
TT,rel
rejection power for events in which the Emiss is gener-
T
ated by a neutrino in a jet or the mismeasurement of an object, since in those events the Emiss tends to point in
T
the direction of the object. After the lepton isolation and
Emiss
requirements that define the pre-selected sample,
T,rel
the multijet background is negligible and the Drell-Yan background is much reduced. The Drell-Yan contribution becomes very small after the topological selections, described below, are applied.
The background rate and composition depend significantly on the jet multiplicity, as does the signal topology. Without accompanying jets, the signal originates almost entirely from the ggF process and the background is dominated by WW events. In contrast, when produced in association with two or more jets, the signal contains a much larger contribution from the VBF process compared to the ggF process, and the background is dominated by tt production. Therefore, to maximise the sensitivity to SM Higgs events, further selection criteria depending on the jet multiplicity are applied to the pre-selected sample. The data are subdivided into 0-jet, 1-jet and 2-jet search channels according to the number of jets in the final state, with the 2-jet channel also including higher jet multiplicities.
Owing to spin correlations in the WW(∗) system arising from the spin-0 nature of the SM Higgs boson and the V-A structure of the W boson decay vertex, the charged leptons tend to emerge from the primary ver tex pointing in the same direction [107]. This kinematic feature is exploited for all jet multiplicities by requiring that |Δφee| < 1.8, and the dilepton invariant mass, mee, be less than 50 GeV for the 0-jet and 1-jet channels. For the 2-jet channel, the mee upper bound is increased to 80 GeV.
In the 0-jet channel, the magnitude pee T of the transverse momentum of the dilepton system, pee = pT e1 +peT2,
T
is required to be greater than 30 GeV. This improves the rejection of the Drell-Yan background.
In the 1-jet channel, backgrounds from top quark production are suppressed by rejecting events containing a b-tagged jet, as determined using a b-tagging algorithm that uses a neural network and exploits the topology of weak decays of b-and c -hadrons [108]. The total trans verse momentum, ptot, defined as the magnitude of the
T T j + Emiss
vector sum ptot = pe1 + pe2 + p, is required
TTT T
to be smaller than 30 GeV to suppress top background events that have jets with pT below the threshold defined for jet counting. In order to reject the background from Z→ ττ, the ττ invariant mass, mττ, is computed under the assumptions that the reconstructed leptons are τ lepton decay products. In addition the neutrinos produced in these decays are assumed to be the only source of
Emiss
T and to be collinear with the leptons [109]. Events with |mττ − mZ | < 25 GeV are rejected if the collinear approximation yields a physical solution.
The 2-jet selection follows the 1-jet selection described above, with the ptot definition modified to in-
T
clude all selected jets. Motivated by the VBF topology, several additional criteria are applied to the tag jets, defined as the two highest-pT jets in the event. These are required to be separated in rapidity by a distance |Δyjj| > 3.8 and to have an invariant mass, mjj, larger than 500 GeV. Events with an additional jet with pT > 20 GeV between the tag jets (yj1 < y < yj2) are rejected.
A transverse mass variable, mT [110], is used to test for the presence of a signal for all jet multiplicities. This variable is defined as:
mT =(Eee + Emiss)2 −|pee + Emiss|2 ,
TT TT
where Eee =|pee |2 + m2 . The statistical analysis of
TT ee
the data uses a fit to the mT distribution in the signal region after the Δφee requirement (see Section 6.4), which results in increased sensitivity compared to the analysis described in Ref. [111].
For a SM Higgs boson with mH = 125 GeV, the cross section times branching ratio to the eνµν final
√
state is 88 fb for s = 7 TeV, increasing to 112 fb at
√
s = 8 TeV. The combined acceptance times efficiency of the 8 TeV 0-jet and 1-jet selection relative to the ggF production cross section times branching ratio is about 7.4%. The acceptance times efficiency of the 8 TeV 2-jet selection relative to the VBF production cross section times branching ratio is about 14%. Both of these figures are based on the number of events selected before the final mT criterion is applied (as described in Sec tion 6.4).
The leading backgrounds from SM processes producing two isolated high-pT leptons are WW and top (in this section, “top” background always includes both tt¯and single top, unless otherwise noted). These are estimated using partially data-driven techniques based on normalising the MC predictions to the data in control regions dominated by the relevant background source. The W+jets background is estimated from data for all jet multiplicities. Only the small backgrounds from Drell-Yan and diboson processes other than WW, as well as the WW background for the 2-jet analysis, are estimated using MC simulation.
The control and validation regions are defined by selections similar to those used for the signal region but with some criteria reversed or modified to obtain signal-depleted samples enriched in a particular background. The term “validation region” distinguishes these regions from the control regions that are used to directly normalise the backgrounds. Some control regions have significant contributions from backgrounds other than the targeted one, which introduces dependencies among the background estimates. These correlations are fully incorporated in the fit to the mT distribution. In the following sections, each background estimate is described after any others on which it depends. Hence, the largest background (WW) is described last.
The W+jets background contribution is estimated using a control sample of events where one of the two leptons satisfies the identification and isolation criteria de scribed in Section 6.1, and the other lepton fails these criteria but satisfies a loosened selection (denoted “antiidentified”). Otherwise, events in this sample are required to pass all the signal selections. The dominant contribution to this sample comes from W+jets events in which a jet produces an object that is reconstructed as a lepton. This object may be either a true electron or muon from the decay of a heavy quark, or else a product of the fragmentation identified as a lepton candidate.
The contamination in the signal region is obtained by scaling the number of events in the data control sample by a transfer factor. The transfer factor is defined here as the ratio of the number of identified lepton candidates passing all selections to the number of anti-identified leptons. It is calculated as a function of the anti-identified lepton pT using a data sample dominated by QCD jet production (dijet sample) after subtracting the residual contributions from leptons produced by leptonic W and Z decays, as estimated from data. The small remaining lepton contamination, which includes Wγ(∗)/WZ(∗) events, is subtracted using MC simulation.
The processes producing the majority of same-charge dilepton events, W+jets, Wγ(∗)/WZ(∗) and Z(∗)Z(∗), are all backgrounds in the opposite-charge signal region. W+jets and Wγ(∗) backgrounds are particularly important in a search optimised for a low Higgs boson mass hypothesis. Therefore, the normalisation and kinematic features of same-charge dilepton events are used to validate the predictions of these backgrounds. The predicted number of same-charge events after the Emiss and
T,rel
zero-jet requirements is 216 ± 7 (stat) ± 42 (syst), while 182 events are observed in the data. Satisfactory agreement between data and simulation is observed in various kinematic distributions, including those of Δφee (see Fig. 5(a)) and the transverse mass.
In the 0-jet channel, the top quark background prediction is first normalised using events satisfying the pre selection criteria described in Section 6.1. This sample is selected without jet multiplicity or b-tagging requirements, and the majority of events contain top quarks. Non-top contributions are subtracted using predictions from simulation, except for W+jets, which is estimated using data. After this normalisation is performed, the fraction of events with zero jets that pass all selections is evaluated. This fraction is small (about 3%), since the top quark decay t→ Wb has a branching ratio of nearly
Figure 5: Validation and control distributions for the H→ WW(∗)→ eνµν analysis. a) Δφee distribution in the same-charge validation region after the Emiss and zero-jet requirements. b)
T,rel
mT distribution in the WW control region for the 0-jet channel. The eµ and µe final states are combined. The hashed area indicates the total uncertainty on the background prediction. The expected signal for mH = 125 GeV is negligible and therefore not visible.
1. Predictions of this fraction from MC simulation are sensitive to theoretical uncertainties such as the modelling of initial-and final-state radiation, as well as experimental uncertainties, especially that on the jet energy scale. To reduce the impact of these uncertainties, the top quark background determination uses data from a b-tagged control region in which the one-to-two jet ra tio is compared to the MC simulation [112]. The result ing correction factor to a purely MC-based background estimate after all selections amounts to 1.11±0.06 (stat).
In the 1-jet and 2-jet analyses, the top quark background predictions are normalised to the data using control samples defined by reversing the b-jet veto and removing the requirements on Δφee and mee. The |Δyjj| and mjj requirements are included in the definition of the 2-jet control region. The resulting samples are dominated by top quark events. The small contributions from other sources are taken into account using MC simulation and the data-driven W+jets estimate. Good agreement between data and MC simulation is observed for the total numbers of events and the shapes of the mT distributions. The resulting normalisation factors are
1.11 ± 0.05 for the 1-jet control region and 1.01 ± 0.26 for the 2-jet control region. Only the statistical uncertainties are quoted.
6.2.3. WW control sample The MC predictions of the WW background in the 0-jet and 1-jet analyses, summed over lepton flavours, are normalised using control regions defined with the same selections as for the signal region except that the Δφee requirement is removed and the upper bound on mee is replaced with a lower bound: mee > 80 GeV. The numbers of events and the shape of the mT distribution in the control regions are in good agreement between data and MC, as shown in Fig. 5(b). WW production contributes about 70% of the events in the 0-jet control region and about 45% in the 1-jet region. Contaminations from sources other than WW are derived as for the signal region, including the data-driven W+jets and top estimates. The resulting normalisation factors with their associated statistical uncertainties are 1.06±0.06 for the
0-jet control region and 0.99 ± 0.15 for the 1-jet control region.
The systematic uncertainties that have the largest impact on the sensitivity of the search are the theoretical uncertainties associated with the signal. These are de scribed in Section 9. The main experimental uncertain ties are associated with the JES, the jet energy resolution (JER), pile-up, Emiss, the b-tagging efficiency, the
T
W+jets transfer factor, and the integrated luminosity. The largest uncertainties on the backgrounds include WW normalisation and modelling, top normalisation, and Wγ(∗) normalisation. The 2-jet systematic uncertainties are dominated by the statistical uncertainties in the data and the MC simulation, and are therefore not discussed further.
Variations of the jet energy scale within the systematic uncertainties can cause events to migrate between the jet bins. The uncertainty on the JES varies from ±2% to ±9% as a function of jet pT and η for jets with pT > 25 GeV and |η| < 4. 5 [101]. The largest impact of this uncertainty on the total signal (background) yield amounts to 7% (4%) in the 0-jet (1-jet) bin. The uncertainty on the JER is estimated from in situ measurements and it impacts mostly the 1-jet channel, where its effect on the total signal and background yields is 4% and 2%, respectively. An additional contribution to the JES uncertainty arises from pile-up, and is estimated to vary between ±1% and ±5% for multiple pp collisions in the same bunch crossing and up to ±10% for neighbouring bunch crossings. This uncertainty affects mainly the 1-jet channel, where its impact on the signal and background yields is 4% and 2%, respectively. JES and lepton momentum scale uncertainties are propagated to the Emiss measurement. Additional contri-
T
butions to the Emiss uncertainties arise from jets with
T
pT < 20 GeV and from low-energy calorimeter deposits not associated with reconstructed physics objects [113]. The impact of the Emiss uncertainty on the total signal
T
and background yields is ∼3%. The efficiency of the b-tagging algorithm is calibrated using samples contain ing muons reconstructed in the vicinity of jets [114]. The uncertainty on the b-jet tagging efficiency varies between ±5% and ±18% as a function of the jet pT, and its impact on the total background yield is 10% for the 1-jet channel. The uncertainty in the W+jets transfer factor is dominated by differences in jet properties between dijet and W+jets events as observed in MC simulations. The total uncertainty on this background is approximately ±40%, resulting in an uncertainty on the total background yield of 5%. The uncertainty on the integrated luminosity is ±3.6%.
A fit to the distribution of mT is performed in order to obtain the signal yield for each mass hypoth esis (see Section 6.4). Most theoretical and experimental uncertainties do not produce statistically significant changes to the mT distribution. The uncertainties that do produce significant changes of the distribution of mT have no appreciable effect on the final results, with the exception of those associated with the WW background. In this case, an uncertainty is included to take into account differences in the distribution of mT and normalisation observed between the MCFM [115], MC@NLO+HERWIG and POWHEG+PYTHIA generators. The potential impact of interference between resonant (Higgs-mediated) and non-resonant gg→ WW diagrams [116] for mT > mH was investigated and found to be negligible. The effect of the WW normalisation, modelling, and shape systematics on the total background yield is 9% for the 0-jet channel and 19% for the 1-jet channel. The uncertainty on the shape of the total background is dominated by the uncertainties on the normalisations of the individual backgrounds. The main uncertainties on the top background in the 0-jet analysis include those associated with interference effects between tt¯and single top, initial state an final state radiation, b-tagging, and JER. The impact on the total background yield in the 0-jet bin is 3%. For the 1-jet analysis, the impact of the top background on the total yield is 14%. Theoretical uncertainties on the Wγ background normalisation are evaluated for each jet bin using the procedure described in Ref. [117]. They are ±11% for the 0-jet bin and ±50% for the 1-jet bin. For Wγ∗ with mee < 7 GeV, a k-factor of 1.3 ±0.3 is applied to the MadGraph LO prediction based on the comparison with the MCFM NLO calculation. The k-factor for Wγ∗/WZ(∗) with mee > 7 GeV is 1.5 ± 0.5. These uncertainties affect mostly the 1-jet channel, where their impact on the total background yield is approximately 4%.
Table 5: The expected numbers of signal (mH = 125 GeV) and background events after all selections, including a cut on the transverse mass of 0.75 mH < mT < mH for mH = 125 GeV. The observed numbers of events in data are also displayed. The eµ and µe channels are combined. The uncertainties shown are the combination of the statistical and all systematic uncertainties, taking into account the constraints from control samples. For the 2-jet analysis, backgrounds with fewer than 0.01 expected events are marked with ‘-’.
0-jet | 1-jet | 2-jet | |
---|---|---|---|
Signal | 20 ± 4 | 5 ± 2 | 0.34 ± 0.07 |
WW | 101 ± 13 | 12 ± 5 | 0.10 ± 0.14 |
WZ(∗)/ZZ/Wγ(∗) | 12 ± 3 | 1.9 ± 1.1 | 0.10 ± 0.10 |
t¯t | 8 ± 2 | 6 ± 2 | 0.15 ± 0.10 |
tW/tb/tqb | 3.4 ± 1.5 | 3.7 ± 1.6 | - |
Z/γ∗ + jets | 1.9 ± 1.3 | 0.10 ± 0.10 | - |
W + jets | 15 ± 7 | 2 ± 1 | - |
Total Background | 142 ± 16 | 26 ± 6 | 0.35 ± 0.18 |
Observed | 185 | 38 | 0 |
Table 5 shows the numbers of events expected from a SM Higgs boson with mH = 125 GeV and from the backgrounds, as well as the numbers of candidates observed in data, after application of all selection criteria plus an additional cut on mT of0.75 mH < mT < mH . The uncertainties shown in Table 5 include the system atic uncertainties discussed in Section 6.3, constrained by the use of the control regions discussed in Sec tion 6.2. An excess of events relative to the background expectation is observed in the data.
Figure 6 shows the distribution of the transverse mass after all selection criteria in the 0-jet and 1-jet channels combined, and for both lepton channels together.
The statistical analysis of the data employs a binned likelihood function constructed as the product of Poisson probability terms for the eµ channel and the µe
Figure 6: Distribution of the transverse mass, mT, in the 0-jet and 1-jet analyses with both eµ and µe channels combined, for events satisfying all selection criteria. The expected signal for mH = 125 GeV is shown stacked on top of the background prediction. The W+jets background is estimated from data, and WW and top background MC predictions are normalised to the data using control regions. The hashed area indicates the total uncertainty on the background prediction.
channel. The mass-dependent cuts on mT described above are not used. Instead, the 0-jet (1-jet) signal regions are subdivided into five (three) mT bins. For the 2-jet signal region, only the results integrated over mT are used, due to the small number of events in the final sample. The statistical interpretation of the observed excess of events is presented in Section 9.
7. Statistical procedure
The statistical procedure used to interpret the data is described in Refs. [17, 118–121]. The parameter of in terest is the global signal strength factor µ, which acts as a scale factor on the total number of events predicted by the Standard Model for the Higgs boson signal. This factor is defined such that µ = 0 corresponds to the background-only hypothesis and µ = 1 corresponds to the SM Higgs boson signal in addition to the background. Hypothesized values of µ are tested with a statistic λ(µ ) based on the profile likelihood ratio [122]. This test statistic extracts the information on the signal strength from a full likelihood fit to the data. The likelihood function includes all the parameters that describe the systematic uncertainties and their correlations.
Exclusion limits are based on the CLs prescrip tion [123]; a value of µ is regarded as excluded at 95% CL when CLs is less than 5%. A SM Higgs boson with mass mH is considered excluded at 95% confidence level (CL) when µ = 1 is excluded at that mass. The significance of an excess in the data is first quantified with the local p0, the probability that the background can produce a fluctuation greater than or equal to the excess observed in data. The equivalent formulation in terms of number of standard deviations, Zl, is referred to as the local significance. The global significance of a local excess anywhere in the search region, correcting for the “look-elsewhere” effect, is estimated with the method described in Ref. [124].
The statistical tests are performed in steps of values of the hypothesized Higgs boson mass mH . The asymp totic approximation [122] upon which the results are based has been validated previously [17].
The combination of individual search sub-channels for a specific Higgs boson decay, and the full combination of all search channels, are based on the global signal strength factor µ and on the identification of the nuisance parameters that correspond to correlated sources of systematic uncertainty.
8. Correlated systematic uncertainties
The individual search channels that enter the combi nation are summarised in Table 6.
The main uncorrelated systematic uncertainties are described in Sections 4–6 for the H→ ZZ(∗)→ 4e, H→ γγ and H→ WW(∗)→ eνeν channels and in Ref. [17] for the other channels. They include the background normalisations or background model parameters from control regions or sidebands, the Monte Carlo simulation statistical uncertainties and the theoretical uncertainties affecting the background processes.
The main sources of correlated systematic uncertainties are the following.
T
cluded, in addition to the JES uncertainty, which is due to low energy jet activity not associated with reconstructed physics objects.
6. Theory uncertainties: Correlated theoretical uncertainties affect mostly the signal predictions. The QCD scale uncertainties for mH=125 GeV amount to
+7%
−8% for the ggF process, ±1% for the VBF and WH/ZH processes, and for the t¯
+4% tH process [52, 53]; the
−9%
small dependence of these uncertainties on mH is taken into account. The uncertainties on the predicted branching ratios amount to ±5%. The uncertainties related to the parton distribution functions amount to ±8% for the predominantly gluon-initiated ggF and t¯
tH processes, and ±4% for the predominantly quark-initiated VBF and WH/ZH processes [78, 134–136]. The theoretical uncertainty associated with the exclusive Higgs boson production process with additional jets in the H→ γγ, H→ WW(∗)→ eνeν and H → τ+τ− channels is estimated using the prescription of Refs. [53, 117, 118], with the noticeable difference that an explicit calculation of the gluon-fusion process at NLO using MCFM [137] in the 2-jet category reduces the uncertainty on this non-negligible contribution to 25 %. An additional theoretical uncertainty on the signal normalisation of ±150%×(mH /TeV)3 (e.g. ±4% for mH = 300 GeV) accounts for effects related to off-shell Higgs boson pro duction and interference with other SM processes [53].
Sources of systematic uncertainty that affect both the 7 TeV and the 8 TeV data are taken as fully correlated. The uncertainties on background estimates based on control samples in the data are considered uncorrelated between the 7 TeV and 8 TeV data.
Table 6: Summary of the individual channels entering the combination. The transition points between separately optimised mH regions are indicated where applicable. In channels sensitive to associated production of the Higgs boson, V indicates a W or Z boson. The symbols ⊗ and ⊕ represent direct products and sums over sets of selection requirements, respectively.
Higgs Boson Decay | Subsequent Decay | Sub-Channels | mH Range[GeV] | L dt [fb−1] | Ref. |
---|---|---|---|---|---|
2011 √ s =7 TeV |
H → ZZ(∗) | 4e eeν¯ν eeq ¯q | {4e, 2e2µ, 2µ2e, 4µ}{ee, µµ} ⊗ {low, high pile-up}{b-tagged, untagged} | 110–600 200–280–600 200–300–600 | 4.8 4.7 4.7 | [87] [125] [126] |
H → γγ | – | 10 categories {pTt ⊗ ηγ ⊗ conversion} ⊕ {2-jet} | 110–150 | 4.8 | [127] |
H → WW(∗) | eνeν eνqq ' | {ee, eµ/µe, µµ} ⊗ {0-jet, 1-jet, 2-jet} ⊗ {low, high pile-up}{e, µ} ⊗ {0-jet, 1-jet, 2-jet} | 110–200–300–600 300–600 | 4.7 4.7 | [106] [128] |
H → ττ | τlepτlep τlepτhad τhadτhad | {eµ} ⊗ {0-jet} ⊕ {ee} ⊗ {1-jet, 2-jet, VH}{e, µ} ⊗ {0-jet} ⊗ {Emiss T < 20 GeV, Emiss T ≥ 20 GeV}⊕ {e, µ} ⊗ {1-jet} ⊕ {e} ⊗ {2-jet}{1-jet} | 110–150 110–150 110–150 | 4.7 4.7 4.7 | [129] |
VH → Vbb | Z → νν W → eν Z → ee | Emiss T ∈ {120 − 160, 160 − 200, ≥ 200 GeV}pW T ∈ {< 50, 50 − 100, 100 − 200, ≥ 200 GeV}pZ T ∈ {< 50, 50 − 100, 100 − 200, ≥ 200 GeV}√ | 110–130 110–130 110–130 | 4.6 4.7 4.7 | [130] |
2012 s =8 TeV
H → ZZ(∗) | 4e | {4e, 2e2µ, 2µ2e, 4µ} | 110–600 | 5.8 | [87] |
---|---|---|---|---|---|
H → γγ | – | 10 categories {pTt ⊗ ηγ ⊗ conversion} ⊕ {2-jet} | 110–150 | 5.9 | [127] |
H → WW(∗) | eνµν | {eµ, µe} ⊗ {0-jet, 1-jet, 2-jet} | 110–200 | 5.8 | [131] |
Table 7: Characterisation of the excess in the H → ZZ(∗) → 4e, H→ γγ and H→ WW(∗)→ eνeν channels and the combination of all channels listed in Table 6. The mass value mmax for which the local significance is maximum, the maximum observed local significance Zl and the expected local significance E(Zl) in the presence of a SM Higgs boson signal at mmax are given. The best fit value of the signal strength parameter ˆµ at mH = 126 GeV is shown with the total uncertainty. The expected and observed mass ranges excluded at 95% CL (99% CL, indicated by a *) are
√√
also given, for the combined s = 7 TeV and s = 8 TeV data.
Search channel | Dataset | mmax [GeV] | Zl [σ] | E(Zl) [σ] | ˆµ(mH = 126 GeV) | Expected exclusion [GeV] | Observed exclusion [GeV] |
---|---|---|---|---|---|---|---|
7 TeV | 125.0 | 2.5 | 1.6 | 1.7 ± 1.1 | |||
H → ZZ(∗) → 4e | 8 TeV | 125.5 | 2.6 | 2.1 | 1.3 ± 0.8 | ||
7 & 8 TeV | 125.0 | 3.6 | 2.7 | 1.4 ± 0.6 | 124–164, 176–500 | 131–162, 170–460 | |
7 TeV | 126.0 | 3.4 | 1.6 | 2.2 ± 0.7 | |||
H→ γγ | 8 TeV | 127.0 | 3.2 | 1.9 | 1.5 ± 0.6 | ||
7 & 8 TeV | 126.5 | 4.5 | 2.5 | 1.8 ± 0.5 | 110–140 | 112–123, 132–143 | |
7 TeV | 135.0 | 1.1 | 3.4 | 0.5 ± 0.6 | |||
H→ WW(∗)→ eνeν | 8 TeV | 120.0 | 3.3 | 1.0 | 1.9 ± 0.7 | ||
7 & 8 TeV | 125.0 | 2.8 | 2.3 | 1.3 ± 0.5 | 124–233 | 137–261 | |
7 TeV | 126.5 | 3.6 | 3.2 | 1.2 ± 0.4 | |||
Combined | 8 TeV 7 & 8 TeV | 126.5 126.5 | 4.9 6.0 | 3.8 4.9 | 1.5 ± 0.4 1.4 ± 0.3 | 110–582 113–532 (*) | 111–122, 131–559 113–114, 117–121, 132–527 (*) |
9. Results
The addition of the 8 TeV data for the H→ ZZ(∗)→ 4e, H→ γγ and H→ WW(∗)→ eνµν channels, as well as the improvements to the analyses of the 7 TeV data in the first two of these channels, bring a significant gain in sensitivity in the low-mass region with respect to the previous combined search [17].
The combined 95% CL exclusion limits on the production of the SM Higgs boson, expressed in terms of the signal strength parameter µ , are shown in Fig. 7(a) as a function of mH . The expected 95% CL exclusion region covers the mH range from 110 GeV to 582 GeV. The observed 95% CL exclusion regions are 111–122 GeV and 131–559 GeV. Three mass regions are excluded at 99% CL, 113–114, 117–121 and 132– 527 GeV, while the expected exclusion range at 99% CL is 113–532 GeV.
An excess of events is observed near mH=126 GeV in the H→ ZZ(∗)→ 4e and H→ γγ channels, both of which
Figure 7: Combined search results: (a) The observed (solid) 95% CL limits on the signal strength as a function of mH and the expectation (dashed) under the background-only hypothesis. The dark and light shaded bands show the ±1σ and ±2σ uncertainties on the background-only expectation. (b) The observed (solid) local p0 asa function of mH and the expectation (dashed) for a SM Higgs boson signal hypothesis (µ = 1) at the given mass. (c) The best-fit signal strength ˆµ as a function of mH . The band indicates the approximate 68% CL interval around the fitted value.
provide fully reconstructed candidates with high reso lution in invariant mass, as shown in Figures 8(a) and 8(b). These excesses are confirmed by the highly sen sitive but low-resolution H→ WW(∗)→ eνeν channel, as shown in Fig. 8(c).
The observed local p0 values from the combination of channels, using the asymptotic approximation, are shown as a function of mH in Fig. 7(b) for the full mass range and in Fig. 9 for the low mass range.
The largest local significance for the combination of the 7 and 8 TeV data is found for a SM Higgs boson mass hypothesis of mH =126.5 GeV, where it reaches
6.0 σ, with an expected value in the presence of a SM Higgs boson signal at that mass of 4.9 σ (see also Ta ble 7). For the 2012 data alone, the maximum local significance for the H→ ZZ(∗)→ 4e, H→ γγ and
Figure 8: The observed local p0 as a function of the hypothesized Higgs boson mass for the (a) H→ ZZ(∗)→ 4e, (b) H→ γγ and (c) H→ WW(∗)→ eνeν channels. The dashed curves show the expected local p0 under the hypothesis of a SM Higgs boson signal at that mass.
√
Results are shown separately for the s = 7 TeV data (dark, blue), the
√
s = 8 TeV data (light, red), and their combination (black).
H→ WW(∗)→ eνµν channels combined is 4.9 σ, and occurs at mH = 126.5 GeV (3.8 σ expected).
The significance of the excess is mildly sensitive to uncertainties in the energy resolutions and energy scale systematic uncertainties for photons and electrons; the effect of the muon energy scale systematic uncertainties is negligible. The presence of these uncertainties, evaluated as described in Ref. [138], reduces the local significance to 5.9 σ.
The global significance of a local 5.9 σ excess anywhere in the mass range 110–600 GeV is estimated to be approximately 5.1 σ, increasing to 5.3 σ in the range 110–150 GeV, which is approximately the mass range not excluded at the 99% CL by the LHC combined SM Higgs boson search [139] and the indirect constraints from the global fit to precision electroweak measure ments [12].
18
Figure 9: The observed (solid) local p0 as a function of mH in the low mass range. The dashed curve shows the expected local p0 under the hypothesis of a SM Higgs boson signal at that mass with its ±1σ band. The horizontal dashed lines indicate the p-values corresponding to significances of 1 to 6 σ.
The mass of the observed new particle is estimated using the profile likelihood ratio λ(mH ) for H→ ZZ(∗)→ 4e and H→ γγ, the two channels with the highest mass resolution. The signal strength is allowed to vary independently in the two channels, although the result is essentially unchanged when restricted to the SM hypothesis µ = 1. The leading sources of systematic uncertainty come from the electron and photon energy scales and resolutions. The resulting estimate for the mass of the observed particle is
The best-fit signal strength ˆµ is shown in Fig. 7(c) as a function of mH . The observed excess corresponds to µˆ= 1.4 ± 0.3 for mH = 126 GeV, which is consistent with the SM Higgs boson hypothesis µ = 1. A summary of the individual and combined best-fit values of the strength parameter for a SM Higgs boson mass hy pothesis of 126 GeV is shown in Fig. 10, while more information about the three main channels is provided in Table 7.
In order to test which values of the strength and mass of a signal hypothesis are simultaneously consistent with the data, the profile likelihood ratio λ(µ, mH ) is used. In the presence of a strong signal, it will produce closed contours around the best-fit point (ˆµ, mˆH ), while in the absence of a signal the contours will be upper limits on µ for all values of mH .
Asymptotically, the test statistic −2 ln λ(µ, mH ) is distributed as a χ2 distribution with two degrees of freedom. The resulting 68% and 95% CL contours for the H→ γγ and H→ WW(∗)→ eνeν channels are shown in
Figure 10: Measurements of the signal strength parameter µ for mH =126 GeV for the individual channels and their combination.
Fig. 11, where the asymptotic approximations have been validated with ensembles of pseudo-experiments. Similar contours for the H→ ZZ(∗)→ 4e channel are also shown in Fig. 11, although they are only approximate confidence intervals due to the smaller number of candidates in this channel. These contours in the (µ, mH) plane take into account uncertainties in the energy scale and resolution.
The probability for a single Higgs boson-like particle to produce resonant mass peaks in the H→ ZZ(∗)→ 4e and H→ γγ channels separated by more than the observed mass difference, allowing the signal strengths to vary independently, is about 20%.
The contributions from the different production modes in the H→ γγ channel have been studied in order to assess any tension between the data and the ratios of the production cross sections predicted in the Standard Model. A new signal strength parameter µi is introduced for each production mode, defined by µi = σi/σi,SM. In order to determine the values of (µi,µ j) that are simultaneously consistent with the data, the profile likelihood ratio λ(µi,µ j) is used with the measured mass treated as a nuisance parameter.
Since there are four Higgs boson production modes at the LHC, two-dimensional contours require either some µi to be fixed, or multiple µi to be related in some way. Here, µggF and µt¯
tH have been grouped together as they scale with the t¯
tH coupling in the SM, and are denoted by the common parameter µggF+t¯
tH . Similarly, µVBF and µVH have been grouped together as they scale with the WWH/ZZH coupling in the SM, and are denoted by the common parameter µVBF+VH . Since the distribution of signal events among the 10 categories of the H→ γγ search is sensitive to these factors, constraints in the
Figure 11: Confidence intervals in the (µ, mH ) plane for the H→ ZZ(∗)→ 4e, H→ γγ, and H→ WW(∗)→ eνeν channels, including all systematic uncertainties. The markers indicate the maximum likelihood estimates (ˆµ, mˆH ) in the corresponding channels (the maximum likelihood estimates for H→ ZZ(∗)→ 4e and H→ WW(∗)→ eνeν coincide).
plane of µggF+t¯
tH × B/BSM and µVBF+VH × B/BSM, where B is the branching ratio for H→ γγ, can be obtained (Fig. 12). Theoretical uncertainties are included so that the consistency with the SM expectation can be quantified. The data are compatible with the SM expectation at the 1.5 σ level.
10. Conclusion
Searches for the Standard Model Higgs boson have been performed in the H→ ZZ(∗)→ 4e, H→ γγ and H→ WW(∗)→ eνµν channels with the ATLAS experiment at the LHC using 5.8–5.9 fb−1 of pp collision data recorded during April to June 2012 at a centre-of-mass energy of 8 TeV. These results are combined with ear lier results [17], which are based on an integrated lu minosity of 4.6–4.8 fb−1 recorded in 2011 at a centreof-mass energy of 7 TeV, except for the H→ ZZ(∗)→ 4e and H→ γγ channels, which have been updated with the improved analyses presented here.
The Standard Model Higgs boson is excluded at 95% CL in the mass range 111–559 GeV, except for the narrow region 122–131 GeV. In this region, an excess of events with significance 5.9 σ, corresponding to p0 = 1.7 × 10−9, is observed. The excess is driven by the two channels with the highest mass resolution, H→ ZZ(∗)→ 4e and H→ γγ, and the equally sensitive but low-resolution H→ WW(∗)→ eνeν channel. Taking into account the entire mass range of the search, 110– 600 GeV, the global significance of the excess is 5.1 σ, which corresponds to p0 = 1.7 × 10−7.
tH (µVBF+VH ) is a common scale factor for the ggF and t¯
tH (VBF and VH) production cross sections. The best fit to the data (+) and 68% (full) and 95% (dashed) CL contours are also indicated, as well as the SM expectation (×).
These results provide conclusive evidence for the discovery of a new particle with mass
126.0 ± 0.4 (stat) ± 0.4 (sys) GeV. The signal strength parameter µ has the value 1.4 ± 0.3 at the fitted mass, which is consistent with the SM Higgs boson hypothesis µ = 1. The decays to pairs of vector bosons whose net electric charge is zero identify the new particle as a neutral boson. The observation in the diphoton channel disfavours the spin-1 hypothe sis [140, 141]. Although these results are compatible with the hypothesis that the new particle is the Standard Model Higgs boson, more data are needed to assess its nature in detail.
Acknowledgements
The results reported in this paper would not have been possible without the outstanding performance of the LHC. We warmly thank CERN and the entire LHC exploitation team, including the operation, technical and infrastructure groups, and all the people who have contributed to the conception, design and construction of this superb accelerator We thank also the support staff at our institutions without whose excellent contributions ATLAS could not have been successfully constructed or operated so efficiently.
We acknowledge the support of ANPCyT, Argentina; YerPhI, Armenia; ARC, Australia; BMWF, Austria; ANAS, Azerbaijan; SSTC, Belarus; CNPq and FAPESP, Brazil; NSERC, NRC and CFI, Canada; CERN; CONICYT, Chile; CAS, MOST and NSFC, China; COLCIENCIAS, Colombia; MSMT CR, MPO CR and VSC CR, Czech Republic; DNRF, DNSRC and Lundbeck Foundation, Denmark; EPLANET and ERC, European Union; IN2P3-CNRS, CEA-DSM/IRFU, France; GNAS, Georgia; BMBF, DFG, HGF, MPG and AvH Foundation, Germany; GSRT, Greece; ISF, MINERVA, GIF, DIP and Benoziyo Center, Israel; INFN, Italy; MEXT and JSPS, Japan; CNRST, Morocco; FOM and NWO, Netherlands; RCN, Norway; MNiSW, Poland; GRICES and FCT, Portugal; MERYS (MECTS), Romania; MES of Russia and ROSATOM, Russian Federation; JINR; MSTD, Serbia; MSSR, Slovakia; ARRS and MVZT, Slovenia; DST/NRF, South Africa; MICINN, Spain; SRC and Wallenberg Foundation, Sweden; SER, SNSF and Cantons of Bern and Geneva, Switzerland; NSC, Taiwan; TAEK, Turkey; STFC, the Royal Society and Leverhulme Trust, United Kingdom; DOE and NSF, United States of America.
The crucial computing support from all WLCG partners is acknowledged gratefully, in particular from CERN and the ATLAS Tier-1 facilities at TRIUMF (Canada), NDGF (Denmark, Norway, Sweden), CCIN2P3 (France), KIT/GridKA (Germany), INFN-CNAF (Italy), NL-T1 (Netherlands), PIC (Spain), ASGC (Taiwan), RAL (UK) and BNL (USA) and in the Tier-2 facilities worldwide.
References
[1] S. L. Glashow, Partial-symmetries of weak interactions, Nucl. Phys. 22 no. 4, (1961) 579.
[2] S. Weinberg, A Model of Leptons , Phys. Rev. Lett. 19 (1967) 1264.
[3] A. Salam, “Weak and electromagnetic interactions, in Elementary particle physics: relativistic groups and analyticity, N. Svartholm, ed., p. 367. Almqvist & Wiskell, 1968. Proceedings of the eighth Nobel symposium.
[4] G. ’t Hooft and M. Veltman, Regularization and Renormalization of Gauge Fields , Nucl. Phys. B44 (1972)
[5] F. Englert and R. Brout, Broken symmetry and the mass of gauge vector mesons , Phys. Rev. Lett. 13 (1964) 321.
[6] P. W. Higgs, Broken symmetries, massless particles and gauge fields , Phys. Lett. 12 (1964) 132.
[7] P. W. Higgs, Broken symmetries and the masses of gauge bosons , Phys. Rev. Lett. 13 (1964) 508.
[8] G. S. Guralnik, C. R. Hagen, and T. W. B. Kibble, Global conservation laws and massless particles , Phys. Rev. Lett. 13 (1964) 585.
[9] P. W. Higgs, Spontaneous symmetry breakdown without massless bosons , Phys. Rev. 145 (1966) 1156.
[10] T. W. B. Kibble, Symmetry breaking in non-Abelian gauge theories , Phys. Rev. 155 (1967) 1554.
[11] L. Evans and P. Bryant (Eds.), LHC Machine , JINST 3 (2008) S08001.
[12] ALEPH, CDF, DØ, DELPHI, L3, OPAL, SLD Collaborations, the LEP Electroweak Working Group, the Tevatron Electroweak Working Group, and the SLD electroweak and
heavy flavour groups, Precision Electroweak Measurements and Constraints on the Standard Model, CERN-PH-EP-2010-095 (2010), arXiv:1012.2367 [hep-ex].
[13] ALEPH, DELPHI, L3 and OPAL Collaborations, The LEP Working Group for Higgs boson searches, , Search for the standard model Higgs boson at LEP , Phys. Lett. B 565 (2003)
[14] CDF Collaboration, T. Aaltonen et al., Combined search for the standard model Higgs boson decaying to a b¯
b pair using
the full CDF data set, submitted to Phys. Rev. Lett. (2012),
[15] DØ Collaboration, V. M. Abazov et al., Combined search for the standard model Higgs boson decaying to b¯
b using the DØ
Run II data set, submitted to Phys. Rev. Lett. (2012),
[16] CDF Collaboration, DØ Collaboration, Evidence for a particle produced in association with weak bosons and decaying to a bottom-antibottom quark pair in Higgs boson searches at the Tevatron, submitted to Phys. Rev. Lett. (2012), arXiv:1207.6436 [hep-ex].
[17] ATLAS Collaboration, Combined search for the Standard
√
Model Higgs boson in pp collisions at s=7 TeV with the ATLAS detector, accepted by Phys. Rev. D (2012), arXiv:1207.0319 [hep-ex].
[18] CMS Collaboration, Combined results of searches for the
√
standard model Higgs boson in pp collisions at s = 7 TeV, Phys. Lett. B 710 (2012) 26.
[19] ATLAS Collaboration, ATLAS: letter of intent for a general-purpose pp experiment at the large hadron collider at CERN , CERN-LHCC-92-004 (1992).
[20] ATLAS Collaboration, The ATLAS Collaboration, ATLAS Technical Proposal for a General-Purpose pp Experiment at the Large Hadron Collider at CERN , CERN-LHCC-94-43 (1994).
[21] ATLAS Collaboration, The ATLAS Experiment at the CERN Large Hadron Collider , JINST 3 (2008) S08003.
[22] H. Georgi, S. Glashow, M. Machacek, and D. V. Nanopoulos, Higgs bosons from two gluon annihilation in proton proton collisions , Phys. Rev. Lett. 40 (1978) 692.
[23] A. Djouadi, M. Spira, and P. M. Zerwas, Production of Higgs bosons in proton colliders: QCD corrections , Phys. Lett. B 264 (1991) 440.
[24] S. Dawson, Radiative corrections to Higgs boson production, Nucl. Phys. B359 (1991) 283.
[25] M. Spira, A. Djouadi, D. Graudenz, and P. M. Zerwas, Higgs boson production at the LHC , Nucl. Phys. B453 (1995) 17.
[26] R. Harlander and W. B. Kilgore, Next-to-next-to-leading order Higgs production at hadron colliders , Phys. Rev. Lett. 88 (2002) 201801.
[27] C. Anastasiou and K. Melnikov, Higgs boson production at hadron colliders in NNLO QCD , Nucl. Phys. B646 (2002)
[28] V. Ravindran, J. Smith, and W. L. van Neerven, NNLO corrections to the total cross section for Higgs boson production in hadron hadron collisions , Nucl. Phys. B665 (2003) 325.
[29] U. Aglietti, R. Bonciani, G. Degrassi, and A. Vicini, Two-loop light fermion contribution to Higgs production and decays, Phys. Lett. B 595 (2004) 432.
[30] S. Actis, G. Passarino, C. Sturm, and S. Uccirati, NLO electroweak corrections to Higgs boson production at hadron colliders , Phys. Lett. B 670 (2008) 12.
[31] S. Catani, D. de Florian, M. Grazzini, and P. Nason,
Soft-gluon re-summation for Higgs boson production at
hadron colliders , JHEP 0307 (2003) 028.
[32] C. Anastasiou, R. Boughezal, and F. Petriello, Mixed QCD-electroweak corrections to Higgs boson production in gluon fusion , JHEP 0904 (2009) 003.
[33] D. de Florian and M. Grazzini, Higgs production at the LHC:
√
updated cross sections at s = 8 TeV, arXiv:1206.4133 [hep-ph].
[34] C. Anastasiou, S. Buehler, F. Herzog, and A. Lazopoulos, Inclusive Higgs boson cross-section for the LHC at 8 TeV, JHEP 1204 (2012) 004.
[35] J. Baglio and A. Djouadi, Higgs production at the lHC , JHEP 1103 (2011) 055.
[36] D. de Florian, G. Ferrera, M. Grazzini, and D. Tommasini, Transverse-momentum resummation: Higgs boson production at the Tevatron and the LHC , JHEP 1111 (2011) 064.
[37] P. S. E. Bagnaschi, G. Degrassi and A. Vicini, Higgs production via gluon fusion in the POWHEG approach in the SM and in the MSSM , JHEP 1202 (2012) 88.
[38] R. Cahn and S. Dawson, Production of very massive Higgs bosons , Phys. Lett. B 136 (1984) 196, Erratum–ibid. B 138 (1984) 464.
[39] M. Ciccolini, A. Denner, and S. Dittmaier, Strong and electroweak corrections to the production of Higgs+2jets via weak interactions at the LHC , Phys. Rev. Lett. 99 (2007) 161803.
[40] M. Ciccolini, A. Denner, and S. Dittmaier, Electroweak and QCD corrections to Higgs production via vector-boson fusion at the LHC , Phys. Rev. D77 (2008) 013002.
[41] K. Arnold, M. Bahr, G. Bozzi, F. Campanario, C. Englert, et al., VBFNLO: A parton level Monte Carlo for processes with electroweak bosons , Comput. Phys. Commun. 180 (2009) 1661.
[42] P. Bolzoni, F. Maltoni, S.-O. Moch, and M. Zaro, Higgs production via vector-boson fusion at NNLO in QCD , Phys. Rev. Lett. 105 (2010) 011801.
[43] S. Glashow, D. V. Nanopoulos, and A. Yildiz, Associated production of Higgs bosons and Z particles , Phys. Rev. D18 (1978) 1724.
[44] T. Han and S. Willenbrock, QCD correction to the pp → WH and ZH total cross sections , Phys. Lett. B 273 (1991) 167.
[45] O. Brein, A. Djouadi, and R. Harlander, NNLO QCD corrections to the Higgs-strahlung processes at hadron colliders , Phys. Lett. B 579 (2004) 149.
[46] M. L. Ciccolini, S. Dittmaier, and M. Kramer, Electroweak radiative corrections to associated WH and ZH production at hadron colliders , Phys. Rev. D68 (2003) 073003.
[47] Z. Kunszt, Associated production of heavy Higgs boson with top quarks , Nucl. Phys. B247 (1984) 339.
[48] W. Beenakker et al., Higgs Radiation Off Top Quarks at the Tevatron and the LHC , Phys. Rev. Lett. 87 (2001) 201805.
[49] W. Beenakker et al., NLO QCD corrections to t¯tH production in hadron collisions , Nucl. Phys. B653 (2003) 151.
[50] S. Dawson, L. H. Orr, L. Reina, and D. Wackeroth, Next-to-leading order QCD corrections to pp → tth at the CERN Large Hadron Collider , Phys. Rev. D67 (2003) 071503.
[51] S. Dawson, C. Jackson, L. H. Orr, L. Reina, and
D. Wackeroth, Associated Higgs production with top quarks at the Large Hadron Collider: NLO QCD corrections , Phys. Rev. D68 (2003) 034022.
[52] LHC Higgs Cross Section Working Group, S. Dittmaier,
C. Mariotti, G. Passarino, and R. Tanaka (Eds.), Handbook of LHC Higgs cross sections: 1. Inclusive observables, CERN-2011-002 (2011), arXiv:1101.0593 [hep-ph].
[53] LHC Higgs Cross Section Working Group, S. Dittmaier, C. Mariotti, G. Passarino, and R. Tanaka (Eds.), Handbook of LHC Higgs Cross Sections: 2. Differential Distributions, CERN-2012-002 (2012), arXiv:1201.3084 [hep-ph].
[54] A. Djouadi, J. Kalinowski, and M. Spira, HDECAY: A program for Higgs boson decays in the standard model and its supersymmetric extension , Comput. Phys. Commun. 108 (1998) 56.
[55] A. Bredenstein, A. Denner, S. Dittmaier, and M. M. Weber,
Precise predictions for the Higgs-boson decay H → WW/ZZ
→ 4 leptons , Phys. Rev. D74 (2006) 013004.
[56] A. Bredenstein, A. Denner, S. Dittmaier, and M. M. Weber, Radiative corrections to the semileptonic and hadronic Higgs-boson decays H → WW/ZZ → 4 fermions , JHEP 0702 (2007) 080.
[57] S. Alioli, P. Nason, C. Oleari, and E. Re, NLO Higgs boson production via gluon fusion matched with shower in POWHEG , JHEP 0904 (2009) 002.
[58] P. Nason and C. Oleari, NLO Higgs boson production via vector-boson fusion matched with shower in POWHEG , JHEP 1002 (2010) 037.
[59] M. L. Mangano et al., ALPGEN, a generator for hard multiparton processes in hadronic collisions , JHEP 0307 (2003) 001.
[60] S. Frixione and B. R. Webber, Matching NLO QCD computations and parton shower simulations , JHEP 0206 (2002) 029.
[61] B. P. Kersevan and E. Richter-Was, The Monte Carlo event generator AcerMC version 2.0 with interfaces to PYTHIA 6.2 and HERWIG 6.5, arXiv:hep-ph/0405247.
[62] T. Binoth, M. Ciccolini, N. Kauer and M. Kr¨amer, Gluon-induced W-boson pair production at the LHC , JHEP 0612 (2006) 046.
[63] T. Melia, P. Nason, R. Rontsch, and G. Zanderighi, W+W− , WZ and ZZ production in the POWHEG BOX , JHEP 1111 (2011) 078.
[64] T. Binoth, N. Kauer and P. Mertsch, Gluon-induced QCD Corrections to pp → ZZ → eee e , arXiv:0807.0024
¯ ' ¯' [hep-ph].
[65] R. C. Gray, C. Kilic, M. Park, S. Somalwar, and S. Thomas, Backgrounds to Higgs Boson Searches from Wγ ∗→ lνl(l) Asymmetric Internal Conversion, arXiv:1110.1368 [hep-ph].
[66] ATLAS Collaborations, New ATLAS event generator tunes to 2010 data, ATL-PHYS-PUB-2011-008 (2011). http://cdsweb.cern.ch/record/1345343.
[67] ATLAS Collaborations, ATLAS tunes of PYTHIA 6 and Pythia 8 for MC11, ATL-PHYS-PUB-2011-009 (2011). http://cdsweb.cern.ch/record/1363300.
[68] ATLAS Collaborations, Further ATLAS tunes of PYTHIA6 and Pythia 8, ATL-PHYS-PUB-2011-014 (2011). http://cdsweb.cern.ch/record/1400677.
[69] T. Sjostrand, S. Mrenna, and P. Z. Skands, PYTHIA 6.4 physics and manual , JHEP 0605 (2006) 026.
[70] J. Alwall et al., MadGraph/MadEvent v4: The New Web Generation , JHEP 0709 (2007) 028.
[71] J. Alwall, M. Herquet, F. Maltoni, O. Mattelaer, and T. Stelzer, MadGraph 5 : Going Beyond , JHEP 1106 (2011) 128.
[72] T. Sjostrand, S. Mrenna, and P. Z. Skands, A Brief Introduction to PYTHIA 8.1 , Comput. Phys. Commun. 178 (2008) 852.
[73] G. Corcella, I. Knowles, G. Marchesini, S. Moretti,
K. Odagiri, et al., HERWIG 6: An Event generator for hadron emission reactions with interfering gluons (including supersymmetric processes) , JHEP 0101 (2001) 010.
[74] T. Gleisberg et al., Event generation with SHERPA 1.1 , JHEP
[75] J. M. Butterworth, J. R. Forshaw, and M. H. Seymour, Multiparton interactions in photoproduction at HERA , Z. Phys. C72 (1996) 637.
[76] S. Jadach, Z. Was, R. Decker, and J. H. Kuhn, The tau decay library TAUOLA: Version 2.4 , Comput. Phys. Commun. 76 (1993) 361.
[77] P. Golonka and Z. Was, PHOTOS Monte Carlo: A Precision tool for QED corrections in Z and W decays , Eur.Phys.J. C45 (2006) 97.
[78] H.-L. Lai et al., New parton distributions for collider physics, Phys. Rev. D82 (2010) 074024.
[79] P. M. Nadolsky et al., Implications of CTEQ global analysis for collider observables , Phys. Rev. D78 (2008) 013004.
[80] A. Sherstnev and R. S. Thorne, Parton distributions for the LHC , Eur. Phys. J C55 (2009) 553, arXiv:0711.2473 [hep-ph].
[81] ATLAS Collaboration, The ATLAS simulation infrastructure, Eur. Phys. J. C70 (2010) 823.
[82] S. Agostinelli et al., Geant4, a simulation toolkit , Nucl. Instrum. Meth. A506 (2003) 250.
[83] ATLAS Collaboration, Search for the Standard Model Higgs boson in the decay channel H → ZZ(∗) → 4e with 4.8 fb−1 of
√
pp collision data at s = 7 TeV with ATLAS , Phys. Lett. B 710 (2012) 383.
[84] ATLAS Collaboration, Muon reconstruction efficiency in reprocessed 2010 LHC proton-proton collision data recorded with the ATLAS detector, ATLAS-CONF-2011-063 (2011). http://cdsweb.cern.ch/record/1345743.
[85] ATLAS Collaboration, Electron performance measurements with the ATLAS detector using the 2010 LHC proton-proton collision data , Eur. Phys. J. C72 (2012) 1909.
[86] ATLAS Collaboration, Improved electron reconstruction in ATLAS using the Gaussian Sum Filter-based model for bremsstrahlung,. http://cdsweb.cern.ch/record/1449796.
[87] ATLAS Collaboration, Observation of an excess of events in the search for the Standard Model Higgs boson in the H → ZZ(∗) → 4e channel with the ATLAS detector, ATLAS-CONF-2012-092 (2012). http://cdsweb.cern.ch/record/1460411.
[88] W. Lampl, S. Laplace, D. Lelas, P. Loch, H. Ma, S. Menke,
S. Rajagopalan, D. Rousseau, S. Snyder, and G. Unal,
Calorimeter Clustering Algorithms: Description and Performance, ATL-LARG-PUB-2008-002 (2008). http://cdsweb.cern.ch/record/1099735.
[89] M. Cacciari and G. P. Salam, Pileup subtraction using jet areas , Phys. Lett. B 659 (2008) 119.
[90] M. Cacciari, G. P. Salam, and G. Soyez, FastJet user manual, Eur. Phys. J. C72 (2012) 1896.
[91] ATLAS Collaboration, Measurements of the electron and muon inclusive cross-sections in proton-proton collisions at
√
s = 7 TeV with the ATLAS detector, Phys. Lett. B 707 (2012) 438, arXiv:1109.0525 [hep-ex].
[92] ATLAS Collaboration, Improved Luminosity Determination in
√
pp Collisions at s = 7 TeV using the ATLAS Detector at the
LHC, ATLAS-CONF-2012-080 (2012).
http://cdsweb.cern.ch/record/1460392.
[93] ATLAS Collaboration, Measurement of the total ZZ production cross section in the four-lepton channel using 5.8 fb−1 of ATLAS data at √ s = 8 TeV, ATLAS-CONF-2012-090 (2012). http://cdsweb.cern.ch/record/1460409.
[94] N. Kauer and G. Passarino, Inadequacy of zero-width approximation for a light Higgs boson signal, arXiv:1206.4803 [hep-ph].
[95] ATLAS Collaboration, Search for the Standard Model Higgs boson in the diphoton decay channel with 4.9 fb−1 of pp
√
collision data at s = 7 TeV with ATLAS , Phys. Rev. Lett. 108 (2012) 111803.
[96] ATLAS Collaboration, Performance of the ATLAS Trigger System in 2010, Eur. Phys. J. C72 (2012) 1849, arXiv:1110.1530 [hep-ex].
[97] ATLAS Collaboration, Measurement of the inclusive isolated
√
prompt photon cross section in pp collisions at s = 7 TeV with the ATLAS detector , Phys. Rev. D83 (2011) 052005.
[98] ATLAS Collaboration, Search for the Standard Model Higgs boson in the two photon decay channel with the ATLAS detector at the LHC , Phys. Lett. B 705 (2011) 452.
[99] OPAL Collaboration Collaboration, K. Ackerstaff et al.,
Search for anomalous production of dilepton events with missing transverse momentum in e+e− collisions at √ s = 161 GeV and 172 GeV , Eur.Phys.J. C4 (1998) 47, arXiv:hep-ex/9710010 [hep-ex].
[100] M. Vesterinen and T. R. Wyatt, A Novel Technique for Studying the Z Boson Transverse Momentum Distribution at Hadron Colliders , Nucl. Instrum. Meth. A602 (2009) 432.
[101] ATLAS Collaboration, Jet energy measurement with the
√
ATLAS detector in proton-proton collisions at s = 7 TeV, submitted to Eur. Phys. J. C (2011), arXiv:1112.6426 [hep-ex].
[102] M. Cacciari, G. P. Salam, and G. Soyez, The anti-kt jet clustering algorithm , JHEP 0804 (2008) 063.
[103] L. J. Dixon and M. S. Siu, Resonance continuum interference in the diphoton Higgs signal at the LHC , Phys. Rev. Lett. 90 (2003) 252001.
[104] J. Gaiser, Charmonium spectroscopy from radiative decays of the J/Ψ and Ψ ', Ph.D. Thesis no. SLAC-R-255, (1982).
[105] S. N. Bernstein, D´eor`emonstration du th´eme de Weierstrass fond´es, Comm. Soc. Math.
ee sur le calcul des probabilit´ Kharkov 13 (1912) 1.
[106] ATLAS Collaboration, Search for the Standard Model Higgs boson in the H → WW(∗) → eνeν decay mode with 4.7 fb−1 of
√
ATLAS data at s = 7 TeV, submitted to Phys. Lett. B (2012), arXiv:1206.0756 [hep-ex].
[107] M. Dittmar and H. Dreiner, How to find a Higgs boson with a mass between 155 and 180 GeV at the CERN LHC , Phys. Rev. D55 (2007) 167.
[108] ATLAS Collaboration, Commissioning of the ATLAS high-performance b-tagging algorithms in the 7 TeV collision data, ATLAS-CONF-2011-102 (2011). http://cdsweb.cern.ch/record/1369219.
[109] R.K. Ellis et al., Higgs decay to τ+τ−: A possible signature of intermediate mass Higgs bosons at the SSC , Nucl. Phys. B297 (1988) 221.
[110] A. J. Barr, B. Gripaios, and C. G. Lester, Measuring the Higgs boson mass in dileptonic W-boson decays at hadron colliders, JHEP 0907 (2009) 072.
[111] ATLAS Collaboration, Search for the Higgs boson in the √H → WW(∗) → eνeν decay channel in pp collisions at
s = 7 TeV with the ATLAS detector , Phys. Rev. Lett. 108 (2012) 111802.
[112] ATLAS Collaboration, Limits on the production of the
√
Standard Model Higgs Boson in pp collisions at s = 7 TeV with the ATLAS detector , Eur. Phys. J. C 71 (2011) 1728.
[113] ATLAS Collaboration, Performance of missing transverse momentum reconstruction in proton-proton collisions at 7 TeV with ATLAS , Eur. Phys. J. C72 (2012) 1844.
[114] ATLAS Collaboration, Measurement of the b-tag Efficiency in a Sample of Jets Containing Muons with 5 fb−1 of Data from the ATLAS Detector, ATLAS-CONF-2012-043 (2012).
http://cdsweb.cern.ch/record/1435197.
[115] J. M. Campbell, R. K. Ellis, and C. Williams, Vector boson pair production at the LHC, JHEP 1107 (2011) 018, arXiv:1105.0020 [hep-ph].
[116] J. Campbell, R. Ellis, and C. Williams, Gluon-gluon contributions to W+W− production and Higgs interference effects, JHEP 1110 (2011) 005, arXiv:1107.5569 [hep-ph].
[117] I. Stewart and F. Tackmann, Theory uncertainties for Higgs mass and other searches using jet bins , Phys. Rev. D85 (2012) 034011.
[118] ATLAS and CMS Collaborations, Procedure for the LHC Higgs boson search combination in Summer 2011, ATL-PHYS-PUB-2011-011, CERN-CMS-NOTE-2011-005 (2011). http://cdsweb.cern.ch/record/1375842.
[119] L. Moneta, K. Belasco, K. S. Cranmer, S. Kreiss, A. Lazzaro, et al., The RooStats Project, PoS ACAT2010 (2010) 057, arXiv:1009.1003 [physics.data-an].
[120] K. Cranmer, G. Lewis, L. Moneta, A. Shibata, and
W. Verkerke, HistFactory: A tool for creating statistical models for use with RooFit and RooStats, CERN-OPEN-2012-016 (2012). http://cdsweb.cern.ch/record/1456844.
[121] W. Verkerke and D. Kirkby, The RooFit toolkit for data modeling, Tech. Rep. physics/0306116, SLAC, Stanford, CA, Jun, 2003. arXiv:physics/0306116v1 [physics.data-an].
[122] G. Cowan, K. Cranmer, E. Gross, and O. Vitells, Asymptotic formulae for likelihood-based tests of new physics , Eur. Phys.
[123] A. L. Read, Presentation of search results: The CLs technique , J. Phys. G28 (2002) 2693.
[124] E. Gross and O. Vitells, Trial factors for the look elsewhere effect in high energy physics , Eur. Phys. J. C 70 (2010) 525.
[125] ATLAS Collaboration, Search for a Standard Model Higgs boson in the H → ZZ → llνν decay channel using 4.7 fb−1 of
√
s = 7 TeV data with the ATLAS detector, submitted to Phys. Lett. B (2012), arXiv:1205.6744 [hep-ex].
[126] ATLAS Collaboration, Search for a Standard Model Higgs boson in the mass range 200-600 GeV in the H → ZZ → llqq decay channel, submitted to Phys. Lett. B (2012), arXiv:1206.2443 [hep-ex].
[127] ATLAS Collaboration, Observation of an excess of events in the search for the Standard Model Higgs boson in the gamma-gamma channel with the ATLAS detector, ATLAS-CONF-2012-091 (2012). http://cdsweb.cern.ch/record/1460410.
[128] ATLAS Collaboration, Search for the Higgs boson in the
√
H → WW → lν j j decay channel at s = 7 TeV with the
ATLAS detector, submitted to Phys. Lett. B (2012),
[129] ATLAS Collaboration, Search for the Standard Model Higgs boson in the H → τ+τ− decay mode in √ s = 7 TeV pp collisions with ATLAS, submitted to JHEP (2012), arXiv:1206.5971 [hep-ex].
[130] ATLAS Collaboration, Search for the Standard Model Higgs boson produced in association with a vector boson and decaying to a b-quark pair with the ATLAS detector, submitted to Phys. Lett. B (2012), arXiv:1207.0210 [hep-ex].
[131] ATLAS Collaboration, Observation of an Excess of Events in the Search for the Standard Model Higgs Boson in the H → WW(∗) → eνeν Channel with the ATLAS Detector, ATLAS-CONF-2012-098 (2012). http://cdsweb.cern.ch/record/1462530.
[132] ATLAS Collaboration, Luminosity Determination in pp
√
Collisions at s = 7 TeV Using the ATLAS Detector at the LHC, Eur. Phys. J. C 71 (2011) 1630, arXiv:1101.2185 [hep-ex].
[133] ATLAS Collaboration, Luminosity Determination in pp
√
Collisions at s = 7 TeV using the ATLAS Detector in 2011, ATLAS-CONF-2011-116 (2011). http://cdsweb.cern.ch/record/1376384.
[134] M. Botje, J. Butterworth, A. Cooper-Sarkar, A. de Roeck,
J. Feltesse, et al., The PDF4LHC Working Group Interim Recommendations, arXiv:1101.0538 [hep-ph].
[135] A. Martin, W. Stirling, R. Thorne, and G. Watt, Parton distributions for the LHC , Eur.Phys.J. C63 (2009) 189.
[136] R. D. Ball et al., Impact of heavy quark masses on parton distributions and LHC phenomenology, Nucl. Phys. B 849 (2011) 296, arXiv:1101.1300 [hep-ph].
[137] J. M. Campbell, R. K. Ellis, and G. Zanderighi, Next-to-Leading order Higgs + 2 jet production via gluon fusion , JHEP 0610 (2006) 028.
[138] ATLAS Collaboration, Observation of an Excess of Events in the Search for the Standard Model Higgs boson with the ATLAS detector at the LHC, ATLAS-CONF-2012-093 (2012). http://cdsweb.cern.ch/record/1460439.
[139] ATLAS and CMS Collaborations, Combined Standard Model Higgs boson searches with up to 2.3 fb−1 of pp collisions at
√
s=7 TeV at the LHC, ATLAS-CONF-2011-157,
CMS-PAS-HIG-11-023 (2011).
http://cdsweb.cern.ch/record/1399599.
[140] L. D. Landau, The moment of a 2-photon system, Dokl. Akad. Nawk. USSR 60 (1948) 207.
[141] C. N. Yang, Selection Rules for the Dematerialization of a Particle Into Two Photons , Phys. Rev. 77 (1950) 242.
G. Aad48, T. Abajyan21, B. Abbott111, J. Abdallah12, S. Abdel Khalek115, A.A. Abdelalim49, O. Abdinov11,
R. Aben105, B. Abi112, M. Abolins88, O.S. AbouZeid158, H. Abramowicz153, H. Abreu136, B.S. Acharya164a,164b,
L. Adamczyk38, D.L. Adams25, T.N. Addy56, J. Adelman176, S. Adomeit98, P. Adragna75, T. Adye129, S. Aefsky23,
I.A. Christidi77, A. Christov48, D. Chromek-Burckhart30, M.L. Chu151, J. Chudoba125, G. Ciapetti132a,132b,
A.K. Ciftci4a, R. Ciftci4a, D. Cinca34, V. Cindro74, C. Ciocca20a,20b, A. Ciocio15, M. Cirilli87, P. Cirkovic13b,
Z.H. Citron172, M. Citterio89a, M. Ciubancan26a, A. Clark49, P.J. Clark46, R.N. Clarke15, W. Cleland123,
124a
C. Collins-Tooth53, J. Collot55, T. Colombo119a,119b, G. Colon84, G. Compostella99, P. Conde Mui˜no,
E. Coniavitis166, M.C. Conidi12, S.M. Consonni89a,89b, V. Consorti48, S. Constantinescu26a, C. Conta119a,119b,
G. Conti57, F. Conventi102a, j, M. Cooke15, B.D. Cooper77, A.M. Cooper-Sarkar118, N.J. Cooper-Smith76, K. Copic15,
T. Cornelissen175, M. Corradi20a, F. Corriveau85,k, A. Cortes-Gonzalez165, G. Cortiana99, G. Costa89a, M.J. Costa167,
D. Costanzo139, D. Cotˆe´30, L. Courneyea169, G. Cowan76, C. Cowden28, B.E. Cox82, K. Cranmer108,
F. Crescioli122a,122b, M. Cristinziani21, G. Crosetti37a,37b, S. Cr´ep´e-Renaudin55, C.-M. Cuciuc26a,
C. Cuenca Almenar176, T. Cuhadar Donszelmann139, M. Curatolo47, C.J. Curtis18, C. Cuthbert150, P. Cwetanski60,
H. Czirr141, P. Czodrowski44, Z. Czyczula176, S. D’Auria53, M. D’Onofrio73, A. D’Orazio132a,132b,
C.N.P. Gee129, D.A.A. Geerts105, Ch. Geich-Gimbel21, K. Gellerstedt146a,146b, C. Gemme50a, A. Gemmell53,
A.R. Gillman129, D.M. Gingrich3,d, J. Ginzburg153, N. Giokaris9, M.P. Giordani164c, R. Giordano102a,102b,
F.M. Giorgi16, P. Giovannini99, P.F. Giraud136, D. Giugni89a, M. Giunta93, P. Giusti20a, B.K. Gjelsten117,
alez de la Hoz167, G. Gonzalez Parra12, M.L. Gonzalez Silva27, S. Gonzalez-Sevilla49
andez Jim´,
L. Hervas30, G.G. Hesketh77, N.P. Hessey105, E. Hig´on-Rodriguez167, J.C. Hill28, K.H. Hiller42, S. Hillert21,
S.J. Hillier18, I. Hinchliffe15, E. Hines120, M. Hirose116, F. Hirsch43, D. Hirschbuehl175, J. Hobbs148, N. Hod153,
83
A.V. Ivashin128, W. Iwanski39, H. Iwasaki65, J.M. Izen41, V. Izzo102a, B. Jackson120, J.N. Jackson73, P. Jackson1,
K.E. Johansson146a, P. Johansson139, S. Johnert42, K.A. Johns7, K. Jon-And146a,146b, G. Jones170, R.W.L. Jones71,
T.J. Jones73, C. Joram30, P.M. Jorge124a, K.D. Joshi82, J. Jovicevic147, T. Jovin13b, X. Ju173, C.A. Jung43,
T.J. Khoo28, G. Khoriauli21, A. Khoroshilov175, V. Khovanskiy95, E. Khramov64, J. Khubua51b, H. Kim146a,146b,
29
P. Kodys126, K. K¨onig104, S. Koenig81, L. K¨,
oneke30, A.C. K¨opke81, F. Koetsveld104, P. Koevesarki21, T. Koffas
E. Koffeman105, L.A. Kogan118, S. Kohlmann175, F. Kohn54, Z. Kohout127, T. Kohriki65, T. Koi143,
S.L. Lloyd75, E. Lobodzinska42, P. Loch7, W.S. Lockman137, T. Loddenkoetter21, F.K. Loebinger82, A. Loginov176,
M.J. Losty159a,∗, X. Lou41, A. Lounis115, K.F. Loureiro162, J. Love6, P.A. Love71, A.J. Lowe143,e, F. Lu33a,
cek74, J. Machado Miguens124a, R. Mackeprang36, R.J. Madaras15, H.J. Maddocks71
F.K. Martens158, L.F. Marti17, S. Marti-Garcia167, B. Martin30, B. Martin88, J.P. Martin93, T.A. Martin18,
V.J. Martin46, B. Martin dit Latour49, S. Martin-Haugh149, M. Martinez12, V. Martinez Outschoorn57,
R.P. Middleton129, S. Migas73, L. Mijovi´c136, G. Mikenberg172, M. Mikestikova125, M. Mikuˇz74, D.W. Miller31,
R.J. Miller88, W.J. Mills168, C. Mills57, A. Milov172, D.A. Milstead146a,146b, D. Milstein172, A.A. Minaenko128,
M.Mi˜,
nano Moya167, I.A. Minashvili64, A.I. Mincer108, B. Mindur38, M. Mineev64, Y. Ming173, L.M. Mir12
G. Mirabelli132a, J. Mitrevski137, V.A. Mitsou167, S. Mitsui65, P.S. Miyagawa139, J.U. Mj¨ornmark79, T. Moa146a,146b,
V. Moeller28, K. M¨oser21, S. Mohapatra148, W. Mohr48, R. Moles-Valls167, A. Molfetas30, J. Monk77
onig42, N. M¨,
E. Monnier83, J. Montejo Berlingen12, F. Monticelli70, S. Monzani20a,20b, R.W. Moore3, G.F. Moorhead86,
C. Mora Herrera49, A. Moraes53, N. Morange136, J. Morel54, G. Morello37a,37b, D. Moreno81, M. Moreno Ll´acer,
P. Morettini50a, M. Morgenstern44, M. Morii57, A.K. Morley30, G. Mornacchi30, J.D. Morris75, L. Morvaj101,
A.A. Nepomuceno24a, M. Nessi30,aa, M.S. Neubauer165, M. Neumann175, A. Neusiedl81, R.M. Neves108, P. Nevski25,
30
A. Policicchio37a,37b, R. Polifka158, A. Polini20a, J. Poll75, V. Polychronakos25, D. Pomeroy23, K. Pomm`,
es
L. Pontecorvo132a, B.G. Pope88, G.A. Popeneciu26a, D.S. Popovic13a, A. Poppleton30, X. Portell Bueso30,
A.S. Randle-Conde40, K. Randrianarivony29, F. Rauscher98, T.C. Rave48, M. Raymond30, A.L. Read117,
R.D. Schamberger148, A.G. Schamov107, V. Scharf58a, V.A. Schegelsky121, D. Scheirich87, M. Schernau163,
L.N. Smirnova97, O. Smirnova79, B.C. Smith57, D. Smith143, K.M. Smith53, M. Smizanska71, K. Smolek127,
G.A. Stewart30, J.A. Stillings21, M.C. Stockton85, K. Stoerig48, G. Stoicea26a, S. Stonjek99, P. Strachota126,
anchez167, D. Ta105, K. Tackmann42, A. Taffard163, R. Tafirout159a
N. Taiblum153, Y. Takahashi101, H. Takai25, R. Takashima68, H. Takeda66, T. Takeshita140, Y. Takubo65, M. Talby83,
A. Talyshev107, f , M.C. Tamsett25, K.G. Tan86, J. Tanaka155, R. Tanaka115, S. Tanaka131, S. Tanaka65,
A.J. Tanasijczuk142, K. Tani66, N. Tannoury83, S. Tapprogge81, D. Tardif158, S. Tarem152, F. Tarrade29,
G.F. Tartarelli89a, P. Tas126, M. Tasevsky125, E. Tassi37a,37b, M. Tatarkhanov15, Y. Tayalati135d, C. Taylor77,
R.J. Teuscher158,k, J. Therhaag21, T. Theveneaux-Pelzer78, S. Thoma48, J.P. Thomas18, E.N. Thompson35, P.D. Thompson18, P.D. Thompson158, A.S. Thompson53, L.A. Thomsen36, E. Thomson120, M. Thomson28,
o Pastor167, J. Toth83,ad, F. Touchard83
D.R. Tovey139, T. Trefzger174, L. Tremblet30, A. Tricoli30, I.M. Trigger159a, G. Trilling15, S. Trincaz-Duvoid78,
O.E. Vickey Boeriu145b, G.H.A. Viehhauser118, S. Viel168, M. Villa20a,20b, M. Villaplana Perez167, E. Vilucchi47,
P.M. Watkins18, A.T. Watson18, I.J. Watson150, M.F. Watson18, G. Watts138, S. Watts82, A.T. Waugh150,
H.G. Wilkens30, J.Z. Will98, E. Williams35, H.H. Williams120, W. Willis35, S. Willocq84, J.A. Wilson18,
M.G. Wilson143, A. Wilson87, I. Wingerter-Seez5, S. Winkelmann48, F. Winklmeier30, M. Wittgen143,
S.J. Wollstadt81, M.W. Wolter39, H. Wolters124a,h, W.C. Wong41, G. Wooden87, B.K. Wosiek39, J. Wotschack30,
Zeniˇ,
D. Zhang33b,ak, H. Zhang88, J. Zhang6, X. Zhang33d, Z. Zhang115, L. Zhao108, Z. Zhao33b, A. Zhemchugov64,
J. Zhong118, B. Zhou87, N. Zhou163, Y. Zhou151, C.G. Zhu33d, H. Zhu42, J. Zhu87, Y. Zhu33b, X. Zhuang98,
V. Zhuravlov99, D. Zieminska60, N.I. Zimin64, R. Zimmermann21, S. Zimmermann21, S. Zimmermann48, M. Ziolkowski141, R. Zitoun5, L. ˇc35, V.V. Zmouchko128,∗, G. Zobernig173, A. Zoccoli20a,20b
Zivkovi´,
M. zur Nedden16, V. Zutshi106, L. Zwalinski30.
1 School of Chemistry and Physics, University of Adelaide, Adelaide, Australia 2 Physics Department, SUNY Albany, Albany NY, United States of America 3 Department of Physics, University of Alberta, Edmonton AB, Canada 4(a)Department of Physics, Ankara University, Ankara; (b)Department of Physics, Dumlupinar University, Kutahya; (c)Department of Physics, Gazi University, Ankara; (d)Division of Physics, TOBB University of Economics and Technology, Ankara; (e)Turkish Atomic Energy Authority, Ankara, Turkey 5 LAPP, CNRS/IN2P3 and Universit´e de Savoie, Annecy-le-Vieux, France 6 High Energy Physics Division, Argonne National Laboratory, Argonne IL, United States of America 7 Department of Physics, University of Arizona, Tucson AZ, United States of America 8 Department of Physics, The University of Texas at Arlington, Arlington TX, United States of America 9 Physics Department, University of Athens, Athens, Greece 10 Physics Department, National Technical University of Athens, Zografou, Greece 11 Institute of Physics, Azerbaijan Academy of Sciences, Baku, Azerbaijan 12 Institut de F´ısica d’Altes Energies and Departament de F´ısica de la Universitat Aut`onoma de Barcelona and ICREA, Barcelona, Spain 13 (a)Institute of Physics, University of Belgrade, Belgrade; (b)Vinca Institute of Nuclear Sciences, University of Belgrade, Belgrade, Serbia 14 Department for Physics and Technology, University of Bergen, Bergen, Norway 15 Physics Division, Lawrence Berkeley National Laboratory and University of California, Berkeley CA, United States of America 16 Department of Physics, Humboldt University, Berlin, Germany 17 Albert Einstein Center for Fundamental Physics and Laboratory for High Energy Physics, University of Bern, Bern, Switzerland 18 School of Physics and Astronomy, University of Birmingham, Birmingham, United Kingdom 19 (a)Department of Physics, Bogazici University, Istanbul; (b)Division of Physics, Dogus University, Istanbul; (c)Department of Physics Engineering, Gaziantep University, Gaziantep; (d)Department of Physics, Istanbul Technical University, Istanbul, Turkey 20 (a)INFN Sezione di Bologna; (b)Dipartimento di Fisica, Universit`a di Bologna, Bologna, Italy 21 Physikalisches Institut, University of Bonn, Bonn, Germany 22 Department of Physics, Boston University, Boston MA, United States of America 23 Department of Physics, Brandeis University, Waltham MA, United States of America 24 (a)Universidade Federal do Rio De Janeiro COPPE/EE/IF, Rio de Janeiro; (b)Federal University of Juiz de Fora (UFJF), Juiz de Fora; (c)Federal University of Sao Joao del Rei (UFSJ), Sao Joao del Rei; (d)Instituto de Fisica, Universidade de Sao Paulo, Sao Paulo, Brazil 25 Physics Department, Brookhaven National Laboratory, Upton NY, United States of America 26 (a)National Institute of Physics and Nuclear Engineering, Bucharest; (b)University Politehnica Bucharest, Bucharest; (c)West University in Timisoara, Timisoara, Romania 27 Departamento de F´ısica, Universidad de Buenos Aires, Buenos Aires, Argentina 28 Cavendish Laboratory, University of Cambridge, Cambridge, United Kingdom 29 Department of Physics, Carleton University, Ottawa ON, Canada 30 CERN, Geneva, Switzerland 31 Enrico Fermi Institute, University of Chicago, Chicago IL, United States of America 32 (a)Departamento de F´ısica, Pontificia Universidad Cat´olica de Chile, Santiago; (b)Departamento de F´ısica, Universidad T´ecnica Federico Santa Mar´ıa, Valpara´ıso, Chile 33 (a)Institute of High Energy Physics, Chinese Academy of Sciences, Beijing; (b)Department of Modern Physics, University of Science and Technology of China, Anhui; (c)Department of Physics, Nanjing University, Jiangsu; (d)School of Physics, Shandong University, Shandong, China 34 Laboratoire de Physique Corpusculaire, Clermont Universit´e Blaise Pascal and CNRS/IN2P3,
e and Universit´Clermont-Ferrand, France 35 Nevis Laboratory, Columbia University, Irvington NY, United States of America 36 Niels Bohr Institute, University of Copenhagen, Kobenhavn, Denmark 37 (a)INFN Gruppo Collegato di Cosenza; (b)Dipartimento di Fisica, Universit`
a della Calabria, Arcavata di Rende, Italy 38 AGH University of Science and Technology, Faculty of Physics and Applied Computer Science, Krakow, Poland 39 The Henryk Niewodniczanski Institute of Nuclear Physics, Polish Academy of Sciences, Krakow, Poland 40 Physics Department, Southern Methodist University, Dallas TX, United States of America 41 Physics Department, University of Texas at Dallas, Richardson TX, United States of America 42 DESY, Hamburg and Zeuthen, Germany 43 Institut f¨ur Experimentelle Physik IV, Technische Universit¨at Dortmund, Dortmund, Germany 44 Institut f¨
ur Kern-und Teilchenphysik, Technical University Dresden, Dresden, Germany 45 Department of Physics, Duke University, Durham NC, United States of America 46 SUPA -School of Physics and Astronomy, University of Edinburgh, Edinburgh, United Kingdom 47 INFN Laboratori Nazionali di Frascati, Frascati, Italy 48 Fakult¨at f¨ur Mathematik und Physik, Albert-Ludwigs-Universit¨at, Freiburg, Germany 49 Section de Physique, Universit´eve, Geneva, Switzerland
e de Gen`50 (a)INFN Sezione di Genova; (b)Dipartimento di Fisica, Universit`a di Genova, Genova, Italy 51 (a)E. Andronikashvili Institute of Physics, Tbilisi State University, Tbilisi; (b)High Energy Physics Institute, Tbilisi State University, Tbilisi, Georgia 52 II Physikalisches Institut, Justus-Liebig-Universit¨at Giessen, Giessen, Germany 53 SUPA -School of Physics and Astronomy, University of Glasgow, Glasgow, United Kingdom 54 II Physikalisches Institut, Georg-August-Universit¨at, G¨ottingen, Germany 55 Laboratoire de Physique Subatomique et de Cosmologie, Universit´e Joseph Fourier and CNRS/IN2P3 and Institut National Polytechnique de Grenoble, Grenoble, France 56 Department of Physics, Hampton University, Hampton VA, United States of America 57 Laboratory for Particle Physics and Cosmology, Harvard University, Cambridge MA, United States of America 58 (a)Kirchhoff-Institut f¨at Heidelberg, Heidelberg; (b)Physikalisches Institut,
ur Physik, Ruprecht-Karls-Universit¨Ruprecht-Karls-Universit¨at Heidelberg, Heidelberg; (c)ZITI Institut f¨ur technische Informatik, Ruprecht-Karls-Universit¨
at Heidelberg, Mannheim, Germany 59 Faculty of Applied Information Science, Hiroshima Institute of Technology, Hiroshima, Japan 60 Department of Physics, Indiana University, Bloomington IN, United States of America 61 Institut f¨at, Innsbruck, Austria
ur Astro-und Teilchenphysik, Leopold-Franzens-Universit¨62 University of Iowa, Iowa City IA, United States of America 63 Department of Physics and Astronomy, Iowa State University, Ames IA, United States of America 64 Joint Institute for Nuclear Research, JINR Dubna, Dubna, Russia 65 KEK, High Energy Accelerator Research Organization, Tsukuba, Japan 66 Graduate School of Science, Kobe University, Kobe, Japan 67 Faculty of Science, Kyoto University, Kyoto, Japan 68 Kyoto University of Education, Kyoto, Japan 69 Department of Physics, Kyushu University, Fukuoka, Japan 70 Instituto de F´ısica La Plata, Universidad Nacional de La Plata and CONICET, La Plata, Argentina 71 Physics Department, Lancaster University, Lancaster, United Kingdom 72 (a)INFN Sezione di Lecce; (b)Dipartimento di Matematica e Fisica, Universit`
a del Salento, Lecce, Italy 73 Oliver Lodge Laboratory, University of Liverpool, Liverpool, United Kingdom 74 Department of Physics, Joˇ
zef Stefan Institute and University of Ljubljana, Ljubljana, Slovenia 75 School of Physics and Astronomy, Queen Mary University of London, London, United Kingdom 76 Department of Physics, Royal Holloway University of London, Surrey, United Kingdom 77 Department of Physics and Astronomy, University College London, London, United Kingdom 78 Laboratoire de Physique Nucl´eaire et de Hautes Energies, UPMC and Universit´e Paris-Diderot and CNRS/IN2P3, Paris, France 79 Fysiska institutionen, Lunds universitet, Lund, Sweden 80 Departamento de Fisica Teorica C-15, Universidad Autonoma de Madrid, Madrid, Spain 81 Institut f¨ur Physik, Universit¨at Mainz, Mainz, Germany 82 School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom 83 CPPM, Aix-Marseille Universit´
e and CNRS/IN2P3, Marseille, France
84 Department of Physics, University of Massachusetts, Amherst MA, United States of America
85 Department of Physics, McGill University, Montreal QC, Canada
86 School of Physics, University of Melbourne, Victoria, Australia
87 Department of Physics, The University of Michigan, Ann Arbor MI, United States of America
88 Department of Physics and Astronomy, Michigan State University, East Lansing MI, United States of America
89 (a)INFN Sezione di Milano; (b)Dipartimento di Fisica, Universit`a di Milano, Milano, Italy
90 B.I. Stepanov Institute of Physics, National Academy of Sciences of Belarus, Minsk, Republic of Belarus
91 National Scientific and Educational Centre for Particle and High Energy Physics, Minsk, Republic of Belarus
92 Department of Physics, Massachusetts Institute of Technology, Cambridge MA, United States of America
93 Group of Particle Physics, University of Montreal, Montreal QC, Canada
94 P.N. Lebedev Institute of Physics, Academy of Sciences, Moscow, Russia
95 Institute for Theoretical and Experimental Physics (ITEP), Moscow, Russia
96 Moscow Engineering and Physics Institute (MEPhI), Moscow, Russia
97 Skobeltsyn Institute of Nuclear Physics, Lomonosov Moscow State University, Moscow, Russia
98 Fakult¨at f¨ur Physik, Ludwig-Maximilians-Universit¨at M¨unchen, M¨unchen, Germany
99 Max-Planck-Institut f¨unchen, Germany
ur Physik (Werner-Heisenberg-Institut), M¨
100 Nagasaki Institute of Applied Science, Nagasaki, Japan
101 Graduate School of Science and Kobayashi-Maskawa Institute, Nagoya University, Nagoya, Japan
102 (a)INFN Sezione di Napoli; (b)Dipartimento di Scienze Fisiche, Universit`a di Napoli, Napoli, Italy
103 Department of Physics and Astronomy, University of New Mexico, Albuquerque NM, United States of America
104 Institute for Mathematics, Astrophysics and Particle Physics, Radboud University Nijmegen/Nikhef, Nijmegen,
Netherlands
105 Nikhef National Institute for Subatomic Physics and University of Amsterdam, Amsterdam, Netherlands
106 Department of Physics, Northern Illinois University, DeKalb IL, United States of America
107 Budker Institute of Nuclear Physics, SB RAS, Novosibirsk, Russia
108 Department of Physics, New York University, New York NY, United States of America
109 Ohio State University, Columbus OH, United States of America
110 Faculty of Science, Okayama University, Okayama, Japan
111 Homer L. Dodge Department of Physics and Astronomy, University of Oklahoma, Norman OK, United States of
America
112 Department of Physics, Oklahoma State University, Stillwater OK, United States of America
113 Palack´y University, RCPTM, Olomouc, Czech Republic
114 Center for High Energy Physics, University of Oregon, Eugene OR, United States of America
115 LAL, Universit´e Paris-Sud and CNRS/IN2P3, Orsay, France
116 Graduate School of Science, Osaka University, Osaka, Japan
117 Department of Physics, University of Oslo, Oslo, Norway
118 Department of Physics, Oxford University, Oxford, United Kingdom
119 (a)INFN Sezione di Pavia; (b)Dipartimento di Fisica, Universit`a di Pavia, Pavia, Italy
120 Department of Physics, University of Pennsylvania, Philadelphia PA, United States of America
121 Petersburg Nuclear Physics Institute, Gatchina, Russia
122 (a)INFN Sezione di Pisa; (b)Dipartimento di Fisica E. Fermi, Universit`a di Pisa, Pisa, Italy
123 Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh PA, United States of America
124 (a)Laboratorio de Instrumentacao e Fisica Experimental de Particulas -LIP, Lisboa, Portugal; (b)Departamento de
Fisica Teorica y del Cosmos and CAFPE, Universidad de Granada, Granada, Spain
125 Institute of Physics, Academy of Sciences of the Czech Republic, Praha, Czech Republic
126 Faculty of Mathematics and Physics, Charles University in Prague, Praha, Czech Republic
127 Czech Technical University in Prague, Praha, Czech Republic
128 State Research Center Institute for High Energy Physics, Protvino, Russia
129 Particle Physics Department, Rutherford Appleton Laboratory, Didcot, United Kingdom 130 Physics Department, University of Regina, Regina SK, Canada 131 Ritsumeikan University, Kusatsu, Shiga, Japan 132 (a)INFN Sezione di Roma I; (b)Dipartimento di Fisica, Universit`a La Sapienza, Roma, Italy
133 (a)INFN Sezione di Roma Tor Vergata; (b)Dipartimento di Fisica, Universit`a di Roma Tor Vergata, Roma, Italy
134 (a)INFN Sezione di Roma Tre; (b)Dipartimento di Fisica, Universit`a Roma Tre, Roma, Italy
135 (a)Facult´e des Sciences Ain Chock, R´eseau Universitaire de Physique des Hautes Energies -Universit´e Hassan II,
Casablanca; (b)Centre National de l’Energie des Sciences Techniques Nucleaires, Rabat; (c)Facult´e des Sciences
Semlalia, Universit´e des Sciences, Universit´e Cadi Ayyad, LPHEA-Marrakech; (d)Facult´e Mohamed Premier and
LPTPM, Oujda; (e)Facult´e des sciences, Universit´e Mohammed V-Agdal, Rabat, Morocco
136 DSM/IRFU (Institut de Recherches sur les Lois Fondamentales de l’Univers), CEA Saclay (Commissariat a
l’Energie Atomique), Gif-sur-Yvette, France
137 Santa Cruz Institute for Particle Physics, University of California Santa Cruz, Santa Cruz CA, United States of
America
138 Department of Physics, University of Washington, Seattle WA, United States of America
139 Department of Physics and Astronomy, University of Sheffield, Sheffield, United Kingdom
140 Department of Physics, Shinshu University, Nagano, Japan
141 Fachbereich Physik, Universit¨at Siegen, Siegen, Germany
142 Department of Physics, Simon Fraser University, Burnaby BC, Canada
143 SLAC National Accelerator Laboratory, Stanford CA, United States of America
144 (a)Faculty of Mathematics, Physics & Informatics, Comenius University, Bratislava; (b)Department of Subnuclear
Physics, Institute of Experimental Physics of the Slovak Academy of Sciences, Kosice, Slovak Republic
145 (a)Department of Physics, University of Johannesburg, Johannesburg; (b)School of Physics, University of the
Witwatersrand, Johannesburg, South Africa
146 (a)Department of Physics, Stockholm University; (b)The Oskar Klein Centre, Stockholm, Sweden
147 Physics Department, Royal Institute of Technology, Stockholm, Sweden
148 Departments of Physics & Astronomy and Chemistry, Stony Brook University, Stony Brook NY, United States of
America
149 Department of Physics and Astronomy, University of Sussex, Brighton, United Kingdom
150 School of Physics, University of Sydney, Sydney, Australia
151 Institute of Physics, Academia Sinica, Taipei, Taiwan
152 Department of Physics, Technion: Israel Institute of Technology, Haifa, Israel
153 Raymond and Beverly Sackler School of Physics and Astronomy, Tel Aviv University, Tel Aviv, Israel
154 Department of Physics, Aristotle University of Thessaloniki, Thessaloniki, Greece
155 International Center for Elementary Particle Physics and Department of Physics, The University of Tokyo, Tokyo,
Japan
156 Graduate School of Science and Technology, Tokyo Metropolitan University, Tokyo, Japan
157 Department of Physics, Tokyo Institute of Technology, Tokyo, Japan
158 Department of Physics, University of Toronto, Toronto ON, Canada
159 (a)TRIUMF, Vancouver BC; (b)Department of Physics and Astronomy, York University, Toronto ON, Canada
160 Faculty of Pure and Applied Sciences, University of Tsukuba, Tsukuba, Japan
161 Department of Physics and Astronomy, Tufts University, Medford MA, United States of America
162 Centro de Investigaciones, Universidad Antonio Narino, Bogota, Colombia
163 Department of Physics and Astronomy, University of California Irvine, Irvine CA, United States of America
164 (a)INFN Gruppo Collegato di Udine; (b)ICTP, Trieste; (c)Dipartimento di Chimica, Fisica e Ambiente, Universit`a
di Udine, Udine, Italy
165 Department of Physics, University of Illinois, Urbana IL, United States of America
166 Department of Physics and Astronomy, University of Uppsala, Uppsala, Sweden
167 Instituto de F´ısica Corpuscular (IFIC) and Departamento de F´omica, Molecular y Nuclear and ısica At´
Departamento de Ingenier´ıa Electr´onica and Instituto de Microelectr´onica de Barcelona (IMB-CNM), University of
Valencia and CSIC, Valencia, Spain
168 Department of Physics, University of British Columbia, Vancouver BC, Canada
169 Department of Physics and Astronomy, University of Victoria, Victoria BC, Canada 170 Department of Physics, University of Warwick, Coventry, United Kingdom 171 Waseda University, Tokyo, Japan 172 Department of Particle Physics, The Weizmann Institute of Science, Rehovot, Israel 173 Department of Physics, University of Wisconsin, Madison WI, United States of America 174 Fakult¨ur Physik und Astronomie, Julius-Maximilians-Universit¨urzburg, Germany
atf¨at,W¨175 Fachbereich C Physik, Bergische Universit¨at Wuppertal, Wuppertal, Germany 176 Department of Physics, Yale University, New Haven CT, United States of America 177 Yerevan Physics Institute, Yerevan, Armenia 178 Centre de Calcul de l’Institut National de Physique Nucl´eaire et de Physique des Particules (IN2P3), Villeurbanne, France a Also at Laboratorio de Instrumentacao e Fisica Experimental de Particulas -LIP, Lisboa, Portugal b Also at Faculdade de Ciencias and CFNUL, Universidade de Lisboa, Lisboa, Portugal c Also at Particle Physics Department, Rutherford Appleton Laboratory, Didcot, United Kingdom d Also at TRIUMF, Vancouver BC, Canada e Also at Department of Physics, California State University, Fresno CA, United States of America f Also at Novosibirsk State University, Novosibirsk, Russia g Also at Fermilab, Batavia IL, United States of America h Also at Department of Physics, University of Coimbra, Coimbra, Portugal i Also at Department of Physics, UASLP, San Luis Potosi, Mexico j Also at Universit`a di Napoli Parthenope, Napoli, Italy k Also at Institute of Particle Physics (IPP), Canada l Also at Department of Physics, Middle East Technical University, Ankara, Turkey m Also at Louisiana Tech University, Ruston LA, United States of America n Also at Dep Fisica and CEFITEC of Faculdade de Ciencias e Tecnologia, Universidade Nova de Lisboa, Caparica, Portugal
o Also at Department of Physics and Astronomy, University College London, London, United Kingdom p Also at Department of Physics, University of Cape Town, Cape Town, South Africa q Also at Institute of Physics, Azerbaijan Academy of Sciences, Baku, Azerbaijan r Also at Institut f¨at Hamburg, Hamburg, Germany
ur Experimentalphysik, Universit¨s Also at Manhattan College, New York NY, United States of America t Also at CPPM, Aix-Marseille Universit´e and CNRS/IN2P3, Marseille, France u Also at School of Physics and Engineering, Sun Yat-sen University, Guanzhou, China v Also at Academia Sinica Grid Computing, Institute of Physics, Academia Sinica, Taipei, Taiwan w Also at Laboratoire de Physique Nucl´eaire et de Hautes Energies, UPMC and Universit´e Paris-Diderot and CNRS/IN2P3, Paris, France x Also at School of Physics, Shandong University, Shandong, China y Also at Dipartimento di Fisica, Universit`
a La Sapienza, Roma, Italy
z Also at DSM/IRFU (Institut de Recherches sur les Lois Fondamentales de l’Univers), CEA Saclay (Commissariat a
l’Energie Atomique), Gif-sur-Yvette, France
aa Also at Section de Physique, Universit´eve, Geneva, Switzerland e de Gen`
ab Also at Departamento de Fisica, Universidade de Minho, Braga, Portugal
ac Also at Department of Physics and Astronomy, University of South Carolina, Columbia SC, United States of
America
ad Also at Institute for Particle and Nuclear Physics, Wigner Research Centre for Physics, Budapest, Hungary
ae Also at California Institute of Technology, Pasadena CA, United States of America
af Also at Institute of Physics, Jagiellonian University, Krakow, Poland
ag Also at LAL, Universit´e Paris-Sud and CNRS/IN2P3, Orsay, France
ah Also at Nevis Laboratory, Columbia University, Irvington NY, United States of America
ai Also at Department of Physics and Astronomy, University of Sheffield, Sheffield, United Kingdom
aj Also at Department of Physics, Oxford University, Oxford, United Kingdom
ak Also at Institute of Physics, Academia Sinica, Taipei, Taiwan
al Also at Department of Physics, The University of Michigan, Ann Arbor MI, United States of America
am Also at Discipline of Physics, University of KwaZulu-Natal, Durban, South Africa
∗ Deceased