Search for third-generation scalar leptoquarks decaying to a top quark and a $\tau$ lepton at $\sqrt{s}=$ 13 TeV

A search for pair production of heavy scalar leptoquarks (LQs), each decaying into a top quark and a $\tau$ lepton, is presented. The search considers final states with an electron or a muon, one or two $\tau$ leptons that decayed to hadrons, and additional jets. The data were collected in 2016 in proton-proton collisions at $\sqrt{s}=$ 13 TeV with the CMS detector at the LHC, and correspond to an integrated luminosity of 35.9 fb$^{-1}$. No evidence for pair production of LQs is found. Assuming a branching fraction of unity for the decay LQ$ \to $t$\tau$, upper limits on the production cross section are set as a function of LQ mass, excluding masses below 900 GeV at 95% confidence level. These results provide the most stringent limits to date on the production of scalar LQs that decay to a top quark and a $\tau$ lepton.


Introduction
Leptoquarks (LQs) are hypothetical particles that carry non-zero baryon and lepton quantum numbers.They are charged under all standard model (SM) gauge groups, and their possible quantum numbers can be restricted by the assumption that their interactions with SM fermions are renormalizable and gauge invariant [1].The spin of an LQ state is either 0 (scalar LQ) or 1 (vector LQ).Leptoquarks appear in theories beyond the SM such as grand unified theories [2][3][4], technicolor models [5,6] and other compositeness scenarios [7,8], and R-parity-violating (RPV) supersymmetric models [9,10].
Third-generation scalar LQs (LQ 3 s) have recently received considerable theoretical interest, as their existence can explain the anomaly in the B → Dτν and B → D * τν decay rates reported by the BaBar [11,12], Belle [13][14][15], and LHCb [16] Collaborations.These decay rates deviate from the SM predictions by about four standard deviations [17], and studies of the flavor structure of LQ couplings reveal that large couplings to third-generation quarks and leptons could explain this anomaly [18][19][20][21].Third-generation LQs can appear in models in which only third-generation quarks and leptons are unified [22,23] and therefore their existence is not constrained by proton decay experiments.All models that predict LQs with masses at the TeV scale and sizable couplings to top quarks and τ leptons can be probed by the CMS experiment at the CERN LHC.In proton-proton (pp) collisions LQs are mainly pair produced through the quantum chromodynamic (QCD) quark-antiquark annihilation and gluon-gluon fusion sand t-channel subprocesses as shown in Fig. 1.There are also lepton-mediated tand u-channel contributions that depend on the unknown lepton-quark-LQ Yukawa coupling, but these contributions to LQ 3 production are negligible at the LHC as they require third-generation quarks in the initial state.Hence, the LQ pair-production cross section can be taken to depend only on the assumed values of the LQ spin and mass, and on the center-of-mass energy.The corresponding pair production cross sections have been calculated up to next-to-leading order (NLO) in perturbative QCD [24].
This paper presents the first search for the production of an LQ 3 decaying into a top quark and a τ lepton at √ s = 13 TeV.The search targets LQ 3 s with electric charges −5/3 e and −1/3 e, where e is the proton charge, and with various possible weak isospin configurations, depending on the model.A previous search for this channel at √ s = 8 TeV by the CMS Collaboration resulted in a lower mass limit of 685 GeV for an LQ 3 with branching fraction B = 1 into a top quark and a τ lepton [25].Other searches for an LQ 3 have targeted the decays LQ 3 → bν and LQ 3 → bτ [26][27][28][29][30][31][32][33][34][35][36][37][38][39].The results of the search presented here are also interpreted in the context of RPV supersymmetric models, where the supersymmetric partner of the bottom quark (bottom squark) decays into a top quark and a τ lepton via the RPV coupling.
We consider events with at least one electron or muon and at least one τ lepton, where the τ lepton undergoes a one-or three-prong hadronic decay, τ h → hadron(s) + ν τ .In LQ 3 LQ 3 events, τ leptons arise directly from LQ 3 decays, as well as from W bosons in the top quark decay chain.Electrons and muons are produced in leptonic decays of W bosons or τ leptons.Two search regions are used in this analysis: a di-τ region with the signature τ h τ h +jets and small background levels from SM processes, which provides high sensitivity for LQ 3 masses below 500 GeV, and a region with a single τ lepton in the final state, τ h +jets, which has higher sensitivity for LQ 3 masses above 500 GeV because of a larger signal efficiency.Here, denotes either an electron or a muon.The dominant backgrounds in this search come from tt+jets and W+jets production, with jets misidentified as hadronically decaying τ leptons.These backgrounds are estimated through measurements in control regions and extrapolated to the signal region.
In this paper, Section 2 describes the CMS detector, while Section 3 discusses the data samples and the properties of simulated events utilized in the analysis.Section 4 outlines the techniques used for event reconstruction and Section 5 describes the selection criteria applied in each analysis channel.The method used for the background estimation is reported in Section 6, and systematic uncertainties are detailed in Section 7. Finally, Section 8 contains the results of the analysis, and Section 9 summarizes this work.

