Revealing timid pseudo-scalars with taus at the LHC

A light pseudo-scalar that is copiously produced at the LHC may still be allowed by present searches. While masses above 65 GeV are effectively covered by di-photon searches, the lower mass window can be tested by a new search for boosted di-tau resonances. We test this strategy on a set of composite Higgs models with top partial compositeness, where most models can be probed with an integrated luminosity below 300 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {fb}^{-1}$$\end{document}fb-1.


Introduction
The search for new resonances is one of the main physics goals at the LHC, with the discovery of a Higgs boson being an illustrious example [1,2]. The efforts continue, mainly focusing on high mass objects typically heavier than the Higgs itself. There are in fact few searches exploring invariant masses of two Standard Model (SM) particles below, say, 100 GeV: one notable case is the search for a di-photon resonance [3,4], mostly motivated by models that feature an extended Higgs sector, like two Higgs doublet models [5] and the next-to-minimal supersymmetric SM [6].
In this article, we focus on the LHC phenomenology of a light new scalar with a mass between 10 and 100 GeV. Generically, light new scalars are strongly constrained from electroweak precision measurements (indirectly) and from direct searches at LEP and Tevatron. At the LHC, besides the above-mentioned di-photon channel, light (pseudo)-scalars are usually searched for in the decays of the 125 GeV Higgs boson. Below roughly 10 GeV, strong bounds arise a e-mail: g.cacciapaglia@ipnl.in2p3.fr b e-mail: ferretti@chalmers.se c e-mail: flacke@ibs.re.kr d e-mail: hugo.serodio@thep.lu.se from searches related to mesons, or in experiments looking for light axion-like particles (ALPs) [7][8][9][10]. Thus, the common lore is that a new scalar, in order to escape detection, needs to be either very heavy or weakly coupled to the SM.
Note, however, that it is enough to have small couplings to electrons and to the electroweak gauge bosons in order to escape direct LEP searches and electroweak precision bounds, as well as small couplings to the Higgs to avoid the Higgs portal constraints. Couplings to gluons (and heavy quarks) are less constrained, and may lead to sizable production rates at the LHC. Candidates of this kind arise naturally in composite Higgs models that enjoy a fermion-gauge underlying description [11][12][13][14][15] providing a partial UV completion. Recent lattice results [16] have started to address the mass spectrum in a specific model [17].
In this article, we will consider this class of models to explore the 10-100 GeV mass window and show that it is, in fact, poorly tested. A timid composite pseudo-scalar (TCP) arises as the pseudo-Nambu-Goldstone boson associated with an anomaly-free U(1) global symmetry in all models of partial compositeness that enjoy a UV completion, as defined in Ref. [12]. All the possible models can be classified, and give precise predictions for the properties of the TCP candidate [15], thus mapping out a complete landscape of possibilities. We show that, while some models are already partly tested by the low mass di-photon searches, others are unconstrained. We point out that searches for boosted di-tau resonances (which could reach a lower invariant mass than the current value of 90 GeV [18,19]) give very promising signals and could be a powerful complementary probe to the di-photon channel, or even be the only way to access this class of TCPs.  [15] used as benchmarks. For the top, several possibilities arise depending on the choice of top partner representation: here, as an illustration, we take the same coupling as for lighter fermions, whose mass arise from bilinear four-fermion interactions. f a / f ψ is an estimate of the ratio between the TCP decay constant f a and the composite Higgs decay constant f ψ

