Discovery reach for wino and higgsino dark matter with a disappearing track signature at a 100 TeV $pp$ collider

Within the theory of supersymmetry, the lightest neutralino is a dark matter candidate and is often assumed to be the lightest supersymmetric particle (LSP) as well. If the neutral wino or higgsino is dark matter, the upper limit of the LSP mass is determined by the observed relic density of dark matter. If the LSP is a pure neutral state of the wino or higgsino, the lightest chargino state is expected to have a significant lifetime due to a tiny mass difference between the LSP and the chargino. This article presents discovery potential of the 100 TeV future circular hadron collider (FCC) for the wino and higgsino dark matter using a disappearing-track signature. The search strategy to extend the discovery reach to the thermal limits of wino/higgsino dark matter is discussed with detailed studies on the background rate and the reference design of the FCC-hadron detector under possible running scenarios of the FCC-hadron machine. A proposal of modifying the detector layout and several ideas to improve the sensitivity further are also discussed.


Introduction
Astrophysical observations of galaxies and the largescale structure in the universe strongly indicate that dark matter predominate the matter contents in the universe. The nature of dark matter is, however, still unknown and the Standard Model (SM) of particle physics has no counterparticles acting as dark matter. Within the theory of supersymmetry (SUSY), a dark matter a e-mail: ryu.sawada@cern.ch b e-mail: koji.terashi@cern.ch candidate is a neutralino, which is a mixed state of neutral supersymmetric partners of the SM U(1)×SU (2) gauge bosons and the two SU(2) Higgs doublets. The lightest neutralino (χ 0 1 ) is often assumed to be the lightest supersymmetric particle (LSP). If we assume that dark matter is the neutral wino (higgsino), which is a supersymmetric partner of the SM SU(2) gauge boson (Higgs boson) and that it was produced in thermal processes, the upper limit of the LSP mass is determined by the observed relic density of dark matter to be 3 (1) TeV for the wino (higgsino) LSP scenario [1,2].
When the pure wino or higgsino state is the LSP, the mass of theχ 0 1 is highly degenerate with the lightest chargino (χ ± 1 ) because the wino/higgsino belongs to the same SU(2) multiplets. The mass difference between theχ ± 1 and theχ 0 1 is typically 160 (355) MeV for the wino (higgsino) LSP scenario, as evaluated from the calculation of radiative correction in [3,4]. Due to the small mass difference, theχ ± 1 has a long lifetime of approximately 0.2 (0.023) ns for the 3 TeV wino (1 TeV higgsino) at the thermal limit [5]. With this mass difference, theχ ± 1 decays into aχ 0 1 and a charged pion predominantly. Theχ 0 1 passes through the detector without any electromagnetic or strong interactions while the pion is not reconstructed as a track in the innertracking detector due to very low transverse momentum (p T ). Therefore, such long-lived charginos are observed in the detector as a short 'disappearing' track, which has hits in several silicon-detector layers near the interaction point and no associated hits after a certain radius, with a typical length of O(1-10) cm in collider experiments [6][7][8].
In this paper, a search strategy for the long-lived charginos with a disappearing-track signature at the 100 TeV FCC-hadron machine (FCC-hh) is discussed with the reference FCC-hadron detector. The expected arXiv:1901.02987v1 [hep-ph] 10 Jan 2019 sensitivity is presented under the scenario of 30 ab −1 of pp data collected under possible running conditions of the FCC-hh. In earlier studies on the prospects of disappearing-track searches in future pp colliders [9][10][11][12][13], effects of multiple pp collisions occurring simultaneously with a signal event (pileup) were assumed not to alter the results significantly by improving the analysis method or the background rate was scaled with a certain factor without detailed studies. In this paper, the background rates were estimated based on the detector layout and simulation. In addition, several ideas that could potentially improve the sensitivity further are discussed.