The CMS detector
The central feature of the CMS apparatus [40] 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 (ECAL), and a brass and scintillator hadron calorimeter (HCAL), each composed of a barrel and two endcap sections.Forward calorimeters extend the pseudorapidity (η) coverage provided by the barrel and endcap detectors.Electron momenta are estimated by combining the energy measurement in the ECAL with the momentum measurement in the tracker.Muons are measured in gas-ionization detectors embedded in the steel flux-return yoke outside the solenoid.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. [40].
Events of interest are selected using a two-tiered trigger system [41], where the first level is composed of custom hardware processors and selects events at a rate of around 100 kHz within a time interval of less than 4 µs.The second level, known as the high-level trigger, uses a version of the full event reconstruction software optimized for fast processing, and reduces the event rate to around 1 kHz before data storage.

Data sample and simulated events
The search for LQ 3 s presented here uses pp collisions at √ s = 13 TeV recorded with the CMS detector in 2016.The data sample corresponds to an integrated luminosity of 35.9 fb −1 [42].
The leading order (LO) Monte Carlo (MC) program PYTHIA 8.205 [43] is used to simulate the LQ 3 pair production signal process.Both LQ 3 s are required to decay into a top quark and a τ lepton, and polarization effects from the chiralities of the top quark and the τ lepton have been neglected.The signal samples are generated for LQ 3 masses ranging from 200 to 2000 GeV.
The principal background processes, top quark pair production (tt) via the strong interaction and electroweak single top quark production in the t-channel and tW processes, are simulated with the NLO generator POWHEG (v1 is used for the single top tW processes and v2 for the single top t-channel and tt processes) [44][45][46][47][48][49].The s-channel process of single top quark production is generated at NLO using the program MADGRAPH5 aMC@NLO (v2.2.2) [50].Other background processes involve W and Z boson production in association with jet radiation.These processes are generated with MADGRAPH5 aMC@NLO (v2.2.2), with W boson production at NLO and Z boson production at LO level.The matrix element generation of W and Z boson production is matched to the parton shower emissions with the Frederix and Frixione [51] and MLM [52] algorithms, respectively.Background processes from QCD multijet production are simulated with PYTHIA 8.205.For all generated events, PYTHIA 8.205 is used for the description of the parton shower and hadronization.In the parton shower, the underlying event tune CUETP8M1 [53,54] has been applied for all samples except for tt and single top quark production in the t-channel, which use the underlying event tune CUETP8M2T4 [53,54].The event generation is performed using the NNPDF 3.0 parton distribution functions (PDFs) [55], for all events.The detector response is modeled with the GEANT4 [56] suite of programs.

