Same-Sign Dilepton Excesses and Vector-like Quarks

Multiple analyses from ATLAS and CMS collaborations, including searches for ttH production, supersymmetric particles and vector-like quarks, observed excesses in the same-sign dilepton channel containing b-jets and missing transverse energy in the LHC Run 1 data. In the context of little Higgs theories with T parity, we explain these excesses using vector-like T-odd quarks decaying into a top quark, a W boson and the lightest T-odd particle (LTP). For heavy vector-like quarks, decay topologies containing the LTP have not been searched for at the LHC. The bounds on the masses of the T-odd quarks can be estimated in a simplified model approach by adapting the search limits for top/bottom squarks in supersymmetry. Assuming a realistic decay branching fraction, a benchmark with a 750 GeV T-odd b-prime quark is proposed. We also comment on the possibility to fit excesses in different analyses in a common framework.


I. INTRODUCTION
It was recently pointed out that both ATLAS and CMS observed in their Run 1 data mild excesses in the same-sign dilepton (SS2l) channel, which contains final states with same-sign dilepton, b-jets, and missing transverse energy (E miss T ) [1]. The excesses were observed in five semi-independent analyses, three from the ATLAS [2][3][4] and two from the CMS [5, 6], using different background subtraction methods. In particular, the ATLAS exotica search [3] and the CMS ttH Higgs search [5] reported about 2σ and 2.5σ significance, respectively, for the excesses. More recently, measurements on the standard model (SM) process ttW from ATLAS [7] and CMS [8] have also reported seeing excesses above the SM expectation in the SS2l channel.
A supersymmetric interpretation of the SS2l excess was put forward in Ref. [1], which proposed using top squarks (stops) or bottom squarks (sbottoms) decaying into two top quarks, two W bosons, and E miss T to explain the excess. However, it is worth recalling that SS2l, b-jets and E miss T are generic signatures of many beyond-the-SM theories and, therefore, not unique to supersymmetry. In fact, it is well-known that many scenarios could "fake" supersymmetry at the LHC, due to the complex environment of pp collision [9]. For example, decays of the stop could be mimicked by that of a fermionic top partner, unless one can measure the spin of the new particles involved in the decay chain, which is a challenging measurement at the LHC [10]. Indeed, should the SS2l excess be confirmed at the LHC Run 2, it would be of paramount importance to determine the specific quantum number of the new particles associated with the excess.
In this study we would like to proceed in an exploratory spirit that is appropriate for this nascent subject, and consider alternative possibilities, other than supersymmetry, for the SS2l excess. The simplest possibility is to invoke a new conserved quantum number at the TeV scale, the T parity, under which the SM particles are neutral and the new particles are charged [11]. Collider phenomenology of T parity is very similar to that of R parity in supersymmetry [12] and KK parity in Universal Extra Dimensions [13], in that all T-odd particles can only be pair-produced and subsequently cascade-decay into SM particles plus the lightest T-odd particle (LTP), which carries away extra E miss T in collider detectors.
This work is organized as follows. In Section II we give a brief overview of T parity and an estimate on the collider bound on the mass of T-odd quarks, followed by Section III which discusses fitting the SS2l excess in the context of SM ttH searches at the LHC Run 1 and beyond. In Section IV we consider whether it is possible to explain using a common benchmark the excesses observed by analyses outside of the SM ttH searches. Then we provide an outlook in Section V.

II. A SIMPLIFIED OVERVIEW OF T PARITY
In this section we provide a "simplified" overview of T parity and consider collider bounds on T-odd quarks.

