Measurement of the differential cross sections for top quark pair production as a function of kinematic event variables in pp collisions at sqrt(s) = 7 and 8 TeV

Measurements are reported of the normalized differential cross sections for top quark pair production with respect to four kinematic event variables: the missing transverse energy; the scalar sum of the jet transverse momentum (pT); the scalar sum of the pT of all objects in the event; and the pT of leptonically decaying W bosons from top quark decays. The data sample, collected using the CMS detector at the LHC, consists of 5.0 inverse femtobarns of proton-proton collisions at sqrt(s) = 7 TeV and 19.7 inverse femtobarns at sqrt(s) = 8 TeV. Top quark pair events containing one electron or muon are selected. The results are presented after correcting for detector effects to allow direct comparison with theoretical predictions. No significant deviations from the predictions of several standard model event simulation generators are observed.


Introduction
The CERN LHC produced millions of top quark pairs (tt) in 2011 and 2012.This allows for a detailed investigation of the kinematic event properties of tt production such as the missing transverse energy (E miss T ), the scalar sum of the jet transverse momenta (H T ), the scalar sum of the transverse momenta of all objects (S T ), and the transverse momentum (p W T ) of leptonically decaying W bosons produced in top quark decays.These measurements can be used to verify current theoretical models, along with their implementation in simulations of tt production, and also to measure rare standard model (SM) processes such as tt production in association with a W, Z, or Higgs boson.Since top quark pair production is a major background for many searches for physics beyond the SM, it is important that the properties of tt events are well understood.
Here, we report measurements carried out using the CMS detector [1] at the LHC at two different proton-proton center-of-mass energies.The data samples used include integrated luminosities of 5.0 fb −1 collected in 2011 at √ s = 7 TeV and 19.7 fb −1 from 2012 at √ s = 8 TeV.The tt production cross section is measured as a function of E miss T , H T , S T , and p W T , corrected for detector effects, and compared with the predictions from different event generators.Differential tt cross sections have previously been measured at the Tevatron [2, 3], and at the LHC [4-9].These previous measurements study the tt production cross section as a function of the top quark kinematics and the kinematics of the tt system.The results presented here are complementary, since the tt production cross section is measured as a function of variables that do not require the reconstruction of the top quarks from their decay products.
Top quarks decay with close to 100% probability into a W boson and a bottom quark.In this article, we consider the channel in which one of the W bosons decays leptonically into a charged lepton (electron or muon) along with its associated neutrino, while the other W boson decays hadronically.This channel has a branching fraction of around 15% for direct decay to each lepton flavor and a relatively clean experimental signature, including an isolated, high-transversemomentum lepton, large E miss T from the undetected neutrino, and multiple hadronic jets.Two jets are expected to contain b hadrons from the hadronization of the b quarks produced directly in the t → bW decay, while other jets (from the hadronic W boson decay or gluon radiation) will typically contain only light and charm quarks.

The CMS detector
The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter, and a brass and scintillator hadron calorimeter, each composed of a barrel and two endcap sections.Muons are measured in gasionisation detectors embedded in the steel flux-return yoke outside the solenoid.Extensive forward calorimetry complements the coverage provided by the barrel and endcap detectors.
A more detailed description of the CMS detector, together with a definition of the coordinate system used and the relevant kinematic variables, can be found in Ref. [1].

