Searches for lepton-ﬂavour-violating decays of the Higgs boson in √ s = 13 TeV pp collisions with the ATLAS detector

This Letter presents direct searches for lepton ﬂavour violation in Higgs boson decays, H → e τ and H → μτ , performed with the ATLAS detector at the LHC. The searches are based on a data sample of proton–proton collisions at a centre-of-mass energy √ s = 13 TeV, corresponding to an integrated luminosity of 36 . 1 fb − 1 . No signiﬁcant excess is observed above the expected background from Standard Model processes. The observed (median expected) 95% conﬁdence-level upper limits on the lepton-ﬂavour-violating branching ratios are 0 . 47% (0 . 34 + 0 . 13 − 0 . 10 %) and 0 . 28% (0 . 37 + 0 . 14 − 0 . 10 %) for H → e τ and H → μτ , respectively. SCOAP 3 .


Introduction
The search for processes beyond the Standard Model (SM) is one of the main goals of the Large Hadron Collider (LHC) programme at CERN.A possible sign of such processes is lepton flavour violation (LFV) in decays of the Higgs boson [1,2].Many beyond-SM theories predict LFV decays of the Higgs boson, such as supersymmetry [3,4], other models with more than one Higgs doublet [5,6], composite Higgs models [7], models with flavour symmetries [8] or warped extra dimensions [9][10][11] models and others [12,13].
In this Letter, searches for LFV decays of the Higgs boson, H → eτ and H → µτ, at the LHC with the ATLAS experiment are presented.Studies are based on proton-proton (pp) collision data recorded in 2015-2016 at a centre-of-mass energy √ s = 13 TeV.The dataset corresponds to an integrated luminosity of 36.1 fb −1 .
The searches presented here involve both leptonic (τ → ν ν 1) and hadronic (τ → hadrons + ν) decays of τ-leptons, denoted τ and τ had respectively.The dilepton final state τ only considers pairs of different-flavour leptons.Same-flavour lepton pairs are rejected due to the large lepton pair-production Drell-Yan background.Two channels are considered for each of the two searches: eτ µ and eτ had for the H → eτ search, µτ e and µτ had for the H → µτ search.The analysis is designed such that any potential LFV signal overlap between the H → eτ and H → µτ searches is negligible.Many methods are reused from the measurement of the Higgs boson cross-section in the H → ττ final state [17].

Simulation samples
Samples of Monte Carlo (MC) simulated events are used to optimize the event selection, and to model the signal and several of the background processes.The samples were produced with the ATLAS simulation infrastructure [21] using the full detector simulation performed by the G 4 [22] toolkit.The Higgs boson mass was set to m H = 125 GeV [23].The four leading Higgs boson production mechanisms are considered: the gluon-gluon fusion (ggF), vector-boson fusion (VBF) and two associated production modes (W H, Z H), while the others give negligible contributions and are ignored.The cross-sections of all Higgs boson production processes were normalized to the SM predictions [24].The LFV Higgs boson decays as well as the H → ττ and H → WW background decays were modelled with Pythia 8 [25].Other background processes involve electroweak production of W/Z bosons via VBF, Drell-Yan production of W/Z in association with jet(s) as well as diboson, single top-quark and top-quark pair (t t) production.The MC generators used for the SM H → ττ cross-section measurement [17] were also employed here for all background components.The generators and parton shower models used to simulate different processes are summarized in Table 1.
Table 1: Generators used to describe the signal and background processes, parton distribution function (PDF) sets for the hard process, and models used for parton showering, hadronization and the underlying event (UEPS).The orders of the total cross-sections used to normalize the events are also given.More details are given in Ref. [17].

