Published for SISSA by

Springer

Received: September 23, 2014
Revised: November 11, 2014
Accepted: November 21, 2014
Published: January 14, 2015

The ATLAS collaboration
E-mail: atlas.publications@cern.ch
Abstract: A search for the bb̄ decay of the Standard Model Higgs boson is performed
with the ATLAS experiment using the full dataset recorded at the LHC in Run 1. The√
integrated luminosities used are 4.7 and 20.3 fb−1 from pp collisions at s = 7 and 8 TeV,
respectively. The processes considered are associated (W/Z)H production, where W →eν/µν, Z → ee/µµ and Z → νν. The observed (expected) deviation from the backgroundonly hypothesis corresponds to a significance of 1.4 (2.6) standard deviations and the ratioof the measured signal yield to the Standard Model expectation is found to be µ = 0.52 ±0.32 (stat.) ± 0.24 (syst.) for a Higgs boson mass of 125.36 GeV. The analysis procedure isvalidated by a measurement of the yield of (W/Z)Z production with Z → bb̄ in the same
final states as for the Higgs boson search, from which the ratio of the observed signal yield
to the Standard Model expectation is found to be 0.74 ± 0.09 (stat.) ± 0.14 (syst.).
Keywords: Hadron-Hadron Scattering
ArXiv ePrint: 1409.6212

Open Access, Copyright CERN,
for the benefit of the ATLAS Collaboration.
Article funded by SCOAP3 .

doi:10.1007/JHEP01(2015)069

JHEP01(2015)069

Search for the bb̄ decay of the Standard Model Higgs
boson in associated (W/Z)H production with the
ATLAS detector

Contents
1

2 The ATLAS detector

3

3 Data and simulated samples

3

4 Object reconstruction

6

5 Event selection

8

6 Multivariate analysis

14

7 Background composition and modelling
7.1 Multijet background
7.2 Corrections to the simulation
7.3 Distributions in the dijet-mass analysis
7.4 Distributions in the multivariate analysis

21
21
23
25
29

8 Systematic uncertainties
8.1 Experimental uncertainties
8.2 Uncertainties on the multijet backgrounds
8.3 Uncertainties on the modelling of the simulated backgrounds
8.4 Uncertainties on the signal modelling

34
34
35
36
40

9 Statistical procedure
9.1 General aspects
9.2 Technical details
9.3 Cross checks using diboson production

42
42
44
46

10 Analysis of the 7 TeV data
10.1 Object reconstruction
10.2 Event selection
10.3 Background composition and modelling10.4 Systematic uncertainties
10.5 Statistical procedure

50
50
50
50
51
51

11 Results
11.1 Nominal results
11.2 Cross-check with the dijet-mass analysis
11.3 Cross-check with the diboson analysis

53
53
55
58

12 Summary

60

–i–

JHEP01(2015)069

1 IntroductionA Tables of event yields

61

The ATLAS collaboration

71

1

Introduction–1–

JHEP01(2015)069

For decades, the Higgs boson [1–4] of the Standard Model (SM) remained an unconfirmed
prediction. In July 2012, the ATLAS and CMS experiments at the LHC reported theobservation of a new particle with a mass of about 125 GeV and with properties consistent
with those expected for the SM Higgs boson [5, 6]. Since then, more precise measurements
have strengthened the hypothesis that the new particle is indeed a Higgs boson [7–9].
These measurements, however, have been mainly performed in the bosonic decay modesof the new particle (H → γγ, H → ZZ, and H → W W ). It is also essential to verify
whether it decays into fermions as predicted by the Standard Model. Recently, the CMS
Collaboration reported evidence for the τ τ decay mode of the Higgs boson at a level ofsignificance of 3.4 standard deviations (σ) for mH = 125 GeV [10].
The H → bb decay mode is predicted in the SM to have a branching ratio of 58%
for mH = 125 GeV [11]. Accessing H → bb decays is therefore crucial for constraining,under fairly general assumptions [12], the overall Higgs boson decay width and, in a global
fit to all accessible combinations of Higgs boson production and decay modes, to allow
for measurements of absolute Higgs boson couplings. An inclusive search for H → bb is
not feasible at hadron colliders because of the overwhelming background from multijet
production. In spite of a cross section more than an order of magnitude lower than thedominant gluon-fusion process, associated production of a Higgs boson with a vector boson,
W or Z [13], offers a viable alternative because leptonic decays of the vector boson, W → `ν,Z → `` (` = e, µ), and Z → νν, can be efficiently used for triggering and background
reduction purposes [14, 15]. The CDF and D0 experiments at the Tevatron reported an
excess of events in their search for associated (W/Z)H production in the H → bb decay
mode at a significance level of 2.8σ for mH = 125 GeV [16]. Recently, the CMS experimentreported an excess of events in the H → bb decay mode with a significance of 2.1σ formH = 125 GeV [17].
In this paper, a search for associated (W/Z)H production of the SM Higgs bosonin the bb decay mode is presented, using the full integrated luminosity accumulated by
ATLAS during Run 1 of the LHC: 4.7 and 20.3 fb−1 from proton-proton (pp) collisions at
centre-of-mass energies of 7 and 8 TeV in 2011 and 2012, respectively. An analysis of the
7 TeV dataset has already been published by ATLAS [18]. In addition to the increase in
the amount of data analysed, the update presented in this paper benefits from numerous
analysis improvements. Some of the improvements to the object reconstruction, however,are available only for the 8 TeV dataset, which leads to separate analysis strategies for thetwo datasets.

–2–

JHEP01(2015)069

The analysis is performed for events containing zero, one, or two charged leptons(electrons or muons), targeting the Z → νν, W → `ν, or Z → `` decay modes of the vectorboson, respectively. In addition to Z → νν decays, the 0-lepton channel has a smaller but
not insignificant contribution from leptonic W decays when the lepton is produced outsideof the detector acceptance or not identified. A b-tagging algorithm is used to identify thejets consistent with originating from an H → bb decay. To improve the sensitivity, the threechannels are each split according to the vector-boson transverse momentum, the number ofjets (two or three), and the number of b-tagged jets. Topological and kinematic selection
criteria are applied within each of the resulting categories.
A binned maximum likelihood fit is used to extract the signal yield and the backgroundnormalisations. Systematic uncertainties on the signal and background modelling are implemented as deviations in their respective models in the form of “nuisance” parametersthat are varied in the fit. Each nuisance parameter is constrained by a penalty term in
the likelihood, associated with its uncertainty. Two versions of the analysis are presented
in this paper: in the first, referred to as the dijet-mass analysis, the mass of the dijet
system of b-tagged jets is the final discriminating variable used in the statistical analysis;in the other, a multivariate analysis (MVA) incorporating various kinematic variables in
addition to the dijet mass, as well as b-tagging information, provides the final discriminating variable. Because the latter information is not available in similar detail for the
7 TeV dataset, the MVA is used only for the 8 TeV dataset. In both analyses, dedicatedcontrol samples, typically with loosened b-tagging requirements, constrain the contributions of the dominant background processes. The most significant background sources are(W/Z)+heavy-flavour-jet production and tt production. The normalisations of these backgrounds are fully determined by the likelihood fit. Other significant background sources
are single-top-quark and diboson (W Z and ZZ) production, with normalisations takenfrom theory, as well as multijet events, normalised using multijet-enriched control samples.
Since the MVA has higher expected sensitivity, it is chosen as the nominal analysis for the
8 TeV dataset to extract the final results. To validate the analysis procedures, both forthe dijet-mass and MVA approaches, a measurement of the yield of (W/Z)Z production is
performed in the same final states and with the same event selection, with H → bb replacedby Z → bb.
This paper is organised as follows. A brief description of the ATLAS detector isgiven in section 2. Details of the data and simulated samples used in this analysis areprovided in section 3. This is followed by sections describing the dijet-mass and multivariate
analyses applied to the 8 TeV data. The reconstruction of physics objects such as leptons
and jets is addressed in section 4. Section 5 details the event selections applied to thedijet-mass and multivariate analyses, while section 6 explains the construction of the final
discriminating variable of the MVA. Section 7 discusses the background composition in thevarious analysis regions, while the systematic uncertainties are addressed in section 8. The
statistical procedure used to extract the results is described in section 9. For the 7 TeV
data, only a dijet-mass analysis is used, and differences with respect to the 8 TeV data
analysis are specified in section 10. The results are presented and discussed in section 11,
and a summary of the paper is given in section 12.

2

The ATLAS detector

3

Data and simulated samples

The datasets used in this analysis include only pp collision data recorded in stable beam conditions and with all relevant sub-detectors providing high-quality data. The corresponding
integrated luminosities are 4.7 and 20.3 fb−1 [20] for the 7 and 8 TeV data, respectively.
miss of
Events in the 0-lepton channel are selected by triggers based on the magnitude ET
missthe missing transverse momentum vector. The ET trigger configuration evolved during
data taking to cope with the increasing luminosity, and the trigger efficiency was improved
miss trigger efficiency on the E miss reconfor the 8 TeV data. The dependence of the ETTstructed offline is measured in W → µν+jets and Z → µµ+jets events collected with
miss calculated without the muon contribution. As
single-muon triggers, with the offline ET
miss triggers were not available forthere was a brief period of data-taking in which the ET
1

ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in thecentre of the detector and the z-axis coinciding with the axis of the beam pipe. The x-axis points from the IP
towards the centre of the LHC ring, and the y-axis points upward. Cylindrical coordinates (r,φ) are used in
the transverse plane, φ being the azimuthal angle around the z-axis. The pseudorapidity
is defined in terms
pof the polar angle θ as η = − ln tan(θ/2). The distance in (η,φ) coordinates, ∆R = (∆φ)2 + (∆η)2 , is also
used to define cone sizes. Transverse momentum and energy are defined as pT = p sin θ and ET = E sin θ,
respectively. For the purpose of object selections, η is calculated relative to the geometric centre of thedetector; otherwise, it is relative to the reconstructed primary vertex of each event.

–3–

JHEP01(2015)069

The ATLAS detector [19] is cylindrically symmetric around the beam axis and is structured
in a barrel and two endcaps. It consists of three main subsystems. The inner tracking detector is immersed in the 2 T axial magnetic field produced by a superconducting solenoid.
Charged-particle position and momentum measurements are made by pixel detectors followed by silicon-strip detectors in the pseudorapidity1 range |η| < 2.5 and by a straw-tube
transition-radiation tracker (TRT) in the range |η| < 2.0. The pixel detectors are crucial
for b-tagging, and the TRT also contributes to electron identification. The calorimeters, located beyond the solenoid, cover the range |η| < 4.9 with a variety of detector technologies.
The liquid-argon electromagnetic calorimeters are divided into barrel (|η| < 1.475), endcap(1.375 < |η| < 3.2), and forward (3.1 < |η| < 4.9) sections. The hadronic calorimeters(using scintillator tiles or liquid argon as active materials) surround the electromagnetic
calorimeters with a coverage of |η| < 4.9. The muon spectrometer measures the deflectionof muon tracks in the field of three large air-core toroidal magnets, each containing eightsuperconducting coils. It is instrumented with separate trigger and high-precision tracking
chambers covering the |η| < 2.4 and |η| < 2.7 ranges, respectively.
The trigger system is organised in three levels. The first level is based on custom-madehardware and uses coarse-granularity calorimeter and muon information. The second andthird levels are implemented as software algorithms and use the full detector granularity.
At the second level, only regions deemed interesting at the first level are analysed, while thethird level, called the event filter, makes use of the full detector read-out to reconstruct andselect events, which are then logged for offline analysis at a rate of up to 400 Hz averaged
over an accelerator fill.

–4–

JHEP01(2015)069

the first bunch crossings of two bunch trains, the integrated luminosity for the 0-lepton
channel in the 7 TeV dataset is reduced to 4.6 fb−1 . Events in the 1-lepton channel are
primarily selected by single-lepton triggers. The ET threshold of the single-electron trigger was raised from 20 to 22 GeV during the 7 TeV data-taking period, and to 24 GeV for
the 8 TeV data. The pT threshold of the single-muon trigger was similarly increased from18 GeV for the 7 TeV data to 24 GeV at 8 TeV. As the single-lepton triggers for the 8 TeV
data include isolation criteria, triggers with higher thresholds (60 GeV for electrons and
36 GeV for muons) but no isolation requirements are used in addition. Single-lepton trigger
efficiencies are measured using a tag-and-probe method applied to Z → ee and Z → µµ
miss triggers are also used to compensate for the limevents. In the 1-muon sub-channel, ET
ited muon trigger-chamber coverage in some regions of the detector. Events in the 2-lepton
channel are selected by a combination of single-lepton, dielectron and dimuon triggers. The
thresholds of the dilepton triggers are 12 GeV for electrons and 13 GeV for muons.
Monte Carlo (MC) simulated samples are produced for signal and background processes using the atlfast-II simulation [21], which includes a full simulation of the ATLAS
detector based on the geant4 program [22], except for the response of the calorimeters
for which a parameterised simulation is used. A list of the generators used for signal andbackground simulations is given in table 1.
The MC generator used for qq-initiated W H and ZH production is pythia8 [23] with
the CTEQ6L1 [24] parton distribution functions (PDFs). The AU2 tune [25, 26] is used
for the parton shower, hadronisation, and multiple parton interactions. The photos program [27] is used for QED final-state radiation. The powheg generator [28–30] is usedwithin the MiNLO approach [31] with the CT10 PDFs [32], interfaced to pythia8 with theAU2 tune, as a cross-check and to evaluate systematic uncertainties on the signal acceptance and kinematic properties. It is also used for the generation of gluon-gluon-initiatedZH production at leading order (LO) in QCD, with results cross-checked by an independent computation [33]. (For the analysis of the 7 TeV data, the pythia8 generator is used
for gg → ZH.) The transverse momentum distributions of the Higgs boson show substantial differences between the two ZH production processes. For qq-initiated W H and ZH
production, the total production cross sections and associated uncertainties are computed
at next-to-next-to-leading order (NNLO) in QCD [34–36], and with electroweak corrections at next-to-leading order (NLO) [37]. Additional normalisation-preserving differential
electroweak NLO corrections are applied as a function of the transverse momentum of thevector boson [38]. For gluon-gluon-initiated ZH production, NLO corrections [39], which
increase the total ZH production cross section by about 5%, are taken into account. The
Higgs boson decay branching ratios are calculated with hdecay [11]. Signal samples are
simulated for Higgs boson masses from 100 to 150 GeV in steps of 5 GeV. All charged-lepton
flavours are simulated in the W and Z decays, as leptonic decays of the τ leptons can also
be selected in the analysis. For the Higgs boson, only the bb decay mode is considered in
the analysis.
The main background processes are (W/Z)+jets and tt production. Version 1.4.1 ofthe sherpa generator [40] is used with the CT10 PDFs to simulate W +jets and Z+jets at
leading-order in QCD, with massive c- and b-quarks. For tt production, the simulation isProcess

Generator

Signal(?)
qq → ZH → ννbb/``bb

pythia8

gg → ZH → ννbb/``bb

powheg+pythia8

qq → W H → `νbb

pythia8

Vector boson + jets
Sherpa 1.4.1

Z/γ∗ → ``Sherpa 1.4.1

Z → νν

Sherpa 1.4.1

Top-quarktt̄

powheg+pythiat-channel

AcerMC+pythias-channel

powheg+pythia

Wtpowheg+pythia

Diboson(?)

powheg+pythia8WWpowheg+pythia8

WZ

powheg+pythia8

ZZ

powheg+pythia8

Table 1. The generators used for the simulation of the signal and background processes. (?) For
the analysis of the 7 TeV data, pythia8 is used for the simulation of the gg → ZH process, andherwig for the simulation of diboson processes.

performed with the powheg generator with the CT10 PDFs, interfaced with pythia6 [41],for which the CTEQ6L1 PDFs and the Perugia2011C tune [25, 26] are used. In thisanalysis, the final normalisations of these dominant backgrounds are constrained by thedata, but theoretical cross sections are used to optimise the selection. The cross sections
are calculated at NNLO for (W/Z)+jets [42] and at NNLO, including resummations of
next-to-next-to-leading logarithmic (NNLL) soft gluon terms, for tt [43].
Additional backgrounds arise from single-top-quark and diboson (W W , W Z, and ZZ)
production. For single-top-quark production, the s-channel exchange process and W t production are simulated with powheg, as for tt, while the t-channel exchange process issimulated with the AcerMC generator [44] interfaced with pythia6, using the CTEQ6L1
PDFs and the Perugia2011C tune. The cross sections are taken from refs. [45–47]. The
powheg generator with the CT10 PDFs, interfaced to pythia8 with the AU2 tune, is used
for diboson processes [48]. (For the analysis of the 7 TeV data, the herwig generator [49]is used instead with the CTEQ6L1 PDFs and the AUET2 tune [25, 26], and the crosssections are obtained at NLO from mcfm [50] with the MSTW2008NLO PDFs [51].)
Events from minimum-bias interactions are simulated with the pythia8 generator with
the MSTW2008LO PDFs [52] and the A2 tune [25, 26]. They are overlaid on the simulated–5–

JHEP01(2015)069

W → `ν

4

Object reconstruction

In this section, the reconstruction of physics objects used in the analysis of the 8 TeV
data is presented. Differences relevant for the analysis of the 7 TeV data are reported in
section 10.
Charged-particle tracks are reconstructed with a pT threshold of 400 MeV. The primary vertex is selected from amongst all reconstructed vertices as the one with the largest
sum of associated-track squared transverse momenta Σp2T and is required to have at least
three associated tracks.
Three categories of electrons [53, 54] and muons [55] are used in the analysis, referred toas loose, medium and tight leptons in order of increasing purity. Loose leptons are selectedwith transverse energy ET > 7 GeV. Loose electrons are required to have |η| < 2.47
and to fulfil the “very loose likelihood” identification criteria defined in ref. [54]. The
likelihood-based electron identification combines shower-shape information, track-quality
criteria, the matching quality between the track and its associated energy cluster in thecalorimeter (direction and momentum/energy), TRT information and a criterion to help
identify electrons originating from photon conversions. The electron energies are calibrated
by making use of reference processes such as Z → ee [56]. Three types of muons are includedin the loose definition to maximise the acceptance: (1) muons reconstructed in both themuon spectrometer and the inner detector (ID); (2) muons with pT > 20 GeV identified in
the calorimeter and associated with an ID track with |η| < 0.1, where there is limited muonchamber coverage; and (3) muons with |η| > 2.5 identified in the muon spectrometer, andwhich do not match full ID tracks due to the limited inner-detector coverage. For muons of–6–

JHEP01(2015)069

signal and background events according to the luminosity profile of the recorded data. The
contributions from these “pile-up” interactions are simulated both within the same bunch
crossing as the hard-scattering process and in neighbouring bunch crossings. The resulting
events are then processed through the same reconstruction programs as the data.
Additional generators are used for the assessment of systematic uncertainties as explained in section 8.
Simulated jets are labelled according to which generated hadrons with pT > 5 GeV
are found within a cone of size ∆R = 0.4 around the reconstructed jet axis. If a b-hadronis found, the jet is labelled as a b-jet. If not and a c-hadron is found, the jet is labelled
as a c-jet. If neither a b- nor a c-hadron is found, the jet is labelled as a light (i.e., u-,d-, or s-quark, or gluon) jet. Simulated V +jet events, where V stands for W or Z, arethen categorised according to the labels of the two jets that are used to reconstruct the
Higgs boson candidate. If one of those jets is labelled as a b-jet, the event belongs to theV b category. If not and one of the jets is labelled as a c-jet, the event belongs to theV c category. Otherwise, the event belongs to the V l category. Further subdivisions aredefined according to the flavour of the other jet from the pair, using the same precedence
order: V bb, V bc, V bl, V cc, V cl. The combination of V bb, V bc, V bl and V cc is denotedV +hf. The V cl final state is not included in V +hf because the main production process is
gs → W c rather than gluon splitting.

–7–

JHEP01(2015)069

the first and second type, the muon-track impact parameters with respect to the primaryvertex must be smaller than 0.1 mm and 10 mm in the transverse plane and along thez-axis, respectively. Finally, the scalar sum of the transverse momenta of tracks within acone of size ∆R = 0.2 centred on the lepton-candidate track, excluding the lepton track,
is required to be less than 10% of the transverse momentum of the lepton.
Medium leptons must meet the loose identification criteria and have ET > 25 GeV.
Medium muons must be reconstructed in both the muon spectrometer and the inner detector and have |η| < 2.5. Tight electrons are required to additionally fulfil the “very tightlikelihood” identification criteria [54]. For both the tight electrons and the tight muons,more stringent isolation criteria must be satisfied: the sum of the calorimeter energy deposits in a cone of size ∆R = 0.3 around the lepton, excluding energy associated with thelepton candidate, must be less than 4% of the lepton energy, and the track-based isolation
requirement is tightened from 10% to 4%.
Jets are reconstructed from noise-suppressed topological clusters of energy in the calorimeters [57] using the anti-kt algorithm [58] with a radius parameter of 0.4. Jet energies arecorrected for the contribution of pile-up interactions using a jet-area-based technique [59]and calibrated using pT - and η-dependent correction factors determined from simulation,with residual corrections from in situ measurements applied to data [60, 61]. Further
adjustments are made based on jet internal properties, which improve the energy resolution
without changing the average calibration (global sequential calibration [60]). To reduce thecontamination by jets from pile-up interactions, the scalar sum of the pT of tracks matchedto the jet and originating from the primary vertex must be at least 50% of the scalar sumof the pT of all tracks matched to the jet. This requirement is only applied to jets withpT < 50 GeV and |η| < 2.4. Jets without any matched track are retained. The jets kept
for the analysis must have pT > 20 GeV and |η| < 4.5.
To avoid double-counting, the following procedure is applied to loose leptons and jets.
First, if a jet and an electron are separated by ∆R < 0.4, the jet is discarded. Next, if
a jet and a muon are separated by ∆R < 0.4, the jet is discarded if it has three or fewer
matched tracks since in this case it is likely to originate from a muon having showered in
the calorimeter; otherwise the muon is discarded. (Such muons are nevertheless included
miss and in the jet energy corrections described in section 5.)
in the computation of the ET
Finally, if an electron and a muon are separated by ∆R < 0.2, the muon is kept unless it
is identified only in the calorimeter, in which case the electron is kept.
The MV1c b-tagging algorithm is used to identify jets originating from b-quark fragmentation. This algorithm combines in a neural network the information from various
algorithms based on track impact-parameter significance or explicit reconstruction of band c-hadron decay vertices. It is an improved version of the MV1 algorithm [62–64] withhigher c-jet rejection. Four b-tagging selection criteria (or operating points) are calibrated
and used in the analysis, corresponding to average efficiencies of 80%, 70%, 60% and 50%
for b-jets with pT > 20 GeV, as measured in simulated tt events. In this analysis, the 80%,
70% and 50% operating points are denoted loose, medium and tight, respectively. For thetight (loose) operating point, the rejection factors are 26 (3) and 1400 (30) against c-jetsand light jets, respectively. For the tight operating point, the c-jet rejection factor is 1.9
times larger than obtained with the MV1 algorithm.

5

Event selection

In this section, the event selection applied in the analysis of the 8 TeV data is presented.
Differences in the analysis of the 7 TeV data are reported in section 10.
The analysis is optimised for a Higgs boson mass of 125 GeV. Events are first categorised according to the numbers of leptons, jets, and b-tagged jets.

–8–

JHEP01(2015)069

The b-tagging efficiencies for b-jets, c-jets and light jets are measured in both data andsimulation using dedicated event samples such as tt events for b-jets, events with identifiedD∗ mesons for c-jets, or multijet events for light jets. The small differences observed areused to correct the simulation by so-called “scale factors” (SFs) within intervals between
two operating points. These SFs are parameterised as a function of the jet pT and, for light
jets, also |η|. The SFs are, however, strictly valid only for the generator used to derivethem. The differences observed when the efficiencies are measured with different generators
are taken into account by additional “MC-to-MC” SFs. Such differences can be caused by,e.g., different production fractions of heavy-flavour hadrons or modelling of their decays.
Because of the large cross sections of V l and V c production, these backgrounds remain
significant despite the powerful rejection of non-b-jets by the b-tagging algorithm. It is impractical to simulate a sufficiently large number of V l and V c events to provide a reliable
description of these backgrounds in the analysis samples for which two b-tagged jets are
required. An alternative procedure, parameterised tagging, is therefore used. Here, insteadof directly tagging the c- and l-labelled jets with the MV1c algorithm, parameterisations
as functions of pT and |η| of their probabilities to be b-tagged are used for the V l, V cand W W processes in all analysis samples in which two b-tagged jets are required. These
parameterisations are, however, integrated over other variables that can affect the c- andlight-jet tagging efficiencies. In particular, a strong dependence of these efficiencies is observed on ∆R, the angular separation from the closest other jet, and a significant differenceis seen between direct and parameterised tagging for V cc events with ∆R < 1. No suchdifference is seen for V cl, V l and W W events. A dedicated correction, depending on ∆R,
is therefore applied to the V cc events.
miss [65, 66] is measured as the negative
The missing transverse momentum vector ET
vector sum of the transverse momenta associated with energy clusters in the calorimeterswith |η| < 4.9. Corrections are applied to the energies of clusters associated with reconstructed objects (jets, electrons, τ leptons, and photons), using the calibrations of theseobjects. The transverse momenta of reconstructed muons are included, with the energy
deposited by these muons in the calorimeters properly removed to avoid double-counting.
In addition, a track-based missing transverse momentum vector, pmiss, is calculated as theTnegative vector sum of the transverse momenta of tracks with |η| < 2.4 associated with theprimary vertex.
Additional corrections are applied to the simulation to account for small differencesfrom data for trigger efficiencies, for lepton reconstruction and identification efficiencies, aswell as for lepton energy and momentum resolutions.

–9–

JHEP01(2015)069

Events containing no loose leptons are assigned to the 0-lepton channel. Events containing one tight lepton and no additional loose leptons are assigned to the 1-lepton channel.
Events containing one medium lepton and one additional loose lepton of the same flavour,and no other loose leptons, are assigned to the 2-lepton channel. In the 1- and 2-lepton
channels, for at least one of the lepton triggers by which the event was selected, the objects
that satisfied the trigger are required to be associated with the selected leptons.
The jets used in this analysis, called “selected jets”, must have pT > 20 GeV and|η| < 2.5, the η range within which b-tagging can be applied. There must be exactly twoor three such selected jets. Events containing an additional jet with pT > 30 GeV and|η| > 2.5 are discarded to reduce the tt background. Only selected jets are consideredfurther, e.g., to define the jet multiplicity, or to calculate kinematic variables.
The b-tagging algorithm is applied to all selected jets. There must be no more than two
such jets loosely b-tagged, and 3-jet events in which the lowest-pT jet is loosely b-tagged
are discarded. At least one of the two b-tagged jets must have pT > 45 GeV. The following
b-tagging categories are then defined as shown in figure 1. Events with two jets satisfying
the tight b-tagging criterion form the TT (or Tight) category; those not classified as TT,
but with two jets satisfying the medium b-tagging criterion, form the MM (or Medium)
category; those not classified as TT or MM, but with two jets satisfying the loose b-tagging
criterion, form the LL (or Loose) category. This categorisation improves the sensitivity
with respect to what would be obtained using a single category, such as TT+MM, with
the LL category providing constraints on the backgrounds not containing two real b-jets.
Events with exactly one jet loosely b-tagged form the 1-tag category, and those with noloosely b-tagged jet form the 0-tag category. In the 3-jet categories, the dijet system isformed by the two b-tagged jets in any of the 2-tag categories, by the b-tagged jet andthe leading (highest-pT ) non-b-tagged jet for events in the 1-tag category, and by the twoleading jets in the 0-tag category.
Additional topological and kinematic criteria are applied to reject background eventsand enhance the sensitivity of the search. They are outlined in table 2 and detailed below.
In general, the selection criteria are looser in the MVA than in the dijet-mass analysis inorder to maximise the information available to the final discriminant.
Further categorisation is performed according to the transverse momentum of thevector boson, pVT , to take advantage of the better signal-to-background ratio at high pVT .
miss in the 0-lepton
The transverse momentum of the vector boson is reconstructed as the ET
W
channel, the magnitude pT of the vector sum of the lepton transverse momentum and themiss in the 1-lepton channel, and the magnitude pZ of the vector sum of the transverseET
T
momenta of the two leptons in the 2-lepton channel. In the dijet-mass analysis, the eventsare categorised in five pVT intervals, with boundaries at 0, 90, 120, 160, and 200 GeV. In the
P jeti0-lepton channel and for events fulfilling the condition on
pT mentioned in table 2, the
miss
miss
miss = 120 GeV, andET trigger is fully efficient for ET > 160 GeV, 97% efficient for ET
miss = 100 GeV, with an efficiency that decreases rapidly for lower E miss .
80% efficient for ET
Tmiss value
Only four intervals are therefore used in the 0-lepton channel, with a minimum ET
miss trigger is used for pW > 120 GeV toof 100 GeV. In the 1-muon sub-channel, the ET
T
recover events not selected by the single-muon trigger, thus increasing the signal acceptancein this channel by 8%. In the MVA, only two intervals are defined, with pVT below or above
120 GeV, but the detailed pVT information is used in the final discriminant.
In the dijet-mass analysis, requirements are applied to the angular separation between
the two jets of the dijet system, ∆R(jet1 , jet2 ), which depend on the pVT interval. The
requirement on the minimum value reduces the background from V +jet production, while
the requirement on the maximum value, which reduces the background from tt production,is tightened with increasing pVT to take advantage of the increasing collimation of the dijet
system for the signal. To increase the signal acceptance, the requirement on the minimum
value is removed in the highest pVT interval, where the amount of background is smallest.
In the MVA, where the ∆R(jet1 , jet2 ) information is used in the final discriminant, only aminimum value is required, a requirement which is also removed for pVT > 200 GeV.
In the 0-lepton channel, the multijet (MJ) background is suppressed by imposing
requirements on the magnitude pmissof the track-based missing transverse momentumT
miss and pmiss , ∆φ(E miss , pmiss ), the azvector pmiss,theazimuthalangle
betweenET
T
T
T
T
miss and the nearest jet, min[∆φ(E miss , jet)], and the azimuthal
imuthal angle between ET
Tmiss and the dijet system, ∆φ(E miss , dijet). In addition, a minimum
angle between ET
T
P jetivalue is required for the scalar sum of the jet transverse momenta,
pT , which depends
on the jet multiplicity. Additional requirements are applied in the lowest pVT interval ofmiss > 100 GeV;
the 0-lepton channel, where the MJ background is largest: Njet = 2; ET∆φ(jet1 , jet2 ) < 2.7; S > 7; and L > 0.5. Here, ∆φ(jet1 , jet2 ) is the azimuthal angle be-– 10 –JHEP01(2015)069

Figure 1. Event classification as a function of the output of the MV1c b-tagging algorithm for thetwo highest pT jets. The bin boundaries denote the operating points (MV1c(jet) OP) as definedin section 4, corresponding to b-tagging efficiencies of 100%, 80%, 70%, 50%, i.e., the b-jet purity
increases from left (bottom) to right (top). The event categories are 0-tag, 1-tag, and TT, MM and
LL for 2-tag, as explained in the text.

Variable

Dijet-mass analysis

Multivariate analysis

Common selection
pVT [GeV]

0–9090(∗) –120

120–160 160–200 > 200 0–120∆R(jet1 , jet2 )0.7–3.4

0.7–3.0

0.7–2.3

0.7–1.8> 120< 1.4 > 0.7 (pVT <200 GeV)

0-lepton selection
pmiss
[GeV]
T

> 30

> 30

> 30

miss , pmiss )∆φ(ETTmiss , jet)]min[∆φ(ETmiss , dijet)
∆φ(ET
NjetP
=2(3)
jet
pT i [GeV]
i=1

< π/2

< π/2

< π/2

–

> 1.5

> 2.2

> 2.8

> 120 (NU)

> 120 (150)

> 1.5

NU

–
> 120 (150)

1-lepton selectionmW
T [GeV]< 120

HT [GeV]> 180

miss [GeV]
ET

–

–
–
> 20

> 50

> 180

–

–

> 20

2-lepton selectionm`` [GeV]

83-99

71-121

miss [GeV]
ET

< 60

–

Table 2. Event topological and kinematic selections. NU stands for ‘Not Used’. (∗) In the 0lepton channel, the lower edge of the second pVT interval is set at 100 GeV instead of 90 GeV. For
the 1-lepton channel, only the 1-muon sub-channel is used in the pVT < 120 GeV intervals.
miss significance, defined as the ratio of E miss to the square
tween the two jets, S is the ET
T
P jetiroot ofpT ; and L is a likelihood ratio constructed to discriminate further against the
MJ background.2
In the 1-lepton channel, a requirement is imposed on the transverse mass3 mW
T in thedijet-mass analysis. This requirement reduces the contamination from the tt background.
miss ) for pV < (>)120 GeV, where H is theRequirements are also imposed on HT (ETT
T
miss
scalar sum of ET and the transverse momenta of the two leading jets and the lepton.
This mainly reduces the MJ background. As discussed in section 7.1, the MJ backgroundis difficult to model and remains substantial in the 1-electron sub-channel in the pVT <120 GeV intervals. Therefore, only the 1-muon sub-channel is used in these intervals.
In the 2-lepton channel, criteria are imposed on the dilepton invariant mass, m`` , which2

miss
The likelihood ratio uses the following inputs: ∆φ(ET, dijet); ∆φ(jet1 , jet2 ); the magnitude of the
P jetimiss
missvector sum of the two jet transverse momenta, HT ; HT divided by
pT ; and the cosine of the helicity
angle in the dijet rest frame as defined in ref. [67]. For the MJ background, the probability density functions
missused in the likelihood ratio are constructed from data events selected with ∆φ(ET, pmiss) > π/2.
T
3
W
The transverse mass mT is calculated from the transverse momentum and the azimuthal angle of the`misscharged lepton, p`T and
p φ , and from the missing transverse momentum magnitude, ET , and azimuthal
miss
W
`
miss
`
miss
angle, φ: mT = 2pT ET (1 − cos(φ − φ)).

– 11 –JHEP01(2015)069

NU

– 12 –JHEP01(2015)069

must be consistent with the mass of the Z boson. In the dijet-mass analysis a requirement
miss ; this variable is used in the final discriminant of the MVA.
is imposed on ET
For events in which two jets are loosely b-tagged, these selection criteria define a set of“2-tag signal regions”, categorised in terms of channel (0, 1, or 2 leptons), pVT interval, andjet multiplicity (2 or 3). In the dijet-mass analysis, a further division is performed into theTT, MM and LL b-tagging categories. In the MVA, where the b-tagging information is usedin the final discriminant, a similar subdivision is performed with the difference that the
TT and MM categories are merged in the 0- and 2-lepton channels. Similarly defined 1-tag
and 0-tag “control regions” are used in the analysis to constrain the main backgrounds.
In the 1-lepton channel, the 2-tag signal regions with a third selected jet act in practice
as control regions because they are largely dominated by tt events. All 2-tag signal and1-tag control regions are used simultaneously in the global fit (described in section 9) used
to extract the results. The 0-tag control regions are used only for background modellingstudies (reported in section 7).
After event selection, the energy calibration of the b-tagged jets is improved as follows.
The energy from muons within a jet is added to the calorimeter-based jet energy after
removing the energy deposited by the muon in the calorimeter (muon-in-jet correction), anda pT -dependent correction is applied to account for biases in the response due to resolution
effects (resolution correction). This latter correction is determined for the pT spectrum ofjets from the decay of a Higgs boson with mH = 125 GeV in simulated (W/Z)H events.
The dijet mass resolution for the signal is improved by 14% after these corrections and is
miss involved
typically 11% (figure 2(a)). In the 2-lepton channel, wherein there is no true ET
except possibly from semileptonic heavy-flavour decays, the energy calibration of the jets
is further improved by a kinematic likelihood fit, which includes a Breit-Wigner constraint
on the dilepton mass, Gaussian constraints on each of the transverse components of the``bb system momentum (with a width of 9 GeV, as determined from ZH simulated events),
dedicated transfer functions relating the true jet transverse momenta to their reconstructedvalues (after the muon-in-jet correction, but without the resolution correction) as well asa prior built from the expected true jet pT spectrum in ZH events (playing a role similar
to the resolution correction). Overall, the bb mass resolution is improved by 30% in the
2-lepton channel (figure 2(b)).
The cross sections times branching ratios for (W/Z)H with W → `ν, Z → ``, Z → νν,and H → bb, as well as the acceptances in the three channels after full selection are given intable 3 for the MVA and the dijet-mass analysis. The acceptance for other production anddecay modes of the Higgs boson is negligible. The 0-lepton channel adds 7% in acceptancefor the W → `ν process with respect to the 1-lepton channel. Similarly, the 1-lepton channel
adds 10% in acceptance for the Z → `` process with respect to the 2-lepton channel.

Resolutions

0.08
0.06

2 lep., 2 jets, 2 b-tags
p V inclusiveT( σGSC-σ)/ σGSC

16.4 GeV--

14.4 GeV

12%

14.1 GeV

14%

Global Sequential Calib. (GSC)

Events / 4.0 GeV

Events / 4.0 GeVATLAS Simulation

0.1 Pythia VH, H → bb MC

ATLAS Simulation

0.04

0.08

--

30%

0.06

Global Sequential Calib. (GSC)+ Muon-in-Jet Correction+ Kinematic Likelihood Fit

0.02

20

40

60

80

100

120

140

160

180 200
mbb [GeV]

0
0

20

40

(a)

60

80

100

120

140

160

180 200
mbb [GeV](b)

Figure 2. Dijet-invariant-mass distribution for the decay products of a Higgs boson with mH =125 GeV in the 2-lepton MVA selection. The distributions are shown (a) using jets after global
sequential calibration (GSC, solid), and after adding muons inside jets (dotted) and after correcting
for resolution effects specific to the kinematics of the decay of a Higgs boson with mH = 125 GeV(dash-dotted); (b) using jets after global sequential calibration (GSC, solid), and after adding
muons inside jets and applying the kinematic fit (dash-dotted). The distributions are fit to the
Bukin function [68] and the parameter representing the width of the core of the distribution isshown in the figures, as well as the relative improvement in the resolution with respect to jets after
the global sequential calibration.

mH = 125 GeV at
Process
qq → (Z → ``)(H → bb)
gg → (Z → ``)(H → bb)
qq → (W → `ν)(H → bb)
qq → (Z → νν)(H → bb)
gg → (Z → νν)(H → bb)√