Simulation samples
Simulated samples of Monte Carlo (MC) events are used to obtain the kinematic distributions of the signal and background processes. The wino mass spectrum is calculated using Softsusy 3.7.3 [14] assuming the minimum Anomaly Mediated SUSY Breaking (AMSB) model with tanβ = 5, positive sign of the higgsino-mass term and the universal scalar mass (m 0 ) of 20 TeV. This setting naturally accommodates the observed Higgs boson mass of 125 GeV without maximal mixing. The higgsino mass spectrum is calculated by Susyhit [15] assuming the general minimal supersymmetric extension of the SM (MSSM) with tanβ = 5. Other SUSY particle masses are set to 100 TeV to prevent the higgsino from being affected from loop effects. Only the pure wino and higgsino LSP scenarios are considered here. The production cross sections for the winos and higgsinos in 100 TeV pp collisions are calculated using Prospino2 [16]. The kinematic distributions are simulated using MadGraph5 aMC@NLO 2.3.3 [17] assuming the wino spectra. The same kinematic distributions are assumed for the wino and higgsino if they have the same masses. In this study, pair-production processes of charginos and neutralinos are considered: are degenerate in the higgsino case, the difference in the cross section or decay processes betweenχ ± 1χ 0 1 and χ ± 1χ 0 2 is negligible. Therefore, only theχ ± 1χ 0 1 process is generated and the cross section is doubled to account forχ ± 1χ 0 2 . Minimum-bias collisions are produced using Pythia 8.230 [18] and are overlaid to evaluate contribution from non-genuine tracks caused by the random combination of hits in the inner-tracking detector, which are referred to as fake tracks hereafter. The details of the minimum-bias event overlay are described in Sect. 5.
SM background processes such as W/Z-boson production process in association with jets, top-quark pair and single top-quark production processes are simulated using MadGraph5 aMC@NLO and the generated samples are used to evaluate the selection efficiency. The kinematic distributions of jets, electrons, muons and the magnitude of missing transverse momentum (E miss T ) are processed using Delphes 3.4.2 [19] with the reference FCC-hh run card to account for the detector reconstruction efficiency and energy/momentum measurement. Jets are reconstructed from particle-flow candidates using the anti-k t algorithm [20] with a distance parameter of 0.4.

