Measurement of the top quark pair production cross section in dilepton final states containing one τ lepton in pp collisions at s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \sqrt{s} $$\end{document} = 13 TeV

The cross section of top quark pair production is measured in the tt¯→ℓνℓτhντbb¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathrm{t}\overline{\mathrm{t}}\to \left(\ell {\nu}_{\ell}\right)\left({\uptau}_{\mathrm{h}}{\nu}_{\uptau}\right)\mathrm{b}\overline{\mathrm{b}} $$\end{document} 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\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \sqrt{s} $$\end{document} = 13 TeV with the CMS detector. The measured cross section is σtt¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\sigma}_{\mathrm{t}\overline{\mathrm{t}}} $$\end{document} = 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¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathrm{t}\overline{\mathrm{t}} $$\end{document} production cross section in proton-proton collisions at s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \sqrt{s} $$\end{document} = 13 TeV that explicitly includes τ leptons. The ratio of the cross sections in the ℓτh and ℓℓ final states yields a value Rℓτh/ll\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {R}_{\ell {\uptau}_{\mathrm{h}}/\mathcal{ll}} $$\end{document} = 0.973 ± 0.009 (stat) ± 0.066 (syst), consistent with lepton universality.


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][10][11][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 -1 -JHEP02(2020)191 (τ ), 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.

JHEP02(2020)191
Signal tt events are simulated with the powheg event generator (v2) [17][18][19][20][21] at next-to-leading-order (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 Mad-Graph5 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) -3 -

JHEP02(2020)191
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 anti-k 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 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 -4 -JHEP02(2020)191 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 [46]. 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. [47]) 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 figure 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 m T = 2| p T || p miss T |(1 − cos ∆ϕ), 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 figure 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 as .5 GeV and m W = 80.385 GeV are, respectively, the masses of the top quark and of the W boson [48], 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. Data / Pred.  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.
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 on the signal strength is determined with this approach.
Simulation CMS (13 TeV) [GeV]  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, which amounts to approximately 5% of all background events and 17% of all signal events.

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,49], 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 -7 -JHEP02(2020)191 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 resolution (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 uncertainties related to b tagging (mistagging) efficiencies are estimated from a variety of control samples enriched in b quarks (c and light-flavour quarks) [47]; the datato-simulation scale factors for b, c, and light-flavour jets are applied to the simulation and the corresponding uncertainties are included in the fit.
The uncertainties in the JES, JER, and p miss T scales are estimated according to the prescription described in ref. [50]. The uncertainty in the JES is evaluated as a function of jet p T and η. The JES and JER uncertainties are propagated to the p miss T scale.

JHEP02(2020)191
The lepton trigger, identification, and isolation efficiencies are measured in data and simulation with a tag-and-probe method in Z → + − events [43][44][45]. The simulated events are corrected with the corresponding data-to-simulation scale factors. The uncertainties in the scale factors are included as systematic uncertainties in the measurement.
The uncertainty in the integrated luminosity is estimated to be 2.5% [51].
The pileup distribution is estimated from the measured luminosity in each bunch crossing multiplied by the average total inelastic cross section. It is used to model the pileup in simulation with an uncertainty obtained by varying the inelastic pp cross section extracted from a control region by its uncertainty of ±4.6% [49].
The measurement includes the uncertainty in the modelling of the b quark fragmentation, which covers e + e − data [52][53][54][55] at the Z pole with the Bowler-Lund [56] and Peterson [57] parametrizations, and the uncertainties in the semileptonic b-flavoured hadron branching fractions according to their measured values [48]. An uncertainty in the modelling of the p T distribution of the top quark in tt processes is included to cover the difference between the predicted and observed spectra [58][59][60][61]. The fit is sensitive to the top quark p T as it affects the shape of the m T distribution. The top quark p T variation also covers the slight trend of the τ h p T distribution.
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 [62]. The underlying event tune is varied within its uncertainties [23,62]. 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. Therefore, the measurement in the fiducial phase space is performed with the nominal acceptance and its uncertainties are only included in the extrapolation to the full phase space. 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 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.

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 figure 3, for both the eτ h and µτ h final states. A good shape agreement is observed between the data and the expected sum of signal and background distributions. 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).    [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 in the τ h and light dilepton final states are correlated, as shown in the last column.
The fiducial cross section for the production of tt events is extracted from the acceptance region of kinematic phase space defined by the selection criteria described earlier.
( 8.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:

JHEP02(2020)191
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 (eτ h ) = 789 ± 11 (stat) ± 71 (syst) ± 20 (lumi) pb, (8.7) σ tt (µτ h ) = 770 ± 8 (stat) ± 63 (syst) ± 20 (lumi) pb, (8.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 figure 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 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 . The ratio is calculated as R Γ = σ tt ( τ h )B(W → τν τ )/σ tt ( ), where 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]. The W boson branching fraction B(W → τν τ ) that is included in the signal acceptance is cancelled out in the multiplication. 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 [48,63,64]. The result is dominated by the systematic uncertainty and it is consistent with the SM value of 0.1083 ± 0.0002 [48]. 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 [65,66]. 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 section 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 [63,64].

JHEP02(2020)191
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. Finally, we acknowledge the enduring support for the construction and operation of the LHC and the CMS detector provided by the following funding agencies: BMBWF and FWF ( , 123959, 124845, 124850, 125105, 128713, 128786, and 129058 (Hungary) Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited. [44] CMS collaboration, Electron and photon performance in CMS with the full 2016 data sample, CMS-DP-2017-004, CERN, Geneva, Switzerland (2017).

References
[45] CMS collaboration, Performance of the CMS muon detector and muon reconstruction with proton-proton collisions at √ s = 13 [54] SLD collaboration, Measurement of the b quark fragmentation function in Z 0 decays, Phys.
Rev. D 65 (2002)  [59] 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. Rev. -20 -