Object reconstruction
The correct identification of H → τ events requires reconstruction of several different objects (electrons, muons, and jets, including those initiated by hadronic decays of τ-leptons) and the missing transverse momentum ì p miss T , whose magnitude is called E miss T .Electrons are reconstructed by matching tracks in the inner detector to clustered energy deposits in the electromagnetic calorimeter [57].Loose likelihood-based identification [58], p T > 15 GeV and fiducial volume requirements (|η| < 2.47, excluding the transition region between the barrel and the endcap calorimeters 1.37 < |η| < 1.52) are applied.Medium identification, corresponding to an efficiency of 87% at p T = 20 GeV, is imposed for the baseline electron selection.
Muons are identified by tracks reconstructed in the inner detector and matched to tracks in the MS.Loose identification [59], p T > 10 GeV and |η| < 2.5 requirements are applied.Medium identification (efficiency of 96.1% for muons with p T > 20 GeV) is imposed for the baseline muon selection.
Isolation criteria exploiting calorimeter and track-based information are applied to both electrons and muons.The gradient working point is used, featuring an efficiency of 90% (99%) obtained for leptons with p T > 25 GeV (60 GeV) originating from the Z → process [58,59].
Jets are reconstructed using the anti-k t algorithm [60] as implemented by the FastJet [61] package.The algorithm is applied to topological clusters of calorimeter cells [62] with a radius parameter R = 0.4.Only jets with p T > 20 GeV and |η| < 4.5 are considered.Jets from other pp interactions in the same and neighbouring bunch crossings (pile-up) are suppressed using jet vertex tagger (JVT) algorithms [63,64].Jets containing b-hadrons (b-jets) are identified by the MV2c20 algorithm [65,66] in the central region (|η| < 2.4).A working point corresponding to 85% average efficiency determined for b-jets in t t simulated events is chosen, rejection factors are 2.8 and 28 against c-jets and light-flavour jets respectively.
The reconstruction of the object formed by the visible products of the τ had decay (τ had-vis ) begins from jets reconstructed by the anti-k t jet algorithm with a radius parameter R = 0.4.Information from the inner detector tracks associated with the energy deposits in the calorimeter is incorporated in the reconstruction.Only τ had-vis candidates with p T > 20 GeV and |η| < 2.5 are considered.3One or three associated tracks with an absolute total charge |q| = 1 are required.An identification algorithm [67,68] based on boosted decision trees (BDT) [69][70][71] is used to reject τ had-vis candidates arising from misidentification of jets or from decays of hadrons with bor c-quark content.Unless otherwise indicated, a tight identification (ID) working point is used for the τ had-vis , corresponding to an efficiency of 60% (45%) for 1-prong (3-prong) candidates.Jets corresponding to identified τ had-vis candidates are removed from the jet collection.The τ had-vis candidates with one track overlapping with an electron candidate with high ID score, as determined by a multivariate (MVA) approach, are rejected.Leptonic τ-decays are reconstructed as electrons or muons.
Events considered in the analysis are triggered with single-electron or single-muon triggers.The p T thresholds depend on the isolation requirement and data-taking period [72,73].The lowest trigger thresholds correspond to 25 − 27 GeV (electrons) and 21 − 27 GeV (muons).

Event selection and categorization
Events selected in the τ channel contain exactly one electron and one muon of opposite-sign (OS) charges.Similarly in the τ had channel, a lepton and a τ had-vis of OS charges are required, and events with more than one baseline lepton are rejected.The selection criteria are summarized in Table 2 for the analysis categories as well as the control regions (CRs), which are described in Section 5.
In the τ channel, 1 and 2 denote the leading and subleading lepton in p T , respectively.Events where the leading lepton is an electron (muon) are used in the search for H → eτ µ (H → µτ e ).A requirement on the dilepton invariant mass, equal to the invariant mass of the lepton and the visible τ-decay products, m vis , reduces backgrounds with top quarks, and the criterion applied to the track-to-cluster p T ratio of the electron reduces the Z → µµ background where a muon deposits a large amount of energy in the electromagnetic calorimeter and is misidentified as an electron in the µτ e channel.The contribution from the H → ττ decay is reduced by the asymmetric p T selection of the two leptons.
In the τ had channel, the criterion based on the azimuthal separations of lepton-E miss T and τ had-vis -E miss T , i= ,τ had-vis cos ∆φ(i, E miss T ), reduces the W+jets background whereas the requirement on |∆η( , τ had-vis )| reduces backgrounds with misidentified τ had-vis candidates.
For both channels of each search, a b-veto requirement reduces the single-top-quark and t t backgrounds.Events are further categorized into VBF (with a focus on the VBF production of the Higgs boson) and non-VBF categories.The VBF selection is based on the kinematics of the two jets with the highest p T , where j 1 and j 2 denote the leading and subleading jet in p T , respectively.The variables m(j 1 , j 2 ) and ∆η(j 1 , j 2 ) stand for the invariant mass and η separation of these two jets.The non-VBF category contains events failing the VBF selection.In the dilepton channel, additional selection criteria are applied to further reject background events in this category.These criteria are also listed in Table 2, where m T stands for the transverse mass4 of the two objects listed in parentheses, and p τ T represents the magnitude of the vector  T and E miss T .The requirement on p τ T /p 1 T reduces the background arising from jets misidentified as leptons.The VBF and non-VBF categories in each of the τ and τ had channels give rise to four signal regions in each search.
The analysis exploits BDT algorithms to enhance the signal separation from the background in the individual searches, channels and categories.The components of the four-momenta of the analysis objects as well as derived event variables (e.g.invariant masses and angular separations) are the input variables of the BDT discriminant.Correlations between these input variables have been carefully checked, highly correlated variables have been removed and the remaining ones are ranked according to their discrimination power [74,75].The list of variables is then optimized, removing the lowest-ranked variables with marginal contribution to the sensitivity.The final list of variables is presented in Table 3 for each channel and category.The invariant mass of the Higgs boson reconstructed under the H → τ decay hypothesis exhibits the highest signal-to-background separation power and it helps to distinguish LFV signal from H → ττ and H → WW backgrounds.For the τ channel the invariant mass is reconstructed with the MMC algorithm [76] and is denoted by m MMC ; for the τ had channel it is reconstructed with the collinear approximation [76] and is denoted by m coll .The analysis CRs are used to validate the level of agreement between data and simulated distributions of the BDT score and input variables, as well as their correlations.
Table 3: BDT input variables used in the analysis.For each channel and category, used input variables are marked with HR (indicating the five variables with the highest rank) or a bullet.Analogous variables between the two channels are listed on the same line.