s = 8TeV

Cross section × BR [fb]
14.9
1.3
131.7
44.2
3.8

Acceptance [%]0-lepton 1-lepton
2-lepton–
1.3 (1.1) 13.4 (10.9)–
0.9 (0.7) 10.5 (8.1)
0.3 (0.3) 4.2 (3.7)–
4.0 (3.8)––
5.5 (5.0)––

Table 3. The cross section times branching ratio (BR) and acceptance for the three channels at
8 TeV. For ZH, the qq- and gg-initiated processes are shown separately. The branching ratios arecalculated considering only decays to muons and electrons for Z → ``, decays to all three lepton
flavours for W → `ν and decays to neutrinos for Z → νν. The acceptance is calculated as thefraction of events remaining in the combined 2-tag signal regions of the MVA (dijet-mass analysis)
after the full event selection.

– 13 –JHEP01(2015)069

0011.4 GeV(σGSC-σ)/ σGSC

T

0.04

0.02

16.4 GeV

2 lep., 2 jets, 2 b-tags
p V inclusive+ Muon-in-Jet Correction+ Resolution Correction

Resolutions

0.1 Pythia VH, H → bb MC

6Multivariate analysis

4

In this and all similar figures, all backgrounds are taken into account, but those contributing less than
1% are omitted from the legend.

– 14 –JHEP01(2015)069

Although the dijet mass is the kinematic variable that provides the best discrimination
between signal and backgrounds, the sensitivity of the search is improved by making use ofadditional kinematic, topological and b-tagging properties of the selected events in a multivariate analysis. The Boosted Decision Tree (BDT) technique [69, 70] is used, which, similarly to other multivariate methods, properly accounts for correlations between variables.
Dedicated BDTs are constructed, trained and evaluated in each of the 0-, 1- and 2lepton channels in the 2-tag regions (with the LL, MM and TT categories combined) andseparately for the events with two and three jets. In the 0-lepton channel, only events with
pVT > 120 GeV are used, whereas for the 1- and 2-lepton channels individual BDTs areused for pVT < 120 GeV and pVT > 120 GeV. Events in the electron and muon sub-channels
are combined since none of the variables used are lepton-flavour specific. In the 0-lepton
channel, the final results are obtained using the MVA for pVT > 120 GeV. For the small
100 < pVT < 120 GeV interval, which has reduced sensitivity, no dedicated BDT is trained
and only the dijet-mass distribution is used.
The BDTs are trained to separate the (V H, H → bb) signal from the sum of the
expected background processes. The input variables used to construct the BDTs are chosen
in order to maximise the separation, while avoiding the use of variables not improving theperformance significantly. Starting from the dijet mass, additional variables are tried one
at a time and the one yielding the best separation gain is kept. This procedure is repeateduntil adding more variables does not result in a significant performance gain. The final sets
of variables for the different channels are listed in table 4. The b-tagged jets belonging to thedijet system (with mass denoted mbb ) are labelled in decreasing pT as b1 and b2 , and theirseparation in pseudorapidity is |∆η(b1 , b2 )|. The b-tagging information is provided by theoutputs of the MV1c neural network, M V 1c(b1 ) and M V 1c(b2 ). The angular separation,in the transverse plane, of the vector boson and the dijet system of b-tagged jets and theirpseudorapidity separation are denoted ∆φ(V, bb) and |∆η(V, bb)|, respectively. In the 0lepton channel, HT is defined as the scalar sum of the transverse momenta of all jets andmiss . In the 1-lepton channel, the angle between the lepton and the closest b-tagged jetETin the transverse plane is denoted min[∆φ(`, b)]. The other variables were defined in theprevious sections. In 3-jet events, the third jet is labelled as jet3 and the mass of the 3-jet
system is denoted mbbj .
The input variables of the BDTs are compared between data and simulation, and good
agreement is found within the assessed uncertainties. Selected input-variable distributions
are shown in figure 3.4 In this figure, as for all figures in this section, the MJ backgroundis estimated as described in section 7.1, corrections to the simulation as explained in section 7.2 are applied, and background normalisations and shapes are adjusted by the global
fit of the MVA as outlined at the beginning of section 7 and presented in more detail in
section 9. A similarly good agreement is found for the correlations between pairs of input
variables, as can be seen in figure 4.

jet

pT 3
mbbj

0-Lepton
×
×
×
×
×
×
×

1-Lepton
×
×
×
×
×
×
×

2-Lepton
×
×
×
×
×
×
×
×
×

×
×
×
×
×

×
×
×

×
×
Only in 3-jet events
×
×
×
×
×
×

Table 4. Variables used in the multivariate analysis for the 0-, 1- and 2-lepton channels.

The Toolkit for Multivariate Data Analysis, TMVA [71], is used to train the BDTs.
The values for the training parameters are found by determining the configuration withthe best separation between signal and background in a coarsely binned multi-dimensional
training parameter space, followed by more finely grained one-dimensional scans of individual training parameters. In order to make use of the complete set of simulated MC events
for the BDT training and evaluation in an unbiased way, the MC events are split into two
samples of equal size, A and B. The performance of the BDTs trained on sample A (B)is evaluated with sample B (A) in order to avoid using identical events for both training
and evaluation of the same BDT. Half of the data are analysed with the BDTs trained on
sample A, and the other half with the BDTs trained on sample B. At the end, the output
distributions of the BDTs trained on samples A and B are merged for both the simulated
and data events.
The values of the BDT outputs do not have a well-defined interpretation. A dedicatedprocedure is applied to transform the BDT-output distributions to obtain a smoother distribution for the background processes and a finer binning in the regions with the largest
signal contribution, while at the same time preserving a sufficiently large number of background events in each bin. Starting from a very fine-binned histogram of the BDT-output
distribution, the procedure merges histogram bins, from high to low BDT-output values,
until a certain requirement, based on the fractions of signal and background events in themerged bin, is satisfied. To limit the number of bins and to reduce the impact of statistical fluctuations, a further condition is that the statistical uncertainty of the expected– 15 –JHEP01(2015)069

