arxiv: v1 [hep-ex] 23 Mar 2018

Transcription

1 EUROPEAN ORGANIZAION FOR NUCLEAR RESEARCH (CERN CERN-EP /03/6 -OP7-00 arxiv: v [ep-ex 3 Mar 08 Measurement of differential cross sections for te production of top quark pairs and of additional jets in lepton+jets events from pp collisions at s = 3 ev e Collaboration Abstract Differential and double-differential cross sections for te production of top quark pairs in proton-proton collisions at s = 3 ev are measured as a function of kinematic variables of te top quarks and te top quark-antiquark (tt system. In addition, kinematic variables and multiplicities of jets associated wit te tt production are measured. is analysis is based on data collected by te experiment at te LHC in 06 corresponding to an integrated luminosity of 35.8 fb. e measurements are performed in te lepton+jets decay cannels wit a single muon or electron and jets in te final state. e differential cross sections are presented at te, witin a pase space close to te experimental acceptance, and at te parton level in te full pase space. e results are compared to several standard model predictions tat use different metods and approximations. e kinematic variables of te top quarks and te tt system are reasonably described in general, toug none predict all te measured distributions. In particular, te transverse momentum distribution of te top quarks is more steeply falling tan predicted. e kinematic distributions and multiplicities of jets are adequately modeled by certain combinations of next-to-leading-order calculations and parton sower models. Submitted to Pysical Review D c 08 CERN for te benefit of te Collaboration. CC-BY-4.0 license See Appendix E for te list of collaboration members

2

3 Introduction Measurements of differential production cross sections of top quark pairs (tt provide important information for testing te standard model and searcing for penomena beyond te standard model. Precise teoretical predictions of tese measurements are callenging since iger-order effects of quantum cromodynamics (QCD and electroweak (EW corrections [ are important. Moreover, te generation of tt events requires a realistic modeling of te parton sower (PS. e measured kinematic properties and multiplicities of jets allow for a detailed comparison of different PS models to te data and provide insigt into teir tuning. In tis paper, differential and double-differential production cross sections as a function of kinematic variables of te top quarks and te tt system are reported. In addition, measurements of multiplicities and kinematic properties of jets in tt events are presented. e measurements are based on proton-proton (pp collision data at a center-of-mass energy of 3 ev corresponding to an integrated luminosity of 35.8 fb [. e data were recorded by te experiment at te CERN LHC in 06. Only tt decays into te l+jets (l = e, µ final state are considered, were, after te decay of eac top quark into a bottom quark and a W boson, one of te W bosons decays adronically and te oter one leptonically. Hence, te experimental signature consists of two jets coming from te adronization of b quarks (b jets, two jets from a adronically decaying W boson, a transverse momentum imbalance associated wit te neutrino from te leptonically decaying W boson, and a single isolated muon or electron. is measurement continues a series of differential tt production cross section measurements in pp collisions at te LHC. Previous measurements of differential cross sections at s = 7 ev [3, 4 and 8 ev [5 ave been performed in various tt decay cannels. First measurements at 3 ev are available [ 4. Previous studies of multiplicities and kinematic properties of jets in tt events can be found in Refs. [5 7. Wit about 5 times more data and an improved understanding of systematic uncertainties, we provide an update and extension to te previous analysis in te l+jets cannel at 3 ev [8. We measure differential cross sections defined in two ways: at te and te parton level. For te particle-level measurement a proxy of te top quark is defined based on experimentally accessible quantities, suc as properties of jets, wic are made up of quasistable particles wit a mean lifetime greater tan 30 ps. ese quantities are described by teoretical predictions tat require modeling of te PS and adronization, in addition to te matrix-element calculations. e kinematic requirements on tese objects are cosen to closely reproduce te experimental acceptance. Muons and electrons stemming from τ lepton decays are not treated separately and can contribute to te particle-level signal. A detailed definition of particle-level objects is given in Section 3. e particle-level approac as te advantage tat it reduces teoretical uncertainties in te experimental results by avoiding teory-based extrapolations from te experimentally accessible portion of te pase space to te full range, and from jets to partons. For te parton-level measurement top quarks in te l+jets decay cannel are defined as signal directly before teir decays into a bottom quark and a W boson. e τ+jets decay cannel is not considered ere as signal even in cases were te τ lepton decays into a muon or electron. No restriction on te pase space is applied for parton-level top quarks. e corrections and extrapolations used in tis measurement are based on a next-to-leading-order (NLO calculation of tt production, combined wit a simulation of te PS. For bot particle- and parton-level measurements te tt system is reconstructed at te detector level wit a likeliood-based approac using te top quark and W boson mass constraints to

4 identify te corresponding top quark decay products. e differential cross sections are measured at te particle and parton levels as a function of te transverse momentum p and te absolute rapidity y of te top quarks, separately for te adronically (labeled t and leptonically (labeled t l decaying W bosons, and te p, y, and invariant mass M of te tt system. In addition, te differential cross sections at te parton level are determined as a function of te lower- and iger-p values of te top quarks in an event. Double-differential cross sections for te following combinations of variables are determined at bot levels: y(t vs. p (t, M(tt vs. y(tt, and p (t vs. M(tt. At, te differential cross sections as a function of p (t, p (tt, and M(tt are measured in bins of jet multiplicity. Using te four jets identified as te tt decay products and te four igest-p additional jets, te cross sections are determined as a function of te jet p and absolute pseudorapidity η, te minimal separation R jt of jets from anoter jet in te tt system, and te separation R t of jets from te closest top quark. Here R = ( φ + ( η, were φ and η are te differences in azimutal angle (in radians and pseudorapidity between te directions of te two objects. Finally, we determine te gap fraction, defined as te fraction of events tat do not contain jets above a given p tresold, and te jet multiplicities for various tresolds of te jet p. is paper is organized as follows: In Section, we provide a description of te signal and background simulations, followed by te definition of te particle-level top quarks in Section 3. After a sort overview of te detector and te particle reconstruction in Section 4, we describe te object and event selections in Section 5. Section 6 contains a detailed description of te reconstruction of te tt system. Details on te background estimation and te unfolding are presented in Sections 7 and 8. After a discussion of systematic uncertainties in Section 9, te differential cross sections as a function of observables of te top quark and te tt system are presented in Section. Finally, observables involving jets are discussed in Section. e results are summarized in Section. Signal and background modeling e Monte Carlo generator POWHEG [9 (v,vq is used to calculate te production of tt events at NLO accuracy in QCD. e renormalization µ r and factorization µ f scales are set to te transverse mass m = m t + p of te top quark, were a top quark mass m t = 7.5 GeV is used in all simulations. e result is combined wit te PS simulations of PYHIA8 [3, 4 (v8.9 using te underlying event tune CUEP8M4 [5, 6, and of HERWIG++ [7 (v.7. using te tune EE5C [8. In addition, MADGRAPH5 amc@nlo [9 (v.. (MG5 amc@nlo is used to produce a simulation of tt events wit additional partons. All processes wit up to two additional partons are calculated at NLO and combined wit te PYHIA8 PS simulation using te FxFx [30 algoritm. e scales are selected as µ r = µ f = (m (t + m ( t. e default parametrization of te parton distribution functions (PDFs used in all simulations is NNPDF30 nlo as 08 [3. e simulations are normalized to an inclusive tt production cross section of pb [3. is value is calculated wit next-to-nlo (NNLO accuracy, including te resummation of next-to-next-to-leading-logaritmic soft-gluon terms. Its uncertainty is evaluated by varying te coice of µ r and µ f and by propagating uncertainties in te PDFs. Distributions tat correspond to variations in te PDFs or te scales µ r and µ f are obtained by applying different event weigts. ese distributions are used for te corresponding uncertainty estimates. For additional uncertainty estimations we use POWHEG+PYHIA8 simulations wit top quark masses of 7.5 and 73.5 GeV, wit initial and final PS scales varied up and down by a factor of two, wit variations of te underlying event tune, and a simulation wit an alternative color-reconnection model.

5 3 e main backgrounds are simulated using te same tecniques. e MG5 amc@nlo generator is used for te simulation of W boson production in association wit jets, t-cannel single top quark production, and Drell Yan (DY production in association wit jets. e generator POWHEG [33 is used for te simulation of single top quark associated production wit a W boson (tw, and PYHIA8 is used for multijet production. In all cases, te PS and te adronization are described by PYHIA8. e W boson and DY backgrounds are normalized to teir NNLO cross sections calculated wit FEWZ [34 (v3.. e t-cannel single top quark production is normalized to te NLO calculation obtained from HAHOR [35 (v.. e production of tw is normalized to te NLO calculation [36, and te multijet simulation is normalized to te LO calculation obtained wit PYHIA8 [4. e detector response is simulated using GEAN4 [37. e simulations include multiple pp interactions per bunc crossing (pileup. e simulated events are weigted, depending on teir number of pileup interactions, to reproduce te measured pileup distribution. Finally, te same reconstruction algoritms tat are applied to te data are used for te simulated events. 3 Particle-level top quark definition e definitions of particle-level objects constructed from quasi-stable simulated particles, obtained from te predictions of tt event generators before any detector simulation, are summarized below. ese particle-level objects are furter used to define te particle-level top quarks. Detailed studies on particle-level definitions can be found in Ref. [38. All simulated muons and electrons are corrected for effects of bremsstralung by adding te poton momenta to te momentum of te closest lepton if teir separation is R < 0.. All potons are considered for te momentum correction. A corrected lepton is selected if it fulfills te isolation requirement tat te p sum of all quasi-stable particles, excluding corrected leptons and neutrinos, witin R = 0.4 is less tan 35% of te corrected lepton p. In addition, we require te lepton to ave p > 5 GeV and η <.4. Simulated potons wit p > 5 GeV and η <.4 tat are not used in te momentum correction of a lepton are considered if teir isolation, defined analogously to te lepton isolation, is below 5%. All neutrinos are selected including tose stemming from decays of adrons. Jets are clustered by te anti-k jet algoritm [39, 40 wit a distance parameter of 0.4. All quasi-stable particles wit te exception of neutrinos are clustered. Jets wit p > 5 GeV and η <.4 are selected if tere is no isolated lepton or poton, as defined above, witin R = 0.4. b jets at te are defined as tose jets tat contain a b adron. As a result of te sort lifetime of b adrons, tese are not quasi-stable particles and only teir decay products sould be considered for te jet clustering. However, to allow teir association wit a jet, te b adrons are also included wit teir momenta scaled down to a negligible value. is preserves te information of teir directions, but removes teir impact on te jet clustering. Based on te invariant masses of tese objects, we construct a pair of particle-level top quarks in te l+jets final state. Events wit exactly one muon or electron wit p > 30 GeV and η <.4 are selected. Simulated events wit an additional muon or electron wit p > 5 GeV and η <.4 are rejected. We take te sum of te four-momenta of all neutrinos as te neutrino

6 4 candidate momentum p ν from te leptonically decaying top quark and find te permutation of jets tat minimizes te quantity [M(p ν + p l + p bl m t + [M(p jw + p jw m W + [M(p jw + p jw + p b m t, ( were p jw, are te four-momenta of two ligt-flavor jet candidates, considered as te decay products of te adronically decaying W boson; p bl, are te four-momenta of two b jet candidates; p l is te four-momentum of te lepton; and m W = 80.4 GeV [4 is te mass of te W boson. All jets wit p > 5 GeV and η <.4 are considered. At least four jets are required, of wic at least two must be b jets. e remaining jets wit p > 30 GeV and η <.4 are defined as additional jets. Events wit a adronically and a leptonically decaying particle-level top quark are not required to be l+jets events at te parton level, e.g., tt dilepton events wit additional jets can be identified as l+jets event at te if one lepton fails to pass te selection. As an example, te comparison between te p (t distributions at te particle and parton levels are sown in Fig. and demonstrates te direct relation between particle-level and parton-level top quarks. o obtain an unambiguous nomenclature for te jets, we define j W to be te jet in te W boson decay wit te iger p. e additional jets j i are sorted by teir transverse momenta were j as te igest p. Fraction Simulation Purity Stability Bin efficiency Non-parton-level signal Particle level p (t [GeV (t Bin at : p Simulation Bin at parton level: p (t Figure : Comparison between te p (t distributions at te particle and parton level, extracted from te POWHEG+PYHIA8 simulation. Left: fraction of parton-level top quarks in te same p bin at te (purity, fraction of particle-level top quarks in te same p bin at te parton level (stability, ratio of te number of particle- to parton-level top quarks (bin efficiency, and fraction of events wit a particle-level top quark pair tat are not considered as signal events at te parton level (non-parton-level signal. Rigt: p -bin migrations between particle and parton level. e p range of te bins can be taken from te left panel. Eac column is normalized suc tat te sum of its entries corresponds to te fraction of particle-level events in tis bin at te parton level in te full pase space Normalized at parton level per bin 4 e detector e central feature of te detector is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8. Witin te solenoid volume are a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL, and a brass and scintillator adron calorimeter (HCAL, eac composed of a barrel and two endcap sections. Forward calorimeters extend te η coverage provided by te barrel and endcap detectors. Muons

7 5 are measured in gas-ionization detectors embedded in te steel flux-return yoke outside te solenoid. A more detailed description of te detector, togeter wit a definition of te coordinate system and relevant kinematic variables, can be found in Ref. [4. e particle-flow (PF event algoritm [43 reconstructs and identifies eac individual particle wit an optimized combination of information from te various elements of te detector. e energy of muons is obtained from te curvature of te corresponding track. e energy of electrons is determined from a combination of te electron momentum at te primary interaction vertex as determined by te tracker, te energy of te corresponding ECAL cluster, and te energy sum of all bremsstralung potons spatially compatible wit originating from te electron track. e energy of potons is directly obtained from te ECAL measurement, corrected for zero-suppression effects. e energy of carged adrons is determined from a combination of teir momentum measured in te tracker and te matcing ECAL and HCAL energy deposits, corrected for zero-suppression effects and for te response function of te calorimeters to adronic sowers. Finally, te energy of neutral adrons is obtained from te corresponding corrected ECAL and HCAL energy. 5 Pysics object reconstruction and event selection e measurements presented in tis paper depend on te reconstruction and identification of muons, electrons, jets, and missing transverse momentum associated wit a neutrino. Muons and electrons are selected if tey are compatible wit originating from te primary vertex, wic, among te reconstructed primary vertices, is te one wit te largest value of summed pysics-object p. e pysics objects are jets, clustered using te jet finding algoritm [39, 40 wit te tracks assigned to te primary vertex as inputs, and te associated missing transverse momentum, taken as te negative vector sum of te p of tose jets. Since leptons from tt decays are typically isolated, a requirement on te lepton isolation is used to reject leptons produced in decays of adrons. e lepton isolation variables are defined as te sum of te p of neutral adrons, carged adrons, and poton PF candidates witin a cone of R = 0.4 for muons and R = 0.3 for electrons. It is required to be less tan 5% (6% of te muon (electron p. Event-by-event corrections are applied to maintain a pileup-independent isolation efficiency. e muon and electron reconstruction and selection efficiencies are measured in te data using tag-and-probe tecniques [ Depending on te p and η, teir product is 75 85% for muons and 50 80% for electrons. Jets are clustered from PF objects using te anti-k jet algoritm wit a distance parameter of 0.4 implemented in te FASJE package [40. Carged particles originating from a pileup interaction vertex are excluded. e total energy of te jets is corrected for energy depositions from pileup. In addition, p - and η-dependent corrections are applied to correct for te detector response effects [47. If an isolated lepton wit p > 5 GeV witin R = 0.4 around a jet exists, te jet is assumed to represent te isolated lepton and is removed from furter consideration. For te identification of b jets, te combined secondary vertex algoritm [48 is used. It provides a discriminant between b and non-b jets based on te combined information of secondary vertices and te impact parameter of tracks at te primary vertex. A jet is identified as a b jet if te associated value of te discriminant exceeds a tresold criterion wit an efficiency of about 63% and a combined carm and ligt-flavor jet rejection probability of 97%. e missing transverse momentum p miss is calculated as te negative of te vectorial sum of transverse momenta of all PF candidates in te event. Jet energy corrections are also propa-

8 6 gated to improve te measurement of p miss. Events considered for tis analysis are selected by single-lepton triggers. ese require p > 4 GeV for muons and p > 7 GeV for electrons, as well as various quality and isolation criteria. o reduce te background contributions and optimize te tt reconstruction, additional requirements are imposed on te recorded events. Events wit exactly one muon or electron wit p > 30 GeV and η <.4 are selected. No additional muons or electrons wit p > 5 GeV and η <.4 are allowed. In addition to te lepton, at least four jets wit p > 30 GeV and η <.4 are required. At least two of tese jets must be identified as b jets. We compare several kinematic distributions in te data to te simulation separately for te muon and electron cannels to verify tat tere are no unexpected differences. e ratios of te measured to te expected event yields in te two cannels agree witin te uncertainty in te lepton reconstruction and selection efficiencies. In te remaining steps of te analysis, te two cannels are combined by adding teir distributions. 6 Reconstruction of te top quark-antiquark system e reconstruction of te tt system follows closely te metods used in Ref. [8. e goal is te correct identification of detector-level objects as parton- or particle-level top quark decay products. In te simulation, a jet or lepton at te can be spatially matced to te corresponding detector-level object. If no one-to-one assignment to a corresponding detector-level object is possible for any of te objects in te particle-level tt system, te event is considered as nonreconstructable in te particle-level measurement. For te parton-level measurement a quark from te tt decay is assigned to te detector-level jet wit te igest p witin R = 0.4 around te parton. If no one-to-one correspondence at detector level is found for any of tese quarks or te leptons, te event is nonreconstructable in te parton-level measurement. In particular, tis includes events wit merged topologies were at least two quarks are matced to te same jet. Based on tese relations between detector level and parton or, te efficiencies of te tt reconstruction are studied. A detailed discussion on te relationsip between quantities at te parton or and detector level is presented in Section 8. For te reconstruction all possible permutations of assigning detector-level jets to te corresponding tt decay products are tested and a likeliood tat a certain permutation is correct is evaluated. Permutations are considered only if te two jets wit te igest b identification probabilities are te two b jet candidates. In eac event, te permutation wit te igest likeliood is selected. e likelioods are evaluated separately for te particle- and te parton-level measurements. For eac tested permutation te neutrino four-momentum p ν is reconstructed using te algoritm of Ref. [49. e idea is to find all possible solutions for te tree components of te neutrino momentum vector using te two mass constraints (p ν + p l = m W and (p ν + p l + p bl = m t. Eac equation describes an ellipsoid in te tree-dimensional momentum space of te neutrino. e intersection of tese two ellipsoids is usually an ellipse. We select p ν as te point on te ellipse for wic te distance D ν,min between te ellipse projection onto te transverse plane and p miss is minimal. is algoritm leads to a unique solution for te longitudinal neutrino momentum and an improved resolution of its transverse component. For te cases were te invariant mass of te lepton and b l candidate is above m t no solution can be found and te corresponding permutation is discarded. e minimum distance D ν,min is also used to identify

9 7 te correct b l, as described below. e value of D ν,min from te neutrino reconstruction and te mass constraints on te adronically decaying top quark are combined in a likeliood function λ, given by log[λ = log[p m (m, m 3 log[p ν (D ν,min, ( were P m is te two-dimensional probability density of te invariant masses of W bosons and top quarks tat are correctly reconstructed, based on te matcing criteria described above. e value of λ is maximized to select te permutation of jets. e probability density P m is calculated as a function of te invariant mass of te two jets, m, tested as te W boson decay products, and te invariant mass of te tree jets, m 3, tested as te decay products of te adronically decaying top quark. e distributions for te correct jet assignments, taken from te POWHEG+PYHIA8 simulation and normalized to unit area, are sown in Fig. (upper for te particle- and parton-level measurements. is part of te likeliood function is sensitive to te correct reconstruction of te adronically decaying top quark. Permutations wit probabilities less tan 0.% of te maximum value of te probability density P m are rejected. is selection criterion discards less tan % of te correctly reconstructed events. Especially in te parton-level measurement, it removes events tat are incompatible wit te ypotesis of a adronically decaying top quark and reduces te background contribution. is is caused by te stringent mass constraints for a parton-level top quark, were, in contrast to te particlelevel top quark, close compatibility wit te top quark and W boson masses are required. e probability density P ν describes te distribution of D ν,min for a correctly selected b l. In Fig. (lower, te normalized distributions of D ν,min for b l and for oter jets are sown. On average, te distance D ν,min for a correctly selected b l is smaller and as a smaller tail compared to te distance obtained for oter jets. Permutations wit values of D ν,min > 50 GeV are rejected since tey are very unlikely to originate from a correct b l association. is part of te likeliood function is sensitive to te correct reconstruction of te leptonically decaying top quark. Since te likeliood function λ combines te probabilities from te reconstruction of te leptonically and adronically decaying top quarks, it provides information on te reconstruction of te wole tt system. e performance of te reconstruction as a function of jet multiplicity is sown for several tt simulations in Fig. 3, were we use te input distributions P m and P ν from te POWHEG+PYHIA8 simulation. e reconstruction efficiency of te algoritm is defined as te probability tat te most-likely permutation, as identified troug te maximization of te likeliood λ, is te correct one, given tat all tt decay products are reconstructed and selected. e performance deteriorates wit te increase in te number of jets, since te number of permutations increases drastically and te probability of selecting a wrong permutation increases. e differences observed between te various simulations are taken into account in te estimation of te systematic uncertainties. We observe a lower reconstruction efficiency for te particle-level measurement. is is caused by te less powerful mass constraints for a particle-level top quark. is can be seen in te mass distributions of Fig. and te likeliood distributions in Fig. 4, were te simulations are normalized to te measured integrated luminosity of te data sample, and te tt simulation is divided into te following categories: correctly reconstructed tt systems (tt rigt reco; events were all decay products are available, but te algoritm failed to identify te correct permutation (tt wrong reco; te nonreconstructable events (tt nonreconstructable; and events tat are according to te parton- or particle-level definitions not tt signal events (tt nonsignal. Only te last category is treated as tt background, wile te oter categories are considered as signal. e lower reconstruction efficiency of te particle-level top quark is compensated by te iger number of reconstructable events.

10 8 M(W [GeV Simulation parton level 4 M(W [GeV Simulation M(t [GeV M(t [GeV Simulation parton level Correct b l Simulation Correct b l 0.03 Incorrect b l 0.03 Incorrect b l D ν,min [GeV D ν,min [GeV Figure : Upper: normalized two-dimensional mass distribution of te correctly reconstructed adronically decaying W bosons M(W and te correctly reconstructed top quarks M(t for te (left parton- and te (rigt particle-level measurements. Lower: normalized distributions of te distance D ν,min for correctly and incorrectly selected b jets from te leptonically decaying top quarks. e distributions are taken from te POWHEG+PYHIA8 tt simulation. reconstruction eff. Simulation parton level FSR down FSR up MG5_aMC@NLO P8 reconstruction eff. Simulation FSR down FSR up MG5_aMC@NLO P8 tt tt Additional jets Additional jets Figure 3: Reconstruction efficiency of te tt system as a function of te number of additional jets for te (left parton- and (rigt particle-level measurements. e efficiencies are calculated based on te simulations wit POWHEG+PYHIA8 (P8 wit scale variations up and down of te final-state PS, POWHEG+HERWIG++ (H++, and MG5 amc@nlo+pyhia8. e vertical bars represent te statistical uncertainties in eac simulation.

11 9 Events / 0.4 Pred parton level 35.8 fb (3 ev tt rigt reco tt wrong reco tt nonreconstructable tt nonsignal Single t Multijet, DY/W+jets Exp. uncertainty log(λ Events / 0.4 Pred fb (3 ev tt rigt reco tt wrong reco tt nonreconstructable tt nonsignal Single t Multijet, DY/W+jets Exp. uncertainty log(λ Figure 4: Distribution of te negative log-likeliood for te selected best permutation in te (left parton- and te (rigt particle-level measurements in data and simulations. Events generated wit POWHEG+PYHIA8 are used to describe te tt production. e contribution of multijet, DY, and W boson plus jets background events is extracted from te data (cf. Section 7. Combined experimental (cf. Section 9 and statistical uncertainties (atced area are sown for te total predicted yields. e data points are sown wit statistical uncertainties. e ratios of data to te sum of te predicted yields are provided at te bottom of eac panel. In Fig. 5, te p of te jets from te tt system, as identified by te reconstruction algoritm, and of te additional jets are presented and compared to te simulation. In Fig. 6, te distributions of p and y of te reconstructed top quarks, and in Fig. 7, te distributions of p (tt, y(tt, and M(tt for te parton- and particle-level measurements are sown. e simulations are normalized according to te measured integrated luminosity of te data. In general, good agreement is observed between te data and te simulation, altoug all measured p spectra are softer tan predicted by te simulation. 7 Background subtraction After te event selection and tt reconstruction, we observe about and events at te parton and s, respectively, were total background contributions of 4.5 and 6.0% from single top quark, DY, W boson, and multijet events are predicted. ese backgrounds ave to be estimated and subtracted from te selected data. In addition, a residual contamination from nonsignal tt events is expected and estimated from te simulation, as detailed below. e predictions of te single top quark background are taken from simulations. Its overall contribution corresponds to about.7 and 3.3% of te selected data in te parton- and particlelevel measurements, respectively. Single top quark production cross sections are calculated wit a precision of a few percent [35, 36. However, tese calculations do not consider te production of additional jets as required by te tt selection. erefore, we use an overall uncertainty of 50%, wic represents a conservative estimate of te PS modeling, scale, and PDF uncertainties. Even wit suc a conservative estimate, its impact on te precision of te final results is negligible. After te full tt selection, te numbers of events in te simulations of multijet, DY, and W boson production are not sufficient to obtain smoot background distributions. erefore, we extract a common distribution for tese backgrounds from a control region in te data. Its selection differs from te signal selection by te requirement of aving no b-tagged jet in te event. In tis control region, te contribution of tt events is estimated to be about 5%, wile te remaining

12 Events / GeV fb (3 ev tt signal tt nonsignal Single t Multijet, DY/W+jets Exp. uncertainty Events / GeV fb (3 ev tt signal tt nonsignal Single t Multijet, DY/W+jets Exp. uncertainty 3 3 Pred. Events / GeV p (b [GeV l 35.8 fb (3 ev tt signal tt nonsignal Single t Multijet, DY/W+jets Exp. uncertainty Pred. Events / GeV p (b (b [GeV [GeV 35.8 fb (3 ev tt signal tt nonsignal Single t Multijet, DY/W+jets Exp. uncertainty 3 3 Pred. Events / GeV p (j (j [GeV [GeV 35.8 fb W W (3 ev tt signal tt nonsignal Single t Multijet, DY/W+jets Exp. uncertainty Pred. Events / GeV p (j (j [GeV [GeV 35.8 fb W W (3 ev tt signal tt nonsignal Single t Multijet, DY/W+jets Exp. uncertainty Pred. Events / GeV p (j (j [GeV [GeV 35.8 fb (3 ev tt signal tt nonsignal Single t Multijet, DY/W+jets Exp. uncertainty Pred. Events / GeV p (j (j [GeV [GeV 35.8 fb (3 ev tt signal tt nonsignal Single t Multijet, DY/W+jets Exp. uncertainty Pred p (j (j [GeV [GeV 3 Pred p (j (j [GeV [GeV Figure 5: Comparisons between data and simulation at te of te reconstructed distributions of te p of jets as identified by te tt reconstruction algoritm. e simulation of POWHEG+PYHIA8 is used to describe te tt production. e contribution of multijet, DY, and W boson plus jets background events is extracted from te data (cf. Section 7. Combined experimental (cf. Section 9 and statistical uncertainties (atced area are sown for te total predicted yields. e data points are sown wit statistical uncertainties. e ratios of data to te predicted yields are given at te bottom of eac panel. 4

13 Events / 40 GeV parton level 35.8 fb (3 ev tt signal tt nonsignal Single t Multijet, DY/W+jets Exp. uncertainty Events / 40 GeV fb (3 ev tt signal tt nonsignal Single t Multijet, DY/W+jets Exp. uncertainty 3 3 Pred. Events / 0.5 Pred p (t [GeV l parton level 35.8 fb (3 ev tt signal tt nonsignal Single t Multijet, DY/W+jets Exp. uncertainty y(t y(t l Pred. Events / 0.5 Pred p (t [GeV fb (3 ev tt signal tt nonsignal Single t Multijet, DY/W+jets Exp. uncertainty y(t y(t Figure 6: Comparisons of te reconstructed p (upper and y (lower in data and simulations for te t l (left at te parton level and te t (rigt at te. e simulation of POWHEG+PYHIA8 is used to describe te tt production. e contribution of multijet, DY, and W boson plus jets background events is extracted from te data (cf. Section 7. Combined experimental (cf. Section 9 and statistical uncertainties (atced area are sown for te total predicted yields. e data points are sown wit statistical uncertainties. e ratios of data to te predicted yields are given at te bottom of eac panel.

14 Events / 0 GeV parton level 35.8 fb (3 ev tt signal tt nonsignal Single t Multijet, DY/W+jets Exp. uncertainty Events / 0 GeV fb (3 ev tt signal tt nonsignal Single t Multijet, DY/W+jets Exp. uncertainty 3 3 Pred. Events / 0.5 Pred. Events / 80 GeV p (tt [GeV fb (3 ev parton level tt signal 0 tt nonsignal Single t 0 Multijet, DY/W+jets 80 Exp. uncertainty y(tt y(tt 35.8 fb (3 ev 6 parton level tt signal tt nonsignal Single t Multijet, DY/W+jets Exp. uncertainty Pred. Events / 0.5 Pred. Events / 80 GeV p (tt (tt [GeV [GeV fb tt signal 60 tt nonsignal 40 Single t 0 Multijet, DY/W+jets 0 Exp. uncertainty y(tt y(tt 35.8 fb (3 ev 6 tt signal tt nonsignal (3 ev Single t Multijet, DY/W+jets Exp. uncertainty Pred M(tt M(tt [GeV [GeV Pred M(tt M(tt [GeV [GeV Figure 7: Comparisons of te reconstructed distributions of p (tt (upper, y(tt (middle, and M(tt (lower for te (left parton- and te (rigt particle-level measurements in data and simulation. e simulation of POWHEG+PYHIA8 is used to describe te tt production. e contribution of multijet, DY, and W boson plus jets background events is extracted from te data (cf. Section 7. Combined experimental (cf. Section 9 and statistical uncertainties (atced area are sown for te total predicted yields. e data points are sown wit statistical uncertainties. e ratios of data to te predicted yields are given at te bottom of eac panel.

15 3 fraction consists of multijet, DY, and W boson events. e background distributions are obtained after applying exactly te same tt reconstruction algoritm as in te signal region. e two b jet candidates still ave te igest value of te b identification discriminant to maintain a similar number of allowed permutations of jets in te control and signal regions. o estimate te sape dependency on te selection of te control region, we vary te selection tresold of te b identification discriminant. is canges te tt signal contribution and te flavor composition. However, we find te observed sape variations to be negligible. For te background subtraction te distributions extracted from te control region are normalized individually in eac bin of jet multiplicity to te yield of multijet, DY, and W boson events predicted by te simulation in te signal region. In te control region, te expected and measured event yields agree witin teir statistical uncertainties. aking into account te statistical uncertainty in te normalization factor and te sape differences between te signal and control regions in te simulation, we derive an overall uncertainty of 0% in tis background estimation. Special care as to be taken wit te contribution of nonsignal tt events. For te parton-level measurement tese are dilepton, all-jets, and τ+jets events. For te particle-level measurement all tt events for wic no pair of particle-level top quarks exists are considered as nonsignal tt events. e corresponding contributions are about.5% for bot te parton- and te particlelevel measurements. e beavior of tese backgrounds depends on te tt cross section, and a subtraction according to te expected value can result in a bias of te measurement, especially if large differences between te simulation and te data are observed. However, te sapes of te distributions from data and simulation are consistent witin teir uncertainties, and we subtract te predicted relative fractions from te remaining event yields. 8 Corrections to particle and parton levels After te subtraction of te backgrounds, an unfolding procedure is used to correct te reconstructed distributions for detector-specific effects, e.g., efficiency and resolutions, and to extrapolate eiter to te parton or. We do not subtract te fractions of wrongly reconstructed or nonreconstructable events, since in many of tese events a rater soft jet is misidentified, wic as little impact on te resolution of te measured quantities. e iterative D Agostini metod [50 is used to unfold te data. e migration matrices, wic relate te quantities at te parton or and at detector level, and te acceptances are needed as te input. However, not only te detector simulation, but also te teoretical description of tt events affects te migration matrix. is dependence is reduced in te particle-level measurement, were no extrapolation to parton-level top quarks is needed. For te central results te migration matrices and te acceptances are taken from te POWHEG+PYHIA8 simulation, and oter simulations are used to estimate te uncertainties. e binning of te distributions is optimized based on te resolution in te simulation. e minimal bin widts are selected suc tat, according to te resolution, at least 50% of te events are reconstructed in te correct bin. As an example, te migration matrices for te parton- and particle-level measurements of p (t are sown in te rigt-and plots of Fig. 8. For te measured parton-level distributions of any quantity we define te purity as te fraction of parton-level top quarks in te same bin at te detector level, te stability as te fraction of detector-level top quarks in te same bin at te parton level, and te bin efficiency as te ratio of te number of detector- to parton-level top quarks in te same bin. Similar parameters are defined for te particle-level distributions. e purity, stability, and bin efficiency are sown for te p (t measurements in te left-and plots of Fig. 8. ese illustrate te improved agreement between te reconstructed and te unfolded quantities, as well as te reduced extrapolation in te particle-level measurement.

16 4 Fraction Fraction Simulation parton level Purity Stability Bin efficiency Detector level p (t [GeV Simulation Purity Stability Bin efficiency Detector level p (t [GeV at detector level (t Bins of p Simulation, parton level Bins of p (t at parton level Simulation, Bins of p (t at Figure 8: Migration studies of te (upper parton- and (lower particle-level measurements of p (t, extracted from te POWHEG+PYHIA8 simulation. Left: purity, stability, and bin efficiency. Rigt: bin migrations between detector and parton (particle level. e p range of te bins can be taken from te left panels. Eac column is normalized suc tat te sum of its entries corresponds to te percentage of reconstructed events in tis bin at te parton (particle level. at detector level (t Bins of p Normalized per unfolded bin [% Normalized per unfolded bin [%

17 5 o control te level of regularization, te iterative D Agostini metod takes te number of iterations as an input parameter. e initial distributions for te D Agostini unfolding are taken from te POWHEG+PYHIA8 simulation. e number of iterations is cosen suc tat te compatibility between a model and te unfolded data at eiter te parton or is te same as te compatibility between te folded model and te data at detector level. e compatibilities are determined by χ tests at eac level tat are based on all te available simulations and on several modified spectra obtained by reweigting te POWHEG+PYHIA8 distributions of p (t, y(t, p (tt, or p (j before te detector simulation. e modified spectra are cosen suc tat te effect of te reweigting corresponds rougly to te observed differences between te data and te unmodified simulation at detector level. We ave found tat te number of iterations needed to fulfill te above criterion is suc tat a second χ test between te detector-level spectrum wit its statistical uncertainty and te refolded spectrum wit zero uncertainty exceeds a probability of 99.9%. e refolded spectrum is obtained by inverting te unfolding step. is consists of a multiplication wit te response matrix and does not need any regularization. e algoritm needs between 4 and 56 iterations depending on te distribution. e numbers of iterations are iger for measurements wit lower purities and stabilities of te migration matrices. is is te case for te measurements of p (t l and y(t l, wose resolutions are significant worse tan tose of p (t and y(t due to te missing neutrino information. For te two-dimensional measurements wit n bins in one quantity and m i, i =... n bins in te oter te D Agostini unfolding can be generalized using a vector of B = n i m i entries of te form: b,, b,... b n,,... b,m, b,m... b n,mn, wit a corresponding B B migration matrix. e number of iterations is optimized in te same way. In te measurements of jet kinematic properties, we do not unfold te measured spectra of eac jet separately, but do correct for te effect of misidentified jets. e response matrix sowing te migration among te identified jets is given in Fig. 9 for te measurements of te jet p spectra. Fraction Simulation, Purity Stability Bin efficiency p (b p (b p (j p (j p (j p (j p (j p (j l W W 3 4 Bins of jet type at detector level Bins of jet type at detector level p (j p (j p (j W p (j W p (b p (b l Simulation p (b l p (b p (j p (j W W p (j, p (j p (j p (j 3 4 Bins of jet type at Figure 9: Migration studies of te particle-level measurement of te jet p spectra, extracted from te POWHEG+PYHIA8 simulation. Left: purity, stability, and bin efficiency. Rigt: bin migrations between detector and. On te axes te p bins for eac jet are sown. Eac column is normalized in te way tat te sum of its entries corresponds to te percentage of reconstructed events in tis bin at te. p (j 4 p (j Normalized per unfolded bin [%

18 6 9 Systematic uncertainties Several sources of experimental and teoretical systematic uncertainty are considered. Uncertainties in te jet and p miss calibrations, pileup modeling, b identification and lepton selection efficiencies, and integrated luminosity fall into te first category. e total uncertainty in te jet energy calibration is te combination of 9 different sources of uncertainty and te jet-flavor-specific uncertainties [47, were te uncertainty for b jets is evaluated separately. For eac uncertainty source te energies of jets in te simulation are sifted up and down. At te same time, p miss is recalculated accordingly to te rescaled jet energies. e recomputed backgrounds, response matrices, and acceptances are used to unfold te data. e observed differences between tese and te original results are taken as an uncertainty in te unfolded event yields. e same tecnique is used to calculate te impact of te uncertainties in te jet energy resolution, te uncertainty in p miss not related to te jet energy calibration, te b identification, te pileup modeling, and te lepton reconstruction and selection. e b identification efficiency in te simulation is corrected using scale factors determined from data [5. ese ave an uncertainty of about 3% depending on te p of te b jet. e effect on te measurement due to te uncertainty in te modeling of pileup in simulation is estimated by varying te average number of pileup events per bunc crossing by 4.6% [5 and reweigting te simulated events accordingly. e trigger, reconstruction, and identification efficiencies of leptons are evaluated wit tagand-probe tecniques using Z boson dilepton decays [45, 46. e uncertainties in te scale factors, wic are used to correct te simulation to matc te data, take into account te different lepton selection efficiencies in events wit ig jet multiplicities as in tt events. e uncertainty in te lepton reconstruction and selection efficiencies depends on p and η and is below % in te relevant pase-space region. e relative uncertainty in te integrated luminosity measurement is.5% [. Uncertainties in te coice of µ r and µ f, te combination of te matrix-element calculation wit te PS, te modeling of te PS and adronization, te top quark mass, and te PDFs fall into te second category of uncertainties. e effects of tese teoretical uncertainties are estimated eiter by using te event weigts introduced in Section, or by using a tt signal simulation wit varied settings. Again, te uncertainties are assessed using te recomputed backgrounds, response matrices, and acceptances to unfold te data. e scales µ r and µ f are varied up and down by a factor of two individually and simultaneously in te same directions. Afterwards, te envelope of te observed variations is quoted as te uncertainty. e uncertainty in te combination of te matrix-element calculation wit te PS is estimated from an 40% variation of te damp parameter in POWHEG, normally set to damp =.58m t. is variation as been found to be compatible wit te modeling of jet multiplicities in previous measurements at s = 8 ev [6. o estimate te uncertainty in te PS, several effects ave been studied and are assessed individually. e scale of te initial- (ISR and final-state (FSR radiation is varied up and down by a factor of and, respectively. ese variations are motivated by te uncertainties in te PS tuning [5. e effect of multiple parton interactions and te parametrization of color reconnection ave been studied in Ref. [6 and are varied accordingly. In addition, we use a simulation wit activated color reconnection of resonant decays. is enables te color recon-