Inner Tracker Layout
The reconstruction of short tracks relies on the performance of the inner-most tracking system. This paper considers the reference design of the FCC-hadron detector [21] as a baseline. For the inner tracker, the most relevant part for the disappearing-track search is the inner-pixel layers, in particular those in the barrel region where the signal acceptance is high. In the present study three different layouts are considered for the barrel region to evaluate impact of the inner-tracker configuration.
The default layout (labelled #1) is the one in the reference design, containing four pixel layers within R = 150 mm in radial distance from the beam line. In the second layout (#2) a pixel layer is added at R = 200 mm to the default layout and all the other layers are unchanged. The third layout (#3) considers an additional pixel layer at R = 150 mm and the inner four pixel layers are moved closer to the beam line. The R-z views of all the three layouts are shown in Fig. 1. 1 The outer macro-pixel layers located at R ≥ 270 mm in the barrel region or the forward tracking system are unchanged from the reference design for all the three layouts. The first five inner layers are most relevant for the shorttrack reconstruction as described below. Their radial positions considered in this study are summarized in Table 1.

Event Selection
The chargino/neutralino-pair production is characterized by the presence of significant E miss T associated with 1 We use a right-handed coordinate system, with the x-, y-, and z-axis are defined by the horizontal, the vertical and the beam directions respectively. Its origin is at the nominal pp interaction point. The pseudorapidity η is defined in terms of the polar angle θ by η = − ln tan(θ/2).  Fig. 1: Three barrel inner-tracker layouts considered in this study, the default layout #1 (top), the alternative layout #2 (middle) and the alternative layout #3 (bottom). The difference between different layouts is restricted within R ≤ 200 mm and |z| ≤ 82 cm, and the detector configuration is identical outside the region. The contour drawn behind the layouts shows the number of chargino decays in 3 TeV wino signal events with 30 ab −1 at a given position. The analysis considers the region |η| < 1, denoted by the dotted lines. theχ 0 1 and one or two disappearing tracks from theχ ± 1 . The SM background is suppressed by requiring an existence of a high-p T disappearing track. The signal event will presumably contain high-p T jet(s) as well, originating from initial state radiation, and they provide the means of triggering on the event in conjunction with E miss T . The analysis therefore selects events which contain at least one high-p T jets, large E miss T and no lepton (electron or muon) with p T > 10 GeV. Figure 2 shows the leading jet p T and E miss T distributions of signal and background events selected with the lepton veto and no requirement on the jet p T or E miss T . The thresholds on the jet p T and E miss T are determined, for each configuration of the inner-tracker layouts, by maximizing the sensitivity based on the signal acceptance and background rate (summarized later in Tables 3 and 4). Although the background yield can be reduced significantly by applying high thresholds to the leading jet p T and E miss T , there will be no discovery power for the wino or higgsino at the thermal limit if the disappearingtrack information is not used in the search.
The chargino-track reconstruction requires the particle to traverse a certain minimum number of innertracking layers before decay. The ATLAS search with 36 fb −1 of Run 2 data at √ s = 13 TeV [22] used a short-track reconstruction with only hits in the four pixel layers. As the chargino lifetime becomes short, the acceptance for the chargino-track reconstruction increases with decreasing radial position of the outermost layer necessary for the track reconstruction. On the other hand, the use of short tracks with less number of layer hits will usually increase fake tracks (discussed in Sect. 5). In order to assess the impact on the sensitivity from different requirements on the number of hits used in the track reconstruction (N hit layer ), the track reconstruction is performed with the requirements of N hit layer = 4 and N hit layer = 5 separately. When the N hit layer = 4 (5), the four (five) inner-most layers with hits are used to reconstruct tracks and the fifth (sixth) layer or above is not considered in the track reconstruction. The contribution from fake tracks increases significantly at high |η| due to higher hit oc- cupancy. Figure 3 shows normalized η distributions of fake tracks (details in Sect. 5.2) and chargino tracks in 3 TeV wino signal events. Given the dramatic increase of fake background with increasing |η|, the analysis requires a candidate track to be within |η| < 1 (further discussion is in Sect. 7).
The signal acceptance for the track requirement is estimated based on the tracker geometry and the chargino lifetime by assuming that the chargino tracks can be reconstructed at 100% of time if the charginos traverse at least four or five inner-most layers before decay. Table 2 summarizes the signal acceptances for the wino and higgsino under the three scenarios of the inner-tracker layouts.

Background estimation
Backgrounds for the disappearing-track signature are categorized into two components [22]. The first one is physical background, which arises from charged particles (mainly electrons or charged pions) scattered by the material in the inner tracker. The second one is unphysical background, which arises from fake tracks. The estimate of both background components is discussed separately below.

Physical background
The physical background due to particle scattering is estimated by assuming that it is solely determined by the amount of detector material passed through by particles. Under this assumption the physical background rate at the FCC-hh is obtained by scaling the background rate measured in ATLAS [22] by the ratio of the material budgets in the FCC-hh and ATLAS inner trackers [23]. A support structure of the FCC-hh tracker has not been fully defined yet, so the expected impact from particle scattering with support material is accounted for by assuming the same relative amounts of tracker layers and support material measured in AT-LAS. The physical background is dominated by W (→ ν)+jets processes, which have large E miss T and a highp T track associated with the lepton from the W boson decay. Therefore, the final estimate of the physical background is obtained by correcting the scaled background rate above for the cross section and kinematic selection efficiency for W +jets processes at √ s = 100 TeV.  Table 2 Signal acceptance for the 3 TeV wino and 1 TeV higgsino with |η| < 1 for the three innertracker layouts. The acceptances are provided separately for the requirements of N hit layer = 4 and 5. The alternative layout #3 has significantly higher acceptance than the others because the relevant layers are located closer to the beam line, particularly for the higgsino case with N hit layer = 5.

Background with fake tracks
The fake-track background is evaluated by counting the number of reconstructed tracks with good quality using simulated minimum-bias events. The samples of minimum-bias events are generated with the average number of pp interactions per bunch crossing ( µ ) of 200 or 500, and simulated for the FCC-hadron detector using Geant4 [24]. The two µ values are assumed to be representative of FCC-pileup conditions. Pileup collisions are produced with two different models of soft QCD processes: the first model considers only non-diffractive processes while the second one the mixture of diffractive and non-diffractive processes with the relative fraction determined according to the Pythia cross sections 2 . The former model predicts more fake tracks than the latter due to the absence of diffractive events, and is used as the nominal pileup model to be more conservative. The track reconstruction is performed as follows. First, the track-seed finding is done by finding a group of inner-tracker hits within a narrow region with a size of about 10 −3 × 10 −3 in η and φ. Second, a track is reconstructed by fitting a line with a given curvature to the collection of the hits. The fit is iterated by removing hits which are far from the line. The following quality criteria are required to select chargino-track candidates: p T > 100 GeV, |η| < 1, χ 2 -probability of the track fit > 0.1, at least one hit in each layer and the association to the hard-scattering vertex by applying the nominal impact-parameter requirement (transverse and longitudinal impact parameters of |d 0 | < 0.05 mm and |z 0 | < 0.5 mm). The final number of fake tracks is obtained by counting the number of tracks with loose impact-2 The non-diffractive processes are set using the SoftQCD:nonDiffractive option while the diffractive ones using SoftQCD:singleDiffractive and SoftQCD:doubleDiffractive options in Pythia.
parameter requirements (|d 0 | < 1 mm, |z 0 | < 250 mm) and scaling that to the nominal impact-parameter requirement assuming a uniform distribution of d 0 and z 0 . Figure 4 shows the estimated probability of finding a fake track per event as a function of the µ . The faketrack finding probability increases significantly with increasing µ . However, the probability can be reduced by 2-3 orders of magnitude by changing the N hit layer requirement from 4 to 5. The number of fake-track backgrounds is estimated by multiplying the number of SM background events satisfying the kinematic selection criteria (in Sect. 4) by these probabilities.
The time resolution of pixel detectors can be better than 50 ps by utilizing, for example, low-gain avalanche detectors [25,26]. Therefore, the time information could be used in the track fit as an additional track parameter to potentially improve the quality of the track. The reduction of fake-track background by requiring a good time quality (i.e, small χ 2 in the estimation of track time obtained from the associated pixel hits) is investigated using the simulated minimum-bias samples with µ = 500. In the simulation, the x, y, z positions and the time of pp collisions are distributed randomly according to Gaussian probability-density functions with the standard deviations of 0.5 mm, 0.5 mm, 50 mm and 160 ps, respectively (without assuming any correlations among them). The χ 2 is computed assuming the constant time resolution of 50 ps for a single layer hit. The fake-track background is found to be reduced by 96% by requiring the χ 2 -probability to be larger than 0.05 for tracks reconstructed with N hit layer = 4 when pp collisions occur in bunch crossings separated by 25 ns. The signal efficiency for this selection is approximately 95%, evaluated using a sample of single muon events.

Result
The discovery sensitivity with the requirement of N hit layer = 5 is higher than that with N hit layer = 4 due to much lower 5-layer-hit fake-track background, as indicated in Fig. 4. The fake-track background rate appears to be rather insensitive to the inner-tracker layout for a fixed N hit layer requirement. Therefore, the alternative layout #2 is expected to have sensitivity in-between the default and alternative layout #3 (hence not shown in the rest of the figures) because the signal acceptance is a factor 2-4 smaller for the #2 than #3 (see Table 2). Figure 5 shows the expected discovery significance for the wino (higgsino) LSP model with the proper lifetime of 0.2 (0.023) ns with 30 ab −1 with the requirement of N hit layer = 5, assigning 30% systematic uncertainty in the background yields. The background uncertainty is based on the estimate in the ATLAS Run 2 analysis [22]. The time information discussed in Sect. 5.2 is not used in Fig. 5. The discovery significance is obtained using an approximate formula to calculate the significance of a signal+background hypothesis against a background-only hypothesis, taking into account the background uncertainty [27,28]. The sensitivity is degraded significantly when the µ increases from 200 to 500, even if five layer hits are used in the track reconstruction. The alternative layout (#3) clearly improves the sensitivity, reaching well above the 5σ discovery for the 3 TeV wino in both pileup scenarios. For the 1 TeV higgsino the 5σ discovery is also feasible in both pileup scenarios but with less margin for the high pileup case. The discovery significance is re-evaluated using the time information and the result is shown in Fig. 6. The sensitivity is restored at the high pileup scenario and the degradation with increasing µ is much more relaxed. The 5σ discovery is also reached with sufficient margin for the 1 TeV higgsino with 30 ab −1 or even less if the time information is used. The discovery sensitivities for the three tracker layouts at the pileup scenarios of µ = 200 and 500 (without time information) are summarized in Table 3 and 4, respectively. The reach of the FCC-hadron machine is shown in Fig. 7 as well as the HL-LHC reach and the current observed limits by the ATLAS experiment. According to the present study the FCC-hadron machine has potential for answering conclusively yes or no to the thermal production of pure wino or higgsino dark matter.

Conclusions and discussions
The discovery potential for the 3 TeV wino and 1 TeV higgsino scenarios has been evaluated using a disappearingtrack signature in pp collisions at √ s = 100 TeV with the FCC-hadron detector. The FCC-hh sensitivity is very promising to both the wino and higgsino scenarios, reaching up to the thermal limits, with inner-tracker layouts optimized for a short-track reconstruction, even in the high pileup environment of µ = 500. If the time resolution of ∼50 ps is achievable at pixel-hit level, wino or higgsino dark matter will be confirmed or refuted up to the thermal limits with an integrated luminosity of 30 ab −1 or smaller.
Finally, we discuss several ideas that could potentially improve the sensitivity further. The |η| range could be expanded from |η| < 1 (used in the present analysis) to e.g, |η| < 4, to improve the overall signal acceptance, then split into smaller |η| regions to optimize the sensitivity. As shown in Fig. 3, the contribution from fake tracks, that grows very rapidly with increasing |η|, will have to be evaluated for the determination of the |η| regions. By tilting the sensor modules and making overlaps between them, the charginos produced at the interaction region may pass through multiple modules in a single layer. This will potentially help improving the acceptance for a short-lifetime signal without decreasing the required number of hits in the track reconstruction, therefore without increasing the fake-track background   Table 4 Signal and background yields as well as the discovery significance for the 3 TeV wino and 1 TeV higgsino with 30 ab −1 with the requirements of N hit layer = 5 under the pileup condition of µ = 500. The fake-track background rate is assumed to be same for both alternative layouts. The kinematic selection requirements on the leading jet p T and E miss T are also given.
level. It would be also beneficial if the barrel pixel layers can be moved even closer to the beam line than the alternative layout #3. This has a large impact on the higgsino signal acceptance. Another interesting avenue to pursue is to use the speed of particles relative to the speed of light (β), which can be measured using the time information of pixel hits, to reject background and characterize the observed disappearing-track signature. With the single muon sample, the β resolution is found to be about 14% (with the pixel-hit time resolution of 50 ps) for the tracks with N hit layer = 5 under the alternative tracker layouts. The estimated β resolution is roughly proportional to the pixel-hit time resolution. The fake-track background will be reduced by requiring the measured β to be compatible with a particle (produced from the primary pp collision vertex) with a certain range of momentum and mass. Because new heavy particles move more slowly in the detector than SM particles, the β could be used as an additional discriminant to separate the signal from scattered SM particles. The grey (red) bands show the significance using the default (alternative) layout #1 (#3). The difference between the solid and hatched bands corresponds to the different pileup conditions of µ = 200 and 500. The band width corresponds to the significance variation due to the two models assumed for soft QCD processes.  Fig. 6: Expected discovery significance for the wino (top) and higgsino (bottom) models with 30 ab −1 with the requirements of N hit layer = 5 and a good time-fit quality. The background reduction rate with the time information is assumed to be the same for both pileup conditions. The grey (red) bands show the significance using the default (alternative) layout #1 (#3). The difference between the solid and hatched bands corresponds to the different pileup conditions of µ = 200 and 500. The band width corresponds to the significance variation due to the two models assumed for soft QCD processes.  Fig. 7: Summary of reach for charginos using a disappearing-track signature. The FCC-hh 5σdiscovery reach assuming the alternative layout, the average number of pp interactions per bunch crossing of 500, only non-diffractive soft-QCD processes and a use of the time information of the pixel detector is shown in blue. The reach for the wino scenario is obtained by extrapolating the sensitivity curve in Fig. 6. The 5σ-discovery reach and the 95% CL exclusion sensitivity are taken from [29]. The observed exclusion in the ATLAS search with 36 fb −1 of Run 2 data [22,30] are shown as well.