Simulation
For the Monte Carlo (MC) simulation of the tt signal sample the leading-order MADGRAPH v5.1.5.11 event generator [10] is used with relevant matrix elements for up to three addi-tional partons implemented.Theoretical production cross section values of 177.3 +4.6  −6.0 (scale) ± 9.0 (PDF+α S ) pb at √ s = 7 TeV, and 252.9 +6.4 −8.6 (scale) ± 11.7 (PDF+α S ) pb at √ s = 8 TeV, are used for the normalization of these samples.These cross sections are calculated with the Top++2.0program to next-to-next-to-leading order (NNLO) in perturbative QCD, including soft-gluon resummation to next-to-next-to-leading-logarithm (NNLL) order [11], and assuming a top quark mass m t = 172.5 GeV.The first uncertainty comes from the independent variation of the renormalization (µ R ) and factorization (µ F ) scales, while the second one is associated with variations in the parton distribution function (PDF) and α S , following the PDF4LHC prescription with the MSTW2008 68% CL NNLO, CT10 NNLO, and NNPDF2.35f FFN PDF sets [12][13][14][15][16].
The generated events are subsequently processed with PYTHIA v6.426 [17] for parton showering and hadronization.The PYTHIA parton shower is matched to the jets from the hard quantum chromodynamics (QCD) matrix element via the MLM prescription [18] with a transverse momentum (p T ) threshold of 20 GeV.The CMS detector response is simulated using GEANT4 [19].
The most significant backgrounds to tt production are events in which a W boson is produced in association with additional jets.Other backgrounds include single top quark production, Z boson production in association with multiple jets, and QCD multijet events where hadronic activity is misidentified as a lepton.The simulation of background from W and Z boson production in association with jets is also performed using the combination of MADGRAPH and PYTHIA, with a p T matching threshold of 10 GeV in this case.These samples are referred to as W+jets and Z+jets, respectively.Single top quark production via tand s-channel W boson exchange [25] and with an associated on-shell W boson [26] are generated using POWHEG.The QCD multijet processes are simulated using PYTHIA.The event yields of the background processes are normalized according to their predicted production cross section values.These are from NNLO calculations for W+jets and Z+jets events [27,28], next-to-leading order calculations with NNLL corrections for single top quark events [29], and leading-order calculations for QCD multijet events [17].
Samples are generated using the CTEQ6L PDFs [30] for MADGRAPH samples, the CT10 PDFs [31] for POWHEG samples, and the CTEQ6M PDFs [30] for MC@NLO.The PYTHIA Z2 tune is used to describe the underlying event in both the MADGRAPH and POWHEG + PYTHIA samples at √ s = 7 TeV, whereas the Z2* tune is used for the corresponding samples at √ s = 8 TeV [32].The underlying event in the POWHEG + HERWIG samples is described by the AUET2 tune [33], whereas the default tune is used in the MC@NLO + HERWIG sample.
The value of the top quark mass is fixed to m t = 172.5 GeV in all samples.In all cases, PYTHIA is used for simulating the gluon radiation and fragmentation, following the prescriptions of Ref. [34].Additional simulated hadronic pp interactions ("pileup"), in the same or nearby beam crossings, are overlaid on each simulated event to match the high-luminosity conditions in actual data taking.
Previous measurements of differential tt production cross sections at the LHC [4, 5, 8] showed that several of the tt event generators considered in this analysis predict a harder top quark p T spectrum than that observed in data.An additional simulated tt sample is considered here, where the sample produced with the MADGRAPH event generator is reweighted to improve the agreement of the top quark p T spectrum with data.