Event reconstruction
Event reconstruction is based on the CMS particle-flow (PF) algorithm [57], which combines information from all subdetectors, including measurements from the tracking system, energy deposits in the ECAL and HCAL, and tracks reconstructed in the muon detectors.Based on this information, all particles in the event are reconstructed as electrons, muons, photons, charged hadrons, or neutral hadrons.
Interaction vertices are reconstructed using a deterministic annealing filtering algorithm [58,59].The reconstructed vertex with the largest value of summed physics-objects p 2 T is taken to be the primary pp interaction vertex.The physics objects are jets, clustered using the jet finding algorithm [60,61] with the tracks assigned to the vertex as inputs, and the associated missing transverse momentum, taken as the negative vector sum of the p T of those jets.Charged particles associated with other interaction vertices are removed from further consideration.
Muons are reconstructed using the information collected in the muon detectors and the inner tracking detectors, and are measured in the range |η| < 2.4.Tracks associated with muon candidates must be consistent with muons originating from the primary vertex, and are required to satisfy a set of identification requirements.Matching muon detector information to tracks measured in the silicon tracker results in a p T resolution for muons with 20 < p T < 100 GeV of 1.3-2.0% in the barrel and better than 6% in the endcaps.The p T resolution in the barrel is better than 10% for muons with p T up to 1 TeV [62].
Electron candidates are reconstructed in the range |η| < 2.5 by combining tracking information with energy deposits in the ECAL.Candidates are identified [63] using information on the spatial distribution of the shower, the track quality and the spatial match between the track and electromagnetic cluster, the fraction of total cluster energy in the HCAL, and the level of activity in the surrounding tracker and calorimeter regions.The transverse momentum p T resolution for electrons with p T ≈ 45 GeV from Z → ee decays ranges from 1.7% for nonshowering electrons in the barrel region to 4.5% for electrons showering in the endcaps [63].
Jets are clustered using PF candidates as inputs to the anti-k T algorithm [60] in the FASTJET 3.0 software package [61], using a distance parameter of 0.4.For all jets, corrections based on the jet area [64] are applied to the energy of the jets to remove the energy contributions from neutral hadrons from additional pp interactions in the same or adjacent bunch crossings (pileup collisions).Subsequent corrections are used to account for the nonlinear calorimetric response in both jet energy and mass, as a function of η and p T [65].The jet energy resolution amounts typ-Table 1: Summary of selection criteria in event categories A ( τ h + jets) and B ( τ h τ h + jets), where = µ, e.In category A, the two subcategories, OS and SS, are defined by the charge of the τ h pair.The fit variable used in each category is also shown.
T in two S T bins number of events ically to 15% at 10 GeV, 8% at 100 GeV, and 4% at 1 TeV [66].Corrections to the jet energy scale and the jet energy resolution are propagated to the determination of the missing transverse momentum [66].Jets associated with b quarks are identified using the combined secondary vertex v2 algorithm [67,68].The working point used for jet b tagging in this analysis has an efficiency of ≈65% (in tt simulated events) and a mistag rate (the rate at which light-flavor jets are incorrectly tagged) of approximately 1% [68].
Hadronically decaying τ leptons are reconstructed with the hadron-plus-strips (HPS) algorithm [69] and are denoted by τ h .The HPS algorithm is based on PF jets and additionally includes photons originating from neutral pion decays.Energy depositions in the ECAL are reconstructed in "strips" elongated in the direction of the azimuthal angle φ, to take account of interactions in the material of the detector and the axial magnetic field.These deposits are associated with one or three charged tracks to reconstruct various hadronic decay modes of τ leptons.To suppress backgrounds from light-quark or gluon jets, a τ h candidate is required to be isolated from other energy deposits in the event.The isolation criterion is based on the scalar p T sum I τ of charged and neutral PF candidates within a cone of radius (∆η) 2 + (∆φ) 2 = 0.5 around the τ h direction, excluding the τ h candidate.The isolation criterion is I τ < 1.5 GeV [70].
The energies and resolutions as well as the selection efficiencies for all reconstructed jets and leptons are studied in data and simulated events [62,63,66,68,70].Based on these studies, the simulation is corrected to match the data.

Event selection and categorization
In the online trigger system, events with an isolated muon (or electron) with p T > 24 (27) GeV and |η| < 2.4 (2.1) are selected in the muon (electron) channel.We select events offline containing exactly one isolated muon (or electron) with p T > 30 GeV and |η| < 2.4 (2.1).For the electron channel, a veto is applied to events with a muon to avoid overlap between the two channels.At least one τ h lepton with p T > 20 GeV and |η| < 2.1 and at least two jets with p T > 50 GeV and |η| < 2.4 are required.Events are selected if a third jet with p T > 30 GeV and |η| < 2.4 is present, and any additional jets are only considered if they have p T > 30 GeV.
The magnitude of the missing transverse momentum, p miss T , is required to be above 50 GeV.Further, the events are divided into two categories corresponding to the number of observed LQ candidates, allowing the sensitivity to be enhanced over a broad range of LQ masses.The event selection was chosen to maximize the expected significance of a possible LQ signal.A summary of the selection criteria for both categories is given in Table 1 and described below.