Variable
pVT
miss
ET
b1
pT
pbT2
mbb
∆R(b1 , b2 )|∆η(b1 , b2 )|∆φ(V, bb)|∆η(V, bb)|HTmin[∆φ(`, b)]mW
Tm``
M V 1c(b1 )M V 1c(b2 )Events / 25 GeV-1Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
W+cl
W+l
Z+hf
Z+cl
Z+l
Uncertainty
Pre-fit background
VH(bb)×10-1

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
W+cl
W+l
Z+hf
Uncertainty
Pre-fit background
VH(bb)×90(d)T

50 100 150 200 250 300 350 400 450 500
pV [GeV]50 100 150 200 250 300 350 400 450 500

T

s = 8 TeV ∫Ldt = 20.3 fb
1 lep., 2 jets, 2 tags
p V >120 GeV

ATLAS

(a)

50 100 150 200 250 300 350 400 450 500
mbb [GeV]

50 100 150 200 250 300 350 400 450 500

T

s = 8 TeV ∫Ldt = 20.3 fb
0 lep., 2 jets, 2 tags
p V >120 GeV

ATLAS
1600

1.5
1
0.5

50

100

150

200

250

300

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
W+cl
W+l
Z+hf
Z+cl
Z+l
Uncertainty
Pre-fit background
VH(bb)×50

T

100

100

150

150

(e)

200

200

-1

(b)

250

250

300

300

T

1

400

350
400
p (b ) [GeV]350

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Z+hf
Z+cl
Z+l
Uncertainty
Pre-fit background
VH(bb)×60

50 100 150 200 250 300 350 400 450 500
Emiss
[GeV]
T

s = 8 TeV ∫Ldt = 20.3 fb
2 lep., 2 jets, 2 tags
p V >120 GeV50-150 100 150 200 250 300 350 400 450 500

ATLAS

50

0.5
0

1

1.50

200

400

600

800

1000

1200

T

s = 8 TeV ∫Ldt = 20.3 fb
0 lep., 2 jets, 2 tags
p V >120 GeV

1400

ATLAS

1800

250

0.5
0

1

1.50

50

100

150

200

20
1.5
1
0.5
0
0

100

200

300

400

500

600

700

800

900

1

0.5

0.5

0.5

T

1.5

1.5

2

2

1.5

1.5

1

1

2

2

-1

2.5

2.5

2.5

2.5

(f)

-1

(c)

s = 8 TeV ∫Ldt = 20.3 fb
2 lep., 2 jets, 2 tags
p V >120 GeV

ATLAS

1

0.5

T

s = 8 TeV ∫Ldt = 20.3 fb
1 lep., 2 jets, 2 tags
p V >120 GeV

ATLAS

3

3

3

3

4

4

1

3.5

3.5

4

4

5

4.5 5
∆ η(V,bb)

4.5

5

4.5 5
∆ R(b ,b2)

4.5

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Z+hf
Z+cl
Z+l
Uncertainty
Pre-fit background
VH(bb)×60

3.5

3.5

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
W+cl
W+l
Z+hf
Uncertainty
Pre-fit background
VH(bb)×50

Figure 3. Examples of variables input to the BDT in the 2-jet 2-tag category (LL, MM and TT combined) for pVT > 120 GeV: (a) 0-lepton channel,b1
miss
dijet mass; (b) 0-lepton channel, ET; (c) 1-lepton channel, ∆R(b1 , b2 ); (d) 1-lepton channel, pW
T ; (e) 2-lepton channel, pT ; (f) 2-lepton channel,|∆η(V, bb)|. The distributions for the 2-lepton channel in (e) and (f) are shown after having applied the kinematic fit as described in section 5.
The background contributions after the global fit of the MVA are shown as filled histograms. The Higgs boson signal (mH = 125 GeV) is shown
as a filled histogram on top of the fitted backgrounds, as expected from the SM (indicated as µ = 1.0), and, unstacked as an unfilled histogram,scaled by the factor indicated in the legend. The dashed histogram shows the total background as expected from the pre-fit MC simulation. The
entries in overflow are included in the last bin. The size of the combined statistical and systematic uncertainty on the sum of the signal and fitted
background is indicated by the hatched band. The ratio of the data to the sum of the signal and fitted background is shown in the lower panel.

0
1.5
1
0.5
0

500

1000

1500

2000

2500

3000

3500

4000

4500

0.5
0

1

1.50

200

400

600

800

1000

Data/Pred

Events / 20 GeV

Data/Pred

Events / 20 GeV
Data/Pred
Events / 20 GeV
Data/Pred

Events / 0.2
Data/Pred
Events / 0.2
Data/Pred

JHEP01(2015)069– 16 –

0-lepton

∫

1

2

(a)

3

MC

5

∆R(b1,b2)

4

0

10

20

30

40

50

60

Eventsmbb [GeV]

0
0

1

2

3

4

5

0
0

100

200

300

400

-1

1.5

V

s = 8 TeV

1

50

100

150

T

(b)

2

(c)

200

0-lepton, 2 tags, 2 jets, p > 120 GeV∫Ldt = 20.3 fb

ATLAS

0.5

T

V

0-lepton, 2 tags, 2 jets, p > 120 GeV

s = 8 TeV

500

-1

∫Ldt = 20.3 fb

ATLAS

600

700

3

250

300

MC:
Data 2012:

2.5

MC:
Data 2012:

4

350

400

mean
mean

3.5

mean
mean

5

500

mbb [GeV]

450

RMS
RMS

4.5

∆R(b1,b2)

RMS
RMS

Figure 4. Correlations between the dijet mass and ∆R(b1 , b2 ) input variables to the BDT in the 2-jet 2-tag category (LL, MM and TT combined)of the 0-lepton channel for pVT > 120 GeV: (a) dijet mass versus ∆R(b1 , b2 ) for the total expected background (shadings indicating the numbers ofevents) and the data (open boxes with the box size being proportional to the number of events), (b) and (c) show the mean values and RMS of the
projections onto the ∆R(b1 , b2 ) and dijet-mass axes, respectively, for the total expected background after the global fit of the MVA and the data.

0
0

50

100

150

200

250

300

-1

Data 2012

Ldt = 20.3 fb

350 2 tags, 2 jets, pV > 120 GeVT400

s = 8 TeV

ATLAS
450

500
mbb [GeV]∆R(b1,b2)

JHEP01(2015)069– 17 –

– 18 –JHEP01(2015)069

total background contribution has to be smaller than 10% in each merged bin. The freeparameters of the transformation algorithm are optimised to maximise the expected signal
sensitivity. For simplicity, these transformed outputs, which are used for the analysis, arecalled “BDTV H discriminants” in the following. An optimisation of the number of bins andbin boundaries is also performed for the mbb distribution used in the dijet-mass analysis
in a similar way, where the free parameters of the transformation algorithm are optimised
separately for the different analysis regions. The effect of the transformation on the BDToutput and dijet-mass distributions can be seen in figure 5 for the 1-lepton channel and one
signal region. The transformation groups into few bins the mbb regions that are far fromthe signal on each of the low and high mass sides, while it expands the region close to thesignal mass, where the signal-to-background ratio is largest. The effect on the BDT output is similar, but simpler to visualise because the signal and the background accumulate
initially on the high and the low sides of the distribution, respectively.
Correlations between input variables and the BDTV H discriminant can provide information on the impact of individual variables on the classification. Figure 6 shows suchcorrelations for the dijet mass, which is the BDT input that provides the best single-variable
discriminating power.

50

Events / 0.13

Events / 25 GeV

60

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
Z+hf
Uncertainty
Pre-fit background
VH(bb)×10

ATLAS
s = 8 TeV ∫Ldt = 20.3 fb
1 lep., 2 jets, 2 Tight tags
160<p V <200 GeV
-1

T

40

70

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
Z+hf
Uncertainty
Pre-fit background
VH(bb)×10

ATLAS
s = 8 TeV ∫Ldt = 20.3 fb
1 lep., 2 jets, 2 Tight tags
160<p V <200 GeV
-1

60

T

50
40

30
30
20

20

20 40 60 80 100 120 140 160 180 200 220

20 40 60 80 100 120 140 160 180 200 220
mbb [GeV]

Data/Pred

1.5
1
0.5

10
20
1.5
1
0.5
0
0

0.1 0.2 0.3

0.4 0.5 0.6

0.7 0.8 0.9

0.1 0.2 0.3

0.4 0.5 0.6

0.7 0.8 0.9
1Transformed mbb

220
200
180

(b)
Events / 0.14

Events / 0.08(a)Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
Z+hf
Uncertainty
Pre-fit background
VH(bb)×20

ATLAS
s = 8 TeV ∫Ldt = 20.3 fb
1 lep., 2 jets, 2 Tight tags
p V >120 GeV
-1

T

160
140

450

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
Z+hf
Uncertainty
Pre-fit background
VH(bb)×20

ATLAS
s = 8 TeV ∫Ldt = 20.3 fb
1 lep., 2 jets, 2 Tight tags
p V >120 GeV
-1

400
350

T

300

120

250

100

200

80

150

60

100

40

50

2-1 -0.8 -0.6 -0.4 -0.2
1.5
1
0.5
0
-1 -0.8 -0.6 -0.4 -0.2

0

0.2

0.4 0.6 0.8

10

0.2 0.4 0.6 0.8
1
Untransformed BDTVH

Data/Pred20
Data/Pred1(c)2-1 -0.8 -0.6 -0.4 -0.2
1.5
1
0.5
0
-1 -0.8 -0.6 -0.4 -0.2

0

0.2

0.4 0.6 0.8

1

0

0.2

0.4 0.6 0.8
1
BDTVH

(d)

Figure 5. Top: the dijet-mass distributions for the expected background and signal contributionsin the 1-lepton channel and the 2-jet 2-tag TT category for 160 GeV < pW
T ≤ 200 GeV (a) before and(b) after applying the transformation of the histogram bins. Bottom: the BDT-output distribution
for the expected background and signal contributions in the 1-lepton channel and the 2-jet 2-tag TT
category for pW
T > 120 GeV (c) before and (d) after applying the transformation of the histogram
bins. The background contributions after the relevant global fit (of the dijet-mass analysis in (a)and (b) and of the MVA in (c) and (d)) are shown as filled histograms. The Higgs boson signal(mH = 125 GeV) is shown as a filled histogram on top of the fitted backgrounds, as expectedfrom the SM (indicated as µ = 1.0), and, unstacked as an unfilled histogram, scaled by the factor
indicated in the legend. The dashed histogram shows the total background as expected from thepre-fit MC simulation. The entries in overflow are included in the last bin. The size of the combined
statistical and systematic uncertainty on the sum of the signal and fitted background is indicatedby the hatched band. The ratio of the data to the sum of the signal and fitted background is shownin the lower panel.

– 19 –JHEP01(2015)069

Data/Pred10

T∫Ldt = 20.3 fb

MC

Data 2012

2 tags, 2 jets, pV > 120 GeV

0-leptons = 8 TeV

ATLAS
-1

0

2

4

6

8

10

12

14

16

18

20

22

24

Eventsmbb [GeV]

-1
0

-0.5

0

0.5

1

0
-1

100

200

300

400

500
-1

s = 8 TeV

-1

-0.4

50

100

150

T

V

s = 8 TeV

-0.6

(c)

200

(b)

-0.2

0-lepton, 2 tags, 2 jets, p > 120 GeV∫Ldt = 20.3 fb

ATLAS

-0.8

T

0-lepton, 2 tags, 2 jets, pV > 120 GeV∫Ldt = 20.3 fb

ATLAS

0.2

250

300

MC:
Data 2012:

0MC:
Data 2012:

0.6

350

400

mean
mean

0.4

mean
mean

1

500

mbb [GeV]

450

RMS
RMS

BDTVH

0.8

RMS
RMS

Figure 6. Correlations between the dijet mass and the BDTV H discriminant in the 2-jet 2-tag category (LL, MM and TT combined) of the
0-lepton channel for pVT > 120 GeV: (a) the dijet mass versus the BDTV H discriminant for the total expected background (shadings indicating thenumbers of events) and the data (open boxes with the box size being proportional to the number of events), (b) and (c) show the mean values and
RMS of the projections onto the axes of the BDTV H discriminant and dijet mass, respectively, for the total expected background after the global
fit of the MVA and the data.

(a)

BDTVH

0
-1 -0.8-0.6 -0.4-0.2 0 0.2 0.4 0.6 0.8 1

50

100

150

200

250

300

350

400

450

500
mbb [GeV]BDTVH

JHEP01(2015)069– 20 –

7

Background composition and modelling

7.1

Multijet background

Multijet events are produced with a huge cross section via the strong interaction, and therefore give rise to potentially large backgrounds. A first class of MJ background arises fromjets or photon conversions misidentified as electrons, or from semileptonic heavy-flavour
decays; the 1- and 2-lepton channels are especially sensitive to this class of background.
Another class, which affects mostly the 0-lepton channel, arises from large fluctuations in
miss . These MJ backjet energy measurements in the calorimeters, which create “fake” ETgrounds cannot be determined reliably by simulation, and are estimated from data in eachof the 0-, 1-, and 2-lepton channels, and in each of the 2- and 3-jet, 0-, 1-, and 2-tag regions.
The MJ background is estimated in the 0-lepton channel using an “ABCD method”,
miss , jet)]within which the data are divided into four regions based on the min[∆φ(ETmiss , pmiss ) variables, such that three of the regions are dominated by backand ∆φ(ETTground. (In the 100–120 GeV pVT interval, the likelihood ratio L designed to suppress themiss , jet)].) For events with real E miss , it
MJ background is used instead of min[∆φ(ETT
is expected that the directions of the calorimeter-based and track-based missing transmiss and pmiss , are similar. In events with fake E miss arising fromverse momenta, ETTTmiss is close to the dia jet energy fluctuation, it is expected that the direction of ET
rection of the poorly measured jet. The signal region (A) is therefore selected with
miss , jet)] > 1.5 and ∆φ(E miss , pmiss ) < π/2. In region C, the requirement onmin[∆φ(ETT
Tmiss , pmiss ) is reversed. In regions B and D, min[∆φ(E miss , jet)] < 0.4 is required,
∆φ(ET
T
T
miss , pmiss ) as in regions A and C, respectively. A comparison
with requirements on ∆φ(ETT
miss , jet)] distributions for ∆φ(E miss , pmiss ) above and below π/2 showsof the min[∆φ(ETT
Tthat these two variables are only weakly correlated, and this observation is confirmed in a

– 21 –JHEP01(2015)069

This section describes the modelling of individual backgrounds. In many cases, the data
are able to constrain the normalisations and shapes better than the a priori estimates. A
likelihood fit (also called “global fit”) is used to simultaneously extract both the signal yieldand constraints on the background normalisations and shapes. The distributions used by
the fit are those of the dijet mass or BDTV H discriminant in the 2-tag signal regions, as
appropriate, as well as those of the M V 1c value of the b-tagged jet in the 1-tag control
regions. More details are provided in section 9.
For the multijet (MJ) backgrounds, the normalisations and shapes provided as inputs
to the fit are estimated from data, as explained below. For the other backgrounds the
inputs are taken from the simulation, except for the normalisations of the V +jets and tt
backgrounds that are left free to float in the fit. The corrections to these two backgrounds,
described below, are applied prior to the fit.
In all distributions presented in this section, unless otherwise specified, the normalisations of the various backgrounds are those extracted from the global fit for the dijet-mass
or multivariate analysis, as appropriate. The fit also adjusts the background shapes inthose distributions within the constraints from the systematic uncertainties discussed in
section 8.

– 22 –JHEP01(2015)069

multijet event sample simulated with pythia8. An MJ template in region A is obtainedusing events in region C after subtracting the contribution of other backgrounds, takenfrom simulation. The template is normalised by the ratio of the number of events in region
B to that in region D, again after subtracting other backgrounds from those regions. The
populations of events in the various regions suffer from low statistical precision after the
2-tag requirement. The b-tagging requirement is therefore dropped in regions B, C andD, and an additional b-tagging normalisation factor is applied to the resulting template,taken as the fraction of 2-tag events in region D. The MJ background in the signal regions
is found to amount to ∼ 1% of the total background.
In the 1-lepton channel, the MJ background is determined separately for the electron
and muon sub-channels. For each signal or control region, an MJ-background template is
obtained in an MJ-dominated region after subtracting the small remaining contribution
from the other backgrounds. The other backgrounds are taken from a simulation improvedby scale factors for the various contributions obtained from a preliminary global fit. The
MJ-dominated region is obtained by modifying the nominal selection to use medium, instead of tight, leptons and loosening both the track and calorimeter-based isolation criteria.
The track-based isolation is changed to the intervals 5%–12% and 7%–50% for electrons
and muons respectively, instead of < 4%; and the calorimeter-based isolation is loosened
to < 7% from < 4%. The sample sizes of the MJ-templates are however rather low in the
2-tag regions. Since it is observed that the kinematic properties of the 1-tag and 2-tag
events in the MJ-dominated regions are similar, 1-tag events are used to enrich the 2-tag
MJ templates. Events in the 1-tag category are promoted to the 2-tag category by assigning to the untagged jet an emulated M V 1c value drawn from the appropriate M V 1c
distribution observed in the corresponding 2-tag MJ template. This distribution depends
on the rank (leading or sub-leading) of the untagged jet and on the M V 1c value of thetagged jet. To cope with residual differences observed in some distributions between thesepseudo-2-tag MJ events and the actual 2-tag MJ events, a reweighting is applied according
to the M V 1c of the tagged jet and, for the electron sub-channel, according to ∆R(jet1 , jet2 )and pW
T . This procedure is applied in each of the 2- and 3-jet, LL, MM and TT categories.
miss
The normalisations of the MJ templates are then obtained from “multijet fits” to the ET
distributions in the 2- and 3-jet, 1- and 2-tag (LL, MM and TT combined) categories,
with floating normalisations for the templates of the other background processes. The
templates for these other background processes are taken from the improved simulation
mentioned above.
The MJ background in the 1-lepton channel is concentrated at low pW
T , and in theW
2-jet 2-tag sample with pT < 120 GeV it ranges from 11% of the total background in the
LL category to 6% in the TT category. The main purpose of including the pW
T < 120 GeV
intervals is to provide constraints on the largest backgrounds (V +jets and tt) in the global
fit. Since the MJ background is twice as large for pW
T < 120 GeV in the 1-electron subchannel than in the 1-muon sub-channel, only the 1-muon sub-channel is kept for pW
T <120 GeV so as to provide the most reliable constraints on the non-MJ backgrounds. The
resulting loss in sensitivity is 0.6%. For pW
T > 120 GeV, the MJ background is muchsmaller: 4% and 2% in the LL and TT categories, respectively, for 2-jet events.

7.2

Corrections to the simulation

The large number of events in the 0-tag samples allows for detailed investigations of themodelling of the V +jet backgrounds by the version of the sherpa generator used in thisanalysis. Given that the search is performed in intervals of pVT , with the higher intervals
providing most of the sensitivity, an accurate modelling of the pVT distribution is important.
Figure 7(a) shows that the pW
T spectrum for W +jets production in the 1-muon subchannel is softer in the data than in the simulation. It is found that this mismodellingis strongly correlated with a mismodelling of the ∆φ(jet1 , jet2 ) distribution,5 shown infigure 8(a).6 In order to address this mismodelling, the W l and W cl simulations are
reweighted based on parameterised fits to the ratio of data to simulation in the ∆φ(jet1 , jet2 )variable in the 0-tag region, where these backgrounds dominate. Four separate functions
are derived: for the 2- and 3-jet categories and for pW
T above and below 120 GeV. The
reweighted ∆φ(jet1 , jet2 ) distributions show good agreement between data and simulation(figure 8(b)). This reweighting increases (reduces) by 0.7% (5.6%) the normalisation of theWpW
T < (>) 120 GeV region. After this reweighting, the modelling of the whole pT distribution is greatly improved, as can be seen in figure 7(b). This reweighting also improves
the modelling of other distributions, most notably the dijet mass. It also improves themodelling in the 1-tag control regions and is therefore applied to the W l and W cl backgrounds in all regions of all channels. The numbers of W cc and W b background events in
the 0- and 1-tag regions are too small to allow conclusive studies of their modelling, so no
reweighting is applied to these backgrounds, but an associated systematic uncertainty isassessed instead, as explained in section 8.
A similar, but not identical, procedure is used for the Z+jet events in the 2-lepton
channel. A ∆φ(jet1 , jet2 ) reweighting is found to improve the modelling of the pZ
T distribution in the 0-tag regions. In the signal-depleted 2-tag regions obtained by the exclusion
5

It has indeed been observed that the shape of the ∆φ(jet1 , jet2 ) distribution in data is better reproduced
by NLO generators than by the baseline sherpa generator used in this analysis [72].
6
The peak around ∆φ(jet1 , jet2 )=0.7 comes from the combination of two effects: a rise towards low∆φ(jet1 , jet2 ) due to gluon splitting, and a drop towards low ∆φ(jet1 , jet2 ) due to the two jets becoming
unresolved.

– 23 –JHEP01(2015)069

A template for the MJ background in the 2-electron sub-channel is obtained in a similarway, by loosening identification and isolation requirements. The normalisation is performed
by a fit to the dilepton-mass distribution, where the Z+jets and MJ components are freeparameters, while the other backgrounds (mostly tt) are taken from the simulation. The
MJ normalisation factors are found to be consistent in the 0-, 1- and 2-tag regions. To copewith the reduced size of the 2-tag MJ event sample, a procedure similar to that used in the
1-lepton channel is used, wherein the pretag MJ sample is weighted by its 2-tag fraction
and combinations of M V 1c values are randomly assigned to the jets according to theirdistribution in the 2-tag MJ template. In the 2-muon sub-channel, the MJ backgroundis found to be negligible from a comparison between data and MC prediction in the sidebands of the Z mass peak. Altogether, the MJ background amounts to <1% of the total
background in the 2-lepton channel.

Data 2012
VH(bb) (µ =1.0)∫ L dt = 20.3 fb

Diboson
tt

-1

60000

Events / 5 GeV

ATLAS
s = 8 TeV

70000

Single top

50000

Multijet
W+hf

80000

ATLAS
s = 8 TeV

70000

Data 2012
VH(bb) (µ =1.0)∫ L dt = 20.3 fb

Diboson
tt

-1

60000

Single top

50000

Multijet
W+hf

40000

W+cl
W+l

40000

W+cl
W+l

30000

Z+hf
Z+cl

30000

Z+hf
Z+cl

Z+l

20000

10000

10000
50

100

150

200

250

Data / Pred

0

1.1

1
0.9
0

50

100

150

Z+l

0

1.1

50

100

150

200

250

50

100

150

200
250
pW [GeV]

1
0.9

200
250
pW [GeV]

0

T

T

(a)

(b)

25000
20000
15000
10000

Events / πx10-2 rad

Events / πx10-2 rad

Figure 7. The pW
T distribution observed in data (points with error bars) and expected (histograms)
for the 2-jet 0-tag control region of the 1-muon sub-channel (MVA selection), (a) before and (b) after∆φ(jet1 , jet2 ) reweighting. The multijet and simulated-background normalisations are provided by
the multijet fits. The size of the statistical uncertainty is indicated by the shaded band. The
data-to-background ratio is shown in the lower panel.

ATLAS
s = 8 TeV

∫ L dt = 20.3 fb-1

VH(bb) (µ=1.0)
tt
Multijet
W+cl
Z+hf
Z+l

Data 2012
DibosonSingle top
W+hf
W+l
Z+cl

ATLAS
s = 8 TeV

∫ L dt = 20.3 fb-1

VH(bb) (µ=1.0)
tt
Multijet
W+cl
Z+hf
Z+l

Data 2012
DibosonSingle top
W+hf
W+l
Z+cl

5000
0.5

1

1.5

2

2.5

3

1
0.9
0

15000

Data / Pred

Data / Pred

0

20000

10000

5000

1.1

25000

0.5

1

1.5

2

2.5
3
∆ φ(jet1,jet2)

0

1.1

1

1.5

2

2.5

3

0.5

1

1.5

2

2.5
3
∆ φ(jet1,jet2)

1
0.9
0

(a)

0.5

(b)

Figure 8. The ∆φ(jet1 , jet2 ) distribution observed in data (points with error bars) and expected(histograms) for the 2-jet 0-tag control region of the 1-muon sub-channel (MVA selection), (a) beforeand (b) after reweighting. All pW
T intervals are combined. The multijet and simulated-backgroundnormalisations are provided by the multijet fits. The size of the statistical uncertainty is indicatedby the shaded band. The data-to-background ratio is shown in the lower panel.

– 24 –JHEP01(2015)069

20000

Data / Pred

Events / 5 GeV

80000

of the 100–150 GeV dijet mass interval, there is no evidence of a need for a ∆φ(jet1 , jet2 )Z
correction, but the pZ
T distribution is mismodelled. A dedicated pT reweighting is therefore
determined in the 2-tag regions. Applying the ∆φ(jet1 , jet2 ) reweighting to the Zl component and the pZ
T reweighting to the Zc and Zb components leads to good modelling also in
the 1-tag regions. This procedure is therefore used in all regions of all channels.
It has been observed in an unfolded measurement of the pT distribution of top quarks
from pair production that the powheg generator interfaced to pythia predicts too harda spectrum [73]. A correction accounting for this discrepancy is therefore applied at thelevel of generated top quarks in the tt production process.
Distributions in the dijet-mass analysis

Distributions of pVT and dijet mass are shown in figure 9 and in figures 10 and 11, respectively, for a selection of 2-tag signal regions of the dijet-mass analysis. It can be seen thatthe background composition in the signal regions varies greatly from channel to channel,
with the pVT interval, with the jet multiplicity, and with the b-tagging category considered.
The signal-to-background ratio is larger in the 2-jet and tighter b-tagging categories, andlower in the 3-jet and loose b-tagging categories.
In the 2-lepton channel, the dominant background is always Zbb. There is also a significant contribution from tt in the lower pZ
T intervals, and the relative diboson contribution
increases with pZ.
TFor the 1-lepton channel and in the 2-jet samples the combination of W bb and tt
accounts for most of the background in the most sensitive MM and TT categories, with therelative contribution of W bb and dibosons being largest in the tighter b-tagging categories
and increasing with pW
T . The flavours of the two selected jets from tt depend on theWreconstructed pT interval. In particular, at high pW
T , when the b-quark and the W from a
top-quark decay are collimated, there is a large bc contribution, where the c-quark comesfrom the W → cs decay. A significant contribution from single-top-quark production
processes is also seen. In the 3-jet category, the tt contribution is in general dominant,but there are significant contributions from single-top-quark production (mostly in the W tchannel) and from W bb, the latter increasing with pW
T . A non-negligible contribution of
MJ background can be seen in the lowest pW
intervalsof the 2-jet category.
T
In the 0-lepton channel, the main backgrounds arise from Zbb and tt, but the W bb
background is also significant. The relative tt contribution is largest in the lowest pVTintervals, and larger in the 3-jet than in the 2-jet category.
The variations in the background composition between categories allow the global fit todisentangle the rates of the various background sources. The non-negligible contributions
from the V cl and, to a lesser extent, the V l backgrounds are constrained in the global fit
by the LL b-tagging categories, and also by the M V 1c distributions of the b-tagged jet in
the 1-tag control regions. The 0-tag control regions are not taken into account in the global
fit, but are mainly used to improve the modelling of the V +jets backgrounds, as explained
in section 7.2.

– 25 –JHEP01(2015)069

7.3

Events / bin

-1

50

50

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
W+cl
W+l
Z+hf
Z+cl
Uncertainty
Pre-fit background
VH(bb)×30

(a)

100

100

0.9
0

1

1.10

102

10

3

104

10

5

6

T

250
pV [GeV]

250

-1

0.8
0

50

50

(d)

100

100

1

1.20

102

3

10

104

105

150

150

s = 8 TeV ∫Ldt = 20.3 fb
1 lep., 3 jets, 2 Medium+Tight tags

ATLAS

200

150

10

200

150

s = 8 TeV ∫Ldt = 20.3 fb
0 lep., 2 jets, 2 Medium+Tight tags

6

10
-1

200

200

250

(b)

100

100

250
pV [GeV]T

Data 2012
VH(bb) (µ=1.0)
ttSingle top
Multijet
W+hf
W+cl
Uncertainty
Pre-fit background
VH(bb)×260

50

50

s = 8 TeV ∫Ldt = 20.3 fb
1 lep., 2 jets, 2 Loose tags

ATLAS

0
1.2
1
0.8
0

10

102

3

10

104

5

T

250
pV [GeV]

250

-1

50

50

(e)

100

100

6

0.8
0

1

1.20

102

3

10

104

5

10

10

150

150

s = 8 TeV ∫Ldt = 20.3 fb
2 lep., 2 jets, 2 Medium+Tight tags

ATLAS

200

150

10

200

150

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
W+cl
W+l
Z+hf
Uncertainty
Pre-fit background
VH(bb)×270
-1

200

200

250

(c)

100

100

250
pV [GeV]T

Data 2012
VH(bb) (µ=1.0)
Diboson
tt
Z+hf
Z+cl
Z+l
Uncertainty
Pre-fit background
VH(bb)×70

50

50

150

150

s = 8 TeV ∫Ldt = 20.3 fb
1 lep., 2 jets, 2 Medium+Tight tags

ATLAS

200

200

250

250
pV [GeV]T

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
W+cl
W+l
Z+hf
Uncertainty
Pre-fit background
VH(bb)×80

channel with the Medium and Tight b-tagging categories (also referred to as MM and TT in the text) combined, (b) the 2-jet signal regions of
the 1-lepton channel for the LL category, (c) the 2-jet signal regions of the 1-lepton channel with the MM and TT categories combined, (d) the
3-jet signal regions of the 1-lepton channel with the MM and TT categories combined, and (e) the 2-jet signal regions of the 2-lepton channel with
the MM and TT categories combined. The background contributions after the global fit of the dijet-mass analysis are shown as filled histograms.
The Higgs boson signal (mH = 125 GeV) is shown as a filled histogram on top of the fitted backgrounds, as expected from the SM (indicated asµ = 1.0), and, unstacked as an unfilled histogram, scaled by the factor indicated in the legend. The dashed histogram shows the total backgroundas expected from the pre-fit MC simulation. Overflow entries are included in the last bin. The size of the combined statistical and systematic
uncertainty on the sum of the signal and fitted background is indicated by the hatched band. The ratio of the data to the sum of the signal andfitted background is shown in the lower panel.

Figure 9. The pVT distribution observed in data (points with error bars) and expected (histograms) for (a) the 2-jet signal regions of the 0-lepton

0.8
0

1

1.20

200

400

600

800

1000

Data/Pred

ATLAS

Events / bin
Data/Pred

Events / bin
Data/Pred

1200

Events / bin
Data/Pred

Events / bin
Data/Pred

JHEP01(2015)069– 26 –

-1

T

s = 8 TeV ∫Ldt = 20.3 fb
0 lep., 2 jets, 2 Medium+Tight tags
100<p V <120 GeV

ATLAS

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
W+cl
W+l
Z+hf
Z+cl
Uncertainty
Pre-fit background
VH(bb)×20(a)

20 40 60 80 100 120 140 160 180 200 220
mbb [GeV]T-1

s = 8 TeV ∫Ldt = 20.3 fb
0 lep., 2 jets, 2 Medium+Tight tagsincl. p V >120 GeV

ATLAS

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
W+cl
Z+hf
Z+cl
Uncertainty
Pre-fit background
VH(bb)×10

0.5

1

(b)

20 40 60 80 100 120 140 160 180 200 220
mbb [GeV]

1.5 20 40 60 80 100 120 140 160 180 200 220

100

200

300

400

500
T

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
W+hf
W+cl
Z+hf
Z+cl
Uncertainty
Pre-fit background
VH(bb)×20

0.5

1

(c)

20 40 60 80 100 120 140 160 180 200 220
mbb [GeV]

1.5 20 40 60 80 100 120 140 160 180 200 220

20

40

60

80

100

120

140

160

s = 8 TeV ∫Ldt = 20.3 fb
0 lep., 3 jets, 2 Medium+Tight tagsincl. p V >120 GeV

180

-1

ATLAS

200

Figure 10. The dijet-mass distribution observed in data (points with error bars) and expected (histograms) in the 0-lepton channel with the Medium
and Tight b-tagging categories (also referred to as MM and TT in the text) combined for (a) the 2-jet signal region in the 100 < pVT < 120 GeV
interval, (b) the 2-jet signal regions in the three intervals with pVT > 120 GeV combined, and (c) the 3-jet signal regions in the three intervals with
pVT > 120 GeV combined. The background contributions after the global fit of the dijet-mass analysis are shown as filled histograms. The Higgs
boson signal (mH = 125 GeV) is shown as a filled histogram on top of the fitted backgrounds, as expected from the SM (indicated as µ = 1.0), and,
unstacked as an unfilled histogram, scaled by the factor indicated in the legend. The dashed histogram shows the total background as expected fromthe pre-fit MC simulation. The entries in overflow are included in the last bin. The size of the combined statistical and systematic uncertainty on
the sum of the signal and fitted background is indicated by the hatched band. The ratio of the data to the sum of the signal and fitted backgroundis shown in the lower panel.

0.51

1.5 20 40 60 80 100 120 140 160 180 200 220

50

100

150

200

250

300

Events / 25 GeV
Data/Pred

Events / 25 GeV

Data/Pred

Events / 25 GeV
Data/Pred

JHEP01(2015)069– 27 –

600
500

s = 8 TeV ∫Ldt = 20.3 fb
1 lep., 2 jets, 2 Medium+Tight tagsincl. p V >120 GeV
-1

T

700
ATLAS
600
500

s = 8 TeV ∫Ldt = 20.3 fb
1 lep., 3 jets, 2 Medium+Tight tagsincl. p V >120 GeV
-1

T

400

300

300

200

200

100

100

1.5 20 40 60 80 100 120 140 160 180 200 220

1.5 20 40 60 80 100 120 140 160 180 200 220

10.5

20 40 60 80 100 120 140 160 180 200 220
mbb [GeV]

1
0.5

20 40 60 80 100 120 140 160 180 200 220
mbb [GeV]80

s = 8 TeV ∫Ldt = 20.3 fb
2 lep., 2 jets, 2 Medium+Tight tagsincl. p V >120 GeV-1T

Data 2012
VH(bb) (µ=1.0)
Diboson
tt
Z+hf
Z+cl
Z+l
Uncertainty
Pre-fit background
VH(bb)×10

Events / 25 GeV

ATLAS
100

(b)

ATLAS
50

40

20

20

10

(c)

-1

T

30

2 20 40 60 80 100 120 140 160 180 200 220
1.5
1
0.5
0
20 40 60 80 100 120 140 160 180 200 220
mbb [GeV]

s = 8 TeV ∫Ldt = 20.3 fb
2 lep., 3 jets, 2 Medium+Tight tagsincl. p V >120 GeV

40

60

Data/Pred

Events / 25 GeV(a)Data/Pred

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
W+hf
W+cl
Uncertainty
Pre-fit background
VH(bb)×40

1.5
1
0.5

Data 2012
VH(bb) (µ=1.0)
Diboson
tt
Z+hf
Z+cl
Z+l
Uncertainty
Pre-fit background
VH(bb)×10

20 40 60 80 100 120 140 160 180 200 220

20 40 60 80 100 120 140 160 180 200 220
mbb [GeV](d)Figure 11. The dijet-mass distribution observed in data (points with error bars) and expected(histograms) with the Medium and Tight b-tagging categories (also referred to as MM and TT inthe text) combined and the three intervals with pVT > 120 GeV combined for (a) the 2-jet signal
regions of the 1-lepton channel, (b) the 3-jet signal regions of the 1-lepton channel, (c) the 2-jet
signal regions of the 2-lepton channel, and (d) the 3-jet signal regions of the 2-lepton channel.
The background contributions after the global fit of the dijet-mass analysis are shown as filled
histograms. The Higgs boson signal (mH = 125 GeV) is shown as a filled histogram on top of thefitted backgrounds, as expected from the SM (indicated as µ = 1.0), and, unstacked as an unfilled
histogram, scaled by the factor indicated in the legend. The dashed histogram shows the total
background as expected from the pre-fit MC simulation. The entries in overflow are included in thelast bin. The size of the combined statistical and systematic uncertainty on the sum of the signal
and fitted background is indicated by the hatched band. The ratio of the data to the sum of thesignal and fitted background is shown in the lower panel.

– 28 –JHEP01(2015)069

Data/Pred

400

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
W+cl
W+l
Z+hf
Uncertainty
Pre-fit background
VH(bb)×20

Events / 25 GeV

ATLAS

Data/Pred

Events / 25 GeV

700

7.4

Distributions in the multivariate analysis– 29 –JHEP01(2015)069

Distributions of the BDTV H discriminants of the MVA are shown in figures 12 to 14 for
2-tag signal regions in the 2- and 3-jet categories of the 0-, 1- and 2-lepton channels.
It can be seen that the backgrounds dominated by light jets and, to a lesser extent, cjets accumulate at lower values of the BDTV H discriminants, due to the inclusion of theMV1c information as inputs to the BDTs. The composition of the dominant backgrounds
accumulating at higher values of the BDTV H discriminant is similar to what was alreadyobserved in the 2-tag signal regions of the dijet-mass analysis, namely V bb and tt, howeverwith a larger contribution of the latter due to the looser requirement on ∆R(jet1 , jet2 ) in
the MVA selection.
Distributions of the output of the MV1c b-tagging algorithm are shown in figure 15
for the b-tagged jet in the 1-tag control regions of the MVA, in the 2-jet category andfor pVT > 120 GeV. In these distributions, the four bins correspond to the four b-tagging
operating points and are ordered from left to right in increasing b-jet purity. It can be seenthat these distributions, which are used in the global fit, provide strong constraints on theV c and V l backgrounds. As in the dijet-mass analysis, the 0-tag control regions are notused in the global fit.

s = 8 TeV ∫Ldt = 20.3 fb
0 lep., 2 jets, 2 Medium+Tight tags
p V >120 GeV
-1

104

T

103
102

Events / 0.22

Events / 0.1

ATLAS

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
W+cl
Z+hf
Z+cl
Z+l
Uncertainty
Pre-fit background
VH(bb)×10

ATLAS
4

10

s = 8 TeV ∫Ldt = 20.3 fb
0 lep., 3 jets, 2 Medium+Tight tags
p V >120 GeV
-1

T

103

102

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
W+hf
W+cl
Z+hf
Z+cl
Uncertainty
Pre-fit background
VH(bb)×20

10
10
1
1

2-1 -0.8 -0.6 -0.4 -0.2
1.5
1
0.5
0
-1 -0.8 -0.6 -0.4 -0.2

0

0.2

0.4 0.6 0.8

1

0

0.2

0.4 0.6 0.8
1
BDTVH

Data/Pred

Data/Pred

10-1

(a)

-1 -0.8 -0.6 -0.4 -0.2
1.5
1
0.5
-1 -0.8 -0.6 -0.4 -0.2

0

0.2

0.4 0.6 0.8

1

0

0.2

0.4 0.6 0.8
1
BDTVH

(b)

Figure 12. The BDTV H -discriminant distribution observed in data (points with error bars) and
expected (histograms) for the 0-lepton channel combining the 2-tag Medium and Tight b-tagging
categories (also referred to as MM and TT in the text) for pVT > 120 GeV for (a) 2-jet events and(b) 3-jet events. The background contributions after the global fit of the MVA are shown as filled
histograms. The Higgs boson signal (mH = 125 GeV) is shown as a filled histogram on top of
the fitted backgrounds, as expected from the SM (indicated as µ = 1.0), and, unstacked as anunfilled histogram, scaled by the factor indicated in the legend. The dashed histogram shows thetotal background as expected from the pre-fit MC simulation. The size of the combined statisticaland systematic uncertainty on the sum of the signal and fitted background is indicated by thehatched band. The ratio of the data to the sum of the signal and fitted background is shown in thelower panel.

– 30 –JHEP01(2015)069

510s = 8 TeV ∫Ldt = 20.3 fb
1 lep., 2 jets, 2 Medium+Tight tags
p V <120 GeV
-1

5

10

T

104

Events / 0.2

Events / 0.2

ATLAS

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
W+cl
W+l
Z+hf
Uncertainty
Pre-fit background
VH(bb)×60

510Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
W+cl
W+l
Z+hf
Uncertainty
Pre-fit background
VH(bb)×40

ATLAS
s = 8 TeV ∫Ldt = 20.3 fb
1 lep., 2 jets, 2 Medium tags
p V >120 GeV
-1

104

T

3

10

102

3

10

10

1.2-1 -0.8 -0.6 -0.4 -0.2

0

0.2

0.4 0.6 0.8

1

0

0.2

0.4 0.6 0.8
1
BDTVH

1
0.8
-1 -0.8 -0.6 -0.4 -0.2

Data/Pred

Data/Pred12-1 -0.8 -0.6 -0.4 -0.2
1.5
1
0.5
0
-1 -0.8 -0.6 -0.4 -0.2

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
Z+hf
Uncertainty
Pre-fit background
VH(bb)×20

ATLAS
104
3

10

0.2

0.4 0.6 0.8

1

0

0.2

0.4 0.6 0.8
1
BDTVH

(b)
Events / 0.25

Events / 0.14(a)0s = 8 TeV ∫Ldt = 20.3 fb
1 lep., 2 jets, 2 Tight tags
p V >120 GeV
-1

T

102

5

10

104

ATLAS
s = 8 TeV ∫Ldt = 20.3 fb
1 lep., 3 jets, 2 Medium+Tight tags
p V >120 GeV
-1

T

3

10

Data 2012
VH(bb) (µ=1.0)
ttSingle top
Multijet
W+hf
W+cl
Uncertainty
Pre-fit background
VH(bb)×80

102
10

1

1

10-1

10-1

2-1 -0.8 -0.6 -0.4 -0.2
1.5
1
0.5
0
-1 -0.8 -0.6 -0.4 -0.2

0

0.2

0.4 0.6 0.8

1

0

0.2

0.4 0.6 0.8
1
BDTVH

Data/Pred

Data/Pred10(c)-1 -0.8 -0.6 -0.4 -0.2
1.5
1
0.5
-1 -0.8 -0.6 -0.4 -0.2

0

0.2

0.4 0.6 0.8

1

0

0.2

0.4 0.6 0.8
1
BDTVH

(d)

Figure 13. The BDTV H -discriminant distribution observed in data (points with error bars) and
expected (histograms) for the 2-tag signal regions of the 1-lepton channel for (a) 2-jet events with
the Medium and Tight b-tagging categories (also referred to as MM and TT in the text) combinedW
and with pW
T ≤ 120 GeV, (b) MM 2-jet events with pT > 120 GeV, (c) TT 2-jet events withW
pT > 120 GeV, and (d) MM and TT combined 3-jet events with pW
T > 120 GeV. The backgroundcontributions after the global fit of the MVA are shown as filled histograms. The Higgs boson
signal (mH = 125 GeV) is shown as a filled histogram on top of the fitted backgrounds, as expectedfrom the SM (indicated as µ = 1.0), and, unstacked as an unfilled histogram, scaled by the factor
indicated in the legend. The dashed histogram shows the total background as expected from thepre-fit MC simulation. The size of the combined statistical and systematic uncertainty on the sumof the signal and fitted background is indicated by the hatched band. The ratio of the data to thesum of the signal and fitted background is shown in the lower panel.

– 31 –JHEP01(2015)069

102

-1

T

0.2

0

1

0.4 0.6 0.8
1
BDTVH

0.4 0.6 0.8

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
Z+hf
Z+cl
Uncertainty
Pre-fit background
VH(bb)×80
T

-1

(b)

0

0

0.2

0.2

s = 8 TeV ∫Ldt = 20.3 fb
2 lep., 2 jets, 2 Medium+Tight tags
p V >120 GeV

ATLAS

2-1 -0.8 -0.6 -0.4 -0.2
1.5
1
0.5
0
-1 -0.8 -0.6 -0.4 -0.2

10-1

1

10

102

10

3

104

1

0.4 0.6 0.8
1
BDTVH

0.4 0.6 0.8

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Z+hf
Z+cl
Z+l
Uncertainty
Pre-fit background
VH(bb)×20
T

-1

0.2

0

(c)

0.2

0

s = 8 TeV ∫Ldt = 20.3 fb
2 lep., 3 jets, 2 Medium+Tight tags
p V >120 GeV

ATLAS

2-1 -0.8 -0.6 -0.4 -0.2
1.5
1
0.5
0
-1 -0.8 -0.6 -0.4 -0.2

10-1

1

10

102

3

10

104

1

0.4 0.6 0.8
1
BDTVH

0.4 0.6 0.8

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Z+hf
Z+cl
Uncertainty
Pre-fit background
VH(bb)×20

Figure 14. The BDTV H -discriminant distribution observed in data (points with error bars) and expected (histograms) for the 2-lepton channel
combining the 2-tag Medium and Tight b-tagging categories (also referred to as MM and TT in the text) for (a) 2-jet events with pZ
T ≤ 120 GeV,Z(b) 2-jet events with pZ
>
120
GeV,
and
(c)
3-jet
events
with
p
>
120
GeV.
The
backgroundcontributions
after
the
globalfitofthe MVA are
T
Tshown as filled histograms. The Higgs boson signal (mH = 125 GeV) is shown as a filled histogram on top of the fitted backgrounds, as expectedfrom the SM (indicated as µ = 1.0), and, unstacked as an unfilled histogram, scaled by the factor indicated in the legend. The dashed histogramshows the total background as expected from the pre-fit MC simulation. The size of the combined statistical and systematic uncertainty on thesum of the signal and fitted background is indicated by the hatched band. The ratio of the data to the sum of the signal and fitted background isshown in the lower panel.

(a)

0.2

0

s = 8 TeV ∫Ldt = 20.3 fb
2 lep., 2 jets, 2 Medium+Tight tags
p V <120 GeV

ATLAS-1 -0.8 -0.6 -0.4 -0.2
1.2
1
0.8
-1 -0.8 -0.6 -0.4 -0.2

10

10

2

10

3

104

5

10

Events / 0.18
Data/Pred

Events / 0.22

Data/Pred

Events / 0.25
Data/Pred

JHEP01(2015)069– 32 –

70

60

(a)

-1

50
0
MV1c(b) OP

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
W+cl
W+l
Z+hf
Z+cl
Z+l
Uncertainty
Pre-fit background
VH(bb)×420

1.05
1
0.95
80

10

20

30

70

-1

60

(b)

50
40

s = 8 TeV ∫Ldt = 20.3 fb
1 lep., 2 jets, 1 tag
p V >120 GeV
T

ATLAS

×103
60

50
0
MV1c(b) OP

Data 2012
VH(bb) (µ=1.0)
Dibosontt
Single top
Multijet
W+hf
W+cl
W+l
Z+l
Uncertainty
Pre-fit background
VH(bb)×620

1.05
1
0.95
80

1000

2000

3000

4000

5000

70

-1

60

(c)

7000
T

s = 8 TeV ∫Ldt = 20.3 fb
2 lep., 2 jets, 1 tag
p V >120 GeV

6000

ATLAS

8000

50
0
MV1c(b) OP

Data 2012
VH(bb) (µ=1.0)
Diboson
tt
Z+hf
Z+cl
Z+l
Uncertainty
Pre-fit background
VH(bb)×320

Figure 15. Distribution of the output of the MV1c b-tagging algorithm observed in data (points with error bars) and expected (histograms) for
the 1-tag control regions of the MVA in the 2-jet category with pVT > 120 GeV (a) in the 0-lepton channel, (b) in the 1-lepton channel, and (c) in the
2-lepton channel. The left bin boundaries denote the operating points (MV1c(b) OP) as defined in section 4, corresponding to b-tagging efficienciesof 80%, 70%, 60%, 50%, i.e., the b-jet purity increases from left to right. The background contributions after the global fit of the MVA are shownas filled histograms. The Higgs boson signal (mH = 125 GeV) is shown as a filled histogram on top of the fitted backgrounds, as expected from theSM (indicated as µ = 1.0), and, unstacked as an unfilled histogram, scaled by the factor indicated in the legend. The dashed histogram shows thetotal background as expected from the pre-fit MC simulation. The size of the combined statistical and systematic uncertainty on the sum of thesignal and fitted background is indicated by the hatched band. The ratio of the data to the sum of the signal and fitted background is shown in
the lower panel.

Ts = 8 TeV ∫Ldt = 20.3 fb
0 lep., 2 jets, 1 tag
p V >120 GeV

ATLAS

×103

1.05
1
0.95
80

5

10

15

20

25

30

35
EventsData/Pred

EventsData/Pred

EventsData/Pred

JHEP01(2015)069– 33 –

8

Systematic uncertainties

The systematic uncertainties discussed in this section are: those of experimental origin;those related to the multijet background estimation; and those associated with the modelling of the simulated backgrounds and Higgs boson signal.
8.1

Experimental uncertainties– 34 –JHEP01(2015)069

All relevant experimental systematic uncertainties are considered, such as those affecting
the trigger selection, the object reconstruction and identification, and the object energy
and momentum calibrations and resolutions. The most relevant ones are discussed in
the following.
miss trigger, an efficiency correction is derived from W → µν+jets and Z →For the ET
miss of 100 GeV,µ+ µ− +jets events. This correction amounts to 4.5% for events with an ET
miss > 120 GeV. The associatedthe threshold required in the analysis, and is below 1% for ET
uncertainties arise from the statistical uncertainties of this method and differences observed
miss (and thus
in the two event classes. They are very small (below 1%) for the high ET
miss interval of the 0-lepton channelhigh pVT ) intervals, and reach about 3% for the low ET(100–120 GeV).
For electrons and muons, uncertainties associated with the corrections for the trigger,
reconstruction, identification and isolation efficiencies are taken into account. Uncertainties
on energy and resolution corrections of the leptons are also considered. The impact of these
uncertainties is very small, typically less than 1%.
Several sources contribute to the uncertainty of the jet energy scale (JES) [61] related
e.g. to uncertainties from in situ calibration analyses, pile-up-dependent corrections andthe flavour composition of jets in different event classes. After being decomposed into uncorrelated components, these are treated as independent sources in the analysis. The total
relative systematic uncertainties on the JES range from about 3% to 1% for central jetswith a pT of 20 GeV and 1 TeV, respectively. An additional specific uncertainty of about
1%–2% affects the energy calibration of b-jets. Small uncertainties on the corrections applied to improve the dijet-mass resolution are also included. Corrections and uncertainties
are also considered for the jet energy resolution (JER) [74], with a separate contribution
for b-jets. The total relative systematic uncertainty on the JER ranges from about 10% to20%, depending on the η range, for jets with pT = 20 GeV to less than 5% for jets withpT > 200 GeV.
miss , as are the much smaller uncerThe JES uncertainties are propagated to the ET
tainties related to the energy and momentum calibration of leptons. An uncertainty on
miss also comes from the uncertainties on the energy calibration (8%) and resolution
the ET(2.5%) of calorimeter energy clusters not associated with any reconstructed object [66].
The b-tagging efficiencies for the different jet flavours are measured in both data andsimulation using dedicated event samples [63, 64]. The b-tagging efficiencies for simulatedjets are corrected within intervals between operating points by MC-to-data SFs, which
depend on the jet kinematics. For b-jets, the precision is driven by an analysis of tt
events in final states containing two leptons. The MC-to-data SFs are close to unity,

8.2

Uncertainties on the multijet backgrounds

In the 0-lepton channel, the robustness of the MJ background estimation is assessed bymiss , jet)] values defining the B and D regions of the ABCD method,varying the min[∆φ(ETand by replacing the b-tagging fractions measured in region D by those measured in region
B. A systematic uncertainty of 100% is assessed for this small (∼ 1%) background, uncorrelated between 2- and 3-jet, 1- and 2-b-tag categories. The MJ background in the 2-lepton
channel is also at the per-cent level, and an uncertainty of 100% is assigned.
In the 1-lepton channel, normalisation uncertainties arise from the statistical uncertainties of the multijet fits and from uncertainties on the non-MJ background subtractions
performed to construct the MJ templates. Normalisation uncertainties are also assessedin the LL, MM and TT categories to cover differences between multijet fits performed
inclusively in the 2-tag regions and in the individual categories. In the 2-jet 2-tag regionof the electron sub-channel, the overall normalisation uncertainties amount to 11%, 14%
and 22% in the LL, MM and TT categories, respectively. In the muon sub-channel, thecorresponding uncertainties are about three times larger because of the smaller size of the
MJ-enriched samples.
In the 1-lepton channel, shape uncertainties are assessed in the various regions bycomparison of evaluations obtained using MJ-enriched samples defined by isolation requirements different from those applied in the nominal selections. In the electron sub-channel,
an alternative template is constructed with a track-based isolation in the 12% to 50% interval, and another alternative template with a calorimeter-based isolation in the 0% to
4% interval. In the muon sub-channel, the results obtained with the nominal MJ template
are compared with those obtained with tighter or looser isolation requirements, definedby track-based isolation intervals of 7%–9.5% and 9.5%–50%, respectively. Furthermore,– 35 –JHEP01(2015)069

with uncertainties at the level of 2–3% over most of the jet pT range, reaching 5% for
pT = 20 GeV and 8% above 200 GeV. The uncertainties, which depend on pT and on
the interval between operating points, are decomposed into uncorrelated components andthe ten most significant ones are kept in the analysis. It was checked that the neglectedcomponents have a negligible impact. The uncertainties on c-jets are decomposed into 15
components, and the uncertainties on light jets, to which the analysis is much less sensitive,
are decomposed into ten components, accounted for in pT and η ranges. For b- and c-jets
further uncertainties are added for the application of the additional MC-to-MC SFs toobtain generator specific MC-to-data SFs as explained in section 4. Half of the correction
is used as systematic uncertainty. As discussed in section 4, a correction to c-jets in the
V cc samples, for which parameterised tagging is used, is applied at low ∆R to the closest
jet. Half of this correction is assigned as a systematic uncertainty.
The uncertainty on the integrated luminosity is 2.8%. It is derived, following the same
methodology as that described in ref. [20], from a preliminary calibration of the luminosity
scale derived from beam-separation scans performed in November 2012. It is applied tothe signal and backgrounds estimates that are taken from simulation.
A 4% uncertainty on the average number of interactions per bunch crossing is taken
into account.

half of the ∆R(jet1 , jet2 ) and pW
T reweightings mentioned in section 7.1 for the electron
sub-channel are taken as systematic uncertainties.
8.3

Uncertainties on the modelling of the simulated backgrounds

Top-quark-pair background. As explained in section 7, the top-quark pT distribution
is reweighted at generator level to bring it into agreement with measurement [73]. A systematic uncertainty amounting to half of this correction is assigned, correlated across channels.
The predictions of the nominal tt generator (powheg+pythia) are compared, focussing on the 1-lepton channel selection, with those obtained using a variety of generators differing by the PDF choice (powheg+pythia with HERAPDF [75]), by the parton
showering and hadronisation scheme (powheg+herwig), by the implementation of the
NLO matrix element and the matching scheme (mc@nlo [76]+herwig), by the amountof initial- and final-state radiation (ISR/FSR) using AcerMC+pythia, or by the implementation of higher-order tree-level matrix elements (alpgen [77]+pythia). It is foundthat, in general, the largest deviations are observed for alpgen, which is therefore used toassess further systematic uncertainties as explained below.
In the global fit, the normalisation of the tt background in the 2-jet category is left
floating freely, independently in each of the lepton channels. An uncertainty of 20% on
the 3-to-2-jet ratio is estimated from the generator comparisons explained above. In the– 36 –JHEP01(2015)069

The physics-modelling systematic uncertainties evaluated focus on the quantities that areused in the global fit, i.e., those affecting the jet multiplicities, the pVT distributions, theflavour composition and the mbb distributions. For the MVA, systematic uncertainties
affecting the other variables used as inputs to the BDTs are also considered. Wheneverpossible, dedicated control regions are used to extract information directly from the data.
This is the case for Z+jets and W +light jets. In other cases, uncertainties are assessed bycomparison of MC predictions based on a variety of generators with the nominal ones.
Details of the assessment of systematic uncertainties are provided below in the context of the MVA. When systematic uncertainties are derived from a comparison between
generators, all relevant variables are considered independently. The variable showing thelargest discrepancy in some generator with respect to the nominal generator is assigned an
uncertainty covering this discrepancy, which is symmetrised. If, once propagated to the
BDTV H discriminant, this uncertainty is sufficient to cover all variations observed with thedifferent generators, it is considered to be sufficient. If not, an uncertainty is consideredin addition on the next most discrepant variable and the procedure is iterated until all
variations of the BDTV H discriminant are covered by the assigned uncertainties.
A given source of systematic uncertainty can affect different analysis regions. Whether
such an uncertainty should be treated as correlated or not depends on whether constraints
resulting from the global fit should be propagated from one region to another. Details of
the procedures leading to such decisions are provided in section 9.2.
A summary of the systematic uncertainties affecting the modelling of the backgrounds
can be found in table 5.

global fit, this uncertainty is treated as correlated between the 0- and 1-lepton channels,
and uncorrelated with the 2-lepton channel.
The shape of the mbb distribution is also studied with the same set of generators,
leading to correlated shape uncertainties for 2- and 3-jet events, and for pVT < 120 GeV
and pVT > 120 GeV. The associated variation is larger in the higher pVT interval: for 2-jet
events, when it increases the distribution by 3% for mbb = 50 GeV, it decreases it by 1%
at 200 GeV; the effect is similar, but of opposite sign, for 3-jet events.

Single-top-quark background. The theoretical uncertainties on the cross sections of
the three processes contributing to single-top production are 4%, 4%, and 7% for thes-channel, t-channel, and W t production, respectively [78].
The predictions of the nominal generators (powheg+pythia for the s-channel andfor W t production; acerMC+pythia for the t-channel) are compared, after the 1-lepton
channel selection, with those obtained using a variety of generators. For the s-channel,
the comparison is made with acerMC and mc@nlo; for W t production with acerMC,
powheg+herwig, and mc@nlo; and for the t-channel with amc@nlo7 [81, 82]+herwig.
For all three processes, the impact of ISR/FSR is evaluated using acerMC. For W t production, there are interference effects with tt production, which need to be considered.
Two methods are available for this: the Diagram Removal (DR) and the Diagram Subtraction (DS) schemes [83]. The former is used in the nominal generation, and the second forcomparison.
Uncertainties on the acceptance for each of the three processes are taken as the largest
deviations observed, separately for pVT < 120 GeV and pVT > 120 GeV, and for 2- and 3-jet
events. They can be as large as 52% for 2-jet events in the t-channel at low pVT , of the orderof 5% for W t production (except for 3-jet events at high pVT : 15%), and typically 20% for
the s-channel.
In addition to the acceptance uncertainties, the effects of the model variations describedabove on variables input to the BDT are evaluated and three shape systematic uncertainties
are found to be needed in W t production. The first uncertainty is on the shape of the mbb
distribution in the high pVT interval for 2-jet events where, when a shift from the nominalmodel increases the rate by 20% for mbb = 50 GeV, it decreases it by 40% at 200 GeV. A
second uncertainty is on the mbb shape for 3-jet events, where the corresponding shifts are
25% and 20%. Finally, a third uncertainty is on the pT distribution of the leading jet in
the low pVT interval for 2-jet events.
7

Event generation with amc@nlo is based on the mc@nlo formalism and the madgraph-5 framework [79, 80].

– 37 –JHEP01(2015)069

The same procedure is used for the pVT distribution, from which a 7.5% uncertainty isassessed on the normalisation of the pVT > 120 GeV interval. Finally, the same approach
miss distribution in the 1-lepton channel,
calls for a shape systematic uncertainty on the ET
different but correlated between pVT < 120 GeV and pVT > 120 GeV. This uncertainty is notapplied in the 0- and 2-lepton channels.

W +jets background. As explained in section 7, a ∆φ(jet1 , jet2 ) reweighting is applied
to the W l and W cl components. Uncorrelated systematic uncertainties amounting to half ofthe correction are assigned to these two components, for each of the 2- and 3-jet categories.
For the W cc + W b component, no reweighting is applied but a systematic uncertainty is
assigned, equal to the full correction applied to the W l and W cl components, uncorrelated
between 2- and 3-jet events.
The normalisation and the 3-to-2-jet ratio for the W l background are taken directlyfrom simulation, both with a 10% uncertainty. This is based on the agreement observed
between data and prediction in the 0-tag sample. The 3-to-2-jet ratio for the W cl background is also assigned an uncertainty of 10%. The normalisations of the W cl and W bb
backgrounds are left free in the global fit.
To assign further uncertainties on the W bb background, for which dedicated control
regions are not available in the data, extensive comparisons are performed at generator level, with kinematic selections mimicking those applied after reconstruction. The
predictions of the sherpa generator are compared to those of powheg+pythia8, of
amc@nlo+herwig++ [84] and of alpgen+herwig. Comparisons are also made between samples generated with amc@nlo with renormalisation (µR ) and factorisation (µF )scales8 independently modified by factors of 2 or 0.5 and also with different PDF sets
8

The nominal scales are taken as µR = µF = [m2W + pT (W )2 + m2b + (pT (b)2 + pT (b)2 )/2]1/2 .

– 38 –JHEP01(2015)069

Z+jets background. As explained in section 7, ∆φ(jet1 , jet2 ) and pZ
T reweightings areapplied to the Zl and Zc+Zb components, respectively. For the ∆φ(jet1 , jet2 ) reweighting,a systematic uncertainty amounting to half of the correction is assigned to the Zl component, while an uncertainty amounting to the full correction is assigned to the Zc + Zb
components. This is done separately for 2- and 3-jet events, and all these uncertainties aretreated as uncorrelated. For the pZ
T reweighting, uncorrelated systematic uncertainties ofhalf the correction are assigned to the Zl and Zc + Zb components. The notation Zc + Zb
is meant to indicate that a systematic uncertainty is treated as correlated between the Zc
and Zb components.
The normalisation and the 3-to-2-jet ratio for the Zl background are determined fromdata in the 0-tag region of the 2-lepton channel, both with an uncertainty of 5%. The
normalisations of the Zcl and Zbb backgrounds are left free in the global fit. The uncertainties on the 3-to-2-jet ratios for the Zcl and Z+hf components are assessed througha comparison of alpgen with the nominal sherpa generator in the 2-tag region of the
2-lepton channel; these are 26% for Zcl and 20% for Z+hf. The same procedure is used toestimate uncertainties on the flavour fractions within Z+hf events, yielding 12% for eachof bl/bb, cc/bb and bc/bb, with bl/bb uncorrelated between 2- and 3-jet samples.
The shape of the mbb distribution is compared between data and simulation in the
2-tag region of the 2-lepton channel, excluding the 100–150 GeV range, from which a shape
uncertainty is derived that, when it increases the dijet-mass distribution by 3% at 50 GeV,
it decreases it by 5% at 200 GeV. This uncertainty is applied uncorrelated to the Zl andZb + Zc components. The differences between alpgen and sherpa are covered by thisuncertainty.

Diboson background. The uncertainties on the cross sections for diboson production(W W , W Z and ZZ) are assessed at parton level using mcfm at NLO in QCD. The sourcesof uncertainty considered are the renormalisation and factorisation scales and the choice of
PDFs. The nominal scales are dynamically set to half of the invariant mass of the diboson
system and the nominal PDFs are the CT10 set.
The scale uncertainties are evaluated by varying simultaneously µR and µF by factorsof 2 or 0.5. Since the analysis is performed in pVT intervals and in exclusive 2- and 3-jet
categories, the uncertainties are evaluated for each channel separately in those intervals andcategories (2 and 3 final-state partons within the nominal selection acceptance) following
the prescription of ref. [86]. This procedure leads, in each pVT interval, to two uncorrelated
uncertainties in the 2-jet category, one for 2+3 jets inclusively and one associated withthe removal of 3-jet events, and to one in the 3-jet category anti-correlated with the latter
uncertainty in the 2-jet category. These uncertainties are largest at high pVT . For pVT >200 GeV, the two uncertainties affecting the 2-jet category can be as large as 29% and 22%
in the W Z channel, roughly half this size in the ZZ channel and intermediate for W W ;and the uncertainty affecting the 3-jet category is about 17% in all channels.
The uncertainties due to the PDF choice are evaluated according to the PDF4LHC recommendation [87], i.e., using the envelope of predictions from the CT10, MSTW2008NLO,
and NNPDF2.3 PDF sets and their associated uncertainties. They range from 2% to 4%,with no pVT dependence observed.
The shape of the reconstructed Z → bb lineshape in V Z production is affected by
the parton-shower and hadronisation model. A shape-only systematic uncertainty is assessed by comparing the lineshapes obtained with the nominal powheg+pythia8 gener-– 39 –JHEP01(2015)069(CT10, MSTW2008NLO and NNPDF2.3 [85]). As a result, a 10% uncertainty is assigned
to the 3-to-2-jet ratio, taken as correlated between all W +hf processes. Shape uncertainties are also assessed for the mbb and pW
T distributions. When the former increases thedijet-mass distribution by 23% at 50 GeV, it decreases it by 28% at 200 GeV. It is taken
as uncorrelated for W l, W cl, W bb + W cc and W bl + W bc. For W bb + W cc, it is furthermore uncorrelated among pW
T intervals (with the three highest intervals correlated for thedijet-mass analysis). When the latter shape uncertainty increases the pW
T distribution by
9% for pW
=
50
GeV,
it
decreases
it
by
23%
at
200
GeV.
It
is
taken
as
correlated for allT
W +hf processes, and uncorrelated between the 2- and 3-jet samples.
Predictions using the inclusive production of all flavours by sherpa and alpgen arecompared after full reconstruction and event selection to assign uncertainties on the flavour
fractions that take properly into account heavy-flavour production at both the matrixelement and parton-shower levels. (For alpgen, the production of light flavours and heavy
flavours are performed separately at the matrix-element level; a dedicated procedure, basedon the ∆R separation between b-partons, is used to remove the overlap between bb pairs
produced at the matrix-element and parton-shower levels.) The following uncertainties are
assigned in the W +hf samples: 35% for bl/bb and 12% for each of bc/bb and cc/bb. The
uncertainty on bl/bb is uncorrelated between pW
T intervals (with the three highest intervals
correlated for the dijet-mass analysis).

ator and with herwig. The relative difference between the shapes is 20% for a dijet massaround 125 GeV.
8.4

Uncertainties on the signal modelling– 40 –JHEP01(2015)069

The qq → W H, qq → ZH, and gg → ZH signal samples are normalised respectively totheir inclusive cross sections as explained in section 3. The uncertainties on these crosssections [88] include those arising from the choice of scales µR and µF and of PDFs.
The scale uncertainty is 1% for W H production. It is larger (3%) for ZH production,
due to the contribution of the gluon-gluon initiated process. Under the assumption that thescale uncertainties are similar (1%) for qq → W H and qq → ZH, a conservative uncertainty
of 50% is inferred for gg → ZH. The same procedure leads to PDF uncertainties of 2.4%
for qq → (W/Z)H and 17% for gg → ZH. The relative uncertainty on the Higgs boson
branching ratio to bb is 3.3% for mH = 125 GeV [11]. The contribution of decays to final
states other than bb is verified to amount to less than 1% after selection.
Acceptance uncertainties due to the choice of scales are determined from signal samples
generated with powheg interfaced to pythia8, with µR and µF varied independently byfactors of 2 or 0.5. The procedure advocated in ref. [86] is used, after kinematic selections
applied at generator level, leading to acceptance uncertainties of 3.0%, 3.4% and 1.5% for
qq → W H, qq → ZH and gg → ZH, respectively, for the 2- and 3-jet categories combined,
and of 4.2%, 3.6% and 3.3% for the 3-jet category. The latter uncertainty is anti-correlatedwith an acceptance uncertainty associated with the removal of 3-jet events from the 2+3-jet
category to form the 2-jet category. In addition, the pVT spectrum is seen to be affected, andshape uncertainties are derived. For the qq → (W/Z)H samples, when they increase the
distribution by 1% for pVT = 50 GeV, they decrease it by 3% at 200 GeV. These variationsare 2% and 8%, respectively, for the gg → ZH samples.
Acceptance uncertainties due to the PDF choice are determined in a similar way,following the PDF4LHC prescription. They range from 2% in the 2-jet gg → ZH samples
to 5% in the 3-jet qq → ZH samples. There is no evidence of a need for pVT shape
uncertainties related to the PDFs.
The applied uncertainties on the shape of the pVT spectrum associated with the NLO
electroweak corrections [38] are typically at the level of 2%, increasing with pVT to reach
2.5% in the highest pVT interval.
The effect of the underlying-event modelling is found to be negligible, using various
pythia tunes. The effect of the parton-shower modelling is examined by comparison ofsimulations by powheg interfaced with pythia8 and with herwig. Acceptance variationsof 8% are seen, except for 3-jet events in the pVT > 120 GeV interval, where the variation isat the level of 13%. These variations are taken as systematic uncertainties.
A summary of the systematic uncertainties affecting the modelling of the Higgs boson
signal is given in table 5.

SignalTable 5. Summary of the systematic uncertainties on the signal and background modelling. An“S” symbol is used when only a shape uncertainty is assessed.

– 41 –JHEP01(2015)069

Cross section (scale)
1% (qq), 50% (gg)
Cross section (PDF)
2.4% (qq), 17% (gg)
Branching ratio
3.3 %Acceptance (scale)
1.5%–3.3%
3-jet acceptance (scale)
3.3%–4.2%
V
pT shape (scale)S
Acceptance (PDF)
2%–5%pVT shape (NLO EW correction)
S
Acceptance (parton shower)
8%–13%
Z+jetsZl normalisation, 3/2-jet ratio5%
Zcl 3/2-jet ratio
26%
Z+hf 3/2-jet ratio
20%
Z+hf/Zbb ratio
12%∆φ(jet1 , jet2 ), pVT , mbb
S
W +jets
W l normalisation, 3/2-jet ratio
10%
W cl, W +hf 3/2-jet ratio
10%
W bl/W bb ratio
35%
W bc/W bb, W cc/W bb ratio
12%∆φ(jet1 , jet2 ), pVT , mbb
Stt
3/2-jet ratio
20%
High/low-pVT ratio
7.5%
miss
Top-quark pT , mbb , ET
S
Single top
Cross section
4% (s-,t-channel), 7% (W t)
Acceptance (generator)
3%–52%b1
mbb , pTS
DibosonCross section and acceptance (scale)
3%–29%
Cross section and acceptance (PDF)
2%–4%
mbb
S
Multijet0-, 2-lepton channels normalisation100%
1-lepton channel normalisation
2%–60%
Template variations, reweightingS9

Statistical procedure

9.1

General aspects

9

While keeping distinct MM and TT categories in the 1-lepton channel improves the sensitivity, this is
not observed for the 0- and 2-lepton channels. Keeping the LL category separated from the others improves
the sensitivity in all lepton channels.

– 42 –JHEP01(2015)069

A statistical fitting procedure based on the Roostats framework [89, 90] is used to extract
the signal strength from the data. The signal strength is a parameter, µ, that multiplies
the SM Higgs boson production cross section times branching ratio into bb. A binned
likelihood function is constructed as the product of Poisson-probability terms over the binsof the input distributions involving the numbers of data events and the expected signal andbackground yields, taking into account the effects of the floating background normalisations
and the systematic uncertainties.
The different regions entering the likelihood fit are summarised in table 6. In the dijetmass analysis, the inputs to the “global fit” are the mbb distributions in the 81 2-tag signal
regions defined by three channels (0, 1 or 2 leptons), up to five pVT intervals, two numberof-jet categories (2 or 3), and three b-tagging categories (LL, MM and TT). Here and in therest of this section, mbb distributions are to be understood as transformed distributions, as
explained in section 6. In the MVA, the inputs are the BDTV H discriminants in the 24 2-tag
signal regions defined by the three lepton channels, up to two pVT intervals, the two numberof-jet categories, and b-tagging categories. In the 1-lepton channel, the b-tagging categories
are LL, MM and TT. In the 0- and 2-lepton channels, they are the LL category and a
combined MM and TT category (MM+TT).9 These BDTV H -discriminant distributions
are supplemented by the three mbb distributions in the 100–120 GeV pVT interval of the
2-jet 2-tag categories (LL, MM, and TT) of the 0-lepton channel. For the MVA, additional
inputs are the M V 1c distributions of the b-tagged jet in the 11 1-tag control regions of
the MVA selection and in the 100–120 GeV pVT interval of the 2-jet 1-tag category of the
0-lepton channel. In the dijet-mass analysis, the M V 1c distributions are combined in eachof the pVT < 120 GeV and pVT > 120 GeV intervals, which also results in 11 1-tag control
regions. Altogether, there are 584 mbb and M V 1c bins in the 92 regions of the dijet-mass
analysis, and 251 BDTV H -discriminant and M V 1c bins in the 38 regions of the MVA, tobe used in the global fits.
The impact of systematic uncertainties on the signal and background expectations is
described by nuisance parameters (NPs), θ, which are constrained by Gaussian or lognormal probability density functions, the latter being used for normalisation uncertainties
to prevent normalisation factors from becoming negative in the fit. The expected numbersof signal and background events in each bin are functions of θ. The parameterisation ofeach NP is chosen such that the predicted signal and background yields in each bin are
log-normally distributed for a normally distributed θ. For each NP, the prior is added as apenalty term to the likelihood, L(µ, θ), which decreases it as soon as θ is shifted away fromits nominal value. The statistical uncertainties of background predictions from simulation
are included through bin-by-bin nuisance parameters.

Channel
1-tag
LL
MM 2-tagTT

Dijet-mass analysis
0-lepton 1-lepton 2-lepton
M V 1cmbb
mbb
mbb

0-lepton
BDT(∗)
BDT(∗)MVA
1-lepton 2-lepton
M V 1c
BDT
BDT
BDT
BDT

The test statistic qµ is then constructed from the profile likelihood ratioˆqµ = −2 ln Λµ with Λµ = L(µ, θ̂ µ )/L(µ̂, θ̂),
where µ̂ and θ̂ are the parameters that maximise the likelihood with the constraint 0 ≤ˆµ̂ ≤ µ, and θ̂ µ are the nuisance parameter values that maximise the likelihood for a givenµ. This test statistic is used for exclusion intervals derived with the CLs method [91,92]. To measure the compatibility of the background-only hypothesis with the observed
data, the test statistic used is q0 = −2 ln Λ0 . The results are presented in terms of: the
95% confidence level (CL) upper limit on the signal strength; the probability p0 of the
background-only hypothesis; and the best-fit signal-strength value µ̂ with its associateduncertainty σµ . The fitted µ̂ value is obtained by maximising the likelihood function withrespect to all parameters. The uncertainty σµ is obtained from the variation of 2 ln Λµ by
one unit, where Λµ is now defined without the constraint 0 ≤ µ̂ ≤ µ. Expected results are
obtained in the same way as the observed results by replacing the data in each input binby the expectation from simulation with all NPs set to their best-fit values, as obtained
from the fit to the data.10
While the analysis is optimised for a Higgs boson of mass 125 GeV, results are also
extracted for other masses. These are obtained without any change to the dijet-mass
analysis, except for the binning of the transformed mbb distribution, which is reoptimised.
For the MVA, it is observed that the performance degrades for masses away from 125 GeV,
for which the BDTs are trained. This is largely due to the fact that mbb is an input tothe BDTs. The MVA results for other masses are therefore obtained using BDTs retrainedfor each of the masses tested at 5 GeV intervals between 100 and 150 GeV. The details
provided in the rest of this section refer to the analysis performed for a Higgs boson massof 125 GeV.
10This type of pseudo-data sample is referred to as an Asimov dataset in ref. [92].

– 43 –JHEP01(2015)069

Table 6. The distributions used in each region by the likelihood fit in the dijet-mass analysis andin the MVA applied to the 8 TeV data. Here, “BDT” stands for “BDTV H discriminant”. For eachentry listed, there are additional divisions into pVT intervals: five in the dijet-mass analysis and twoin the MVA, as shown in table 2. These distributions are input to the fit for the 2-jet and 3-jet
categories separately, except in the low pVT interval (100–120 GeV) of the 0-lepton channel where
only the 2-jet category is used. In the 0- and 2-lepton channels, the MM and TT 2-tag categories
are combined in the MVA. (∗) In the low pVT interval of the 0-lepton channel, the MVA uses the
mbb distributions in the LL, MM and TT 2-tag categories as well as the M V 1c distribution in the
1-tag category.

Process
tt 0-lepton
tt 1-lepton
tt 2-lepton
W bb
W cl
Zbb
Zcl

Scale factor
1.36 ± 0.14
1.12 ± 0.09
0.99 ± 0.04
0.83 ± 0.15
1.14 ± 0.10
1.09 ± 0.05
0.88 ± 0.12

9.2

Technical details

The data have sufficient statistical power to constrain the largest background-normalisation
NPs, which are left free to float in the fit. This applies to the tt, W bb, W cl, Zbb and Zcl
processes. The corresponding factors applied to the nominal background normalisations
as resulting from the global fit of the MVA to the 8 TeV data, are shown in table 7. As
stated in section 8, the tt background is normalised in the 2-jet category independently ineach of the lepton channels. The reason for uncorrelating the normalisations in the three
lepton channels is that the regions of phase space probed in the 2-jet category are very
different between the three channels. In the 2-lepton channel, the tt background is almost
entirely due to events in which both top quarks decay into (W → `ν)b (fully leptonic
decays) with all final-state objects detected (apart from the neutrinos). In the 1-lepton
channel, it is in part due to fully leptonic decays with one of the leptons (often a τ lepton)
undetected, and in part to cases where one of the top quarks decays as above and the otherinto (W → qq 0 )b (semileptonic decays) with a missed light-quark jet. Finally, in the 0lepton channel, the main contributions are from fully leptonic decays with the two leptons
undetected and from semileptonic decays with a missed lepton and a missed light-quark
jet; here again, the missed leptons are often τ leptons. Futhermore, the pVT range probed
is different in the 0-lepton channel: pVT > 100 GeV in contrast to being inclusive in the 1and 2-lepton channels.
As described in detail in section 8, a large number of sources of systematic uncertainty
are considered. The number of nuisance parameters is even larger because care is taken to
appropriately uncorrelate the impact of the same source of systematic uncertainty acrossbackground processes or across regions accessing very different parts of phase space. Thisavoids unduly propagating constraints. For instance, the tt background contributes quite
differently in the 2-tag 3-jet regions of the 0- and 1-lepton channels on one side, and of
the 2-lepton channel on the other. In the 0- and 1-lepton channels, it is likely that a jetfrom a t → b(W → qq) decay is missed, while in the 2-lepton channel it is likely that an
ISR or FSR jet is selected. This is the reason for not correlating, between these two setsof lepton channels, the systematic uncertainty attached to the 3-to-2 jet ratio for the tt– 44 –JHEP01(2015)069

Table 7. Factors applied to the nominal normalisations of the tt, W bb, W cl, Zbb, and Zcl backgrounds, as obtained from the global MVA fit to the 8 TeV data. The tt background is normalisedin the 2-jet category independently in each of the lepton channels. The errors include the statisticaland systematic uncertainties.

– 45 –JHEP01(2015)069

background. Another example is the ∆φ reweighting in the W +jets processes, which isderived in the 0-tag sample and applied to the W cl and W l backgrounds. As explained
in section 7, this reweighting is not applied to the W cc and W b backgrounds but, in theabsence of further information, an uncertainty is assessed for the ∆φ distributions of theW cc and W b backgrounds, uncorrelated with the uncertainty applied to the W cl and W lbackgrounds. Altogether, the fit has to handle almost 170 NPs, with roughly half of thosebeing of experimental origin.
The fit uses templates constructed from the predicted yields for the signal and thevarious backgrounds in the bins of the input distribution in each region. The systematic
uncertainties are encoded in templates of variations relative to the nominal template foreach up-and-down (±1σ) variation. The limited size of the MC samples for some simulated background processes in some regions can cause large local fluctuations in templatesof systematic variations. When the impact of a systematic variation translates into a
reweighting of the nominal template, no statistical fluctuations are expected beyond thosealready present in the nominal template. This is the case, for instance, for the b-tagging
uncertainties. For those, no specific action is taken. On the other hand, when a systematic
variation may introduce changes in the events selected, as is the case for instance with
the JES uncertainties, additional statistical fluctuations may be introduced, which affectthe templates of systematic variations. In such cases, a smoothing procedure is applied
to each systematic-variation template in each region. Bins are merged based on the constraints that the statistical uncertainty in each bin should be less than 5% and that theshapes of the systematic-variation templates remain physical: monotonous for a BDTV H
discriminant, and with at most one local extremum for a dijet mass.
Altogether, given the number of regions and NPs, the number of systematic-variation
template pairs (+1σ and −1σ) is close to twenty thousand, which renders the fits highly
time consuming. To address this issue, systematic uncertainties that have a negligible
impact on the final results are pruned away, region by region. A normalisation (shape)
uncertainty is dropped if the associated template variation is below 0.5% (below 0.5% in
all bins). Additional pruning criteria are applied to regions where the signal contribution isless than 2% of the total background and where the systematic variations impact the total
background prediction by less than 0.5%. Furthermore, shape uncertainties are dropped
if the up- and down-varied shapes are more similar to each other than to the nominal
shape. This is only done for those systematic uncertainties where opposite-sign variations
are expected. This procedure reduces the number of systematic-variation templates by afactor of two.
The behaviour of the global fit is evaluated by a number of checks, including how mucheach NP is pulled away from its nominal value, how much its uncertainty is reduced withrespect to its nominal uncertainty, and which correlations develop between initially uncorrelated systematic uncertainties. To assess these effects, comparisons are made between
the expectations from simulation and the observations in the data. When differences arise,
their source is investigated, and this leads in a number of cases to uncorrelating further
systematic uncertainties by means of additional NPs. This is to prevent a constraint frombeing propagated from one kinematic region to another if this is not considered well moti-

9.3

Cross checks using diboson production

Diboson production with a Z boson decaying to a pair of b-quarks and produced in association with either a W or Z boson has a signature very similar to the one consideredin this analysis, but with a softer pbb
T spectrum and with a mbb distribution peaking at
lower values. The cross section is about five times larger than for the SM Higgs boson witha mass of 125 GeV. Diboson production is therefore used as a validation of the analysis
procedure. For the dijet-mass analysis, the binning of the transformed mbb distribution isreoptimised for the Z boson mass. For the MVA, the BDTs are retrained to discriminate– 46 –JHEP01(2015)069

vated. The stability of the results is also tested by performing fits for each lepton channel
independently, which can also help to identify from which region each constraint originates.
It is particularly useful to understand which systematic uncertainties have the largest
impact on the final results, and therefore should be considered with greater care. For thispurpose, a so-called ranking of the NPs is established. For each systematic uncertainty, thefit is performed again with the corresponding NP fixed to its fitted value, θ̂, shifted up or
down by its fitted uncertainty, with all the other parameters allowed to vary so as to takeproperly into account the correlations between systematic uncertainties. The magnitudeof the shift in the fitted signal strength µ̂ is a measure of the observed impact of theconsidered NP. The same procedure is repeated, using the nominal values of the NP andof its associated uncertainty to provide its expected impact. To reduce the computation
time and therefore to enable more detailed fit studies, some of the NPs which have anegligible effect on the expected fitted uncertainty on µ̂ are dropped: those associatedwith the muon momentum scale and resolution and with the electron energy resolution;
one of those associated with the jet energy scale; and those associated with the quarkgluon composition of the backgrounds, which turn out to be fully correlated with those
associated with the difference in energy response between quark and gluon jets. The
ranking of the systematic uncertainties obtained with the MVA applied to the 8 TeV datais shown in figure 16 with the NPs ordered by decreasing post-fit impact on µ̂. The fivesystematic uncertainties with the largest impact are, in descending order, those: on thedijet-mass shape for the W bb and W cc backgrounds for pW
T > 120 GeV; on the W bl/W bb
Wnormalisation ratio for pT > 120 GeV; on the W bb background normalisation; on the pWTshape in the 3-jet category for the W +hf background; and on the signal acceptance due tothe parton-shower modelling.
Since the same data sample is used for both the dijet-mass analysis and the MVA, theconsistency of the two final results, i.e., the two fitted signal strengths, is assessed using
the “bootstrap” method [93]. A large number of event samples are randomly extracted
from the simulated samples, with the signal strength µ set to unity, the SM value. Each ofthem is representative of the integrated luminosity used for the data analysis in terms of
expected yields as well as of associated Poisson fluctuations. Each of these event samples
is subjected to both the dijet-mass analysis and the MVA, thus allowing the two fittedµ̂ values to be compared and their statistical correlation to be extracted. At the sametime, the expected distributions of µ̂ and of its uncertainty are determined for both thedijet-mass analysis and the MVA.

∆µ-0.15 -0.1 -0.05

00.05 0.1 0.15 0.2

W+bb, W+cc mjj shape(pV > 120 GeV)T

W+bl to W+bb normalisation(pV > 120 GeV)
T

W+bb normalisation
W+HF pV shape (3-jet)
Signal acceptance (parton shower)
Z+bl to Z+bb normalisation (2-jet)b-jet energy resolution
Z+bb, Z+cc mjj shape
Jet energy resolution
Dilepton tt normalisation
W+HF pV shape (2-jet)T

Z+bb normalisation
Jet energy scale 1
b-jet tagging efficiency 4
ttbar high pV normalisation
T

ATLAS

Pull: (θ - θ0)/∆θ

s = 8 TeV,

∫ Ldt = 20.3 fb-1

mH=125 GeV-2 -1.5-4-1 -2-0.50Normalisation+1σ Postfit Impact on µ
-1σ Postfit Impact on µ

0.5 2 1

1.5
4

2 6

Figure 16. Impact of systematic uncertainties on the fitted signal-strength parameter µ̂ for the
MVA applied to the 8 TeV data. The systematic uncertainties are listed in decreasing order of theirimpact on µ̂ on the y-axis. The boxes show the variations of µ̂, referring to the top x-axis, whenfixing the corresponding individual nuisance parameter θ to its post-fit value θ̂ modified upwards or
downwards by its post-fit uncertainty, and repeating the fit as explained in the text. The hatchedand open areas correspond to the upwards and downwards variations, respectively. The filled circles, referring to the bottom x-axis, show the deviations of the fitted nuisance parameters θ̂ fromtheir nominal values θ0 , expressed in terms of standard deviations with respect to their nominal
uncertainties ∆θ. The associated error bars show the post-fit uncertainties of the nuisance parameters, relative to their nominal uncertainties. The open circles with their error bars, also referringto the bottom x-axis, show the fitted values and uncertainties of the normalisation parameters thatare freely floating in the fit. The normalisation parameters have a pre-fit value of one. As explained
in section 8.1, the jet energy scale and b-tagging uncertainties are decomposed into uncorrelated
components; the labels 1 and 4 refer to such components.

– 47 –JHEP01(2015)069

T– 48 –JHEP01(2015)069

the diboson signal from all backgrounds (including the Higgs boson). So-called “V Z fits”
are performed, where the normalisation of the diboson contributions is allowed to varywith a multiplicative scale factor µV Z with respect to the SM expectation, except for thesmall contribution from W W production, which is treated as a background and constrained
within its uncertainty. A SM Higgs boson with mH = 125 GeV is included as a background,with a production cross section at the SM value with an uncertainty of 50%. Distributionsof the BDTV Z discriminants of the MVA are shown in figure 17 for 2-tag signal regionswith pVT > 120 GeV in the 2-jet category of the 0-, 1- and 2-lepton channels.
As an additional check, fits are also performed with both the diboson and Higgs boson
signal-strength parameters µV Z and µ left freely floating, to study the correlation between
the two strength parameters. The fits in the dijet-mass analysis use the mbb distributions
with binning optimised for a Higgs boson mass of 125 GeV. The fits in the MVA use BDTs
trained for that same mass, as well as the associated optimised binnings.

5

-1

T

0

0.4 0.6 0.8

1

0.2 0.4 0.6 0.8
1
BDT discriminant, VZ

0.2

Data 2012
VH(bb) (µ=1.0)
VZ
tt
Single top
W+hf
W+cl
Z+hf
Z+cl
Z+l
Uncertainty
Pre-fit background
VH(bb)× 40
VZ× 10
T

-1

(b)

0

0

0.4 0.6 0.8

Data 2012
VH(bb) (µ=1.0)
VZ
tt
Single top
Multijet
W+hf
W+cl
W+l
Z+hf
Uncertainty
Pre-fit background
VH(bb)× 40
VZ× 301

0.2 0.4 0.6 0.8
1
BDT discriminant, VZ

0.2

s = 8 TeV ∫Ldt = 20.3 fb
1 lep., 2 jets, 2 Medium+Tight tags
p V >120 GeV

ATLAS

1.2-1 -0.8 -0.6 -0.4 -0.2
10.8
-1 -0.8 -0.6 -0.4 -0.2

10

102

3

10

104

105

T

-1

(c)

0

0

0.4 0.6 0.8

Data 2012
VH(bb) (µ=1.0)
VZ
tt
Single top
Z+hf
Z+cl
Z+l
Uncertainty
Pre-fit background
VH(bb)× 40
VZ× 201

0.2 0.4 0.6 0.8
1
BDT discriminant, VZ

0.2

s = 8 TeV ∫Ldt = 20.3 fb
2 lep., 2 jets, 2 Medium+Tight tags
p V >120 GeV

ATLAS

1.5-1 -0.8 -0.6 -0.4 -0.2
10.5
-1 -0.8 -0.6 -0.4 -0.2

10-1

1

10

102

3

10

104

Figure 17. Distribution of the BDTV Z discriminant used in the V Z fit observed in data (points with error bars) and expected (histograms) for
the 2-jet 2-tag signal regions with pVT > 120 GeV in the Medium (MM) and Tight (TT) categories combined of (a) the 0-lepton channel, (b) the
1-lepton channel, (c) the 2-lepton channel. The V Z signal and background contributions are shown as filled histograms after the global V Z fit of
the MVA, except for the Higgs boson contribution (mH = 125 GeV). The latter is shown as expected from the SM (indicated as µ = 1.0) both asa filled histogram on top of the sum of the fitted V Z signal and other backgrounds, and unstacked as an unfilled histogram scaled by the factor
indicated in the legend. The expected V Z contribution is shown in a similar way as an unfilled dotted histogram. The dashed histogram shows
the sum of the total V Z signal and background as expected from the pre-fit MC simulation. The size of the combined statistical and systematic
uncertainty on the fitted V Z signal and background is indicated by the hatched band. The ratio of the data to the sum of the fitted V Z signal
and background is shown in the lower panel.

(a)

1
0.5
-1 -0.8 -0.6 -0.4 -0.2

0s = 8 TeV ∫Ldt = 20.3 fb
0 lep., 2 jets, 2 Medium+Tight tags
p V >120 GeV

ATLAS

1
1.5-1 -0.8 -0.6 -0.4 -0.2

10

102

10

3

104

10
Events / 0.25
Data/Pred

Events / 0.15

Data/Pred

Events / 0.25
Data/Pred

JHEP01(2015)069– 49 –

10

Analysis of the 7 TeV data

For the 7 TeV dataset, only a dijet-mass analysis is performed. It is similar but not identicalto the corresponding analysis for the 8 TeV data, since some of the object reconstruction
tools, such as the simultaneous use of multiple b-tagging operating points, are not available
for the 7 TeV data. In this section, the main differences between the two analyses are
summarised.
10.1

Object reconstruction

10.2

Event selection

The selection criteria are those used in the dijet-mass analysis of the 8 TeV data, withthe following differences. With only one b-tagging operating point, a single 2-tag category
is defined. In the 0-lepton channel, the 100–120 GeV pVT interval is not used, and the
P jetimiss trigger is usedcriterion for
pT is not applied. In the 1-muon sub-channel, the ET
only in the 2-jet 2-tag category for pW
T > 160 GeV, and the events selected only by themiss trigger constitute distinct signal regions. In the 1-lepton channel, mW > 40 GeV isETTmiss > 25 GeV is imposed
required for pW
<
160
GeV;
there
is
no
requirement
on
H
,
but
E
T
T
T
for pW
<
200
GeV.
In
the
2-leptonchannel,nokinematicfit
isperformed.
Different leptonT
flavour events are used to define a 2-tag tt-dominated e–µ control region in the 2-lepton
channel; the region is defined to be inclusive in jet multiplicity (≥ 2).
10.3

Background composition and modelling

The templates used to model the MJ background in the 1-lepton channel are obtained byinversion of the track-based isolation criterion, and the normalisations are performed onmiss distributions in the electron and muon sub-channels, respectively.
the mW
T and ET
Corrections to the simulation of the V +jet backgrounds are determined in the 1- and2-lepton 0-tag samples inclusively in pVT , and applied as ∆φ(jet1 , jet2 ) reweightings to the
W +jet and Z+jet components in all channels. Selected dijet-mass distributions showing
the background composition in various analysis regions are shown in figure 18.

– 50 –JHEP01(2015)069

The three categories of electrons are selected according to the loose, medium, and tightcriteria defined in ref. [53]. The transverse energy threshold for loose electrons is set
at 10 GeV, instead of 7 GeV. For tight electrons and muons, the calorimeter isolation
requirement is loosened from 4% to 7%. The procedure used to avoid double-countingof reconstructed muon and jet objects removes muons separated by ∆R < 0.4 from any
jet, irrespective of the multiplicity of tracks associated with the jet. For jets, the global
sequential calibration is not used and the requirement on the fraction of track pT carriedby tracks originating from the primary vertex is raised from 50% to 75%. The b-tagging
algorithm used is MV1 [94–97] instead of MV1c, with a single operating point to defineb-tagged jets corresponding to an efficiency of 70%.

10.4

Systematic uncertainties

10.5

Statistical procedure

The inputs to the likelihood fits are the mbb distributions (not transformed) in the 28 pVT
intervals of the 2-tag signal regions. Additional inputs are the event yields in the five pVTintervals of the 2-tag e–µ control region and the 26 pVT intervals of the 1-tag control regions.
For the tt background, a single floating normalisation is determined by the global fit,instead of one in each of the 0-, 1-, and 2-lepton channels. In addition to the other floating
normalisations mentioned for the 8 TeV data analysis, the MJ background normalisation is
also left freely floating in all regions of the 1-lepton channel, except in the 2-tag 3-jet regions
where the statistical power of the data is not sufficient to provide a reliable constraint. In
these regions, an uncertainty of 30% is assigned to the MJ background normalisation, using
a method similar to what is done for the analysis of the 8 TeV data.

– 51 –JHEP01(2015)069

The differences with respect to the 8 TeV data analysis arise mainly from experimental
systematic uncertainties. Many of them are evaluated using independent data samples
miss trigger efficiencies or JES. Others refer to different(7 TeV data vs. 8 TeV data), e.g., ET
identification algorithms, e.g., electron identification or b-tagging. The uncertainty on the
integrated luminosity is 1.8% for the 2011 dataset [20].
The uncertainties affecting the signal and background simulation are estimated ina similar way as for the 8 TeV data, i.e., from comparisons between the baseline andalternative generators. For V +jets, the V bc and V bb backgrounds are merged into a single
component. For dibosons, the baseline generator is herwig instead of pythia8; systematic
uncertainties on the 3-to-2-jet ratios and on the pVT distributions are estimated at generator
level for the different diboson processes by comparison with mcfm at NLO. For the signal,
the gg → ZH samples are generated with pythia8 instead of powheg; for all processes,
the alternative generators used are pythia6 and herwig.
Due to these differences, and because the phase space within which the systematic
uncertainties are evaluated is more restricted than for the MVA applied to the 8 TeV
data, all systematic uncertainties, except for the theoretical uncertainties on the signal, aretreated as uncorrelated between the analyses of the 7 TeV and 8 TeV data in the global fit
to the combined dataset, in which the MVA is used for the 8 TeV data.

(a)

150

T

50

(b)

-1

100

s = 7 TeV ∫Ldt = 4.7 fb
1 lep., 2 jets, 2 tags
p V >120 GeV

ATLAS

150

200
250
mbb [GeV]

Data 2011
VH(bb) ( µ =1.0)
Dibosontt
Single top
Multijet
W+hf
W+cl
Uncertainty
Pre-fit background
VH(bb)×10

0
2
1.5
1
0.5
0

2

4

6

8

10

12

14

16

18

T

50

(c)

-1

100

s = 7 TeV ∫Ldt = 4.7 fb
2 lep., 2 jets, 2 tags
p V >120 GeV

ATLAS

150

200
250
mbb [GeV]

Data 2011
VH(bb) ( µ =1.0)
Diboson
tt
Z+hf
Z+cl
Uncertainty
Pre-fit background
VH(bb)×10

Figure 18. The dijet-mass distributions observed in the 7 TeV data (points with error bars) and expected (histograms) for the 2-jet 2-tag
signal regions with the pVT > 120 GeV intervals combined: (a) 0-lepton channel, (b) 1-lepton channel, and (c) 2-lepton channel. The backgroundcontributions after the global fit of the dijet-mass analysis are shown as filled histograms. The Higgs boson signal (mH = 125 GeV) is shown asa filled histogram on top of the fitted backgrounds, as expected from the SM (indicated as µ = 1.0), and, unstacked as an unfilled histogram,scaled by the factor indicated in the legend. The dashed histogram shows the total background as expected from the pre-fit MC simulation. The
entries in overflow are included in the last bin. The size of the combined statistical and systematic uncertainty on the sum of the signal and fitted
background is indicated by the hatched band. The ratio of the data to the sum of the signal and fitted background is shown in the lower panel.

200
250
mbb [GeV]

10

10
2
1.5
1
0.5
0

20

20

1.5
1
0.5

30

50

60

70

80

30

Data 2011
VH(bb) ( µ =1.0)
Dibosontt
Single top
Multijet
W+hf
W+cl
Z+hf
Z+cl
Uncertainty
Pre-fit background
VH(bb)×10

40

100

-1

40

50

50

70

T

s = 7 TeV ∫Ldt = 4.7 fb
0 lep., 2 jets, 2 tags
p V >120 GeV

60

ATLAS

80

Events / 20 GeV
Data/Pred

Events / 20 GeV

Data/Pred

Events / 20 GeV
Data/Pred

JHEP01(2015)069– 52 –

11

Results

11.1

Nominal results

The nominal results are obtained from global fits using the MVA for the 8 TeV data andthe dijet-mass analysis of the 7 TeV data.
Distributions of the BDTV H discriminant and of M V 1c, with background normalisations and nuisance parameters adjusted by the global fit to the 8 TeV data were alreadypresented in section 7.4. Dijet-mass distributions in the 7 TeV data analysis were shownin section 10. Agreement between data and estimated background is observed within theuncertainties shown by the hatched bands.
Figure 19 shows the 95% CL upper limits on the cross section times branching ratio
for pp → (W/Z)(H → bb) in the Higgs boson mass range 110–140 GeV. The observed limit
for mH = 125 GeV is 1.2 times the SM value, to be compared to an expected limit, in theabsence of signal, of 0.8. For the 8 TeV (7 TeV) data only, the observed and expected limits
are 1.4 (2.3) and 0.8 (3.2), respectively.
The probability p0 of obtaining from background alone a result at least as signallike as the observation is 8% for a tested Higgs boson mass of 125 GeV; in the presenceof a Higgs boson with that mass and the SM signal strength, the expected p0 value is0.5%. This corresponds to an excess observed with a significance of 1.4σ, to be compared
to an expectation of 2.6σ. For the 8 TeV data alone, the observed and expected levelsof significance are 1.7σ and 2.5σ, respectively. For the 7 TeV data alone, the expectedsignificance is 0.7σ and there is a deficit rather than an excess in the data, as can be seenin figure 18. Figure 20 shows the p0 values in the mass range 110–140 GeV, as obtained forthe 7 TeV and 8 TeV combined dataset.
The fitted µ values for mH = 125 GeV are shown in figure 21 for the 7 TeV, 8 TeV andcombined datasets. With all lepton channels and data-taking periods combined, the fitted
value of the signal-strength parameter is11 µ = 0.51±0.31(stat.)±0.24(syst.). For the 8 TeV
11The uncertainties of the normalisations of the floating backgrounds are included in the systematic
uncertainties; their contribution is 0.07.

– 53 –JHEP01(2015)069

As explained in section 9, the results are obtained from maximum-likelihood fits to thedata, where the inputs are the distributions of final discriminants in the 2-tag signal regions and the M V 1c distributions of the b-tagged jet in the 1-tag control regions, withnuisance parameters either floating or constrained by priors. The final discriminants are
the transformed mbb for the dijet-mass analysis and the BDTV H discriminants for the MVA.
Results are extracted independently for the dijet-mass and multivariate analyses. Since the
MVA has better expected sensitivity to a Higgs boson signal, it is used for the nominal
results, while the dijet-mass analysis provides a cross-check (cf. section 11.2). For the
7 TeV data, however, only a dijet-mass analysis is performed. Unless otherwise specified,
all results refer to a Higgs boson mass of 125 GeV.
In the following, the fitted signal-strength parameters are simply denoted µ and µV Z ,rather than µ̂ and µ̂V Z .

95% C.L. limit on σ/ σSM

4

3.5

ATLAS
Observed (CLs)
Expected (no Higgs)
Expected (m = 125 GeV)
H
± 1σ
± 2σ

3

2.5

∫ Ldt = 4.7 fb
s = 8 TeV, ∫ Ldt = 20.3 fb-1

s = 7 TeV,

-1

2

1
0.5
0

110

115

120

125

130

135

140
mH [GeV]

Local p0

Figure 19. Observed (solid) and expected 95% CL cross-section upper limits, normalised to the SM
Higgs boson production cross section, as a function of mH for all channels and data-taking periods
combined, as obtained using the dijet-mass analysis for the 7 TeV dataset and BDTs trained ateach individual mass for the 8 TeV dataset. The expected upper limit is given for the backgroundonly hypothesis (dashed) and with the injection of a SM Higgs boson signal at a mass of 125 GeV(dotted). The dark and light shaded bands represent the 1σ and 2σ ranges of the expectation in
the absence of a signal. For all curves shown, the results obtained at the tested masses are linearlyinterpolated.
10

1

ATLAS

s=7 TeV ∫ Ldt = 4.7 fb-1

Observed
Expected
Injected 125 GeV

s=8 TeV ∫ Ldt = 20.3 fb

-1

1σ

10-1
2σ

10-2
3σ

10-3

10-4
110

115

120

125

130

135

140
mH [GeV]Figure 20. Observed (solid) and expected p0 values as a function of mH for all channels anddata-taking periods combined, as obtained using the dijet-mass analysis for the 7 TeV dataset and
BDTs trained at each individual mass for the 8 TeV dataset. The expected p0 values are given forthe background-only hypothesis in the presence of a SM Higgs boson: for the dashed curve the
Higgs boson mass corresponds to each tested mass point in turn; for the dotted curve the Higgs
boson mass is 125 GeV. For all curves shown, the results obtained at the tested masses are linearlyinterpolated.

– 54 –JHEP01(2015)069

1.5

ATLAS

s=7 TeV, ∫ Ldt=4.7 fb-1; s=8 TeV, ∫ Ldt=20.3 fb-1
tot.

tot ( stat

stat.

syst )+ 0.86 )-1.61 +− 1.50
( +− 1.22
1.13 − 0.92
1.46

8 TeV

+ 0.28 )0.65 +− 0.43
( +− 0.33
0.32 − 0.24
0.40

Combination+ 0.25 )0.51 +− 0.40
( +− 0.31
0.30 − 0.22
0.37

-4

-2

0

2

4

best fit µ=σ/ σ

SM

6

8

for mH=125 GeV

Figure 21. The fitted values of the Higgs boson signal-strength parameter µ for mH = 125 GeV
for the 7 TeV and 8 TeV datasets and the combination of the 7 TeV and 8 TeV datasets.

data, the fitted value of the signal-strength parameter is µ = 0.65±0.32(stat.)±0.26(syst.).
For the 7 TeV data, it is µ = −1.6 ± 1.2(stat.) ± 0.9(syst.).
For a Higgs boson with a mass of 125.36 GeV, as measured by ATLAS [98], the signalstrength parameter is µ = 0.52 ± 0.32(stat.) ± 0.24(syst.).
Fits are also performed where the signal strengths are floated independently for (i)the W H and ZH production processes, or (ii) the three lepton channels. The results ofthese fits are shown in figures 22 and 23 respectively. The consistency of the fitted signal
strengths in the W H and ZH processes is at the level of 20%. For the lepton channels,
the consistency between the three fitted signal strengths is at the level of 72% for the
7 TeV data, and of 8% for the 8 TeV data. The low values of the fitted signal strengths for
the ZH process and in the 0-lepton channel are associated with the data deficit observedin the most sensitive bins of the BDTV H discriminant in the 0-lepton channel, shown infigure 12(a).
Figure 24 shows the data, background and signal yields, where the final-discriminantbins in all signal regions are combined into bins of log(S/B), separately for the 7 and 8 TeV
datasets. Here, S is the expected signal yield and B is the fitted background yield. Detailsof the fitted values of the signal and of the various background components are provided
in table 8.
11.2

Cross-check with the dijet-mass analysis

The distributions of mbb in the dijet-mass analysis, with background normalisations andnuisance parameters adjusted by the global fit to the 8 TeV data were already presented
in section 7.3. Agreement between data and estimated background is observed within theuncertainties shown by the hatched bands.

– 55 –JHEP01(2015)069

7 TeV

ATLAS

s=7 TeV, ∫ Ldt=4.7 fb-1; s=8 TeV, ∫ Ldt=20.3 fb-1
tot.

tot ( stat

stat.

syst )+ 0.27 )0.05 +− 0.52
( +− 0.44
0.42 − 0.25
0.49

WH

+ 0.42 )1.11 +− 0.65
( +− 0.50
0.48 − 0.38
0.61

Combination+ 0.25 )0.51 +− 0.40
( +− 0.31
0.30 − 0.22
0.37

-1

0

1

2

3

4

best fit µ=σ/ σ

SM

5

6

7

for mH=125 GeV

Figure 22. The fitted values of the Higgs boson signal-strength parameter µ for mH = 125 GeV
for the W H and ZH processes and the combination of the W H and ZH processes, with the 7
and 8 TeV datasets combined. The individual µ values for the (W/Z)H processes are obtained
from a simultaneous fit with the signal strength for each of the W H and ZH processes floating
independently.

ATLAS

s=7 TeV, ∫ Ldt=4.7 fb-1; s=8 TeV, ∫ Ldt=20.3 fb-1
tot.

tot ( stat

stat.

syst )

2 lepton+ 0.50 )0.94 +− 0.88
( +− 0.72
0.68 − 0.41
0.79

1 lepton+ 0.43 )1.17 +− 0.66
( +− 0.50
0.48 − 0.37
0.60

0 lepton+ 0.26 )-0.35 +− 0.55
( +− 0.49
0.44 − 0.27
0.52

Combination+ 0.25 )0.51 +− 0.40
( +− 0.31
0.30 − 0.22
0.37

-1

0

1

2

3

4

best fit µ=σ/ σ

SM

5

6

7

for mH=125 GeV

Figure 23. The fitted values of the Higgs boson signal-strength parameter µ for mH = 125 GeV
for the 0-, 1- and 2-lepton channels and the combination of the three channels, with the 7 and 8 TeV
datasets combined.The individual µ values for the lepton channels are obtained from a simultaneous
fit with the signal strength for each of the lepton channels floating independently.

– 56 –JHEP01(2015)069

ZH

ATLAS

Events / 0.5

Events / 0.5

107

Data 2012
VH(bb) (µ =1.0)
Dibosontt
Single top
Multijet
W+hf
W+cl
W+l
Z+hf
Z+cl
Z+l

s = 8 TeV ∫ Ldt = 20.3 fb-1

6

10

5

10

104
3

10

Data 2011
VH(bb) (µ =1.0)
Dibosontt
Single top
Multijet
W+hf
W+cl
W+l
Z+hf
Z+cl
Z+l

s = 7 TeV ∫ Ldt = 4.7 fb-1

5

10

3

10

102
10
1
10-1

-4

-3.5

-3

-2.5

-2

-1.5

-1 -0.5

0

0.5

-4

-3.5

-3

-2.5

-2

-1.5

-1 -0.50 0.5log (S/B)

Pull (stat.)

1

2
0
-2

-4

-3.5

-3

-2.5

-2

-1.5

-1 -0.50 0.5log (S/B)
10

10

(a)

(b)

Figure 24. Event yields as a function of log(S/B) for data, background and Higgs boson signal
with mH = 125 GeV for the (a) 8 TeV data and (b) 7 TeV data. Final-discriminant bins in all signal
regions are combined into bins of log(S/B). The signal S and background B yields are expected
and fitted, respectively. The Higgs boson signal contribution is shown as expected for the SM crosssection (indicated as µ = 1.0). The pull of the data with respect to the background-only predictionis also shown with statistical uncertainties only. The full line indicates the pull of the prediction
for signal (µ = 1.0) and background with respect to the background-only prediction.

Process
Bin 1
Bin 2
Bin 3 Bin 4 Bin 5 Bin 6 Bin 7 Bin 8
Data368550 141166 111865 20740 5538 2245
382
41
Signal
29
43
96
57
58
62
32
10.7
Background 368802 140846 111831 20722 5467 2189
364
37.9
−5
S/B
8 × 100.0003 0.0009 0.003 0.01
0.03
0.09
0.3
W +hf
14584
10626 15297 1948
618
250
45
8.2
W cl
96282
30184 15227 1286
239
47
4.2
0.2
Wl
125676
14961
3722
588
107
16
1.3
0.03
Z+hf
10758
14167 21684 7458 1178
577
130
14.8
Zcl
13876
11048
4419
941
61
22
2.1
0.1
Zl
49750
18061
3044
537
48
15
1
0.05
tt
30539
24824 26729 5595 2238
922
137
10
Single top
10356
9492
14279 1494
688
252
31
2.7
Diboson
4378
1831
1247
474
186
62
9.7
1
Multijet
12603
5650
6184
400
103
26
3
0.9

Bin 9
4
2.3
3.4
0.7
0.7
0.005
0.001
2.2
0.008
0.004
0.3
0.1
0.2
0

Table 8. The numbers of expected signal and fitted background events and the observed numbersof events after MVA selection in the bins of figure 24(a). These numbers are for both the 1-tag and2-tag events in the 8 TeV dataset, corresponding to an integrated luminosity of 20.3 fb−1 .

– 57 –JHEP01(2015)069

10Pull (stat.)

ATLAS

104

102

2
0
-2

6

10

s = 8 TeV ∫Ldt = 20.3 fb -1
0+1+2 lep., 2+3 jets, 2 tags
Weighted by Higgs S/B

6
42

Weighted events after subtraction / 20.0 GeV

8

Data 2012
VH(bb) (µ=1.0)
DibosonUncertainty

ATLAS

1.2
1
0.8

Data 2011
VH(bb) (µ=1.0)
DibosonUncertainty

ATLAS
s = 7 TeV ∫Ldt = 4.7 fb -1
0+1+2 lep., 2+3 jets, 2 tags
Weighted by Higgs S/B

0.6
0.4
0.2
0

-0.2
-0.4

0

-0.6

50

100

150

200

250

mbb [GeV]
(a)

50

100

150

200
250
mbb [GeV](b)

Figure 25. The distribution of mbb in data after subtraction of all backgrounds except for thediboson processes, as obtained with the dijet-mass analysis for the (a) 8 TeV and (b) 7 TeV data. The
contributions from all lepton channels, pVT intervals, number-of-jets and 2-tag b-tagging categories
are summed weighted by their respective values of the ratio of expected Higgs boson signal to fitted
background. The contribution of the associated W H and ZH production of a SM Higgs bosonwith mH = 125 GeV is shown as expected for the SM cross section (indicated as µ = 1.0). The sizeof the combined statistical and systematic uncertainty on the fitted background is indicated by thehatched band.

In the dijet-mass analysis, a µ value of 1.23 ± 0.44(stat.) ± 0.41(syst.) is obtained forthe 8 TeV dataset. The consistency of the results of the three lepton channels is at the levelof 8%. Using the “bootstrap” method mentioned in section 9.2, the results for the 8 TeV
data with the dijet-mass analysis and with the MVA are expected to be 67% correlated,
and the observed results are found to be statistically consistent at the level of 8%. The
observed significance in the dijet-mass analysis is 2.2σ. The expected significance is 1.9σ,
to be compared to 2.5σ for the MVA, which is the reason for choosing the MVA for thenominal results.
Figure 25 shows the mbb distribution in data after subtraction of all backgrounds except
for diboson production for the 7 and 8 TeV data, as obtained with the dijet-mass analysis.
In this figure, the contributions of all 2-tag signal regions in all channels are summed
weighted by their respective ratios of expected Higgs boson signal to fitted background.
The V Z contribution is clearly seen, located at the expected Z mass. The Higgs boson
signal contribution is shown as expected for the SM cross section.
11.3

Cross-check with the diboson analysis

To validate the analysis procedures, V Z fits are performed, the technical details of whichwere discussed in section 9.3.

– 58 –JHEP01(2015)069

Weighted events after subtraction / 20.0 GeV

10

– 59 –JHEP01(2015)069

The measured signal strength for the 8 TeV dataset with the MVA is µV Z = 0.77 ±0.10(stat.) ± 0.15(syst.). This result is consistent with the observations already made on
figure 25. The signal strengths obtained for the three lepton channels are consistent at the
85% level. In the dijet-mass analysis at 8 TeV, a µV Z value of 0.79±0.11(stat.)±0.16(syst.)
is obtained. The correlation of the systematic uncertainties on µV Z and µ is 35% in the
MVA and 67% in the dijet-mass analysis.
Fits are performed with the same final discriminants as used to obtain the results for
the Higgs boson based on the 8 TeV dataset, but with both the V Z and Higgs boson signalstrength parameters µV Z and µ left freely floating. The result for the Higgs boson signal
strength is unchanged from the nominal result, and the statistical correlation between
the two signal-strength parameters is found to be −3% in the MVA and 9% in the dijetmass analysis. The main reason for these low correlations is the different shape of the pVT
distributions for V Z and for the Higgs boson signal, the pVT variable being used by both
the MVA and the dijet-mass analysis. The yield tables in the appendix show that the ratioof the diboson contribution to that of the Higgs boson is indeed smaller in the higher pVTinterval than in the lower one. The additional variables input to the BDT provide further
separation in the MVA, leading to a very small diboson contribution in the most significant
bins of the BDTV H discriminant, as can be seen in table 8.
A value of µV Z = 0.50 ± 0.30(stat.) ± 0.38(syst.) is obtained for the 7 TeV dataset.
The signal strength obtained for the combined 7 and 8 TeV dataset is 0.74 ± 0.09(stat.) ±0.14(syst.) The V Z signal is observed with a significance of 4.9σ, to be compared to an
expected significance of 6.3σ. The fitted µV Z values are shown in figure 26 for the 7 TeV,
8 TeV and combined datasets, and for the three lepton channels separately for the combined
dataset, all with the MVA used for the 8 TeV data. A measurement of V Z production in√pp collisions at s = 8 TeV in final states with b-tagged jets was recently reported by the
CMS Collaboration [99].

ATLASATLAS

s=7 TeV, ∫ Ldt=4.7 fb-1; s=8 TeV, ∫ Ldt=20.3 fb-1

VZ fittot.

tot ( stat

stat.

stat.

+ 0.15 )0.74 +− 0.17
( +− 0.09
0.09 − 0.13
0.16

Combination

0.5

1

1.5

2

2.5

syst )+ 0.21 )0.67 +− 0.25
( +− 0.14
0.14 − 0.20
0.24

1 lepton+ 0.33 )0.80 +− 0.39
( +− 0.20
0.20 − 0.29
0.35

0 lepton+ 0.16 )0.77 +− 0.22
( +− 0.15
0.15 − 0.14
0.21

Combination+ 0.15 )0.74 +− 0.17
( +− 0.09
0.09 − 0.13
0.16

3

best fit µ =σ/ σSM for mH=125 GeV
VZ

(a)

0

0.5

1

1.5

2

2.5

3

best fit µ =σ/ σSM for mH=125 GeV
VZ(b)Figure 26. The fitted values of the diboson signal strength µV Z for (a) the 7 TeV, 8 TeV andcombined datasets, and (b) for the three lepton channels separately and combined, for the combined
dataset. The MVA is used for the 8 TeV data. The individual µV Z values for the lepton channels
are obtained from a simultaneous fit with the signal strength for each floating independently.

12

Summary

A search for the Standard Model Higgs boson produced in association with a W or Z
boson and decaying into bb has been presented. The (W/Z) decay channels considered are
W → `ν, Z → `` and Z → νν. The dataset corresponds to integrated luminosities of
4.7 fb−1 and 20.3 fb−1 from pp collisions at 7 TeV and 8 TeV, respectively, recorded by the
ATLAS experiment during Run 1 of the LHC.
The analysis is carried out in event categories based on the numbers of leptons, jets, andjets tagged as originating from b-quark fragmentation, and on the transverse momentumof the vector-boson candidate. A multivariate analysis provides the nominal results. An
alternative analysis using invariant-mass distributions of the Higgs boson candidates leads
to consistent results.
For a Higgs boson mass of 125.36 GeV, the observed (expected) deviation from thebackground-only hypothesis corresponds to a significance of 1.4 (2.6) standard deviations
and the ratio of the measured signal yield to the Standard Model expectation is found to beµ = 0.52 ± 0.32(stat.) ± 0.24(syst.). The analysis procedure is validated by a measurementof the yield of (W/Z)Z production with Z → bb, from which the ratio of the observed signal
yield to the Standard Model expectation is found to be 0.74 ± 0.09(stat.) ± 0.14(syst.).

Acknowledgments
We thank CERN for the very successful operation of the LHC, as well as the support staff
from our institutions without whom ATLAS could not be operated efficiently.
We acknowledge the support of ANPCyT, Argentina; YerPhI, Armenia; ARC, Australia; BMWFW and FWF, Austria; ANAS, Azerbaijan; SSTC, Belarus; CNPq and
FAPESP, Brazil; NSERC, NRC and CFI, Canada; CERN; CONICYT, Chile; CAS, MOST

– 60 –JHEP01(2015)069

0

tot ( stat

2 lepton+ 0.16 )0.77 +− 0.19
( +− 0.10
0.09 − 0.14
0.17

8 TeV

VZ fittot.

syst )+ 0.39 )0.50 +− 0.49
( +− 0.30
0.29 − 0.37
0.47

7 TeV

s=7 TeV, ∫ Ldt=4.7 fb-1; s=8 TeV, ∫ Ldt=20.3 fb-1

A

Tables of event yields

The event yields in each category for the multivariate analysis are shown in tables 9–11.

– 61 –JHEP01(2015)069

and NSFC, China; COLCIENCIAS, Colombia; MSMT CR, MPO CR and VSC CR,
Czech Republic; DNRF, DNSRC and Lundbeck Foundation, Denmark; EPLANET, ERC
and NSRF, European Union; IN2P3-CNRS, CEA-DSM/IRFU, France; GNSF, Georgia;
BMBF, DFG, HGF, MPG and AvH Foundation, Germany; GSRT and NSRF, Greece;
ISF, MINERVA, GIF, I-CORE and Benoziyo Center, Israel; INFN, Italy; MEXT and
JSPS, Japan; CNRST, Morocco; FOM and NWO, Netherlands; BRF and RCN, Norway; MNiSW and NCN, Poland; GRICES and FCT, Portugal; MNE/IFA, Romania; MES
of Russia and ROSATOM, Russian Federation; JINR; MSTD, Serbia; MSSR, Slovakia;
ARRS and MIZŠ, Slovenia; DST/NRF, South Africa; MINECO, 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), CC-IN2P3 (France), KIT/GridKA (Germany), INFN-CNAF(Italy), NL-T1 (Netherlands), PIC (Spain), ASGC (Taiwan), RAL (U.K.) and BNL(U.S.A.) and in the Tier-2 facilities worldwide.

7.9
235
840
531
5470
2230
762
3890
1590
2550
18340 ± 150
18343
–
–
–
–
–
–
–
–
–
–
–
–

VH
W Z, ZZ, W W
tt̄
Single top
Wl
W cl
W +hf
Zl
Zcl
Z+hf
Total
Data

VH
W Z, ZZ, W W
tt̄
Single top
Wl
W cl
W +hf
Zl
Zcl
Z+hf
Total
Data

LL-tag100 − 120 GeV > 120 GeV
2-jet
2.2
6.9
18
49
114
183
40
56
159
206
106
159
54
124
96
272
59
162
225
607
886 ± 171841 ± 25
887
1860
3-jet
–
2
–
17
–
186
–
25
–
59
–
53
–
60
–
79
–
65
–
170
–
718 ± 12
–
657
–
–
–
–
–
–
–
–
–
–
–
–

3.5
14.7
129
32.6
14.4
23.4
33.6
6.1
9.4
186
458.5 ± 9.8477

MM-tag100 − 120 GeV
4.0
13.3
116
22.8
0.2
1.1
21.8
0.1
0.3
151
333.2 ± 8.9
306
–
–
–
–
–
–
–
–
–
–
–
–

7
20
315
30
4.3
8.9
62
4.5
9.8
259
719 ± 17710

TT-tag100 − 120 GeV

23
81
332
66
16
27
128
17
24
876
1599 ± 28
1592

MM+TT-tag> 120 GeV

Table 9. The expected signal and fitted background yields for each category of the 0-lepton channel after the full selection of the multivariate
analysis. The mbb distribution is used in the 100 < pVT < 120 GeV interval. The background yields are normalised by the results of the global
likelihood fit. All systematic uncertainties are included in the indicated uncertainties.

8
260
1670
318
2280
1240
750
3190
1620
1890
13310 ± 100
13344

23
635
1520
704
7100
3710
1520
10750
3990
6510
36890 ± 200
36903

1-tag
100 − 120 GeV> 120 GeV

SamplepVT

JHEP01(2015)069– 62 –

15
1100
18660
7390
24980
25900
6530
1780
690
1024
–
5300
93350 ± 320
93359

VH
WZ,ZZ,WWtt̄
Single top
W +l
W +cl
W +hf
Z+l
Z+cl
Z+hf
MJ (e)MJ (µ)
Total
Data

14
689
10490
2815
11320
10080
4740
572
193
272
1290
91
42570 ± 200
42557

28
927
6641
3730
23702
21650
6112
1223
333
475
2618
164
67600 ± 290
67603

> 120

LL-tagMM-tag< 120 GeV > 120 GeV < 120 GeV > 120 GeV
2-jet
11
9.3
17
14
103.6
62
64
36
1954
1051
2426
1080
1150
398
975
307
1603
697
124
48
2630
966
465
153
1470
648
954
402
101
37
7.4
2.8
53
18
10
3.1
161
45
126
30
–
162
–
89
721
16.0
329
4.8
9953 ± 91 4106 ± 50 5492 ± 66 2161 ± 33
9941
4072
5499
2199
3-jet
3.2
3.8
4.8
5.8
50
39.6
22.6
18
3240
1622
4119
1670
66
318
619
261
588
322
42
20
952
454
164
72
576
490
353
297
43
18.1
2.8
1.4
27
9.8
4.5
1.6
77
25.9
54
18.8
–
68.6
–
36
227
4.9
117
3.2
6447 ± 57 3376 ± 43 5501 ± 50 2405 ± 33
6336
3472
5551
2356
16
29.5
890
219
0.9
6.5
227
0.0
0.2
24
40
1.3
1448 ± 27
1405
6.5
14
1388
188
0.4
3.2
168
0.0
0.1
14
15
0.8
1796 ± 30
1838

19
52
2290
739
3.1
23
506
0.1
0.5
85
–
178
3889 ± 55
3923
5.4
14
4181
503
1.1
7.7
187
0.0
0.2
40
–
58
4995 ± 55
4977

TT-tag< 120 GeV > 120 GeV

Table 10. The expected signal and fitted background yields for each category of the 1-lepton channel after the full selection of the multivariate
analysis. The background yields are normalised by the results of the global likelihood fit. All systematic uncertainties are included in the indicateduncertainties.

39
1950
11380
13680
65980
71930
16030
3940
1350
2080
–
10230
198540 ± 500
1985441-tag< 120 GeV

VH
W Z, ZZ, W W
tt̄
Single top
Wl
W cl
W +hf
Zl
Zcl
Z+hf
MJ (e)MJ (µ)
Total
Data

SamplepVT

JHEP01(2015)069– 63 –

20
800
3490
385
245
166
90
29410
12130
21090
69360 ± 280
69313
7.4
347
2170
144
76
57
40
8870
4650
5790
22720 ± 150
22662

VH
WZ,ZZ,WWtt̄
Single top
W +l
W +cl
W +hf
Z+l
Z+cl
Z+hf
Total
Data

VH
WZ,ZZ,WWtt̄
Single top
W +l
W +cl
W +hf
Z+l
Z+cl
Z+hf
Total
Data

LL-tag< 120 GeV > 120 GeV
2-jet
5.5
2.0
73
12.2
1023
49
40
5.8
3.7
0.4
5.8
0.6
7.5
0.7
648
101
421
64
1916
287.7
4271 ± 50 524.6 ± 7.9
4313
517
3-jet
1.5
0.9
22
6.9
486
29
18
3.0
1.4
0.2
2.3
0.4
3.1
0.8
190
46
164
38
497
121
1403 ± 20 245.2 ± 5.2
1428
253
4.9
29
1476
41
0.1
0.4
0.8
12
25
766
2365 ± 43
2394

19
129
3180
96
0.3
1.0
4.9
41
71
2858
6473 ± 76
6501
2.8
8.7
74
6.5
0.0
0.1
0.4
2.7
5.7
180
280.9 ± 8.4
283

7.0
22
137
12
0.0
0.1
0.6
6.2
9.7
414
607 ± 13570

MM+TT-tag< 120 GeV > 120 GeV

Table 11. The expected signal and fitted background yields for each category of the 2-lepton channel after the full selection of the multivariate
analysis. The background yields are normalised by the results of the global likelihood fit. All systematic uncertainties are included in the indicateduncertainties.

3.9
100
182
22
7.1
8.9
6.4
1913
949
1223
4430 ± 56
4436

7.4
135
213
44
12
18
13
4180
1622
2853
9123 ± 86
9150

1-tag
< 120 GeV
> 120 GeV

SamplepVT

JHEP01(2015)069– 64 –

Open Access. This article is distributed under the terms of the Creative Commons
Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction inany medium, provided the original author(s) and source are credited.

References[1] F. Englert and R. Brout, Broken Symmetry and the Mass of Gauge Vector Mesons, Phys.
Rev. Lett. 13 (1964) 321 [INSPIRE].

[3] P.W. Higgs, Broken symmetries, massless particles and gauge fields, Phys. Lett. 12 (1964)
132 [INSPIRE].
[4] G.S. Guralnik, C.R. Hagen and T.W.B. Kibble, Global Conservation Laws and MasslessParticles, Phys. Rev. Lett. 13 (1964) 585 [INSPIRE].
[5] ATLAS collaboration, Observation of a new particle in the search for the Standard Model
Higgs boson with the ATLAS detector at the LHC, Phys. Lett. B 716 (2012) 1[arXiv:1207.7214] [INSPIRE].
[6] CMS collaboration, Observation of a new boson at a mass of 125 GeV with the CMS
experiment at the LHC, Phys. Lett. B 716 (2012) 30 [arXiv:1207.7235] [INSPIRE].
[7] ATLAS collaboration, Measurements of Higgs boson production and couplings in diboson
final states with the ATLAS detector at the LHC, Phys. Lett. B 726 (2013) 88[arXiv:1307.1427] [INSPIRE].
[8] ATLAS collaboration, Evidence for the spin-0 nature of the Higgs boson using ATLAS data,
Phys. Lett. B 726 (2013) 120 [arXiv:1307.1432] [INSPIRE].
[9] CMS collaboration, Observation of a new boson with mass near 125 GeV in pp collisions at√
s = 7 and 8 TeV, JHEP 06 (2013) 081 [arXiv:1303.4571] [INSPIRE].
[10] CMS collaboration, Evidence for the 125 GeV Higgs boson decaying to a pair of τ leptons,
JHEP 05 (2014) 104 [arXiv:1401.5041] [INSPIRE].
[11] A. Djouadi, J. Kalinowski and M. Spira, HDECAY: A program for Higgs boson decays in thestandard model and its supersymmetric extension, Comput. Phys. Commun. 108 (1998) 56[hep-ph/9704448] [INSPIRE].
[12] R. Lafaye, T. Plehn, M. Rauch, D. Zerwas and M. Dührssen, Measuring the Higgs Sector,
JHEP 08 (2009) 009 [arXiv:0904.3866] [INSPIRE].
[13] S.L. Glashow, D.V. Nanopoulos and A. Yildiz, Associated Production of Higgs Bosons and Z
Particles, Phys. Rev. D 18 (1978) 1724 [INSPIRE].
[14] A. Stange, W.J. Marciano and S. Willenbrock, Higgs bosons at the Fermilab Tevatron, Phys.
Rev. D 49 (1994) 1354 [hep-ph/9309294] [INSPIRE].
[15] A. Stange, W.J. Marciano and S. Willenbrock, Associated production of Higgs and weakbosons, with H → bb, at hadron colliders, Phys. Rev. D 50 (1994) 4491 [hep-ph/9404247][INSPIRE].
[16] CDF and D0 collaborations, T. Aaltonen et al., Evidence for a particle produced inassociation with weak bosons and decaying to a bottom-antibottom quark pair in Higgs boson
searches at the Tevatron, Phys. Rev. Lett. 109 (2012) 071804 [arXiv:1207.6436] [INSPIRE].

– 65 –JHEP01(2015)069[2] P.W. Higgs, Broken Symmetries and the Masses of Gauge Bosons, Phys. Rev. Lett. 13(1964) 508 [INSPIRE].

[17] CMS collaboration, Search for the standard model Higgs boson produced in association witha W or a Z boson and decaying to bottom quarks, Phys. Rev. D 89 (2014) 012003[arXiv:1310.3687] [INSPIRE].
[18] 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, Phys. Lett. B
718 (2012) 369 [arXiv:1207.0210] [INSPIRE].

[21] ATLAS collaboration, The ATLAS Simulation Infrastructure, Eur. Phys. J. C 70 (2010)
823 [arXiv:1005.4568] [INSPIRE].
[22] GEANT4 collaboration, S. Agostinelli et al., GEANT4: A simulation toolkit, Nucl. Instrum.
Meth. A 506 (2003) 250 [INSPIRE].
[23] T. Sjöstrand, S. Mrenna and P.Z. Skands, A Brief Introduction to PYTHIA 8.1, Comput.
Phys. Commun. 178 (2008) 852 [arXiv:0710.3820] [INSPIRE].
[24] J. Pumplin et al., New generation of parton distributions with uncertainties from global QCD
analysis, JHEP 07 (2002) 012 [hep-ph/0201195] [INSPIRE].
[25] ATLAS collaboration, ATLAS tunes of PYTHIA6 and PYTHIA8 for MC11,
ATL-PHYS-PUB-2011-009 (2011).
[26] ATLAS collaboration, New ATLAS event generator tunes to 2010 data,
ATL-PHYS-PUB-2011-008 (2011).
[27] P. Golonka and Z. Was, PHOTOS Monte Carlo: A Precision tool for QED corrections in Z
and W decays, Eur. Phys. J. C 45 (2006) 97 [hep-ph/0506026] [INSPIRE].
[28] P. Nason, A new method for combining NLO QCD with shower Monte Carlo algorithms,
JHEP 11 (2004) 040 [hep-ph/0409146] [INSPIRE].
[29] S. Frixione, P. Nason and C. Oleari, Matching NLO QCD computations with Parton Shower
simulations: the POWHEG method, JHEP 11 (2007) 070 [arXiv:0709.2092] [INSPIRE].
[30] S. Alioli, P. Nason, C. Oleari and E. Re, A general framework for implementing NLO
calculations in shower Monte Carlo programs: the POWHEG BOX, JHEP 06 (2010) 043[arXiv:1002.2581] [INSPIRE].
[31] G. Luisoni, P. Nason, C. Oleari and F. Tramontano, HW ± /HZ + 0 and 1 jet at NLO with
the POWHEG BOX interfaced to GoSam and their merging within MiNLO, JHEP 10 (2013)083 [arXiv:1306.2542] [INSPIRE].
[32] H.-L. Lai et al., New parton distributions for collider physics, Phys. Rev. D 82 (2010) 074024[arXiv:1007.2241] [INSPIRE].
[33] C. Englert, M. McCullough and M. Spannowsky, Gluon-initiated associated production boosts
Higgs physics, Phys. Rev. D 89 (2014) 013013 [arXiv:1310.4828] [INSPIRE].
[34] J. Ohnemus and W.J. Stirling, Order αs corrections to the differential cross-section for theW H intermediate mass Higgs signal, Phys. Rev. D 47 (1993) 2722 [INSPIRE].

– 66 –JHEP01(2015)069[19] ATLAS collaboration, The ATLAS Experiment at the CERN Large Hadron Collider, 2008
JINST 3 S08003 [INSPIRE].
√[20] ATLAS collaboration, Improved luminosity determination in pp collisions at s = 7 TeV
using the ATLAS detector at the LHC, Eur. Phys. J. C 73 (2013) 2518 [arXiv:1302.4393][INSPIRE].

[35] H. Baer, B. Bailey and J.F. Owens, O(αs ) Monte Carlo approach to W+ Higgs associatedproduction at hadron supercolliders, Phys. Rev. D 47 (1993) 2730 [INSPIRE].
[36] O. Brein, A. Djouadi and R. Harlander, NNLO QCD corrections to the Higgs-strahlung
processes at hadron colliders, Phys. Lett. B 579 (2004) 149 [hep-ph/0307206] [INSPIRE].
[37] M.L. Ciccolini, S. Dittmaier and M. Krämer, Electroweak radiative corrections to associatedWH and ZH production at hadron colliders, Phys. Rev. D 68 (2003) 073003[hep-ph/0306234] [INSPIRE].
[38] A. Denner, S. Dittmaier, S. Kallweit and A. Mück, EW corrections to Higgs strahlung at theTevatron and the LHC with HAWK, PoS(EPS-HEP2011)235 [arXiv:1112.5258] [INSPIRE].

[40] T. Gleisberg et al., Event generation with SHERPA 1.1, JHEP 02 (2009) 007[arXiv:0811.4622] [INSPIRE].
[41] T. Sjöstrand, S. Mrenna and P.Z. Skands, PYTHIA 6.4 Physics and Manual, JHEP 05(2006) 026 [hep-ph/0603175] [INSPIRE].
[42] K. Melnikov and F. Petriello, Electroweak gauge boson production at hadron colliders throughO(αs2 ), Phys. Rev. D 74 (2006) 114017 [hep-ph/0609070] [INSPIRE].
[43] M. Czakon, P. Fiedler and A. Mitov, Total Top-Quark Pair-Production Cross section at
Hadron Colliders Through O(αS4 ), Phys. Rev. Lett. 110 (2013) 252004 [arXiv:1303.6254][INSPIRE].
[44] B.P. Kersevan and E. Richter-Was, The Monte Carlo event generator AcerMC versions 2.0
to 3.8 with interfaces to PYTHIA 6.4, HERWIG 6.5 and ARIADNE 4.1, Comput. Phys.
Commun. 184 (2013) 919 [hep-ph/0405247] [INSPIRE].
[45] N. Kidonakis, Next-to-next-to-leading-order collinear and soft gluon corrections for t-channel
single top quark production, Phys. Rev. D 83 (2011) 091503 [arXiv:1103.2792] [INSPIRE].
[46] N. Kidonakis, NNLL resummation for s-channel single top quark production, Phys. Rev. D
81 (2010) 054028 [arXiv:1001.5034] [INSPIRE].
[47] N. Kidonakis, Two-loop soft anomalous dimensions for single top quark associated production
with a W− or H− , Phys. Rev. D 82 (2010) 054018 [arXiv:1005.4451] [INSPIRE].
[48] P. Nason and G. Zanderighi, W + W − , W Z and ZZ production in the POWHEG-BOX-V2,
Eur. Phys. J. C 74 (2014) 2702 [arXiv:1311.1365] [INSPIRE].
[49] G. Corcella et al., HERWIG 6: An Event generator for hadron emission reactions withinterfering gluons (including supersymmetric processes), JHEP 01 (2001) 010[hep-ph/0011363] [INSPIRE].
[50] J.M. Campbell and R.K. Ellis, MCFM for the Tevatron and the LHC, Nucl. Phys. Proc.
Suppl. 205-206 (2010) 10 [arXiv:1007.3492] [INSPIRE].
[51] A.D. Martin, W.J. Stirling, R.S. Thorne and G. Watt, Parton distributions for the LHC,
Eur. Phys. J. C 63 (2009) 189 [arXiv:0901.0002] [INSPIRE].
[52] A. Sherstnev and R.S. Thorne, Parton Distributions for LO Generators, Eur. Phys. J. C 55
(2008) 553 [arXiv:0711.2473] [INSPIRE].

– 67 –JHEP01(2015)069[39] L. Altenkamp, S. Dittmaier, R.V. Harlander, H. Rzehak and T.J.E. Zirke, Gluon-induced
Higgs-strahlung at next-to-leading order QCD, JHEP 02 (2013) 078 [arXiv:1211.5015][INSPIRE].

[53] ATLAS collaboration, Electron reconstruction and identification efficiency measurementswith the ATLAS detector using the 2011 LHC proton-proton collision data, Eur. Phys. J. C
74 (2014) 2941 [arXiv:1404.2240] [INSPIRE].
[54] ATLAS collaboration, Electron efficiency measurements with the ATLAS detector using the2012 LHC proton-proton collision data, ATLAS-CONF-2014-032 (2014).
[55] ATLAS collaboration, Muon reconstruction efficiency and momentum resolution of the√
ATLAS experiment in proton-proton collisions at s = 7 TeV in 2010, Eur. Phys. J. C 74
(2014) 3034 [arXiv:1404.4562] [INSPIRE].

[57] W. Lampl et al., Calorimeter clustering algorithms: Description and performance,
ATL-LARG-PUB-2008-002 (2008).
[58] M. Cacciari, G.P. Salam and G. Soyez, The Anti-kt jet clustering algorithm, JHEP 04 (2008)
063 [arXiv:0802.1189] [INSPIRE].
[59] M. Cacciari and G.P. Salam, Pileup subtraction using jet areas, Phys. Lett. B 659 (2008) 119[arXiv:0707.1378] [INSPIRE].
[60] ATLAS collaboration, Jet energy measurement with the ATLAS detector in proton-proton√
collisions at s = 7 TeV, Eur. Phys. J. C 73 (2013) 2304 [arXiv:1112.6426] [INSPIRE].
[61] ATLAS collaboration, Jet energy measurement and its systematic uncertainty in√proton-proton collisions at s = 7 TeV with the ATLAS detector, arXiv:1406.0076[INSPIRE].
[62] ATLAS collaboration, Commissioning of the ATLAS high-performance b-tagging algorithms
in the 7 TeV collision data, ATLAS-CONF-2011-102 (2011).
[63] ATLAS collaboration, Calibration of b-tagging using dileptonic top pair events in a
combinatorial likelihood approach with the ATLAS experiment, ATLAS-CONF-2014-004(2014).
[64] ATLAS collaboration, Calibration of the performance of b-tagging for c and light-flavour jetsin the 2012 ATLAS data, ATLAS-CONF-2014-046 (2014).
[65] ATLAS collaboration, Performance of Missing Transverse Momentum Reconstruction in
Proton-Proton Collisions at 7 TeV with ATLAS, Eur. Phys. J. C 72 (2012) 1844
[arXiv:1108.5602] [INSPIRE].
[66] ATLAS collaboration, Performance of Missing Transverse Momentum Reconstruction in
ATLAS studied in Proton-Proton Collisions recorded in 2012 at 8 TeV,ATLAS-CONF-2013-082 (2013).
[67] J. Gallicchio et al., Multivariate discrimination and the Higgs + W/Z search, JHEP 04(2011) 069 [arXiv:1010.3698] [INSPIRE].
[68] A. Bukin, Fitting function for asymmetric peaks, arXiv:0711.4449.
[69] L. Breiman, J. Friedman, R. Olshen and C. Stone, Classification and Regression Trees,
Wadsworth and Brooks, Monterey, CA U.S.A. (1984).
[70] Y. Freund and R. Schapire, Experiments with a new boosting algorithm, in Machine Learning,
Proceedings of the Thirteenth International Conference (ICML), M. Kaufmann, (1996).

– 68 –JHEP01(2015)069[56] ATLAS collaboration, Electron and photon energy calibration with the ATLAS detector
using LHC Run 1 data, Eur. Phys. J. C 74 (2014) 3071 [arXiv:1407.5063] [INSPIRE].

[71] A. Höcker et al., TMVA: Toolkit for Multivariate Data Analysis, PoS(ACAT)040[physics/0703039] [INSPIRE].
[72] S. Höche, F. Krauss, M. Schönherr and F. Siegert, QCD matrix elements + parton showers:
The NLO case, JHEP 04 (2013) 027 [arXiv:1207.5030] [INSPIRE].

[75] H1 and ZEUS collaborations, F.D. Aaron et al., Combined Measurement and QCD Analysisof the Inclusive e± p Scattering Cross sections at HERA, JHEP 01 (2010) 109[arXiv:0911.0884] [INSPIRE].
[76] S. Frixione and B.R. Webber, Matching NLO QCD computations and parton shower
simulations, JHEP 06 (2002) 029 [hep-ph/0204244] [INSPIRE].
[77] M.L. Mangano, M. Moretti, F. Piccinini, R. Pittau and A.D. Polosa, ALPGEN, a generator
for hard multiparton processes in hadronic collisions, JHEP 07 (2003) 001 [hep-ph/0206293][INSPIRE].
[78] N. Kidonakis, Differential and total cross sections for top pair and single top production,
arXiv:1205.3453 [INSPIRE].
[79] J. Alwall, M. Herquet, F. Maltoni, O. Mattelaer and T. Stelzer, MadGraph 5: Going Beyond,
JHEP 06 (2011) 128 [arXiv:1106.0522] [INSPIRE].
[80] J. Alwall et al., The automated computation of tree-level and next-to-leading orderdifferential cross sections and their matching to parton shower simulations, JHEP 07 (2014)079 [arXiv:1405.0301] [INSPIRE].
[81] V. Hirschi et al., Automation of one-loop QCD corrections, JHEP 05 (2011) 044[arXiv:1103.0621] [INSPIRE].
[82] R. Frederix et al., W and Z/γ ∗ boson production in association with a bottom-antibottom
pair, JHEP 09 (2011) 061 [arXiv:1106.6019] [INSPIRE].
[83] S. Frixione, E. Laenen, P. Motylinski, B.R. Webber and C.D. White, Single-top
hadroproduction in association with a W boson, JHEP 07 (2008) 029 [arXiv:0805.3067][INSPIRE].
[84] K. Arnold et al., HERWIG++ 2.6 Release Note, arXiv:1205.4902 [INSPIRE].
[85] R.D. Ball et al., Parton distributions with LHC data, Nucl. Phys. B 867 (2013) 244[arXiv:1207.1303] [INSPIRE].
[86] I.W. Stewart and F.J. Tackmann, Theory Uncertainties for Higgs and Other Searches Using
Jet Bins, Phys. Rev. D 85 (2012) 034011 [arXiv:1107.2117] [INSPIRE].
[87] M. Botje et al., The PDF4LHC Working Group Interim Recommendations,
arXiv:1101.0538 [INSPIRE].
[88] LHC Higgs Cross section Working Group collaboration, S. Dittmaier et al., Handbookof LHC Higgs Cross sections: 1. Inclusive Observables, arXiv:1101.0593 [INSPIRE].
[89] L. Moneta et al., The RooStats Project, PoS(ACAT)057 [arXiv:1009.1003] [INSPIRE].

– 69 –JHEP01(2015)069[73] ATLAS collaboration, Measurements of normalized differential cross sections for tt̄
√production in pp collisions at s = 7 TeV using the ATLAS detector, Phys. Rev. D 90
(2014) 072004 [arXiv:1407.0371] [INSPIRE].
√[74] ATLAS collaboration, Jet energy resolution in proton-proton collisions at s = 7 TeV
recorded in 2010 with the ATLAS detector, Eur. Phys. J. C 73 (2013) 2306
[arXiv:1210.6210] [INSPIRE].

[90] W. Verkerke and D.P. Kirkby, The RooFit toolkit for data modeling, eConf C 0303241(2003) MOLT007 [physics/0306116] [INSPIRE].
[91] A.L. Read, Presentation of search results: The CLs technique, J. Phys. G 28 (2002) 2693
[INSPIRE].
[92] G. Cowan, K. Cranmer, E. Gross and O. Vitells, Asymptotic formulae for likelihood-based
tests of new physics, Eur. Phys. J. C 71 (2011) 1554 [Erratum ibid. C 73 (2013) 2501][arXiv:1007.1727] [INSPIRE].
[93] G. Bohm and G. Zech, Introduction to statistics and data analysis for physicts, Verlag
Deutsches Elektronen-Synchrotron, Hamburg Germany (2010).

[95] 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).
[96] ATLAS collaboration, b-jet tagging calibration on c-jets containing D∗+ mesons,
ATLAS-CONF-2012-039 (2012).
[97] ATLAS collaboration, Measurement of the Mistag Rate with 5 fb−1 of Data Collected by the
ATLAS Detector, ATLAS-CONF-2012-040 (2012).
[98] ATLAS collaboration, Measurement of the Higgs boson mass from the H → γγ and
H → ZZ ∗ → 4` channels with the ATLAS detector using 25 fb−1 of pp collision data, Phys.
Rev. D 90 (2014) 052004 [arXiv:1406.3827] [INSPIRE].
√[99] CMS collaboration, Measurement of WZ and ZZ production in pp collisions at s = 8 TeV
in final states with b-tagged jets, Eur. Phys. J. C 74 (2014) 2973 [arXiv:1403.3047][INSPIRE].

– 70 –JHEP01(2015)069[94] ATLAS collaboration, Measuring the b-tag efficiency in a top-pair sample with 4.7 fb−1 ofdata from the ATLAS detector, ATLAS-CONF-2012-097 (2012).

The ATLAS collaboration– 71 –JHEP01(2015)069

G. Aad85 , B. Abbott113 , J. Abdallah152 , S. Abdel Khalek117 , O. Abdinov11 , R. Aben107 ,B. Abi114 , M. Abolins90 , O.S. AbouZeid159 , H. Abramowicz154 , H. Abreu153 , R. Abreu30 ,Y. Abulaiti147a,147b , B.S. Acharya165a,165b,a , L. Adamczyk38a , D.L. Adams25 , J. Adelman177 ,S. Adomeit100 , T. Adye131 , T. Agatonovic-Jovin13a , J.A. Aguilar-Saavedra126a,126f ,M. Agustoni17 , S.P. Ahlen22 , F. Ahmadov65,b , G. Aielli134a,134b , H. Akerstedt147a,147b ,T.P.A. Åkesson81 , G. Akimoto156 , A.V. Akimov96 , G.L. Alberghi20a,20b , J. Albert170 ,S. Albrand55 , M.J. Alconada Verzini71 , M. Aleksa30 , I.N. Aleksandrov65 , C. Alexa26a ,G. Alexander154 , G. Alexandre49 , T. Alexopoulos10 , M. Alhroob113 , G. Alimonti91a , L. Alio85 ,J. Alison31 , B.M.M. Allbrooke18 , L.J. Allison72 , P.P. Allport74 , A. Aloisio104a,104b , A. Alonso36 ,F. Alonso71 , C. Alpigiani76 , A. Altheimer35 , B. Alvarez Gonzalez90 , M.G. Alviggi104a,104b ,K. Amako66 , Y. Amaral Coutinho24a , C. Amelung23 , D. Amidei89 , S.P. Amor Dos Santos126a,126c ,A. Amorim126a,126b , S. Amoroso48 , N. Amram154 , G. Amundsen23 , C. Anastopoulos140 ,L.S. Ancu49 , N. Andari30 , T. Andeen35 , C.F. Anders58b , G. Anders30 , K.J. Anderson31 ,A. Andreazza91a,91b , V. Andrei58a , X.S. Anduaga71 , S. Angelidakis9 , I. Angelozzi107 , P. Anger44 ,A. Angerami35 , F. Anghinolfi30 , A.V. Anisenkov109,c , N. Anjos12 , A. Annovi47 , A. Antonaki9 ,M. Antonelli47 , A. Antonov98 , J. Antos145b , F. Anulli133a , M. Aoki66 , L. Aperio Bella18 ,R. Apolle120,d , G. Arabidze90 , I. Aracena144 , Y. Arai66 , J.P. Araque126a , A.T.H. Arce45 ,F.A. Arduh71 , J-F. Arguin95 , S. Argyropoulos42 , M. Arik19a , A.J. Armbruster30 , O. Arnaez30 ,V. Arnal82 , H. Arnold48 , M. Arratia28 , O. Arslan21 , A. Artamonov97 , G. Artoni23 , S. Asai156 ,N. Asbah42 , A. Ashkenazi154 , B. Åsman147a,147b , L. Asquith6 , K. Assamagan25 , R. Astalos145a ,M. Atkinson166 , N.B. Atlay142 , B. Auerbach6 , K. Augsten128 , M. Aurousseau146b , G. Avolio30 ,B. Axen15 , G. Azuelos95,e , Y. Azuma156 , M.A. Baak30 , A.E. Baas58a , C. Bacci135a,135b ,H. Bachacou137 , K. Bachas155 , M. Backes30 , M. Backhaus30 , J. Backus Mayes144 , E. Badescu26a ,P. Bagiacchi133a,133b , P. Bagnaia133a,133b , Y. Bai33a , T. Bain35 , J.T. Baines131 , O.K. Baker177 ,P. Balek129 , F. Balli137 , E. Banas39 , Sw. Banerjee174 , A.A.E. Bannoura176 , V. Bansal170 ,H.S. Bansil18 , L. Barak173 , S.P. Baranov96 , E.L. Barberio88 , D. Barberis50a,50b , M. Barbero85 ,T. Barillari101 , M. Barisonzi176 , T. Barklow144 , N. Barlow28 , S.L. Barnes84 , B.M. Barnett131 ,R.M. Barnett15 , Z. Barnovska5 , A. Baroncelli135a , G. Barone49 , A.J. Barr120 , F. Barreiro82 ,J. Barreiro Guimarães da Costa57 , R. Bartoldus144 , A.E. Barton72 , P. Bartos145a , V. Bartsch150 ,A. Bassalat117 , A. Basye166 , R.L. Bates53 , J.R. Batley28 , M. Battaglia138 , M. Battistin30 ,F. Bauer137 , H.S. Bawa144,f , M.D. Beattie72 , T. Beau80 , P.H. Beauchemin162 ,R. Beccherle124a,124b , P. Bechtle21 , H.P. Beck17 , K. Becker176 , S. Becker100 , M. Beckingham171 ,C. Becot117 , A.J. Beddall19c , A. Beddall19c , S. Bedikian177 , V.A. Bednyakov65 , C.P. Bee149 ,L.J. Beemster107 , T.A. Beermann176 , M. Begel25 , K. Behr120 , C. Belanger-Champagne87 ,P.J. Bell49 , W.H. Bell49 , G. Bella154 , L. Bellagamba20a , A. Bellerive29 , M. Bellomo86 ,K. Belotskiy98 , O. Beltramello30 , O. Benary154 , D. Benchekroun136a , K. Bendtz147a,147b ,N. Benekos166 , Y. Benhammou154 , E. Benhar Noccioli49 , J.A. Benitez Garcia160b ,D.P. Benjamin45 , J.R. Bensinger23 , S. Bentvelsen107 , D. Berge107 , E. Bergeaas Kuutmann167 ,N. Berger5 , F. Berghaus170 , J. Beringer15 , C. Bernard22 , P. Bernat78 , C. Bernius79 ,F.U. Bernlochner170 , T. Berry77 , P. Berta129 , C. Bertella85 , G. Bertoli147a,147b ,F. Bertolucci124a,124b , C. Bertsche113 , D. Bertsche113 , M.I. Besana91a , G.J. Besjes106 ,O. Bessidskaia147a,147b , M. Bessner42 , N. Besson137 , C. Betancourt48 , S. Bethke101 , W. Bhimji46 ,R.M. Bianchi125 , L. Bianchini23 , M. Bianco30 , O. Biebel100 , S.P. Bieniek78 , K. Bierwagen54 ,J. Biesiada15 , M. Biglietti135a , J. Bilbao De Mendizabal49 , H. Bilokon47 , M. Bindi54 , S. Binet117 ,A. Bingul19c , C. Bini133a,133b , C.W. Black151 , J.E. Black144 , K.M. Black22 , D. Blackburn139 ,R.E. Blair6 , J.-B. Blanchard137 , T. Blazek145a , I. Bloch42 , C. Blocker23 , W. Blum83,∗ ,– 72 –JHEP01(2015)069

U. Blumenschein54 , G.J. Bobbink107 , V.S. Bobrovnikov109,c , S.S. Bocchetta81 , A. Bocci45 ,C. Bock100 , C.R. Boddy120 , M. Boehler48 , T.T. Boek176 , J.A. Bogaerts30 , A.G. Bogdanchikov109 ,A. Bogouch92,∗ , C. Bohm147a , J. Bohm127 , V. Boisvert77 , T. Bold38a , V. Boldea26a ,A.S. Boldyrev99 , M. Bomben80 , M. Bona76 , M. Boonekamp137 , A. Borisov130 , G. Borissov72 ,M. Borri84 , S. Borroni42 , J. Bortfeldt100 , V. Bortolotto60a , K. Bos107 , D. Boscherini20a ,M. Bosman12 , H. Boterenbrood107 , J. Boudreau125 , J. Bouffard2 , E.V. Bouhova-Thacker72 ,D. Boumediene34 , C. Bourdarios117 , N. Bousson114 , S. Boutouil136d , A. Boveia31 , J. Boyd30 ,I.R. Boyko65 , I. Bozic13a , J. Bracinik18 , A. Brandt8 , G. Brandt15 , O. Brandt58a , U. Bratzler157 ,B. Brau86 , J.E. Brau116 , H.M. Braun176,∗ , S.F. Brazzale165a,165c , B. Brelier159 , K. Brendlinger122 ,A.J. Brennan88 , R. Brenner167 , S. Bressler173 , K. Bristow146c , T.M. Bristow46 , D. Britton53 ,F.M. Brochu28 , I. Brock21 , R. Brock90 , C. Bromberg90 , J. Bronner101 , G. Brooijmans35 ,T. Brooks77 , W.K. Brooks32b , J. Brosamer15 , E. Brost116 , J. Brown55 ,P.A. Bruckman de Renstrom39 , D. Bruncko145b , R. Bruneliere48 , S. Brunet61 , A. Bruni20a ,G. Bruni20a , M. Bruschi20a , L. Bryngemark81 , T. Buanes14 , Q. Buat143 , F. Bucci49 ,P. Buchholz142 , R.M. Buckingham120 , A.G. Buckley53 , S.I. Buda26a , I.A. Budagov65 ,F. Buehrer48 , L. Bugge119 , M.K. Bugge119 , O. Bulekov98 , A.C. Bundock74 , H. Burckhart30 ,S. Burdin74 , B. Burghgrave108 , S. Burke131 , I. Burmeister43 , E. Busato34 , D. Büscher48 ,V. Büscher83 , P. Bussey53 , C.P. Buszello167 , B. Butler57 , J.M. Butler22 , A.I. Butt3 ,C.M. Buttar53 , J.M. Butterworth78 , P. Butti107 , W. Buttinger28 , A. Buzatu53 , M. Byszewski10 ,S. Cabrera Urbán168 , D. Caforio20a,20b , O. Cakir4a , P. Calafiura15 , A. Calandri137 , G. Calderini80 ,P. Calfayan100 , R. Calkins108 , L.P. Caloba24a , D. Calvet34 , S. Calvet34 , R. Camacho Toro49 ,S. Camarda42 , D. Cameron119 , L.M. Caminada15 , R. Caminal Armadans12 , S. Campana30 ,M. Campanelli78 , A. Campoverde149 , V. Canale104a,104b , A. Canepa160a , M. Cano Bret76 ,J. Cantero82 , R. Cantrill126a , T. Cao40 , M.D.M. Capeans Garrido30 , I. Caprini26a , M. Caprini26a ,M. Capua37a,37b , R. Caputo83 , R. Cardarelli134a , T. Carli30 , G. Carlino104a , L. Carminati91a,91b ,S. Caron106 , E. Carquin32a , G.D. Carrillo-Montoya146c , J.R. Carter28 , J. Carvalho126a,126c ,D. Casadei78 , M.P. Casado12 , M. Casolino12 , E. Castaneda-Miranda146b , A. Castelli107 ,V. Castillo Gimenez168 , N.F. Castro126a , P. Catastini57 , A. Catinaccio30 , J.R. Catmore119 ,A. Cattai30 , G. Cattani134a,134b , J. Caudron83 , V. Cavaliere166 , D. Cavalli91a , M. Cavalli-Sforza12 ,V. Cavasinni124a,124b , F. Ceradini135a,135b , B.C. Cerio45 , K. Cerny129 , A.S. Cerqueira24b ,A. Cerri150 , L. Cerrito76 , F. Cerutti15 , M. Cerv30 , A. Cervelli17 , S.A. Cetin19b , A. Chafaq136a ,D. Chakraborty108 , I. Chalupkova129 , P. Chang166 , B. Chapleau87 , J.D. Chapman28 ,D. Charfeddine117 , D.G. Charlton18 , C.C. Chau159 , C.A. Chavez Barajas150 , S. Cheatham87 ,A. Chegwidden90 , S. Chekanov6 , S.V. Chekulaev160a , G.A. Chelkov65,g , M.A. Chelstowska89 ,C. Chen64 , H. Chen25 , K. Chen149 , L. Chen33d,h , S. Chen33c , X. Chen33f , Y. Chen67 , Y. Chen35 ,H.C. Cheng89 , Y. Cheng31 , A. Cheplakov65 , R. Cherkaoui El Moursli136e , V. Chernyatin25,∗ ,E. Cheu7 , L. Chevalier137 , V. Chiarella47 , G. Chiefari104a,104b , J.T. Childers6 , A. Chilingarov72 ,G. Chiodini73a , A.S. Chisholm18 , R.T. Chislett78 , A. Chitan26a , M.V. Chizhov65 , S. Chouridou9 ,B.K.B. Chow100 , D. Chromek-Burckhart30 , M.L. Chu152 , J. Chudoba127 , J.J. Chwastowski39 ,L. Chytka115 , G. Ciapetti133a,133b , A.K. Ciftci4a , R. Ciftci4a , D. Cinca53 , V. Cindro75 , A. Ciocio15 ,Z.H. Citron173 , M. Citterio91a , M. Ciubancan26a , A. Clark49 , P.J. Clark46 , R.N. Clarke15 ,W. Cleland125 , J.C. Clemens85 , C. Clement147a,147b , Y. Coadou85 , M. Cobal165a,165c ,A. Coccaro139 , J. Cochran64 , L. Coffey23 , J.G. Cogan144 , B. Cole35 , S. Cole108 , A.P. Colijn107 ,J. Collot55 , T. Colombo58c , G. Compostella101 , P. Conde Muiño126a,126b , E. Coniavitis48 ,S.H. Connell146b , I.A. Connelly77 , S.M. Consonni91a,91b , V. Consorti48 , S. Constantinescu26a ,C. Conta121a,121b , G. Conti57 , F. Conventi104a,i , M. Cooke15 , B.D. Cooper78 ,A.M. Cooper-Sarkar120 , N.J. Cooper-Smith77 , K. Copic15 , T. Cornelissen176 , M. Corradi20a ,F. Corriveau87,j , A. Corso-Radu164 , A. Cortes-Gonzalez12 , G. Cortiana101 , G. Costa91a ,– 73 –JHEP01(2015)069

M.J. Costa168 , D. Costanzo140 , D. Côté8 , G. Cottin28 , G. Cowan77 , B.E. Cox84 , K. Cranmer110 ,G. Cree29 , S. Crépé-Renaudin55 , F. Crescioli80 , W.A. Cribbs147a,147b , M. Crispin Ortuzar120 ,M. Cristinziani21 , V. Croft106 , G. Crosetti37a,37b , C.-M. Cuciuc26a , T. Cuhadar Donszelmann140 ,J. Cummings177 , M. Curatolo47 , C. Cuthbert151 , H. Czirr142 , P. Czodrowski3 , S. D’Auria53 ,M. D’Onofrio74 , M.J. Da Cunha Sargedas De Sousa126a,126b , C. Da Via84 , W. Dabrowski38a ,A. Dafinca120 , T. Dai89 , O. Dale14 , F. Dallaire95 , C. Dallapiccola86 , M. Dam36 , A.C. Daniells18 ,M. Dano Hoffmann137 , V. Dao48 , G. Darbo50a , S. Darmora8 , J.A. Dassoulas42 , A. Dattagupta61 ,W. Davey21 , C. David170 , T. Davidek129 , E. Davies120,d , M. Davies154 , O. Davignon80 ,A.R. Davison78 , P. Davison78 , Y. Davygora58a , E. Dawe143 , I. Dawson140 ,R.K. Daya-Ishmukhametova86 , K. De8 , R. de Asmundis104a , S. De Castro20a,20b , S. De Cecco80 ,N. De Groot106 , P. de Jong107 , H. De la Torre82 , F. De Lorenzi64 , L. De Nooij107 , D. De Pedis133a ,A. De Salvo133a , U. De Sanctis150 , A. De Santo150 , J.B. De Vivie De Regie117 , W.J. Dearnaley72 ,R. Debbe25 , C. Debenedetti138 , B. Dechenaux55 , D.V. Dedovich65 , I. Deigaard107 , J. Del Peso82 ,T. Del Prete124a,124b , F. Deliot137 , C.M. Delitzsch49 , M. Deliyergiyev75 , A. Dell’Acqua30 ,L. Dell’Asta22 , M. Dell’Orso124a,124b , M. Della Pietra104a,i , D. della Volpe49 , M. Delmastro5 ,P.A. Delsart55 , C. Deluca107 , S. Demers177 , M. Demichev65 , A. Demilly80 , S.P. Denisov130 ,D. Derendarz39 , J.E. Derkaoui136d , F. Derue80 , P. Dervan74 , K. Desch21 , C. Deterre42 ,P.O. Deviveiros107 , A. Dewhurst131 , S. Dhaliwal107 , A. Di Ciaccio134a,134b , L. Di Ciaccio5 ,A. Di Domenico133a,133b , C. Di Donato104a,104b , A. Di Girolamo30 , B. Di Girolamo30 ,A. Di Mattia153 , B. Di Micco135a,135b , R. Di Nardo47 , A. Di Simone48 , R. Di Sipio20a,20b ,D. Di Valentino29 , F.A. Dias46 , M.A. Diaz32a , E.B. Diehl89 , J. Dietrich42 , T.A. Dietzsch58a ,S. Diglio85 , A. Dimitrievska13a , J. Dingfelder21 , P. Dita26a , S. Dita26a , F. Dittus30 , F. Djama85 ,T. Djobava51b , J.I. Djuvsland58a , M.A.B. do Vale24c , A. Do Valle Wemans126a,126g , D. Dobos30 ,C. Doglioni49 , T. Doherty53 , T. Dohmae156 , J. Dolejsi129 , Z. Dolezal129 , B.A. Dolgoshein98,∗ ,M. Donadelli24d , S. Donati124a,124b , P. Dondero121a,121b , J. Donini34 , J. Dopke131 , A. Doria104a ,M.T. Dova71 , A.T. Doyle53 , M. Dris10 , J. Dubbert89 , S. Dube15 , E. Dubreuil34 , E. Duchovni173 ,G. Duckeck100 , O.A. Ducu26a , D. Duda176 , A. Dudarev30 , F. Dudziak64 , L. Duflot117 ,L. Duguid77 , M. Dührssen30 , M. Dunford58a , H. Duran Yildiz4a , M. Düren52 , A. Durglishvili51b ,M. Dwuznik38a , M. Dyndal38a , J. Ebke100 , W. Edson2 , N.C. Edwards46 , W. Ehrenfeld21 ,T. Eifert144 , G. Eigen14 , K. Einsweiler15 , T. Ekelof167 , M. El Kacimi136c , M. Ellert167 , S. Elles5 ,F. Ellinghaus83 , N. Ellis30 , J. Elmsheuser100 , M. Elsing30 , D. Emeliyanov131 , Y. Enari156 ,O.C. Endner83 , M. Endo118 , R. Engelmann149 , J. Erdmann177 , A. Ereditato17 , D. Eriksson147a ,G. Ernis176 , J. Ernst2 , M. Ernst25 , J. Ernwein137 , D. Errede166 , S. Errede166 , E. Ertel83 ,M. Escalier117 , H. Esch43 , C. Escobar125 , B. Esposito47 , A.I. Etienvre137 , E. Etzion154 ,H. Evans61 , A. Ezhilov123 , L. Fabbri20a,20b , G. Facini31 , R.M. Fakhrutdinov130 , S. Falciano133a ,R.J. Falla78 , J. Faltova129 , Y. Fang33a , M. Fanti91a,91b , A. Farbin8 , A. Farilla135a , T. Farooque12 ,S. Farrell15 , S.M. Farrington171 , P. Farthouat30 , F. Fassi136e , P. Fassnacht30 , D. Fassouliotis9 ,A. Favareto50a,50b , L. Fayard117 , P. Federic145a , O.L. Fedin123,k , W. Fedorko169 ,M. Fehling-Kaschek48 , S. Feigl30 , L. Feligioni85 , C. Feng33d , E.J. Feng6 , H. Feng89 ,A.B. Fenyuk130 , S. Fernandez Perez30 , S. Ferrag53 , J. Ferrando53 , A. Ferrari167 , P. Ferrari107 ,R. Ferrari121a , D.E. Ferreira de Lima53 , A. Ferrer168 , D. Ferrere49 , C. Ferretti89 ,A. Ferretto Parodi50a,50b , M. Fiascaris31 , F. Fiedler83 , A. Filipčič75 , M. Filipuzzi42 , F. Filthaut106 ,M. Fincke-Keeler170 , K.D. Finelli151 , M.C.N. Fiolhais126a,126c , L. Fiorini168 , A. Firan40 ,A. Fischer2 , J. Fischer176 , W.C. Fisher90 , E.A. Fitzgerald23 , M. Flechl48 , I. Fleck142 ,P. Fleischmann89 , S. Fleischmann176 , G.T. Fletcher140 , G. Fletcher76 , T. Flick176 , A. Floderus81 ,L.R. Flores Castillo60a , A.C. Florez Bustos160b , M.J. Flowerdew101 , A. Formica137 , A. Forti84 ,D. Fortin160a , D. Fournier117 , H. Fox72 , S. Fracchia12 , P. Francavilla80 , M. Franchini20a,20b ,S. Franchino30 , D. Francis30 , L. Franconi119 , M. Franklin57 , S. Franz62 , M. Fraternali121a,121b ,– 74 –JHEP01(2015)069

S.T. French28 , C. Friedrich42 , F. Friedrich44 , D. Froidevaux30 , J.A. Frost28 , C. Fukunaga157 ,E. Fullana Torregrosa83 , B.G. Fulsom144 , J. Fuster168 , C. Gabaldon55 , O. Gabizon176 ,A. Gabrielli20a,20b , A. Gabrielli133a,133b , S. Gadatsch107 , S. Gadomski49 , G. Gagliardi50a,50b ,P. Gagnon61 , C. Galea106 , B. Galhardo126a,126c , E.J. Gallas120 , V. Gallo17 , B.J. Gallop131 ,P. Gallus128 , G. Galster36 , K.K. Gan111 , J. Gao33b,h , Y.S. Gao144,f , F.M. Garay Walls46 ,F. Garberson177 , C. Garcı́a168 , J.E. Garcı́a Navarro168 , M. Garcia-Sciveres15 , R.W. Gardner31 ,N. Garelli144 , V. Garonne30 , C. Gatti47 , G. Gaudio121a , B. Gaur142 , L. Gauthier95 ,P. Gauzzi133a,133b , I.L. Gavrilenko96 , C. Gay169 , G. Gaycken21 , E.N. Gazis10 , P. Ge33d ,Z. Gecse169 , C.N.P. Gee131 , D.A.A. Geerts107 , Ch. Geich-Gimbel21 , K. Gellerstedt147a,147b ,C. Gemme50a , A. Gemmell53 , M.H. Genest55 , S. Gentile133a,133b , M. George54 , S. George77 ,D. Gerbaudo164 , A. Gershon154 , H. Ghazlane136b , N. Ghodbane34 , B. Giacobbe20a ,S. Giagu133a,133b , V. Giangiobbe12 , P. Giannetti124a,124b , F. Gianotti30 , B. Gibbard25 ,S.M. Gibson77 , M. Gilchriese15 , T.P.S. Gillam28 , D. Gillberg30 , G. Gilles34 , D.M. Gingrich3,e ,N. Giokaris9 , M.P. Giordani165a,165c , R. Giordano104a,104b , F.M. Giorgi20a , F.M. Giorgi16 ,P.F. Giraud137 , D. Giugni91a , C. Giuliani48 , M. Giulini58b , B.K. Gjelsten119 , S. Gkaitatzis155 ,I. Gkialas155,l , E.L. Gkougkousis117 , L.K. Gladilin99 , C. Glasman82 , J. Glatzer30 ,P.C.F. Glaysher46 , A. Glazov42 , G.L. Glonti65 , M. Goblirsch-Kolb101 , J.R. Goddard76 ,J. Godlewski30 , C. Goeringer83 , S. Goldfarb89 , T. Golling177 , D. Golubkov130 ,A. Gomes126a,126b,126d , L.S. Gomez Fajardo42 , R. Gonçalo126a ,J. Goncalves Pinto Firmino Da Costa137 , L. Gonella21 , S. González de la Hoz168 ,G. Gonzalez Parra12 , S. Gonzalez-Sevilla49 , L. Goossens30 , P.A. Gorbounov97 , H.A. Gordon25 ,I. Gorelov105 , B. Gorini30 , E. Gorini73a,73b , A. Gorišek75 , E. Gornicki39 , A.T. Goshaw6 ,C. Gössling43 , M.I. Gostkin65 , M. Gouighri136a , D. Goujdami136c , M.P. Goulette49 ,A.G. Goussiou139 , C. Goy5 , S. Gozpinar23 , H.M.X. Grabas138 , L. Graber54 , I. Grabowska-Bold38a ,P. Grafström20a,20b , K-J. Grahn42 , J. Gramling49 , E. Gramstad119 , S. Grancagnolo16 ,V. Grassi149 , V. Gratchev123 , H.M. Gray30 , E. Graziani135a , O.G. Grebenyuk123 ,Z.D. Greenwood79,m , K. Gregersen78 , I.M. Gregor42 , P. Grenier144 , J. Griffiths8 , A.A. Grillo138 ,K. Grimm72 , S. Grinstein12,n , Ph. Gris34 , Y.V. Grishkevich99 , J.-F. Grivaz117 , J.P. Grohs44 ,A. Grohsjean42 , E. Gross173 , J. Grosse-Knetter54 , G.C. Grossi134a,134b , J. Groth-Jensen173 ,Z.J. Grout150 , L. Guan33b , J. Guenther128 , F. Guescini49 , D. Guest177 , O. Gueta154 ,C. Guicheney34 , E. Guido50a,50b , T. Guillemin117 , S. Guindon2 , U. Gul53 , C. Gumpert44 ,J. Guo35 , S. Gupta120 , P. Gutierrez113 , N.G. Gutierrez Ortiz53 , C. Gutschow78 , N. Guttman154 ,C. Guyot137 , C. Gwenlan120 , C.B. Gwilliam74 , A. Haas110 , C. Haber15 , H.K. Hadavand8 ,N. Haddad136e , P. Haefner21 , S. Hageböck21 , Z. Hajduk39 , H. Hakobyan178 , M. Haleem42 ,D. Hall120 , G. Halladjian90 , K. Hamacher176 , P. Hamal115 , K. Hamano170 , M. Hamer54 ,A. Hamilton146a , S. Hamilton162 , G.N. Hamity146c , P.G. Hamnett42 , L. Han33b , K. Hanagaki118 ,K. Hanawa156 , M. Hance15 , P. Hanke58a , R. Hanna137 , J.B. Hansen36 , J.D. Hansen36 ,P.H. Hansen36 , K. Hara161 , A.S. Hard174 , T. Harenberg176 , F. Hariri117 , S. Harkusha92 ,D. Harper89 , R.D. Harrington46 , O.M. Harris139 , P.F. Harrison171 , F. Hartjes107 , M. Hasegawa67 ,S. Hasegawa103 , Y. Hasegawa141 , A. Hasib113 , S. Hassani137 , S. Haug17 , M. Hauschild30 ,R. Hauser90 , M. Havranek127 , C.M. Hawkes18 , R.J. Hawkings30 , A.D. Hawkins81 , T. Hayashi161 ,D. Hayden90 , C.P. Hays120 , H.S. Hayward74 , S.J. Haywood131 , S.J. Head18 , T. Heck83 ,V. Hedberg81 , L. Heelan8 , S. Heim122 , T. Heim176 , B. Heinemann15 , L. Heinrich110 , J. Hejbal127 ,L. Helary22 , C. Heller100 , M. Heller30 , S. Hellman147a,147b , D. Hellmich21 , C. Helsens30 ,J. Henderson120 , R.C.W. Henderson72 , Y. Heng174 , C. Hengler42 , A. Henrichs177 ,A.M. Henriques Correia30 , S. Henrot-Versille117 , G.H. Herbert16 , Y. Hernández Jiménez168 ,R. Herrberg-Schubert16 , G. Herten48 , R. Hertenberger100 , L. Hervas30 , G.G. Hesketh78 ,N.P. Hessey107 , R. Hickling76 , E. Higón-Rodriguez168 , E. Hill170 , J.C. Hill28 , K.H. Hiller42 ,– 75 –JHEP01(2015)069

S. Hillert21 , S.J. Hillier18 , I. Hinchliffe15 , E. Hines122 , M. Hirose158 , D. Hirschbuehl176 ,J. Hobbs149 , N. Hod107 , M.C. Hodgkinson140 , P. Hodgson140 , A. Hoecker30 , M.R. Hoeferkamp105 ,F. Hoenig100 , J. Hoffman40 , D. Hoffmann85 , M. Hohlfeld83 , T.R. Holmes15 , T.M. Hong122 ,L. Hooft van Huysduynen110 , W.H. Hopkins116 , Y. Horii103 , J-Y. Hostachy55 , S. Hou152 ,A. Hoummada136a , J. Howard120 , J. Howarth42 , M. Hrabovsky115 , I. Hristova16 , J. Hrivnac117 ,T. Hryn’ova5 , C. Hsu146c , P.J. Hsu83 , S.-C. Hsu139 , D. Hu35 , X. Hu89 , Y. Huang42 , Z. Hubacek30 ,F. Hubaut85 , F. Huegging21 , T.B. Huffman120 , E.W. Hughes35 , G. Hughes72 , M. Huhtinen30 ,T.A. Hülsing83 , M. Hurwitz15 , N. Huseynov65,b , J. Huston90 , J. Huth57 , G. Iacobucci49 ,G. Iakovidis10 , I. Ibragimov142 , L. Iconomidou-Fayard117 , E. Ideal177 , Z. Idrissi136e , P. Iengo104a ,O. Igonkina107 , T. Iizawa172 , Y. Ikegami66 , K. Ikematsu142 , M. Ikeno66 , Y. Ilchenko31,o ,D. Iliadis155 , N. Ilic159 , Y. Inamaru67 , T. Ince101 , P. Ioannou9 , M. Iodice135a , K. Iordanidou9 ,V. Ippolito57 , A. Irles Quiles168 , C. Isaksson167 , M. Ishino68 , M. Ishitsuka158 ,R. Ishmukhametov111 , C. Issever120 , S. Istin19a , J.M. Iturbe Ponce84 , R. Iuppa134a,134b ,J. Ivarsson81 , W. Iwanski39 , H. Iwasaki66 , J.M. Izen41 , V. Izzo104a , B. Jackson122 , M. Jackson74 ,P. Jackson1 , M.R. Jaekel30 , V. Jain2 , K. Jakobs48 , S. Jakobsen30 , T. Jakoubek127 , J. Jakubek128 ,D.O. Jamin152 , D.K. Jana79 , E. Jansen78 , H. Jansen30 , J. Janssen21 , M. Janus171 , G. Jarlskog81 ,N. Javadov65,b , T. Javůrek48 , L. Jeanty15 , J. Jejelava51a,p , G.-Y. Jeng151 , D. Jennens88 ,P. Jenni48,q , J. Jentzsch43 , C. Jeske171 , S. Jézéquel5 , H. Ji174 , J. Jia149 , Y. Jiang33b ,M. Jimenez Belenguer42 , S. Jin33a , A. Jinaru26a , O. Jinnouchi158 , M.D. Joergensen36 ,K.E. Johansson147a,147b , P. Johansson140 , K.A. Johns7 , K. Jon-And147a,147b , G. Jones171 ,R.W.L. Jones72 , T.J. Jones74 , J. Jongmanns58a , P.M. Jorge126a,126b , K.D. Joshi84 , J. Jovicevic148 ,X. Ju174 , C.A. Jung43 , R.M. Jungst30 , P. Jussel62 , A. Juste Rozas12,n , M. Kaci168 ,A. Kaczmarska39 , M. Kado117 , H. Kagan111 , M. Kagan144 , E. Kajomovitz45 , C.W. Kalderon120 ,S. Kama40 , A. Kamenshchikov130 , N. Kanaya156 , M. Kaneda30 , S. Kaneti28 , V.A. Kantserov98 ,J. Kanzaki66 , B. Kaplan110 , A. Kapliy31 , D. Kar53 , K. Karakostas10 , N. Karastathis10 ,M.J. Kareem54 , M. Karnevskiy83 , S.N. Karpov65 , Z.M. Karpova65 , K. Karthik110 ,V. Kartvelishvili72 , A.N. Karyukhin130 , L. Kashif174 , G. Kasieczka58b , R.D. Kass111 ,A. Kastanas14 , Y. Kataoka156 , A. Katre49 , J. Katzy42 , V. Kaushik7 , K. Kawagoe70 ,T. Kawamoto156 , G. Kawamura54 , S. Kazama156 , V.F. Kazanin109 , M.Y. Kazarinov65 ,R. Keeler170 , R. Kehoe40 , M. Keil54 , J.S. Keller42 , J.J. Kempster77 , H. Keoshkerian5 ,O. Kepka127 , B.P. Kerševan75 , S. Kersten176 , K. Kessoku156 , J. Keung159 , F. Khalil-zada11 ,H. Khandanyan147a,147b , A. Khanov114 , A. Khodinov98 , A. Khomich58a , T.J. Khoo28 ,G. Khoriauli21 , A. Khoroshilov176 , V. Khovanskiy97 , E. Khramov65 , J. Khubua51b , H.Y. Kim8 ,H. Kim147a,147b , S.H. Kim161 , N. Kimura172 , O. Kind16 , B.T. King74 , M. King168 , R.S.B. King120 ,S.B. King169 , J. Kirk131 , A.E. Kiryunin101 , T. Kishimoto67 , D. Kisielewska38a , F. Kiss48 ,K. Kiuchi161 , E. Kladiva145b , M. Klein74 , U. Klein74 , K. Kleinknecht83 , P. Klimek147a,147b ,A. Klimentov25 , R. Klingenberg43 , J.A. Klinger84 , T. Klioutchnikova30 , P.F. Klok106 ,E.-E. Kluge58a , P. Kluit107 , S. Kluth101 , E. Kneringer62 , E.B.F.G. Knoops85 , A. Knue53 ,D. Kobayashi158 , T. Kobayashi156 , M. Kobel44 , M. Kocian144 , P. Kodys129 , T. Koffas29 ,E. Koffeman107 , L.A. Kogan120 , S. Kohlmann176 , Z. Kohout128 , T. Kohriki66 , T. Koi144 ,H. Kolanoski16 , I. Koletsou5 , J. Koll90 , A.A. Komar96,∗ , Y. Komori156 , T. Kondo66 ,N. Kondrashova42 , K. Köneke48 , A.C. König106 , S. König83 , T. Kono66,r , R. Konoplich110,s ,N. Konstantinidis78 , R. Kopeliansky153 , S. Koperny38a , L. Köpke83 , A.K. Kopp48 , K. Korcyl39 ,K. Kordas155 , A. Korn78 , A.A. Korol109,c , I. Korolkov12 , E.V. Korolkova140 , V.A. Korotkov130 ,O. Kortner101 , S. Kortner101 , V.V. Kostyukhin21 , V.M. Kotov65 , A. Kotwal45 , C. Kourkoumelis9 ,V. Kouskoura155 , A. Koutsman160a , R. Kowalewski170 , T.Z. Kowalski38a , W. Kozanecki137 ,A.S. Kozhin130 , V.A. Kramarenko99 , G. Kramberger75 , D. Krasnopevtsev98 , M.W. Krasny80 ,A. Krasznahorkay30 , J.K. Kraus21 , A. Kravchenko25 , S. Kreiss110 , M. Kretz58c , J. Kretzschmar74 ,– 76 –JHEP01(2015)069

K. Kreutzfeldt52 , P. Krieger159 , K. Kroeninger54 , H. Kroha101 , J. Kroll122 , J. Kroseberg21 ,J. Krstic13a , U. Kruchonak65 , H. Krüger21 , T. Kruker17 , N. Krumnack64 , Z.V. Krumshteyn65 ,A. Kruse174 , M.C. Kruse45 , M. Kruskal22 , T. Kubota88 , H. Kucuk78 , S. Kuday4c , S. Kuehn48 ,A. Kugel58c , A. Kuhl138 , T. Kuhl42 , V. Kukhtin65 , Y. Kulchitsky92 , S. Kuleshov32b ,M. Kuna133a,133b , J. Kunkle122 , A. Kupco127 , H. Kurashige67 , Y.A. Kurochkin92 , R. Kurumida67 ,V. Kus127 , E.S. Kuwertz148 , M. Kuze158 , J. Kvita115 , A. La Rosa49 , L. La Rotonda37a,37b ,C. Lacasta168 , F. Lacava133a,133b , J. Lacey29 , H. Lacker16 , D. Lacour80 , V.R. Lacuesta168 ,E. Ladygin65 , R. Lafaye5 , B. Laforge80 , T. Lagouri177 , S. Lai48 , H. Laier58a , L. Lambourne78 ,S. Lammers61 , C.L. Lampen7 , W. Lampl7 , E. Lançon137 , U. Landgraf48 , M.P.J. Landon76 ,V.S. Lang58a , A.J. Lankford164 , F. Lanni25 , K. Lantzsch30 , S. Laplace80 , C. Lapoire21 ,J.F. Laporte137 , T. Lari91a , F. Lasagni Manghi20a,20b , M. Lassnig30 , P. Laurelli47 , W. Lavrijsen15 ,A.T. Law138 , P. Laycock74 , O. Le Dortz80 , E. Le Guirriec85 , E. Le Menedeu12 , T. LeCompte6 ,F. Ledroit-Guillon55 , C.A. Lee146b , H. Lee107 , J.S.H. Lee118 , S.C. Lee152 , L. Lee1 , G. Lefebvre80 ,M. Lefebvre170 , F. Legger100 , C. Leggett15 , A. Lehan74 , G. Lehmann Miotto30 , X. Lei7 ,W.A. Leight29 , A. Leisos155 , A.G. Leister177 , M.A.L. Leite24d , R. Leitner129 , D. Lellouch173 ,B. Lemmer54 , K.J.C. Leney78 , T. Lenz21 , G. Lenzen176 , B. Lenzi30 , R. Leone7 , S. Leone124a,124b ,C. Leonidopoulos46 , S. Leontsinis10 , C. Leroy95 , C.G. Lester28 , C.M. Lester122 , M. Levchenko123 ,J. Levêque5 , D. Levin89 , L.J. Levinson173 , M. Levy18 , A. Lewis120 , G.H. Lewis110 , A.M. Leyko21 ,M. Leyton41 , B. Li33b,t , B. Li85 , H. Li149 , H.L. Li31 , L. Li45 , L. Li33e , S. Li45 , Y. Li33c,u ,Z. Liang138 , H. Liao34 , B. Liberti134a , P. Lichard30 , K. Lie166 , J. Liebal21 , W. Liebig14 ,C. Limbach21 , A. Limosani151 , S.C. Lin152,v , T.H. Lin83 , F. Linde107 , B.E. Lindquist149 ,J.T. Linnemann90 , E. Lipeles122 , A. Lipniacka14 , M. Lisovyi42 , T.M. Liss166 , D. Lissauer25 ,A. Lister169 , A.M. Litke138 , B. Liu152 , D. Liu152 , J.B. Liu33b , K. Liu33b,w , L. Liu89 , M. Liu45 ,M. Liu33b , Y. Liu33b , M. Livan121a,121b , A. Lleres55 , J. Llorente Merino82 , S.L. Lloyd76 ,F. Lo Sterzo152 , E. Lobodzinska42 , P. Loch7 , W.S. Lockman138 , T. Loddenkoetter21 ,F.K. Loebinger84 , A.E. Loevschall-Jensen36 , A. Loginov177 , T. Lohse16 , K. Lohwasser42 ,M. Lokajicek127 , V.P. Lombardo5 , B.A. Long22 , J.D. Long89 , R.E. Long72 , L. Lopes126a ,D. Lopez Mateos57 , B. Lopez Paredes140 , I. Lopez Paz12 , J. Lorenz100 , N. Lorenzo Martinez61 ,M. Losada163 , P. Loscutoff15 , X. Lou41 , A. Lounis117 , J. Love6 , P.A. Love72 , A.J. Lowe144,f ,F. Lu33a , N. Lu89 , H.J. Lubatti139 , C. Luci133a,133b , A. Lucotte55 , F. Luehring61 , W. Lukas62 ,L. Luminari133a , O. Lundberg147a,147b , B. Lund-Jensen148 , M. Lungwitz83 , D. Lynn25 ,R. Lysak127 , E. Lytken81 , H. Ma25 , L.L. Ma33d , G. Maccarrone47 , A. Macchiolo101 ,J. Machado Miguens126a,126b , D. Macina30 , D. Madaffari85 , R. Madar48 , H.J. Maddocks72 ,W.F. Mader44 , A. Madsen167 , M. Maeno8 , T. Maeno25 , A. Maevskiy99 , E. Magradze54 ,K. Mahboubi48 , J. Mahlstedt107 , S. Mahmoud74 , C. Maiani137 , C. Maidantchik24a , A.A. Maier101 ,A. Maio126a,126b,126d , S. Majewski116 , Y. Makida66 , N. Makovec117 , P. Mal137,x , B. Malaescu80 ,Pa. Malecki39 , V.P. Maleev123 , F. Malek55 , U. Mallik63 , D. Malon6 , C. Malone144 , S. Maltezos10 ,V.M. Malyshev109 , S. Malyukov30 , J. Mamuzic13b , B. Mandelli30 , L. Mandelli91a , I. Mandić75 ,R. Mandrysch63 , J. Maneira126a,126b , A. Manfredini101 , L. Manhaes de Andrade Filho24b ,J.A. Manjarres Ramos160b , A. Mann100 , P.M. Manning138 , A. Manousakis-Katsikakis9 ,B. Mansoulie137 , R. Mantifel87 , L. Mapelli30 , L. March146c , J.F. Marchand29 , G. Marchiori80 ,M. Marcisovsky127 , C.P. Marino170 , M. Marjanovic13a , C.N. Marques126a , F. Marroquim24a ,S.P. Marsden84 , Z. Marshall15 , L.F. Marti17 , S. Marti-Garcia168 , B. Martin30 , B. Martin90 ,T.A. Martin171 , V.J. Martin46 , B. Martin dit Latour14 , H. Martinez137 , M. Martinez12,n ,S. Martin-Haugh131 , A.C. Martyniuk78 , M. Marx139 , F. Marzano133a , A. Marzin30 , L. Masetti83 ,T. Mashimo156 , R. Mashinistov96 , J. Masik84 , A.L. Maslennikov109,c , I. Massa20a,20b ,L. Massa20a,20b , N. Massol5 , P. Mastrandrea149 , A. Mastroberardino37a,37b , T. Masubuchi156 ,P. Mättig176 , J. Mattmann83 , J. Maurer26a , S.J. Maxfield74 , D.A. Maximov109,c , R. Mazini152 ,– 77 –JHEP01(2015)069

L. Mazzaferro134a,134b , G. Mc Goldrick159 , S.P. Mc Kee89 , A. McCarn89 , R.L. McCarthy149 ,T.G. McCarthy29 , N.A. McCubbin131 , K.W. McFarlane56,∗ , J.A. Mcfayden78 , G. Mchedlidze54 ,S.J. McMahon131 , R.A. McPherson170,j , J. Mechnich107 , M. Medinnis42 , S. Meehan31 ,S. Mehlhase100 , A. Mehta74 , K. Meier58a , C. Meineck100 , B. Meirose81 , C. Melachrinos31 ,B.R. Mellado Garcia146c , F. Meloni17 , A. Mengarelli20a,20b , S. Menke101 , E. Meoni162 ,K.M. Mercurio57 , S. Mergelmeyer21 , N. Meric137 , P. Mermod49 , L. Merola104a,104b , C. Meroni91a ,F.S. Merritt31 , H. Merritt111 , A. Messina30,y , J. Metcalfe25 , A.S. Mete164 , C. Meyer83 ,C. Meyer122 , J-P. Meyer137 , J. Meyer30 , R.P. Middleton131 , S. Migas74 , L. Mijović21 ,G. Mikenberg173 , M. Mikestikova127 , M. Mikuž75 , A. Milic30 , D.W. Miller31 , C. Mills46 ,A. Milov173 , D.A. Milstead147a,147b , D. Milstein173 , A.A. Minaenko130 , Y. Minami156 ,I.A. Minashvili65 , A.I. Mincer110 , B. Mindur38a , M. Mineev65 , Y. Ming174 , L.M. Mir12 ,G. Mirabelli133a , T. Mitani172 , J. Mitrevski100 , V.A. Mitsou168 , A. Miucci49 , P.S. Miyagawa140 ,J.U. Mjörnmark81 , T. Moa147a,147b , K. Mochizuki85 , S. Mohapatra35 , W. Mohr48 ,S. Molander147a,147b , R. Moles-Valls168 , K. Mönig42 , C. Monini55 , J. Monk36 , E. Monnier85 ,J. Montejo Berlingen12 , F. Monticelli71 , S. Monzani133a,133b , R.W. Moore3 , N. Morange63 ,D. Moreno83 , M. Moreno Llácer54 , P. Morettini50a , M. Morgenstern44 , M. Morii57 ,V. Morisbak119 , S. Moritz83 , A.K. Morley148 , G. Mornacchi30 , J.D. Morris76 , L. Morvaj103 ,H.G. Moser101 , M. Mosidze51b , J. Moss111 , K. Motohashi158 , R. Mount144 , E. Mountricha25 ,S.V. Mouraviev96,∗ , E.J.W. Moyse86 , S. Muanza85 , R.D. Mudd18 , F. Mueller58a , J. Mueller125 ,K. Mueller21 , T. Mueller28 , T. Mueller83 , D. Muenstermann49 , Y. Munwes154 ,J.A. Murillo Quijada18 , W.J. Murray171,131 , H. Musheghyan54 , E. Musto153 , A.G. Myagkov130,z ,M. Myska128 , O. Nackenhorst54 , J. Nadal54 , K. Nagai120 , R. Nagai158 , Y. Nagai85 , K. Nagano66 ,A. Nagarkar111 , Y. Nagasaka59 , K. Nagata161 , M. Nagel101 , A.M. Nairz30 , Y. Nakahama30 ,K. Nakamura66 , T. Nakamura156 , I. Nakano112 , H. Namasivayam41 , G. Nanava21 ,R.F. Naranjo Garcia42 , R. Narayan58b , T. Nattermann21 , T. Naumann42 , G. Navarro163 ,R. Nayyar7 , H.A. Neal89 , P.Yu. Nechaeva96 , T.J. Neep84 , P.D. Nef144 , A. Negri121a,121b ,G. Negri30 , M. Negrini20a , S. Nektarijevic49 , C. Nellist117 , A. Nelson164 , T.K. Nelson144 ,S. Nemecek127 , P. Nemethy110 , A.A. Nepomuceno24a , M. Nessi30,aa , M.S. Neubauer166 ,M. Neumann176 , R.M. Neves110 , P. Nevski25 , P.R. Newman18 , D.H. Nguyen6 , R.B. Nickerson120 ,R. Nicolaidou137 , B. Nicquevert30 , J. Nielsen138 , N. Nikiforou35 , A. Nikiforov16 ,V. Nikolaenko130,z , I. Nikolic-Audit80 , K. Nikolics49 , K. Nikolopoulos18 , P. Nilsson25 ,Y. Ninomiya156 , A. Nisati133a , R. Nisius101 , T. Nobe158 , L. Nodulman6 , M. Nomachi118 ,I. Nomidis29 , S. Norberg113 , M. Nordberg30 , O. Novgorodova44 , S. Nowak101 , M. Nozaki66 ,L. Nozka115 , K. Ntekas10 , G. Nunes Hanninger88 , T. Nunnemann100 , E. Nurse78 , F. Nuti88 ,B.J. O’Brien46 , F. O’grady7 , D.C. O’Neil143 , V. O’Shea53 , F.G. Oakham29,e , H. Oberlack101 ,T. Obermann21 , J. Ocariz80 , A. Ochi67 , M.I. Ochoa78 , S. Oda70 , S. Odaka66 , H. Ogren61 , A. Oh84 ,S.H. Oh45 , C.C. Ohm15 , H. Ohman167 , H. Oide30 , W. Okamura118 , H. Okawa25 , Y. Okumura31 ,T. Okuyama156 , A. Olariu26a , A.G. Olchevski65 , S.A. Olivares Pino46 , D. Oliveira Damazio25 ,E. Oliver Garcia168 , A. Olszewski39 , J. Olszowska39 , A. Onofre126a,126e , P.U.E. Onyisi31,o ,C.J. Oram160a , M.J. Oreglia31 , Y. Oren154 , D. Orestano135a,135b , N. Orlando73a,73b ,C. Oropeza Barrera53 , R.S. Orr159 , B. Osculati50a,50b , R. Ospanov122 , G. Otero y Garzon27 ,H. Otono70 , M. Ouchrif136d , E.A. Ouellette170 , F. Ould-Saada119 , A. Ouraou137 ,K.P. Oussoren107 , Q. Ouyang33a , A. Ovcharova15 , M. Owen84 , V.E. Ozcan19a , N. Ozturk8 ,K. Pachal120 , A. Pacheco Pages12 , C. Padilla Aranda12 , M. Pagáčová48 , S. Pagan Griso15 ,E. Paganis140 , C. Pahl101 , F. Paige25 , P. Pais86 , K. Pajchel119 , G. Palacino160b , S. Palestini30 ,M. Palka38b , D. Pallin34 , A. Palma126a,126b , J.D. Palmer18 , Y.B. Pan174 , E. Panagiotopoulou10 ,C.E. Pandini80 , J.G. Panduro Vazquez77 , P. Pani107 , N. Panikashvili89 , S. Panitkin25 ,D. Pantea26a , L. Paolozzi134a,134b , Th.D. Papadopoulou10 , K. Papageorgiou155,l , A. Paramonov6 ,– 78 –JHEP01(2015)069

D. Paredes Hernandez155 , M.A. Parker28 , F. Parodi50a,50b , J.A. Parsons35 , U. Parzefall48 ,E. Pasqualucci133a , S. Passaggio50a , A. Passeri135a , F. Pastore135a,135b,∗ , Fr. Pastore77 ,G. Pásztor29 , S. Pataraia176 , N.D. Patel151 , J.R. Pater84 , S. Patricelli104a,104b , T. Pauly30 ,J. Pearce170 , L.E. Pedersen36 , M. Pedersen119 , S. Pedraza Lopez168 , R. Pedro126a,126b ,S.V. Peleganchuk109 , D. Pelikan167 , H. Peng33b , B. Penning31 , J. Penwell61 , D.V. Perepelitsa25 ,E. Perez Codina160a , M.T. Pérez Garcı́a-Estañ168 , L. Perini91a,91b , H. Pernegger30 ,S. Perrella104a,104b , R. Perrino73a , R. Peschke42 , V.D. Peshekhonov65 , K. Peters30 ,R.F.Y. Peters84 , B.A. Petersen30 , T.C. Petersen36 , E. Petit42 , A. Petridis147a,147b , C. Petridou155 ,E. Petrolo133a , F. Petrucci135a,135b , N.E. Pettersson158 , R. Pezoa32b , P.W. Phillips131 ,G. Piacquadio144 , E. Pianori171 , A. Picazio49 , E. Piccaro76 , M. Piccinini20a,20b , R. Piegaia27 ,D.T. Pignotti111 , J.E. Pilcher31 , A.D. Pilkington78 , J. Pina126a,126b,126d , M. Pinamonti165a,165c,ab ,A. Pinder120 , J.L. Pinfold3 , A. Pingel36 , B. Pinto126a , S. Pires80 , M. Pitt173 , C. Pizio91a,91b ,L. Plazak145a , M.-A. Pleier25 , V. Pleskot129 , E. Plotnikova65 , P. Plucinski147a,147b , D. Pluth64 ,S. Poddar58a , F. Podlyski34 , R. Poettgen83 , L. Poggioli117 , D. Pohl21 , M. Pohl49 , G. Polesello121a ,A. Policicchio37a,37b , R. Polifka159 , A. Polini20a , C.S. Pollard45 , V. Polychronakos25 ,K. Pommès30 , L. Pontecorvo133a , B.G. Pope90 , G.A. Popeneciu26b , D.S. Popovic13a ,A. Poppleton30 , X. Portell Bueso12 , S. Pospisil128 , K. Potamianos15 , I.N. Potrap65 , C.J. Potter150 ,C.T. Potter116 , G. Poulard30 , J. Poveda61 , V. Pozdnyakov65 , P. Pralavorio85 , A. Pranko15 ,S. Prasad30 , R. Pravahan8 , S. Prell64 , D. Price84 , J. Price74 , L.E. Price6 , D. Prieur125 ,M. Primavera73a , M. Proissl46 , K. Prokofiev47 , F. Prokoshin32b , E. Protopapadaki137 ,S. Protopopescu25 , J. Proudfoot6 , M. Przybycien38a , H. Przysiezniak5 , E. Ptacek116 ,D. Puddu135a,135b , E. Pueschel86 , D. Puldon149 , M. Purohit25,ac , P. Puzo117 , J. Qian89 , G. Qin53 ,Y. Qin84 , A. Quadt54 , D.R. Quarrie15 , W.B. Quayle165a,165b , M. Queitsch-Maitland84 ,D. Quilty53 , A. Qureshi160b , V. Radeka25 , V. Radescu42 , S.K. Radhakrishnan149 , P. Radloff116 ,P. Rados88 , F. Ragusa91a,91b , G. Rahal179 , S. Rajagopalan25 , M. Rammensee30 ,A.S. Randle-Conde40 , C. Rangel-Smith167 , K. Rao164 , F. Rauscher100 , T.C. Rave48 ,T. Ravenscroft53 , M. Raymond30 , A.L. Read119 , N.P. Readioff74 , D.M. Rebuzzi121a,121b ,A. Redelbach175 , G. Redlinger25 , R. Reece138 , K. Reeves41 , L. Rehnisch16 , H. Reisin27 ,M. Relich164 , C. Rembser30 , H. Ren33a , Z.L. Ren152 , A. Renaud117 , M. Rescigno133a ,S. Resconi91a , O.L. Rezanova109,c , P. Reznicek129 , R. Rezvani95 , R. Richter101 , M. Ridel80 ,P. Rieck16 , J. Rieger54 , M. Rijssenbeek149 , A. Rimoldi121a,121b , L. Rinaldi20a , E. Ritsch62 , I. Riu12 ,F. Rizatdinova114 , E. Rizvi76 , S.H. Robertson87,j , A. Robichaud-Veronneau87 , D. Robinson28 ,J.E.M. Robinson84 , A. Robson53 , C. Roda124a,124b , L. Rodrigues30 , S. Roe30 , O. Røhne119 ,S. Rolli162 , A. Romaniouk98 , M. Romano20a,20b , E. Romero Adam168 , N. Rompotis139 ,M. Ronzani48 , L. Roos80 , E. Ros168 , S. Rosati133a , K. Rosbach49 , M. Rose77 , P. Rose138 ,P.L. Rosendahl14 , O. Rosenthal142 , V. Rossetti147a,147b , E. Rossi104a,104b , L.P. Rossi50a ,R. Rosten139 , M. Rotaru26a , I. Roth173 , J. Rothberg139 , D. Rousseau117 , C.R. Royon137 ,A. Rozanov85 , Y. Rozen153 , X. Ruan146c , F. Rubbo12 , I. Rubinskiy42 , V.I. Rud99 , C. Rudolph44 ,M.S. Rudolph159 , F. Rühr48 , A. Ruiz-Martinez30 , Z. Rurikova48 , N.A. Rusakovich65 ,A. Ruschke100 , J.P. Rutherfoord7 , N. Ruthmann48 , Y.F. Ryabov123 , M. Rybar129 , G. Rybkin117 ,N.C. Ryder120 , A.F. Saavedra151 , G. Sabato107 , S. Sacerdoti27 , A. Saddique3 , I. Sadeh154 ,H.F-W. Sadrozinski138 , R. Sadykov65 , F. Safai Tehrani133a , H. Sakamoto156 , Y. Sakurai172 ,G. Salamanna135a,135b , A. Salamon134a , M. Saleem113 , D. Salek107 , P.H. Sales De Bruin139 ,D. Salihagic101 , A. Salnikov144 , J. Salt168 , D. Salvatore37a,37b , F. Salvatore150 , A. Salvucci106 ,A. Salzburger30 , D. Sampsonidis155 , A. Sanchez104a,104b , J. Sánchez168 , V. Sanchez Martinez168 ,H. Sandaker14 , R.L. Sandbach76 , H.G. Sander83 , M.P. Sanders100 , M. Sandhoff176 , T. Sandoval28 ,C. Sandoval163 , R. Sandstroem101 , D.P.C. Sankey131 , A. Sansoni47 , C. Santoni34 ,R. Santonico134a,134b , H. Santos126a , I. Santoyo Castillo150 , K. Sapp125 , A. Sapronov65 ,– 79 –JHEP01(2015)069

J.G. Saraiva126a,126d , B. Sarrazin21 , G. Sartisohn176 , O. Sasaki66 , Y. Sasaki156 , G. Sauvage5,∗ ,E. Sauvan5 , P. Savard159,e , D.O. Savu30 , C. Sawyer120 , L. Sawyer79,m , D.H. Saxon53 , J. Saxon122 ,C. Sbarra20a , A. Sbrizzi20a,20b , T. Scanlon78 , D.A. Scannicchio164 , M. Scarcella151 ,V. Scarfone37a,37b , J. Schaarschmidt173 , P. Schacht101 , D. Schaefer30 , R. Schaefer42 , S. Schaepe21 ,S. Schaetzel58b , U. Schäfer83 , A.C. Schaffer117 , D. Schaile100 , R.D. Schamberger149 , V. Scharf58a ,V.A. Schegelsky123 , D. Scheirich129 , M. Schernau164 , M.I. Scherzer35 , C. Schiavi50a,50b ,J. Schieck100 , C. Schillo48 , M. Schioppa37a,37b , S. Schlenker30 , E. Schmidt48 , K. Schmieden30 ,C. Schmitt83 , S. Schmitt58b , B. Schneider17 , Y.J. Schnellbach74 , U. Schnoor44 , L. Schoeffel137 ,A. Schoening58b , B.D. Schoenrock90 , A.L.S. Schorlemmer54 , M. Schott83 , D. Schouten160a ,J. Schovancova25 , S. Schramm159 , M. Schreyer175 , C. Schroeder83 , N. Schuh83 , M.J. Schultens21 ,H.-C. Schultz-Coulon58a , H. Schulz16 , M. Schumacher48 , B.A. Schumm138 , Ph. Schune137 ,C. Schwanenberger84 , A. Schwartzman144 , T.A. Schwarz89 , Ph. Schwegler101 , Ph. Schwemling137 ,R. Schwienhorst90 , J. Schwindling137 , T. Schwindt21 , M. Schwoerer5 , F.G. Sciacca17 , E. Scifo117 ,G. Sciolla23 , W.G. Scott131 , F. Scuri124a,124b , F. Scutti21 , J. Searcy89 , G. Sedov42 , E. Sedykh123 ,S.C. Seidel105 , A. Seiden138 , F. Seifert128 , J.M. Seixas24a , G. Sekhniaidze104a , S.J. Sekula40 ,K.E. Selbach46 , D.M. Seliverstov123,∗ , G. Sellers74 , N. Semprini-Cesari20a,20b , C. Serfon30 ,L. Serin117 , L. Serkin54 , T. Serre85 , R. Seuster160a , H. Severini113 , T. Sfiligoj75 , F. Sforza101 ,A. Sfyrla30 , E. Shabalina54 , M. Shamim116 , L.Y. Shan33a , R. Shang166 , J.T. Shank22 ,M. Shapiro15 , P.B. Shatalov97 , K. Shaw165a,165b , C.Y. Shehu150 , P. Sherwood78 , L. Shi152,ad ,S. Shimizu67 , C.O. Shimmin164 , M. Shimojima102 , M. Shiyakova65 , A. Shmeleva96 ,M.J. Shochet31 , D. Short120 , S. Shrestha64 , E. Shulga98 , M.A. Shupe7 , S. Shushkevich42 ,P. Sicho127 , O. Sidiropoulou155 , D. Sidorov114 , A. Sidoti133a , F. Siegert44 , Dj. Sijacki13a ,J. Silva126a,126d , Y. Silver154 , D. Silverstein144 , S.B. Silverstein147a , V. Simak128 , O. Simard5 ,Lj. Simic13a , S. Simion117 , E. Simioni83 , B. Simmons78 , R. Simoniello91a,91b , P. Sinervo159 ,N.B. Sinev116 , G. Siragusa175 , A. Sircar79 , A.N. Sisakyan65,∗ , S.Yu. Sivoklokov99 ,J. Sjölin147a,147b , T.B. Sjursen14 , H.P. Skottowe57 , K.Yu. Skovpen109 , P. Skubic113 , M. Slater18 ,T. Slavicek128 , M. Slawinska107 , K. Sliwa162 , V. Smakhtin173 , B.H. Smart46 , L. Smestad14 ,S.Yu. Smirnov98 , Y. Smirnov98 , L.N. Smirnova99,ae , O. Smirnova81 , K.M. Smith53 ,M. Smizanska72 , K. Smolek128 , A.A. Snesarev96 , G. Snidero76 , S. Snyder25 , R. Sobie170,j ,F. Socher44 , A. Soffer154 , D.A. Soh152,ad , C.A. Solans30 , M. Solar128 , J. Solc128 , E.Yu. Soldatov98 ,U. Soldevila168 , A.A. Solodkov130 , A. Soloshenko65 , O.V. Solovyanov130 , V. Solovyev123 ,P. Sommer48 , H.Y. Song33b , N. Soni1 , A. Sood15 , A. Sopczak128 , B. Sopko128 , V. Sopko128 ,V. Sorin12 , M. Sosebee8 , R. Soualah165a,165c , P. Soueid95 , A.M. Soukharev109,c , D. South42 ,S. Spagnolo73a,73b , F. Spanò77 , W.R. Spearman57 , F. Spettel101 , R. Spighi20a , G. Spigo30 ,L.A. Spiller88 , M. Spousta129 , T. Spreitzer159 , B. Spurlock8 , R.D. St. Denis53,∗ , S. Staerz44 ,J. Stahlman122 , R. Stamen58a , S. Stamm16 , E. Stanecka39 , R.W. Stanek6 , C. Stanescu135a ,M. Stanescu-Bellu42 , M.M. Stanitzki42 , S. Stapnes119 , E.A. Starchenko130 , J. Stark55 ,P. Staroba127 , P. Starovoitov42 , R. Staszewski39 , P. Stavina145a,∗ , P. Steinberg25 , B. Stelzer143 ,H.J. Stelzer30 , O. Stelzer-Chilton160a , H. Stenzel52 , S. Stern101 , G.A. Stewart53 , J.A. Stillings21 ,M.C. Stockton87 , M. Stoebe87 , G. Stoicea26a , P. Stolte54 , S. Stonjek101 , A.R. Stradling8 ,A. Straessner44 , M.E. Stramaglia17 , J. Strandberg148 , S. Strandberg147a,147b , A. Strandlie119 ,E. Strauss144 , M. Strauss113 , P. Strizenec145b , R. Ströhmer175 , D.M. Strom116 , R. Stroynowski40 ,A. Strubig106 , S.A. Stucci17 , B. Stugu14 , N.A. Styles42 , D. Su144 , J. Su125 , R. Subramaniam79 ,A. Succurro12 , Y. Sugaya118 , C. Suhr108 , M. Suk128 , V.V. Sulin96 , S. Sultansoy4d , T. Sumida68 ,S. Sun57 , X. Sun33a , J.E. Sundermann48 , K. Suruliz140 , G. Susinno37a,37b , M.R. Sutton150 ,Y. Suzuki66 , M. Svatos127 , S. Swedish169 , M. Swiatlowski144 , I. Sykora145a , T. Sykora129 , D. Ta90 ,C. Taccini135a,135b , K. Tackmann42 , J. Taenzer159 , A. Taffard164 , R. Tafirout160a , N. Taiblum154 ,H. Takai25 , R. Takashima69 , H. Takeda67 , T. Takeshita141 , Y. Takubo66 , M. Talby85 ,– 80 –JHEP01(2015)069

A.A. Talyshev109,c , J.Y.C. Tam175 , K.G. Tan88 , J. Tanaka156 , R. Tanaka117 , S. Tanaka132 ,S. Tanaka66 , A.J. Tanasijczuk143 , B.B. Tannenwald111 , N. Tannoury21 , S. Tapprogge83 ,S. Tarem153 , F. Tarrade29 , G.F. Tartarelli91a , P. Tas129 , M. Tasevsky127 , T. Tashiro68 ,E. Tassi37a,37b , A. Tavares Delgado126a,126b , Y. Tayalati136d , F.E. Taylor94 , G.N. Taylor88 ,W. Taylor160b , F.A. Teischinger30 , M. Teixeira Dias Castanheira76 , P. Teixeira-Dias77 ,K.K. Temming48 , H. Ten Kate30 , P.K. Teng152 , J.J. Teoh118 , S. Terada66 , K. Terashi156 ,J. Terron82 , S. Terzo101 , M. Testa47 , R.J. Teuscher159,j , J. Therhaag21 , T. Theveneaux-Pelzer34 ,J.P. Thomas18 , J. Thomas-Wilsker77 , E.N. Thompson35 , P.D. Thompson18 , P.D. Thompson159 ,R.J. Thompson84 , A.S. Thompson53 , L.A. Thomsen36 , E. Thomson122 , M. Thomson28 ,W.M. Thong88 , R.P. Thun89,∗ , F. Tian35 , M.J. Tibbetts15 , V.O. Tikhomirov96,af ,Yu.A. Tikhonov109,c , S. Timoshenko98 , E. Tiouchichine85 , P. Tipton177 , S. Tisserant85 ,T. Todorov5 , S. Todorova-Nova129 , J. Tojo70 , S. Tokár145a , K. Tokushuku66 , K. Tollefson90 ,E. Tolley57 , L. Tomlinson84 , M. Tomoto103 , L. Tompkins31 , K. Toms105 , N.D. Topilin65 ,E. Torrence116 , H. Torres143 , E. Torró Pastor168 , J. Toth85,ag , F. Touchard85 , D.R. Tovey140 ,H.L. Tran117 , T. Trefzger175 , L. Tremblet30 , A. Tricoli30 , I.M. Trigger160a , S. Trincaz-Duvoid80 ,M.F. Tripiana12 , W. Trischuk159 , B. Trocmé55 , C. Troncon91a , M. Trottier-McDonald15 ,M. Trovatelli135a,135b , P. True90 , M. Trzebinski39 , A. Trzupek39 , C. Tsarouchas30 , J.C-L. Tseng120 ,P.V. Tsiareshka92 , D. Tsionou137 , G. Tsipolitis10 , N. Tsirintanis9 , S. Tsiskaridze12 ,V. Tsiskaridze48 , E.G. Tskhadadze51a , I.I. Tsukerman97 , V. Tsulaia15 , S. Tsuno66 ,D. Tsybychev149 , A. Tudorache26a , V. Tudorache26a , A.N. Tuna122 , S.A. Tupputi20a,20b ,S. Turchikhin99,ae , D. Turecek128 , I. Turk Cakir4c , R. Turra91a,91b , A.J. Turvey40 , P.M. Tuts35 ,A. Tykhonov49 , M. Tylmad147a,147b , M. Tyndel131 , K. Uchida21 , I. Ueda156 , R. Ueno29 ,M. Ughetto85 , M. Ugland14 , M. Uhlenbrock21 , F. Ukegawa161 , G. Unal30 , A. Undrus25 ,G. Unel164 , F.C. Ungaro48 , Y. Unno66 , C. Unverdorben100 , D. Urbaniec35 , P. Urquijo88 , G. Usai8 ,A. Usanova62 , L. Vacavant85 , V. Vacek128 , B. Vachon87 , N. Valencic107 , S. Valentinetti20a,20b ,A. Valero168 , L. Valery34 , S. Valkar129 , E. Valladolid Gallego168 , S. Vallecorsa49 ,J.A. Valls Ferrer168 , W. Van Den Wollenberg107 , P.C. Van Der Deijl107 , R. van der Geer107 ,H. van der Graaf107 , R. Van Der Leeuw107 , D. van der Ster30 , N. van Eldik30 , P. van Gemmeren6 ,J. Van Nieuwkoop143 , I. van Vulpen107 , M.C. van Woerden30 , M. Vanadia133a,133b , W. Vandelli30 ,R. Vanguri122 , A. Vaniachine6 , P. Vankov42 , F. Vannucci80 , G. Vardanyan178 , R. Vari133a ,E.W. Varnes7 , T. Varol86 , D. Varouchas80 , A. Vartapetian8 , K.E. Varvell151 , F. Vazeille34 ,T. Vazquez Schroeder54 , J. Veatch7 , F. Veloso126a,126c , S. Veneziano133a , A. Ventura73a,73b ,D. Ventura86 , M. Venturi170 , N. Venturi159 , A. Venturini23 , V. Vercesi121a , M. Verducci133a,133b ,W. Verkerke107 , J.C. Vermeulen107 , A. Vest44 , M.C. Vetterli143,e , O. Viazlo81 , I. Vichou166 ,T. Vickey146c,ah , O.E. Vickey Boeriu146c , G.H.A. Viehhauser120 , S. Viel169 , R. Vigne30 ,M. Villa20a,20b , M. Villaplana Perez91a,91b , E. Vilucchi47 , M.G. Vincter29 , V.B. Vinogradov65 ,J. Virzi15 , I. Vivarelli150 , F. Vives Vaque3 , S. Vlachos10 , D. Vladoiu100 , M. Vlasak128 , A. Vogel21 ,M. Vogel32a , P. Vokac128 , G. Volpi124a,124b , M. Volpi88 , H. von der Schmitt101 ,H. von Radziewski48 , E. von Toerne21 , V. Vorobel129 , K. Vorobev98 , M. Vos168 , R. Voss30 ,J.H. Vossebeld74 , N. Vranjes137 , M. Vranjes Milosavljevic13a , V. Vrba127 , M. Vreeswijk107 ,T. Vu Anh48 , R. Vuillermet30 , I. Vukotic31 , Z. Vykydal128 , P. Wagner21 , W. Wagner176 ,H. Wahlberg71 , S. Wahrmund44 , J. Wakabayashi103 , J. Walder72 , R. Walker100 , W. Walkowiak142 ,R. Wall177 , P. Waller74 , B. Walsh177 , C. Wang152,ai , C. Wang45 , F. Wang174 , H. Wang15 ,H. Wang40 , J. Wang42 , J. Wang33a , K. Wang87 , R. Wang105 , S.M. Wang152 , T. Wang21 ,X. Wang177 , C. Wanotayaroj116 , A. Warburton87 , C.P. Ward28 , D.R. Wardrope78 ,M. Warsinsky48 , A. Washbrook46 , C. Wasicki42 , P.M. Watkins18 , A.T. Watson18 , I.J. Watson151 ,M.F. Watson18 , G. Watts139 , S. Watts84 , B.M. Waugh78 , S. Webb84 , M.S. Weber17 ,S.W. Weber175 , J.S. Webster31 , A.R. Weidberg120 , B. Weinert61 , J. Weingarten54 , C. Weiser48 ,1

Department of Physics, University of Adelaide, Adelaide, Australia
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, Gazi University,
Ankara; (c) Istanbul Aydin University, Istanbul; (d) Division of Physics, TOBB University of
Economics and Technology, Ankara, Turkey
5
LAPP, CNRS/IN2P3 and Université 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 ofAmerica9
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 Autònoma deBarcelona, 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
2– 81 –JHEP01(2015)069

H. Weits107 , P.S. Wells30 , T. Wenaus25 , D. Wendland16 , Z. Weng152,ad , T. Wengler30 , S. Wenig30 ,N. Wermes21 , M. Werner48 , P. Werner30 , M. Wessels58a , J. Wetter162 , K. Whalen29 , A. White8 ,M.J. White1 , R. White32b , S. White124a,124b , D. Whiteson164 , D. Wicke176 , F.J. Wickens131 ,W. Wiedenmann174 , M. Wielers131 , P. Wienemann21 , C. Wiglesworth36 , L.A.M. Wiik-Fuchs21 ,P.A. Wijeratne78 , A. Wildauer101 , M.A. Wildt42,aj , H.G. Wilkens30 , H.H. Williams122 ,S. Williams28 , C. Willis90 , S. Willocq86 , A. Wilson89 , J.A. Wilson18 , I. Wingerter-Seez5 ,F. Winklmeier116 , B.T. Winter21 , M. Wittgen144 , T. Wittig43 , J. Wittkowski100 , S.J. Wollstadt83 ,M.W. Wolter39 , H. Wolters126a,126c , B.K. Wosiek39 , J. Wotschack30 , M.J. Woudstra84 ,K.W. Wozniak39 , M. Wright53 , M. Wu55 , S.L. Wu174 , X. Wu49 , Y. Wu89 , E. Wulf35 ,T.R. Wyatt84 , B.M. Wynne46 , S. Xella36 , M. Xiao137 , D. Xu33a , L. Xu33b,ak , B. Yabsley151 ,S. Yacoob146b,al , R. Yakabe67 , M. Yamada66 , H. Yamaguchi156 , Y. Yamaguchi118 ,A. Yamamoto66 , K. Yamamoto64 , S. Yamamoto156 , T. Yamamura156 , T. Yamanaka156 ,K. Yamauchi103 , Y. Yamazaki67 , Z. Yan22 , H. Yang33e , H. Yang174 , U.K. Yang84 , Y. Yang111 ,S. Yanush93 , L. Yao33a , W-M. Yao15 , Y. Yasu66 , E. Yatsenko42 , K.H. Yau Wong21 , J. Ye40 ,S. Ye25 , I. Yeletskikh65 , A.L. Yen57 , E. Yildirim42 , M. Yilmaz4b , R. Yoosoofmiya125 , K. Yorita172 ,R. Yoshida6 , K. Yoshihara156 , C. Young144 , C.J.S. Young30 , S. Youssef22 , D.R. Yu15 , J. Yu8 ,J.M. Yu89 , J. Yu114 , L. Yuan67 , A. Yurkewicz108 , I. Yusuff28,am , B. Zabinski39 , R. Zaidan63 ,A.M. Zaitsev130,z , A. Zaman149 , S. Zambito23 , L. Zanello133a,133b , D. Zanzi88 , C. Zeitnitz176 ,M. Zeman128 , A. Zemla38a , K. Zengel23 , O. Zenin130 , T. Ženiš145a , D. Zerwas117 ,G. Zevi della Porta57 , D. Zhang89 , F. Zhang174 , H. Zhang90 , J. Zhang6 , L. Zhang152 , X. Zhang33d ,Z. Zhang117 , Y. Zhao33d , Z. Zhao33b , A. Zhemchugov65 , J. Zhong120 , B. Zhou89 , L. Zhou35 ,N. Zhou164 , C.G. Zhu33d , H. Zhu33a , J. Zhu89 , Y. Zhu33b , X. Zhuang33a , K. Zhukov96 ,A. Zibell175 , D. Zieminska61 , N.I. Zimine65 , C. Zimmermann83 , R. Zimmermann21 ,S. Zimmermann21 , S. Zimmermann48 , Z. Zinonos54 , M. Ziolkowski142 , G. Zobernig174 ,A. Zoccoli20a,20b , M. zur Nedden16 , G. Zurzolo104a,104b , V. Zutshi108 , L. Zwalinski30 .

18

School of Physics and Astronomy, University of Birmingham, Birmingham, United Kingdom
Department of Physics, Bogazici University, Istanbul; (b) Department of Physics, Dogus
University, Istanbul; (c) Department of Physics Engineering, Gaziantep University, Gaziantep,
Turkey20 (a)INFN Sezione di Bologna; (b) Dipartimento di Fisica e Astronomia, Università 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 Universityof Juiz de Fora (UFJF), Juiz de Fora; (c) Federal University of Sao Joao del Rei (UFSJ), Sao Joaodel 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) National Institute for
Research and Development of Isotopic and Molecular Technologies, Physics Department, Cluj
Napoca; (c) University Politehnica Bucharest, Bucharest; (d) 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 Católica de Chile, Santiago; (b) Departamento deFı́sica, Universidad Técnica Federico Santa Marı́a, Valparaı́so, Chile
33 (a)
Institute of High Energy Physics, Chinese Academy of Sciences, Beijing; (b) Department ofModern Physics, University of Science and Technology of China, Anhui; (c) Department of Physics,
Nanjing University, Jiangsu; (d) School of Physics, Shandong University, Shandong; (e) Physics
Department, Shanghai Jiao Tong University, Shanghai; (f ) Physics Department, Tsinghua
University, Beijing 100084, China
34
Laboratoire de Physique Corpusculaire, Clermont Université and Université Blaise Pascal and
CNRS/IN2P3, 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, Laboratori Nazionali di Frascati; (b) Dipartimento di Fisica,
Università della Calabria, Rende, Italy
38 (a)AGH University of Science and Technology, Faculty of Physics and Applied Computer Science,
Krakow; (b) Marian Smoluchowski Institute of Physics, Jagiellonian University, Krakow, Poland
39
The Henryk Niewodniczanski Institute of Nuclear Physics, Polish Academy of Sciences, Krakow,
Poland40
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 für Experimentelle Physik IV, Technische Universität Dortmund, Dortmund, Germany
44
Institut für Kern- und Teilchenphysik, Technische Universität 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
Fakultät für Mathematik und Physik, Albert-Ludwigs-Universität, Freiburg, Germany
49
Section de Physique, Université de Genève, Geneva, Switzerland
50 (a)INFN Sezione di Genova; (b) Dipartimento di Fisica, Università di Genova, Genova, Italy
51 (a)
E. Andronikashvili Institute of Physics, Iv. Javakhishvili Tbilisi State University, Tbilisi; (b) HighEnergy Physics Institute, Tbilisi State University, Tbilisi, Georgia
19 (a)JHEP01(2015)069– 82 –