Event reconstruction and selection
Parallel selection paths for the two lepton types are implemented, resulting in samples classified as electron+jets and muon+jets.The trigger for the electron+jets channel during the √ s = 7 TeV data taking selects events containing an electron candidate with p T > 25 GeV and at least three reconstructed hadronic jets with p T > 30 GeV.In the √ s = 8 TeV data, at least one electron candidate with p T > 27 GeV is required, with no additional requirement for jets.In the muon+jets channel, at least one isolated muon candidate with p T > 24 GeV is required at the trigger level.Each candidate event is required to contain at least one well-measured vertex [35], located within the pp luminous region in the center of CMS.
Events are reconstructed using a particle-flow (PF) technique [36,37], which combines information from all subdetectors to optimize the reconstruction and identification of individual long-lived particles.
Electron candidates are selected with a multivariate technique using calorimetry and tracking information [38].Inputs to the discriminant include information about the calorimeter shower shape, track quality, track-shower matching, and a possible photon conversion veto.Electron candidates are required to have E T > 30 GeV and pseudorapidity in the range |η| < 2.5.The low-efficiency region 1.44 < |η| < 1.57 between the barrel and endcap sections of the detector is excluded.Muon candidates are selected with tight requirements on track and vertex quality, and on hit multiplicity in the tracker and muon detectors [39].These requirements suppress cosmic rays, misidentified muons, and nonprompt muons from decay of hadrons in flight.Muon candidates are required to have p T > 26 GeV and |η| < 2.1.
For the lepton isolation requirement, a cone of size ∆R = √ (∆η) 2 + (∆φ) 2 is constructed around the lepton direction, where ∆η and ∆φ are the differences in pseudorapidity and azimuthal angle (in radians), respectively, between the directions of the lepton and another particle.The p T values of charged and neutral particles found in this cone are summed, excluding the lepton itself and correcting for the effects of pileup [38].The relative isolation variable I(∆R) is defined as the ratio of this sum to the lepton p T .Lepton candidates are selected if they satisfy I(0.3) < 0.1 for electrons, and I(0.4) < 0.12 for muons.
Reconstructed particles are clustered into jets using the anti-k T algorithm [40] with a distance parameter of 0.5.The measured p T of each jet is corrected [41] for known variations in the jet energy response as a function of the measured jet η and p T .The jet energy is also corrected for the extra energy deposition from pileup interactions [42,43].Jets are required to pass loose identification requirements to remove calorimeter noise [44].Any such jet whose direction is less than ∆R = 0.3 from the identified lepton direction is removed.For the identification of b quark jets ("b tagging"), a "combined secondary vertex" algorithm [45] is used, taking into account the reconstructed secondary vertices and track-based lifetime information.The b tagging threshold is chosen to give an acceptance of 1% for light-quark and gluon jets with a tagging efficiency of 65% for b quark jets.
The final selection requires exactly one high-p T , isolated electron or muon.Events are vetoed if they contain an additional lepton candidate satisfying either of the following criteria: an electron with p T > 20 GeV, |η| < 2.5, and I(0.3) < 0.15; or a muon, with looser requirements on hit multiplicity, and with p T > 10 GeV, |η| < 2.5, and I(0.4) < 0.2.The event must have at least four jets with p T > 30 GeV, of which at least two are tagged as containing b hadrons.
After the final selection, 26 290 data events are found at √ s = 7 TeV, and 153 223 at √ s = 8 TeV.The tt contribution to these event samples, as estimated from simulation, is about 92%.The fraction of true signal events in the samples is 78%.Misidentified all-hadronic or dileptonic tt events, and events containing tau leptons among the tt decay products, comprise 14% of the samples.The remaining events are approximately 4% single top quark events, 2% W/Z+jets events, and 2% QCD multijet events.The efficiency for signal events to satisfy the final selection criteria is about 8%, as determined from simulation.

