Search for anomalous couplings in the W tb vertex from the measurement of double differential angular decay rates of single top quarks produced in the t-channel with the ATLAS detector

The electroweak production and subsequent decay of single top quarks is determined by the properties of the Wtb vertex. This vertex can be described by the complex parameters of an effective Lagrangian. An analysis of angular distributions of the decay products of single top quarks produced in the t -channel constrains these parameters simultaneously. The analysis described in this paper uses 4.6 fb−1 of proton-proton collision data at s=7\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \sqrt{s}=7 $$\end{document} TeV collected with the ATLAS detector at the LHC. Two parameters are measured simultaneously in this analysis. The fraction f1 of decays containing transversely polarised W bosons is measured to be 0.37 ± 0.07 (stat.⊕syst.). The phase δ− between amplitudes for transversely and longitudinally polarised W bosons recoiling against left-handed b-quarks is measured to be −0.014π ± 0.036π (stat.⊕syst.). The correlation in the measurement of these parameters is 0.15. These values result in two-dimensional limits at the 95% confidence level on the ratio of the complex coupling parameters gR and VL, yielding Re[gR/VL] ∈ [−0.36, 0.10] and Im[gR/VL] ∈ [−0.17, 0.23] with a correlation of 0.11. The results are in good agreement with the predictions of the Standard Model.


Introduction
The top quark is the heaviest known fundamental particle, making the measurement of its production and decay kinematics a unique probe of physical processes beyond the Standard Model (SM).The top quark was first observed in 1995 at the Tevatron [1,2] through its dominant production channel, top-quark pair (t t) production via the flavour-conserving strong interaction.An alternative process produces single top quarks through the weak interaction, first observed at the Tevatron in 2009 [3,4].
The t-channel exchange of a virtual W boson is the dominant process contributing to single top-quark production (see Figure 1).The production cross-section is calculated for proton-proton (pp) collisions at √ s = 7 TeV using a next-to-leading-order (NLO) QCD prediction with resummed next-to-next-to-leading logarithmic (NNLL) accuracy, referred to as approximate next-to-next-to-leading order (NNLO).With a top-quark mass of 172.5 GeV and MSTW2008NNLO [5] parton distribution function (PDF) sets, the cross-section for the t-channel process is calculated to be 64.6 +2.6  −1.7 pb [6].The uncertainties correspond to the sum in quadrature of the error obtained from the MSTW PDF set at the 90% confidence level (C.L.) and the factorisation and renormalisation scale uncertainties.
Figure 1: Representative Feynman diagrams for t-channel single top-quark production and decay.Here q represents a u or d quark, and q represents a d or ū quark, respectively.The initial b-quark arises from (a) a sea b-quark in the 2 → 2 process, or (b) a gluon splitting into a b b pair in the 2 → 3 process.
New physics resulting in corrections to the Wtb vertex would affect t-channel single top-quark production and decay, and can be probed through processes which affect variables sensitive to the angular distributions of the final-state particles from the top-quark decay.Within the effective field theory framework, the Lagrangian can be expressed in full generality as [7,8]: where g is the weak coupling constant, and m W and q ν are the mass and the four-momentum of the W boson, respectively.P L,R ≡ (1 ∓ γ 5 )/2 are the left-and right-handed projection operators and σ µν = i[γ µ , γ ν ]/2.V L,R and g L,R are the complex left-and right-handed vector and tensor couplings, respectively.The couplings could also be expressed using the Wilson coefficients [9] of the relevant dimension-six operators 1 , described in Refs.[7,10].In the SM at tree level, V L = V tb , which is an element in the Cabibbo-Kobayashi-Maskawa (CKM) matrix, while the anomalous couplings V R = g L,R = 0. Deviations from these values would provide hints of physics beyond the SM, and complex values could imply that top-quark decay has a CP-violating component [11].
Indirect constraints on V R , g L , and g R were obtained [12,13] from precision measurements of B-meson decays.These result in very tight constraints on V R and g L , but much looser constraints on g R .Calculations of the anomalous couplings in specific models predict a much larger contribution to g R than to V R or g L [14].However, within the models studied so far for tensor couplings, the magnitude of g R is expected to be < 0.01 [14,15].Measurements of the W boson polarisation fractions [16][17][18][19] are sensitive to the magnitude of combinations of anomalous couplings.In order to obtain constraints on Re g R in those measurements, the anomalous couplings are taken to be purely real, V R = g L = 0 is assumed and V L is fixed to 1.In this analysis, the anomalous couplings V L and g R are allowed to be complex, and the ratios Re g R /V L and Im g R /V L are measured.
This paper is organised as follows.Section 2 describes the coordinate system and parameterisation used in the measurement.Section 3 gives a short description of the ATLAS detector, then Section 4 describes the simulated samples used to predict properties of the t-channel signal and the background processes, and the data samples used.Section 5 describes the definitions of the objects used to identify t-channel events, and the set of criteria applied to reduce the number of background events.The procedures for modelling background processes are described in Section 6, along with comparisons between the predictions and data.Section 7 describes the efficiency, resolution, and background models used to translate the distribution of true t-channel signal events to the distribution of reconstructed signal and background events, and how the parameters of the model are estimated.Section 8 quantifies the sources of uncertainty important in this measurement, Section 9 presents the resulting central value and covariance matrix for the model parameters and the ratios Re g R /V L and Im g R /V L , and the conclusions are given in Section 10.