Background modelling
The most significant backgrounds in the search are from events with Z → ττ decays or with (single or pair-produced) top quarks, especially in the τ channel, as well as from events with misidentified objects, which are estimated using data-driven (d.d.) techniques.The relative contribution from misidentified objects to the total background yield is 5-25% in the τ channel and 25-45% in the τ had channel, depending on the search and the analysis category.The shapes of distributions from the Z → ττ and top-quark (single-top-quark and t t) processes are modelled by simulation in both the τ and τ had decay channels.In the τ channel, the relative contributions of Z → ττ and top-quark production processes are 20-35% and 20-55%, respectively; the top-quark background dominates in the VBF category.In the τ had channel, the top-quark background fraction is 1-10%, while the Z → ττ process contributes to 45-55% of the total background.The individual contributions are listed in Tables 4 and 5. Smaller background components are also modelled by simulation and are grouped together: Z → µµ, diboson production, H → ττ and H → WW.
Good modelling of the background is demonstrated in Figure 1 for a selection of important BDT input variables.Details of the background estimation techniques are given below.

τ channel
Two sets of CRs, as defined in Table 2, are used to constrain the normalization of Z → ττ and top-quark background components.These CRs inherit their definitions from the corresponding analysis category but invert one requirement to ensure orthogonality with the nominal selection.The normalization factors are determined during the statistical analysis by fitting the event yields in all signal and control regions simultaneously.For each search, separate Z → ττ normalization factors are used for the VBF and non-VBF categories.In the case of the top-quark background, in which leading jets are produced at a lower order of the perturbative expansion of the scattering process, a combined normalization factor across the two categories is used in the τ channel.
Top-quark CRs are almost exclusively composed of top-quark backgrounds: the purity is 95% across both searches and categories, with t t process accounting for more than 90% of the top-quark backgrounds.The Z → ττ CRs achieved a purity of ∼80% in the non-VBF categories, while a lower purity of ∼60% is observed in the VBF categories.The contributions of all other background components are normalized to their SM predictions when the likelihood fit (Section 7) is applied.
The shape and normalization of diboson and Z → µµ background distributions are validated with data in dedicated regions where their contributions are enhanced.The latter process only contributes sizeably in the µτ e channel, where it represents up to 10% of the total background.Another source of background comes from W+jets, top-quark and multi-jet events, where jets are misidentified as leptons.This background is estimated directly from OS data events where an inverted isolation requirement is imposed on the subleading lepton [17].Normalization factors are applied to correct for the inverted isolation requirement.The normalization factors are derived in a dedicated region where the leptons are required to have same-sign (SS) charges.Additional corrections are made by reweighting the MC distributions of ∆φ( 1 , E miss T ) and ∆φ( 2 , E miss T ) to data in the SS region, which improves the modelling of azimuthal angles between leptons and the E miss T direction as well as the modelling of m T ( 2 , E miss T ).A similar improvement is observed in the nominal OS region.In most of the cases, the misidentified jet mimics the lepton of lower p T , 2 , while the fraction of events where both leptons are misidentified varies between 2% to 8% across categories.The systematic uncertainties of the estimation of the misidentified lepton background include contributions from closure tests in SS and OS regions enriched with misidentified leptons, from the corrections made to the ∆φ distributions, and from the composition of the misidentified lepton background.