Cross section measurements
We study the normalized tt differential production cross section as a function of four kinematic event variables: E miss T , H T , S T , and p W T .The variable E miss T is the magnitude of the missing transverse momentum vector p miss T , which is defined as the projection on the plane perpendicular to the beams of the negative vector sum of the momenta of all PF candidates in the event: , where p i x and p i y are the x and y momentum components of the ith candidate, and the sums extend over all PF candidates.The measured E miss T is corrected for pileup and nonuniformities in response as a function of φ [46].
The variable H T is defined as the scalar sum of the transverse momenta of all jets in the event, where the sum extends over all jets having p T > 20 GeV and |η| < 2.5.
The variable S T is the scalar sum of H T , E miss T , and the p T of the identified lepton, T is the magnitude of the transverse momentum of the leptonically decaying W boson, which is derived from the momentum of the isolated lepton and p miss    The points are the data histograms, with the vertical bars showing the statistical uncertainty, and the predictions from the simulation are the solid histograms.The shaded region shows the uncertainty in the values from simulation.These include contributions from the statistical uncertainty and the uncertainty in the tt cross section.The lower plots show the ratio of the number of events from data and the prediction from the MC simulation.
For simulated tt signal events, these four kinematic variables are also calculated using the momenta of particles in the event, before the simulation of the detector response.We refer to the quantities calculated in this way as the generated variables.The generated value of E miss T is the magnitude of the vector sum of the p T of all neutrinos in the event.The long-lived particles in the event are clustered into jets in the same way as the reconstructed particles.The generated value of H T is the sum of the magnitudes of the p T of these jets with p T > 20 GeV and |η| < 2.5.The generated values of S T and p W T are calculated in the same way as the corresponding reconstructed variables, using the p T of the charged lepton from the leptonic decay of a W boson coming from t → bW decay.
The choice of bin widths for this measurement is optimized separately for each kinematic event variable to minimize the migration between bins.This optimization is based on three criteria: (i) of the simulated signal events for which the value of the generated variable falls in the bin, at least 50% are required to have the reconstructed variable in the same bin (this is sensitive to migration of events out of the bin); (ii) of the simulated signal events for which the value of the reconstructed variable falls in the bin, at least 50% are required to have the generated variable in the same bin (this is sensitive to migration of events into the bin); (iii) the number of reconstructed simulation events in a bin is required to be more than 100.These criteria ensure that bin-to-bin migrations are kept small, while allowing a differential cross section measurement with reasonable granularity.
The number of tt events in each bin of each kinematic event variable, and in each channel, is obtained by subtracting the expected contributions of background processes from data.The contributions of single top quark, and W or Z boson plus jet events are estimated from simulation.
In the case of the QCD multijet background, the contribution is estimated from data using a control region where the selection criteria are modified to enrich the contribution of QCD multijet events.In the electron+jets channel, the control region is obtained by inverting the photon conversion veto on the electron.In addition to this, the number of b-tagged jets is required to be exactly zero.The small contamination of tt, single top, W+jets, and Z+jets events in this control region, as estimated from simulation, is subtracted from the data.Then, the ratio of simulated QCD multijet events in the control region and the signal region is used to scale the normalization of the data-driven QCD multijet estimate from the control region to the signal region in the data.The control region in the muon+jets channel is obtained by inverting the isolation criterion on the muon in the selected events, and by requiring exactly zero b-tagged jets.The jet selection criterion is also modified, requiring at least three jets.The same procedure is then followed to estimate the contribution of QCD mulitjet events in the muon+jets signal region.
The number of tt events from data in each bin is then corrected for the small fractions of dileptonic, all-hadronic, and tau tt events in the final sample, as determined from simulation, and for experimental effects, such as detector resolution, acceptance, and efficiency.This correction is performed by constructing a response matrix that maps the generated values to the reconstructed values for the four kinematic variables in the simulated tt signal events.The response matrix is constructed using the MADGRAPH tt sample.This matrix is then inverted, using regularized singular-value decomposition [47] in the ROOUNFOLD [48] software framework.Since we impose no requirements on the generated events, the procedure corrects to the full signal phase space.
The fully-corrected numbers of tt events in the electron+jets and muon+jets channels yield consistent results.These are then added and used to calculate the normalized tt differential production cross section with respect to each kinematic event variable, X, using: where x j represents the number of unfolded signal events in bin j, ∆ X j is the width of bin j; σ is the total tt production cross section, and N = ∑ i x i is the total number of unfolded signal events.