A. Simplified T Parity
The case for a new symmetry at the TeV scale, under which all SM particles are neutral, is motivated by two considerations [11]: • Precision Electroweak Constraints: If we parameterize flavor-conserving new physics in terms of higher-dimensional operators, precision electroweak constraints indicate that the mass scale suppressing these operators tend to be at around 5 -10 TeV, assuming all dimensionless coefficients to be order unity [14]. This is referred to as the little hierarchy problem, since naturalness principle expects new physics at around 1 TeV scale to stabilize the Higgs mass. If there exists a new symmetry at the TeV scale such that all SM particles are neutral while all new particles are charged under the new symmetry, then new particles would enter into precision electroweak observables only through loop-induced effects, thereby introducing a loop factor of about 1/16π 2 in front of the higher dimensional operators. This allows one to lower the scale of new physics down to below 1 TeV.
• Dark Matter: The existence of dark matter calls for (at least) a new particle that is electrically neutral and stable at the cosmological time scale. A simple possibility is to make the dark matter particle absolutely stable. This can be achieved again by postulating a new symmetry such that all SM particles are neutral. Then the lightest particle charged under the new symmetry will be stable, assuming the symmetry is exact, or almost exact so that the particle is stable cosmologically. If such a particle is also electrically neutral, it can be a dark matter candidate.
There are many realizations of such a new TeV symmetry in explicit models. The most popular models use the simplest symmetry: a new parity (Z 2 ) symmetry. R parity in supersymmetry [12], KK parity in Universal Extra Dimensions [13] and warped extra dimension [15,16], and T parity in little Higgs models [17,18] all fall into this simple category. Nevertheless, larger symmetry groups such as the Z 3 group have also been employed [19].
It is possible to adopt the "simplified model" approach [20], by postulating the existence of a new parity, without committing to a specific model realization. In particular, one could assume there is a "simplified T parity," under which all SM particles are neutral and all new particles are charged. Then at colliders any new particles must be pair-produced, in order to conserve the simplified T parity, and eventually cascade-decay into the LTP. The LTP would serve as the dark matter candidate, assuming it is electrically neutral, and manifest itself as E miss T in colliders. Phenomenology of the simplified T parity at the LHC is very similar to that of supersymmetry with R parity, which is characterized by leptons, jets and can be impersonated in simplified T parity by a T-odd vector-like quark with the same quantum number as the bottom quark: where W ± H is a pair of heavy charged vector boson and A H is a heavy neutral vector boson that is the LTP. For the SS2l excess, the bottom squark decay chain Eq. (1) was discussed to some extent in Ref. [1], which focused on supersymmetric theories. In this work we consider the alternative decay chain in Eq. (2) based on the (simplified) T parity. While both models can produce SS2l events, the kinematic distributions in general are different, and might be used to distinguish different models. Of course, to fully identify the decay chain requires 1 For a related approach to the simplified T parity, see Ref. [21]. measuring the spin of the mother and intermediate particles [22], which would be a top priority should the excess be confirmed at the LHC.
Explicit constructions of the T-odd top and bottom partners in little Higgs models were given in Ref. [18,23,24]. It was shown that in order to cut off the contributions to the µ . In this case, the decay branching fractions of the T-odd b and t are more or less determined by the Goldstone equivalence theorem [25]: The other decay channels, b → bA H and t → tA H , are subdominant. These observations will be taken into account when we consider model-specific T-odd quark decays. 2

B. Collider Bounds on Third Generation T-odd Quarks
Although there are many searches for stops and sbottoms in SUSY at the LHC Run 1, there hasn't been many dedicated searches on the closely related decay chains involving T-odd quarks. 3 In particular, the stop and sbottom searches usually adopt the simplified model approach, by assuming 100% decay branching ratio (BR) into the final states being searched for. Current limits from the LHC Run 1 data on the stop and bottoms are typically around 500 -700 GeV [27].
Since the decays of T-odd quarks often give the same final state signatures as the squarks, in this subsection we will provide a rough estimate on the bounds on the masses of T-odd quarks, based on the experimental searches for third generation squarks. Recasting the search limits on squarks is not a straightforward task, as the signal selection efficiencies depend non-trivially on the kinematics of the decay products [28], which in turn is determined by the mass spectrum of the mother and daughter particles. For example, the strong limits on the sbottom mass in the tripleton channel of the decay chain in Eq. (1) varies significantly when the lightest neutralino [which is assumed to be the lightest supersymmetric particle (LSP)] and chargino masses change [29]. Moreover, the bound disappears completely when the LSP becomes heavier than 170 -240 GeV, depending on the chargino mass. It is clear that a full-fledge study on the experimental constraints of T-odd quarks requires dedicated efforts and is beyond the scope of this paper.
In this work we will settle for a naïve estimate on the experimental bounds on third generation T-odd quarks by translating the strongest limits on the third generation squark masses into upper bounds on the production cross-sections, and then computing the b and t masses that give rise to the same production cross-section as the upper bounds. The outcome based on the simplified model assumption is shown in Fig. 1, where the squark and T-odd quark cross sections are quoted from the LHC SUSY Cross-section Working Group [30,31] and from HATHOR [32,33], respectively. We see that the lower bound on the b and t quarks are about 800 GeV and 825 GeV, respectively. The bounds in Fig. 1 are conservative in the sense that they correspond to the most stringent limits on the m LSP -mb 1 /t 1 plane, which would loosen when the m LSP becomes larger. Beyond the simplified model approach, the exclusion limit degrades quickly as the decay BR decreases from the assumed 100%, because the signal strength is usually proportional to the square of the BR due to the assumption of pair-production of the mother particles.
For example, in the littlest Higgs with T parity model (LHT) [18,25], the decay BR of T-odd b into tW H final state is typically about 55% at around 800 GeV, as shown in Fig. 2.
Consequently, the collider bounds on the T-odd quarks are considerably weaker in a full model than in the simplified model approach. This comparison is shown in Fig. 3, where we compared the bound assuming 100% BR versus 55% for both T-odd b and t . The bounds on the b (t ) mass is only about 680 (700) GeV. Ref. [29], which gives the most stringent limits, while the bounds on the stop are from Ref. [27].
In this plot the decay of squarks and T-odd quarks into the desired final states are assumed to be 100%.

