Constraining Higgsino Kink Tracks from Existing LHC Searches

Considering the supersymmetric model with the long-lived charged Higgsino, we discuss how Higgsino kink signals can show up in the latest LHC Disappearing Track (DT) and stable chargino searches. We derive constraints on the Higgsino kink signal, and characterize it in comparison to the Wino DT and the slepton kink track. We also discuss how to infer Higgsino model kinematics.


I. INTRODUCTION
The "kink" signal is often predicted in models beyond the Standard Model (SM). In the gauge-mediated supersymmetric models, for example, the charged slepton next-to-lightest superparticle (NLSP) becomes long-lived due to very weak interactions of the gravitino lightest superparticle (LSP), and it can subsequently decay to a charged lepton within a detector. Since the momentum directions of these charged particles are different in general, tracks will suddenly be deflected (or kinked) at the decay position, and their tracks form a kink track [1][2][3]. 1 The kink has been searched for at the LEP [4,5], but no direct constraints on it have been reported from LHC experiments. A recent study by theorists [6] claims that the slepton kink can be constrained by Disappearing Track (DT) searches at the LHC [7,8]. The DT search is aimed to look for a different signal -the charged Wino decaying inside the LHC tracker to the almost degenerate neutral Wino, producing only an invisibly soft pion [10][11][12][13] and leaving no energetic tracks afterwards. But the kink track resembles the DT in the sense that the original track does not extrapolate to the end of detectors.
In this paper, we consider another supersymmetric model where the charged Higgsino can be the kink track. The charged Higgsino is a long-lived NLSP, eventually decaying to a weakly interacting gravitino LSP and a charged lepton, producing a kink signal. We consider this model because this model realized in the MSSM contains light Higgsinos that might be related to the weak-scale hierarchy problem. This model has been studied in only a small number of theoretical papers [14], and our study derives the latest new constraints on this model. The Higgsino kink also has several features that can be compared with often discussed signals of the Wino DT and the slepton kink. Although the Wino DT produces only invisibly soft pions, the Higgsino kink can accompany energetic leptons and/or quarks (from intermediate W boson decays) moving in different directions from original Higgsino momentum. Also, although the slepton kink decays to only a single lepton (for simplicity, we do not consider the tau decay), the Higgsino kink can decay to multiple quarks as well as to neutrinos. In addition to these differences in decay kinematics, the kinematics of the fermionic Higgsino pair production is different from that of the scalar slepton pair, as will be discussed. We will study how such differences can become apparent and can be utilized in DT searches to search for Higgsino kinks and to characterize them.
The rest of this paper is organized as follow. In Sec. II, we introduce the model of the long-lived charged Higgsino NLSP. In Sec. III, we present our results for how Higgsino kink would show up in the LHC DT observables and discuss how they can characterize similarities and differences among the aforementioned kink and DT models. Then we derive constraints on the Higgsino kink in the later part of Sec. III, and we conclude in Sec. IV.