52

II Physikalisches Institut, Justus-Liebig-Universität Giessen, Giessen, Germany
SUPA - School of Physics and Astronomy, University of Glasgow, Glasgow, United Kingdom
54
II Physikalisches Institut, Georg-August-Universität, Göttingen, Germany
55
Laboratoire de Physique Subatomique et de Cosmologie, Université Grenoble-Alpes, CNRS/IN2P3,
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 Statesof America
58 (a)
Kirchhoff-Institut für Physik, Ruprecht-Karls-Universität Heidelberg, Heidelberg;(b)Physikalisches Institut, Ruprecht-Karls-Universität Heidelberg, Heidelberg; (c) ZITI Institut für
technische Informatik, Ruprecht-Karls-Universität Heidelberg, Mannheim, Germany
59
Faculty of Applied Information Science, Hiroshima Institute of Technology, Hiroshima, Japan
60 (a)
Department of Physics, The Chinese University of Hong Kong, Shatin, N.T., Hong Kong;(b)
Department of Physics, The University of Hong Kong, Hong Kong; (c) Department of Physics,
The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong,
China
61
Department of Physics, Indiana University, Bloomington IN, United States of America
62
Institut für Astro- und Teilchenphysik, Leopold-Franzens-Universität, Innsbruck, Austria
63
University of Iowa, Iowa City IA, United States of America
64
Department of Physics and Astronomy, Iowa State University, Ames IA, United States of America
65
Joint Institute for Nuclear Research, JINR Dubna, Dubna, Russia
66
KEK, High Energy Accelerator Research Organization, Tsukuba, Japan
67
Graduate School of Science, Kobe University, Kobe, Japan
68
Faculty of Science, Kyoto University, Kyoto, Japan
69
Kyoto University of Education, Kyoto, Japan70
Department of Physics, Kyushu University, Fukuoka, Japan
71
Instituto de Fı́sica La Plata, Universidad Nacional de La Plata and CONICET, La Plata, Argentina
72
Physics Department, Lancaster University, Lancaster, United Kingdom
73 (a)INFN Sezione di Lecce; (b) Dipartimento di Matematica e Fisica, Università del Salento, Lecce,
Italy
74
Oliver Lodge Laboratory, University of Liverpool, Liverpool, United Kingdom
75
Department of Physics, Jožef Stefan Institute and University of Ljubljana, Ljubljana, Slovenia
76
School of Physics and Astronomy, Queen Mary University of London, London, United Kingdom
77
Department of Physics, Royal Holloway University of London, Surrey, United Kingdom
78
Department of Physics and Astronomy, University College London, London, United Kingdom
79
Louisiana Tech University, Ruston LA, United States of America
80
Laboratoire de Physique Nucléaire et de Hautes Energies, UPMC and Université Paris-Diderot and
CNRS/IN2P3, Paris, France
81
Fysiska institutionen, Lunds universitet, Lund, Sweden
82
Departamento de Fisica Teorica C-15, Universidad Autonoma de Madrid, Madrid, Spain
83
Institut für Physik, Universität Mainz, Mainz, Germany
84
School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom
85
CPPM, Aix-Marseille Université and CNRS/IN2P3, Marseille, France
86
Department of Physics, University of Massachusetts, Amherst MA, United States of America
87
Department of Physics, McGill University, Montreal QC, Canada
88
School of Physics, University of Melbourne, Victoria, Australia
89
Department of Physics, The University of Michigan, Ann Arbor MI, United States of America
90
Department of Physics and Astronomy, Michigan State University, East Lansing MI, United Statesof America
91 (a)INFN Sezione di Milano; (b) Dipartimento di Fisica, Università di Milano, Milano, Italy
92
B.I. Stepanov Institute of Physics, National Academy of Sciences of Belarus, Minsk, Republic of
Belarus53