III. SAME-SIGN DILEPTON SIGNALS AT THE LHC
Having estimated the LHC bounds on third generation T-odd quarks, in this Section we consider an interpretation of the SS2l excess from the pair-production of T-odd quarks b .
In particular, the decay chain we concentrate on is  [18,25]: Ref.
[5] we impose further the following kinematic cuts: and SM ttH are then used to calculate the ratio of signal strength µ b /µ ttH . In the end, the total signal strength is which is to be compared with the ATLAS result µ = 2.8 +2.1 −1.9 [4] and the CMS fit µ = 5.3 +2.1 −1.8 [5]. If we consider the "simplified T parity model," namely the 100% branching fraction of b → tW H , a b with a mass of 850 GeV and production cross-section 12.85 fb will generate the similar signal strength.
It is interesting to compare with the SUSY benchmark considered in Ref. [1]. There the spectrum was taken to be: The pair production cross section of 550 GeVt 1 is 45.2 fb at 8 TeV and the branching fraction oft 1 → tW ±χ∓ 1 is close to 100%. The total signal strength is µ = µt 1 + µ ttH = 2.83. This is because the final state particles in our benchmark model are much more energetic due to the heavy mother particle mass as well as the large splitting in the spectrum, and hence are easier to pass the cuts.
At the LHC Run 2 with the same cuts imposed in Eq. (5), the total signal strength for the b benchmark model increases to µ = 3.2. The reason of the increase is similar to the situation in the stop interpretation: the production cross-section for the heavy particles grows at a faster pace than that of SM ttH [1]. The increase is also more significant than that of the SUSY benchmark model [which has µ(13 TeV) = 3.69]. In Fig. 4  GeV.
It will also be interesting to contemplate further kinematic cuts at the LHC Run 2 that could help disentangle the b signal from that coming from the SM ttH or the stop/sbottom SUSY signals. In this regard, we show in Fig. 6 some kinematic distributions of events from the SM ttH, the T-odd b and the stop in SUSY. One sees that b has the hardest spectra among the three benchmarks, which is due to the heaviness of the b in the benchmark, resulting in more energetic decay products. This feature is quite generic, since the fermion has a significantly larger cross section than the scalar at the same mass. Therefore, given a particular signal strength, one can always fit it with a fermion mass that is heavier than that of a scalar. As long as the spectrum is not degenerate, it will result in harder distributions of the decay products.
to enhance the contribution of b to SS2l signal region in the context of the CMS SM ttH analysis. With these further cuts, the total signal strength, in unit of the SM ttH strength, raises dramatically to In other words, the signal would come almost exclusively from T-odd b .

IV. A BROADER PICTURE
So far we have focused on fitting the SS2l excess in the context of the ttH multilepton analyses [4,5], by normalizing to the SM ttH signal strength. Given that several other analyses