Measurement definition
An event-specific coordinate system is defined for analysing the decay of the top quark in its rest frame using the directions of the jets, lepton, and W boson in the final state, depicted in Figure 2. The ẑ-axis is chosen along the direction of the W boson boosted into the top-quark rest frame, ẑ ≡ q = q/| q|.The angle θ * between q and the momentum of the lepton in the W boson rest frame, p , is the same angle used to measure the W boson helicity fractions in top-quark decays [17].The spin of single top quarks produced in the t-channel, s t , is predicted to be polarised along the direction of the light quark (q in Figure 1) in the top-quark rest frame [20], called the spectator-quark direction, ps .If this jet defines the spin analysing direction, the degree of polarisation is shown in Refs.[21,22] to be P ≡ ps • s t /| s t | ≈ 0.9 at √ s = 7 TeV for SM couplings.A three-dimensional coordinate system is defined from the q-ps plane and the perpendicular direction, with ŷ = ps × q and x = ŷ × q.In this coordinate system, the direction of the lepton ( ) in the W boson rest frame p is specified by θ * and the complementary azimuthal angle φ * .These angles arise as a natural choice for measuring a normalised double differential distribution, in the helicity basis (top-quark rest frame), of the leptonic decay of the W boson arising from the top quark.In this basis, the t → Wb transition is determined by four helicity amplitudes, where the W boson and b-quark are either both right-handed, both left-handed, or the W boson is longitudinal and the b-quark is either left-or right-handed.The resulting angular distribution of final-state leptons in θ * and φ * is Figure 2: Definition of the coordinate system with x, ŷ, and ẑ defined as shown from the momentum directions of the W boson, q ≡ ẑ, and the spectator jet, ps , in the top-quark rest frame.The angles θ * and φ * indicate the lepton direction p in this coordinate system.expressed [23] as a series in spherical harmonics, Y m l (θ * , φ * ), parameterised by ( The restriction to l ≤ 2 in Equation ( 2) is caused by the limited spin states of the initial and final state fermions and the vector boson at the weak vertex.
The parameters α describing this distribution can be expressed at tree level in terms of four linear combinations of the anomalous couplings.The fractions, f 1 , f + 1 , and f + 0 , depend on the magnitude of these four combinations, while δ + and δ − depend on the phases between pairs of combinations.In addition there is a small correction due to the finite mass of the b-quark, m b .Defining x W = m W /m t and x b = m b /m t , where m t is the top-quark mass, the parameters are given by: • f 1 , the fraction of decays containing transversely polarised W bosons, • f + 1 , the fraction of transversely polarised W boson decays that are right-handed, • f + 0 , in events with longitudinally polarised W bosons, the fraction of b-quarks that are right-handed, • δ + , the phase between amplitudes for longitudinally polarised and transversely polarised W bosons recoiling against right-handed b-quarks, • δ − , the phase between amplitudes for longitudinally polarised and transversely polarised W bosons recoiling against left-handed b-quarks, • P, which is considered separately from α because it depends on the production of the top quark, rather than the decay.There is no analytical expression for P in terms of anomalous couplings, but a parameterisation is determined in Ref. [24] by fitting simulated samples produced with the leading-order (LO) Protos [25] generator 2 with different values for the various couplings.
Through these expressions, measurements of the parameterisation variables, ( α, P), can be used to set limits on the values of the couplings V L,R and g L,R .
The parameters to which this analysis is most sensitive are the fraction f 1 ∈ [0, 1] and the phase δ − ∈ [−π, π].The parameter f 1 can be related to the W boson helicity fractions via . Since f + 1 and f + 0 are small near the SM point, the term in Equation (2) that is proportional to e iδ + is nearly zero.Thus δ + cannot be constrained and f + 0 cannot be separated from P. Constraints on f + 1 are currently better determined from limits on F R [17,18].The variations that the parameters f 1 and δ − induce in the model are demonstrated in Figure 3 for t-channel signal events generated with Protos.For regions of the parameter space ( f 1 , δ − ) close to the SM, only these parameters can be significantly constrained.The dependence of δ − on V R and g L is suppressed by x b , while both f 1 and δ − are dependent on V L and g R directly or through x W . Thus, to simplify the analysis, only variations in V L and g R are considered, while V R and g L are assumed to be zero.With this assumption, the values of f + 1 and f + 0 are very small.The value of P is also determined from the values of V L and g R .The highest-order dependence of f 1 and δ − on the couplings V L and g R appear as the ratio g R V L , where the real and imaginary parts of this ratio are measured separately.This motivates quoting the results in both the parameter space ( f 1 , δ − ) and the coupling space (Re g R /V L , Im g R /V L ).

The ATLAS detector
The ATLAS detector [26] consists of a set of sub-detector systems, cylindrical in the central region and planar in the two end-cap regions, that cover almost the full solid angle around the interaction point 3 .ATLAS is composed of an inner tracking detector (ID) close to the interaction point, surrounded by a superconducting solenoid providing a 2 T axial magnetic field, electromagnetic and hadronic calorimeters, and a muon spectrometer (MS).The ID consists of a silicon pixel detector, a silicon microstrip detector (SCT), providing tracking information within pseudorapidity |η| < 2.5, and a straw-tube transition radiation tracker (TRT) that covers |η| < 2.0.The central electromagnetic calorimeter is a lead and liquid-argon (LAr) sampling calorimeter with high granularity, and is divided into a barrel region that covers |η| < 1.475 and end-cap regions that cover 1.375 < |η| < 3.2.An iron/scintillator tile calorimeter provides hadronic energy measurements in the central pseudorapidity range.The forward regions are 3 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 upward.Cylindrical coordinates (r, φ) are used in the transverse plane, φ being the azimuthal angle around the beam pipe.The pseudorapidity is defined in terms of the polar angle θ as η = − ln tan(θ/2).
instrumented with LAr calorimeters covering 3.1 < |η| < 4.9 for both the electromagnetic and hadronic energy measurements.The MS covers |η| < 2.7 and consists of three large superconducting toroid magnets with eight coils each, a system of trigger chambers, and precision tracking chambers.The ATLAS detector has a three-level trigger system [27], used to select events to be recorded for offline analysis.The first-level trigger is hardware-based and uses a subset of the detector information to reduce the physical event rate from 40 MHz to at most 75 kHz.The second and third levels are software-based and together reduce the event rate to about 300 Hz.

Data and simulated samples
This analysis is performed using pp collision data delivered by the LHC [28] in 2011 at √ s = 7 TeV and recorded by the ATLAS experiment.Stringent detector and data quality requirements are applied, resulting in a data sample corresponding to a total integrated luminosity of 4.59 ± 0.08 fb −1 [29].The events are selected by un-prescaled single-lepton triggers [27,30] that require a minimum transverse energy, E T , of 20 GeV or 22 GeV for electrons and a minimum transverse momentum, p T , of 18 GeV for muons, depending on the data-taking conditions.Samples of events generated using Monte Carlo (MC) simulations are produced for t-channel signal and the background processes, and are used to evaluate models of efficiency and resolution, and to estimate systematic uncertainties.
Samples of simulated t-channel single top-quark events are produced with the AcerMC multi-leg LO generator [31] (version 3.8) using the LO CTEQ6L1 [32] PDF sets.AcerMC incorporates both the 2 → 2 and 2 → 3 processes (see Figure 1) and features an automated procedure to remove the overlap in phase space between them [33].The factorisation and renormalisation scales are set to µ F = µ R = 172.5 GeV.Additional t-channel samples with different anomalous couplings are produced with Protos [25] (version 2.2) using the CTEQ6L1 PDF sets.Events are generated using Protos with the four-flavour scheme 4 , incorporating only the 2 → 3 process, and anomalous couplings are enabled in both the production and the decay vertices, varying Re [V L ] and Im g R simultaneously to keep the top-quark width invariant.The factorisation scale is set to µ 2 F = −p 2 W for the light quark and µ 2 F = p 2 b + m 2 b for the gluon.These samples are used to evaluate t-channel generator modelling uncertainties as described in Section 8.2.They are also used to derive folding models with non-SM values of the couplings, used for performing validation tests in MC simulation and measurements in real data, as described in Section 7. The parton showering (PS), hadronisation and underlying-event (UE) modelling in these samples are simulated with Pythia [34] (version 6.426) using the Perugia 2011C set of tuned parameters (P2011C tune) [35] with the CTEQ6L1 [32] PDF sets.Samples of t t events, s-channel single top-quark events, and associated production of a W boson and top quark (Wt) are produced using the Powheg-box NLO generator (version 1.0) coupled with the NLO CT10 [36] set of PDFs and interfaced with Pythia for showering and hadronisation using the P2011C tune with CTEQ6L1 PDF sets.Additional t t samples are produced with Protos (version 2.2) using the CTEQ6L1 PDF sets.These samples are used to evaluate the background model with non-SM values of the couplings.The PS, hadronisation and UE in these samples are simulated with Pythia (version 6.426) using the AUET2B tune [37] with the MRST LO** [38] PDF sets.All processes involving top quarks are produced assuming m t = 172.5 GeV.
Vector-boson production in association with jets (V+jets) is simulated using the multi-leg LO Alpgen generator [39] (version 2.13) using CTEQ6L1 PDF sets and interfaced to Herwig [40] (version 6.5.20) together with the Jimmy UE model [41] (version 4.31).The contributions of W+light-jets and W+heavyjets (W+bb, W+cc, W+c) are simulated separately.To remove overlaps between the n and n + 1 parton samples the MLM matching scheme [39] is used.Double counting between the inclusive W + n parton samples and samples with associated heavy-quark pair production is removed using an overlap-removal algorithm based on parton-jet ∆R matching [42].Diboson processes (WW, WZ and ZZ) are produced using the Herwig generator with the MRST LO** PDF sets.
After the event-generation step, the samples are passed through the full simulation of the ATLAS detector [43] based on Geant4 [44].Additional Protos samples are passed through the ATLFAST2 simulation [43,45] of the ATLAS detector, which uses a fast simulation for the calorimeters and their response.They are then reconstructed using the same procedure as for data.The simulation includes the effect of multiple pp collisions per bunch crossing (pile-up).The events are weighted such that the distribution of the average number of collisions per bunch crossing is the same as in data.

Object definitions and event selection
The definitions of objects used in the analysis, including electrons, muons, jets and b-tagged jets, and missing transverse momentum, as well as the basic event selection are chosen to be identical to those used for the t-channel cross-section measurements in Ref. [46].This analysis requires exactly one isolated charged light lepton (electron or muon) with transverse momentum p T > 25 GeV and pseudorapidity |η| < 2.5.Exactly one b-tagged jet with |η| < 2.5 and exactly one untagged jet with |η| < 4.5 are required, both with p T > 30 GeV.The second b-quark coming from gluon splitting as shown in Figure 1(b) can result in an additional b-tagged jet.This second b-tagged jet generally has a softer p T spectrum and a broader η distribution compared to the b-tagged jet produced in the top-quark decay.It is often not detected in the experiment and is thus not required in the event selection.To reject jets from pile-up collisions, the jet-vertex fraction JVF , defined as the ratio of the p T of all tracks in the jet originating at the primary vertex to the p T of all tracks associated with the jet, is required to satisfy JVF > 0.75.The magnitude of the missing transverse momentum must be E miss T > 30 GeV.Two additional multijet background rejection criteria are applied.The transverse mass of the lepton-E miss T system, is required to be larger than 30 GeV.Also, a more stringent cut on the lepton p T is applied to events in which the lepton and leading jet, j 1 , are back-to-back, where ∆φ( j 1 , ) is the difference in azimuthal angle between the lepton momentum and the leading jet.
The W boson coming from the decay of the top quark can be reconstructed from the momenta of the lepton and the neutrino by using four-momentum conservation.The missing transverse momentum vector, with magnitude E miss T , is used to represent the transverse component of the neutrino momentum and the longitudinal component p z ν is chosen such that the resulting W boson is on its mass shell.A quadratic expression is found for p z ν .If there are two possible real values, the value closer to zero is taken.If the solutions are complex, the assumption that the neutrino is the only contributor to the E miss T is not valid.Therefore, the E miss T is rescaled, preserving its direction, in order to have physical (real) solutions for p z ν .If two solutions for E miss T are found, the one resulting in the smaller value of |p z ν | is taken.The top-quark momentum is then reconstructed from this W boson momentum and the momentum of the b-tagged jet.Finally, the momenta of the W boson and spectator jet are boosted into the top-quark rest frame to obtain q and p s , used to generate the coordinate system in Figure 2, and the lepton is boosted into the W boson rest frame to obtain p .
In addition to this basic event selection, further discrimination between t-channel single top-quark events and background events is achieved by applying additional criteria: • The pseudorapidity of the untagged jet must satisfy |η| > 2, since the spectator jet tends to be forward in the t-channel signature.
• The sum of the p T of all final-state objects, H T , must be larger than 210 GeV, since the H T distributions of the backgrounds have their peaks just below this value.
• The mass of the top quark reconstructed from its decay products is required to be within 150-190 GeV, to reject background events from processes not involving top quarks.
• The distance in η between the untagged jet and the b-tagged jet must be larger than 1, to further reduce t t contributions.
These criteria and the basic event selection together define the signal region of the analysis.
Two control regions are defined, enhanced in each of the two dominant backgrounds, t t and W+jets, to validate the modelling of the backgrounds by simulated events.In the control region where the t t background is enriched, events with exactly four jets, exactly one of which is b-tagged, passing the basic event selection together with the top-quark mass and H T requirements are selected.In the control region where W+jets backgrounds are enriched, events with exactly two jets passing the basic event selection and an inversion of the top-quark mass criterion are required.

Background estimation and event yields
The largest background contributions to single top-quark t-channel production in the signal region arise from t t production and from W boson production in association with jets (W+jets).Events containing t t production are difficult to distinguish from single top-quark events, as they also have real top quarks in the final state.Also, the production of a W boson with two jets, where one is identified as containing a b-hadron, can mimic the t-channel final state.Multijet production via the strong interaction can enter the signal region as well, and data-driven techniques are required to accurately model it, as explained at the end of this section.Other backgrounds originate from Wt production and s-channel single top-quark, diboson (WW, WZ, and ZZ) and Z+jets production.Most of these background models are derived directly from MC simulation, but specialised procedures are implemented for W+jets and multijet production.The t t cross-section is calculated at NNLO in QCD including resummation of NNLL soft-gluon terms with Top++ 2.0 [47][48][49][50][51][52][53] and is found to be 177 +10 −11 pb.The uncertainties due to the PDFs at 68% C.L. and α S are calculated following the PDF4LHC [54] prescription for the MSTW2008NNLO, CT10, and NNPDF2.3[55] error PDF sets, and are added in quadrature to the scale uncertainty to yield a total uncertainty of 6%.Single top-quark production cross-sections for Wt and s-channel production are calculated at approximate NNLO and are found to be 15.7±1.2pb [56] and 4.63 +0.20  −0.18 pb [57], respectively.The uncertainties correspond to the sum in quadrature of the uncertainty derived from the MSTW2008NNLO error PDF sets at 90% C.L. and the scale uncertainties, yielding a final uncertainty of about 8% for Wt production and about 4% for s-channel production.The Z+jets inclusive production cross-section is estimated with NNLO precision using the FEWZ program [58].Diboson events (WW, WZ, and ZZ) are normalised to the NLO cross-section prediction calculated with MCFM [59].The uncertainty on the combined Z+jets and diboson background is estimated to be 60% as in Ref. [46].
For the production of a W boson in association with jets, the shapes of the relevant kinematic distributions are predicted from the Alpgen sample.They are initially normalised to make the inclusive W cross-section correspond to the NNLO prediction using the FEWZ program [58], with the same scaling factor being applied to the Alpgen prediction for the W+b b, W+cc, and W+light-jets samples.The Alpgen prediction for the W+c process is scaled by a factor that is obtained from a study based on NLO calculations using MCFM [59].Data-driven techniques are then used to estimate the flavour composition and the overall normalisation as described in Ref. [60], and an uncertainty of 18% (15%) for electron (muon) selections is applied.
Multijet background events pass the signal region selection if a jet is mis-identified as an isolated lepton or if the event has a non-prompt lepton from the decay of a hadron that appears to be isolated.Since it is neither possible to simulate a sufficient number of those events nor to calculate the rate precisely, dedicated techniques are used to model multijet events and to estimate their production rate, employing both collision data and simulated events.In the electron channel, mis-identified jets are the main source of multijet background events.This is modelled by the jet-lepton method [46], in which an electron-like jet is selected in Pythia di-jet MC with special requirements and redefined as a lepton.In the muon channel, a data-driven matrix method [46,61] based on the muon impact parameter significance is used.An uncertainty of 50% is assigned to the estimated yield based on comparisons of the rates obtained by using alternative methods [46], i.e. the matrix method in the electron channel and the jet-lepton method in the muon channel, and using an alternative variable, i.e. m T ( , E miss T ) instead of E miss T .Table 1 provides the event yields for the electron and muon channels after the event selection.The predictions are normalised to 4.6 fb −1 and their uncertainties are derived taking into account the statistical uncertainties of the simulated samples.For the multijet background, the uncertainties on the prediction are the statistical uncertainties of the data and simulated samples of the muon and electron channel, respectively.Figure 4 shows the angular distributions cos θ * and φ * after the event selection for the electron and muon channels in data, compared to predictions using SM couplings.The uncertainties are derived taking into account the statistical uncertainties of the data samples for data and simulated samples on the prediction and the 50% systematic uncertainty on the normalisation of the multijet background.Since the number of events in W+light-jets and multijet samples are low, statistical fluctuations and events with large generator weights affect the sample shapes.Therefore, shape templates from events selected without the b-tagging requirement are used for these two backgrounds to smooth the simulated models.In these events, the hardest jet in the event is chosen to take the place of the b-tagged-like jet for the purposes of the event selection and reconstruction.Overall, good agreement between the observed and expected distributions is found.

(d)
Figure 4: Angular distributions cos θ * (upper row) and φ * (lower row) in the signal region for electrons (left) and muons (right), comparing data, shown as the black points with statistical uncertainties, to SM signal and background predictions.The uncertainties shown on the prediction take into account MC sample sizes and the 50% systematic uncertainty on the normalisation of the multijet background.The W+jets backgrounds are normalised to the observed number of events in the control region, and the muon multijet backgrounds are predicted from observed data.All other samples are scaled to their theoretical predictions.The lower plots show the ratio of data to prediction in each bin.
Table 1: Event yields for the electron and muon channels in the signal region.Individual predictions are rounded to integers while "Total expected" corresponds to the rounding of the sum of full precision individual predictions.The uncertainties shown are statistical only.Uncertainties of less than 0.5 events appear as zero.

Analysis method
The model describing the t-channel signal in Section 2 is connected to the angles measured in reconstructed events via an analytic folding procedure.Models of the selection efficiency and detector resolution are derived from simulated t-channel signal events, and a model of the reconstructed background is derived from the sum of the combined background processes described in Section 4. Since the values of the parameters f 1 and δ − depend on top-quark couplings, both the t-channel signal and the t t background are sensitive to their values.Only the angular distributions for the leptonic decay of the W boson are measured, but the efficiency and resolution may also depend on other unmeasured distributions, such as those associated with the η distribution of the top or spectator quark in the laboratory frame.The efficiency, resolution, and background models are constructed such that they account for any dependence on the values of f 1 and δ − , which is required to remove biases dependent on the values of these parameters (see Appendix A).A probability density is derived for all events in the signal region, as a function of the reconstructed angles cos θ * and φ * , and conditional on the parameters.This density is then used to construct a likelihood, from which f 1 and δ − are measured in data.
The signal model described by Equation ( 2) is a series in spherical harmonics, only containing terms up to a maximum value of l, l sig max = 2. Describing the efficiency and resolution models similarly allows the use of the orthogonality properties of spherical harmonics to construct analytically-folded distributions [23].An efficiency function, given by (θ * T , φ * T ; α, P) = describes the probability that a t-channel signal event with the given true angles θ * T and φ * T will be reconstructed in the signal region.The series contains all allowed values of l and m up to a maximum spherical harmonic degree, l eff max .The selected signal density, ρ s , is then defined as the product of the efficiency function and the signal model, ρ, normalised to the overall rate, The Gaunt coefficients G are defined in terms of Clebsch-Gordan coefficients C m 1 m 2 M l 1 l 2 L , given by, The efficiency function is determined from a likelihood fit to simulated t-channel events.A value is chosen for l eff max by comparing the associated values of the corrected Akaike Information Criteria, AICc [62,63], a likelihood-ratio test with an additional penalty term which increases with the number of parameters.From this test, l eff max = 6 is chosen.The deviation in the results from the selection of efficiency models with l eff max varied by ±2 are included as systematic uncertainties.Other criteria, such as the Schwarz Criteria [64], select values of l eff max within the range chosen for this uncertainty.When reconstructing events in the signal region, the finite resolution of the detector results in a migration from the true angles to measured angles θ * and φ * .This migration is modelled by the resolution function, This series contains all terms with allowed values of λ, µ, L , and M up to the maximum degree parameters l true max , associated with the dependence on θ * T and φ * T , and l reco max , associated with the dependence on θ * and φ * .The reconstructed signal density, ρ r , is defined as the convolution of this function with the selected signal density, The resolution function is determined using a spherical Fourier technique [65].In principle, describing a very narrow resolution function could require large values of l true max and l reco max .In practice, the mathematics of angular-momentum addition guarantee that there are no terms in Equation ( 9) with L > l sig max + l eff max .Thus, l true max = l sig max + l eff max .The parameter l reco max is determined using the mean integrated squared error, MISE = (ρ(x) − ρ(x)) 2 dx, where ρ(x) is the true probability density and ρ(x) is a distribution estimating that density.Calculating the MISE with ρ r yields a broad minimum for l reco max > 4, and for this analysis l reco max = 6 is chosen.The deviation in the results from the selection of resolution models with l reco max varied by ±2 are included as systematic uncertainties.This density of reconstructed t-channel signal events is a series in spherical harmonics, with coefficients which are functions of the physics parameters describing the production (P) and decay ( α) of the top quark.A background model, is derived which is also a series in spherical harmonics, and is added to the reconstructed signal density via the signal fraction, f s .This background function is determined from a likelihood fit to all simulated background events in the signal region, and contains all terms with allowed values of λ and µ up to the maximum degree parameter l bkg max .This parameter is determined like the equivalent parameter in the efficiency function, using the AICc.The variations in the results from the selection of background models with l bkg max ± 2 are included as systematic uncertainties.The probability density of all events in the signal region, ρ t , is the sum of the reconstructed signal and background densities with signal fraction f s , Figure 5 shows the data summed over the electron and muon channels for the signal region.Three representative models are also shown: one for the SM expectations, one that is near the expected 95% C.L. sensitivity in f 1 , and one near the expected sensitivity in δ − .These values for f 1 and δ − are the same as those used in Figure 3, but the curves now include the effects of efficiency, resolution and background.The main differences between the two figures are due to the isolation requirements placed on the leptons.For cos θ * = −1, the lepton overlaps the b-tagged jet, while for φ * = π, the lepton may overlap the spectator jet.The region where the efficiency is less than 0.05% (cos θ * < −0.95, and −0.95 < cos θ * < −0.85 for π − 0.1 < φ * < π + 0.1), is excluded from further analysis.For other regions, the efficiency varies, as a function of φ * and cos θ * , from 0.1% up to 1.7% for larger cos θ * and away from the plane containing the spectator jet.As the features of the angular distributions are broad, the effects of resolution are small.The main effect is to smooth out some of the effects of the acceptance and to increase the number of events at large cos θ * where the SM contribution is small.As seen in Figure 4, the background is also somewhat larger for large cos θ * .
An unbinned likelihood function is constructed from the probability density in Equation ( 13) and from a set of events where θ * i and φ * i can be taken either from data or MC simulation, and the event weights w i are 1 for measured data.The likelihood is evaluated on a grid with spacing 0.01 in f 1 and 0.01π in δ − .The    efficiency and resolution models at each point are derived from Protos t-channel simulated events, and the background model uses Protos t t simulated events.Interpreting these parameters as varying the coupling ratio g R /V L , the polarisation P is varied simultaneously according to the parameterisation in Ref. [24].
A central value is measured at the maximum value of L on this grid, determined from a Gaussian fit to the points at which L is evaluated.The 68% and 95% confidence limits are defined by the region where the likelihood ratio, L/L max , is larger than the value that would yield the same likelihood for a twodimensional Gaussian distribution.They are refined to increase the accessible precision by interpolating between points on either side of the contours determined from these evaluated points.The statistical uncertainty is estimated from the symmetrised 68% C.L. interval of each parameter.This model is checked for closure both with MC samples, and by defining toy-MC samples based on Equation (13) with the same number of events as in data.These samples were generated at various points in ( f 1 , δ − ) space.In all cases they are found to reproduce the expected values of f 1 and δ − .These toy-MC samples were also used to derive pull distributions for f 1 and δ − to check the statistical uncertainties returned by the fit.

Sources of systematic uncertainty
Systematic uncertainties are evaluated in the full ( f 1 , δ − ) parameter space.Efficiency, resolution, and background models are determined from MC samples with a parameter varied by its uncertainty, or a subset where appropriate.A likelihood is constructed from the resulting model, using events generated with nominal values of the varied parameters.The difference between the central values estimated at the nominal value of a parameter and at the value varied by its uncertainty, or half the difference between central values estimated with the parameter varied up and down by its uncertainty, is used to construct a covariance matrix for each source of systematic uncertainty.The total covariance matrix for the systematic uncertainties and its correlation matrix are found from the sum of the covariance matrices determined for individual uncertainties.

Object modelling
The uncertainties in the reconstruction of jets, leptons, and E miss T are propagated to the analysis, following the same procedures as described in Ref. [46].The main source of uncertainty from these physics objects is the jet energy scale (JES) [66], which is the largest systematic uncertainty on the measurement of f 1 .
To estimate the impact of the JES uncertainty on the result, the jet energy is scaled up and down by its uncertainty [66], which ranges from 2.5% in the central region with high-p T jets to 14% in the far forward region with low-p T jets.Uncertainties are also estimated for jet energy resolution and reconstruction efficiency; the impact of varying the jet-vertex fraction requirement; E miss T reconstruction and the effect of pile-up collisions on the E miss T ; b-tagging efficiency and mistagging rate; lepton trigger, identification, and reconstruction efficiencies; and lepton momentum, energy scale, and resolution.

MC generators and PDFs
Multiple MC event generators are used to model the t-channel and t t processes in this analysis, and the differences between these generators are included as systematic uncertainties.Comparing the AcerMC generator used for t-channel events and the Powheg-box generator used for t t events to the Protos generator used for both the t-channel single top-quark events and t t events yields the largest uncertainty in δ − .Additional t-channel comparisons are performed between the AcerMC and Powheg-box generators and between AcerMC events with showering by Pythia and Herwig.The renormalisation and factorisation scales in the Powheg-box t-channel sample are also varied independently by factors of 0.5 and 2.0.Additional t t comparisons are performed between the Powheg-box generator and the MC@NLO [67] generator with showering by Herwig, and between Powheg-box events with showering by Pythia and Herwig.The variations between PDF sets and within individual sets are used to estimate a systematic uncertainty following the PDF4LHC prescription.

Signal and background normalisation
Cross-sections for background processes given in Section 6 are varied up and down by their uncertainties.The multijet normalisation is varied by 50% as discussed in Section 6, and the remaining backgrounds are varied simultaneously to produce a conservative estimate.The impact of object modelling and other uncertainties on the W+jets background normalisation are considered in parallel with the variations made in t-channel and t t samples.A separate shape uncertainty is assigned to the W+jets samples by varying the matching and factorisation scales in the Alpgen generator.The integrated luminosity is varied up and down by its uncertainty, ±1.8%, derived as detailed in Ref. [29].

Detector correction and background parameterisation
An uncertainty due to the limited size of the MC samples used to estimate the efficiency and resolution models is estimated from the statistical uncertainties derived from a measurement of f 1 and δ − in the t-channel signal MC sample alone.The background statistical uncertainty is estimated by varying the background model according to the eigenvectors of its covariance matrix.The background statistical uncertainty dominates the total "MC statistics" uncertainty listed in Table 2. Its significance reflects the small size of some background samples in the signal region, and the resulting disparate values of the sample weights.The effect of varying the cutoff degree l max in the determination of these models is also estimated, and is found to be small compared to MC statistical uncertainty.

Uncertainty combination
Table 2 shows the contribution of each source of uncertainty to the measurement of the parameters f 1 and δ − and their correlation, ρ( f 1 , δ − ).The total systematic uncertainty and correlation is obtained from the sum of the covariance matrices determined for each source.It is combined with the covariance matrix of the statistical uncertainty via the following method.At each point α i in the ( f 1 , δ − ) space, a multivariate normal distribution N i is constructed with the covariance matrix representing the systematic uncertainties, Σ syst , and the mean, α i .The resulting distribution is evaluated at a point α j and multiplied by the likelihood at this point, L stat j .The maximum modified likelihood value, over all possible points α j , is kept, and the resulting broadened likelihood distribution L stat+syst i is used to represent the measurement with both statistical and systematic variation incorporated: The resulting confidence limits are taken as the full uncertainty in the measurement.

Results
The result for ( f 1 , δ − ) and the coupling ratios (Re g R /V L , Im g R /V L ) is shown in Figure 6.The 68% contour represents the total uncertainty on the measurement.
The parameters f 1 and δ − and their uncertainties are measured to be f 1 = 0.37 ± 0.05 (stat.)± 0.05 (syst.),δ − = − 0.014π ± 0.023π (stat.)± 0.028π (syst.).( 16) The correlation in the measurement of these parameters is ρ ( f 1 , δ − ) = 0.15.The results are compatible with the SM expectations at LO, derived from expressions in Refs.[11,68] with m t = 172.5 GeV, m W = 80.399 GeV, and m b = 4.95 GeV: f 1 = 0.304 and δ − = 0.The dependence of the parameters f 1 and δ − on the top-quark mass is evaluated using t-channel and t t simulation samples with a range of different top-quark masses.A linear dependence is found, resulting from changes in acceptance at different masses, with a slope of −0.019 GeV −1 for f 1 and a negligible slope for δ − .The uncertainty due to the top-quark mass dependence is not included in the total systematic uncertainty since it has no significant impact on the results.

The propagation of the uncertainties to the (Re
The correlation in the measurement of these coupling ratios is ρ Re g R /V L , Im g R /V L = 0.11.The effect on the propagation due to the current uncertainty in the top-quark, W boson and b-quark masses [69] is < 0.01 in Re g R /V L , and < 0.0001 in Im g R /V L . Limits are placed simultaneously on the possible complex values of the ratio of the anomalous couplings   13) and data events into (a) f 1 , (b) δ − , (c) f 1 vs. δ − , and (d) Re g R /V L vs. Im g R /V L , with systematic uncertainties incorporated.The black points indicate the largest evaluated likelihood in each bin of the projected variable.Gaussian fits to the one-dimensional projections were performed, displayed as the red curve.Regions shown in green and yellow represent the 68% and 95% confidence level regions, respectively.A black line or cross indicates the observed value, and the grey line or point indicates the SM expectation.
g R and V L at 95% C.L.,

Re
The best constraints on Re g R come from W boson helicity fractions in top-quark decays, with Re g R of [−0.08, 0.04] and [−0.08, 0.07], both at 95% C.L., from ATLAS [17] and from CMS [18], respectively.However, these limits use the measured single top-quark production cross-section [46,70] along with the assumption that V L = 1 and Im g R = 0. Without these assumptions no value within the range 0.0 Re g R /V L 0.8 can be excluded.The limits presented in this paper remove these assumptions and extend the knowledge of g R to the whole complex plane by simultaneously measuring information about Re g R /V L and Im g R /V L ; the latter is previously unmeasured.

Conclusion
The analysis described in this publication is a two-dimensional measurement of anomalous couplings in the Wtb vertex in t-channel single top-quarks events, performed on 4.6 fb −1 of pp collisions at a centre-of-mass energy of 7 TeV collected by the ATLAS detector at the LHC.An analytic folding model and likelihood maximisation techniques are used to extract values of the parameters f 1 and δ − in a parameterisation of the decay of the top quark, which translate to limits on the couplings Re g R /V L and Im g R /V L with V R = g L = 0.The coupling parameterisation in terms of the spherical angles θ * and φ * defines the angular distribution of the t-channel signal process.Efficiency and resolution functions are fit in this space, which are then folded into the signal model describing the underlying physics analytically.A background function is fit and added to the signal model.The full model is used to construct a likelihood and its characterisation provides estimators of the best-fit central value, uncertainties, and correlations between f 1 and δ − , or Re g R /V L and Im g R /V L .The result is combined with estimated sources of systematic uncertainty to give the measurement of the parameters, f 1 = 0.37 ± 0.05 (stat.)± 0.05 (syst.),δ − = −0.014π± 0.023π (stat.)± 0.028π (syst.), and ρ( f 1 , δ − ) = 0.15, or Re g R /V L = −0.13± 0.07 (stat.)± 0.10 (syst.),Im g R /V L = 0.03 ± 0.06 (stat.)± 0.07 (syst.), and ρ(Re g R /V L , Im g R /V L ) = 0.11 in terms of the couplings.The measurement sets limits at the 95% confidence level on Re g R /V L ∈ [−0.36, 0.10] and Im g R /V L ∈ [−0.17, 0.23].The limits on Re g R /V L are complementary to limits from B-meson decay and W helicity measurements, and contain previously unmeasured information about Im g R /V L .The measured values are in agreement with SM expectations.

Appendix
A. Parameter dependence of the folding and background models The derivation of the efficiency, resolution, and background models, described in Equations (8)(9)(10)(11)(12)(13), is based on the form of these distributions in MC simulation.For the background model, constructed from the sum of all predicted backgrounds with an appreciable effect on the distribution, this includes events containing top quarks, primarily from t t production, the distribution of which is affected by changing the values of the couplings V L and g R .The efficiency and resolution models are averages over all unmeasured distributions in the signal.Variations in the values of V L and g R alter those unmeasured distributions, which could lead to a dependence on these couplings for the efficiency and resolution models.For instance, t-channel single top-quark production depends on anomalous coupling in both the top-quark production and decay vertices, so varying the couplings alters production-side distributions, such as the p T and η distributions of the top or spectator quark.
To test how these potential sources of dependence affect the models, closure tests are performed with input MC events generated at different values of the parameters f 1 and δ − using efficiency, resolution, and background models derived from MC events generated with SM parameter values.The t-channel signal and t t background events are generated with Protos for both the input MC events and the events used to derive the models.The results of these closure tests can then be compared to the values with which the input events were generated; any trend connecting them other than a line with a slope of one and an intercept of zero indicates a dependence which could bias the measurement.The results of these tests, shown in Figure 7, indicate biases in the extracted parameters.The dashed line with a slope of one and intercept of zero, representing no bias in the measurement, is added to guide the eye.A non-negligible bias is evident, and is represented at the lowest reasonable order with a polynomial fit.For δ − a linear fit is sufficient, while for f 1 a quadratic fit gives a better description.
If the biases shown are due to the mismatch between the values of f 1 and δ − in the efficiency, resolution, and background models, they can be removed by varying these parameters in the functions describing these models to match their values in the signal model.This is achieved by deriving these models from tchannel signal MC events representing each point in ( f 1 , δ − ) independently, and calculating the likelihood using the same f 1 and δ − in the signal model and in the t-channel signal MC events.The results of these tests, shown in Figure 8, indicate that this procedure eliminates the biases shown in Figure 7.This procedure for deriving the efficiency, resolution, and background models is used in the measurement.The dashed line with a slope of one and intercept of zero, representing no bias in the measurement, is added to guide the eye.All estimates are consistent with this line; linear fits produce a slope consistent with one and an intercept consistent with zero.

Figure 3 :
Figure 3: Projections into (a) φ * in bins of cos θ * , (b) cos θ * , and (c) φ * in Equation (2), for different values of f 1 and δ − .The black points represent the Protos t-channel signal generated with SM parameters, and the curves shown represent the signal model.For the three curves shown, the parameters f 1 and δ − are set to their values in the SM, f 1 = 0.3, δ − = 0 (solid red), and to two sets of beyond-the-SM values, f 1 = 0.1, δ − = 0 (dashed blue) and f 1 = 0.3, δ − = 0.1π (dotted green).

Figure 6 :
Figure6: Projections of the likelihood function constructed from the signal region probability density Equation(13) and data events into (a) f 1 , (b) δ − , (c) f 1 vs. δ − , and (d) Re g R /V L vs. Im g R /V L , with systematic uncertainties incorporated.The black points indicate the largest evaluated likelihood in each bin of the projected variable.Gaussian fits to the one-dimensional projections were performed, displayed as the red curve.Regions shown in green and yellow represent the 68% and 95% confidence level regions, respectively.A black line or cross indicates the observed value, and the grey line or point indicates the SM expectation.

Figure 7 :
Figure 7: Biases in the estimation of the values of the parameters (a) f 1 and (b) δ − which arise when the efficiency, resolution, and background models are determined from SM simulated events, where f 1 = 0.3 and δ − = 0.The points represent the central values obtained from the likelihood fit, with the expected statistical uncertainty in 4.6 fb −1 of data at √ s = 7 TeV shown by the error bars.The dashed line with a slope of one and intercept of zero, representing no bias in the measurement, is added to guide the eye.A non-negligible bias is evident, and is represented at the lowest reasonable order with a polynomial fit.For δ − a linear fit is sufficient, while for f 1 a quadratic fit gives a better description.

Figure 8 :
Figure 8: Comparison of the estimation of the values of the parameters (a) f 1 and (b) δ − derived by allowing their values in the efficiency, resolution, and background models to vary together with the values in the signal model.The points represent the central values obtained from the likelihood fit, with the expected statistical uncertainty in 4.6 fb −1 of data at √ s = 7 TeV shown by the error bars.The dashed line with a slope of one and intercept of zero, representing no bias in the measurement, is added to guide the eye.All estimates are consistent with this line; linear fits produce a slope consistent with one and an intercept consistent with zero. 1

Table 2 :
Sources of systematic uncertainty on the measurement of f 1 and δ − using AcerMC t-channel single topquark simulated events and backgrounds estimated from both MC simulation and data, including Powheg-box t t simulation.Individual sources are evaluated separately for shifts up and down, and symmetrised uncertainties σ( f 1 ), σ(δ − ) and correlation coefficients ρ( f 1 , δ − ) are given.