JHEP01(2015)069– 83 –

93– 84 –JHEP01(2015)069

National Scientific and Educational Centre for Particle and High Energy Physics, Minsk, Republicof Belarus
94
Department of Physics, Massachusetts Institute of Technology, Cambridge MA, United States ofAmerica
95
Group of Particle Physics, University of Montreal, Montreal QC, Canada
96
P.N. Lebedev Institute of Physics, Academy of Sciences, Moscow, Russia
97
Institute for Theoretical and Experimental Physics (ITEP), Moscow, Russia
98
National Research Nuclear University MEPhI, Moscow, Russia
99
D.V.Skobeltsyn Institute of Nuclear Physics, M.V.Lomonosov Moscow State University, Moscow,
Russia100
Fakultät für Physik, Ludwig-Maximilians-Universität München, München, Germany
101
Max-Planck-Institut für Physik (Werner-Heisenberg-Institut), München, Germany
102
Nagasaki Institute of Applied Science, Nagasaki, Japan103
Graduate School of Science and Kobayashi-Maskawa Institute, Nagoya University, Nagoya, Japan
104 (a)INFN Sezione di Napoli; (b) Dipartimento di Fisica, Università di Napoli, Napoli, Italy
105
Department of Physics and Astronomy, University of New Mexico, Albuquerque NM, United Statesof America
106
Institute for Mathematics, Astrophysics and Particle Physics, Radboud University
Nijmegen/Nikhef, Nijmegen, Netherlands
107
Nikhef National Institute for Subatomic Physics and University of Amsterdam, Amsterdam,
Netherlands
108
Department of Physics, Northern Illinois University, DeKalb IL, United States of America
109
Budker Institute of Nuclear Physics, SB RAS, Novosibirsk, Russia
110
Department of Physics, New York University, New York NY, United States of America
111
Ohio State University, Columbus OH, United States of America
112
Faculty of Science, Okayama University, Okayama, Japan
113
Homer L. Dodge Department of Physics and Astronomy, University of Oklahoma, Norman OK,
United States of America
114
Department of Physics, Oklahoma State University, Stillwater OK, United States of America
115
Palacký University, RCPTM, Olomouc, Czech Republic
116
Center for High Energy Physics, University of Oregon, Eugene OR, United States of America
117
LAL, Université Paris-Sud and CNRS/IN2P3, Orsay, France
118
Graduate School of Science, Osaka University, Osaka, Japan
119
Department of Physics, University of Oslo, Oslo, Norway
120
Department of Physics, Oxford University, Oxford, United Kingdom
121 (a)INFN Sezione di Pavia; (b) Dipartimento di Fisica, Università di Pavia, Pavia, Italy
122
Department of Physics, University of Pennsylvania, Philadelphia PA, United States of America
123
Petersburg Nuclear Physics Institute, Gatchina, Russia
124 (a)INFN Sezione di Pisa; (b) Dipartimento di Fisica E. Fermi, Università di Pisa, Pisa, Italy
125
Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh PA, United States ofAmerica
126 (a)Laboratorio de Instrumentacao e Fisica Experimental de Particulas - LIP, Lisboa; (b) Faculdadede Ciências, Universidade de Lisboa, Lisboa; (c) Department of Physics, University of Coimbra,
Coimbra; (d) Centro de Fı́sica Nuclear da Universidade de Lisboa, Lisboa; (e) Departamento de
Fisica, Universidade do Minho, Braga; (f ) Departamento de Fisica Teorica y del Cosmos and
CAFPE, Universidad de Granada, Granada (Spain); (g) Dep Fisica and CEFITEC of Faculdade deCiencias e Tecnologia, Universidade Nova de Lisboa, Caparica, Portugal
127
Institute of Physics, Academy of Sciences of the Czech Republic, Praha, Czech Republic
128
Czech Technical University in Prague, Praha, Czech Republic
129
Faculty of Mathematics and Physics, Charles University in Prague, Praha, Czech Republic
130State Research Center Institute for High Energy Physics, Protvino, Russia
131
Particle Physics Department, Rutherford Appleton Laboratory, Didcot, United Kingdom

