J H E P 1 0 ( 2 0 1 3 ) 1 6 7 Published for SISSA by Springer Received: August 18, 2013 Accepted: September 16, 2013 Published: October 24, 2013 Measurement of the W-boson helicity in top-quark decays from tt production in lepton+jets events in pp collisions at √ s = 7TeV The CMS collaboration E-mail: cms-publication-committee-chair@cern.ch Abstract: The W-boson helicity fractions in top-quark decays are measured with tt events in the lepton+jets final state, using proton-proton collisions at a centre-of-mass energy of 7 TeV, collected in 2011 with the CMS detector at the LHC. The data sample corresponds to an integrated luminosity of 5.0 fb−1. The measured fractions of longitudinal, left-, and right-handed helicity are F0 = 0.682± 0.030 (stat.)± 0.033 (syst.), FL = 0.310± 0.022 (stat.) ± 0.022 (syst.), and FR = 0.008 ± 0.012 (stat.) ± 0.014 (syst.), consistent with the standard model predictions. The measured fractions are used to probe the existence of anomalous Wtb couplings. Exclusion limits on the real components of the anomalous couplings gL, gR are also derived. Keywords: Hadron-Hadron Scattering, Top physics Open Access, Copyright CERN, for the benefit of the CMS collaboration doi:10.1007/JHEP10(2013)167 mailto:cms-publication-committee-chair@cern.ch http://dx.doi.org/10.1007/JHEP10(2013)167 J H E P 1 0 ( 2 0 1 3 ) 1 6 7 Contents 1 Introduction 1 2 The CMS detector 3 3 Data and simulated samples 4 4 Event reconstruction and selection 4 5 Top-quark reconstruction 7 6 Background estimation 8 6.1 Normalisation and shape of DY production with jets 8 6.2 Normalisation and shape of W-boson production with jets 9 7 Determination of helicity fractions 9 8 Comparison between data and SM predictions 11 9 Systematic uncertainties 12 10 Results 17 11 Limits on anomalous couplings 18 12 Summary 22 The CMS collaboration 27 1 Introduction Following the discovery of the top quark in proton-antiproton collisions at the Tevatron collider [1, 2], measurements of W-boson helicity fractions in top-quark decays have been an important subject of investigation, because of their relationship to the V−A structure of the weak charged current and their sensitivity to physics beyond the standard model (SM). With the large samples of tt events produced in proton-proton collisions at the Large Hadron Collider (LHC), the W-boson helicity fraction measurements can be improved considerably, enhancing the search for anomalous Wtb couplings, i.e. those that do not arise from the SM. Previous measurements of W-boson helicity fractions in top-quark decays have been performed by the CDF, D0 [3], and ATLAS [4] Collaborations. The helicity fractions of the W boson produced in a t→Wb decay are defined as the partial rate for a given helicity state divided by the total decay rate: FL,R,0 ≡ ΓL,R,0/Γ, – 1 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 where FL, FR, and F0 are the left-handed, right-handed, and longitudinal helicity fractions, respectively. For SM couplings and unpolarised top-quark production, the helicity fractions are approximately 70% longitudinal and 30% left-handed. At leading order (LO) and in the limit mb = 0 (where mb is the b-quark mass), the right-handed helicity fraction is zero due to helicity suppression. For finite mb, the helicity fractions are [5]: F0 = (1− y2)2 − x2(1 + y2) (1− y2)2 + x2(1− 2x2 + y2) , (1.1) FL = x2(1− x2 + y2 + √ λ) (1− y2)2 + x2(1− 2x2 + y2) , (1.2) FR = x2(1− x2 + y2 − √ λ) (1− y2)2 + x2(1− 2x2 + y2) , (1.3) where x = MW/mt, y = mb/mt, λ = 1 +x4 + y4− 2x2y2− 2x2− 2y2, and MW, mt are the masses of the W boson and top quark, respectively. These fractions are minimally modified by higher-order corrections. Recent next-to-next-to-leading-order (NNLO) calculations [6] predict F0 = 0.687± 0.005, FL = 0.311± 0.005, and FR = 0.0017± 0.0001 for a top-quark mass of mt = 172.8± 1.3 GeV/c2. Experimentally, the W-boson helicity components can be extracted through the study of angular distributions of top-quark decay products in tt final states. The helicity angle θ∗ is defined as the angle between the W-boson momentum in the top-quark rest frame and the momentum of the down-type decay fermion in the rest frame of the W boson. The distribu- tion of the cosine of the helicity angle has a dependence on the helicity fractions given by: 1 Γ dΓ d cos θ∗ = 3 8 FL (1− cos θ∗)2 + 3 4 F0(sin θ ∗)2 + 3 8 FR (1 + cos θ∗)2 . (1.4) Deviations of the measured helicity fractions from the SM predictions can be inter- preted in terms of anomalous Wtb couplings [7, 8], using the most general dimension-six Lagrangian: LWtb = − g√ 2 b̄γµ(VLPL + VRPR)tW−µ − g√ 2 b̄ iσµνqν MW (gLPL + gRPR) t W−µ + h.c., (1.5) where VL, VR, gL, gR are dimensionless complex constants, q = pt − pb, where pt (pb) is the four-momentum of the top quark (b quark), PL (PR) are the left (right) projec- tor operators, and h.c. denotes the Hermitian conjugate. Hermiticity conditions on the possible dimension-six Lagrangians also impose Im(VL) = 0 [8]. In the SM and at tree level, VL = Vtb, where Vtb is the Cabibbo-Kobayashi-Maskawa matrix element Vtb ' 1, and VR = gL = gR = 0. The relation between helicity fractions and anomalous couplings, including dependencies on the b-quark mass, are given in ref. [9]. We report on a study of the W-boson helicity fractions in top-quark decays using a sample of tt events where one of the top quarks decays semileptonically (e.g. t→W+b→ `+ν`b, where ` is either an electron or a muon) and the other decays hadronically (e.g. t→W−b→ qq′b). A kinematic fit is used to determine the best association of b jets, other jets, and lepton candidates to the top quark and antiquark decay hypotheses, interpreting – 2 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 the measured momentum imbalance as due to the presence of a neutrino. In this kinematic fit, top-quark and W-boson mass constraints are employed to improve the resolution of the measured jet and lepton energies, resulting in an improved reconstruction of the W-boson rest frame and the helicity angles in the weak decays of top quarks. In this article, the cosine of the helicity angles in the semileptonic and hadronic top-quark decays will be referred to as cos θ∗ and coshad θ∗, respectively. In the leptonic branch, the down-type fermion corresponds to the charged lepton. In the hadronic branch, since the down-type quark is not experimentally identified, only the absolute value of coshad θ∗ is used, providing information on the relative proportion between longitudinal (F0) and total transverse (FL+FR = 1−F0) fractions. The resulting helicity angle distributions are fitted to measure the W-boson helicity fractions and to determine possible anomalous Wtb couplings. 2 The CMS detector The central feature of the Compact Muon Solenoid (CMS) detector is a superconducting solenoid, 13 m in length and 6 m in diameter, which provides a uniform axial magnetic field of 3.8 T. The CMS experiment uses a right-handed coordinate system, with the origin at the nominal interaction point, the x axis pointing to the centre of the LHC, the y axis pointing up (perpendicular to the LHC plane), and the z axis along the anticlockwise-beam direction. The polar angle θ is measured from the positive z axis and the azimuthal angle φ is measured in the x-y plane. The bore of the solenoid is instrumented with various particle detection systems. Charged-particle trajectories are measured with a silicon pixel tracker with three barrel layers at radii between 4.4 and 10.2 cm, and a silicon strip tracker with 10 barrel de- tection layers that extend outward reaching a radius of 1.1 m. Each system is completed by two endcaps, and provides angular coverage of 0 < φ < 2π in azimuth and |η| < 2.5, where the pseudorapidity η is defined as η = − ln[tan(θ/2)], with θ the polar angle of the trajectory of the particle with respect to the z axis. A lead-tungstate crystal electromagnetic calorime- ter (ECAL) and a brass/scintillator hadron calorimeter (HCAL) surround the tracking vol- ume and cover the region |η| < 3. The forward calorimeter further extends the HCAL cov- erage in the region 3 < |η| < 5, improving the determination of the momentum imbalance. Muons are measured and identified in gas-ionisation detectors embedded in the steel flux return yoke outside the solenoid in the range |η| < 2.4. The barrel region is covered by drift-tube chambers and the endcap region by cathode strip chambers, each comple- mented by resistive plate chambers. The detector is nearly hermetic, allowing for energy balance measurements in the plane transverse to the beam. A two-level trigger system selects the most interesting pp collision events for use in physics analysis. The first level of the CMS trigger system, composed of custom hardware processors, uses information from the calorimeters and muon detectors to select the most interesting events in a fixed time interval of less than about 3.2µs. The second level, known as the high-level trigger (HLT), is a processor farm that further decreases the event rate from around 100 kHz to around 300 Hz, before data storage. A more detailed description of the CMS detector can be found elsewhere [10]. – 3 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 3 Data and simulated samples The measurements presented in this paper are performed using events from proton-proton collisions at a centre-of-mass energy of 7 TeV, collected by the CMS detector in 2011. The data correspond to an integrated luminosity of 5.0 ± 0.1 fb−1 [11]. In order to account for effects of detector resolution and acceptance, as well as to estimate the contribution from background processes that can satisfy the tt event selection criteria, simulated event samples based on Monte Carlo (MC) event generator programs are used. A tt signal sample was generated using the MadGraph generator version 5.1.1 [12] with matrix elements having up to three extra partons in the final state and an assumed top-quark mass of 172.5 GeV/c2. MadGraph is interfaced with pythia generator version 6.424 [13] to simulate hadronization and parton fragmentation and with tauola program version 27.121.5 [14] to simulate τ -lepton decays. The helicity fractions used as a SM reference are the LO predictions for mt = 172.5 GeV/c2: F0 = 0.6902, FL = 0.3089, FR = 0.0009, consistent with next-to-leading-order (NLO) and NNLO expectations [5, 6]. Single-top quark events in t and tW (or W-boson associated) channels were generated using powheg program version 1.0 [15] and pythia interfaced with tauola. Other rel- evant background processes, such as W boson and Drell-Yan production accompanied by multiple jets, were simulated using MadGraph. In all LO simulations, the parton distri- bution function (PDF) set CTEQ6L1 [16] is used. The powheg simulations use the NLO set CT10 [17]. The effect of multiple proton-proton collisions occurring within the same bunch crossing (pileup events) is taken into account in the simulation and matches the pileup distribution observed in data. 4 Event reconstruction and selection A set of requirements is applied to all samples, selecting candidate events compatible with the topology of tt production. Almost all top quarks decay into a W boson and a b quark. In the decay modes considered for this study, one of the W bosons decays hadronically into two jets and the other W boson decays leptonically into an electron or muon and a neutrino. Hence, final states containing a muon or electron and at least four jets are selected for further consideration. Top-quark decay products are reconstructed using the particle-flow (PF) algorithm de- scribed in detail in refs. [18] and [19]. The particle-flow event reconstruction identifies and measures the properties of each particle, using an optimised combination of the information from all subdetectors. The energy of photons is directly obtained from the ECAL measure- ment, corrected for zero-suppression effects. The energy of electrons is determined from a combination of the momentum of the track originated at the main interaction vertex, the corresponding ECAL cluster energy, and the energy sum of all bremsstrahlung pho- tons attached to the track. The energy of muons is obtained from the corresponding track momentum. The energy of charged hadrons is determined from a combination of the track momentum and the corresponding ECAL and HCAL energy, corrected for zero-suppression effects, and calibrated for the nonlinear response of the calorimeters. Finally, the energy of – 4 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 neutral hadrons is obtained from the corresponding calibrated ECAL and HCAL energy. The missing transverse energy Emiss T is defined as the magnitude of the transverse momen- tum imbalance, ~pmiss T , which is the negative of the vectorial sum of the transverse momenta pT of all the particles reconstructed with the particle-flow algorithm. Tracks belonging to the primary or secondary vertices of the most energetic pp interaction are retained, while particles identified as coming from pileup interactions are removed from the event. Events in the muon+jets channel were selected by a trigger that required at least one isolated, high-momentum muon with a HLT pT threshold varying between 17 and 24 GeV/c for the running period used in this analysis. Events containing a muon, well measured in the silicon tracker and identified in the muon chambers with pT > 25 GeV/c and pseudorapidity |η| < 2.1, are then selected offline. The track associated with the muon candidate is required to have a minimum number of hits in the silicon tracker, to be consistent with originating from the beam spot, and to have a high-quality global fit including a minimum number of hits in the muon detector. More details on the muon quality requirements are given in [20]. The trigger used to collect events in the electron+jets channel required at least one isolated electron with pT > 25 GeV/c, accompanied by at least three jets with pT > 30 GeV/c. Electron candidates are reconstructed from clusters of energy deposits in the electromagnetic calorimeter, which are then matched to hits in the silicon tracker [21]. Electrons are identified by using shower shape and track-cluster matching variables. Offline, events with exactly one electron candidate with pT > 30 GeV/c and |η| < 2.5 are selected. Events having electron candidates in the transition region between the barrel and endcap calorimeters, 1.44 < |η| < 1.56, are excluded because reconstruction in this region is degraded due to additional material there. The electron track must lie within 0.02 cm of the primary vertex in the plane transverse to the beam axis. Additionally, the background due to electrons from photon conversions is reduced by rejecting tracks with missing hits in the inner tracker layers or that are near a track with an opposite charge and a similar polar angle. Events with an additional muon with pT > 10 GeV/c or and additional electron with pT > 15 GeV/c are vetoed in order to reject backgrounds from dileptonic tt and Drell-Yan events. To reduce backgrounds further, muons and electrons are required to be prompt and are therefore typically well isolated from the rest of the event. This is achieved via an offline particle-flow-based relative isolation (PFIso) algorithm, which is defined as the sum of the transverse momenta over all charged hadrons, neutral hadrons, and photons reconstructed inside a cone of radius ∆R = √ (∆η)2 + (∆φ)2 = 0.3, centered around the lepton (muon or electron) and divided by the lepton transverse momentum. In top-quark decays, where the b jet coincidentally overlaps with the prompt lepton, cos θ∗ takes values close to −1. The offline lepton isolation requirement therefore has the effect of reducing the signal selection efficiency near cos θ∗ ≈ −1. In addition, both the muon and electron online triggers impose loose isolation criteria. Hence, selected muons (electrons) are required to have PFIso < 0.125 (0.100). These values are chosen to be tight enough to provide a high trigger efficiency, yet loose enough to maintain reasonable signal efficiency near cos θ∗ ≈ −1. In order to mitigate effects of pileup, a correction based on the average energy density in the event is applied to the electron isolation. – 5 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 Jets are reconstructed using the anti-kT clustering algorithm [22, 23], with a distance parameter of 0.5, applied to the entire list of reconstructed particles that are not identified as isolated muons or electrons in the event. The resulting uncalibrated jet momenta are found to be within 5% to 10% of the true momenta over the whole pT spectrum and detector acceptance. Charged particles not associated with the primary vertex are explicitly removed from the jet, as stated above. A residual correction is applied to account for the energy of any extra neutral particles arising from pileup interactions. Further residual jet energy calibrations are derived from simulation and corrected for any discrepancies with data using in situ measurements of object balancing in both dijet events (to render the jet response in pseudorapidity uniform) as well as photon+jet events (to provide the absolute jet energy scale) [24]. Additional selection criteria are applied to remove spurious events with identified noise patterns in the HCAL. Calibrated jets with pT > 10 GeV/c are used to correct the scale of Emiss T [25]. Jet candidates from the top-quark decay are required to have calibrated pT > 30 GeV/c and |η| < 2.4. Events with less than four jets passing the above mentioned top-quark decay product criteria are not used in the analysis. To reduce the QCD multijet background, the transverse mass, MT, of the lep- tonically decaying W boson, is required to be greater than 30 GeV/c2, where MT =√ 2pTE miss T (1− cos(∆φ))× 1/c3 and ∆φ is the angle in the x-y plane between the direction of the lepton and ~pmiss T . In events where the top-quark pair decays dileptonically and one lepton escapes detection, the MT variable can assume very high values. Background events from this process are rejected by requiring MT < 200 GeV/c2. Due to its relatively high rate and similar final state topology, the main remaining background source for this analysis is the production of several jets in association with a W boson that decays leptonically. This background source can be reduced, and the QCD multijet background even further suppressed, by requiring that at least two of the selected jets be identified as b jets. A high-efficiency tagging algorithm, known as the combined secondary-vertex (CSV) algorithm [26], is used to separate jets originating from light quarks (or gluons) and heavy quarks, i.e. charm or bottom quarks. Jets are first divided into cat- egories according to the probability of reconstructing a secondary vertex and its quality. Then, within each category, several variables including the three-dimensional signed im- pact parameter significance, secondary vertex mass, fractional charge, and charged particle multiplicity are used to form a likelihood that discriminates light-flavour jets from heavy- flavour jets. A selection working point is chosen so that the efficiency to identify a b jet is high (nearly 70%), while the probability that a light-flavour jet is mistaken as a b jet is small (about 1%). Requiring that there be two b-tagged jets in the event reduces the remaining QCD multijet background to negligible levels (less than 0.4%). Trigger, lepton identification, and lepton isolation efficiencies are estimated with a tag- and-probe method [27] using leptons from a sample of events containing Z→ `+`− decays. The efficiencies are computed for the Z-boson events in both data and simulation, as a function of the lepton pT and η. The overall efficiencies in the typical pT and η ranges of the selected leptons are 80% for muons and 70% for electrons. The ratio of the efficiencies in data and simulation, εDATA ` /εMC ` , is used as a scale factor to correct the simulated samples. Likewise, simulated samples are scaled to account for any differences in b-tagging – 6 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 efficiencies between data and simulation according to the ratio of efficiencies εDATA b-tag /ε MC b-tag, which is determined as a function of the pT of the b jet [26]. The number of data events reconstructed with these selection criteria is 9268 in the muon channel and 6526 in the electron channel. Comparisons between data observations and the SM expectations are presented in section 8. 5 Top-quark reconstruction Once events are selected according to the criteria described in section 4, the reconstruction of each top-quark candidate proceeds by testing all selected jets, ~pmiss T , and the lepton for their compatibility with decay products of the hadronic branch (t→ bW→ bqq̄′) and the leptonic branch (t → bW → b`ν). The initial value for the neutrino momentum is set to ~p ν = (~pmiss T , pνz), where pνz is determined by requiring that the invariant mass of the neutrino and lepton be equal to the W-boson mass, which is assumed to be 80.4 GeV/c2 [28]. For each possible neutrino solution and jet assignment to either the leptonic branch or hadronic branch, a χ2 is built, containing the following terms: χ2 comb = ( mt −mref t σmt )2 + ( mt −mref t σmt )2 + ( M lep W − 80.4 σ M lep W )2 + ( Mhad W − 80.4 σMhad W )2 (5.1) − ∑ i=1,4 2 ln pi(disc|f), where mt, mt, M lep W and Mhad W are the reconstructed invariant masses for a given combina- torial assignment of four jets to the final-state particles in the tt decay. The reference value for the top-quark mass mref t is taken to be 173.3 GeV/c2 [29] for data and 172.5 GeV/c2 for the simulated samples. The term pi(disc|f) is the probability for the ith jet to have flavour f (either a b jet from the direct top-quark decay or a jet from the hadronically decaying W boson), given its measured value from the CSV tagger discriminant (see section 4). Since the top-quark and W-boson reconstructed masses are dominated by experimental resolu- tion effects, the parameters σmt,t and σ M lep, had W in eq. (5.1) are approximated as Gaussian widths, which are determined from simulation. In the process, a constrained kinematic fit is performed, which leads to an improved determination of the unmeasured neutrino momentum component pνz , or in some cases a valid physical solution to be found when the analytical solution is imaginary due to detector mismeasurements, and to a more accurate reconstruction of the tt system. The momenta of the measured jets and lepton are allowed to vary within their resolutions and are required to comply with the same kinematic constraints used in eq. (5.1). The t and t candidates that are chosen for further analysis correspond to the particular configuration of lepton, neutrino, and jets that minimise χ2 comb. Any additional jets passing the event selection criteria that are not chosen as one of the four jets belonging to the tt system are discarded and no longer used in the analysis. Events for which the kinematic fit fails to find a real solution complying with the constraints are discarded. Following the full event selection in simulation studies, including the requirement of at least two identified b jets, this top-quark reconstruction algorithm associates the correct – 7 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 jet to the leptonic top-quark decay branch in 71% of all cases. If instead one loosens the requirement to at least one b jet in the event, the fraction of correctly assigned jets to the leptonic branch decreases to 63%. 6 Background estimation The main source of backgrounds in the analysis comes from top-quark pairs that decay into either fully leptonic or fully hadronic modes, or into semileptonic tt decays involving taus, and that pass the `+jets selections. Other backgrounds are, in order of decreasing importance: single top-quark events, events from processes involving W-boson production with jets (W+jets) and Drell-Yan production with jets (DY+jets). The normalisation of the tt processes is determined by the fitting procedure described in section 7, and the predictions for single top-quark processes are determined from simulation (see section 3). The cross section for inclusive W+jets and DY+jets production could be poorly predicted in the specific phase space region where those events become background for top-quark production, corresponding to high jet multiplicities and events containing heavy-flavour jets. For this reason, the normalisation of the background coming from W+jets and DY+jets is determined using an approach partially based on the data. Muons and electrons are treated separately, since important requirements on background rejection, such as lepton isolation, are different. 6.1 Normalisation and shape of DY production with jets The normalisation and shape of distributions from DY+jets production, are studied using a data control sample defined by applying the same selection criteria described in section 4, except that events are required to contain an additional lepton of the same flavour and opposite charge. In this case, the MT variable is computed with the event Emiss T and the highest-pT lepton. Using a reference cross section of 3048 pb [30] for the DY production decaying into pairs of muons, electrons, and taus with invariant mass above 50 GeV/c2, the normalisation of simulated samples for DY+jets is found to agree with the data within the statistical uncertainty of about 10%. The following systematic uncertainties are then considered in the estimation of the normalisation scale factor Ndata/Nsimulation, between the amount of data and simulated events: • the PDFs [31, 32] used to generate the samples: ±3.6%; • the uncertainties in the jet energy scale: ±12.7%; • the uncertainties in the jet energy resolution: ±2.6%; • the uncertainty in the integrated luminosity [11]: ±2.2%; • the difference observed on Ndata/Nsimulation for events inside and outside a window of 30±20 GeV/c2 around the Z-boson mass, where the DY+jets background is probed: ±20.0%. – 8 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 Including all systematic sources, this leads to a total estimated uncertainty of 30.0% for the DY+jets normalisation. The shape of the cos θ∗ distribution from DY+jets background events is verified using this same control sample. The top-quark reconstruction algorithm described in section 5 is applied, except that the highest-pT lepton is used in the kinematic fit. The shapes observed in data and simulation are found to be in agreement. 6.2 Normalisation and shape of W-boson production with jets Due to the charge of the valence quarks in the colliding protons, more positively (`+) than negatively (`−) charged leptons are produced in pp collisions for single top-quark production and W+jets events. On the other hand, the amount of `+ and `− produced in DY+jets and tt processes is, to a very good approximation, the same in pp collisions. Hence, by keeping the predicted cross section of single-top-quark production fixed, the background contribution from W+jets production can be determined from the charge asymmetry N+ − N−, where N+ (N−) is the number of `+ (`−). The total number of W+jets events in the data sample is thus predicted to be (N+ +N−)W+jets predicted = RW ±(MC) × (N+ −N−)data, (6.1) where RW ±(MC) = (N+ + N−)W+jets MC /(N+ − N−)W+jets MC is estimated using simulated events and (N+−N−)data is the measured asymmetry in data. The fixed contribution of single-top quark from simulation is subtracted from the total charge asymmetry in eq. (6.1) so that (N+ −N−)data = (N+ −N−)totaldata − (N+ −N−)single-topMC . The predicted normalisation is found to be consistent with the expectations from the simulation within relatively large statistical uncertainties. A total systematic uncertainty of 100% is assigned to the normalisation of the predicted W+jets background. To test the effect of possible biases due to the assumed background shape from simulation, the analysis is repeated by dividing the data sample in bins of cos θ∗. To increase the statistical power of the test, the jet selection criteria are partially relaxed. Within the precision of the test, the shape of the cos θ∗ distribution predicted from the simulation is found to be in agreement with the data, and any possible bias is negligible compared with the normalisation uncertainty considered. 7 Determination of helicity fractions A fitting procedure based on a MC simulation reweighting technique is used to simul- taneously account for experimental resolution effects and for the dependencies on the W-boson helicity fractions. Any new helicity configuration can be obtained from the original configuration used in the simulation via an algorithm that (re)weights each event according to the generated, matrix-element level cos θ∗ values for each W → `ν and W→ qq decay in the tt simulated sample. Because of the QCD production mechanism for tt events, top quarks can be considered as unpolarised on average, to a high degree of precision. Spin correlations between the two – 9 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 top quarks in the event do not modify this picture. This is because the average spin proper- ties of the top quark associated with the leptonic W-boson decay branch do not change after averaging over the phase space variables of the other top quark (from the hadronic W-boson decay branch) in the event. This scenario, which has been assumed in past and recent W- boson helicity studies [4, 33, 34], implies that the phase space density for the reconstructed cos θ∗ variable of each decay branch, ρ(cos θ∗rec), can be decoupled from the rest of the event. Therefore, at matrix-element (generator) level, we assume the following dependence: ρ(cos θ∗gen) ≡ 1 N dN d cos θ∗gen = 3 8 FL(1− cos θ∗gen)2 + 3 4 F0 sin2 θ∗gen + 3 8 FR(1 + cos θ∗gen)2, (7.1) where (FL, F0, FR) ≡ ~F are the helicity fractions to be measured, and θ∗gen is the matrix-element level quantity for the helicity angle. For normalisation reasons, the helicity fractions are constrained to satisfy FL + F0 + FR = 1. A new cos θ∗rec distribution for a particular configuration of helicity fractions FL, F0, FR is then obtained by reweighting each fully simulated event by the weight W (cos θ∗gen; ~F ) = Wlep(cos θ∗gen; ~F )×Whad(cos θ∗gen; ~F ), (7.2) where Wlep, had is defined for the leptonic and hadronic branches respectively as: Wlep, had(cos θ∗gen; ~F ) ≡ ρ(cos θ∗gen) ρSM(cos θ∗gen) = = 3 8 FL(1−cos θ∗gen)2+ 3 4 F0 sin2 θ∗gen+ 3 8 FR(1+cos θ∗gen)2 3 8 F SM L (1−cos θ∗gen)2+ 3 4 F SM 0 sin2 θ∗gen+ 3 8 F SM R (1+cos θ∗gen)2 . (7.3) In the expression above, F SM L , F SM 0 , F SM R are the helicity fractions that correspond to the expected SM values and which are used to generate the reference simulated sample. The new reweighted distribution automatically takes into account, by construction, all detector resolution and acceptance effects, as described by the simulation. The helicity fractions are measured by maximising a binned Poisson likelihood function, L(~F ) = ∏ bin i NMC(i; ~F )Ndata(i) (Ndata(i))! exp (−NMC(i; ~F )), (7.4) which is constructed using the numbers of observed Ndata(i) and expected NMC(i; ~F ) events in each cos θ∗rec bin i. The numbers of expected events are given by: NMC(i, ~F ) = NBKG(i) +Ntt(i; ~F ), (7.5) Ntt(i; ~F ) = Ftt  ∑ tt events W (cos θ∗gen(i); ~F )  , (7.6) NBKG(i) = NW+jets(i) +NDY+jets(i) +Nsingle-top(i) + Ftt ×Ntt non-`+jets(i). (7.7) – 10 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 The normalisation parameter for the tt component Ftt is not sensitive to the helicity fractions, but it does absorb a large fraction of the experimental and theoretical systematic uncertainties in the predicted rates. Uncertainties on the normalisation of backgrounds are considered as a separate source of systematics. Two types of fits are performed. In the fits denoted by 3D, the fractions F0, FL, and the normalisation factor Ftt (see eq. (7.6)) are treated as free parameters, and FR is determined by the constraint FR = 1−F0−FL. In the fits denoted by 2D, F0 and Ftt are taken as free parameters, setting FR = 0, and solving for FL from the constraint FL = 1 − F0. Within the expected experimental uncertainties, the FR = 0 condition is satisfied both in the SM and in anomalous coupling scenarios where only gR is different from zero. We perform the measurements by fitting either the cos θ∗ distribution from the leptonic branch or the |coshad θ∗| distribution from the hadronic branch. The two variables are fitted separately and not simultaneously, so as to avoid any possible biases due to top-quark spin correlations. At matrix-element level, |coshad θ∗| is only sensitive to the F0 fraction (or, alternatively, to the combination FL + FR), and is therefore only used in the context of the 2D fits. 8 Comparison between data and SM predictions The agreement between data and expectation from the SM is extensively investigated. First, the normalisations of simulated samples involving W+jets and DY+jets are determined and applied using control data samples discussed in section 6. Next, ref- erence cross sections 64.6 pb [35] for single-top-quark (t-channel), and 15.74 pb [36] for W-boson-associated single-top-quark processes, which correspond to approximated NNLO predictions, and NLO tt cross sections are used to scale the respective simulated samples to an integrated luminosity of 5.0 fb−1. Finally, the simulated samples involving top quarks are corrected for both the presence of pileup and the efficiency correction factors εDATA ` /εMC ` and εDATA b-tag /ε MC b-tag. Table 1 presents the number of data events observed in the muon+jets and elec- tron+jets channels which is compared to the predictions for the signal (tt) and background processes. In the table, columns two and four display the number of remaining events after applying the selection criteria described in section 4. Columns three and five display the subsample of these events where a tt pair candidate is found, reconstructed as described in the previous section. The yields observed in data are in agreement with the predicted yields expected from the SM. A comparison of shapes between data and simulation, after having applied the kine- matic fit, for the variables relevant to this analysis is presented in figures 1 to 3. Figure 1 shows the kinematic distributions for the leptons in the event: transverse momentum and pseudorapidity of muons and electrons, and the neutrino transverse momentum. Figure 2 displays the pT distributions for jets related to the top-quark reconstruction, including the b jets (from leptonic and hadronic tops) and the jets from the hadronic W-boson decay. The shapes of all important kinematic distributions in the data, which describe the tt sys- – 11 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 muon+jets electron+jets Process Selected KF Selected KF Single top 372 319 265 223 W+jets 351 282 234 203 DY+jets 43 38 43 34 tt non `+jets 928 828 629 547 Total bkg. 1694 1467 1171 1007 tt signal 7597 7321 5391 5179 Total expected 9291 8788 6562 6186 Data 9268 8772 6526 6135 Table 1. Number of events expected from signal and background processes, together with the observed number of events in data for both the muon+jets and electron+jets channels. The columns labelled as “Selected” represent the number of events passing the analysis selection criteria; the columns labelled as “KF” represent the fraction of those events containing a reconstructed top- quark pair via a kinematic fit that has converged. tem, are well reproduced by the simulation, including those that do not depend strongly on the W-boson helicity (i.e. global properties of top-quark and W-boson systems). The most relevant distributions for this analysis, the cosine distributions of the helicity angles computed according to the definitions discussed in section 1, are shown in figure 3 for both the muon+jets and electron+jets channels. The agreement between data and simulation for cos θ∗ and |coshad θ∗| is observed to be satisfactory. This suggests that, even before any attempt to measure the W-boson helicity fractions is made, the data prefer values close to the SM predictions. 9 Systematic uncertainties Several sources of systematic effects, which could possibly bias the measurement of the W-boson helicity fractions, have been investigated and their corresponding uncertainties in the measurement determined. The scale of the jet energy (JES) calibration is determined from data, which is then applied as a pT- and η-dependent correction to the simulation; the JES calibration has an uncertainty that typically varies between 2% and 4% [24]. To estimate the effect that the JES uncertainty has on the W-boson helicity measurement, the pT of all jets are system- atically shifted together, either up or down, by their corresponding pT- and η-dependent uncertainty. Because the missing transverse energy is corrected due to the presence of JES calibrated jets, the systematic shifts of the jet pTs in the event are propagated and lead to systematic shifts in the momentum imbalance. The full analysis, including the reconstruc- tion of top quarks and the resulting measurements of the W-boson helicity fractions, is then repeated. The jet energy resolution (JER) is observed to be different in data compared with the simulation. Jet energy resolutions are typically 5–13% larger in data than in – 12 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 40 60 80 100 120 140 160 180 200 E n tr ie s 100 200 300 400 500 600 700 muon+jetstt bkgtt W+jets & DY+jets Single Top Muon+jets -1 = 7 TeV, 5.0 fbsCMS [GeV/c] T Muon p 40 60 80 100 120 140 160 180 200M C /d at a 0.8 1 1.2 40 60 80 100 120 140 160 180 200 E n tr ie s 100 200 300 400 500 electron+jetstt bkgtt W+jets & DY+jets Single Top Electron+jets -1 = 7 TeV, 5.0 fbsCMS [GeV/c] T Electron p 40 60 80 100 120 140 160 180 200M C /d at a 0.8 1 1.2 -3 -2 -1 0 1 2 3 E n tr ie s 100 200 300 400 500 muon+jetstt bkgtt W+jets & DY+jets Single Top Muon+jets -1 = 7 TeV, 5.0 fbsCMS ηMuon -3 -2 -1 0 1 2 3M C /d at a 0.8 1 1.2 -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 E n tr ie s 50 100 150 200 250 300 electron+jetstt bkgtt W+jets & DY+jets Single Top Electron+jets -1 = 7 TeV, 5.0 fbsCMS ηElectron -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5M C /d at a 0.8 1 1.2 0 20 40 60 80 100 120 140 160 180 200 E n tr ie s 100 200 300 400 500 600 700 muon+jetstt bkgtt W+jets & DY+jets Single Top Muon+jets -1 = 7 TeV, 5.0 fbsCMS [GeV/c]TNeutrino P 0 20 40 60 80 100 120 140 160 180 200M C /d at a 0.8 1 1.2 0 20 40 60 80 100 120 140 160 180 200 E n tr ie s 100 200 300 400 500 electron+jetstt bkgtt W+jets & DY+jets Single Top Electron+jets -1 = 7 TeV, 5.0 fbsCMS [GeV/c]TNeutrino P 0 20 40 60 80 100 120 140 160 180 200M C /d at a 0.8 1 1.2 Figure 1. Distributions of data compared to SM predictions for signal and expected backgrounds: charged lepton pT (top), η (centre) and neutrino pT (bottom) for the muon+jets (left) and elec- tron+jets (right) channels. Data are displayed as solid points, simulated tt signal distributions as red histograms, and the contribution from other background processes as coloured histograms. Overflows are displayed in the last bin of each histogram. At the bottom, the ratio between pre- diction and data is displayed. Only statistical uncertainties are shown. simulation, with uncertainties smaller than 5% [24]. The systematic uncertainties related to JER are estimated by over-smearing the reconstructed jets in simulated events, so that their transverse momentum resolution is the same as measured in the data. The effect is propagated to the missing transverse energy, and the full analysis is repeated, similar to the procedure used to estimate the JES uncertainty. – 13 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 50 100 150 200 250 E n tr ie s 100 200 300 400 500 600 muon+jetstt bkgtt W+jets & DY+jets Single Top Muon+jets -1 = 7 TeV, 5.0 fbsCMS [GeV/c] T Leptonic b-jet p 50 100 150 200 250M C /d at a 0.8 1 1.2 50 100 150 200 250 E n tr ie s 100 200 300 400 500 electron+jetstt bkgtt W+jets & DY+jets Single Top Electron+jets -1 = 7 TeV, 5.0 fbsCMS [GeV/c] T Leptonic b-jet p 50 100 150 200 250M C /d at a 0.8 1 1.2 50 100 150 200 250 E n tr ie s 100 200 300 400 500 600 700 muon+jetstt bkgtt W+jets & DY+jets Single Top Muon+jets -1 = 7 TeV, 5.0 fbsCMS [GeV/c] T Hadronic b-jet p 50 100 150 200 250M C /d at a 0.8 1 1.2 50 100 150 200 250 E n tr ie s 100 200 300 400 500 electron+jetstt bkgtt W+jets & DY+jets Single Top Electron+jets -1 = 7 TeV, 5.0 fbsCMS [GeV/c] T Hadronic b-jet p 50 100 150 200 250M C /d at a 0.8 1 1.2 0 20 40 60 80 100 120 140 160 180 200 E n tr ie s 200 400 600 800 1000 1200 1400 1600 1800 2000 muon+jetstt bkgtt W+jets & DY+jets Single Top Muon+jets -1 = 7 TeV, 5.0 fbsCMS [GeV/c] T Hadronic jets p 0 20 40 60 80 100 120 140 160 180 200M C /d at a 0.8 1 1.2 0 20 40 60 80 100 120 140 160 180 200 E n tr ie s 200 400 600 800 1000 1200 electron+jetstt bkgtt W+jets & DY+jets Single Top Electron+jets -1 = 7 TeV, 5.0 fbsCMS [GeV/c] T Hadronic jets p 0 20 40 60 80 100 120 140 160 180 200M C /d at a 0.8 1 1.2 Figure 2. Distributions of transverse momenta in data compared to SM predictions for signal and expected backgrounds: the b jets identified in the leptonic (top row) and hadronic (second row) branches, and the jets from the hadronic W decay with two entries per event (bottom row) for the muon+jets (left) and electron+jets (right) channels. Data are displayed as solid points, simulated tt signal distributions as red histograms, and the contribution from other background processes as coloured histograms. Overflows are displayed in the last bin of each histogram. At the bottom, the ratio between prediction and data is displayed. Only statistical uncertainties are shown. Uncertainties in the lepton identification efficiency are investigated by varying the efficiency correction factor εDATA ` /εMC ` . In the case of muons, the efficiency correction factor depends on the η position of the muon in the detector. Since the measurement of the W-boson helicity is mainly affected by shape-dependent effects, uncertainties due to the muon efficiency correction factor are estimated by repeating the full analysis replacing – 14 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 E n tr ie s 50 100 150 200 250 300 muon+jetstt bkgtt W+jets & DY+jets Single Top Muon+jets -1 = 7 TeV, 5.0 fbsCMS *)θcos( -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1M C /d at a 0.8 1 1.2 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 E n tr ie s 20 40 60 80 100 120 140 160 180 200 220 240 electron+jetstt bkgtt W+jets & DY+jets Single Top Electron+jets -1 = 7 TeV, 5.0 fbsCMS *)θcos( -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1M C /d at a 0.8 1 1.2 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 E n tr ie s 50 100 150 200 250 300 muon+jetstt bkgtt W+jets & DY+jets Single Top Muon+jets -1 = 7 TeV, 5.0 fbsCMS *)|θ(had|cos 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1M C /d at a 0.8 1 1.2 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 E n tr ie s 20 40 60 80 100 120 140 160 180 200 220 electron+jetstt bkgtt W+jets & DY+jets Single Top Electron+jets -1 = 7 TeV, 5.0 fbsCMS *)|θ(had|cos 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1M C /d at a 0.8 1 1.2 Figure 3. Cosine of the helicity angles cos θ∗ (top) and |coshad θ∗| (bottom) for the muon+jets (left) and electron+jets (right) channels. Data are displayed as solid points, simulated tt signal distributions as red histograms, and the contribution from other background processes as coloured histograms. At the bottom, the ratio between prediction and data is displayed. Systematic uncer- tainties are shown as hatched histograms. the η-dependent factors by an uniform correction. In the case of electrons, only a very mild dependency on η is observed and the corresponding uncertainty, derived by assuming a flat η-dependence, has very little impact on the W-boson helicity measurement. Therefore, the efficiency correction factor for electron identification is shifted together with a shift in the scale factor for the jet component of the trigger according to the pT and η position of the electrons and jets. The combination of electron and jet scale factors that lead to the maximum possible η-dependent effect is then applied and the full analysis is repeated. The efficiencies for b tagging are measured using a control sample of multijet events [26], in data and simulation. The correction factors εDATA b-tag /ε MC b-tag, which are ap- plied to the simulated samples, are functions of the pT and η of the jets, the number of the required b tags, as well as the number of heavy-flavour and light-flavour jets in the event. The scale factors are relatively uniform in the typical pT range of the selected jets, resulting in only a very small dependence of the W-boson helicity results on the b-tagging efficiency. The scale factors are varied by their uncertainties and the resulting differences are taken as a systematic uncertainty in the measured fractions. – 15 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 The effect of pileup is estimated by varying the number of minimum bias collisions superimposed on each simulated signal event, according to an uncertainty of 5% which includes the uncertainties in the inelastic pp cross section, luminosity and other modelling uncertainties. To account for a bias on the W-boson helicity measurement due to uncertainties in the normalisation from simulated background samples involving single top quarks, rela- tive to the signal, the assumed reference cross sections are varied. While CMS has mea- sured the single-top-quark production cross section in the t-channel with a 9% precision (67.2 ± 6.1 pb [37]), a systematic variation of ±15% is applied to cover the case of single- top-quark events produced with additional jets from radiation, which comprise the main contribution to this background component. Likewise, while CMS measures a cross section of 16+5 −4 pb [38] for the associated tW production case, the reference cross section used for the normalisation, 15.74 pb, is shifted by ±40%. Finally, since the tt sample normalisation is a free parameter in the helicity fits, the uncertainty associated with its initially assumed reference cross section does not affect the measurement. The normalisation of the W+jets and DY+jets samples is estimated using the method described in section 6 with an uncertainty of 100% and 30%, respectively. In both the W+jets and DY+jets events cases, the shapes of all relevant distributions in the simula- tion agree with the data and any systematic effects from variations in normalisations are much larger than any variations arising from shape; hence, systematic uncertainties due to possible differences in shape are negligible. While the sample size for the reference simulated tt dataset is chosen to be five times that of the processed data, it is possible for the reweighting method to introduce a sys- tematic bias by degrading the statistical power of the simulated sample due to weights which can be larger than unity. Hence, special care must be taken to ensure that the statistical uncertainty of the MC prediction for each cos θ∗rec bin is substantially smaller than the corresponding statistical uncertainty of the data. The uncertainties are estimated by repeating the analysis using a subsample of events that correspond to a fraction 1/N of the entire sample. The procedure is repeated many times, and the uncertainties on the W-boson helicity fractions taken as σ/ √ N where σ is the spread observed on the fraction. Several values of N , between 2 and 10, are tested, and result in very similar uncertainties. Since the W-boson helicity fractions depend directly on the top-quark mass, uncertain- ties in the latter could bias the measurement. This systematic effect is studied and taken into account, via tt samples simulated using MadGraph for different mt hypotheses. A variation of ±1.4 GeV/c2 [29] about the assumed central value of 172.5 GeV/c2 is assumed. Uncertainties on the helicity measurement from the choice of renormalisation and fac- torisation scales for the simulated signal samples are estimated using dedicated tt simulated samples that vary the renormalisation and factorisation scales and the scale of the first emis- sion in the parton shower, in a consistent manner, by factors of 0.5 and 2 with respect to a central value of Q, with Q2 = m2 t c 2 + ( ∑ pjetT )2. The kinematic scale used to match jets to partons in the signal simulation is estimated using dedicated tt samples where that match- ing parameter is varied by factors of 0.75 and 1.5 with respect to its central value of 40 GeV. – 16 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 The systematic effects due to the PDFs used to simulate the signal and background samples are estimated using two different methods [32], according to a reweighting tech- nique. Firstly, events are reweighted using 100 members [39] of the NNPDF21 set [40], and the W-boson helicity fractions are remeasured for each of them. For a given helicity fraction, the RMS of the distribution of measurements from the different PDF members pro- vides an uncertainty estimate that corresponds to 68% confidence level (CL). Secondly, the difference between the central values for CTEQ6L1 [16] (used in the analysis simulations) and MSTW2008lo68cl [31] is estimated. The systematic uncertainties in the measurements for the W-boson helicity fractions, due to intrinsic PDF uncertainties, are then taken as the largest difference between two different estimates. The impact of all of the above systematic effects on the W-boson helicity fractions is detailed in table 2, for the measurements for the leptonic side of the events (cos θ∗). The table shows results for the 3D and 2D fits, obtained by fitting two of the fractions or set- ting FR = 0, respectively. Measurements using final states containing either a muon or an electron are presented in columns two through seven. The last three columns display sys- tematic uncertainties for the combined measurements of the muon+jets and electron+jets channels, using the measurements from columns two through seven as inputs. The mea- surements are combined taking into account both statistical and systematic uncertainties. Common sources of uncertainties between the different measurements are assumed to be fully correlated. The dominant sources of systematic uncertainties in the leptonic side are the W+jets background normalisation, the signal modelling (tt renormalisation and factorization scales, and top-quark mass), and the statistics of the simulated samples. Systematic uncertainties in the measurements for the hadronic branch, using the |coshad θ∗|, are presented in table 3. Since |coshad θ∗| has no sensitivity to a measurement of FL−FR, only the 2D fits are performed. Uncertainties on the individual measurements in the electron and muon channels are presented in the first two columns; in the last column, the combination of muons+jets and electrons+jets channels is shown. The hadronic branch is seen to have larger systematic uncertainties, compared with the leptonic branch, due, in part, to the dominant W+jets background and the importance of uncertainties from the JES and JER, as well as PDFs. 10 Results The W-boson helicity fractions are measured according to the fits described in section 7. The unitary condition F0+FL+FR = 1 is used to determine either (a) the right-handed frac- tion FR from measurements of the free parameters, F0 and FL, in the 3D fits or (b) the left- handed fraction FL from the measurement of the free parameter F0 in the 2D fits assuming FR = 0. Table 4 presents the fit measurement of each helicity parameter, one decay channel at a time, together with the statistical and systematic uncertainties. The statistical corre- lation factor ρstat0L between F0 and FL is presented in the last column; the measurements are seen to be highly correlated. All measurements, from either the muon or electron channels, using either the leptonic or hadronic branches, are observed to be compatible within uncer- tainties. The measurements using the leptonic branch cos θ∗ are more precise, as expected. – 17 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 µ+jets (cos θ∗) e+jets (cos θ∗) `+jets (cos θ∗) Systematic 3D fit 2D fit 3D fit 2D fit 3D fit 2D fit Uncertainties ± ∆F0 ± ∆FL ± ∆F0 ± ∆F0 ± ∆FL ± ∆F0 ± ∆F0 ± ∆FL ± ∆F0 JES 0.005 0.003 0.001 0.006 0.002 0.003 0.006 0.003 0.001 JER 0.009 0.005 0.001 0.014 0.009 0.003 0.011 0.007 0.001 Lepton eff. 0.001 0.001 0.001 0.009 0.012 0.015 0.001 0.002 0.002 b-tag eff. 0.001 0.001 < 10−3 < 10−3 < 10−3 0.001 0.001 < 10−3 < 10−3 Pileup 0.013 0.011 0.008 0.008 0.007 0.005 0.002 < 10−3 0.008 Single-t bkg. 0.004 < 10−3 0.003 0.004 < 10−3 0.004 0.004 0.001 0.003 W+jets bkg. 0.019 0.007 0.006 0.009 0.006 0.022 0.013 0.004 0.006 DY+jets bkg. 0.002 0.001 0.001 0.001 < 10−3 0.001 0.001 < 10−3 0.001 MC statistics 0.016 0.012 0.009 0.019 0.015 0.012 0.016 0.012 0.010 Top-quark mass 0.011 0.008 0.007 0.025 0.018 0.014 0.016 0.011 0.019 tt scales 0.013 0.009 0.007 0.015 0.018 0.030 0.009 0.009 0.011 tt match. scale 0.004 0.001 0.006 0.010 0.013 0.016 0.011 0.010 0.008 PDF 0.002 0.001 0.003 0.004 0.002 0.002 0.002 < 10−3 0.003 Table 2. Summary of the systematic uncertainties for the analysis using only the leptonic branch of the event, for the 3D fit, fitting F0, FL, and Ftt (columns 2–3 for muon+jets analysis, 5–6 for electron+jets analysis, and 8–9 for the combination of both decay modes); and the 2D fit, fitting F0 and Ftt only (column 4 for muon+jets analysis, 7 for electron+jets analysis, and 10 for the combination of both decay modes). The numbers given correspond to the absolute uncertainty with respect to the central analysis: ∆F = (F central − F check). Table 5 presents various combinations of the results presented in table 4. Firstly, the muon+jets and electron+jets channels are combined using the leptonic branch measure- ments from the 3D fits. The χ2 per degree of freedom for that combination is 0.109/2, corresponding to a χ2-probability of 94.7%. Secondly, the 2D fit measurements of the F0 helicity fraction from the leptonic (cos θ∗) and hadronic (coshad θ∗) branches are combined, separately for each decay channel. While the leptonic branch dominates with a weight of about 90%, the total uncertainty of the combination nevertheless decreases. Finally, the most precise measurement of F0 is obtained by subsequently combining the 2D fit mea- surements across the muon+jets and electron+jets channels, following the combination of the 2D fit measurements from the leptonic and hadronic branches. Summaries of all measurements and their various combinations are presented in fig- ures 4 and 5 for the 3D and 2D types of fits, respectively. All measurements are compatible with each other, and also compatible with the expectations from the SM [6]. 11 Limits on anomalous couplings The measured helicity fractions can be used to set limits on anomalous Wtb couplings. We assume the minimal parametrisation of the Wtb vertex suggested in refs. [7, 8, 41] and as described in the introduction. We consider two specific scenarios. First, we assume VL = 1, VR = gL = 0 and leave Re(gR) as a free parameter. This CP-conserving scenario – 18 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 µ+jets (|coshad θ∗|) e+jets (|coshad θ∗|) `+jets (|coshad θ∗|) Systematic 2D fit 2D fit 2D fit Uncertainties ± ∆F0 ± ∆F0 ± ∆F0 JES 0.010 0.008 0.002 JER 0.042 0.032 0.038 Lepton eff. 0.002 0.002 0.001 b-tag eff. 0.003 < 10−3 0.002 Pileup 0.018 0.006 0.015 Single-t bkg. 0.005 0.007 0.006 W+jets bkg. 0.060 0.050 0.040 DY+jets bkg. 0.002 0.005 0.002 MC statistics 0.023 0.028 0.025 Top-quark mass 0.008 0.041 0.014 tt scales 0.022 0.033 0.027 tt match. scale 0.002 0.035 0.013 PDF 0.013 0.014 0.014 Table 3. Systematic uncertainties for the 2D fits using the hadronic branch of the tt system, and for the muon channel, electron channel, as well as the combination of both decay channels. The numbers given correspond to the absolute uncertainty with respect to the central analysis: ∆F = (F central − F check). Leptonic branch: cos θ∗ Fit Channel F0 ± (stat.) ± (syst.) FL ± (stat.) ± (syst.) FR ± (stat.) ± (syst.) ρstat0L 3D µ+jets 0.674 ±0.039±0.035 0.314 ±0.028±0.022 0.012 ±0.016±0.020 −0.95 3D e+jets 0.688 ±0.045±0.042 0.310 ±0.033±0.037 0.002 ±0.017±0.023 −0.95 2D µ+jets 0.698 ±0.021±0.019 0.302 ±0.021±0.019 fixed at 0 −1 2D e+jets 0.691 ±0.025±0.047 0.309 ±0.025±0.047 fixed at 0 −1 Hadronic branch: |coshad θ∗| Fit Channel F0 ± (stat.) ± (syst.) FL ± (stat.) ± (syst.) FR ± (stat.) ± (syst.) ρ0L 2D µ+jets 0.651 ±0.060±0.084 0.349 ±0.060±0.084 fixed at 0 −1 2D e+jets 0.629 ±0.060±0.093 0.371 ±0.060±0.093 fixed at 0 −1 Table 4. Measurements of the W-boson helicity fractions from the cos θ∗ (leptonic branch) and |coshad θ∗| (hadronic branch) distributions. The columns show the fit type, the decay channel, and the measurement of each helicity parameter, together with the statistical and systematic uncertain- ties. For the 3D fits, the last column presents the statistical correlation between F0 and FL, while for the 2D fit, total anticorrelation (FL = 1− F0) is assumed. – 19 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 Fit Channel(s) Branch Fraction ± (stat.) ± (syst.) [total] ρtotal0L F0 0.682 ±0.030±0.033 [0.045] 3D `+jets l FL 0.310 ±0.022±0.022 [0.032] −0.95 FR 0.008 ±0.012±0.014 [0.018] 2D µ+jets l+h F0 0.694 ±0.020±0.025 [0.032] FL 0.306 ±0.020±0.025 [0.032] −1 2D e+jets l+h F0 0.674 ±0.025±0.028 [0.037] FL 0.326 ±0.025±0.028 [0.037] −1 2D `+jets l+h F0 0.685 ±0.017±0.021 [0.027] FL 0.315 ±0.017±0.021 [0.027] −1 Table 5. The combined helicity fractions and their uncertainties, including the type of fit performed, the channels (` = e, µ combination) and branches of the tt system (“l” for leptonic, cos θ∗, and “h” for hadronic, |coshad θ∗|) used in the combination, as well as the total correlation between F0 and FL. 0 200 400 600 800 10000 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 -1=7 TeV, 5.0 fbsCMS 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0 0.5 1 1.5 2 2.5 3 +jets; leptonic branchµ )+jets; leptonic branchµcombined (e, e+jets; leptonic branch 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0 0.5 1 1.5 2 2.5 3 0F 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0 0.5 1 1.5 2 2.5 3 LF -0.15 -0.1 -0.05 0 0.05 0.1 0.15 0 0.5 1 1.5 2 2.5 3 RF Figure 4. Summary of the W-boson helicity measurements in semileptonic decays of top-quark pairs with 2011 data for 3D fits. The inner error bars represent the statistical uncertainties and the outer error bars the statistical and systematic uncertainties, added in quadrature. NNLO predictions from ref. [6] with their theoretical uncertainties are represented as hatched bands. is particularly interesting because indirect constraints to gR from radiative B-meson decay measurements are currently poor, Re(gR) ∈ [−0.15,+0.57] [42]. A specific feature of this scenario is that it does not provide any contribution to the right-handed helicity of the W boson, FR. The grand combination of the longitudinal helicity fraction F0 measurements, across both the leptonic and hadronic branches including both the muon and electron channels, and assuming FR = 0, is reinterpreted in terms of Re(gR), yielding Re(gR) = −0.008± 0.024 (stat.)+0.029 −0.030 (syst.), – 20 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 0 100 200 300 400 500 600 700 800 9000 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 -1=7 TeV, 5.0 fbsCMS 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0 1 2 3 4 5 6 7 8 9 +jets; combined leptonic+hadronic branchesµ +jets; hadronic branchµ +jets; leptonic branchµ )+jets and leptonic+hadronic branchesµcombined (e, )+jets; hadronic branchµcombined (e, )+jets; leptonic branchµcombined (e, e+jets; combined leptonic+hadronic branches e+jets; hadronic branch e+jets; leptonic branch 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0 1 2 3 4 5 6 7 8 9 0F 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0 1 2 3 4 5 6 7 8 9 LF Figure 5. Summary of the W-boson helicity measurements in semileptonic decays of top-quark pairs with 2011 data, for 2D fits assuming FR = 0. The inner error bars represent the statistical uncertainties and the outer error bars the statistical and systematic uncertainties, added in quadrature. NNLO predictions from ref. [6] with their theoretical uncertainties are represented as hatched bands. which is consistent with the SM expectations within the quoted uncertainties. In quoting this result we have omitted another minimum of the fit closer to the Re(gR) ≈ 0.8 region, which would lead to an increase of almost a factor of three in the single-top-quark cross section [43], which would not be consistent with the recent CMS measurement [37]. In terms of the effective dimension-six Lagrangian O33 uW defined in refs. [7, 8] we obtain the equivalent result, Re(C33 uW )/Λ2 = −0.088± 0.280 (stat.)+0.339 −0.352 (syst.) TeV−2, where Λ is the scale of new physics and Re(C33 uW ) the effective operator coefficient. In the second scenario, again assuming CP is conserved, we choose Re(gL) and Re(gR) as free parameters of the fit. Limits on those parameters are determined using the combined measurements of the muon+jets and electron+jets channels from the 3D fit of the leptonic branch, cos θ∗. The results of the likelihood fit for the parameters F0 and FL can be reinterpreted in terms of the parameters Re(gL) and Re(gR). Figure 6 shows the regions of the Re(gL), Re(gR) plane allowed at 68% and 95% CL. As in the first scenario, a region near Re(gL) = 0 and Re(gR) � 0, allowed by the fit but excluded by the CMS single-top quark measurement, is not shown. The result obtained from the first scenario represents an improvement of about 50% on the precision of Re(gR) with respect to previous measurements [4], while the limits from the second scenario are similar to those from ref. [4]. – 21 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 ) L Re(g -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 ) R R e( g -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 = 0 R = 1, VLV Muon and electron channels, combined -1 = 7 TeV, 5.0 fbsCMS 68% CL 95% CL SM Figure 6. Limits on the real components of the anomalous couplings gL, gR at 68% and 95% CL, for VL = 1 and VR = 0. The SM prediction (gR = 0 and gL = 0) is also shown. 12 Summary The W-boson helicity has been measured in top-quark-pair events decaying semileptoni- cally, both in the muon+jets and in the electron+jets channels, using proton-proton col- lisions data at √ s = 7 TeV, corresponding to an integrated luminosity of 5.0 fb−1. Both the leptonic and the hadronic branches of the decay have been studied. The most precise measurement, not constraining FR to the SM, corresponding to the combination of muon and electron channels, using only the leptonic branch, yields: F0 = 0.682± 0.030 (stat.)± 0.033 (syst.), FL = 0.310± 0.022 (stat.)± 0.022 (syst.), FR = 0.008± 0.012 (stat.)± 0.014 (syst.), with a correlation coefficient of −0.95 between F0 and FL. The measured W-boson helicity fractions are in agreement with the predictions from the standard model. Assuming a minimal parametrisation of the Wtb vertex, stringent limits on the real components of the anomalous couplings gL and gR are also derived. Acknowledgments We congratulate our colleagues in the CERN accelerator departments for the excellent per- formance of the LHC and thank the technical and administrative staffs at CERN and at – 22 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centres and personnel of the Worldwide LHC Computing Grid for delivering so effectively the computing infrastructure essential to our analyses. Finally, we acknowledge the enduring support for the construction and operation of the LHC and the CMS detector provided by the following funding agencies: the Austrian Federal Ministry of Science and Research and the Austrian Science Fund; the Belgian Fonds de la Recherche Scientifique, and Fonds voor Wetenschappelijk Onderzoek; the Brazilian Funding Agencies (CNPq, CAPES, FAPERJ, and FAPESP); the Bulgarian Ministry of Education and Science; CERN; the Chinese Academy of Sciences, Ministry of Science and Technology, and National Natural Science Foundation of China; the Colombian Funding Agency (COLCIENCIAS); the Croatian Ministry of Science, Education and Sport; the Re- search Promotion Foundation, Cyprus; the Ministry of Education and Research, Recurrent financing contract SF0690030s09 and European Regional Development Fund, Estonia; the Academy of Finland, Finnish Ministry of Education and Culture, and Helsinki Institute of Physics; the Institut National de Physique Nucléaire et de Physique des Particules / CNRS, and Commissariat à l’Énergie Atomique et aux Énergies Alternatives / CEA, France; the Bundesministerium für Bildung und Forschung, Deutsche Forschungsgemeinschaft, and Helmholtz-Gemeinschaft Deutscher Forschungszentren, Germany; the General Secretariat for Research and Technology, Greece; the National Scientific Research Foundation, and National Office for Research and Technology, Hungary; the Department of Atomic En- ergy and the Department of Science and Technology, India; the Institute for Studies in Theoretical Physics and Mathematics, Iran; the Science Foundation, Ireland; the Istituto Nazionale di Fisica Nucleare, Italy; the Korean Ministry of Education, Science and Tech- nology and the World Class University program of NRF, Republic of Korea; the Lithuanian Academy of Sciences; the Mexican Funding Agencies (CINVESTAV, CONACYT, SEP, and UASLP-FAI); the Ministry of Science and Innovation, New Zealand; the Pakistan Atomic Energy Commission; the Ministry of Science and Higher Education and the National Sci- ence Centre, Poland; the Fundação para a Ciência e a Tecnologia, Portugal; JINR, Dubna; the Ministry of Education and Science of the Russian Federation, the Federal Agency of Atomic Energy of the Russian Federation, Russian Academy of Sciences, and the Russian Foundation for Basic Research; the Ministry of Education, Science and Technological De- velopment of Serbia; the Secretaŕıa de Estado de Investigación, Desarrollo e Innovación and Programa Consolider-Ingenio 2010, Spain; the Swiss Funding Agencies (ETH Board, ETH Zurich, PSI, SNF, UniZH, Canton Zurich, and SER); the National Science Council, Taipei; the Thailand Center of Excellence in Physics, the Institute for the Promotion of Teaching Science and Technology of Thailand, Special Task Force for Activating Research and the National Science and Technology Development Agency of Thailand; the Scientific and Technical Research Council of Turkey, and Turkish Atomic Energy Authority; the Science and Technology Facilities Council, U.K.; the US Department of Energy, and the US National Science Foundation. Individuals have received support from the Marie-Curie programme and the European Research Council and EPLANET (European Union); the Leventis Foundation; the A. P. Sloan Foundation; the Alexander von Humboldt Founda- tion; the Belgian Federal Science Policy Office; the Fonds pour la Formation à la Recherche – 23 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 dans l’Industrie et dans l’Agriculture (FRIA-Belgium); the Agentschap voor Innovatie door Wetenschap en Technologie (IWT-Belgium); the Ministry of Education, Youth and Sports (MEYS) of Czech Republic; the Council of Science and Industrial Research, India; the Compagnia di San Paolo (Torino); the HOMING PLUS programme of Foundation for Pol- ish Science, cofinanced by EU, Regional Development Fund; and the Thalis and Aristeia programmes cofinanced by EU-ESF and the Greek NSRF. Open Access. This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited. References [1] CDF collaboration, F. Abe et al., Observation of top quark production in p̄p collisions, Phys. Rev. Lett. 74 (1995) 2626 [hep-ex/9503002] [INSPIRE]. [2] D0 collaboration, S. Abachi et al., Observation of the top quark, Phys. Rev. Lett. 74 (1995) 2632 [hep-ex/9503003] [INSPIRE]. [3] CDF and D0 collaborations, T. Aaltonen et al., Combination of CDF and D0 measurements of the W boson helicity in top quark decays, Phys. Rev. D 85 (2012) 071106 [arXiv:1202.5272] [INSPIRE]. [4] ATLAS collaboration, Measurement of the W boson polarization in top quark decays with the ATLAS detector, JHEP 06 (2012) 088 [arXiv:1205.2484] [INSPIRE]. [5] M. Fischer, S. Groote, J. Korner and M. Mauser, Longitudinal, transverse plus and transverse minus W bosons in unpolarized top quark decays at O(αs), Phys. Rev. D 63 (2001) 031501 [hep-ph/0011075] [INSPIRE]. [6] A. Czarnecki, J.G. Korner and J.H. Piclum, Helicity fractions of W bosons from top quark decays at NNLO in QCD, Phys. Rev. D 81 (2010) 111503 [arXiv:1005.2625] [INSPIRE]. [7] J. Aguilar-Saavedra, A minimal set of top anomalous couplings, Nucl. Phys. B 812 (2009) 181 [arXiv:0811.3842] [INSPIRE]. [8] J. Aguilar-Saavedra, A minimal set of top-Higgs anomalous couplings, Nucl. Phys. B 821 (2009) 215 [arXiv:0904.2387] [INSPIRE]. [9] J. Aguilar-Saavedra and J. Bernabeu, W polarisation beyond helicity fractions in top quark decays, Nucl. Phys. B 840 (2010) 349 [arXiv:1005.5382] [INSPIRE]. [10] CMS collaboration, The CMS experiment at the CERN LHC, 2008 JINST 3 S08004 [INSPIRE]. [11] CMS collaboration, Absolute calibration of the luminosity measurement at CMS: winter 2012 update, CMS-PAS-SMP-12-008, CERN, Geneva Switzerland (2012). [12] J. Alwall, M. Herquet, F. Maltoni, O. Mattelaer and T. Stelzer, MadGraph 5: going beyond, JHEP 06 (2011) 128 [arXiv:1106.0522] [INSPIRE]. [13] T. Sjöstrand, S. Mrenna and P.Z. Skands, PYTHIA 6.4 physics and manual, JHEP 05 (2006) 026 [hep-ph/0603175] [INSPIRE]. – 24 – http://dx.doi.org/10.1103/PhysRevLett.74.2626 http://dx.doi.org/10.1103/PhysRevLett.74.2626 http://arxiv.org/abs/hep-ex/9503002 http://inspirehep.net/search?p=find+EPRINT+hep-ex/9503002 http://dx.doi.org/10.1103/PhysRevLett.74.2632 http://dx.doi.org/10.1103/PhysRevLett.74.2632 http://arxiv.org/abs/hep-ex/9503003 http://inspirehep.net/search?p=find+EPRINT+hep-ex/9503003 http://dx.doi.org/10.1103/PhysRevD.85.071106 http://arxiv.org/abs/1202.5272 http://inspirehep.net/search?p=find+EPRINT+arXiv:1202.5272 http://dx.doi.org/10.1007/JHEP06(2012)088 http://arxiv.org/abs/1205.2484 http://inspirehep.net/search?p=find+EPRINT+arXiv:1205.2484 http://dx.doi.org/10.1103/PhysRevD.63.031501 http://dx.doi.org/10.1103/PhysRevD.63.031501 http://arxiv.org/abs/hep-ph/0011075 http://inspirehep.net/search?p=find+EPRINT+hep-ph/0011075 http://dx.doi.org/10.1103/PhysRevD.81.111503 http://arxiv.org/abs/1005.2625 http://inspirehep.net/search?p=find+EPRINT+arXiv:1005.2625 http://dx.doi.org/10.1016/j.nuclphysb.2008.12.012 http://dx.doi.org/10.1016/j.nuclphysb.2008.12.012 http://arxiv.org/abs/0811.3842 http://inspirehep.net/search?p=find+EPRINT+arXiv:0811.3842 http://dx.doi.org/10.1016/j.nuclphysb.2009.06.022 http://dx.doi.org/10.1016/j.nuclphysb.2009.06.022 http://arxiv.org/abs/0904.2387 http://inspirehep.net/search?p=find+EPRINT+arXiv:0904.2387 http://dx.doi.org/10.1016/j.nuclphysb.2010.07.012 http://arxiv.org/abs/1005.5382 http://inspirehep.net/search?p=find+EPRINT+arXiv:1005.5382 http://dx.doi.org/10.1088/1748-0221/3/08/S08004 http://inspirehep.net/search?p=find+J+JINST,3,S08004 http://cds.cern.ch/record/1434360 http://dx.doi.org/10.1007/JHEP06(2011)128 http://arxiv.org/abs/1106.0522 http://inspirehep.net/search?p=find+EPRINT+arXiv:1106.0522 http://dx.doi.org/10.1088/1126-6708/2006/05/026 http://dx.doi.org/10.1088/1126-6708/2006/05/026 http://arxiv.org/abs/hep-ph/0603175 http://inspirehep.net/search?p=find+EPRINT+hep-ph/0603175 J H E P 1 0 ( 2 0 1 3 ) 1 6 7 [14] N. Davidson, G. Nanava, T. Przedzinski, E. Richter-Was and Z. Was, Universal interface of TAUOLA technical and physics documentation, Comput. Phys. Commun. 183 (2012) 821 [arXiv:1002.0543] [INSPIRE]. [15] S. Alioli, P. Nason, C. Oleari and E. Re, NLO vector-boson production matched with shower in POWHEG, JHEP 07 (2008) 060 [arXiv:0805.4802] [INSPIRE]. [16] J. Pumplin et al., New generation of parton distributions with uncertainties from global QCD analysis, JHEP 07 (2002) 012 [hep-ph/0201195] [INSPIRE]. [17] H.-L. Lai et al., New parton distributions for collider physics, Phys. Rev. D 82 (2010) 074024 [arXiv:1007.2241] [INSPIRE]. [18] CMS collaboration, Particle-flow event reconstruction in CMS and performance for jets, taus and MET, CMS-PAS-PFT-09-001, CERN, Geneva Switzerland (2009). [19] CMS collaboration, Jet performance in pp collisions at 7 TeV, CMS-PAS-JME-10-003, CERN, Geneva Switzerland (2010). [20] CMS collaboration, Measurement of the top-antitop production cross section in pp collisions at √ s = 7 TeV using the kinematic properties of events with leptons and jets, Eur. Phys. J. C 71 (2011) 1721 [arXiv:1106.0902] [INSPIRE]. [21] CMS collaboration, Electron reconstruction and identification at √ s = 7 TeV, CMS-PAS-EGM-10-004, CERN, Geneva Switzerland (2010). [22] M. Cacciari, G.P. Salam and G. Soyez, The anti-kT jet clustering algorithm, JHEP 04 (2008) 063 [arXiv:0802.1189] [INSPIRE]. [23] M. Cacciari, G.P. Salam and G. Soyez, FastJet user manual, Eur. Phys. J. C 72 (2012) 1896 [arXiv:1111.6097] [INSPIRE]. [24] CMS collaboration, Determination of jet energy calibration and transverse momentum resolution in CMS, 2011 JINST 6 P11002 [arXiv:1107.4277] [INSPIRE]. [25] CMS collaboration, Missing transverse energy performance of the CMS detector, 2011 JINST 6 P09001 [arXiv:1106.5048] [INSPIRE]. [26] CMS collaboration, Identification of b-quark jets with the CMS experiment, 2013 JINST 8 P04013 [arXiv:1211.4462] [INSPIRE]. [27] CMS collaboration, Measurements of inclusive W and Z cross sections in pp collisions at√ s = 7 TeV, JHEP 01 (2011) 080 [arXiv:1012.2466] [INSPIRE]. [28] Particle Data Group collaboration, J. Beringer et al., Review of particle physics (RPP), Phys. Rev. D 86 (2012) 010001 [INSPIRE]. [29] ATLAS, D0, CDF and CMS collaborations, F. Déliot, Combination of the top-quark mass measurements from the Tevatron and from the LHC colliders, PoS(ICHEP2012)218 [arXiv:1302.0830] [INSPIRE]. [30] R. Gavin, Y. Li, F. Petriello and S. Quackenbush, FEWZ 2.0: a code for hadronic Z production at next-to-next-to-leading order, Comput. Phys. Commun. 182 (2011) 2388 [arXiv:1011.3540] [INSPIRE]. [31] A. Martin, W. Stirling, R. Thorne and G. Watt, Parton distributions for the LHC, Eur. Phys. J. C 63 (2009) 189 [arXiv:0901.0002] [INSPIRE]. [32] M. Botje et al., The PDF4LHC working group interim recommendations, arXiv:1101.0538 [INSPIRE]. – 25 – http://dx.doi.org/10.1016/j.cpc.2011.12.009 http://arxiv.org/abs/1002.0543 http://inspirehep.net/search?p=find+EPRINT+arXiv:1002.0543 http://dx.doi.org/10.1088/1126-6708/2008/07/060 http://arxiv.org/abs/0805.4802 http://inspirehep.net/search?p=find+EPRINT+arXiv:0805.4802 http://dx.doi.org/10.1088/1126-6708/2002/07/012 http://arxiv.org/abs/hep-ph/0201195 http://inspirehep.net/search?p=find+J+JHEP,0207,012 http://dx.doi.org/10.1103/PhysRevD.82.074024 http://arxiv.org/abs/1007.2241 http://inspirehep.net/search?p=find+EPRINT+arXiv:1007.2241 http://cds.cern.ch/record/1194487 http://cds.cern.ch/record/1279362 http://dx.doi.org/10.1140/epjc/s10052-011-1721-3 http://dx.doi.org/10.1140/epjc/s10052-011-1721-3 http://arxiv.org/abs/1106.0902 http://inspirehep.net/search?p=find+J+Eur.Phys.J.,C71,1721 http://cds.cern.ch/record/1299116 http://dx.doi.org/10.1088/1126-6708/2008/04/063 http://dx.doi.org/10.1088/1126-6708/2008/04/063 http://arxiv.org/abs/0802.1189 http://inspirehep.net/search?p=find+EPRINT+arXiv:0802.1189 http://dx.doi.org/10.1140/epjc/s10052-012-1896-2 http://arxiv.org/abs/1111.6097 http://inspirehep.net/search?p=find+EPRINT+arXiv:1111.6097 http://dx.doi.org/10.1088/1748-0221/6/11/P11002 http://arxiv.org/abs/1107.4277 http://inspirehep.net/search?p=find+J+JINST,6,P11002 http://dx.doi.org/10.1088/1748-0221/6/09/P09001 http://dx.doi.org/10.1088/1748-0221/6/09/P09001 http://arxiv.org/abs/1106.5048 http://inspirehep.net/search?p=find+EPRINT+arXiv:1106.5048 http://dx.doi.org/10.1088/1748-0221/8/04/P04013 http://dx.doi.org/10.1088/1748-0221/8/04/P04013 http://arxiv.org/abs/1211.4462 http://inspirehep.net/search?p=find+EPRINT+arXiv:1211.4462 http://dx.doi.org/10.1007/JHEP01(2011)080 http://arxiv.org/abs/1012.2466 http://inspirehep.net/search?p=find+J+JHEP,1101,080 http://dx.doi.org/10.1103/PhysRevD.86.010001 http://inspirehep.net/search?p=find+J+Phys.Rev.,D86,010001 http://pos.sissa.it/cgi-bin/reader/contribution.cgi?id=PoS(ICHEP2012)218 http://arxiv.org/abs/1302.0830 http://inspirehep.net/search?p=find+EPRINT+arXiv:1302.0830 http://dx.doi.org/10.1016/j.cpc.2011.06.008 http://arxiv.org/abs/1011.3540 http://inspirehep.net/search?p=find+EPRINT+arXiv:1011.3540 http://dx.doi.org/10.1140/epjc/s10052-009-1072-5 http://dx.doi.org/10.1140/epjc/s10052-009-1072-5 http://arxiv.org/abs/0901.0002 http://inspirehep.net/search?p=find+EPRINT+arXiv:0901.0002 http://arxiv.org/abs/1101.0538 http://inspirehep.net/search?p=find+EPRINT+arXiv:1101.0538 J H E P 1 0 ( 2 0 1 3 ) 1 6 7 [33] CDF collaboration, T. Aaltonen et al., Measurement of W -boson helicity fractions in top-quark decays using cos(θ∗), Phys. Lett. B 674 (2009) 160 [arXiv:0811.0344] [INSPIRE]. [34] D0 collaboration, V.M. Abazov et al., Measurement of the W boson helicity in top quark decays using 5.4 fb−1 of pp̄ collision data, Phys. Rev. D 83 (2011) 032009 [arXiv:1011.6549] [INSPIRE]. [35] 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]. [36] 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]. [37] CMS collaboration, Measurement of the single-top-quark t-channel cross section in pp collisions at √ s = 7 TeV, JHEP 12 (2012) 035 [arXiv:1209.4533] [INSPIRE]. [38] CMS collaboration, Evidence for associated production of a single top quark and W boson in pp collisions at √ s = 7 TeV, Phys. Rev. Lett. 110 (2013) 022003 [arXiv:1209.3489] [INSPIRE]. [39] M. Whalley, D. Bourilkov and R. Group, The Les Houches accord PDFs (LHAPDF) and LHAGLUE, hep-ph/0508110 [INSPIRE]. [40] R.D. Ball et al., Impact of heavy quark masses on parton distributions and LHC phenomenology, Nucl. Phys. B 849 (2011) 296 [arXiv:1101.1300] [INSPIRE]. [41] J. Aguilar-Saavedra, J. Carvalho, N.F. Castro, F. Veloso and A. Onofre, Probing anomalous Wtb couplings in top pair decays, Eur. Phys. J. C 50 (2007) 519 [hep-ph/0605190] [INSPIRE]. [42] B. Grzadkowski and M. Misiak, Anomalous Wtb coupling effects in the weak radiative B-meson decay, Phys. Rev. D 78 (2008) 077501 [Erratum ibid. D 84 (2011) 059903] [arXiv:0802.1413] [INSPIRE]. [43] J. Aguilar-Saavedra, Single top quark production at LHC with anomalous Wtb couplings, Nucl. Phys. B 804 (2008) 160 [arXiv:0803.3810] [INSPIRE]. – 26 – http://dx.doi.org/10.1016/j.physletb.2009.02.040 http://arxiv.org/abs/0811.0344 http://inspirehep.net/search?p=find+EPRINT+arXiv:0811.0344 http://dx.doi.org/10.1103/PhysRevD.83.032009 http://arxiv.org/abs/1011.6549 http://inspirehep.net/search?p=find+EPRINT+arXiv:1011.6549 http://dx.doi.org/10.1103/PhysRevD.83.091503 http://arxiv.org/abs/1103.2792 http://inspirehep.net/search?p=find+EPRINT+arXiv:1103.2792 http://dx.doi.org/10.1103/PhysRevD.82.054018 http://arxiv.org/abs/1005.4451 http://inspirehep.net/search?p=find+EPRINT+arXiv:1005.4451 http://dx.doi.org/10.1007/JHEP12(2012)035 http://arxiv.org/abs/1209.4533 http://inspirehep.net/search?p=find+EPRINT+arXiv:1209.4533 http://dx.doi.org/10.1103/PhysRevLett.110.022003 http://arxiv.org/abs/1209.3489 http://inspirehep.net/search?p=find+EPRINT+arXiv:1209.3489 http://arxiv.org/abs/hep-ph/0508110 http://inspirehep.net/search?p=find+EPRINT+hep-ph/0508110 http://dx.doi.org/10.1016/j.nuclphysb.2011.03.021 http://arxiv.org/abs/1101.1300 http://inspirehep.net/search?p=find+EPRINT+arXiv:1101.1300 http://dx.doi.org/10.1140/epjc/s10052-007-0289-4 http://arxiv.org/abs/hep-ph/0605190 http://inspirehep.net/search?p=find+EPRINT+hep-ph/0605190 http://dx.doi.org/10.1103/PhysRevD.84.059903 http://dx.doi.org/10.1103/PhysRevD.84.059903 http://arxiv.org/abs/0802.1413 http://inspirehep.net/search?p=find+EPRINT+arXiv:0802.1413 http://dx.doi.org/10.1016/j.nuclphysb.2008.06.013 http://arxiv.org/abs/0803.3810 http://inspirehep.net/search?p=find+EPRINT+arXiv:0803.3810 J H E P 1 0 ( 2 0 1 3 ) 1 6 7 The CMS collaboration Yerevan Physics Institute, Yerevan, Armenia S. Chatrchyan, V. Khachatryan, A.M. Sirunyan, A. Tumasyan Institut für Hochenergiephysik der OeAW, Wien, Austria W. Adam, T. Bergauer, M. Dragicevic, J. Erö, C. Fabjan1, M. Friedl, R. Frühwirth1, V.M. Ghete, N. Hörmann, J. Hrubec, M. Jeitler1, W. Kiesenhofer, V. Knünz, M. Krammer1, I. Krätschmer, D. Liko, I. Mikulec, D. Rabady2, B. Rahbaran, C. Rohringer, H. Rohringer, R. Schöfbeck, J. Strauss, A. Taurok, W. Treberer-Treberspurg, W. Wal- tenberger, C.-E. Wulz1 National Centre for Particle and High Energy Physics, Minsk, Belarus V. Mossolov, N. Shumeiko, J. Suarez Gonzalez Universiteit Antwerpen, Antwerpen, Belgium S. Alderweireldt, M. Bansal, S. Bansal, T. Cornelis, E.A. De Wolf, X. Janssen, A. Knutsson, S. Luyckx, L. Mucibello, S. Ochesanu, B. Roland, R. Rougny, Z. Staykova, H. Van Haevermaet, P. Van Mechelen, N. Van Remortel, A. Van Spilbeeck Vrije Universiteit Brussel, Brussel, Belgium F. Blekman, S. Blyweert, J. D’Hondt, A. Kalogeropoulos, J. Keaveney, M. Maes, A. Ol- brechts, S. Tavernier, W. Van Doninck, P. Van Mulders, G.P. Van Onsem, I. Villella Université Libre de Bruxelles, Bruxelles, Belgium C. Caillol, B. Clerbaux, G. De Lentdecker, L. Favart, A.P.R. Gay, T. Hreus, A. Léonard, P.E. Marage, A. Mohammadi, L. Perniè, T. Reis, T. Seva, L. Thomas, C. Vander Velde, P. Vanlaer, J. Wang Ghent University, Ghent, Belgium V. Adler, K. Beernaert, L. Benucci, A. Cimmino, S. Costantini, S. Dildick, G. Garcia, B. Klein, J. Lellouch, A. Marinov, J. Mccartin, A.A. Ocampo Rios, D. Ryckbosch, M. Sigamani, N. Strobbe, F. Thyssen, M. Tytgat, S. Walsh, E. Yazgan, N. Zaganidis Université Catholique de Louvain, Louvain-la-Neuve, Belgium S. Basegmez, C. Beluffi3, G. Bruno, R. Castello, A. Caudron, L. Ceard, G.G. Da Silveira, C. Delaere, T. du Pree, D. Favart, L. Forthomme, A. Giammanco4, J. Hollar, P. Jez, V. Lemaitre, J. Liao, O. Militaru, C. Nuttens, D. Pagano, A. Pin, K. Piotrzkowski, A. Popov5, M. Selvaggi, J.M. Vizan Garcia Université de Mons, Mons, Belgium N. Beliy, T. Caebergs, E. Daubie, G.H. Hammad Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro, Brazil G.A. Alves, M. Correa Martins Junior, T. Martins, M.E. Pol, M.H.G. Souza Universidade do Estado do Rio de Janeiro, Rio de Janeiro, Brazil W.L. Aldá Júnior, W. Carvalho, J. Chinellato6, A. Custódio, E.M. Da Costa, D. De Jesus Damiao, C. De Oliveira Martins, S. Fonseca De Souza, H. Malbouisson, M. Malek, D. Matos – 27 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 Figueiredo, L. Mundim, H. Nogima, W.L. Prado Da Silva, A. Santoro, A. Sznajder, E.J. Tonelli Manganote6, A. Vilela Pereira Universidade Estadual Paulista a, Universidade Federal do ABC b, São Paulo, Brazil C.A. Bernardesb, F.A. Diasa,7, T.R. Fernandez Perez Tomeia, E.M. Gregoresb, C. Laganaa, P.G. Mercadanteb, S.F. Novaesa, Sandra S. Padulaa Institute for Nuclear Research and Nuclear Energy, Sofia, Bulgaria V. Genchev2, P. Iaydjiev2, S. Piperov, M. Rodozov, G. Sultanov, M. Vutova University of Sofia, Sofia, Bulgaria A. Dimitrov, R. Hadjiiska, V. Kozhuharov, L. Litov, B. Pavlov, P. Petkov Institute of High Energy Physics, Beijing, China J.G. Bian, G.M. Chen, H.S. Chen, C.H. Jiang, D. Liang, S. Liang, X. Meng, J. Tao, X. Wang, Z. Wang, H. Xiao State Key Laboratory of Nuclear Physics and Technology, Peking University, Beijing, China C. Asawatangtrakuldee, Y. Ban, Y. Guo, Q. Li, W. Li, S. Liu, Y. Mao, S.J. Qian, D. Wang, L. Zhang, W. Zou Universidad de Los Andes, Bogota, Colombia C. Avila, C.A. Carrillo Montoya, L.F. Chaparro Sierra, J.P. Gomez, B. Gomez Moreno, J.C. Sanabria Technical University of Split, Split, Croatia N. Godinovic, D. Lelas, R. Plestina8, D. Polic, I. Puljak University of Split, Split, Croatia Z. Antunovic, M. Kovac Institute Rudjer Boskovic, Zagreb, Croatia V. Brigljevic, K. Kadija, J. Luetic, D. Mekterovic, S. Morovic, L. Tikvica University of Cyprus, Nicosia, Cyprus A. Attikis, G. Mavromanolakis, J. Mousa, C. Nicolaou, F. Ptochos, P.A. Razis Charles University, Prague, Czech Republic M. Finger, M. Finger Jr. Academy of Scientific Research and Technology of the Arab Republic of Egypt, Egyptian Network of High Energy Physics, Cairo, Egypt Y. Assran9, S. Elgammal10, A. Ellithi Kamel11, A.M. Kuotb Awad12, M.A. Mahmoud12, A. Radi13,14 National Institute of Chemical Physics and Biophysics, Tallinn, Estonia M. Kadastik, M. Müntel, M. Murumaa, M. Raidal, L. Rebane, A. Tiko – 28 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 Department of Physics, University of Helsinki, Helsinki, Finland P. Eerola, G. Fedi, M. Voutilainen Helsinki Institute of Physics, Helsinki, Finland J. Härkönen, V. Karimäki, R. Kinnunen, M.J. Kortelainen, T. Lampén, K. Lassila-Perini, S. Lehti, T. Lindén, P. Luukka, T. Mäenpää, T. Peltola, E. Tuominen, J. Tuominiemi, E. Tuovinen, L. Wendland Lappeenranta University of Technology, Lappeenranta, Finland T. Tuuva DSM/IRFU, CEA/Saclay, Gif-sur-Yvette, France M. Besancon, F. Couderc, M. Dejardin, D. Denegri, B. Fabbro, J.L. Faure, F. Ferri, S. Ganjour, A. Givernaud, P. Gras, G. Hamel de Monchenault, P. Jarry, E. Locci, J. Malcles, L. Millischer, A. Nayak, J. Rander, A. Rosowsky, M. Titov Laboratoire Leprince-Ringuet, Ecole Polytechnique, IN2P3-CNRS, Palaiseau, France S. Baffioni, F. Beaudette, L. Benhabib, M. Bluj15, P. Busson, C. Charlot, N. Daci, T. Dahms, M. Dalchenko, L. Dobrzynski, A. Florent, R. Granier de Cassagnac, M. Hague- nauer, P. Miné, C. Mironov, I.N. Naranjo, M. Nguyen, C. Ochando, P. Paganini, D. Sabes, R. Salerno, Y. Sirois, C. Veelken, A. Zabi Institut Pluridisciplinaire Hubert Curien, Université de Strasbourg, Univer- sité de Haute Alsace Mulhouse, CNRS/IN2P3, Strasbourg, France J.-L. Agram16, J. Andrea, D. Bloch, J.-M. Brom, E.C. Chabert, C. Collard, E. Conte16, F. Drouhin16, J.-C. Fontaine16, D. Gelé, U. Goerlach, C. Goetzmann, P. Juillot, A.-C. Le Bihan, P. Van Hove Centre de Calcul de l’Institut National de Physique Nucleaire et de Physique des Particules, CNRS/IN2P3, Villeurbanne, France S. Gadrat Université de Lyon, Université Claude Bernard Lyon 1, CNRS-IN2P3, Institut de Physique Nucléaire de Lyon, Villeurbanne, France S. Beauceron, N. Beaupere, G. Boudoul, S. Brochet, J. Chasserat, R. Chierici, D. Contardo, P. Depasse, H. El Mamouni, J. Fay, S. Gascon, M. Gouzevitch, B. Ille, T. Kurca, M. Lethuillier, L. Mirabito, S. Perries, L. Sgandurra, V. Sordini, M. Vander Donckt, P. Verdier, S. Viret Institute of High Energy Physics and Informatization, Tbilisi State University, Tbilisi, Georgia Z. Tsamalaidze17 RWTH Aachen University, I. Physikalisches Institut, Aachen, Germany C. Autermann, S. Beranek, B. Calpas, M. Edelhoff, L. Feld, N. Heracleous, O. Hindrichs, K. Klein, A. Ostapchuk, A. Perieanu, F. Raupach, J. Sammet, S. Schael, D. Sprenger, H. Weber, B. Wittmer, V. Zhukov5 – 29 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 RWTH Aachen University, III. Physikalisches Institut A, Aachen, Germany M. Ata, J. Caudron, E. Dietz-Laursonn, D. Duchardt, M. Erdmann, R. Fischer, A. Güth, T. Hebbeker, C. Heidemann, K. Hoepfner, D. Klingebiel, S. Knutzen, P. Kreuzer, M. Merschmeyer, A. Meyer, M. Olschewski, K. Padeken, P. Papacz, H. Pieta, H. Reithler, S.A. Schmitz, L. Sonnenschein, J. Steggemann, D. Teyssier, S. Thüer, M. Weber RWTH Aachen University, III. Physikalisches Institut B, Aachen, Germany V. Cherepanov, Y. Erdogan, G. Flügge, H. Geenen, M. Geisler, W. Haj Ahmad, F. Hoehle, B. Kargoll, T. Kress, Y. Kuessel, J. Lingemann2, A. Nowack, I.M. Nugent, L. Perchalla, O. Pooth, A. Stahl Deutsches Elektronen-Synchrotron, Hamburg, Germany I. Asin, N. Bartosik, J. Behr, W. Behrenhoff, U. Behrens, A.J. Bell, M. Bergholz18, A. Bethani, K. Borras, A. Burgmeier, A. Cakir, L. Calligaris, A. Campbell, S. Choudhury, F. Costanza, C. Diez Pardos, S. Dooling, T. Dorland, G. Eckerlin, D. Eckstein, G. Flucke, A. Geiser, I. Glushkov, A. Grebenyuk, P. Gunnellini, S. Habib, J. Hauk, G. Hellwig, D. Horton, H. Jung, M. Kasemann, P. Katsas, C. Kleinwort, H. Kluge, M. Krämer, D. Krücker, E. Kuznetsova, W. Lange, J. Leonard, K. Lipka, W. Lohmann18, B. Lutz, R. Mankel, I. Marfin, I.-A. Melzer-Pellmann, A.B. Meyer, J. Mnich, A. Mussgiller, S. Naumann-Emme, O. Novgorodova, F. Nowak, J. Olzem, H. Perrey, A. Petrukhin, D. Pitzl, R. Placakyte, A. Raspereza, P.M. Ribeiro Cipriano, C. Riedl, E. Ron, M.Ö. Sahin, J. Salfeld-Nebgen, R. Schmidt18, T. Schoerner-Sadenius, N. Sen, M. Stein, R. Walsh, C. Wissing University of Hamburg, Hamburg, Germany M. Aldaya Martin, V. Blobel, H. Enderle, J. Erfle, E. Garutti, U. Gebbert, M. Görner, M. Gosselink, J. Haller, K. Heine, R.S. Höing, G. Kaussen, H. Kirschenmann, R. Klanner, R. Kogler, J. Lange, I. Marchesini, T. Peiffer, N. Pietsch, D. Rathjens, C. Sander, H. Schet- tler, P. Schleper, E. Schlieckau, A. Schmidt, M. Schröder, T. Schum, M. Seidel, J. Sibille19, V. Sola, H. Stadie, G. Steinbrück, J. Thomsen, D. Troendle, E. Usai, L. Vanelderen Institut für Experimentelle Kernphysik, Karlsruhe, Germany C. Barth, C. Baus, J. Berger, C. Böser, E. Butz, T. Chwalek, W. De Boer, A. De- scroix, A. Dierlamm, M. Feindt, M. Guthoff2, F. Hartmann2, T. Hauth2, H. Held, K.H. Hoffmann, U. Husemann, I. Katkov5, J.R. Komaragiri, A. Kornmayer2, P. Lobelle Pardo, D. Martschei, Th. Müller, M. Niegel, A. Nürnberg, O. Oberst, J. Ott, G. Quast, K. Rabbertz, F. Ratnikov, S. Röcker, F.-P. Schilling, G. Schott, H.J. Simonis, F.M. Stober, R. Ulrich, J. Wagner-Kuhr, S. Wayand, T. Weiler, M. Zeise Institute of Nuclear and Particle Physics (INPP), NCSR Demokritos, Aghia Paraskevi, Greece G. Anagnostou, G. Daskalakis, T. Geralis, S. Kesisoglou, A. Kyriakis, D. Loukas, A. Markou, C. Markou, E. Ntomari, I. Topsis-giotis University of Athens, Athens, Greece L. Gouskos, A. Panagiotou, N. Saoulidou, E. Stiliaris – 30 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 University of Ioánnina, Ioánnina, Greece X. Aslanoglou, I. Evangelou, G. Flouris, C. Foudas, P. Kokkas, N. Manthos, I. Papadopou- los, E. Paradas KFKI Research Institute for Particle and Nuclear Physics, Budapest, Hungary G. Bencze, C. Hajdu, P. Hidas, D. Horvath20, F. Sikler, V. Veszpremi, G. Vesztergombi21, A.J. Zsigmond Institute of Nuclear Research ATOMKI, Debrecen, Hungary N. Beni, S. Czellar, J. Molnar, J. Palinkas, Z. Szillasi University of Debrecen, Debrecen, Hungary J. Karancsi, P. Raics, Z.L. Trocsanyi, B. Ujvari National Institute of Science Education and Research, Bhubaneswar, India S.K. Swain22 Panjab University, Chandigarh, India S.B. Beri, V. Bhatnagar, N. Dhingra, R. Gupta, M. Kaur, M.Z. Mehta, M. Mittal, N. Nishu, A. Sharma, J.B. Singh University of Delhi, Delhi, India Ashok Kumar, Arun Kumar, S. Ahuja, A. Bhardwaj, B.C. Choudhary, S. Malhotra, M. Naimuddin, K. Ranjan, P. Saxena, V. Sharma, R.K. Shivpuri Saha Institute of Nuclear Physics, Kolkata, India S. Banerjee, S. Bhattacharya, K. Chatterjee, S. Dutta, B. Gomber, Sa. Jain, Sh. Jain, R. Khurana, A. Modak, S. Mukherjee, D. Roy, S. Sarkar, M. Sharan, A.P. Singh Bhabha Atomic Research Centre, Mumbai, India A. Abdulsalam, D. Dutta, S. Kailas, V. Kumar, A.K. Mohanty2, L.M. Pant, P. Shukla, A. Topkar Tata Institute of Fundamental Research - EHEP, Mumbai, India T. Aziz, R.M. Chatterjee, S. Ganguly, S. Ghosh, M. Guchait23, A. Gurtu24, G. Kole, S. Ku- mar, M. Maity25, G. Majumder, K. Mazumdar, G.B. Mohanty, B. Parida, K. Sudhakar, N. Wickramage26 Tata Institute of Fundamental Research - HECR, Mumbai, India S. Banerjee, S. Dugad Institute for Research in Fundamental Sciences (IPM), Tehran, Iran H. Arfaei, H. Bakhshiansohi, S.M. Etesami27, A. Fahim28, A. Jafari, M. Khakzad, M. Mohammadi Najafabadi, S. Paktinat Mehdiabadi, B. Safarzadeh29, M. Zeinali University College Dublin, Dublin, Ireland M. Grunewald – 31 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 INFN Sezione di Bari a, Università di Bari b, Politecnico di Bari c, Bari, Italy M. Abbresciaa,b, L. Barbonea,b, C. Calabriaa,b, S.S. Chhibraa,b, A. Colaleoa, D. Creanzaa,c, N. De Filippisa,c, M. De Palmaa,b, L. Fiorea, G. Iasellia,c, G. Maggia,c, M. Maggia, B. Marangellia,b, S. Mya,c, S. Nuzzoa,b, N. Pacificoa, A. Pompilia,b, G. Pugliesea,c, G. Selvaggia,b, L. Silvestrisa, G. Singha,b, R. Vendittia,b, P. Verwilligena, G. Zitoa INFN Sezione di Bologna a, Università di Bologna b, Bologna, Italy G. Abbiendia, A.C. Benvenutia, D. Bonacorsia,b, S. Braibant-Giacomellia,b, L. Brigliadoria,b, R. Campaninia,b, P. Capiluppia,b, A. Castroa,b, F.R. Cavalloa, G. Codispotia,b, M. Cuffiania,b, G.M. Dallavallea, F. Fabbria, A. Fanfania,b, D. Fasanellaa,b, P. Giacomellia, C. Grandia, L. Guiduccia,b, S. Marcellinia, G. Masettia, M. Meneghellia,b, A. Montanaria, F.L. Navarriaa,b, F. Odoricia, A. Perrottaa, F. Primaveraa,b, A.M. Rossia,b, T. Rovellia,b, G.P. Sirolia,b, N. Tosia,b, R. Travaglinia,b INFN Sezione di Catania a, Università di Catania b, Catania, Italy S. Albergoa,b, M. Chiorbolia,b, S. Costaa,b, F. Giordanoa,2, R. Potenzaa,b, A. Tricomia,b, C. Tuvea,b INFN Sezione di Firenze a, Università di Firenze b, Firenze, Italy G. Barbaglia, V. Ciullia,b, C. Civininia, R. D’Alessandroa,b, E. Focardia,b, S. Frosalia,b, E. Galloa, S. Gonzia,b, V. Goria,b, P. Lenzia,b, M. Meschinia, S. Paolettia, G. Sguazzonia, A. Tropianoa,b INFN Laboratori Nazionali di Frascati, Frascati, Italy L. Benussi, S. Bianco, F. Fabbri, D. Piccolo INFN Sezione di Genova a, Università di Genova b, Genova, Italy P. Fabbricatorea, R. Ferrettia,b, F. Ferroa, M. Lo Veterea,b, R. Musenicha, E. Robuttia, S. Tosia,b INFN Sezione di Milano-Bicocca a, Università di Milano-Bicocca b, Milano, Italy A. Benagliaa, M.E. Dinardoa,b, S. Fiorendia,b, S. Gennaia, A. Ghezzia,b, P. Govonia,b, M.T. Lucchinia,b,2, S. Malvezzia, R.A. Manzonia,b,2, A. Martellia,b,2, D. Menascea, L. Moronia, M. Paganonia,b, D. Pedrinia, S. Ragazzia,b, N. Redaellia, T. Tabarelli de Fatisa,b INFN Sezione di Napoli a, Università di Napoli ’Federico II’ b, Università della Basilicata (Potenza) c, Università G. Marconi (Roma) d, Napoli, Italy S. Buontempoa, N. Cavalloa,c, A. De Cosaa,b, F. Fabozzia,c, A.O.M. Iorioa,b, L. Listaa, S. Meolaa,d,2, M. Merolaa, P. Paoluccia,2 INFN Sezione di Padova a, Università di Padova b, Università di Trento (Trento) c, Padova, Italy P. Azzia, N. Bacchettaa, M. Bellatoa, D. Biselloa,b, A. Brancaa,b, R. Carlina,b, P. Checchiaa, T. Dorigoa, U. Dossellia, F. Fanzagoa, M. Galantia,b,2, F. Gasparinia,b, U. Gasparinia,b, P. Giubilatoa,b, A. Gozzelinoa, K. Kanishcheva,c, S. Lacapraraa, I. Lazzizzeraa,c, – 32 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 M. Margonia,b, A.T. Meneguzzoa,b, J. Pazzinia,b, M. Pegoraroa, N. Pozzobona,b, P. Ronchesea,b, F. Simonettoa,b, E. Torassaa, M. Tosia,b, A. Triossia, P. Zottoa,b, A. Zucchettaa,b, G. Zumerlea,b INFN Sezione di Pavia a, Università di Pavia b, Pavia, Italy M. Gabusia,b, S.P. Rattia,b, C. Riccardia,b, P. Vituloa,b INFN Sezione di Perugia a, Università di Perugia b, Perugia, Italy M. Biasinia,b, G.M. Bileia, L. Fanòa,b, P. Laricciaa,b, G. Mantovania,b, M. Menichellia, A. Nappia,b†, F. Romeoa,b, A. Sahaa, A. Santocchiaa,b, A. Spieziaa,b INFN Sezione di Pisa a, Università di Pisa b, Scuola Normale Superiore di Pisa c, Pisa, Italy K. Androsova,30, P. Azzurria, G. Bagliesia, T. Boccalia, G. Broccoloa,c, R. Castaldia, M.A. Cioccia, R.T. D’Agnoloa,c,2, R. Dell’Orsoa, F. Fioria,c, L. Foàa,c, A. Giassia, M.T. Grippoa,30, A. Kraana, F. Ligabuea,c, T. Lomtadzea, L. Martinia,30, A. Messineoa,b, C.S. Moona, F. Pallaa, A. Rizzia,b, A. Savoy-Navarroa,31, A.T. Serbana, P. Spagnoloa, P. Squillaciotia, R. Tenchinia, G. Tonellia,b, A. Venturia, P.G. Verdinia, C. Vernieria,c INFN Sezione di Roma a, Università di Roma b, Roma, Italy L. Baronea,b, F. Cavallaria, D. Del Rea,b, M. Diemoza, M. Grassia,b, E. Longoa,b, F. Margarolia,b, P. Meridiania, F. Michelia,b, S. Nourbakhsha,b, G. Organtinia,b, R. Paramattia, S. Rahatloua,b, C. Rovellia, L. Soffia,b INFN Sezione di Torino a, Università di Torino b, Università del Piemonte Orientale (Novara) c, Torino, Italy N. Amapanea,b, R. Arcidiaconoa,c, S. Argiroa,b, M. Arneodoa,c, R. Bellana,b, C. Biinoa, N. Cartigliaa, S. Casassoa,b, M. Costaa,b, A. Deganoa,b, N. Demariaa, C. Mariottia, S. Masellia, E. Migliorea,b, V. Monacoa,b, M. Musicha, M.M. Obertinoa,c, N. Pastronea, M. Pelliccionia,2, A. Potenzaa,b, A. Romeroa,b, M. Ruspaa,c, R. Sacchia,b, A. Solanoa,b, A. Staianoa, U. Tamponia INFN Sezione di Trieste a, Università di Trieste b, Trieste, Italy S. Belfortea, V. Candelisea,b, M. Casarsaa, F. Cossuttia,2, G. Della Riccaa,b, B. Gobboa, C. La Licataa,b, M. Maronea,b, D. Montaninoa,b, A. Penzoa, A. Schizzia,b, A. Zanettia Kangwon National University, Chunchon, Korea S. Chang, T.Y. Kim, S.K. Nam Kyungpook National University, Daegu, Korea D.H. Kim, G.N. Kim, J.E. Kim, D.J. Kong, S. Lee, Y.D. Oh, H. Park, D.C. Son Chonnam National University, Institute for Universe and Elementary Particles, Kwangju, Korea J.Y. Kim, Zero J. Kim, S. Song Korea University, Seoul, Korea S. Choi, D. Gyun, B. Hong, M. Jo, H. Kim, T.J. Kim, K.S. Lee, S.K. Park, Y. Roh – 33 – J H E P 1 0 ( 2 0 1 3 ) 1 6 7 University of Seoul, Seoul, Korea M. Choi, J.H. Kim, C. Park, I.C. Park, S. Park, G. Ryu Sungkyunkwan University, Suwon, Korea Y. Choi, Y.K. Choi, J. Goh, M.S. Kim, E. Kwon, B. Lee, J. Lee, S. Lee, H. Seo, I. Yu Vilnius University, Vilnius, Lithuania I. Grigelionis, A. Juodagalvis Centro de Investigacion y de Estudios Avanzados de