Abstract
This paper presents a measurement of the polarisation of W bosons from \(t\bar{t}\) decays, reconstructed in events with one high-\(p_{\text{ T }}\) lepton and at least four jets. Data from pp collisions at the LHC were collected at \(\sqrt{s}\) = 8 TeV and correspond to an integrated luminosity of 20.2 fb\(^{-1}\). The angle \(\theta ^{*}\) between the b-quark from the top quark decay and a direct W boson decay product in the W boson rest frame is sensitive to the W boson polarisation. Two different W decay products are used as polarisation analysers: the charged lepton and the down-type quark for the leptonically and hadronically decaying W boson, respectively. The most precise measurement of the W boson polarisation via the distribution of \(\cos {\theta ^{*}}\) is obtained using the leptonic analyser and events in which at least two of the jets are tagged as b-quark jets. The fitted fractions of longitudinal, left- and right-handed polarisation states are \(F_{\mathrm {0}}=~0.709~\pm ~{0.019}\), \(F_{\mathrm {L}}=~0.299~\pm ~{0.015}\) and \(F_{\mathrm {R}}=~-0.008~\pm ~{0.014}\), and are the most precisely measured W boson polarisation fractions to date. Limits on anomalous couplings of the Wtb vertex are set.
Similar content being viewed by others
1 Introduction
The top quark, discovered in 1995 by the CDF and D0 collaborations [1, 2] is the heaviest known elementary particle. It decays almost exclusively into a \(W\) boson and a \(b\text {-quark}\). The properties of the top decay vertex Wtb are determined by the structure of the weak interaction. In the Standard Model (SM) this interaction has a (\(V - A\)) structure, where V and A refer to the vector and axial vector components of the weak coupling. The \(W\) boson, which is produced as a real particle in the decay of top quarks, possesses a polarisation which can be left-handed, right-handed or longitudinal. The corresponding fractions, referred to as helicity fractions, are determined by the Wtb vertex structure and the masses of the particles involved. Calculations at next-to-next-to-leading order (NNLO) in QCD predict the fractions to be \(F_{\mathrm {L}} = 0.311 \pm 0.005\), \(F_{\mathrm {R}} = 0.0017 \pm 0.0001\), \(F_{\mathrm {0}} = 0.687 \pm 0.005\) [3].
By measuring the polarisation of the \(W\) boson with high precision, the SM prediction can be tested, and new physics processes which modify the structure of the Wtb vertex can be probed. The structure of the Wtb vertex can be expressed in a general form using left- and right-handed vector (\(V_{\mathrm {L/R}}\)) and tensor (\(g_{\mathrm {L/R}}\)) couplings:
Here, \(P_{\mathrm {L/R}}\) refer to the left- and right-handed chirality projection operators, \(m_W\) to the \(W\) boson mass, and g to the weak coupling constant. At tree level, all of the vector and tensor couplings vanish in the SM, except \(V_{\mathrm {L}}\), which corresponds to the CKM matrix element \(V_{tb}\) and has a value of approximately one. Dimension-six operators, introduced in effective field theories, can lead to anomalous couplings, represented by non-vanishing values of \(V_{\mathrm {R}}\), \(g_{\mathrm {L}}\) and \(g_{\mathrm {R}}\) [4,5,6].
The \(W\) boson helicity fractions can be accessed via angular distributions of polarisation analysers. Such analysers are \(W\) boson decay products whose angular distribution is sensitive to the \(W\) polarisation and determined by the Wtb vertex structure. In case of a leptonic decay of the \(W\) boson (\(W\rightarrow \ell \nu \)), the charged lepton serves as an ideal analyser: its reconstruction efficiency is very high and the sensitivity of its angular distribution to the \(W\) boson polarisation is maximal due to its weak isospin component \(T_3 = -\frac{1}{2}\). If the \(W\) boson decays hadronically (\(W\rightarrow q\bar{q}^{\prime }\)), the down-type quark is used, as it carries the same weak isospin as the charged lepton. This provides it with the same analysing power as the charged lepton, which is only degraded by the lower reconstruction efficiency and resolution of jets compared to charged leptons. The reconstruction of the down-type quark is in particular difficult as the two decay products of a hadronically decaying \(W\) boson are experimentally hard to separate. In the \(W\) boson rest frame, the differential cross-section of the analyser follows the distribution
which directly relates the \(W\) boson helicity fractions \(F_{i}\) to the angle \(\theta ^{*}\) between the analyser and the reversed direction of flight of the \(b\text {-quark}\) from the top quark decay in the W boson rest frame. Previous measurements of the \(W\) boson helicity fractions from the ATLAS, CDF, CMS and D0 collaborations show agreement with the SM within the uncertainties [7,8,9,10,11].
In this paper, the \(W\) boson helicity fractions are measured in top quark pair (\(t\bar{t}\)) events. Data corresponding to an integrated luminosity of 20.2 \(\text{ fb }^{-1}\)of proton–proton (pp) collisions, produced at the LHC with a centre-of-mass energy of \(\sqrt{s}\) = 8 \(\text {TeV}\), and recorded with the ATLAS [12] detector, are analysed. The final state of the \(t\bar{t}\) events is characterised by the decay of the \(W\) bosons. This analysis considers the lepton+jets channel in which one of the \(W\) bosons decays leptonically and the other decays hadronically. Both \(W\) boson decay modes are utilised for the measurement of \(\cos {\theta ^{*}}\). The signal selection and reconstruction includes direct decays of the \(W\) boson into an electron or muon as well as \(W\) boson decays into a \(\tau \)-lepton which subsequently decays leptonically.
2 The ATLAS detector
The ATLAS experiment at the LHC is a multi-purpose particle detector with a forward-backward symmetric cylindrical geometry and a near \(4\pi \) coverage in solid angle.Footnote 1 It consists of an inner tracking detector surrounded by a thin superconducting solenoid providing a 2 T axial magnetic field, electromagnetic and hadron calorimeters, and a muon spectrometer. The inner tracking detector covers the pseudorapidity range \(|\eta | < 2.5\). It consists of silicon pixel, silicon microstrip, and transition-radiation tracking detectors. Lead/liquid-argon (LAr) sampling calorimeters provide electromagnetic energy measurements with high granularity. A hadron (steel/scintillator-tile) calorimeter covers the central pseudorapidity range (\(|\eta | < 1.7\)). The end-cap and forward regions are instrumented with LAr calorimeters for electromagnetic and hadronic energy measurements up to \(|\eta | = 4.9\). The muon spectrometer surrounds the calorimeters and is based on three large air-core toroid superconducting magnets with eight coils each. Its bending power ranges from 2.0 to 7.5 T m. It includes a system of precision tracking chambers and fast detectors for triggering. A three-level trigger system is used to select events. The first-level trigger is implemented in hardware and uses a subset of the detector information to reduce the accepted rate to at most 75 kHz. This is followed by the high-level trigger, two software-based trigger levels that together reduce the accepted event rate to 400 Hz on average depending on the data-taking conditions.
3 Data and simulated samples
The data set consists of pp collisions, recorded at the LHC with \(\sqrt{s}\) = 8 \(\text {TeV}\), and corresponds to an integrated luminosity of \(20.2~\text{ fb }^{-1}\). Single-lepton triggers with a threshold of 24 \(\text {GeV}\) of transverse momentum (energy) for isolated muons (electrons) and 36 (60) \(\text {GeV}\) for muons (electrons) without an isolation criterion are used to select \(t\bar{t}\) candidate events. The lower trigger thresholds include isolation requirements on the candidate lepton, resulting in inefficiencies at high \(p_{\text {T}}\) that are recovered by the triggers with higher \(p_{\text {T}}\) thresholds.
Samples obtained from Monte Carlo (MC) simulations are used to characterise the detector response and reconstruction efficiency of \(t\bar{t}\) events, estimate systematic uncertainties and predict the background contributions from various processes. The response of the full ATLAS detector is simulated [13] using Geant 4 [14]. For the estimation of some systematic uncertainties, generated samples are passed through a faster simulation with parameterised showers in the calorimeters [15], while still using the full simulation of the tracking systems. Simulated events include the effect of multiple pp collisions from the same and nearby bunch-crossings (in-time and out-of-time pile-up) and are reweighted to match the number of collisions observed in data. All simulated samples are normalised using the most precise cross-section calculations available.
Signal \(t\bar{t}\) events are generated using the next-to-leading-order (NLO) QCD MC event generator Powheg-Box [16,17,18,19] using the CT10 parton distribution function (PDF) set [20]. Powheg-Box is interfaced to Pythia 6.425 [21] (referred to as the Powheg+Pythia sample), which is used to model the showering and hadronisation, with the CTEQ6L1 PDF set [22] and a set of tuned parameters called the Perugia2011C tune [23] for the modelling of the underlying event. The model parameter \(h_{\text {damp}}\) is set to \(m_{t}\) and controls matrix element to parton shower matching in Powheg-Box and effectively regulates the amount of high-\(p_{\text {T}}\) radiation.
The \(t\bar{t}\) cross-section is \(\sigma (t\bar{t}) = 253^{+13}_{-15}\) pb. This value is the result of a NNLO QCD calculation that includes resummation of next-to-next-to-leading logarithmic soft gluon terms with top++2.0 [24,25,26,27,28,29,30].
A sample generated with Powheg-Box interfaced with Herwig 6.520 [31] using Jimmy 4.31 [32] to simulate the underlying event (referred to as the Powheg+Herwig sample) is compared to a Powheg+Pythia sample to assess the impact of the different parton shower models. For both the Powheg+Herwig sample and this alternate Powheg+Pythia sample, the \(h_{\text {damp}}\) parameter is set to infinity.
To estimate the uncertainty due to the choice of the MC event generator, an alternate \(t\bar{t}\) MC sample is produced with MC@NLO [33, 34] with the CT10 PDF set interfaced to Herwig 6.520 using the AUET2 tune [35] and the CT10 PDF set for showering and hadronisation. In addition, samples generated with Powheg-Box interfaced to Pythia with variations in the amount of QCD initial- and final-state radiation (ISR/FSR) are used to estimate the effect of such uncertainty. The factorisation and renormalisation scales and the \(h_{\text {damp}}\) parameter in Powheg-Box as well as the transverse momentum scale of the space-like parton-shower evolution in Pythia are varied within the constraints obtained from an ATLAS measurement of \(t\bar{t}\) production in association with jets [36].
Single-top-quark-processes for the t-channel, s-channel and Wt associated production are also simulated with Powheg-Box [37, 38] using the CT10 PDF set. The samples are interfaced to Pythia 6.425 with the CTEQ6L1 PDF set and the Perugia2011C underlying event tune. Overlaps between the \(t\bar{t}\) and Wt final states are removed [39]. The single-top-quark samples are normalised using the approximate NNLO theoretical cross-sections [40,41,42] calculated with the MSTW2008 NNLO PDF set [43, 44]. All \(t\bar{t}\) and single-top samples are generated assuming a top quark mass of 172.5 \(\text {GeV}\), compatible with the ATLAS measurement of \(m_{t} = 172.84 \pm 0.70\,\text {GeV}\) [45].
Events with a \(W\) or \(Z\) boson produced in association with jets are generated using the leading-order (LO) event generator Alpgen 2.14 [46] with up to five additional partons and the CTEQ6L1 PDF set, interfaced to Pythia 6.425 for the parton showering and hadronisation. Separate samples for W / Z+light-jets, \(W/Zb\bar{b}\)+jets, \(W/Zc\bar{c}\)+jets and Wc+jets were generated. A parton–jet matching scheme (“MLM matching”) [47] is employed to avoid double-counting of jets generated from the matrix element and the parton shower. Overlap between the \(W/ZQ\bar{Q} \left( Q=b,c \right) \) events generated at the matrix element level and those generated by the parton shower evolution of the W / Z+light-jets sample are removed with an angular separation algorithm. If the angular distance \(\Delta R \) between the heavy-quark pair is larger than 0.4, the matrix element prediction is used instead of the parton shower prediction. Event yields from the Z+jets background are normalised using their inclusive NNLO theoretical cross-sections [48]. The predictions of normalisation and flavour composition of the \(W\)+jets background are affected by large uncertainties. Hence, a data-driven technique is used to determine both the inclusive normalisation and the heavy-flavour fractions of this process. The approach followed exploits the fact that the \(W^{\pm }\) boson production is charge-asymmetric at a pp collider. The W boson charge asymmetry depends on the flavour composition of the sample. Thus, correction factors estimated from data are used to rescale the fractions of \(Wb\bar{b}/c\bar{c}\)+jets, Wc+jets and \(W+\)light-jets events in the MC simulation: \(K_{bb}\) = \(K_{cc}\) = 1.50 ± 0.11 (stat. + syst.), \(K_{c}\) = 1.07 ± 0.27 (stat. + syst.) and \(K_{\text {light}}\) = 0.80 ± 0.04 (stat. + syst.) [49].
Diboson samples (WW, ZZ, WZ) are generated using the Sherpa 1.4.1 [50] event generator with the CT10 PDF set, with massive b- and c-quarks and with up to three additional partons in the LO matrix elements. The yields of these backgrounds are normalised using their NLO QCD theoretical cross-sections [51].
Multijet events can contain jets misidentified as leptons or non-prompt leptons from hadron decays and hence satisfy the selection criteria of the lepton+jets topology. This source of background events is referred to as fake-lepton background and is estimated using a data-driven approach (“matrix method”) which is based on the measurement of lepton selection efficiencies using different identification and isolation criteria [52].
4 Event selection and \(t\bar{t}\) reconstruction
4.1 Object reconstruction
The final state contains electrons, muons, jets with some of them originating from b-quarks, as well as missing transverse momentum.
Electrons are reconstructed from energy depositions in the electromagnetic calorimeter matching tracks in the inner detector. The transverse component of the energy deposition has to exceed 25 \(\text {GeV}\) and the pseudorapidity of the energy cluster, \(\eta _{\text {cluster}}\), has to fullfil \(|\eta _{\text {cluster}}| < 2.47\), excluding the transition region between the barrel and end-cap sections of the electromagnetic calorimeter at \(1.37< |\eta _{\text {cluster}}|\ < 1.52\). Electrons are further required to have a longitudinal impact parameter with respect to the hard-scattering vertex of less than 2 mm.
To reduce the background from non-prompt electrons (i.e. electrons produced within jets), electron candidates are also required to be isolated. Two \(\eta \)-dependent isolation criteria are applied. The first one considers the energy deposited in the calorimeter cells within a cone of size \(\Delta R = 0.2\) around the electron direction. The second one sums the transverse momenta (\(p_{\text {T}}\)) of all tracks with \(p_{\text {T}}\) > 400 \(\text {MeV}\) within a cone of size \(\Delta R = 0.3\) around the electron track. For each quantity, the transverse energy or momentum of the electron are subtracted. The isolation requirement is applied in such a way as to retain 90% of signal electrons, independent of their \(p_{\text {T}}\) value. This constant efficiency is verified in a data sample of \(Z\rightarrow ee\) decays [53].
For the reconstruction of muons, information from the muon spectrometer and the inner detector is combined. The combined muon track must satisfy \(p_{\text {T}} > 25\,\text {GeV}\) and \(|\eta | < 2.5\). The longitudinal impact parameter with respect to the hard-scattering vertex (defined in next section) is required to be less than 2 mm. Furthermore, muons are required to satisfy a \(p_{\text {T}} \)-dependent track-based isolation requirement. The scalar sum of the track \(p_{\text {T}} \) in a cone of variable size \(\Delta R < 10\,\text {GeV}/p_{\text {T}} ^\mu \) around the muon (excluding the muon track itself) has to be less than 5% of the muon \(p_{\text {T}} \).
Jets are reconstructed from topological clusters [12] built from energy depositions in the calorimeters using the anti-\(k_{t}\) algorithm [54, 55] with a radius parameter of 0.4. Before being processed by the jet-finding algorithm, the topological cluster energies are corrected using a local calibration scheme [56, 57] to account for inactive detector material, out-of-cluster leakage and the noncompensating calorimeter response. After energy calibration [58], the jets are required to have \(p_{\text {T}} > 25\,\text {GeV}\) and \(|\eta | < 2.5\). To suppress jets from pile-up, the jet vertex fractionFootnote 2 is required to be above 0.5 for all jets with \(p_{\text {T}}\) < 50 \(\text {GeV}\) and \(|\eta |\) < 2.4. As all electron candidates are also reconstructed as jets, the closest jet within a cone of size \(\Delta R \) = 0.2 around an electron candidate is discarded to avoid double-counting of electrons as jets. After this removal procedure, electrons within \(\Delta R \) = 0.4 of any remaining jet are removed.
Jets are identified as originating from the hadronisation of a \(b\text {-quark}\) (b-tagged) via a multivariate algorithm [59]. It makes use of the lifetime and mass of b-hadrons and accounts for displaced tracks and topological properties of the jets. A working point with 70% efficiency to tag a \(b\text {-quark}\) jet (\(b\)-jet) is used. The rejection factor for light-quark and gluon jets (light jets) is around 130 and about 5 for charm jets, as determined for b-tagged jets with \(p_{\text {T}}\) > 20 \(\text {GeV}\) and \(|\eta |<2.5\) in simulated \(t\bar{t}\) events. The simulated b-tagging efficiency is corrected to that measured in data using calibrations from statistically independent event samples of \(t\bar{t}\) pairs decaying into a \(b\bar{b}\ell ^{+}\ell ^{-}\nu _{\ell }\bar{\nu _{\ell }}\) final state [60].
The reconstruction of the transverse momentum of the neutrino from the leptonically decaying \(W\) boson is based on the negative vector sum of all energy deposits and momenta of reconstructed and calibrated objects in the transverse plane (missing transverse momentum with magnitude \(E_{\text {T}}^{\text {miss}}\)) as well as unassociated energy depositions [61].
4.2 Event selection
Events are selected from data taken in stable beam conditions with all relevant detector components being functional. At least one primary collision vertex is required with at least five associated tracks with \(p_{\text {T}}\) > 400 \(\text {MeV}\). If more than one primary vertex is reconstructed, the one with the largest scalar sum of transverse momenta is selected as the hard-scattering vertex. If the event contains at least one jet with \(p_{\text {T}} \)>\( 20\,\text {GeV}\) that is identified as out-of-time activity from a previous pp collision or as calorimeter noise [62], the event is rejected.
In order to select events from \(t\bar{t}\) decays in the lepton+jets channel, exactly one reconstructed electron or muon with \(p_{\text {T}}\) \(>25\) \(\text {GeV}\) and at least four jets, of which at least one is b-tagged, are required. A match (\(\Delta R < 0.15\)) between the offline reconstructed electron or muon and the lepton reconstructed by the high-level trigger is required. The selected events are separated into two orthogonal b-tag regions: one region with exactly one b-tag and a second region with two or more b-tags. Thus, the data sample is split into four channels depending on the lepton flavour and the b-jet multiplicity: “e+jets, 1 b-tag”, “e+jets, \(\ge \)2 b-tags”, “\(\mu \)+jets, 1 b-tag” and “\(\mu \)+jets, \(\ge \)2 b-tags”.
For events with one \(b\)-tag, \(E_{\text {T}}^{\text {miss}}\) is required to be greater than 20 \(\text {GeV}\) and the sum of \(E_{\text {T}}^{\text {miss}}\) and transverse mass of the leptonically decaying \(W\) boson, \(m_{\text {T}}(W)\), is required to be greater than 60 \(\text {GeV}\) in order to suppress multijet background. In the case of two \(b\)-tags, no further requirement on the \(E_{\text {T}}^{\text {miss}}\) and transverse mass of the \(W\) boson is applied.
After this selection, the \(t\bar{t}\) candidate events are reconstructed using a kinematic likelihood fit as described next.
4.3 Reconstruction of the \(t\bar{t}\) system
The measurement of the \(W\) boson polarisation in \(t\bar{t}\) events requires the reconstruction and identification of all \(t\bar{t}\) decay products. For this, a kinematic likelihood fitter (KLFitter) [63] is utilised. It maps the four model partons (two b-quarks and the \(q\bar{q}^{\prime }\) pair from a \(W\) boson decay) to four reconstructed jets. The numbers of jets used as input for KLFitter can be larger than four. The two jets with the largest output of the b-tagging algorithm together with two (three) remaining jets with the highest \(p_{\text {T}}\) were chosen as KLFitter input as this selection leads to the highest reconstruction efficiency for events with four (at least five) jets. For each of the \(4! = 24\) (\(5!=120\) for events with at least five jets) possible jet-to-parton permutations, it maximises a likelihood, \(\mathscr {L}\), that incorporates Breit–Wigner distributions for the \(W\) boson and top quark masses as well as transfer functions mapping the reconstructed jet and lepton energies to parton level or true lepton level, respectively. The expression for the likelihood is given by
where the \(BW(m_{ij(k)}| m_{t/W}\Gamma _{t/W})\) terms are the Breit-Wigner functions used to evaluate the mass of composite reconstructed particles (\(W\) bosons and top quarks) and \(W(E_{i}^{\text {meas}}|E_{j})\) are the transfer functions, with \(E_{i}^{\text {meas}}\) being the measured energy of object i and \(E_{j}\) the “true” energy of the reconstructed parton j or true lepton \(\ell \). The transverse components \(p^{x/y}_{\nu }\) of the neutrino momentum are mapped to the missing transverse momentum \(E^{\text {miss,} x/y}\) via transfer functions \(W(E^{\text {miss,} x/y}|p^{x/y}_{\nu })\). Individual transfer functions for electrons, muons, \(b\)-jets, light jets (including c-jets) and missing transverse momentum are used. These transfer functions are obtained from \(t\bar{t}\) events simulated with MC@NLO. The top quark decay products are uniquely matched to reconstructed objects to obtain a continuous function describing the relative energy difference between parton and reconstructed level as a function of the parton-level energy. Individual parameterisations are derived for different regions of \(|\eta |\). The measurement of the \(W\) boson polarisation in the lepton+jets channel is performed for both the top and the anti-top quarks in each event. The anti-down-type quark from the top quark decay (down-type quark from the anti-top quark decay) is used as the hadronic analyser and the charged lepton from the decay of the anti-top quark (charged anti-lepton from the top quark decay) as the leptonic analyser.
Since the likelihood defined in Eq. (3) is invariant under exchange of the \(W\) decay products, it needs further extensions to incorporate information related to down-type quarks. This is achieved by multiplying the likelihood by probability distributions of the \(b\)-tagging algorithm output as a function of the transverse momentum of the jets. These probability distributions are obtained from MC@NLO for b-quark jets as well as u / c- and d / s-quark jets. Since the \(W\) boson decays into a pair of charm and strange quarks in 50% of decays into hadrons, the higher values of the \(b\)-tagging algorithm output for the charm quark allows for a separation of the two. This increases the fraction of events with correct matching of the two jets originating from a \(W\) boson decay to the corresponding up- and down-quark type jet to 60%, compared to 50% for the case of no separation power. The extended likelihood is normalised with respect to the sum of the extended likelihoods for all 120 (24) permutations and this quantity is called the “event probability”. This up- versus down-type quark separation method was established in an ATLAS measurement of the \(t\bar{t}\) spin correlation in the lepton+jets channel [64].
The permutation with the largest event probability is chosen. Figure 1a shows the distributions of the logarithm of the likelihood value for the permutation with the highest event probability for simulated \(t\bar{t}\) events. Correctly reconstructed events (“\(t\bar{t}\) right”) peak at high values of the likelihood. Other contributions come from incorrect assignments of jets (i.e. choosing the wrong permutation, “\(t\bar{t}\) wrong”), non-reconstructable events where for example a quark is out of the acceptance (“\(t\bar{t}\) non-reco”) and \(t\bar{t}\) events which do not have a lepton+jets topology (such as dileptonic \(t\bar{t}\) events, “\(t\bar{t}\) background”). In Fig. 1b the corresponding distribution of the event probability is shown. The peak at 0.5 corresponds to events where no separation between up- and down-type quarks is achieved, leading to two permutations with similar event probabilities. High event probability indicates a correct down-type quark reconstruction.
To select the final data sample, the event probability is used to obtain the best jet-to-parton permutation per event. Events are required to have a reconstruction likelihood of \(\log {\mathscr {L}} > -48\) to reject poorly reconstructed \(t\bar{t}\) events. The value of \(\log {\mathscr {L}} > -48\) was selected to minimise the expected statistical uncertainty. The fraction of events where all jets were correctly assigned to the corresponding partons out of all events that have the corresponding jets present varies between 45 and 50%. The event yields after the final event selection are presented in Table 1.
Figure 2 shows the likelihood and the event probability as well as the reconstructed \(\cos {\theta ^{*}}\) distribution after the final event selection. Good agreement between data and prediction is achieved.
5 Measurement of the W boson helicity fractions
The W boson helicity fractions \(F_{i}\) are defined as the fraction of produced \(t\bar{t}\) events \(N_i\) in a given polarisation state divided by all produced \(t\bar{t}\) events:
The selection efficiency \(\epsilon _{i}^{\text {sel}}\) is different for each polarisation state and determines the number of selected events \(n_i\):
Dedicated \(t\bar{t}\) signal templates for a specific \(F_{i}\) are created by reweighting the simulated SM \(t\bar{t}\) events. These are produced by fitting the \(\cos {\theta ^{*}}\) distribution for the full phase space and calculating per-event weights for each helicity fraction using the functional forms in Eq. (2). Individual templates are created for each lepton flavour and b-tag channel. Figure 3 shows the templates for the \(\mu \) + jets channel with \(\ge \) 2 b-tags.
In addition to these signal templates, templates are derived for each source j of background events. These are independent of the helicity fractions \(F_{i}\). Five different background templates are included: three W+jets templates (W+light-jets, Wc+jets and \(Wc\bar{c}/b\bar{b}\)+jets), a fake-lepton template, and one template for all remaining backgrounds, including contributions from electroweak processes (single top, diboson and Z+jets). The total number of expected events \(n_{\text {exp}}\) in each channel is then given by
The signal and background templates are used to perform a likelihood fit with the number of background events \(n_{\text {bkg,}j}\) and the efficiency corrected signal events \(N_i\) as free parameters:
Here, \(n_{\text {data},k}\) represents the number of events in each bin k. The expected number of background events \(\hat{n}_{\text {bkg,}j}\) of each background source j and their normalisation uncertainties \(\sigma _{\text {bkg,}j}\) are used to constrain the fit. The fit parameters scaling the background contributions are treated as correlated across all channels except for the fake-lepton background, which is uncorrelated across lepton flavours and \(b\)-tag regions. The size of the background normalisation uncertainties \(\sigma _{\text {bkg,}j}\) is described in Sect. 6.
Combined fits of the \(\cos {\theta ^{*}}\) distributions using up to four different channels (\(e\) + jets and \(\mu \) + jets, both with 1 b-tag or \(\ge \) 2 b-tags) are performed for the leptonic and hadronic analyser individually. For each channel, individual templates of the signal and backgrounds are utilised. The combination leading to the lowest total uncertainty is used to quote the result. The helicity fractions are obtained from the fitted values of \(n_{i}\) using Eqs. 4–6. The fit method is validated using pseudo-experiments varying \(F_{\mathrm {0}}\) over the range [0.4, 1.0], \(F_{\mathrm {L}}\) over the range [0.15, 0.45] and \(F_{\mathrm {R}}\) over the range [−0.15, 0.15]. For each set, the unitarity constraint (\(F_{\mathrm {0}}\) + \(F_{\mathrm {L}}\) + \(F_{\mathrm {R}}\) = 1) is imposed. No bias is observed.
The uncertainties in the helicity fractions obtained from the fit include both the statistical uncertainty of the data and the systematic uncertainty of the background normalisations. For the leptonic analyser, the most sensitive results are obtained for the two-channel combination (electron + muon) in the \(\ge \)2 b-tags region. Adding further channels increases the total systematic uncertainty, in particular due to uncertainties in the b-tagging, which do not compensate with the decrease in the statistical uncertainty. For the hadronic analyser, the four-channel combination (including both the 1 b-tag and \(\ge \)2 b-tags regions) improves the sensitivity compared to the two-channel combination. For each source of systematic uncertainty, modified pseudo-data templates are created and evaluated via ensemble testing. The differences between the mean helicity fractions measured using the nominal templates and those varied to reflect systematic errors are quoted as systematic uncertainty. Systematic uncertainties from different sources, described in the following section, are treated as uncorrelated.
6 Systematic uncertainties
Systematic uncertainties from several sources can affect the normalisation of the signal and background and/or the shape of the \(\cos {\theta ^{*}}\) distribution. Correlations of a given systematic uncertainty are maintained across processes and channels, unless otherwise stated. The impact of uncertainties from the various sources is determined using a frequentist method based on the generation of pseudo-experiments.
6.1 Uncertainties associated with reconstructed objects
Different sources of systematic uncertainty affect the reconstructed objects used in this analyses. All these sources, described in the following, are propagated to changes in the shape of the \(\cos {\theta ^{*}}\) distributions.
Uncertainties associated with the lepton selection arise from the trigger, reconstruction, identification and isolation efficiencies, as well as the lepton momentum scale and resolution. They are estimated from \(Z \rightarrow \ell ^{+} \ell ^{-} (\ell = e,\,\mu )\), \(J/\psi \rightarrow \ell ^{+} \ell ^{-}\) and \(W \rightarrow e\nu \) processes in data and in simulated samples using tag-and-probe techniques described in Refs. [65,66,67,68,69]. Since small differences are observed between data and simulation, correction factors and their related uncertainties are considered to account for these differences. The effect of these uncertainties is propagated through the analysis and represent a minor source of uncertainty in this measurement.
Uncertainties associated with the jet selection arise from the jet energy scale, jet energy resolution, jet vertex fraction requirement and jet reconstruction efficiency. The jet energy scale and its uncertainty are derived combining information from test-beam data, LHC collision data, and simulation [58]. The jet energy scale uncertainty is split into 22 uncorrelated sources that have different jet \(p_{\text {T}}\) and \(\eta \) dependencies and are treated independently in this analysis. The uncertainty related to the jet energy resolution is estimated by smearing the energy of jets in simulation by the difference between the jet energy resolutions for data and simulation [70]. The efficiency for each jet to satisfy the jet vertex fraction requirement is measured in \(Z \rightarrow \ell ^{+} \ell ^{-}+\text {1-jet}\) events in data and simulation [71]. The corresponding uncertainty is evaluated in the analysis by changing the nominal jet vertex fraction cut value and repeating the analysis using the modified cut value [72]. The jet reconstruction efficiency is found to be about \(0.2\%\) lower in simulation than in data for jets below 30 \(\text {GeV}\) and consistent with data for higher jet \(p_T\). All jet-related kinematic variables (including the missing transverse momentum) are recomputed by removing randomly \(0.2\%\) of the jets with \(p_{\text {T}}\) below 30 \(\text {GeV}\) and the event selection is repeated.
Since the b-tagging efficiencies and misidentification rates are not modelled satisfactorily in MC simulation, all jets are assigned a specific \(p_{\text {T}}\)- and \(\eta \)-dependent scale factor to account for this difference. The uncertainties in these scale factors are propagated to the measured value.
An additional uncertainty is assigned due to the extrapolation of the b-tagging efficiency measurement to the high-\(p_{\text {T}}\) region. Twelve uncertainties are considered for the light-jet tagging, all depending on jet \(p_{\text {T}}\) and \(\eta \). These systematic uncertainties are taken as uncorrelated.
The uncertainties from the energy scale and resolution corrections for leptons and jets are propagated into the \(E_{\text {T}} ^{\text {miss}}\) calculation. Additional uncertainties are added to account for contributions from energy deposits not associated with any jet and due to soft-jets (7 GeV\(~<p_{\text {T}} <~20\) GeV), and are treated as fully correlated with each other. The uncertainty in the description of extra energy deposited due to pile-up interactions is treated as a separate \(E_{\text {T}} ^{\text {miss}}\) scale uncertainty. This uncertainty has a negligible effect on the measured \(W\) boson helicity fractions.
6.2 Uncertainties in signal modelling
The uncertainties in the signal modelling affect the kinematic properties of simulated \(t\bar{t}\) events and thus the acceptance and the shape of the reconstructed \(\cos {\theta ^{*}}\) distribution.
To assess the impact of the different parton shower and hadronisation models, the Powheg+Herwig sample is compared to a Powheg+Pythia sample and the symmetrised difference is taken as a systematic uncertainty. Similarly, an uncertainty due to the matrix element (ME) MC event generator choice for the hard process is estimated by comparing events produced by Powheg-Box and MC@NLO, both interfaced to Herwig for showering and hadronisation. The uncertainties due to QCD initial- and final-state radiation (ISR/FSR) modelling are estimated using two Powheg+Pythia samples with varied parameters producing more and less radiation. The larger of the changes due to the two variations is taken and symmetrised.
The uncertainty in the \(t\bar{t}\) signal due to the PDF choice is estimated following the PDF4LHC recommendations [73]. It takes into account the differences between three PDF sets: CT10 NLO, MSTW2008 68% CL NLO and NNPDF 2.3 NLO [74]. The final PDF uncertainty is an envelope of an intra-PDF uncertainty, which evaluates the changes due to the variation of different PDF parameters within a single PDF error set, and an inter-PDF uncertainty, which evaluates differences between different PDF sets. Each PDF set has a prescription to evaluate an overall uncertainty using its error sets: symmetric Hessian in the case of CT10, asymmetric Hessian for MSTW and sample standard deviation in the NNPDF case. Half the width of the envelope of the three estimates is taken as the PDF systematic uncertainty.
The effect of the uncertainty in the top quark mass is estimated using MC samples with different input top masses for the signal process. The dependence of the obtained helicity fractions on the top quark mass is fitted with a linear function. The uncertainties in the helicity fractions are obtained from the slopes multiplied by the uncertainty in the top quark mass of \(172.84 \pm 0.70\,\text {GeV}\) [45] measured by ATLAS at \(\sqrt{s}\) = 8 \(\text {TeV}\).
6.3 Uncertainties in background modelling
The different flavour samples of the W+jets background are scaled by data-driven calibration factors [49] as explained in Sect. 3. All sources of uncertainty on the correction factors other than normalisation (e.g. associated with the objects identification, reconstruction and calibration, etc.) are propagated to the W+jets estimation. Their normalisation uncertainty (5% for W+light-jets, 25% for W+c-jets and 7% for W+bb/cc) is taken into account in the likelihood fit as explained in Sect. 5.
A relative uncertainty of 30%, estimated using various control regions in the matrix method calculation [52], is used for the fake-lepton contribution.
For single top quark production, a normalisation uncertainty of 17% is assumed, which takes into account the weighted average of the theoretical uncertainties in s-, t- and Wt-channel production (+5/\(-4\)%) as well as additional uncertainties due to variations in the amount of initial- and final-state radiation and the extrapolation to high jet multiplicity. The uncertainty in the single-top background shape is assessed by comparing Wt-channel Monte Carlo samples generated using alternative methods to take into account Wt and \(t\bar{t}\) diagrams interference: diagram removal and diagram subtraction [39].
An overall normalisation uncertainty of 48% is applied to Z+jets and diboson contributions. It takes into account a 5% uncertainty in the theoretical (N)NLO cross-section as well as the uncertainty associated with the extrapolation to high jet multiplicity (24% per jet).
All normalisation uncertainties are included in the fit of the \(W\) boson helicity fractions via priors for the background yields. While the W+jets and fake-lepton uncertainties are included directly, the uncertainty in the total remaining background from other sources is combined to 16% (17%) in the \(\ge \)2 b-tags regions (1 b-tag + \(\ge \)2 b-tags regions) by adding the uncertainties in the theoretical cross-sections of the single top quark, diboson and Z+jets contributions in quadrature. The uncertainty in the shape of the W+jets background is considered by jet flavour decomposition. Further background shape uncertainties were evaluated and found to be negligible.
6.4 Other uncertainties
The uncertainty associated with the limited number of MC events in the signal and background templates is evaluated by performing pseudo-experiments on MC events.
The impact of the 1.9% luminosity uncertainty [75] is found to be negligible since the background normalisations are constrained in the fit.
7 Results
The measured \(W\) boson helicity fractions obtained using the leptonic analyser in semileptonic \(t\bar{t}\) events with \(\ge \)2 b-tags are presented in Table 2.
By construction, the individual fractions sum up to one. The \(F_{\mathrm {0}}\) value is anti-correlated with both \(F_{\mathrm {L}}\) and \(F_{\mathrm {R}}\) (\(\rho _{F_{\mathrm {0}}, F_{\mathrm {L}}}=-0.55\), \(\rho _{F_{\mathrm {0}}, F_{\mathrm {R}}}=-0.75)\), and \(F_{\mathrm {L}}\) and \(F_{\mathrm {R}}\) are positively correlated (\(\rho _{F_{\mathrm {L}}, F_{\mathrm {R}}}=+0.16)\). The quoted values correspond to the total correlation coefficient, considering statistical and systematic uncertainties. These results are the most precise \(W\) boson helicity fractions measured so far and are consistent with the SM predictions given at NNLO accuracy [3]. The inclusion of single b-tag regions does not improve the sensitivity, due to larger systematic uncertainties.
The \(W\) boson helicity fractions obtained using the hadronic analyser of semileptonic \(t\bar{t}\) events with 1 b-tag and \(\ge \)2b-tags are given in Table 3. Using the hadronic analyser, the correlations between the helicity fraction are \(\rho _{F_{\mathrm {0}}, F_{\mathrm {L}}}=0.56\), \(\rho _{F_{\mathrm {0}}, F_{\mathrm {R}}}=-0.91\) and \(\rho _{F_{\mathrm {L}}, F_{\mathrm {R}}}=-0.92\). The large anticorrelation between \(F_{\mathrm {L}}\) and \(F_{\mathrm {R}}\) is a consequence of the low separation power between the up- and down-type quark from the W decay and the resulting similar shapes of the templates of \(F_{\mathrm {L}}\) and \(F_{\mathrm {R}}\) (see Fig. 3). The results obtained with the two analysers agree well. The combination of leptonic and hadronic analysers has been tested and, despite the improvement in the statistical uncertainty, it does not improve the total uncertainty.
Figure 4 shows, separately for the e+jets and \(\mu \)+jets channels, the distributions of \(\cos {\theta ^{*}}\) from the leptonic analyser. The distributions for the hadronic analyser are presented in Fig. 5. The uncertainty band in the data-to-best-fit ratio represents the statistical and background normalisation uncertainty. The deviations observed in the ratio are covered by the systematic uncertainties. The peak at \(\cos {\theta ^{*}} \approx -0.7\) as seen in the single b-tag channels in Fig. 5 is caused by misreconstructed events. A missing second b-tag increases the probability of swapping the b-quark jet from the top quark decay with the up-type quark jet from the W decay.
The contributions of the various systematic uncertainties are quoted in Table 4. In the case of the leptonic analyser, the dominant contributions come from the jet energy scale and resolution and the statistical error in the MC templates. For the hadronic analyser, the systematic uncertainties are larger. Including the 1 b-tag region aids in reducing the error. One of the main contributions is the \(b\)-tagging uncertainty, affecting both the event selection and b-tag categorisation, as well as the up- vs down-type quark separation. Other major contributions come from the jet energy resolution and the modelling of \(t\bar{t}\) events (initial- and final-state radiation, parton showering and hadronisation, and Monte Carlo event generator choice for the matrix elements).
Within the effective field theory framework [76], the Wtb decay vertex can be parameterised in terms of anomalous couplings as shown in Eq. (1). Limits on these anomalous left- and right-handed vector and tensor couplings are set using the EFTfitter tool [77] and the model of [76]. The anomalous couplings are assumed to be real, corresponding to the CP-conserving case. As the \(W\) helicity fractions only allow the ratios of couplings to be constrained, the value of \(V_{\mathrm {L}}\) is fixed to the Standard Model prediction of one. The correlations of systematic uncertainties are taken into account. Figure 6 shows the limits on \(g_{\mathrm {L}}\) and \(g_{\mathrm {R}}\) couplings while \(V_{\mathrm {L}}\) and \(V_{\mathrm {R}}\) are fixed to their SM values, as well as \(V_{\mathrm {R}}\) and \(g_{\mathrm {R}}\) limits, where the other couplings are fixed to their SM values. The intervals are obtained using the leptonic analyser since it provides the most sensitive results. Table 5 shows the 95% confidence level (CL) intervals for each anomalous coupling while fixing all others to their SM value. These limits correspond to the set of smallest intervals containing 95% of the marginalised posterior distribution for the corresponding parameter.
Similar limits on the anomalous couplings were derived by both the ATLAS and CMS experiments using the measured helicity fractions of \(W\) bosons [10, 11]. Complementary limits can be set by other measurements: the allowed region of \(g_{\mathrm {R}} \approx 0.75\) is excluded by measurements of the t-channel single top quark production [77,78,79,80] which also constrains \(V_{\mathrm {L}}\). The branching fraction of \(\bar{B} \rightarrow X_s \gamma \) allow more stringent limits to be set on \(g_{\mathrm {L}}\) and \(V_{\mathrm {R}}\) [81].
8 Conclusion
The longitudinal, left- and right-handed \(W\) boson helicity fractions are measured using the angle between the charged lepton (down-type quark) and the reversed \(b\text {-quark}\) direction in the \(W\) boson rest frame for leptonically (hadronically) decaying \(W\) bosons from \(t\bar{t}\) decays. A data set corresponding to 20.2 \(\text{ fb }^{-1}\) of pp collisions at the LHC with a centre-of-mass energy of \(\sqrt{s}\) = 8 \(\text {TeV}\), recorded by the ATLAS experiment, is analysed. Events are required to include one isolated electron or muon and at least four jets, with at least one of them tagged as a \(b\)-jet. Events are reconstructed using a kinematic likelihood fit based on mass constraints for the top quarks and \(W\) bosons. It utilises the weight of the \(b\)-jet tagging algorithm to further separate the up- and down-type quarks from the hadronically decaying \(W\) bosons. The fractions for left-handed, right-handed and longitudinally polarised \(W\) bosons are found to be \(F_{\mathrm {0}}\) = \(0.709\) ± \(0.012\) (stat.+bkg. norm.) \({}\pm {0.015}\) (syst.), \(F_{\mathrm {L}}\) = \(0.299\) ± \(0.008\) (stat.+bkg. norm.) \({}\pm {0.013}\) (syst.) and \(F_{\mathrm {R}}\) = \(-0.008\) ± \(0.006\) (stat.+bkg. norm.) \({}\pm {0.012}\) (syst.). These results constitute the most precise measurement of the \(W\) helicity fractions in \(t\bar{t}\) events to date and are in good agreement with the Standard Model predictions within uncertainties. Using these results, limits on anomalous couplings of the Wtb vertex are set.
Change history
09 January 2019
This change does not have any impact on the measured helicity fractions, but it changes the obtained limits on the anomalous couplings.
Notes
ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y-axis points upwards. Cylindrical coordinates \((r,\phi )\) are used in the transverse plane, \(\phi \) being the azimuthal angle around the z-axis. The pseudorapidity is defined in terms of the polar angle \(\theta \) as \(\eta = -\ln \tan (\theta /2)\). The angular distance is measured in units of \(\Delta R \equiv \sqrt{(\Delta \eta )^{2} + (\Delta \phi )^{2}}\).
The jet vertex fraction is defined as the scalar sum of the transverse momenta of a jet’s tracks stemming from the primary collision vertex divided by the scalar sum of the transverse momenta of all tracks in a jet.
References
CDF Collaboration, F. Abe et al., Observation of top quark production in \(\bar{p}p\) collisions. Phys. Rev. Lett. 74, 2626 (1995). doi:10.1103/PhysRevLett. 74.2626. arXiv: hep-ex/9503002
D0 Collaboration, S. Abachi et al., Search for high mass top quark production in \(p\bar{p}\) collisions at \(\sqrt{s} = 1.8\) TeV, Phys. Rev. Lett. 74, 2422 (1995). doi:10.1103/PhysRevLett.74.2422. arXiv: hep-ex/9411001
A. Czarnecki, J.G. Korner, J.H. Piclum, Helicity fractions of W bosons from top quark decays at NNLO in QCD. Phys. Rev. D 81, 111503 (2010). doi:10.1103/PhysRevD.81.111503. arXiv:1005.2625 [hep-ph]
W. Buchmuller, D. Wyler, Effective Lagrangian analysis of new interactions and flavor conservation. Nucl. Phys. B 268, 621 (1986). doi:10.1016/0550-3213(86)90262-2
J.A. Aguilar-Saavedra, A minimal set of top anomalous couplings. Nucl. Phys. B 812, 181 (2009). doi:10.1016/j.nuclphysb.2008.12.012. arXiv:0811.3842 [hep-ph]
C. Zhang, S. Willenbrock, Effective-field-theory approach to top-quark production and decay. Phys. Rev. D 83, 034006 (2011). doi:10.1103/PhysRevD.83.034006. arXiv:1008.3869 [hep-ph]
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, 071106 (2012). doi:10.1103/PhysRevD.85.071106. arXiv:1202.5272 [hep-ex]
CDF Collaboration, T. Aaltonen et al., Measurement of W-Boson Polarization in Top-quark Decay in \(p\bar{p}\) Collisions at \(\sqrt{s}\) = 1.96 TeV. Phys. Rev. Lett 105, 042002 (2010). doi:10.1103/PhysRevLett.105.042002. arXiv:1003.0224 [hep-ex]
D0 Collaboration, V.M. Abazov et al., Measurement of the W boson helicity in top quark decays using 5.4 fb\(^{-1}\) of \(p\bar{p}\) collision data. Phys. Rev. D 83, 032009 (2011). doi:10.1103/PhysRevD.83.032009. arXiv:1011.6549 [hep-ex]
ATLAS Collaboration, Measurement of the \(W\) boson polarization in top quark decays with the ATLAS detector. JHEP 06, 088 (2012). doi:10.1007/JHEP06(2012)088. arXiv:1205.2484 [hep-ex]
CMS Collaboration, Measurement of the W boson helicity fractions in the decays of top quark pairs to lepton \(+\) jets final states produced in pp collisions at \(\sqrt{s}=\) 8TeV. Phys. Lett. B 762, 512 (2016). doi:10.1016/j.physletb.2016.10.007. arXiv:1605.09047 [hep-ex]
ATLAS Collaboration, The ATLAS experiment at the CERN large hadron collider. JINST 3, S08003 (2008). doi:10.1088/1748-0221/3/08/S08003
ATLAS Collaboration, The ATLAS simulation infrastructure. Eur. Phys. J. C 70, 823 (2010). doi:10.1140/epjc/s10052-010-1429-9. arXiv:1005.4568 [hep-ex]
S. Agostinelli et al., GEANT4: a simulation toolkit. Nucl. Instrum. Methods A 506, 250 (2003). doi:10.1016/S0168-9002(03)01368-8
ATLAS Collaboration, The simulation principle and performance of the ATLAS fast calorimeter simulation FastCaloSim. ATL-PHYS-PUB-2010-013,2010, http://cds.cern.ch/record/1300517
P. Nason, A new method for combining NLO QCD with shower Monte Carlo algorithms. JHEP 11, 040 (2004). doi:10.1088/1126-6708/2004/11/040. arXiv:hep-ph/0409146
S. Frixione, P. Nason, G. Ridolfi, A positive-weight next-to-leading-order Monte Carlo for heavy flavour hadroproduction. JHEP 09, 126 (2007). doi:10.1088/1126-6708/2007/09/126. arXiv:0707.3088 [hep-ph]
S. Frixione, P. Nason, C. Oleari, Matching NLO QCD computations with Parton shower simulations: the POWHEG method. JHEP 11, 070 (2007). doi:10.1088/1126-6708/2007/11/070. arXiv:0709.2092 [hep-ph]
S. Alioli, P. Nason, C. Oleari, E. Re, A general framework for implementing NLO calculations in shower Monte Carlo programs: the POWHEG BOX. JHEP 06, 043 (2010). doi:10.1007/JHEP06(2010)043. arXiv:1002.2581 [hep-ph]
H.-L. Lai et al., New parton distributions for collider physics. Phys. Rev. D 82, 074024 (2010). doi:10.1103/PhysRevD.82.074024. arXiv:1007.2241 [hep-ph]
T. Sjöstrand, S. Mrenna, P.Z. Skands, PYTHIA 6.4 physics and manual. JHEP 05, 026 (2006). doi:10.1088/1126-6708/2006/05/026. arXiv:hep-ph/0603175
J. Pumplin et al., New generation of parton distributions with uncertainties from global QCD analysis. JHEP 07, 012 (2002). doi:10.1088/1126-6708/2002/07/012. arXiv:hep-ph/0201195
P.Z. Skands, Tuning Monte Carlo generators: the perugia tunes. Phys. Rev. D 82, 074018 (2010). doi:10.1103/PhysRevD.82.074018. arXiv:1005.3457 [hep-ph]
M. Cacciari, M. Czakon, M. Mangano, A. Mitov, P. Nason, Top-pair production at hadron colliders with next-to-next-to-leading logarithmic soft-gluon resummation. Phys. Lett. B 710, 612 (2012). doi:10.1016/j.physletb.2012.03.013. arXiv:1111.5869 [hep-ph]
M. Beneke, P. Falgari, S. Klein, C. Schwinn, Hadronic top-quark pair production with NNLL threshold resummation. Nucl. Phys. B 855, 695 (2012). doi:10.1016/j.nuclphysb.2011.10.021. arXiv:1109.1536 [hep-ph]
P. Bärnreuther, M. Czakon, A. Mitov, Percent level precision physics at the tevatron: first genuine NNLO QCD corrections to \(q \bar{q} \rightarrow t \bar{t} + X\). Phys. Rev. Lett. 109, 132001 (2012). doi:10.1103/PhysRevLett.109.132001. arXiv:1204.5201 [hep-ph]
M. Czakon, A. Mitov, NNLO corrections to top pair production at hadron colliders: the quark-gluon reaction. JHEP 01, 080 (2013). doi:10.1007/JHEP01(2013)080. arXiv:1210.6832 [hep-ph]
M. Czakon, A. Mitov, NNLO corrections to top-pair production at hadron colliders: the all-fermionic scattering channels. JHEP 12, 054 (2012). doi:10.1007/JHEP12(2012)054. arXiv:1207.0236 [hep-ph]
M. Czakon, P. Fiedler, A. Mitov, Total top-quark pair-production cross section at hadron colliders through \(O(\alpha ^{4}_{S})\). Phys. Rev. Lett. 110, 252004 (2013). doi:10.1103/PhysRevLett.110.252004. arXiv:1303.6254 [hep-ph]
M. Czakon, A. Mitov, Top++: a program for the calculation of the top-pair cross-section at hadron colliders. Comput. Phys. Commun. 185, 2930 (2014). doi:10.1016/j.cpc.2014.06.021. arXiv:1112.5675 [hep-ph]
G. Corcella et al., HERWIG 6: an event generator for hadron emission reactions with interfering gluons (including supersymmetric processes). JHEP 01, 010 (2001). doi:10.1088/1126-6708/2001/01/010. arXiv:hep-ph/0011363
J.M. Butterworth, J.R. Forshaw, M.H. Seymour, Multiparton interactions in photoproduction at HERA. Z. Phys. C 72, 637 (1996). doi:10.1007/BF02909195, doi:DOIurl10.1007/s002880050286. arXiv:hep-ph/9601371
S. Frixione, B.R. Webber, Matching NLO QCD computations and parton shower simulations. JHEP 06, 029 (2002). doi:10.1088/1126-6708/2002/06/029. arXiv:hep-ph/0204244
S. Frixione, P. Nason, B.R. Webber, Matching NLO QCD and parton showers in heavy flavor production. JHEP 08, 007 (2003). doi:10.1088/1126-6708/2003/08/007. arXiv:hep-ph/0305252
ATLAS Collaboration, New ATLAS event generator tunes to 2010 data. ATL-PHYS-PUB-2011-008, 2011. http://cds.cern.ch/record/1345343
ATLAS Collaboration, Measurement of \(t \bar{t}\) production with a veto on additional central jet activity in \(pp\) collisions at \(\sqrt{s} = 7\;\text{TeV}\)using the ATLAS detector. Eur. Phys. J. C 72, 2043 (2012). doi:10.1140/epjc/s10052-012-2043-9. arXiv:1203.5015 [hep-ex]
S. Alioli, P. Nason, C. Oleari, E. Re, NLO single-top production matched with shower in POWHEG: s- and t-channel contributions. JHEP 09, 111 (2009). [Erratum: JHEP 02, 011 (2010)]. doi:10.1007/JHEP02(2010)011, doi:10.1088/1126-6708/2009/09/111. arXiv:0907.4076 [hep-ph]
E. Re, Single-top Wt-channel production matched with parton showers using the POWHEG method. Eur. Phys. J. C 71, 1547 (2011). doi:10.1140/epjc/s10052-011-1547-z. arXiv:1009.2450 [hep-ph]
S. Frixione, E. Laenen, P. Motylinski, B.R. Webber, C.D. White, Single-top hadroproduction in association with a W boson. JHEP 07, 029 (2008). doi:10.1088/1126-6708/2008/07/029. arXiv:0805.3067 [hep-ph]
N. Kidonakis, Next-to-next-to-leading-order collinear and soft gluon corrections for t-channel single top quark production. Phys. Rev. D 83, 091503 (2011). doi:10.1103/PhysRevD.83.091503. arXiv:1103.2792 [hep-ph]
N. Kidonakis, NNLL resummation for s-channel single top quark production. Phys. Rev. D 81, 054028 (2010). doi:10.1103/PhysRevD.81.054028. arXiv:1001.5034 [hep-ph]
N. Kidonakis, Two-loop soft anomalous dimensions for single top quark associated production with a W- or H-. Phys. Rev. D 82, 054018 (2010). doi:10.1103/PhysRevD.82.054018. arXiv:1005.4451 [hep-ph]
A.D. Martin, W.J. Stirling, R.S. Thorne, G. Watt, Parton distributions for the LHC. Eur. Phys. J. C 63, 189 (2009). doi:10.1140/epjc/s10052-009-1072-5. arXiv:0901.0002 [hep-ph]
A.D. Martin, W.J. Stirling, R.S. Thorne, G. Watt, Uncertainties on alpha(S) in global PDF analyses and implications for predicted hadronic cross sections. Eur. Phys. J. C 64, 653 (2009). doi:10.1140/epjc/s10052-009-1164-2. arXiv:0905.3531 [hep-ph]
ATLAS Collaboration, Measurement of the top quark mass in the \(t\bar{t}\rightarrow \) dilepton channel from \(\sqrt{s} = 8\;\text{ TeV }\) ATLAS data. Phys. Lett. B 761, 350 (2016). doi:10.1016/j.physletb.2016.08.042. arXiv:1606.02179 [hep-ex]
M.L. Mangano, M. Moretti, F. Piccinini, R. Pittau, A.D. Polosa, ALPGEN, a generator for hard multiparton processes in hadronic collisions. JHEP 07, 001 (2003). doi:10.1088/1126-6708/2003/07/001. arXiv:hep-ph/0206293
M.L. Mangano, M. Moretti, R. Pittau, Multijet matrix elements and shower evolution in hadronic collisions: \(W b \bar{b}\) + \(n\) jets as a case study. Nucl. Phys. B 632, 343 (2002). doi:10.1016/S0550-3213(02)00249-3. arXiv:hep-ph/0108069
K. Melnikov, F. Petriello, Electroweak gauge boson production at hadron colliders through O(alpha(s)**2). Phys. Rev. D 74, 114017 (2006). doi:10.1103/PhysRevD.74.114017. arXiv:hep-ph/0609070
ATLAS Collaboration, Measurement of the charge asymmetry in top-quark pair production in the lepton-plus-jets final state in \(pp\) collision data at \(\sqrt{s} = 8\;\text{ TeV }\) with the ATLAS detector. Eur. Phys. J. C 76, 87 (2016). doi:10.1140/epjc/s10052-016-3910-6. arXiv:1509.02358 [hep-ex]
T. Gleisberg et al., Event generation with SHERPA 1.1. JHEP 02, 007 (2009). doi:10.1088/1126-6708/2009/02/007. arXiv:0811.4622 [hep-ph]
J.M. Campbell, R.K. Ellis, An Update on vector boson pair production at hadron colliders. Phys. Rev. D 60, 113006 (1999). doi:10.1103/PhysRevD.60.113006. arXiv:hep-ph/9905386
ATLAS Collaboration, Estimation of non-prompt and fake lepton backgrounds in final states with top quarks produced in proton–proton collisions at \(\sqrt{s} = 8\;\text{ TeV }\) with the ATLAS Detector. ATLAS-CONF-2014-058, 2014. http://cds.cern.ch/record/1951336
ATLAS Collaboration, Electron efficiency measurements with the ATLAS detector using the 2012 LHC proton–proton collision data. ATLAS-CONF-2014-032, 2014. http://cds.cern.ch/record/1706245
M. Cacciari, G.P. Salam, G. Soyez, The anti-\(k_t\) jet clustering algorithm. JHEP 04, 063 (2008). doi:10.1088/1126-6708/2008/04/063. arXiv:0802.1189 [hep-ph]
M. Cacciari, G.P. Salam, Dispelling the \(N^{3}\) myth for the \(k_t\) jet-finder. Phys. Lett. B 641, 57 (2006). doi:10.1016/j.physletb.2006.08.037. arXiv:hep-ph/0512210
W. Lampl et al., Calorimeter Clustering Algorithms: Description and Performance. ATL-LARG-PUB-2008-002, 2008. http://cds.cern.ch/record/1099735
C. Cojocaru et al., Hadronic calibration of the ATLAS liquid argon end-cap calorimeter in the pseudorapidity region \(1.6 < |\eta | < 1.8\) in beam tests. Nucl. Instrum. Methods A 531, 481 (2004). doi:10.1016/j.nima.2004.05.133. arXiv:physics/0407009
ATLAS Collaboration, Jet energy measurement and its systematic uncertainty in proton–proton collisions at \(\sqrt{s} = 7\;\text{ TeV }\) with the ATLAS detector. Eur. Phys. J. C 75, 17 (2015). doi:10.1140/epjc/s10052-014-3190-y. arXiv:1406.0076 [hep-ex]
ATLAS Collaboration, Performance of \(b\)-jet identification in the ATLAS experiment. JINST 11, P04008 (2016). doi:10.1088/1748-0221/11/04/P04008. arXiv:1512.01094 [hep-ex]
ATLAS Collaboration, Calibration of \(b\)-tagging using dileptonic top pair events in a combinatorial likelihood approach with the ATLAS experiment. ATLAS-CONF-2014-004, 2014. http://cds.cern.ch/record/1664335
ATLAS Collaboration, Performance of missing transverse momentum reconstruction in proton–proton collisions at \(\sqrt{s} = 7\;\text{ TeV }\) with ATLAS. Eur. Phys. J. C 72, 1844 (2012). doi:10.1140/epjc/s10052-011-1844-6. arXiv:1108.5602 [hep-ex]
ATLAS Collaboration, Selection of jets produced in proton–proton collisions with the ATLAS detector using 2011 data. ATLAS-CONF-2012-020, 2012. http://cds.cern.ch/record/1430034
J. Erdmann et al., A likelihood-based reconstruction algorithm for top-quark pairs and the KLFitter framework. Nucl. Instrum. Methods A 748, 18 (2014). doi:10.1016/j.nima.2014.02.029. arXiv:1312.5595 [hep-ex]
ATLAS Collaboration, Measurements of spin correlation in top–antitop quark events from proton–proton collisions at \(\sqrt{s} = 7\;\text{ TeV }\) using the ATLAS detector. Phys. Rev. D 90, 112016 (2014). doi:10.1103/PhysRevD.90.112016. arXiv:1407.4314 [hep-ex]
ATLAS Collaboration, Electron reconstruction and identification efficiency measurements with the ATLAS detector using the 2011 LHC proton–proton collision data. Eur. Phys. J. C 74, 2941 (2014). doi:10.1140/epjc/s10052-014-2941-0. arXiv:1404.2240 [hep-ex]
ATLAS Collaboration, Muon reconstruction efficiency and momentum resolution of the ATLAS experiment in proton–proton collisions at \(\sqrt{s} = 7\;\text{ TeV }\) in 2010. Eur. Phys. J. C 74, 3034 (2014). doi:10.1140/epjc/s10052-014-3034-9. arXiv:1404.4562 [hep-ex]
ATLAS Collaboration, Electron performance measurements with the ATLAS detector using the 2010 LHC proton–proton collision data. Eur. Phys. J. C 72, 1909 (2012). doi:10.1140/epjc/s10052-012-1909-1. arXiv:1110.3174 [hep-ex]
ATLAS Collaboration, Measurement of the muon reconstruction performance of the ATLAS detector using 2011 and 2012 LHC proton–proton collision data. Eur. Phys. J. C 74, 3130 (2014). doi:10.1140/epjc/s10052-014-3130-x. arXiv:1407.3935 [hep-ex]
ATLAS Collaboration, Electron and photon energy calibration with the ATLAS detector using LHC run 1 data. Eur. Phys. J. C 74, 3071 (2014). doi:10.1140/epjc/s10052-014-3071-4. arXiv:1407.5063 [hep-ex]
ATLAS Collaboration, Jet energy resolution in proton–proton collisions at \(\sqrt{s} = 7\;\text{ TeV }\) recorded in 2010 with the ATLAS detector. Eur. Phys. J. C 73, 2306 (2013). doi:10.1140/epjc/s10052-013-2306-0. arXiv:1210.6210 [hep-ex]
ATLAS Collaboration, Performance of pile-up mitigation techniques for jets in \(pp\) collisions at \(\sqrt{s}=8\) TeV using the ATLAS detector. Eur. Phys. J. C 76, 581 (2016). doi:10.1140/epjc/s10052-016-4395-z. arXiv: 1510.03823 [hep-ex]
ATLAS Collaboration, Pile-up subtraction and suppression for jets in ATLAS. ATLAS-CONF-2013-083, 2013. http://cds.cern.ch/record/1570994
M. Botje et al., The PDF4LHC Working Group Interim Recommendations (2011). arXiv:1101.0538 [hep-ph]
R.D. Ball et al., Parton distributions with LHC data. Nucl. Phys. B 867, 244 (2013). doi:10.1016/j.nuclphysb.2012.10.003. arXiv:1207.1303 [hep-ph]
ATLAS Collaboration, Luminosity determination in pp collisions at \(\sqrt{s}\) = 8 TeV using the ATLAS detector at the LHC. Eur. Phys. J. C 76, 653 (2016). doi:10.1140/epjc/s10052-016-4466-1. arXiv:1608.03953 [hep-ex]
J.A. Aguilar-Saavedra, J. Carvalho, N.F. Castro, F. Veloso, A. Onofre, Probing anomalous Wtb couplings in top pair decays. Eur. Phys. J. C 50, 519 (2007). doi:10.1140/epjc/s10052-007-0289-4. arXiv:hep-ph/0605190 [hep-ph]
N. Castro, J. Erdmann, C. Grunwald, K. Kröninger, N.-A. Rosien, EFTfitter–a tool for interpreting measurements in the context of effective field theories. Eur. Phys. J. C 76, 432 (2016). doi:10.1140/epjc/s10052-016-4280-9. arXiv:1605.05585 [hep-ex]
ATLAS Collaboration, Search for anomalous couplings in the \(Wtb\) vertex from the measurement of double differential angular decay rates of single top quarks produced in the \(t\)-channel with the ATLAS detector. JHEP 04, 023 (2016). doi:10.1007/JHEP04(2016)023. arXiv:1510.03764 [hep-ex]
CMS Collaboration, Measurement of the \(t\)-channel single-top-quark production cross section and of the \(|V_{tb}|\) CKM matrix element in \(pp\) collisions at \(\sqrt{s} = 8\;\text{ TeV }\). JHEP 06, 090 (2014). doi:10.1007/JHEP06(2014)090. arXiv:1403.7366 [hep-ex]
ATLAS Collaboration, Comprehensive measurements of \(t\)-channel single top-quark production cross sections at \(\sqrt{s} = 7\;\text{ TeV }\) with the ATLAS detector. Phys. Rev. D 90, 112006 (2014). doi:10.1103/PhysRevD.90.112006. arXiv:1406.7844 [hep-ex]
B. Grzadkowski, M. Misiak, Anomalous Wtb coupling effects in the weak radiative B-meson decay. Phys. Rev. D 78, 077501 (2008). [Erratum: Phys. Rev. D 84, 059903 (2011)]. doi:10.1103/PhysRevD.84.059903, doi:10.1103/PhysRevD.78.077501. arXiv:0802.1413 [hep-ph]
ATLAS Collaboration, ATLAS Computing Acknowledgements 2016–2017. ATL-GEN-PUB-2016-002, 2016. http://cds.cern.ch/record/2202407
Acknowledgements
We thank CERN for the very successful operation of the LHC, as well as the support staff from our institutions without whom ATLAS could not be operated efficiently. We acknowledge the support of ANPCyT, Argentina; YerPhI, Armenia; ARC, Australia; BMWFW and FWF, Austria; ANAS, Azerbaijan; SSTC, Belarus; CNPq and FAPESP, Brazil; NSERC, NRC and CFI, Canada; CERN; CONICYT, Chile; CAS, MOST and NSFC, China; COLCIENCIAS, Colombia; MSMT CR, MPO CR and VSC CR, Czech Republic; DNRF and DNSRC, Denmark; IN2P3-CNRS, CEA-DSM/IRFU, France; SRNSF, Georgia; BMBF, HGF, and MPG, Germany; GSRT, Greece; RGC, Hong Kong SAR, China; ISF, I-CORE and Benoziyo Center, Israel; INFN, Italy; MEXT and JSPS, Japan; CNRST, Morocco; NWO, Netherlands; RCN, Norway; MNiSW and NCN, Poland; FCT, Portugal; MNE/IFA, Romania; MES of Russia and NRC KI, Russian Federation; JINR; MESTD, Serbia; MSSR, Slovakia; ARRS and MIZŠ, Slovenia; DST/NRF, South Africa; MINECO, Spain; SRC and Wallenberg Foundation, Sweden; SERI, SNSF and Cantons of Bern and Geneva, Switzerland; MOST, Taiwan; TAEK, Turkey; STFC, United Kingdom; DOE and NSF, United States of America. In addition, individual groups and members have received support from BCKDF, the Canada Council, CANARIE, CRC, Compute Canada, FQRNT, and the Ontario Innovation Trust, Canada; EPLANET, ERC, ERDF, FP7, Horizon 2020 and Marie Skłodowska-Curie Actions, European Union; Investissements d’Avenir Labex and Idex, ANR, Région Auvergne and Fondation Partager le Savoir, France; DFG and AvH Foundation, Germany; Herakleitos, Thales and Aristeia programmes co-financed by EU-ESF and the Greek NSRF; BSF, GIF and Minerva, Israel; BRF, Norway; CERCA Programme Generalitat de Catalunya, Generalitat Valenciana, Spain; the Royal Society and Leverhulme Trust, United Kingdom. The crucial computing support from all WLCG partners is acknowledged gratefully, in particular from CERN, the ATLAS Tier-1 facilities at TRIUMF (Canada), NDGF (Denmark, Norway, Sweden), CC-IN2P3 (France), KIT/GridKA (Germany), INFN-CNAF (Italy), NL-T1 (Netherlands), PIC (Spain), ASGC (Taiwan), RAL (UK) and BNL (USA), the Tier-2 facilities worldwide and large non-WLCG resource providers. Major contributors of computing resources are listed in Ref. [82].
Author information
Authors and Affiliations
Consortia
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Funded by SCOAP3
About this article
Cite this article
Aaboud, M., Aad, G., Abbott, B. et al. Measurement of the W boson polarisation in \(t\bar{t}\) events from pp collisions at \(\sqrt{s}\) = 8 TeV in the lepton + jets channel with ATLAS. Eur. Phys. J. C 77, 264 (2017). https://doi.org/10.1140/epjc/s10052-017-4819-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-017-4819-4