Systematic uncertainties
The systematic uncertainties in the experimental and theoretical input quantities are evaluated and propagated to the final results, taking correlations into account.Since the final result is normalized to the total number of events, the effect of uncertainties that are correlated across all bins is negligible.As such, only uncertainties that affect the shape of the measured distributions are significant.
The uncertainty coming from the choice of renormalization and factorization scales in the physics modeling of tt events is determined by producing two additional simulated event samples.These samples are generated with both scales simultaneously varied by a factor of two up or down from their default values equal to the Q of the hard process in the event; where the sum is over all additional final-state partons in the matrix element.The effect of varying the renormalization and factorization scales in the W+jets and Z+jets samples is also considered to determine the uncertainty in the shape of this background.The uncertainty arising from the choice of parton shower matching threshold in the event generation is determined in a similar fashion, using additional samples in which the threshold is varied up or down.Uncertainties from the modeling of the hadronization are evaluated by comparing POWHEG v1 simulated samples with two different hadron shower generators (PYTHIA and HERWIG).The uncertainty owing to the choice of the PDF is determined by reweighting the simulated events and repeating the analysis using the 44 CTEQ6L PDF error sets [30].The maximum variation is taken as the uncertainty.Simulated samples with the top quark mass varied by ±1 GeV, which corresponds to the precision of the measured top quark mass [49], are generated to evaluate the effect of the uncertainty in this parameter.The effect of reweighting the top quark p T spectrum in simulation, as described in Section 3, is found to have a negligible effect for low values of the kinematic event variables, and increases to 3-7% for the highest values.
Other uncertainties are associated with imperfect understanding of efficiencies, resolutions, and scales describing the detector response.The uncertainty arising from each source is estimated, and the analysis repeated with each corresponding parameter varied within its uncertainty.
The efficiencies and associated uncertainties for triggering and lepton identification are determined from data by a tag-and-probe method [50].The probabilities for identifying and misidentifying b jets in the simulation are compared to those measured in data, and the resulting correction factors and their uncertainties are determined as a function of jet energy and quark flavor.The uncertainties in the correction factors are typically 2%.
The uncertainty in the jet energy scale (JES) is determined as a function of the jet p T and η [41], and an uncertainty of 10% is included in the jet energy resolution (JER) [41].The effect of this limited knowledge of the JES and JER is determined by varying the JES and JER in the simulated samples within their uncertainties.The uncertainty in the JES and JER, as well as uncertainties in the electron, photon, tau, and muon energy scale, are propagated into the calculation of E miss T .The uncertainty in the electron and photon energy scale is 0.6% in the barrel, and 1.5% in the endcap [38].The uncertainty in the tau lepton energy scale is estimated to be ±3% [51], while the effect of the uncertainty in the muon momentum measurement is found to be negligible.A 10% uncertainty is assigned to the estimate of the nonclustered energy used in the calculation of E miss T [46].
The effect of the uncertainty in the level of pileup is estimated by varying the inelastic pp cross section used in the simulation by ±5% [52].
The uncertainty in the normalization of the background is determined by varying the normalization of the single top, W+jets, and Z+jets processes by ±30%, and the QCD multijet processes by ±100%.The uncertainty in the shape of the QCD multijet distribution in the electron channel is estimated by using an alternative control region in data to determine the contribution of QCD multijet events.This uncertainty is found to have a negligible effect.
The dominant systematic effects are caused by the uncertainties in the modeling of the hadronization and the tt signal.For illustrative purposes, typical systematic uncertainties in the √ s = 8 TeV results coming from each of the sources described above are presented in Table 1.The values shown for each kinematic event variable are the median uncertainties over all of the bins for that variable.
Table 1: Typical relative systematic uncertainties in percent (median values) in the normalized tt differential cross section measurement as a function of the four kinematic event variables at a center-of-mass energy of 8 TeV (combination of electron and muon channels).Typical values of the total systematic uncertainty are also shown.

Uncertainty source
Relative