Description of the models
The effective Lagrangian we consider is the SM Lagrangian augmented by the following terms, up to dimension five operators (counting powers of f a ): A pseudo-scalar a described by this general Lagrangian arises, for example, in UV completions of composite Higgs models, which were classified and studied in Refs [12,15]. In this section, we briefly summarize the main results relevant for our phenomenological study. Further background information on the models is provided in Appendix A. Within this class of models, the coupling to the SM fermions in Eq. (1) is only the first term of the expansion of the spurion coupling − m f (h) e iC f a/ f aΨ f L Ψ f R + h.c. (generating the fermions masses), which breaks explicitly the U (1) shift symmetry. A derivative coupling of the TCP to fermions of the form (∂ μ a/ f a )Ψ f γ 5 γ μ Ψ f is absent in these models since the SM fermions are neutral under the TCP U (1) charge. Although such a coupling can be obtained by using the fermion equations of motion on the leading term given in Eq. (1), the two couplings are of genuinely different origin [20], as manifested in the higher-order expansion of the spurion coupling. Starting from the complete spurion term, couplings of the Higgs to two TCPs, as well as to one TCP and Z boson, arise at loop level and are given by (see Appendix A.1 for the derivation) where we list only the effect of the log-divergence (Λ ∼ 4π f a ), g A = −g/(4 cos θ W ) is the axial coupling of the Z to tops, and κ V,t are the corrections from compositeness to the couplings of the Higgs to vectors and tops, respectively. As , our result agrees with the fact that the only non-zero contribution to the h Za coupling arises from a dimension 7 operator [21].
The couplings to gauge bosons in Eq. (1) arise as anomalous couplings if the TCP is a (SM singlet) bound state of underlying SM charged fermions. In this case, the anomaly coefficients K g,W,B are fully determined by the charges of the hyper-fermions. We refer to [15] for an extensive description of a classification of UV completions giving rise to this TCP, which yields twelve models. For the purpose of this article, the TCP dynamics in the twelve models is fully specified 1 by the numerical couplings in Table 1. Note that, due to the small TCP mass, top loops also give additional sizable contributions to the couplings to gauge bosons (not included in the table, but included in our analysis).
Our goal is to confront the TCP with the existing searches and to propose a new, more sensitive search for such object. We treat the mass m a and the decay constant f a of the TCP as free parameters. In composite Higgs UV completions, f a is related to the composite Higgs decay constant f ψ , entering in the usual alignment parameter ξ = v 2 / f 2 ψ , by a relative coefficient that was estimated in [15] and is summarized in Table 1. Since bounds on composite Higgs models require f ψ 800 GeV, f a is expected to be naturally of the order of 1 ÷ 2 TeV.

Bounds from existing searches
Since the TCP is a gauge singlet, its couplings to Z and W are induced by the anomaly and by top loops, thus they are always much smaller than those of a SM Higgs boson. Hence, bounds from all LEP searches for a light Higgs, which are based on Z associated production, are evaded. At hadron colliders the TCP is copiously produced via gluon fusion. However, only very few Tevatron or LHC two-body reso-  Table 1. The bounds arise from di-muon searches [22,23] in the low mass range, di-photon searches [3,4] in the higher one, and from the BSM decay width of the Higgs [24] below 65 GeV. We have also indicated the current bounds obtained by adapting the results in [10] in the region between 20 and 65 GeV for the two models (M9 and M10) where they are the strongest nant searches reach down to resonance masses below ∼ 100 GeV. Relevant bounds arise from Run-I ATLAS [3] and Run-II CMS [4] di-photon searches, which reach down to masses of 65 and 70 GeV, respectively, as well as ATLAS and CMS low-mass di-muon searches [22,23], reaching up to 14 GeV. The bounds on f a from these searches are shown in Fig. 1, 2 for our models. The bounds are obtained by calculating the leading order TCP production cross section following from the Lagrangian (1) with PDF set NNPDF23_nnlo_as _0119_qed [25] and a conservative k-factor of 3.3 applied [26], and using branching ratios into γ γ and μμ (computed at NLO following Ref. [9]) for the models listed in Table 1. Figures on the production cross sections and branching ratios in the twelve sample models are provided in Appendix A.2. Resonant di-tau searches reach values of the mass as low as 90 GeV [18,19], however the current bounds are never competitive with the di-photon ones in that range, mainly due to the presence of the Z-peak background. As noted in Ref. [10], a CMS search looking for boosted Z in di-jet [27] may give additional bounds above 50 GeV.
For other processes, a recent comprehensive review of the existing bounds on ALPs [9] can be directly used to confront TCPs. Firstly, the one-loop suppression (1/16π 2 ) of couplings to vector bosons in the TCP Lagrangian (1) renders bounds from vector-boson-fusion or photon-fusion production very weak. This includes Z → aγ processes and production by photon fusion in Pb-Pb ultra-peripheral collisions [28]. The up-to-date most constraining searches in the mass window between 14 and 65 GeV rely on the indirect production via Higgs portal, h → aa. As compared to the generic ALP model discussed in Ref. [9], the bounds from direct searches are weakened due to the smallness of the h → aa branching fraction following from Eq. (2), and due to the smallness of the a → γ γ and a → μμ branching fractions. Nevertheless, indirect constraints arise from the bounds on the BSM decay width of the Higgs, which is currently below 34% [24]: as shown in Fig. 1, the lower bound on f a always falls short of 1 TeV for the models under consideration. 3 For m a < 34 GeV there is also a bound from h → Za (following from Eq.(3)), however it turns out even weaker than the Higgs portal one. Associated tta production may yield a bound on TCPs: using the results of the feasibility study [30] at √ s = 14 TeV with 3 ab −1 , which focuses on a → bb in the mass range between 20 and 100 GeV, significant bounds on f a can be found only for a few models in the low mass end. Associated bba production yields weaker bounds [31]. Lastly, we should mention that the contribution of the TCP to the anomalous magnetic moment of the muon [32,33] is also small. For m a = 10 GeV and f a = 1 TeV it varies from Δa μ = −5.7 × 10 −11 for M9 to Δa μ = 2.7 × 10 −10 for M7, the current discrepancy being a exp.
μ − a SM μ = (29.3 ± 7.6) × 10 −10 . As shown, the TCP represents an example of a light pseudo-scalar which would evade all existing bounds, while being copiously produced at the LHC in gluon fusion. Searches in final states from which current bounds arise can be extended in mass range. The low-mass di-muon search [23] (performed at √ s = 7 TeV) terminated at m a = 14 GeV, but the first severe physical barrier at higher mass is the dimuon background from Drell-Yan Z production. However, a dedicated low-mass di-muon trigger and a very high invariant mass resolution would be required. 4 A recent study on inclusive di-photon cross section measurements [10] has shown how to extend the low-mass reach of di-photon searches for a generic ALP. Applying their projected reach to our models we find a nice complementarity between the di-photon channel and our proposal to use the di-tau channel to be discussed below. 5