τ had channel
The main background contributions come from the Z → ττ process and events where either a jet or an electron is misidentified as τ had-vis .The shape of the Z → ττ background distribution is modelled by simulation, and the corresponding normalization factors are determined from the simultaneous fit of the event yields in all signal and control regions.The Z → ττ normalization factors are fully correlated with those of the τ channel, in each VBF and non-VBF category.Top-quark production represents less than 1% of the total background in the τ had channel and is determined by simulation, including its normalization, which is kept fixed in the fit.
The main contributions to jets misidentified as τ had-vis come from multi-jet events and W-boson production in association with jets, and a fake-factor method is used to estimate the contribution of each component separately.A fake factor is defined as the ratio of the number of events where the highest-p T jet is identified as a tight τ had-vis candidate to the number of events where the highest-p T jet fails to satisfy this τ-ID criterion but satisfies a looser criterion.The procedure, including systematic uncertainties, is described in Ref. [17].Since a different τ-ID working point is considered in this analysis, fake factors are re-derived as a function of p T and track multiplicity of the τ had-vis candidate.
Electrons misidentified as τ had-vis , denoted by "Z → ee (d.d.)" in the following figures and tables, represent another background component in the eτ had channel, with a contribution about five times smaller than that of jets misidentified as τ had-vis .While the rate of electrons misidentified as 3-prong τ had-vis makes a negligible contribution and is modelled by simulation, the rate of electrons misidentified as 1-prong τ had-vis is determined with a fake-factor method.This time, the fake factor is defined as the ratio of the number of events with tight τ-ID to the number of events with anti-identified τ had-vis (such a candidate satisfies all criteria but the requirement on the high electron ID score is inverted).These fake factors are derived in a dedicated Z → ee enriched region defined by |m vis − m Z | < 5 GeV, m T ( , E miss T ) < 40 GeV, and m T (τ had-vis , E miss T ) < 60 GeV, where the τ had-vis candidate satisfies the medium τ-ID (corresponding to an efficiency of 55% and 40% for 1-prong and 3-prong candidates, respectively) but not the tight τ-ID criterion to avoid overlap with the τ had signal region.These fake factors are applied to signal-like events with the anti-identified τ had-vis to determine the background contribution in the categories of the analysis.The systematic uncertainties include the statistical uncertainty of the fake factors and account for looser τ-ID in the Z → ee enriched region as well as for the subtraction of the not misidentified components in this region.