132

Ritsumeikan University, Kusatsu, Shiga, JapanINFN Sezione di Roma; (b) Dipartimento di Fisica, Sapienza Università di Roma, Roma, Italy
134 (a)INFN Sezione di Roma Tor Vergata; (b) Dipartimento di Fisica, Università di Roma Tor Vergata,
Roma, Italy
135 (a)INFN Sezione di Roma Tre; (b) Dipartimento di Matematica e Fisica, Università Roma Tre,
Roma, Italy
136 (a)Faculté des Sciences Ain Chock, Réseau Universitaire de Physique des Hautes Energies Université Hassan II, Casablanca; (b) Centre National de l’Energie des Sciences Techniques
Nucleaires, Rabat; (c) Faculté des Sciences Semlalia, Université Cadi Ayyad, LPHEA-Marrakech;(d)Faculté des Sciences, Université Mohamed Premier and LPTPM, Oujda; (e) Faculté des sciences,
Université Mohammed V-Agdal, Rabat, Morocco137
DSM/IRFU (Institut de Recherches sur les Lois Fondamentales de l’Univers), CEA Saclay(Commissariat à l’Energie Atomique et aux Energies Alternatives), Gif-sur-Yvette, France
138
Santa Cruz Institute for Particle Physics, University of California Santa Cruz, Santa Cruz CA,
United States of America
139
Department of Physics, University of Washington, Seattle WA, United States of America
140
Department of Physics and Astronomy, University of Sheffield, Sheffield, United Kingdom
141
Department of Physics, Shinshu University, Nagano, Japan142
Fachbereich Physik, Universität Siegen, Siegen, Germany
143
Department of Physics, Simon Fraser University, Burnaby BC, Canada
144
SLAC National Accelerator Laboratory, Stanford CA, United States of America
145 (a)Faculty of Mathematics, Physics & Informatics, Comenius University, Bratislava; (b) Departmentof Subnuclear Physics, Institute of Experimental Physics of the Slovak Academy of Sciences,
Kosice, Slovak Republic
146 (a)Department of Physics, University of Cape Town, Cape Town; (b) Department of Physics,
University of Johannesburg, Johannesburg; (c) School of Physics, University of the Witwatersrand,
Johannesburg, South Africa
147 (a)Department of Physics, Stockholm University; (b) The Oskar Klein Centre, Stockholm, Sweden
148
Physics Department, Royal Institute of Technology, Stockholm, Sweden
149
Departments of Physics & Astronomy and Chemistry, Stony Brook University, Stony Brook NY,
United States of America
150
Department of Physics and Astronomy, University of Sussex, Brighton, United Kingdom
151
School of Physics, University of Sydney, Sydney, Australia
152
Institute of Physics, Academia Sinica, Taipei, Taiwan
153
Department of Physics, Technion: Israel Institute of Technology, Haifa, Israel
154
Raymond and Beverly Sackler School of Physics and Astronomy, Tel Aviv University, Tel Aviv,
Israel155
Department of Physics, Aristotle University of Thessaloniki, Thessaloniki, Greece
156
International Center for Elementary Particle Physics and Department of Physics, The Universityof Tokyo, Tokyo, Japan
157
Graduate School of Science and Technology, Tokyo Metropolitan University, Tokyo, Japan
158
Department of Physics, Tokyo Institute of Technology, Tokyo, Japan
159
Department of Physics, University of Toronto, Toronto ON, Canada
160 (a)TRIUMF, Vancouver BC; (b) Department of Physics and Astronomy, York University, Toronto
ON, Canada
161
Faculty of Pure and Applied Sciences, University of Tsukuba, Tsukuba, Japan
162
Department of Physics and Astronomy, Tufts University, Medford MA, United States of America
163
Centro de Investigaciones, Universidad Antonio Narino, Bogota, Colombia
164
Department of Physics and Astronomy, University of California Irvine, Irvine CA, United States ofAmerica165 (a)INFN Gruppo Collegato di Udine, Sezione di Trieste, Udine; (b) ICTP, Trieste; (c) Dipartimento
di Chimica, Fisica e Ambiente, Università di Udine, Udine, Italy
133 (a)JHEP01(2015)069– 85 –