Results
The normalized differential tt cross sections as a function of each of the kinematic event variables are shown in Figs. 3 and 4  The data distributions in the figures are compared with the predictions from the event generators in the left-hand plots: MADGRAPH and POWHEG V2 with two different hadron shower generators, PYTHIA and HERWIG.For the √ s = 8 TeV results, the predictions from the MC@NLO and POWHEG V1 generators are also shown.The effect on the predicted distributions from varying the modeling parameters (the matching threshold and renormalization scale Q 2 ) up and down by a factor of two for the MADGRAPH event generator is shown in the right-hand plots for the two MADGRAPH simulations.The uncertainties shown by the vertical bars on the points in the figures and given in the tables include both the statistical uncertainties and those resulting from the unfolding procedure.
The measurements at √ s = 7 TeV are well described by all the event generators in the distribution of E miss T .For S T , p W T , and H T , the event generators predict a somewhat harder spectrum than seen in data.However, the POWHEG V2 + PYTHIA event generator provides a reasonable description of the H T and S T differential cross sections.
The results at √ s = 8 TeV are generally well described by the MC@NLO and the POWHEG V2 + PYTHIA event generators.The POWHEG V2 + HERWIG event generator describes the E miss T and p W T distributions well.However, for H T and S T this event generator predicts a harder spectrum than seen in data, at both center-of-mass energies.
The MADGRAPH event generator generally predicts a harder spectrum than seen in data for all variables.The variations in matching threshold and Q 2 in the MADGRAPH event generator are not sufficient to explain this difference between the prediction and data.However, the MADGRAPH event generator provides a good description of the data after reweighting the top quark p T spectrum, as described in Section 3. The prediction obtained from the MADGRAPH event generator after the reweighting is shown on all the plots.

Summary
A measurement of the normalized differential cross section of top quark pair production with respect to the four kinematic event variables E miss T , H T , S T , and p W T has been performed in pp collisions at a center-of-mass energy of 7 TeV using 5.0 fb −1 and at 8 TeV using 19.7 fb −1 of data collected by the CMS experiment.
This study confirms previous CMS findings that the observed top quark p T spectrum is softer than predicted by the MADGRAPH, POWHEG, and MC@NLO event generators, but otherwise there is broad consistency between the MC event generators and observation.This result provides confidence in the description of tt production in the SM and its implementation in the most frequently used simulation packages.acknowledge the computing centers and personnel of the Worldwide LHC Computing Grid for delivering so effectively the computing infrastructure essential to our analyses.Finally, we acknowledge the enduring support for the construction and operation of the LHC and the CMS detector provided by the following funding agencies: BMWFW and FWF (Austria); FNRS and FWO (Belgium); CNPq, CAPES, FAPERJ, and FAPESP (Brazil); MES (Bulgaria); CERN;       [5] CMS Collaboration, "Measurement of the differential cross section for top quark pair production in pp collisions at √ s = 8 TeV", Eur.Phys.J. C 75 (2015) 542, doi:10.1140/epjc/s10052-015-3709-x,arXiv:1505.04480.
[7] CMS Collaboration, "Measurement of the differential tt production cross section for high-p T top quarks in e/µ+jets final states at 8 TeV", CMS Physics Analysis Summary CMS-PAS-TOP-14-012, CERN, 2015.

A Additional tables
The measured values of the tt differential cross sections as a function of E miss T , H T , S T , and p W T for √ s = 7 TeV and √ s = 8 TeV are given in the tables below, along with their statistical, systematic, and total uncertainties.variable at a center-of-mass energy of 8 TeV (combination of electron and muon channels).The rightmost three columns show the relative uncertainties on the measured values, in percent.The statistical and systematic uncertainties are listed separately, and are combined in quadrature to give the overall relative uncertainty.
Figures 1 and 2 show the observed distributions of E miss T , H T , S T , and p W T , in the √ s = 8 TeV data samples, compared to the sum of the corresponding signal and background distributions from simulation.

Figure 1 :
Figure 1: The observed distributions of E miss T (top) and H T (bottom) in the √ s = 8 TeV elec-tron+jets (left) and muon+jets (right) data samples, compared to predictions from simulation.The points are the data histograms, with the vertical bars showing the statistical uncertainty, and the predictions from the simulation are the solid histograms.The shaded region shows the uncertainty in the values from simulation.These include contributions from the statistical uncertainty and the uncertainty in the tt cross section.The lower plots show the ratio of the number of events from data and the prediction from the MC simulation.