Boosted di-tau searches as a chance to explore the TCP
As TCP decays to muons and photons have small rates, it is of interest to also look at other final states. The dom- inant TCP decay channels are gg and bb, but both have very large irreducible QCD background. 6 The next-most frequent final state is τ + τ − : sizable rates of a few % are possible and the models with the lowest rates are the ones with better di-photon reach (see Fig. 1). Compared to the dimuon channel, the branching ratios are larger by a factor of ∼ m 2 τ /m 2 μ ∼ 280. One of the main challenges for low-mass di-tau resonant searches is the trigger. The topology that we find most promising is that of a boosted TCP recoiling against an initial state radiation (ISR) jet, and then decaying into τ + τ − . The boost needs to be sufficient to allow the event to pass the highlevel trigger requirement in at least one category (jet, tau or lepton p T ) and yet leave enough observable signal. Boosted di-tau pairs have already been considered by CMS [37] in searches for heavy resonances which decay to hh, h Z, or hW . In our case, the mass of the TCP is not known, thus it may not be necessary to require a full reconstruction of the taus with subsequent increase of the systematic uncertainty associated to the procedure. Furthermore, we are interested in light TCPs with a large boost, and thus smaller separation angles between the di-tau decay products can be expected. All decay modes of the di-tau system -fully hadronic, semileptonic, and leptonic -are potentially interesting. However, for the reasons mentioned above, we focus on the opposite flavor leptonic channel, in which one τ decays to an electron and the other to a muon. One crucial issue, to be discussed more extensively in Appendix B, is related to the minimum angular separation ΔR eμ between the leptons, since the boosted tau pairs have a small separation angle. We generate the signal sample pp → a → τ + τ − for m a = 10, 20, · · · 100 GeV with up to two jets at the partonic level using MadGraph [38]. We shower and hadronize with Pythia [39] and run the fast detector simulation of Delphes [40] using the standard CMS card after removing the isolation requirement between electrons and muons. Table 2 shows the value of the signal cross section σ prod. ×BR τ τ times the efficiency expected for each of the benchmark models with f a = 1 TeV after imposing the following requirements: 7 p T μ > 52 GeV, p T e > 10 GeV, The upper cut on the separation ΔR eμ < 1 is essential in reducing the background from W 's [41], while we do not impose any minimum value yet. (This last issue is discussed below and in Appendix B.) The leading irreducible SM backgrounds are tt+ single top, γ * /Z and V V , the last one being mainly W + W − . Our simulation of these backgrounds yields σ × = 59, 25 and 11 fb, respectively, after imposing the same cuts as above. As discussed in Appendix B, we expect the reducible backgrounds of single vector boson+fakes and QCD to be subleading (of the order of a few fb) in the (μ, e) channel. Note that we do not require full reconstruction of the tau momentum.
To be able to estimate the reach of Run II + III at the LHC, it is essential to know the systematic uncertainties because we are in a situation where the signal over background ratio is small, S/B 1, as can be seen from Table 2. Since a search of this type has not been done by the experimental collaborations, we cannot reliably quantify the systematic uncertainties yet. Data driven techniques can certainly be used to reduce the systematic uncertainties on the different backgrounds. For the lepton identification, the systematic uncertainties typically amount to 1% [41], while we do not require tau identification, which would increase the systematics to 10-20%. To assess the feasibility of the analysis, we will include the systematic uncer- in the proposed di-tau search after an integrated Luminosity of 300 fb −1 . We assume a systematic error on the background δ = 1%. Shown in grey are the current bounds as of Fig. 1 tainty δ in the significance Z according the approximate formula where S and B are the number of signal and background events at a given integrated luminosity and δ the relative systematic error on the background. Eq. (4) works quite well in the regime of interest for this work when compared to the more general treatment in [42]. A projection of the bound on f a for the various models after 300 fb −1 integrated luminosity is shown in Fig. 2, including an estimated systematic uncertainty of 1%. We can see that for all models with exception of M9 and M10, the boosted di-tau search we propose can probe the mass range of 10-70 GeV with integrated luminosity below 300 fb −1 . n Fig. 3 we show the relative change in the projected bound on f a if a minimum ΔR eμ cut of 0.1 or 0.2 is imposed as well as its dependence on a change in systematic uncertainty from 1 to 0, 0.5, and 2%. These changes apply to all models in a universal way. We can see a loss of sensitivity for masses below ∼ 30 − 40 GeV when raising ΔR eμ , while above this mass range the search is barely affected. Thus, being able to remove or reduce the minimum separation angle is important for the lowest mass region, as long as it does not imply an increase in the systematic errors.
The plot also clearly shows the importance of controlling the systematic uncertainties to a level close to 1%. The latter values are what CMS and ATLAS typically require for opposite sign leptons [41] in current searches.
It is possible to improve sensitivity by imposing variable cuts on the invariant mass m eμ and particularly on the angular separation ΔR eμ of the lepton pair depending on the mass range of interest. For guidance we show in Fig. 4 the kine- in the projected bounds on f a with 300 fb −1 of data. We plot the relative change against the baseline presented in Fig. 2 for different values of systematic uncertainties δ = 0, 0.5, 1, and 2% (green, blue, black, and red) and choosing three different separation cuts ΔR eμ > 0 (solid) 0.1 (dashed) and 0.2 (dotted) respectively As mentioned above, fully or semi-hadronic decays of the di-tau system may also be testable by designing appropriate di-tau jet algorithms. For the semi-hadronic case, a sophisticated isolation procedure has been used by CMS for boosted Higgs tagging in the di-tau channel. However, large systematic uncertainties, or the order of 20-30% [37], are introduced due to the modified isolation and tau-identification procedures. Furthermore, the signal we are interested in features smaller separation between the two taus, thus a better performance may be achieved by a dedicated identification procedure. For instance, the technique of "mini-isolation" proposed in Ref. [43] may be adapted to this case, although this is beyond the scope of this paper. For the fully hadronic case, preliminary studies in Refs. [44,45] show that a good discrimination between di-tau jets and single tau or QCD jets can be achieved using sub-jet variables. However, a correct estimate of the background (especially from QCD) can only be done with data driven techniques, thus we do not attempt to quantify the sensitivity of these channels.