Category A: τ h + jets
In this category, exactly one τ h lepton is required in addition to the presence of one electron or muon.High p T requirements are applied to maximize the sensitivity at high LQ masses.The leading jet is required to have p T > 150 GeV.In addition we define two subcategories based on the electric charges of the particles in the τ h pair: opposite-sign (OS) and same-sign (SS).Events passing the OS τ h pair requirement must contain at least four jets and have p miss T > 100 GeV.For both subcategories, we require that the leading tau lepton has p T > 100 GeV and that there is at least one b-tagged jet.Finally the events are divided into two regions of S T , where S T is the scalar p T sum of all selected jets, leptons, and p miss T .In the low (high)-S T search regions, events must satisfy S T < 1200 (≥ 1200) GeV.This division adds sensitivity for LQ 3 masses of 600 GeV and higher.
The top quarks originating from the decay of a heavy LQ 3 are expected to be produced with larger p T than the top quarks produced in background processes.Therefore, the transverse momentum distribution of the top quark candidate decaying into hadronic jets (p t T ) gives discrimination power between background and signal events, and a measurement of the p t T spectrum is performed in category A.
A kinematic reconstruction of the top quark candidate is performed by building top quark hypotheses using between one and five jets.Because of the presence of multiple hypotheses in each event, we choose the hypothesis in which the reconstructed top quark mass is closest to the value of 172.5 GeV.
The statistical evaluation in this category is performed through a template-based fit to the measured p t T distribution.

Category B: τ h τ h + jets
In this category events are required to have at least two τ h leptons and one electron or muon.This requirement of two τ h leptons removes a large fraction of the SM background processes.The exception to this exclusion of SM backgrounds are diboson production events that contain one or more τ h leptons, but the cross sections for these processes are small.The selection criteria in this category are adapted to provide good sensitivity for low LQ masses.
Each event is required to contain an OS τ h τ h pair.If the event contains more than one τ h τ h pair, the OS pair with the largest scalar p T sum is selected.Moreover, the leading and subleading τ h must satisfy p T > 65 GeV and 35 GeV, respectively.
In this category a counting experiment is performed, as the number of expected background events is too small for results to benefit from a shape-based analysis.

Background estimation
The background in this analysis consists of samples of events that are selected because of jets misidentified as τ h leptons and events with one electron or muon together with one or more τ h leptons.
In the following, events from tt and W+jets production that contain at least one misidentified τ h lepton are obtained from control regions (CRs) separately defined for the two search regions (SRs) A and B. We consider the following contributions: the tt background that consists of only misidentified τ h leptons (or exactly one misidentified τ h lepton as in category A), denoted by tt f , the tt background that consists of (at least) one τ h lepton and (at least) one misidentified τ h  An extrapolation method is used to derive the background due to misidentified τ h leptons.The normalization, and in category A also the shape, of the tt background is estimated using where N is the total number of events for the respective process in the signal region or control region and where "other" denotes all non-tt background processes that are estimated from simulation.The contribution to the background from events with τ h leptons only is estimated from simulated events.

Backgrounds in category A
In each subcategory of category A, the largest fraction of background events originates from tt production.The second largest source of background events arises from W+jets production, while minor contributions come from single top quark and Z +jets production.
The tt f background and the W+jets background that contain a misidentified τ h lepton are derived from a single control region (CR A ), which is defined through the same selection requirements as for the SR, but with an inverted isolation requirement for the τ h lepton.
The shape of the p t T distribution is compared between the CR A and SR in simulated tt and W+jets events.Since the inversion of the τ h isolation criterion introduces kinematic differences between the SRs and CRs, the jet multiplicity and p t T are corrected in order to reproduce the shape of the tt and W+jets backgrounds in the SRs [71], as shown in Fig. 2.
Once the kinematic distributions in the CR A are corrected, we use Eq. ( 1) to extrapolate the tt and W+jets background yields to the SR.In this equation, we replace N tt with N tt, W+jets for category A.