Figure 2 :
Figure 2: The observed distributions of S T (top) and p W T (bottom) in the √ s = 8 TeV elec-tron+jets (left) and muon+jets (right) data samples, compared to predictions from simulation.The points are the data histograms, with the vertical bars showing the statistical uncertainty, and the predictions from the simulation are the solid histograms.The shaded region shows the uncertainty in the values from simulation.These include contributions from the statistical uncertainty and the uncertainty in the tt cross section.The lower plots show the ratio of the number of events from data and the prediction from the MC simulation.

for the √ s = 7
TeV data, and in Figs. 5 and 6 for the √ s = 8 TeV data.The results are also presented in Tables A.1-A.8 of Appendix A.

Figure 3 :
Figure 3: Normalized E miss T (top) and H T (bottom) differential tt cross sections from the combined electron and muon data at √ s = 7 TeV.The vertical bars on the data points represent the statistical and systematic uncertainties added in quadrature.The inner section of the vertical bars, denoted by the tick marks, show the statistical uncertainty.Left: comparison with different simulation event generators: MADGRAPH + PYTHIA (both the default and after reweighting the top quark p T spectrum), POWHEG V2 + HERWIG, and POWHEG V2 + PYTHIA.Right: comparison with predictions from the MADGRAPH + PYTHIA event generator found by varying the matching threshold and renormalization scales (µ R , µ F ) up and down by a factor of two.The lower plots show the ratio of the predictions to the data, with the statistical and total uncertainties in the ratios indicated by the two shaded bands.

Figure 4 :
Figure 4: Normalized S T (top) and p W T (bottom) differential tt cross sections from combined electron and muon data at √ s = 7 TeV.The vertical bars on the data points represent the statistical and systematic uncertainties added in quadrature.The inner section of the vertical bars, denoted by the tick marks, show the statistical uncertainty.Left: comparison with different simulation event generators: MADGRAPH + PYTHIA (both the default and after reweighting the top quark p T spectrum), POWHEG V2 + HERWIG, and POWHEG V2 + PYTHIA.Right: comparison with predictions from the MADGRAPH + PYTHIA event generator found by varying the matching threshold and renormalization scales (µ R , µ F ) up and down by a factor of two.The lower plots show the ratio of the predictions to the data, with the statistical and total uncertainties in the ratios indicated by the two shaded bands.

Figure 5 :
Figure 5: Normalized E miss T (top) and H T (bottom) differential tt cross sections from combined electron and muon data at √ s = 8 TeV.The vertical bars on the data points represent the statistical and systematic uncertainties added in quadrature.The inner section of the vertical bars, denoted by the tick marks, show the statistical uncertainty.Left: comparison with different simulation event generators: MADGRAPH + PYTHIA (both the default and after reweighting the top quark p T spectrum), MC@NLO + HERWIG, POWHEG V1 + HERWIG, POWHEG V1 + PYTHIA, POWHEG V2 + HERWIG, and POWHEG V2 + PYTHIA.Right: comparison with predictions from the PYTHIA event generator found by varying the matching threshold and renormalization scales (µ R , µ F ) up and down by a factor of two.The lower plots show the ratio of the predictions to the data, with the statistical and total uncertainties in the ratios indicated by the two shaded bands.