166

Department of Physics, University of Illinois, Urbana IL, United States of America
Department of Physics and Astronomy, University of Uppsala, Uppsala, Sweden
168
Instituto de Fı́sica Corpuscular (IFIC) and Departamento de Fı́sica Atómica, Molecular y Nuclear
and Departamento de Ingenierı́a Electrónica and Instituto de Microelectrónica de Barcelona(IMB-CNM), University of Valencia and CSIC, Valencia, Spain
169
Department of Physics, University of British Columbia, Vancouver BC, Canada
170
Department of Physics and Astronomy, University of Victoria, Victoria BC, Canada
171
Department of Physics, University of Warwick, Coventry, United Kingdom
172
Waseda University, Tokyo, Japan173
Department of Particle Physics, The Weizmann Institute of Science, Rehovot, Israel
174
Department of Physics, University of Wisconsin, Madison WI, United States of America
175
Fakultät für Physik und Astronomie, Julius-Maximilians-Universität, Würzburg, Germany
176
Fachbereich C Physik, Bergische Universität Wuppertal, Wuppertal, Germany
177
Department of Physics, Yale University, New Haven CT, United States of America
178
Yerevan Physics Institute, Yerevan, Armenia
179
Centre de Calcul de l’Institut National de Physique Nucléaire et de Physique des Particules(IN2P3), Villeurbanne, FranceaAlso at Department of Physics, King’s College London, London, United Kingdomb
Also at Institute of Physics, Azerbaijan Academy of Sciences, Baku, AzerbaijancAlso at Novosibirsk State University, Novosibirsk, RussiadAlso at Particle Physics Department, Rutherford Appleton Laboratory, Didcot, United KingdomeAlso at TRIUMF, Vancouver BC, CanadafAlso at Department of Physics, California State University, Fresno CA, United States of AmericagAlso at Tomsk State University, Tomsk, RussiahAlso at CPPM, Aix-Marseille Université and CNRS/IN2P3, Marseille, Francei
Also at Università di Napoli Parthenope, Napoli, Italyj
Also at Institute of Particle Physics (IPP), CanadakAlso at Department of Physics, St. Petersburg State Polytechnical University, St. Petersburg,
Russial
Also at Department of Financial and Management Engineering, University of the Aegean, Chios,
GreecemAlso at Louisiana Tech University, Ruston LA, United States of AmericanAlso at Institucio Catalana de Recerca i Estudis Avancats, ICREA, Barcelona, Spaino
Also at Department of Physics, The University of Texas at Austin, Austin TX, United States ofAmericap
Also at Institute of Theoretical Physics, Ilia State University, Tbilisi, Georgiaq
Also at CERN, Geneva, Switzerlandr
Also at Ochadai Academic Production, Ochanomizu University, Tokyo, Japans
Also at Manhattan College, New York NY, United States of AmericatAlso at Institute of Physics, Academia Sinica, Taipei, TaiwanuAlso at LAL, Université Paris-Sud and CNRS/IN2P3, Orsay, Francev
Also at Academia Sinica Grid Computing, Institute of Physics, Academia Sinica, Taipei, Taiwan
w
Also at Laboratoire de Physique Nucléaire et de Hautes Energies, UPMC and Université
Paris-Diderot and CNRS/IN2P3, Paris, Francex
Also at School of Physical Sciences, National Institute of Science Education and Research,
Bhubaneswar, Indiay
Also at Dipartimento di Fisica, Sapienza Università di Roma, Roma, Italyz
Also at Moscow Institute of Physics and Technology State University, Dolgoprudny, Russia
aaAlso at section de Physique, Université de Genève, Geneva, Switzerlandab
Also at International School for Advanced Studies (SISSA), Trieste, Italyac
Also at Department of Physics and Astronomy, University of South Carolina, Columbia SC, UnitedStates of America
167

JHEP01(2015)069– 86 –

ad

Also at School of Physics and Engineering, Sun Yat-sen University, Guangzhou, China
Also at Faculty of Physics, M.V.Lomonosov Moscow State University, Moscow, Russia
af
Also at National Research Nuclear University MEPhI, Moscow, Russiaag
Also at Institute for Particle and Nuclear Physics, Wigner Research Centre for Physics, Budapest,
HungaryahAlso at Department of Physics, Oxford University, Oxford, United KingdomaiAlso at Department of Physics, Nanjing University, Jiangsu, China
aj
Also at Institut für Experimentalphysik, Universität Hamburg, Hamburg, GermanyakAlso at Department of Physics, The University of Michigan, Ann Arbor MI, United States ofAmericaal
Also at Discipline of Physics, University of KwaZulu-Natal, Durban, South AfricaamAlso at University of Malaya, Department of Physics, Kuala Lumpur, Malaysia ∗ Deceasedae

JHEP01(2015)069– 87 –