Backgrounds in category B
In category B, the dominant background also originates from tt production.As the fraction of misidentified electrons and muons was found to be negligible in this analysis, at least one of the two τ h leptons is mimicked by a jet.Thus, background events from tt production consist either of only misidentified τ h leptons or one τ h lepton and one misidentified τ h lepton, plus an electron or a muon.A separate CR is defined for each component.The strategy for determining this background in category B is shown in Fig. 3.
The first control region (CR B1 ) is defined by inverting the isolation criterion for all τ h leptons with respect to the isolation criterion applied in the SR.The region CR B1 is used to extrapolate the tt f background to the SR.In contrast to the SR, the charge criterion on the τ h lepton is removed and the leading τ h lepton must have p T < 100 GeV to avoid overlap between the control region CR B1 and control region CR A .The tt f background normalization is then derived as in Eq. (1).
A second control region (CR B2 ) to estimate the tt p+f background is defined, in which at least one isolated and at least one nonisolated τ h lepton are required.In contrast to the SR, the charge criterion on the τ h lepton is removed and the leading τ h lepton must have p T < 45 GeV.The event must have an opposite-sign τ h pair.For this requirement, the pair with the largest summed p T is chosen.In addition, the events must satisfy M T ( , p miss T ) > 100 GeV, where M T ( , p miss T ) is the transverse mass of the lepton-p miss T system and defined as The largest non-tt p+f fraction in control region CR B2 arises from the tt f events.The estimate of this background is derived from the control region CR B1 and extrapolated to the control region CR B2 by using the extrapolation method as in Eq. ( 1).Once the tt f background is estimated from CR B1 , it is subtracted from CR B2 .The tt p+f background is extrapolated to the SR by using the extrapolation method as in Eq. (1).

Systematic uncertainties
Systematic uncertainties can affect both the overall normalization of background components, and the shapes of the p t T distributions for signal and background processes.Uncertainties in the MC simulation are applied to all simulated events used in the signal and in the various control regions.For each systematic uncertainty, the background estimation procedure described in Section 6 is repeated to study the impact of the respective systematic variation on the final result of the analysis.In the following, the systematic uncertainties applied to the analysis are summarized.
• The estimation of pileup effects is based on the total inelastic cross section.This cross section is determined to be 69.2 mb.The uncertainty is taken into account by varying the total inelastic cross section by 5% [80].
• Simulated events are corrected for lepton identification, trigger, and isolation efficiencies.The corresponding scale factors are applied as functions of |η| and p T .The systematic uncertainties due to these corrections are taken into account by varying each scale factor within its uncertainty.
• The scale factors for the jet energy scale and the jet energy resolution are determined as functions of |η| and p T [66].The effect of the uncertainties in these scale factors are considered by varying the scale factors within their uncertainties.These variations are propagated to the measurement of the p miss T .
• Scale factors for the b tagging efficiencies are applied.These scale factors are measured as a function of the jet p T [68].The corresponding uncertainty is taken into account by varying the scale factors within their uncertainties.
• Various uncertainties in the τ lepton reconstruction are considered.An uncertainty of 5% in the τ lepton identification is applied, with an additional uncertainty of 0.2 p T /(1 TeV).An uncertainty of 3% in the τ lepton energy scale is taken into account, and an uncertainty in the charge misidentification rate of 2% is applied [70].
• Parton distribution functions from the NNPDF 3.0 set are used to generate simulated events for both background and signal samples.The uncertainties in the PDFs are determined according to the procedure described in Ref. [81].The associated PDF uncertainties in the signal acceptance are estimated following the prescription for the LHC [81].
• We consider uncertainties in the renormalization (µ R ) and factorization (µ F ) scales by varying the respective scales, both simultaneously and independently, by factors between 0.5 and 2.
• We apply an uncertainty in the background estimation method by varying the extrapolation factors for background processes without τ leptons within their uncertainties.An additional uncertainty due to the correction factors used to reweight events in control region CR A is applied.
The systematic uncertainties with the largest effects on the most important background processes and on the signal are summarized in Table 2.The most important background processes are the tt f , tt f and W+jets, and tt p+f backgrounds derived from data, and the tt p background taken from simulation.Also shown is the systematic uncertainty associated with the signal produced by an LQ 3 whose mass is 700 GeV.The impact of the different sources of uncertainty varies for different processes.The uncertainty due to the variation in the scales µ R and µ F has a large impact on the tt p background, and is derived from simulation.The uncertainty in the τ lepton identification has the largest effect on the signal sample.For the backgrounds derived from several CRs, the uncertainty in the extrapolation factor has the largest impact.