A. CMS SUSY Analysis
The SS2l analysis in the CMS SUSY working group in Ref. [6] presented exclusion limits on the sbottom decay chain that are degraded from the expected limits, implying more events were observed than expected. The most significant excess appeared in the signal region SR24, which requires SS2l and N b−jets ≥ 2, N jets ≥ 4, 50 GeV ≤ E miss T ≤ 120 GeV, H T ≥ 400 GeV.
The expected number of events in SR24 is 4.4±1.7 (2.8±1.2) in the low (high) p T region and the observed numbers are 11 (7).
We simulated the contribution to SR24 from both the b and the stop benchmarks. The b (stop) would give rise to 0.4 (3) and 0.3 (2.3) events in the low and high p T region, respectively. It is clear that T-odd b has a lot more difficulty fitting the SR24 excess than the stop benchmark. This is due to the fact that SR24 only selects events with a relatively small E miss T , while the E miss T distribution from b benchmark is much harder than that from the stop, as can be seen from Fig. 6.
It is also interesting to note that the b-tagging selection efficiency is typically 70%. As a result, events contributing to N b−jets = 2 region will also contribute to N b−jet = 1 region.
Therefore  We implemented the above cuts in both the b and the stop benchmarks, which contributed 3.7 and 3.1 events, respectively, to the SR1b signal region. Therefore, the excess observed by the ATLAS SUSY search could be explained by both the b and the stop benchmarks. The selection efficiency for the b benchmark is much better than that from the stop because of the hard cuts on E miss T and H T , which can be seen from the distributions shown in Fig. 6.

C. ATLAS Exotica Analysis
The ATLAS search for heavy vector-like quarks in the SS2l + 2 b-jets region reported an excess in the SR4t3 signal region, which is defined by requiring SS2l and In this signal region 4.4 ± 1.1 ± 1.1 events are expected and 12 are observed. There is also a separate excess in SR4t4, the 3 b-jets category, which sees 6 events while expecting only 1.1 ± 0.9 ± 0.4 events. The p-values for both excesses are about 2σ. However, recall that the ATLAS SUSY search in Ref. [2] sees no excess in SR3b, the three b-jets region.
Implementing the above selection cuts in our simulation, we find that the b benchmark contributes about 2.3 events to SR4t3, while the stop benchmark contributes about 1.8 events. Again b contributes more to the signal region because it has a harder decay spectrum.

V. SUMMARY AND OUTLOOK
In this work we studied the possibility of fitting the SS2l excess observed in the LHC Run 1 data by vector-like quarks that are odd under a new parity at the TeV scale, the T parity.
Phenomenology of T-odd quarks is quite different from those heavy quark decay chains that are being searched for at the LHC so far, since the T-odd quarks decays into SM particles plus the LTP, which carries away extra missing transverse energy at collider detectors. The current bounds on the masses of the third generation T-odd quarks can be estimated by recasting the exclusion limits for stops and sbottoms in supersymmetry, which have similar decays chains to the third generation T-odd quarks,.
We proposed a vector-like quark benchmark containing a 750 GeV b quark, decaying 55% of the time into a SM top quark and a 320 GeV W H . The W H subsequently decays into the SM W boson and the LTP, whose mass is at 66 GeV. Using the ttH searches as a starting point, we normalized the b signal strength to that of the SM ttH expectation and obtained a total signal strength of µ = 2 in the context of ttH searches.
At the LHC Run 2, kinematic cuts are suggested to further enhance the b signal over the SM ttH signal. Moreover, we also studied differences in the kinematic distributions between the b benchmark and the stop benchmark proposed in Ref. [1], and observed that the decay spectra of b are generically more energetic than those from the stop decays. Should the SS2l excess be confirmed at the LHC Run 2, it would be a top priority to determine whether the excess is due to the production of new colored particles, as well as the possible quantum numbers of such new particles.
Looking beyond the excess in the ttH searches, we also considered whether the SS2l excesses in searches for new particles can be explained in a common framework. As it stands there is some difficulty in explaining all the excesses using the b or the stop benchmarks.
However, the number of expected signal events is quite small at the LHC Run 1 and obviously more data from Run 2 is needed to clarify the nature of the excess. In particular, if the excess is due to the pair-production of new colored particles, a small amount of data, of the order of 5 fb −1 , at the 13 TeV LHC is sufficient to reach the same sensitivity as the entire Run 1 dataset.
Last but not least, even if the SS2l excess disappears as more data is collected, it is clear that searches for vector-like quarks at the LHC need to be extended to the well-motivated scenario of T-odd quarks, whose decay phenomenology requires dedicated analyses that are not covered by current searches for vector-like quarks.