Systematic uncertainties
The systematic uncertainties affect the normalization of signal and background, and/or the shape of their corresponding final discriminant distributions.Each source of systematic uncertainty is considered to be uncorrelated with the other sources.The effect of each systematic uncertainty is fully considered in each category, including control regions.Correlations of each systematic uncertainty are maintained across processes, channels, categories and regions.The size of the systematic uncertainties and their impact on the fitted branching ratio are discussed in Section 7. The main sources of systematic uncertainties are related to the estimation of the backgrounds originating from mis-identified leptons/jets and to the jet energy scale uncertainties.
Experimental uncertainties include those originating from the reconstruction, identification, tagging and triggering efficiencies of all physics objects as well as their momentum scale and resolution.These include effects from leptons [57][58][59], τ had-vis [68], jets [63,64,77] and E miss T [78].Uncertainties affecting the kinematics of the physics objects are propagated to the BDT input variables.The corresponding shape and normalization variations of the BDT discriminant are considered in the statistical analysis.Uncertainties of the luminosity measurement [79], pile-up modelling and uncertainties specific to mis-identified background estimation techniques mentioned in Section 5 are included.
The procedures to estimate the uncertainty of the Higgs boson production cross-sections follow the recommendations of the LHC Higgs Cross-Section Working Group [80].Theoretical uncertainties affecting the ggF signal originate from nine sources [24].Two sources account for yield uncertainties, which are evaluated by an overall variation of all relevant scales and are correlated across all bins of the BDT discriminant distribution [81].Another two sources account for migration uncertainties of zero to one jet and one to at least two jets in the event [81][82][83], two for Higgs boson p T shape uncertainties, one for the treatment of the top-quark mass in the loop corrections, and two for the acceptance uncertainties of ggF production in the VBF phase space from selecting exactly two and at least three jets, respectively [84,85].For VBF and W H, Z H production cross-sections, the uncertainties due to missing higher-order QCD corrections are estimated by varying the factorization and renormalization scales up and down by factors of two around the nominal scale.For all signal samples, PDF uncertainties are estimated using 30 eigenvector variations and two α s variations using the default PDF set PDF4LHC15 [32].Uncertainties related to the simulation of the underlying event, hadronization and parton shower are estimated by comparing the acceptances when using Pythia 8.212 [25] or Herwig 7.0.3[86,87].
The sources of modelling uncertainties considered for the Z → ττ process are the same as in Ref. [17] and their effect on the event migrations between categories and on the shape of the BDT discriminant are considered, since the overall normalizations are determined from data in the statistical analysis.These systematic uncertainties include variations of PDF sets, factorization and renormalization scales, CKKW matching [88], resummation scale and parton shower modelling.The other background processes are either normalized using data (processes with top-quarks and mis-identified leptons and τ had-vis candidates) or their cross-section uncertainties have negligible impact and therefore are not included.The shape uncertainties of these backgrounds originate from experimental uncertainties only.

Statistical analysis
The searches for H → eτ and H → µτ are treated independently.For each search, the analysis exploits the four signal regions and the two control regions specified in Table 2.The BDT score distributions of all signal regions are analysed to test the presence of a signal, simultaneously with the event yields in control regions, which are included to constrain the normalizations of the major backgrounds estimated from simulation.The statistical analysis uses a binned likelihood function L(µ, θ), constructed as a product of Poisson probability terms over all bins considered in the search.This function depends on the parameter µ, defined as the branching ratio B(H → τ), and a set of nuisance parameters θ that encode the effect of systematic uncertainties in the signal and background expectations.All nuisance parameters are implemented in the likelihood function as Gaussian or log-normal constraints.The normalization factors of the single-top-quark and t t backgrounds in the τ channel and of the Z → ττ background component are unconstrained parameters of the fit.Estimates of the parameters of interest are calculated with the profile-likelihood-ratio test statistic qµ [89], and the upper limits on the branching ratios are derived by using qµ and the CL S method [90].
The discriminant distributions after the fit in each channel are shown in Figures 2 and 3. Good agreement between data and the background expectation is observed.The event yields after the background-only fit are summarized in Tables 4 and 5.In the non-VBF category, the yields in the τ had channel are larger than in the τ channel due to the looser selection criteria defined for the former channel (Section 4).Table 6 shows a summary of the uncertainties of B(H → τ).The uncertainties associated with misidentified leptons and τ had-vis candidated and those related to the jet energy scale and resolution exhibit the highest impact on the best-fit branching ratios in both searches.The combined impact from all systematic uncertainties and the data statistics ranges from 0.17% to 0.19%.Table 4: Event yields and predictions as determined by the background-only fit in different signal regions of the H → eτ analysis.Uncertainties include both the statistical and systematic contributions."Other" contains diboson, Z → , H → ττ and H → WW background processes.For the eτ had channel the "Z → ee (d.d.)" component corresponds to electrons misidentified as τ had-vis .This contribution is summed with "Other" since there are few events in the VBF category.The uncertainty of the total background includes all correlations between channels.The normalizations of top-quark ( τ channel only) and Z → ττ background components are determined by the fit, while the expected signal event yields are given for B(H → eτ) = 1%.