Conclusions
The search for new resonances at the LHC continues, and many searches for high-mass resonances are being performed. Nevertheless, complementary searches for lowermass resonances which have evaded current constraints must not be forgotten. We observe that light pseudo-scalars in the mass regime between 14 and 65 GeV can be copiously produced at the LHC while avoiding current experimental constraints.
We propose to search for boosted di-tau resonances, produced via gluon fusion, that can effectively cover this open window. We test this strategy on a set of twelve benchmark models of composite Higgs with top partial compositeness, which have a simple gauge-fermion underlying description. Low mass di-photon searches effectively cover masses above 65 GeV. Extending the di-photon search to lower masses is challenging due to triggers (but potential solutions have been presented [10]), while resuming low-mass di-muon resonant searches and extending them to higher masses is possible but challenging due to increased muon p T trigger thresholds. The boosted di-tau search we propose allows to access the open window below 65 GeV, and, for some models, it can be competitive with the di-photon channel at higher masses. IBS CTPU, where part of this work was performed. TF was supported by IBS under the project code IBS-R018-D1. GF is supported in part by a grant from the Wallenberg foundation. GC acknowledges partial support from the Institut Franco-Suedois (project Tör) and the Labex Lyon Institute of the Origins -LIO. HS has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (Grant agreement No 668679).
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .

Appendices A Details of the models
A light Higgs may emerge as the pseudo-Nambu-Goldstone boson (pNGB) of a spontaneously broken global symmetry of a new strongly interacting sector. Most of the developments in these models have been achieved by use of effective descriptions. However, as shown in [15], the knowledge of the underlying theory (a partial UV completion) based on gauge-fermion models allows for specific, and testable, predictions.
A particular minimalist class of UV completions was found in [12]. In such models, dangerous leptoquark pNGBs are avoided by having two distinct cosets, one associated with color and the other with electroweak quantum numbers. The presence of two cosets also requires the presence of two species of underlying fermions in the theory, ψ and χ . Imposing further dynamical assumptions, one is left with twelve models, whose properties are summarized in Table 3. The details of the models have been extensively explored in Ref. [15], so here we will simply recall their main features.
In order to study the low energy degrees of freedom of the present models we use the formulation of chiral perturbation theory. We parameterize relevant degrees of freedom as where r = ψ or χ . The meson matrices Σ r are associated with electroweak and color cosets. The number of fields π a r , associated with the broken generators of the non-Abelian symmetries, is model dependent. The matrix Σ 0,r is the gauge-preserving vacuum. The meson matrix Σ ψ contains a Higgs in the custodial representation on top of extra scalar multiplets. The coefficient c 5 is √ 2 when the ψ irrep is real and 1 otherwise.
While the meson matrices Σ r have a different structure for each coset, the Abelian terms Φ r characterize a universal feature of all the models, i.e. the presence of two pseudo-scalars in the spectrum. They are associated with the freedom to preform a chiral rotation in the underlying fermions ψ and χ , or, in other words, with the presence of the Abelian symmetries U (1) ψ,χ . These symmetries are spontaneously broken by the respective condensates, and explicitly broken by the underlying fermion masses and the gauging of the SM symmetries (via anomalies). On top of this, one combination has an anomaly with the new strong dynamics gauge bosons. The charges for the non-anomalous U (1) are given in Table 3. For each model we can, therefore, find the linear combinations of a ψ,χ that are anomaly-free and anomalous, i.e. Table 3 The first column shows the EW and color cosets, respectively. The −q χ /q ψ column indicates the ratio of charges of the fermions under the non anomalous U (1) combination. HC is the confining hyper-color gauge group. F and A 2 denote the fundamental and anti-symmetric representation of HC respectively. These combinations are, in general, not the physical mass eigenstates. Ref. [15] studied the phenomenology of these two states in the case where they are both heavy (above 500 GeV). In this article, we are only interested in the limit where one of them is light: this can only happen for the pNGB associated to the non-anomalous U(1), i.e.ã, while the other one is massive and decouples. Thus, the light mass eigenstate a coincides with the non-anomalousã. The couplings of the light pNGB to gauge bosons and fermions can be expressed as (we match the notation used in (1) with the one used in Ref. [15]): with A = g, W, B, and we have normalized the coefficients in such a way to render them independent on the normalization of the charges. The numerical values for each model are given in Table 1 To quantify the amount of fine-tuning, it is necessary to give an estimate of this decay constant with respect to the one entering the Higgs sector. To do so, we relate the decay constants in the abelian and non-abelian sectors by use of large-N estimates: f a ψ = N ψ f ψ and f a χ = N χ f χ . This leads to: The ratio f ψ / f χ can be estimated based on the MAC hypothesis [15], leading to for models M1, …, M12.
We finally want to comment on the couplings to fermions. For the light quarks and leptons, we assume that they couple to the strong dynamics via bi-linear four fermion interactions only involving ψ's. This is mainly due to the fact that it is impossible to generate enough partners to make all SM fermions partially composite. Thus, the couplings are obtained by setting (n ψ , n χ ) = (2, 0) in the above expressions. For the top, the coupling crucially depends on the charges of the top partners: if the bound state is made of ψψχ, the possible assignments are (n ψ , n χ ) = (±2, 0), (0, ± 2), (± 4, 2), while for ψχχ one has (n ψ , n χ ) = (± 2, 0), (0, ± 2), (2, ± 4). The values of the couplings C t for all models and all assignments are reported in Table 4. In the main text, we present results for the case (2, 0), which gives the same coupling to all SM fermions. This is only a representative case. Note that the bounds from other searches also depend on this choice.