II. MODEL -CHARGED HIGGSINO NLSP AND GRAVITINO LSP
We consider the charged Higgsino as NLSP in the minimal supersymmetric standard model (MSSM) extended with the weakly interacting gravitino LSP. The chargino NLSP is long-lived and eventually decays to the gravitino LSP via H ± → W ± G, leaving kink tracks [14].
The tree-level mass splitting between the lightest neutral Higgsino and the charged Higgsino is given by [14] (for M 2 > 0) where M 1,2 are gaugino mass parameters and µ is the Higgsino mass parameter. t W ≡ tan θ W is the weak mixing parameter, and tan β ≡ v u /v d is the ratio of the vacuum expectation values of the two Higgs doublets. The ∆m tree is negative when M 2 > |M 1 |( µ) > 0 and M 1 < 0 with small tan β, so that the charged Higgsino can be the lightest Higgsino components; see Refs. [14,15] for more discussions.
The charged Higgsino decay width is [16] Γ( where the SUSY breaking scale F is related to the gravitino mass m 3/2 by m 3/2 =  [17]. We take the Higgsino mass and the lifetime cτ as free parameters. We assume that slightly heavier neutral Higgsino components promptly decay to the charged Higgsino (the decay to the gravitino LSP is slower); this is realized for |M 1 |, |M 2 | 2 TeV and small tan β. Thus, the effective production rate of the charged Higgsino pair is a sum of all possible pair production of neutral and charged Higgsinos (and each event eventually contains two long-lived charginos). We have used cross-sections from the LHC Working Group results [18].

III. ANALYSIS AND RESULTS
The (Higgsino) kink track can be characterized by two types of features. First, a Higgsino kink track is usually isolated from other leptons or quarks. Second, the kink track from a new physics beyond the SM is likely a heavy particle while SM backgrounds will be from light particles. Meanwhile, the Wino DT is also characterized by similar features: no energetic particles along the DT direction and the heaviness of a DT. Thus, we consider the following observables, similarly to the latest LHC DT search observables: (1) E ∆R<0.5 , total energy of all charged leptons and quarks within ∆R < 0.  2. ∆Rmin, the minimum ∆R between the 300 GeV charged Higgsino and charged leptons and quarks. The dilepton channel (blue) where both W bosons decay to leptons and the dijet channel (red) where both W decay to quarks are compared. Higgsino lifetime cτ =10 cm (solid), 100 cm (dashed) and the massless LSP. The fraction of events decaying within the appropriate part of detector is shown. CMS pre-selection cuts are applied. For comparison, slepton kinks (gray) are also shown.
the charged Higgsino, and (2) p T (Higgsino). In addition, we look at (3) ∆R min , the minimum ∆R between the charged Higgsino and any charged leptons or quarks, which can characterize charged activities along the original Higgsino direction.
We restrict to the events where the charged Higgsino decays within the relevant part of detectors: 30cm ≤ L 90cm [7,8], where the L is the decay length in the transverse plane. For this, each event generated by MadGraph5 v2.3.3 MSSM model [9] is reweighted by the probability for at least one Higgsino to decay within this part of detector. For the given transverse momentum γβ and the lifetime cτ of the charged Higgsino, the probability for the Higgsino to decay within the part is The probability for at least one Higgsino to decay properly can be obtained from this. For simplicity, we ignore longitudinal decay lengths which would not significantly change our results. We show the resulting observables for the 300 GeV Higgsino and the massless LSP in Fig. 2, 3, and 4. These are obtained after applying CMS pre-selection cuts (triggering the pair production of charginos plus initial-state radiation) [8]: at least one well-separated prompt jet and missing transverse energy (MET) with p T > 110 GeV, MET > 110 GeV, ∆φ(jet, MET)> 0.5, and no isolated prompt leptons. Although we base on the CMS analysis for our convenience, similar discussions can be made with the ATLAS analysis [7]. Slepton kink events are also shown for comparison, where the long-lived slepton decays to the gravitino LSP and a charged lepton.
"∆R min " shown in Fig. 2   along the original direction of the Higgsino kink. It is usually larger in the dilepton channel (where both W bosons decay to leptons) than in the dijet channel (where both W decay to quarks) because a smaller number of charged leptons and quarks are produced. Charged activities near the Wino DT is also small. Thus, the latest CMS DT search requires charged activities within ∆R < 0.03 of the DT to be very small: i∈∆R<0.03 p i T /p T (DT) < 0.05 [8] (the scalar sum of all track p T compared to the p T (DT)). To approximately mimic this cut, we require ∆R min > 0.03 in our final analysis; we can see that more than 95% of Higgsino kink events can satisfy this cut.
∆R min of the slepton kink may be expected to be similar to that of the dilepton channel since the same number of final-state charged leptons are produced. However, the slepton kink has somewhat smaller ∆R min as shown in Fig. 2 mainly due to the threshold suppression of the slepton scalar pair production as will be discussed.
Higgsino kinks have high p T ∼ m(Higgsino) ∼ O(100) GeV as shown in Fig. 3. The p T (Higgsino) distribution for the appropriate part of detector does not depend on the decay modes of the W boson, but it de-  [7,8] and stable chargino searches [19] with 20/fb. The region left to each curve is constrained. The uncertainty band is obtained by varying selection efficiencies by 50%.
pends sensitively on the chargino proper lifetime cτ as its decay length is (probabilistically) determined by cτ and its boost. The Higgsino kink and the Wino DT have the same p T distribution because both are Drell-Yan pair produced electroweak fermions; whereas, the slepton kink has harder p T distribution because a scalar pair Drell-Yan production suffers a threshold suppression due to angular momentum conservation. Thus, the p T distribution can be used to distinguish slepton kinks from others. The threshold suppression also explains why slepton kinks tend to have smaller ∆R min in Fig. 2; harder sleptons with higher boosts lead to more collimated decay products. Following the latest CMS DT search, we apply p T (Higgsino)> 50 GeV in our final analysis; again, the majority of the Higgsino kink events satisfy this. "E ∆R<0.5 " shown in Fig. 4 is the total energy of all charged leptons and quarks within ∆R = 0.5 of the Higgsino. Since the decay products of the chargino are energetic (p T ∼ O(100) GeV), either very small (if no charged decay products are nearby) or large (if at least one is nearby) value is predicted. The dilepton channel, having a smaller number of charged leptons, tends to have more events at small E ∆R<0.5 than the dijet channel. On the other hand, the Wino DT always has small value of E ∆R<0.5 . Following the latest CMS DT search, we require E ∆R<0.5 < 10 GeV in our final analysis; still more than a half of Higgsino kink events would satisfy this.
From these considerations, we expect that a large fraction of Higgsino kink events would be captured by LHC Wino DT searches; in other words, our kink track that suddenly changes its direction (due to charged Higgsino's decays) can mimic DTs that actually disappear suddenly. We derive constraints on the Higgsino kink by applying the final dicsovery cuts mentioned above. We first obtain the ratio of Higgsino kink efficiency to that of the Wino DT. 2 For the given Higgsino mass, we then find the cτ of the Higgsino that can compensate the efficiency ratio and different effective chargino production rates. It is the constraint shown in Fig. 5. 3 Constraints are strongest near cτ 50 cm, excluding up to the 360 GeV Higgsino. There is no constraint for cτ 3 cm as the Higgsino has to decay after 30 cm; while a longer lifetime cτ 1 m is already constrained by stable chargino searches [19]. The blue band is obtained by varying our efficiency by 50%. We have assumed the massless LSP in this figure. But current CMS (and AT-LAS too) selection efficiencies do not depend strongly on the mass difference between the NLSP and LSP as long as the W boson can be on-shell; the massless LSP case has only 10% larger efficiency than that of the compressed case. Thus, similar results will be obtained for other LSP masses.
In comparison with the Wino DT, the efficiency of the Higgsino kink is usually about the half of the Wino DT's with the same mass; somewhat smaller (larger) for the lighter (heavier) Higgsino. E ∆R<0.5 in Fig. 4 is a main difference for this. The total effective charged Higgsino pair production rate (sum of all possible neutral and charged Higgsino pair productions in our model; see Sec. II) is about 20% smaller than the effective Wino DT production rate. In all, the maximum reach is reduced from about 520 GeV for the Wino DT [8] to 360 GeV for the Higssino kink. But the strongest constraints are similarly obtained for cτ ∼ O(10) cm due to similar p T (Higgsino/Wino) distributions.
The NLSP-LSP mass difference of the Higgsino model can be apparent in the p T of the charged lepton or the quark nearest to the Higgsino. This reasonably measures the energy of chargino decay products, and it directly probes the mass difference. It is because the smaller the mass difference, the softer the visible decay products. This is clearly shown in Fig. 6. This result is obtained after all CMS discovery cuts, so presumably backgrounds are already small there. Thus, this distribution can be useful in inferring mass parameters or in distinguishing Higgsino kink versus Wino DT that does not accompany energetic charged particles nearby.

IV. CONCLUSION
We have constrained the supersymmeric model of the long-lived charged Higgsino NLSP and the weakly interacting gravitino LSP from LHC searches. For cτ ∼ 50 cm, up to the 360 GeV Higgsino can be excluded by the latest 8 TeV LHC DT searches. But we encourage a more dedicated study by including various realistic effects that our parton-level results do not take into account. We have also discussed similarities and differences among the Higgsino kink versus the much-discussed Wino DT and the slepton kink. The differences are most apparent in observables measuring the activity nearby DT/kink (for Higgsino kink versus Wino DT) and in the p T of DT/kink (for Higgsino kink versus slepton kink). Moreover, we have discussed that the p T of the charged lepton or the quark nearest to the Higgsino can be used to infer the NLSP-LSP mass difference. Similar discussions can also be applied to various other models including the weakly interacting axino LSP [20][21][22] and R-parity violating slow decays (through LLE c and LQD c operators) of the charged LSP.
Note added: After the submission of the first draft to arXiv, more experimental information became available, and we have revised our study significantly. Some conclusions and focus are changed accordingly.
Acknowledgements. We thank Brian Batell, Suyong Choi, Andy Haas, Teruki Kamon, and Jianming Qian for discussions. SJ thanks the Aspen Center for Physics where part of the work was done. HL appreciates the hospitality during his visit to KIAS. The work of SJ was supported in part by National Research Foundation of Korea (Grant 2013R1A1A2058449) and US Department of Energy under contract DE-AC02-76SF00515. The work of HL was supported in part by CERN-Korea fellowship and IBS (Project Code IBS-R018-D1). 2 The number of background events reported in Ref. [8] is less than 10 for m H + 200 GeV, so our approximation in this way is somewhat reliable. 3 We encourage a more dedicated study improving our parton-level results by including realistic detector effects, parton showering and reconstruction algorithms; our treatment of particle momentum is ideal and accurate efficiencies of our model can only be obtained by dedicated simulation studies.