Results
The best-fit branching ratios and upper limits are computed while assuming B(H → µτ) = 0 for the H → eτ search and B(H → eτ) = 0 for the H → µτ search.The best-fit values of the LFV Higgs boson branching ratios are equal to (0.15 +0.18 −0.17 )% and (−0.22 ± 0.19)% for the H → eτ and H → µτ search, respectively.In the absence of a significant excess, upper limits on the LFV branching ratios are set for a Higgs boson with m H = 125 GeV.The observed (median expected) 95% CL upper limits are 0.47% (0.34 +0.13  −0.10 %) and 0.28% (0.37 +0.14 −0.10 %) for the H → eτ and H → µτ searches, respectively.These limits are significantly lower than the corresponding Run 1 limits of Refs.[14,15].The breakdown of contributions from different signal regions is shown in Figure 4.The branching ratio of the LFV Higgs boson decay is related to the non-diagonal Yukawa coupling matrix elements [91] by the formula where Γ H (SM) = 4.07 MeV [92] stands for the Higgs boson width as predicted by the Standard Model.Thus, the observed limits on the branching ratio correspond to the following limits on the coupling matrix elements: |Y τe | 2 + |Y eτ | 2 < 0.0020, and |Y τµ | 2 + |Y µτ | 2 < 0.0015.Figure 5 shows the limits on the individual coupling matrix elements Y τ and Y τ together with the limits from the ATLAS Run 1 analysis and from τ → γ searches [91,93].Also indicated are limits corresponding to different branching ratios (0.01%, 0.1%, 1%, 10% and 50%) and the naturalness limit (denoted n.l.) |Y τ Y τ | m τ m v 2 [91] where v is the vacuum expectation value of the Higgs field.

Conclusions
Direct searches for the decays H → eτ and H → µτ are performed with proton-proton collisions recorded by the ATLAS detector at the LHC corresponding to an integrated luminosity of 36.1 fb −1 at a centre-of-mass energy of √ s = 13 TeV.No significant excess is observed above the expected background from Standard Model processes.The observed (expected) upper limits at 95% confidence level on the branching ratios of H → eτ and H → µτ are 0.47% (0.34 +0.13  −0.10 %) and 0.28% (0.37 +0.14 −0.10 %), respectively.These limits are more stringent by a factor of 2 (5) quark CR inverted b-veto: VBF and non-VBF ≥ 1 b-tagged jet (p T > 25 GeV and |η| < 2.4) Z → ττ CR inverted p 1 T requirement: VBF and non-VBF 35 GeV < p 1 T < 45 GeV sum of p 2

Figure 1 :
Figure 1: Distributions of representative kinematic quantities for different searches, channels and categories, before the fit as described in Section 7 is applied.Top row: transverse mass m T ( 1 , E miss T ) (eτ µ non-VBF), collinear mass m coll (eτ had non-VBF) and m MMC (eτ µ VBF).Bottom row: m MMC (µτ e non-VBF), muon p T (µτ had non-VBF) and m coll (µτ had VBF).Entries with values that would exceed the x-axis range are included in the last bin of each distribution.The size of the combined statistical, experimental and theoretical uncertainties in the background is indicated by the hatched bands.The H → eτ (H → µτ) signal overlaid in top (bottom) plots assumes B(H → τ) = 1% and is enhanced by a factor 10. In the data/background prediction ratio plots, points outside the displayed y-axis range are shown by arrows.

Figure 2 :
Figure 2: Distributions of the BDT score after the background+signal fit in each signal region of the eτ search, with the LFV signal overlaid, normalized with B(H → eτ) = 1% and enhanced by a factor 10 for visibility.The top and bottom plots display eτ µ and eτ had BDT scores respectively, the left (right) column corresponds to the non-VBF (VBF) category.The size of the combined statistical, experimental and theoretical uncertainties of the background is indicated by the hatched bands.The binning is shown as in the statistical analysis.

Figure 3 :
Figure 3: Distributions of the BDT score after the background+signal fit in each signal region of the µτ search, with the LFV signal overlaid, normalized with B(H → µτ) = 1% and enhanced by a factor 10 for visibility.The top and bottom plots display µτ e and µτ had BDT scores respectively, the left (right) column corresponds to the non-VBF (VBF) category.The size of the combined statistical, experimental and theoretical uncertainties of the background is indicated by the hatched bands.The binning is shown as in the statistical analysis.In the data/background prediction ratio plots, points outside the displayed y-axis range are shown by arrows.

