Measurement of the top quark pair production cross section in dilepton final states containing one tau lepton in pp collisions at 13 TeV

The cross section of top quark pair production is measured in the tt → (`ν`)(τhντ )bb final state, where τh refers to the hadronic decays of the τ lepton, and ` is either an electron or a muon. The data sample corresponds to an integrated luminosity of 35.9 fb−1 collected in proton-proton collisions at √ s = 13 TeV with the CMS detector. The measured cross section is σtt = 781± 7 (stat)± 62 (syst)± 20 (lumi) pb, and the ratio of the partial width Γ(t → τντ b) to the total decay width of the top quark is measured to be 0.1050 ± 0.0009 (stat) ± 0.0071 (syst). This is the first measurement of the tt production cross section in proton-proton collisions at √ s = 13 TeV that explicitly includes τ leptons. The ratio of the cross sections in the `τh and `` final states yields a value R`τh/`` = 0.973± 0.009 (stat)± 0.066 (syst), consistent with lepton universality. Submitted to the Journal of High Energy Physics c © 2019 CERN for the benefit of the CMS Collaboration. CC-BY-4.0 license ∗See Appendix A for the list of collaboration members ar X iv :1 91 1. 13 20 4v 1 [ he pex ] 2 9 N ov 2 01 9


Introduction
In proton-proton (pp) collisions at the CERN LHC, top quarks are produced mainly in pairs (tt) and subsequently decay to b quarks and W bosons: pp → tt → W + bW − b. The decay modes of the two W bosons determine the observed event signature. The dilepton decay channel denotes the case where both W bosons decay leptonically. In this paper, we consider the process tt → ( ν )(τν τ )bb, where one W boson decays into ν where is either an electron (e) or a muon (µ), and the other into a tau lepton and a neutrino (τν τ ). The expected fraction of events in this final state corresponds to ≈4/81 (≈5%) of all tt decays, i.e. equivalent to the fraction of all light dilepton channels (ee, µµ, eµ).
Recent checks of lepton flavour universality violation [1-8] sparked a renewed interest towards measurements involving τ leptons, owing to a potential disagreement with standard model (SM) predictions. The t → (τν τ )b decay exclusively involves third-generation leptons and quarks which, owing to their large masses, may be particularly sensitive to beyond SM contributions. For example, the existence of a charged Higgs boson [9-12] may give rise to anomalous τ lepton production that could be observed in this decay channel. This is the first measurement of the tt production cross section in pp collisions at √ s = 13 TeV that explicitly includes τ leptons. The data sample was collected in 2016 with the CMS detector at the LHC and corresponds to an integrated luminosity of 35.9 fb −1 . The τ lepton is identified through its visible decay products, either hadrons (τ h ) or leptons (τ ), with the corresponding branching fractions B(τ h → hadrons + ν τ ) ≈ 65% and B(τ → ν ν τ ) ≈ 35%. In the first case, the τ h decays into a narrow jet with a distinct signature, whereas the leptonic decays are difficult to distinguish from prompt electron or muon production. In this measurement, the signal includes only τ leptons that decay hadronically, and does not include leptons from τ decays. The dominant background contribution comes from events where a jet is misidentified as a τ h , mostly from tt lepton+jets events, i.e. tt → ( ν )(qq )bb. The cross section is measured by performing a profile likelihood ratio (PLR) fit [13] to the transverse mass of the system containing the lepton (e or µ) and the missing transverse momentum, in two kinematic categories of the selected events for each of the eτ h and µτ h final states. The cross section is measured in the fiducial phase space of the detector and also extrapolated to the full phase space. The ratio of the cross sections in the τ and light dilepton [14] final states σ tt ( τ)/σ tt ( ), and the ratio of the partial to the total decay width of the top quark Γ(t → τν τ b)/Γ total are evaluated. This paper is organized as follows: the CMS detector layout is briefly described in Section 2; details about the simulated event samples used in the data analysis are provided in Section 3; Section 4 covers the event reconstruction and the event selection; the event categorization and the fit procedures are described in Section 5; the background determination procedure is given in Section 6; the description of the systematic uncertainties is presented in Section 7; measurements of the cross sections, and the ratio of the partial to the total tt decay width are discussed in Section 8; and the results are summarized in Section 9.

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, covering 0 < ϕ < 2π in azimuth and |η| < 2.5 in pseudorapidity. Forward calorimeters extend the pseudorapidity coverage provided by the barrel and endcap detectors. Muons are detected in gas-ionization chambers embedded in the steel flux-return yoke outside the solenoid. The detector is nearly hermetic, providing reliable measurement of the momentum imbalance in the plane transverse to the beams. A two-level trigger system [15] selects the most interesting pp collision events for use in physics analysis. 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. [16].

Event simulation
The analysis makes use of simulated samples of tt events, as well as other processes that result in reconstructed τ leptons in the final state. These samples are used to design the event selection, to calculate the acceptance for tt events, and to estimate most of the backgrounds in the analysis.
Signal tt events are simulated with the POWHEG event generator (v2) [17][18][19][20][21] at next-to-leadingorder (NLO) accuracy in quantum chromodynamics (QCD). The parton showers are modelled using PYTHIA (v8.2) [22] with the CUETP8M2T4 underlying event (UE) tune [23]. The background samples used in the measurement of the cross section are simulated with POWHEG and MADGRAPH5 aMC@NLO (v2.2.2) [24]. The MADGRAPH5 aMC@NLO generator with MLM matching [25] is used for the simulation of W boson production in association with jets (W+jets), and Drell-Yan (DY) production in association with jets at leading-order (LO) accuracy. Here, only the leptonic decays of DY events and W bosons are simulated, and up to four additional jets are included. The diboson processes are produced with NLO accuracy: WW with POWHEG, WZ and ZZ with MADGRAPH5 aMC@NLO with FxFx matching [26]. The POWHEG generator is used for the simulation of t-channel single top quark production and single top quark production associated with a W boson (tW) [27,28]. The single top quark s-channel sample is produced with MADGRAPH5 aMC@NLO at NLO accuracy with FxFx matching scheme. The simulated events are produced with a top quark mass of m t = 172.5 GeV. The generated events are subsequently processed with PYTHIA using the underlying event tune CUETP8M1 to provide the showering of the partons, and to perform the matching of the soft radiation with the contributions from direct emissions included in the matrix-element (ME) calculations. The default parton distribution functions (PDFs) are the NNPDF3.0 [29]. The τ decays are simulated with PYTHIA, which correctly accounts for the τ lepton polarization in describing the kinematic properties of the decay. The CMS detector response is simulated with GEANT4 [30]. Additional pp interactions in the same or nearby bunch crossings (pileup, PU) are superimposed on the hard collision. Simulated events are reweighted to match the distribution of the number of pileup collisions per event in data. This distribution is derived from the instantaneous luminosity and the inelastic cross section [31].
The next-to-next-to-leading-order (NNLO) expected SM tt pair production cross section of 832 +20 −29 (scale) ± 35 (PDF+α S ) pb [32] (m t = 172.5 GeV) is used for the normalization of the number of tt events in the simulation. The first uncertainty includes the uncertainties in the factorization and renormalization scales, while the second is associated with possible choices of PDFs and the value of the strong coupling constant (α S ). The proton structure is described by the CT14 (NNLO) PDF set with the corresponding PDF and α S uncertainties [33]. The W+jets and DY+jets backgrounds are normalized to their NNLO cross sections calculated with FEWZ (v3.1) [34]. The t-channel and the s-channel single top quark production are normalized to the NLO calculations obtained from HATHOR (v2.1) [35,36]. The production of tW is normalized to the NNLO calculation [37,38]. Finally, the production of diboson pairs is normalized to the NLO cross section prediction calculated with MCFM [39,40] (v7.0).

Event reconstruction and selection
The signal event topology is defined by the presence of two b quark jets from the top quark decays, one W boson decaying leptonically into eν or µν, and a second W boson decaying into τ h ν. In each event, all objects are reconstructed with a particle-flow (PF) algorithm [41]. This algorithm combines the information from all subdetectors to identify and reconstruct all types of particles in the event, namely charged and neutral hadrons, photons, muons, and electrons, together referred to as PF objects. These objects are used to construct a variety of higher-level objects and observables, including jets and missing transverse momentum ( p miss T ), which is the negative vector sum of transverse momenta of all reconstructed PF objects. Parameters of jets and the tracks associated with jets provide input variables for b tagging discriminators. The reconstructed vertex with the largest value of summed physics-object p 2 T is taken to be the primary pp interaction vertex. Jets are reconstructed by clustering PF objects with the antik T [42] jet algorithm with a distance parameter R = 0.4.
Electron or muon candidates are required to originate from the primary vertex, pass quality selection criteria, and be isolated relative to other activity in the event. The relative isolation is based on PF objects within a cone of ∆R = √ (∆η) 2 + (∆ϕ) 2 = 0.4 around the electron or muon track, and defined as I rel = (E ch + E nh + E ph − 0.5 E PU ch )/p T , where E ch is the transverse energy deposited by charged hadrons from the primary vertex, E nh and E ph are the respective transverse energies of the neutral hadrons and photons, and 0.5 E PU ch is the estimation of the contribution of neutral particles from pileup vertices, calculated as half of the energy of the charged particles from pileup; p T is the electron or muon transverse momentum. Electron candidates with I rel < 0.0588 in the barrel or I rel < 0.0571 in the endcaps are considered isolated. The muon candidate is isolated if I rel < 0.15 in either the barrel or the endcaps. The lepton isolation requirements are used to suppress backgrounds from to multijet production. The charge misidentification probability for electrons and muons is less than 0.5% and 0.1%, respectively, and is measured from Z boson decays and simulation [43][44][45].
Hadronic τ lepton decays are reconstructed with the hadron-plus-strips (HPS) algorithm [46], which starts from reconstructed jets. In each jet, a charged hadron is combined with other nearby charged hadrons or photons to identify the decay modes. The identification of π 0 mesons is enhanced by clustering electrons and photons in "strips" along the track bending direction to take into account possible broadening of calorimeter signatures by early showering photons. The τ h candidates are selected from the following combinations of charged hadrons and strips that correspond to the τ decay modes: single hadron, hadron plus a strip, hadron plus two strips, and three hadrons. A multivariate analysis of these HPS τ h candidates is used to reduce the contamination from quark and gluon jets. A boosted decision tree is trained using a sample of DY events with τ h decays as signal and a sample of QCD multijet events as background, both from simulation. Input variables include the multiplicity and the transverse momenta of electron and photon candidates in the vicinity of the τ h , the kinematic properties of hadrons and strips, and the τ h lifetime information, such as the impact parameter of the leading track and the significance of the length of flight to the secondary vertex of the τ h candidates with three charged hadrons. Additional requirements are applied to discriminate genuine τ h leptons from prompt electrons and muons. The τ h charge is taken as the sum of the charges of the corresponding charged hadrons. The misidentification probability for the charge is less than 1% and it is estimated from Z → ττ → µτ h data events with same-charge µ and τ h . The τ h identification efficiency of this algorithm is estimated to be approximately 60% for p T > 20 GeV, and it is measured in a sample enriched in Z → ττ → µτ h data events with a "tag-and-probe" technique [47]. The corresponding probability for generic hadronic jets to be misidentified as τ h is less than 1% [46].
For the eτ h (µτ h ) final state, data are collected with a trigger requiring at least one isolated electron (muon) with a threshold of p T > 27 (24) GeV.
Events are selected by requiring one isolated electron (muon) with transverse momentum p T > 30 (26) GeV and |η| < 2.4, at least two jets with p T > 30 GeV and |η| < 2.5, and exactly one τ h candidate with p T > 30 GeV and |η| < 2.4. The τ h candidate and the selected lepton are required to have opposite electric charges (OC). Electrons or muons are required to be separated from any jet and from the τ h candidate in the η-ϕ plane by a distance ∆R > 0.4. Events with any additional loosely isolated electron (muon) of p T > 15 (10) GeV are rejected. An electron is considered loosely isolated if I rel < 0.0994 in the barrel or I rel < 0.107 in the endcaps. A muon is loosely isolated if I rel < 0.25 in either the barrel or the endcaps. At least one jet is required to be identified as originating from b quark hadronization ("b tagged"). The b tagging algorithm used ("CSVv2" in Ref. [48]) combines the information of displaced tracks and secondary vertices associated with the jet in a multivariate technique. The working point selected provides a b tagging efficiency of about 66% with a corresponding light-flavour misidentification rate of 1%. The selected events exhibit good agreement between the observed data and the expectation, as shown in Fig. 1 for the p T distribution of the τ h candidate. The dominant background contribution comes from other tt decays, mostly from lepton+jets final states where a jet is misidentified as a τ h candidate.

Event categories and fit procedure
The tt production cross section is extracted from a PLR fit of the binned distribution of the transverse mass of the lepton and p miss T in two kinematic event categories, for each of the eτ h and µτ h final states. The transverse mass is defined as where ∆ϕ is the azimuthal angle difference between the lepton transverse momentum vector, p T , and p miss T . The m T distribution provides separation between signal and background processes (as shown in Fig. 2) and does not significantly depend on p T and η of the τ candidate, or other jet characteristics in the kinematic ranges of this study. Two event categories are defined according to the kinematic properties of jets in the event. In order to discriminate against the main background of misidentified τ h from the tt lepton+jets process, the constraints from top quark and W boson masses in the decay t → bW → b(qq ) are used. Jet triplets are constructed for each combination of one b-tagged jet and two untagged jets, chosen from all jets in the event, including the τ h candidate. The distance parameter for each triplet is calculated .5 GeV and m W = 80.385 GeV are, respectively, the masses of the top quark and of the W boson [49], m jj is the invariant mass of the two untagged jets, and m jjb is the invariant mass of the jet triplet. The event is assigned to the "signal-like" category if there is only one untagged jet, or if the minimum parameter value D min jjb is larger than 60 GeV. Otherwise, it is assigned to the "background-like" event category. The threshold of 60 GeV provides an optimal separation of signal and background event categories, together with a maximization of the yields in each of the two categories in order to reduce the statistical uncertainties. In the fit, the two event categories provide an additional constraint on the background processes independent from the details of the m T distribution.   In the m T distribution, the signal may extend beyond the W boson mass endpoint because of the two-neutrino final state, whereas the background process cannot. The last bin in both distributions includes overflow events. In the D min jjb distribution, the downward arrow points at the threshold of the cut used (D min jjb >60 GeV), and the panel on the right shows the fraction of events in the "signal-like" category where there is only one untagged jet.
The cross section is derived from the signal strength measured in the fit, i.e. its ratio to the value expected in the SM. It is estimated for both event categories, in each of the eτ h and µτ h final states. The expected number of events in a given bin of the m T distribution is parametrized as a function of signal strength and nuisance parameters. The nuisance parameters encode the effects of systematic uncertainties. The signal strength is a free parameter in the fit. The fitted variables do not significantly depend on the kinematic properties of the τ lepton in the specific tt signal model considered here, i.e. tt → ( ν )(τ h ν τ )bb. The likelihood function is defined as a product of Poisson distributions of the expected number of events in bins of the m T distribution and nuisance constraints. Based on the likelihood function, the PLR test statistic is defined as the ratio between the maximum of the likelihood for a given value of signal strength and the global maximum of the likelihood function. The effect of the systematic uncertainties in the signal strength is determined with this approach.

Background estimate
The main background contribution comes from events with one lepton, significant p miss T , and three or more jets, dominated by the lepton+jets tt process, where one of the jets is falsely identified as a τ h . Misidentified τ h candidates also come from multijet and W+jet background processes. There is a small contribution from processes with genuine hadronic τ h : tW single top quark production, τ τ h from DY decays, tt → τ τ h bb, and diboson processes. All processes, except multijet, are estimated from simulation after applying appropriate corrections. The pileup, trigger efficiencies, lepton identification, jet energy corrections, and b tagging efficiencies in the simulation are corrected with scale factors measured in separate publications [43, 45, 50], as described in Section 7.
The τ h misidentification contribution is determined by constraining the falsely identified τ h in the overall fit to the data in the m T distribution. In the fit, the event yields of the background processes with a misidentified τ h are determined by adjusting the normalization of the shapes of the m T distributions. The normalization factors are introduced as nuisance parameters with constraints determined from studies in other processes [46]. The corresponding uncertainties are discussed in Section 7.
The background from the multijet processes is determined from data as it provides a more accurate description with a smaller statistical uncertainty. The shape of the m T distribution is obtained from a sample of events containing lepton and τ h candidates of the same charge (SC). It is estimated by subtracting from the data all other processes taken from simulation, including the fully hadronic final states in tt, single top quark, and dibosons. The m T shapes for SC and OC events are the same within the uncertainties in a control region with a relaxed τ h identification requirement, and in agreement with the simulation. The normalization is corrected by multiplying the SC m T distribution by the OC-to-SC ratio, f OC/SC , as determined in a control region from events with a relaxed τ h identification and an inverted lepton isolation requirement, where the multijet contribution is dominant. All other event selection requirements remain the same as in the main selection. The ratio is measured to be f OC/SC = 1.05 ± 0.05 (stat + syst), in agreement with simulation. As one of the processes with misidentified τ h , the normalization of the multijet contribution is varied in the fit as a separate nuisance parameter, as described in Section 7.

Systematic uncertainties
The main sources of systematic uncertainty are from τ h identification and misidentification, b tagging, estimation of pileup in the pp collisions, jet energy scale (JES), and jet energy res-olution (JER). Other sources of uncertainty are from lepton identification, trigger efficiency, and the calibration of the integrated luminosity. Theoretical uncertainties are also included in the event simulation. Uncertainties are applied in a coherent way to signal and background processes. The corresponding corrections and their uncertainties are measured in dedicated studies, which are described below.
The uncertainty in the efficiency of τ h identification is 5% for all τ h with p T > 20 GeV and is applied to all processes with a genuine τ h . It is measured with a tag-and-probe technique in samples enriched in Z → τ τ h events [46]. The τ h charge confusion probability, estimated to be less than 1%, is considered a part of the τ h identification efficiency uncertainty. The correction to the reconstructed energy of the τ h jet (τ energy scale) and the corresponding uncertainty is estimated in a fit of the data in distributions sensitive to the τ energy, such as the τ h visible mass [46]. The dominant background contribution arises from processes where a jet is misidentified as τ h , mainly lepton+jets tt, W+jets, and multijet production. The τ h misidentification probability and its uncertainty in these processes are directly measured in the fit. The misidentification probability is varied within ±50% of the expected values in all processes with a jet falsely identified as the τ h candidate. The variation covers the differences between expected and observed misidentification probabilities and the possible dependence on other kinematic properties of the τ h candidate [46]. The misidentification probability is significantly constrained in the fit and is not the dominant source of the uncertainty in the final result. The cross section is measured by the fit in the fiducial phase space of the detector. The fiducial cross section is extrapolated to the full phase space by correcting for the acceptance of the tt signal process. The fit and the acceptance include the following modelling uncertainties: the renormalization and factorization scales, and PDFs including α S . The uncertainty in the PDF is estimated by using the CT14 (NNLO) set as alternative PDFs. The renormalization and factorization scales in the ME calculations are varied independently by factors of 0.5 and 2.0 from their nominal values, and the envelope of the variations is included in the measurement. The scale is varied by factors of 0.5 and 2.0 in the parton shower (PS) simulation of final-state and initial-state radiation, FSR and ISR. The h damp parameter regulating the real emissions in POWHEG (ME-PS matching) is varied from its central value of 1.58 m t using samples with h damp set to 0.99 m t and 2.24 m t (m t = 172.5 GeV), as obtained from tuning this parameter to tt data at √ s = 8 TeV [61]. The underlying event tune is varied within its uncertainties [23,61]. The effect of these uncertainties on the final state objects is included in the fit in the fiducial phase space by adding the corresponding systematic variations normalized to the nominal acceptance. The uncertainties in the fit are not correlated with the acceptance uncertainty in the extrapolation to the full phase space.
The theoretical uncertainties are implemented by reweighting the simulated events with corresponding scale factors. The differences between weighted and unweighted distributions are taken as the uncertainties in the modelling. Separate data sets with varied parameters are used for determining FSR, ISR, ME-PS matching, and underlying event uncertainties.
The impact of the systematic uncertainties on the measurement is given in Table 2.

Results
The event yields expected from the signal and background processes, as well as the observed event yields are summarized in Table 1, for the signal-like and the background-like event categories (described in Section 5) in each of the eτ h and µτ h final states. The observed event yields in data show good agreement with the prediction. The m T distributions in the two categories of the selected events are shown in Fig. 3, for both the eτ h and µτ h final states. These distributions show good agreement between the observed event yields and the expected numbers of signal and background events. Table 1: Expected and observed event yields in the τ h ( = e, µ) final state for signal and background processes for an integrated luminosity of 35.9 fb −1 . Statistical and systematic uncertainties are shown. The expected prefit contributions of all processes are presented separately for background-like and signal-like event categories. The statistical uncertainties of the modelling are shown for the processes estimated from the simulation. The multijet contribution and the corresponding statistical uncertainties are estimated using data, as described in Section 6. Events / GeV  Table 2 lists the systematic uncertainties in the signal strength after the fit. The effect of the uncertainties on the signal strength is estimated by a likelihood scan where only one nuisance parameter (or a group of them) is varied at once while the others are fixed to their nominal postfit values. The largest experimental uncertainties are from τ h identification and misidentification, and pileup estimation. The largest theoretical uncertainties are due to the modelling of top quark p T in tt processes, b quark fragmentation, and PS modelling (ISR and FSR).
The fiducial cross section for the production of tt events is defined in the acceptance region of kinematic phase space defined by the selection criteria described earlier. The estimate of the fiducial cross section includes the branching fractions of the final states, trigger, lepton identification and isolation, and the overall reconstruction efficiency. The cross sections in the fiducial phase space for the individual eτ h and µτ h final states, as well as the τ h combined final state, are measured from the PLR fit to be: σ fid tt (µτ h ) = 135.2 ± 1.5 (stat) ± 9.9 (syst) ± 3.4 (lumi) pb, σ fid tt ( τ h ) = 134.5 ± 1.2 (stat) ± 9.5 (syst) ± 3.4 (lumi) pb. (3) The acceptance A tt is the fraction of signal events in the fiducial phase space, and it is determined with respect to all signal events in the nominal tt simulation. It includes kinematic selection cuts and is evaluated for the different signal final states as: A tt ( τ h ) = 0.1722 ± 0.0003 (stat) ± 0.0062 (syst), where the systematic uncertainties include the uncertainties of the modelling as described in Section 7 and listed as "Extrapolation uncertainties" in Table 2.
The cross section values in the full phase space are obtained from the extrapolation of the fiducial cross sections using the acceptances A tt estimated from the simulation: σ tt (µτ h ) = 770 ± 8 (stat) ± 63 (syst) ± 20 (lumi) pb, (8) The expected and observed dependence of the likelihood on the cross section in the full phase space in the τ h combined final state are shown in Fig. 4. The result of the fit is consistent with the predicted SM tt production cross section of 832 +20 −29 (scale) ± 35 (PDF+α S ) pb [32]. Using simulated tt samples with different m t values, we find that the cross section changes by 1.5% per ∆m t = 1 GeV.
The ratio of the cross section in the τ h final state divided by the cross section measured in the dilepton final state in the same data-taking period [14] yields a value of R τ h / = 0.973 ± 0.009 (stat) ± 0.066 (syst), consistent with unity as expected from lepton flavour universality. The relative systematic uncertainty in the ratio is 6.8%. About 5% comes from the uncertainties in the τ h identification (4.5%) and misidentification probability in tt events (2.3%). The rest comes from the other uncorrelated uncertainties in the ratio and the treatment of the correlated uncertainties in the calculation of the ratio. In particular, a small contribution comes from the uncertainties in the extrapolation to the full phase space that are considered uncorrelated because the two measurements extrapolate from different fiducial phase Table 2: Systematic and statistical uncertainties determined from the fit to the data in the eτ h and µτ h final states, and their combination. Uncertainties are grouped by their origin: experimental, theoretical, and extrapolation. The uncertainties in the measurement in the dilepton final state [14] used in the partial width ratio estimate are also quoted (column "dileptons"), where the asymmetric extrapolation uncertainties are symmetrized by adding them in quadrature. As both measurements use the same data, some uncertainties are correlated, as shown in the last column.

Source
Uncertainty spaces. Also, the triggers are not the same. The measurement also provides an estimate of the ratio of the partial to the total width of the top quark decay, R Γ = Γ(t → τν τ b)/Γ total . In the ratio, the cross section measured in the τ h final state is multiplied by the branching fraction B(W → τν τ ) and divided by the inclusive tt cross section measured in the dilepton final state [14]. Since both measurements are performed in the same data-taking period with the same reconstruction algorithms, the uncertainty in the ratio includes the correlations between common sources of uncertainties as indicated in Table 2. The estimate yields the value R Γ = 0.1050 ± 0.0009 (stat) ± 0.0071 (syst), improving over the previous measurements [49,62,63]. The result is dominated by the systematic uncertainty and it is consistent with the SM value of 0.1083 ± 0.0002 [49]. While in Ref.

Summary
A measurement of the top quark pair production cross section in the tt → ( ν )(τ h ν τ )bb channel, where is either an electron or a muon, is performed by CMS in proton-proton collisions at LHC, using a data sample corresponding to an integrated luminosity of 35.9 fb −1 obtained at √ s = 13 TeV. Events are selected by requiring the presence of an electron or a muon, and at least three jets, of which at least one is b tagged and one is identified as a τ lepton decaying to hadrons (τ h ). The largest background contribution arises from tt lepton+jets events, tt → ( ν )(qq )bb, where one jet is misidentified as the τ h . The background contribution is constrained in a fit to the distribution of the transverse mass of the light lepton and missing transverse momentum system in two event categories, constructed according to the kinematic properties of the jets in the tt lepton+jets final state. The signal enters as a free parameter without constraining the kinematic properties of the τ lepton. Assuming a top quark mass of 172.5 GeV, the measured total tt cross section σ tt ( τ h ) = 781 ± 7 (stat) ± 62 (syst) ± 20 (lumi) pb is in agreement with the standard model expectation. This is the first measurement of the tt production cross section in proton-proton collisions at √ s = 13 TeV that explicitly includes hadronically decaying τ leptons, and it improves the relative precision with respect to the 7 and 8 TeV results [64,65]. The higher precision is achieved through a shape fit to the kinematic distributions of the events, thus better constraining the backgrounds. The measurement of the ratio of the cross sections in the τ h final state to the light dilepton cross section [14] yields a value of R τ h / = 0.973 ± 0.009 (stat) ± 0.066 (syst), consistent with lepton universality. The ratio of the partial to the total width of the top quark Γ(t → τν τ b)/Γ total = 0.1050 ± 0.0009 (stat) ± 0.0071 (syst) is measured with respect to the tt inclusive cross section extrapolated from the light dilepton final state, improving the precision over the previous measurements [62,63].

Acknowledgments
We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centres and personnel of the Worldwide LHC Computing Grid for delivering so effectively the computing infrastructure essential to our analyses.
[60] CMS Collaboration, "Measurement of differential cross sections for top quark pair production using the lepton+jets final state in proton-proton collisions at 13 TeV", Phys.