Results
The results of all search categories in the electron and muon channels are combined in a binnedlikelihood fit.A statistical template-based analysis, using the measured p t T distributions in category A and a counting experiment with the events measured in category B, is performed by using the THETA software package [82].Each systematic uncertainty discussed in Sec.7 is accounted for by a nuisance parameter in the likelihood formation.
The post-fit p t T distributions in the electron and muon channels in category A are shown in Figs. 4 and 5, respectively.Contributions from tt and W+jets production with a misidentified τ h lepton are derived from control region CR A , whereas SM backgrounds with a τ h lepton and other small backgrounds are taken from simulation.
In Table 3, the total number of events from background processes and signal processes in category B is summarized.No significant deviation from the SM prediction is observed in the data in either category A or category B.
A Bayesian statistical method [82,83] is used to derive 95% confidence level (CL) upper limits on the product of the cross section and the branching fraction squared for LQ 3 pair production.Pseudo-experiments are performed to extract expected limits under a background-only hypothesis.For the signal cross section parameter, we use a uniform prior distribution.For the nuisance parameters, log-normal prior distributions are used.These are randomly varied within their ranges of validity to estimate the 68 and 95% CL expected limits.Correlations between the systematic uncertainties across all channels are taken into account.The statistical uncertainties of simulated samples are treated as an additional Poisson nuisance parameter in each bin of the p t T distribution.The 95% CL upper limits on the product of the cross section and the branching fraction squared B 2 as a function of LQ 3 mass and the 95% CL upper limits on the LQ 3 mass as a function of B are shown in Fig. 6 (left).The cross section for pair production of scalar LQs at NLO accuracy [24] is shown as the dashed line.The dotted lines indicate the uncertainty due to the PDFs and to variations of the renormalization and factorization scales by factors of 0.5 and 2.
Production cross sections of 0.6 pb for LQ 3 masses of 300 GeV and of about 0.01 pb for masses up to 1.5 TeV are excluded at 95% CL under the assumption of B = 1 for LQ 3 decays to a top  T for events in the electron channel passing the full selection in category A. The events are separated into OS (upper), SS (lower), low S T (left) and high S T (right) categories.The hatched areas represent the total uncertainties of the SM background.In the bottom panel, the ratio of data to SM background is shown together with statistical (dark gray) and total (light gray) uncertainties of the total SM background.T for events in the muon channel passing the full selection in category A. The events are separated into OS (upper), SS (lower), low S T (left) and high S T (right) categories.The hatched areas represent the total uncertainties of the SM background.In the bottom panel, the ratio of data to SM background is shown together with statistical (dark gray) and total (light gray) uncertainties of the total SM background.quark and τ lepton.Comparing these limits with the NLO cross sections, LQ 3 masses up to 900 GeV (930 GeV expected) can be excluded.
Exclusion limits with varying branching fractions B are presented in Fig. 6 (right), where limits on the complementary LQ 3 → bν (B = 0) decay channel are also included.The results for B = 0 are obtained from a search for pair-produced bottom squarks [38] with subsequent decays into b quark and neutralino pairs, in the limit of vanishing neutralino masses.Scalar LQ 3 s can be excluded for masses below 1150 GeV for B = 0 and for masses below 700 GeV over the full B range.For the assumptions of a LQ with symmetric couplings under the SM gauge symmetry and with decays to only bν and tτ, B can only take values of 1 or 0.5.When these assumptions are lifted, B can take all possible values between 0 and 1.Note that if upper limits on B are to be used to constrain the lepton-quark-LQ 3 Yukawa couplings, λ bν and λ tτ , kinematic suppression factors that favor bν decay over the tτ decay have to be considered as well [26,27].
The results presented here can be directly reinterpreted in the context of pair produced downtype squarks decaying into top quark and τ lepton pairs.Such squarks appear in RPV SUSY scenarios and correspond to LQs with B = 0.5.These squarks are excluded up to a mass of 810 GeV, and the decay mode is dominated by the RPV coupling λ 333 [84].