Figure 6 :
Figure 6: Normalized S T (top) and p W T (bottom) differential tt cross sections from combined electron and muon data at √ s = 8 TeV.The vertical bars on the data points represent the statistical and systematic uncertainties added in quadrature.The inner section of the vertical bars, denoted by the tick marks, show the statistical uncertainty.Left: comparison with different simulation event generators: MADGRAPH + PYTHIA (both the default and after reweighting the top quark p T spectrum), MC@NLO + HERWIG, POWHEG V1 + HERWIG, POWHEG V1 + PYTHIA, POWHEG V2 + HERWIG, and POWHEG V2 + PYTHIA.Right: comparison with predictions from the MADGRAPH + PYTHIA event generator found by varying the matching threshold and renormalization scales (µ R , µ F ) up and down by a factor of two.The lower plots show the ratio of the predictions to the data, with the statistical and total uncertainties in the ratios indicated by the two shaded bands.
Individuals have received support from the Marie-Curie program and the European Research Council and EPLANET (European Union); the Leventis Foundation; the A. P. Sloan Founda-tion; the Alexander von Humboldt Foundation; the Belgian Federal Science Policy Office; the Fonds pour la Formation à la Recherche dans l'Industrie et dans l'Agriculture (FRIA-Belgium); the Agentschap voor Innovatie door Wetenschap en Technologie (IWT-Belgium); the Ministry of Education, Youth and Sports (MEYS) of the Czech Republic; the Council of Science and Industrial Research, India; the HOMING PLUS program of the Foundation for Polish Science, cofinanced from European Union, Regional Development Fund, the Mobility Plus program of the Ministry of Science and Higher Education, the OPUS program contract 2014/13/B/ST2/02543 and contract Sonata-bis DEC-2012/07/E/ST2/01406 of the National Science Center (Poland); the Thalis and Aristeia programmes cofinanced by EU-ESF and the Greek NSRF; the National Priorities Research Program by Qatar National Research Fund; the Programa Clarín-COFUND del Principado de Asturias; the Rachadapisek Sompot Fund for Postdoctoral Fellowship, Chulalongkorn University and the Chulalongkorn Academic into Its 2nd Century Project Advancement Project (Thailand); and the Welch Foundation, contract C-1845.

Table A .
1: Normalized tt differential cross section measurements with respect to the E miss T variable at a center-of-mass energy of 7 TeV (combination of electron and muon channels).The rightmost three columns show the relative uncertainties on the measured values, in percent.The statistical and systematic uncertainties are listed separately, and are combined in quadrature to give the overall relative uncertainty.

Table A .
2: Normalized tt differential cross section measurements with respect to the H T variable at a center-of-mass energy of 7 TeV (combination of electron and muon channels).The rightmost three columns show the relative uncertainties on the measured values, in percent.The statistical and systematic uncertainties are listed separately, and are combined in quadrature to give the overall relative uncertainty.

Table A .
3: Normalized tt differential cross section measurements with respect to the S T variable at a center-of-mass energy of 7 TeV (combination of electron and muon channels).The rightmost three columns show the relative uncertainties on the measured values, in percent.The statistical and systematic uncertainties are listed separately, and are combined in quadrature to give the overall relative uncertainty.

Table A .
4: Normalized tt differential cross section measurements with respect to the p W T variable at a center-of-mass energy of 7 TeV (combination of electron and muon channels).The rightmost three columns show the relative uncertainties on the measured values, in percent.The statistical and systematic uncertainties are listed separately, and are combined in quadrature to give the overall relative uncertainty.

Table A .
5: Normalized tt differential cross section measurements with respect to the E miss T

Table A .
6: Normalized tt differential cross section measurements with respect to the H T variable at a center-of-mass energy of 8 TeV (combination of electron and muon channels).The rightmost three columns show the relative uncertainties on the measured values, in percent.The statistical and systematic uncertainties are listed separately, and are combined in quadrature to give the overall relative uncertainty.

Table A .
7: Normalized tt differential cross section measurements with respect to the S T variable at a center-of-mass energy of 8 TeV (combination of electron and muon channels).The rightmost three columns show the relative uncertainties on the measured values, in percent.The statistical and systematic uncertainties are listed separately, and are combined in quadrature to give the overall relative uncertainty.TableA.8:Normalized tt differential cross section measurements with respect to the p W T variable at a center-of-mass energy of 8 TeV (combination of electron and muon channels).The rightmost three columns show the relative uncertainties on the measured values, in percent.The statistical and systematic uncertainties are listed separately, and are combined in quadrature to give the overall relative uncertainty.