A.1 Loop calculation for the h → aa and h → Za decays
The coupling of the Higgs to two TCPs is generated mainly by loops of top quarks (the contribution of lighter fermions being suppressed by the mass, while the gauge contribution is suppressed by the small anomalous couplings). The relevant vertices can be read off by expanding the spurion term −m t (h) e iC t a/ f aΨ t L Ψ t R + h.c. as follows: where the Higgs coupling is defined as so that κ t encodes the deviations from the SM coupling m t /v. The Lagrangian above allows for four diagrams, depicted in Fig. 5. The last three contain a quadratic divergence, that vanishes once they are summed. Thus, we are left with a log divergence that contributes to the amplitude as: where p h and p ai are the four-momenta of the Higgs and of the two TCPs respectively. The operator in Eq. (2) matches the divergent part of the above amplitude. Note also that the result differs from the one in Ref. [9] by a factor of 1/4 while having the same form: the two calculations indeed refer to two different models, as in Ref. [9] only a derivative coupling is considered. The calculation of the h → Za vertex is the same as in Ref. [21], except for the modifications of the Higgs couplings κ t and κ V that prevent the cancellation of the log divergences.

A.2 Cross sections, widths and branching ratios
In Fig. 6 we plot the TCP production cross section for gluon fusion at the LHC with a center-of-mass energy of 7, 8 and 13 TeV. We also show, in Fig. 7, the branching ratios in the main decay channels.