Summary
A search has been conducted for pair production of third-generation scalar leptoquarks (LQ 3 s) decaying into a top quark and a τ lepton.Proton-proton collision data recorded in 2016 at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 35.9 fb −1 , has 500 1000 1500 [GeV] been analyzed.The search has been carried out in the τ h +jets and τ h τ h +jets channels, where is either an electron or muon and τ h indicates a tau lepton decaying to hadrons.Standard model backgrounds due to misidentified τ h leptons are derived from control regions.The measured transverse momentum distributions for the reconstructed top quark candidate are analyzed in four search regions in the τ h +jets channel.The observed number of events are found to be in agreement with the background predictions.
Upper limits on the production cross section of LQ 3 pairs are set between 0.6 and 0.01 pb at 95% confidence level for LQ 3 masses between 300 and 1700 GeV, assuming a branching fraction of B = 1.The scalar LQ 3 s are excluded with masses below 900 GeV, for B = 1.This result represents the most stringent limits to date on LQ 3 s coupled to τ leptons and top quarks and constrains models explaining flavor anomalies in the b quark sector through contributions from scalar LQs.

Figure 1 :
Figure 1: Dominant leading order Feynman diagrams for the production of leptoquark pairs in proton-proton collisions.

Figure 2 :
Figure2: Shape comparison between the category A signal region and the corresponding control region, as a function of p t T , for simulated tt and W+jets events.Events with an opposite-sign µτ h pair are shown on the left, while those with a same-sign µτ h pair are shown on the right.The full selection is applied and the S T categories are combined.All histograms are normalized to the total number of entries.Uncertainties of the signal region and control region are indicated by red error bars and gray hatched areas, respectively.The gray band in the ratio plot corresponds to the statistical uncertainty in the simulated samples.lepton (only used in category B), denoted by tt p+f , and the tt background that consists of one τ h lepton, denoted by tt p .

Figure 3 :
Figure 3: Strategy for the background estimation in category B. The tt f background in the signal region is derived from the control region CR B1 .The tt p+f background in the signal region is derived from the control region CR B2 .To obtain an estimate of the tt f background in the control region CR B2 , the control region CR B1 is used.

Figure 4 :
Figure 4: Distributions of p tT for events in the electron channel passing the full selection in category A. The events are separated into OS (upper), SS (lower), low S T (left) and high S T (right) categories.The hatched areas represent the total uncertainties of the SM background.In the bottom panel, the ratio of data to SM background is shown together with statistical (dark gray) and total (light gray) uncertainties of the total SM background.

Figure 5 :
Figure 5: Distributions of p tT for events in the muon channel passing the full selection in category A. The events are separated into OS (upper), SS (lower), low S T (left) and high S T (right) categories.The hatched areas represent the total uncertainties of the SM background.In the bottom panel, the ratio of data to SM background is shown together with statistical (dark gray) and total (light gray) uncertainties of the total SM background.

Figure 6 :
Figure6: Upper limits at 95% confidence level on the product of the cross section and the branching fraction squared (left), and on the leptoquark mass as a function of the branching fraction (right), for the pair production of scalar LQs decaying to a top quark and a τ lepton.In the left plot, the theoretical curve corresponds to the NLO cross section with uncertainties from PDF and scale variations[24], shown by the dotted lines.The right plot additionally includes results from a search for pair-produced bottom squarks[38].

Table 2 :
Summary of largest systematic uncertainties for the tt f (and W+jets) and tt p+f backgrounds derived from data, for the tt p background obtained from simulation and for a leptoquark signal with a mass of 700 GeV.Shown are the ranges of uncertainties, which are dependent on the search regions and the lepton channel type.

Table 3 :
Final event yield in category B in the muon and electron channels for different leptoquark mass hypotheses, the background processes, and data.The total uncertainties for the signal and the background processes are shown.