Figure 4 :
Figure4: Upper limits at 95% CL on the LFV branching ratios of the Higgs boson, H → eτ (left) and H → µτ (right), indicated by solid and dashed lines.Best-fit values of the branching ratios ( μ) are also given, in %.The limits are computed while assuming that either B(H → µτ) = 0 (left) or B(H → eτ) = 0 (right).First, the results of the fits are shown, when only the data of an individual channel or of an individual category are used; in these cases the signal and control regions from all other channels/categories are removed from the fit.These results are finally compared with the full fit displayed in the last row.

Figure 5 :
Figure5: Upper limits on the absolute value of the couplings Y τ and Y τ together with the limits from the ATLAS Run 1 analysis (light grey line) and the most stringent indirect limits from τ → γ searches (dark purple region).Also indicated are limits corresponding to different branching ratios (0.01%, 0.1%, 1%, 10% and 50%) and the naturalness limit (denoted n.l.) |Y τ Y τ | m τ m v 2[91] where v is the vacuum expectation value of the Higgs field.

Table 2 :
Baseline event selection and further categorization for the τ and τ had channels.The same criteria are also used for the control region (CR) definitions in the τ channel (Section 5), but one requirement of the baseline selection is inverted to achieve orthogonal event selection.There is no CR in the τ had channel.

Table 5 :
Event yields and predictions as determined by the background-only fit in different signal regions of the H → µτ analysis.Uncertainties include both the statistical and systematic contributions."Other" contains diboson, Z → , H → ττ and H → WW background processes.The uncertainty of the total background includes all correlations between channels.The normalizations of top-quark ( τ channel only) and Z → ττ background components are determined by the fit, while the expected signal event yields are given for B(H → µτ) = 1%.

Table 6 :
Summary of the systematic uncertainties and their impact on the best-fit value of B in the H → eτ and H → µτ searches.The measured values are obtained by the fit to data, while the expected values are determined by the fit to a background-only sample.
than the corresponding limits for the H → eτ (H → µτ) decay determined by ATLAS at √ s = 8 TeV.STFC, United Kingdom; DOE and NSF, United States of America.In addition, individual groups and members have received support from BCKDF, CANARIE, CRC and Compute Canada, Canada; COST, ERC, ERDF, Horizon 2020, and Marie Skłodowska-Curie Actions, European Union; Investissements d' Avenir Labex and Idex, ANR, France; DFG and AvH Foundation, Germany; Herakleitos, Thales and Aristeia programmes co-financed by EU-ESF and the Greek NSRF, Greece; BSF-NSF and GIF, Israel; CERCA Programme Generalitat de Catalunya, Spain; The Royal Society and Leverhulme Trust, United Kingdom.The crucial computing support from all WLCG partners is acknowledged gratefully, in particular from CERN, the ATLAS Tier-1 facilities at TRIUMF (Canada), NDGF (Denmark, Norway, Sweden), CC-IN2P3 (France), KIT/GridKA (Germany), INFN-CNAF (Italy), NL-T1 (Netherlands), PIC (Spain), ASGC (Taiwan), RAL (UK) and BNL (USA), the Tier-2 facilities worldwide and large non-WLCG resource providers.Major contributors of computing resources are listed in Ref. [94].Also at Joint Institute for Nuclear Research, Dubna; Russia.ai Also at LAL, Université Paris-Sud, CNRS/IN2P3, Université Paris-Saclay, Orsay; France.a j Also at Louisiana Tech University, Ruston LA; United States of America.ak Also at LPNHE, Sorbonne Université, Paris Diderot Sorbonne Paris Cité, CNRS/IN2P3, Paris; France.al Also at Manhattan College, New York NY; United States of America.am Also at Moscow Institute of Physics and Technology State University, Dolgoprudny; Russia.an Also at National Research Nuclear University MEPhI, Moscow; Russia.ao Also at Physics Department, An-Najah National University, Nablus; Palestine.ap Also at Physics Dept, University of South Africa, Pretoria; South Africa.aq Also at Physikalisches Institut, Albert-Ludwigs-Universität Freiburg, Freiburg; Germany.ar Also at School of Physics, Sun Yat-sen University, Guangzhou; China.as Also at The City College of New York, New York NY; United States of America.at Also at The Collaborative Innovation Center of Quantum Matter (CICQM), Beijing; China.au Also at Tomsk State University, Tomsk, and Moscow Institute of Physics and Technology State University, Dolgoprudny; Russia.
ah av Also at TRIUMF, Vancouver BC; Canada.aw Also at Universita di Napoli Parthenope, Napoli; Italy.* Deceased