B Details of the simulation
As already stated in the main text, we generate a sample of signal events pp → a → τ + τ − for m a = 10, 20, · · · 100 GeV with up to two jets at the partonic level using MadGraph/ MadEvent. We shower and hadronize with Pythia and pass the resulting events through the fast detector simulation of Delphes using the standard CMS card with the modification on electron and muon isolation to be discussed below.
For the signal sample we set a p T -cut of 100 GeV on the first jet at partonic level in order to increase the efficiency of the analysis at the reconstructed level. We use MLM matching and nn23lo1 PDFs for both signal and background generation.
First, we generate the signal samples with 8 K g,eff = C τ = 1 and f a = 1 TeV with the same pNGB width of 1 GeV for all masses and re-scale the production cross-section multiplying by GeV/Γ τ τ for each model.
We then perform the cuts discussed in the main paper and below on these samples. This gives the value of the TCP production cross-section after the cuts for K g,eff = 1 and f a = 1 TeV that we denote byσ prod. × , where is the efficiency of the cut flow.
The true value of σ prod. ×BR τ τ × displayed in Table 2 for each model is obtained by multiplyingσ prod. × by K 2 g,eff × BR τ τ shown in Table 5. We do not include a k-factor for this analysis. The efficiencies of the cuts depend on m a but are independent on the type of model. Thus the expected signal S is obtained as S = σ prod. × BR τ τ × × L × (TeV 2 / f 2 a ), where L is the total luminosity.
The Z /γ * → τ τ background cross-section is estimated from the Monte Carlo. We generate this sample in exactly the same way as the signal sample, i.e. with up to two jets MLM matched and a p T cut of 100 GeV on the leading jet at partonic level. We find the cross-section after matching to be 49 pb.
For the remaining background processes we generate all the fully leptonic channels. We use the total cross-sections published by ATLAS, multiplied by the appropriate leptonic branching ratios for the W and Z (BR(W → lν) = 0.326, BR(Z → l + l − ) = 0.101, where l = e, μ, τ , assuming BR(t → W b) = 100%. We find σ tt,lep. = 82. 9  For the (e, μ) channel the QCD background and the single vector boson production+fakes background are expected to be sub-leading with respect to the irreducible backgrounds above. These can only be reliably computed by the experi-ment using data driven techniques. An order of magnitude estimate using a fake-rate of 10 −3 for j → e and 10 −4 for j → μ and (conservatively) efficiencies similar to those of the signal sample ≈ 0.003 leads us to estimate their total contribution after cuts to be at most a few fb. The tt and W t backgrounds are further reduced by imposing a b-jet veto, while the di-boson backgrounds containing a Z boson are further reduced by vetoing on a third lepton. Neither of these last two cuts has any significant effect on the signal sample.
We now would like to pick a set of cuts that maximizes the figure of merit Z defined in Eq. (4) for our channel (opposite sign, opposite flavor di-lepton channel). We need to satisfy both trigger and isolation requirements for the leptons. From the trigger menus discussed in [46] we chose to retain only muons with p T μ > 52 GeV (off-line selection). This allows us to go rather low in the selection of p T j of the leading jet, on which we are not triggering. We find that p T j > 150 GeV gives an acceptable compromise between the signal rate an the reach in Z. The p T e of the electron is chosen to be above 10 GeV and we use the standard isolation requirements  between leptons and jets ΔR l j > 0.5. Furthermore, we put a third lepton veto and a b-jet veto in order to reduce the top and di-boson background. The separation between the electron and the muon requires a more detailed discussion. Figure 4 shows the ΔR eμ distribution of the signal (for two reference masses) and the main backgrounds before imposing ΔR eμ cuts. The background distributions of di-tops and di-bosons are rather flat while the signal -especially for light TCPs -is characterized by small ΔR eμ . Thus, the sensitivity of the search is increased by imposing an upper cut which we choose at ΔR eμ < 1.0.
In current searches, the ATLAS and CMS collaborations typically also use a lower cut of ΔR eμ > 0.1 or 0.2 as lepton isolation requirement. For these values, efficiencies and systematic errors are known and determined from data driven methods. Imposing an analogous isolation cut thus gives a reliable estimate of systematics and the reach of the study, but it comes at a price because for a light TCP, the signal efficiency is reduced. To study the impact of lepton isolation in more detail we explore three cases: without μ-e isolation, with ΔR eμ > 0.1, and 0.2. The resulting values for σ prod. × BR τ τ × for the case of no isolation are given in Table 2. The analogous values with an isolation criterion ΔR eμ > 0.1 and 0.2 are shown in Tables 6 and 7, respectively.
The background cross sections are barely affected by the modified isolation cuts. They are given in the caption of the respective table. The signal cross sections for the different isolations are comparable for a high mass TCP, but get reduced for m a 20 GeV (30 GeV) for ΔR eμ > 0.1 (0.2). Thus removing or at least reducing the isolation cut is advantageous for the TCP low mass regime if systematic errors can be kept under control.
Additionally, one could consider other channels such as (τ h , μ) and (τ h , e). Perhaps even the (μ, μ) channel could be relevant in spite of the large Z /γ * background. We refrain from doing this since the systematic error and the backgrounds become harder to estimate with our tools.

C Complementarity of τ τ and γ γ searches
To conclude, we show in Fig. 8 a comparison of the reach of the di-photon analysis proposed in [10] to the di-tau search proposed in this paper. We see that there is a nice complementarity between the two approaches, the di-tau being more sensitive in the low mass region and the di-photon nicely covering the high mass regions left uncovered. The complementarity also extends to the models themselves, namely some models, like M5 and M8, are much more sensitive to the di-tau signal, while others like M9 and M10 are covered by the di-photon analysis. A combination of the two approaches would essentially allow to test all of these models in the full mass range.