Indirect Studies of Electroweakly Interacting Particles at 100 TeV Hadron Colliders

There are many extensions of the standard model that predict the existence of electroweakly interacting massive particles (EWIMPs), in particular in the context of the dark matter. In this paper, we provide a way for indirectly studying EWIMPs through the precise study of the pair production processes of charged leptons or that of a charged lepton and a neutrino at future 100 TeV collider experiments. It is revealed that this search method is suitable in particular for Higgsino, providing us the $5\sigma$ discovery reach of Higgsino in supersymmetric model with mass up to 850 GeV. We also discuss how accurately one can extract the mass, gauge charge, and spin of EWIMPs in our method.


Introduction
ElectroWeakly Interacting Massive Particles (EWIMPs) are theoretically well-motivated particles that appear in many models beyond the standard model (SM). They are widely discussed in the context of the dark matter (DM), with identifying an electrically neutral (or milli-charged) component as the DM. An attractive feature of this scenario is that the vanilla thermal freeze-out scenario predicts the correct amount of the relic abundance for the EWIMP mass range of O(1 -10) TeV, and this mass range is within the scope of current and future experiments. Well-known examples of the EWIMPs are Higgsino and Wino that arise within the supersymmetric extension of the SM. Assuming that Higgsino (Wino) is the lightest supersymmetric particle and its stability is assured by the R-parity, its thermal relic abundance becomes consistent with the DM abundance if the mass is 1.1 TeV [1,2] (2.9 TeV [2][3][4][5]). Another example is the minimal dark matter [1,6,7], where a particle with a large SU (2) L charge is identified as the DM. The stability is automatically assured since operators that cause its decay are suppressed by the cut-off scale of the theory thanks to the large SU (2) L charge, provided that one chooses a correct combination of the charge and spin. A 5-plet Majorana fermion with a mass of O(10) TeV is the most popular in this context, but there are also other possibilities, including both scalar and fermionic particles.
EWIMPs are extensively searched for by many experiments, including DM direct, indirect detections and collider searches (in particular, the mono-X search and the disappearing charged track search). While EWIMPs with relatively large SU (2) L charges such as Wino and the 5-plet fermion are promising for these searches, Higgsino is typically more challenging to probe [8]. Given this situation, another search strategy is proposed [9][10][11][12][13][14][15][16] that probes EWIMPs via the electroweak precision measurement at colliders. It utilizes a pair production of charged leptons or that of a charged lepton and a neutrino, where EWIMPs affect the pair production processes through the vacuum polarizations of the electroweak gauge bosons. It is an indirect search method in the sense that it does not produce on-shell EWIMPs as final states. The current status and future prospects have been analyzed for LHC, ILC, CLIC, and 100 TeV colliders [17][18][19], indicating that it provides a promising way to probe Higgsino as well as the other EWIMPs. A virtue of this method is that it is robust against the change of the lifetime and the decay modes of EWIMPs and whether an EWIMP constitutes a sizable portion of the DM or not. Another important point is that, due to EWIMPs, the invariant mass distributions of the final state particles show sharp dip-like behavior at the invariant mass close to twice the EWIMP mass. It helps us to distinguish the EWIMP effects from backgrounds and systematic errors.
In this paper, we pursue this indirect search method further. In particular, we demonstrate that the indirect search method can be applied not only to discover EWIMPs but also to investigate their properties, such as charges, masses, and spins. To be more specific, in this paper we focus on the future prospect of the indirect studies of EWIMPs at 100 TeV colliders such as FCC-hh [20] and SppC [21,22]. We update our previous analysis [14] that has considered only the neutral current (NC) processes (mediated by photon and Z-boson) by including the charged current (CC) processes (mediated by W -boson) as well, as in Refs. [15,16]. It is crucial not only to improve the sensitivity but also to break some degeneracy among different EWIMP charge assignments; the NC and CC processes depend on different combinations of the SU (2) L and U (1) Y charges, and hence the inclusion of both processes allows us to extract these charges separately.
The rest of this paper is organized as follows. In Sec. 2, we discuss how EWIMPs affect the production processes of a charged lepton pair and those of a charged lepton and a neutrino. There we see that the EWIMP correction to the cross section, as a function of the lepton pair invariant mass, develops a dip-like structure when the invariant mass is around twice the EWIMP mass. This feature is essential in distinguishing the EWIMP effect from backgrounds and systematic errors, as discussed in detail in Sec. 3. Although we have to rely on the transverse mass instead of the invariant mass for the CC process, a similar dip-like structure appears in the transverse mass distribution. Sec. 3 is divided into three parts. First, we explain our fitting based statistical approach, in which we absorb various sources of systematic errors into a choice of nuisance parameters. Next, we study the result of the EWIMP detection reach, updating our previous results [14] by taking into account the CC processes. We then move to our main focus of this paper, namely the future prospect of the mass, charge, and spin determination of the EWIMP. Finally Sec. 4 is devoted to conclusions.

EWIMP effect on the lepton production processes
We investigate contributions of the EWIMPs to the Drell-Yan process through the vacuum polarization of the electroweak gauge bosons at the loop level. Throughout the paper, we assume that all the other beyond the SM particles are heavy enough so that they do not affect the following discussion. After integrating out the EWIMPs, the effective lagrangian is expressed as where L SM is the SM Lagrangian, D is a covariant derivative, m is the EWIMP mass, 1 g and g are the SU (2) L and U (1) Y gauge coupling constants, and W a µν and B µν are the field strength associated with the SU (2) L and U (1) Y gauge group, respectively. The function f (x) is defined as 1 Here we neglect a small mass splitting among the SU (2) L multiplet. Table 1: Coefficients of the weak interaction defined as Γ Here, e = gs W and g Z = g/c W , where s W ≡ sin θ W and c W ≡ cos θ W with θ W being the weak mixing angle.
where the first (second) line corresponds to a fermionic (scalar) EWIMP, respectively. The coefficients C 1 and C 2 for an SU (2) L n-plet EWIMP with hypercharge Y are given by where κ = 1, 2, 8, 16 for a real scalar, a complex scalar, a Weyl or Majorana fermion, and a Dirac fermion, respectively. The Dynkin index I(n) for the n dimensional representation of SU (2) L is given by which is normalized so that I(2) = 1/2. The coefficients are uniquely determined by the representation of the EWIMPs. For example, (C 1 , C 2 ) = (1, 1) for Higgsino, and (C 1 , C 2 ) = (0, 2) for Wino. We emphasize that, contrary to the usual effective field theory, our prescription is equally applied when the typical scale of the gauge boson four-momentum, q, is larger than the EWIMP mass scale m since we do not perform a derivative expansion of f in Eq. (1). It is important because, as we see soon, the effect of the EWIMPs are maximized when q 2 ∼ m, where the derivative expansion is not applicable. At the leading order (LO), we are interested in u(p)ū(p ) → − (k) + (k ) and d(p)d(p ) → − (k) + (k ) as the NC processes and u(p)d(p ) → ν(k) + (k ) and d(p)ū(p ) → − (k)ν(k ) as the CC processes. Here, u and d collectively denote up-type and down-type quarks, respectively, and p, p , k, and k are initial and final state momenta. In the SM, the amplitudes for both the NC and CC processes at the LO are expressed as where √ s is the invariant mass of the final state leptons, which is denoted as m for the NC processes and m ν for the CC processes. The relevant gauge bosons are V = γ, Z for the NC processes and V = W ± for the CC processes, with m V being the corresponding gauge boson mass. In addition, with v (V ) f and a (V ) f given in Tab. 1. The EWIMP contribution is given by where W ), and C W W = 4C 2 g 2 . Again V, V = γ, Z for the NC processes and V, V = W ± for the CC processes.
We use dΠ LIPS for a Lorentz invariant phase space factor for the two particles final state. Then, using Eqs. (6) and (8), we define where we take the average and summation over spins. Here, dL ab /d √ s is the luminosity function for a fixed √ s : where a and b denote species of initial partons, √ s is the center of mass energy of the proton collision ( √ s = 100 TeV in our case), and f a (x) is a parton distribution function (PDF) of the given parton a. Eq. (9) represents the SM cross section, while Eq. (10) the EWIMP contribution to the cross section. For the statistical treatment in the next section, we introduce a parameter µ that parametrizes the strength of the EWMP effect, and express the cross section with µ as Obviously, µ = 0 corresponds to the pure SM, while µ = 1 corresponds to the SM+EWIMP model. Hereafter, we use to denote the correction from the EWIMP. In Fig. 1, we plot δ σ for the CC processes as a function of √ s . The purple, blue, and red lines correspond to Higgsino, Wino, and 5-plet scalar, respectively. There is a dip around √ s = 2m for all the cases of the EWIMPs which originates from the loop function f in Eq. (2). The EWIMP contributions to the NC processes show a similar dip structure that  again comes from f . This dip is crucial not only for the discovery of the EWIMP signal (see Sec. 3.3) but also for the determination of the properties of the EWIMPs (see Sec. 3.4). In particular, the EWIMP mass can be extracted from the dip position, while the EWIMP charges (n and Y ) can be determined from the depth of the dip.
For the NC processes, the momenta of two final state charged leptons are measurable and we can use the invariant mass distribution of the number of events for the study of the EWIMPs. For the CC processes, on the contrary, we cannot measure the momentum of the neutrino in real experiments, and hence we instead use the missing transverse momentum p T,miss . We use the transverse mass defined as where p T, denotes the transverse momentum of the charged lepton and φ T, ,miss ≡ φ − φ miss is the difference between the azimuth angles of p T, and p T,miss . The important property of m T is that the distribution of m T peaks at m T = m ν . Because of this property, the characteristic shape of δ σ remains in the m T distribution in the CC events. To see this, we plot in Fig. 2 the EWIMP effect on the number of events as a function of m T . Here, the vertical axis is the ratio of the EWIMP correction to the number of events ∆N to the number of events in the SM N for each bin with the bin width of 100 GeV. 2 We find that the dip structure remains in the m T distribution, though the depth of the dip is smaller compared to the m ν distribution.

Event generation
Now we discuss how well we can extract information about EWIMPs from the invariant mass and transverse mass distributions for the processes of our concern at future 100 TeV pp collider experiments. We take into account the effects of the next-to-leading order QCD corrections in the events as well as detector effects through Monte-Carlo simulations.
In our analysis, we first generate the SM event sets for the NC processes pp → e − e + /µ − µ + and for the CC processes pp → e ± ν e /µ ± ν µ . We use MadGraph5 aMC@NLO (v2.6.3.2) [23,24] for the event generation with the successive use of Pythia8 [25] for the parton shower and the hadronization and Delphes (v3.4.1) [26] for the detector simulation. We use NNPDF2.3QED with α s (M Z ) = 0.118 [27] as a canonical set of PDFs. For the renormalization and factorization scales, we use the default values of MadGraph5 aMC@NLO, i.e., the central m 2 T scale after k T -clustering of the event (which we denote by Q). The events are binned by the characteristic mass m char for each process: we use the lepton invariant mass m char = m for the NC processes, and the transverse mass m char = m T for the CC processes, respectively. In both cases, we generated events with the characteristic mass within the range of 500 GeV < m char < 7.5 TeV and divide them into 70 bins with the equal width of 100 GeV.
As for the event selection by a trigger, we may have to impose some cut on the lepton transverse momentum p T . As we will see, we concentrate on events with high p T charged lepton(s) with which we expect the event may be triggered. For the NC processes, we use events with at least two high p T leptons. For our analysis, we use events with m > 500 GeV; we assume that such events are triggered by using two energetic charged leptons so that we do not impose extra kinematical requirements. On the contrary, the CC events are characterized only by a lepton and a missing transverse momentum. For such events, we require that the p T of the charged lepton should be larger than 500 GeV. 3 For the CC events, the cut reduces the number of events in particular for the bins with the low transverse mass m T ∼ 500 GeV, and thus affects the sensitivity of the CC processes to relatively light EWIMPs. We will come back to this point later.
The EWIMP effect is incorporated by rescaling the SM event by δ σ defined in Eq. (13). With the parameter µ defined in Eq. (12), the number of events corresponding to the SM+EWIMP hypothesis in i-th bin, characterized by m i,min < m char < m i,max , is where the sum runs over all the events of the final state f whose characteristic mass m char (after taking into account the detector effects) falls into the bin. Note that the true value of √ s should be used for each event for the computation of δ σ : we extract it from the hard process information. 4

Statistical treatment
We now explain the statistical method we will adopt in our analysis. We collectively denote our theoretical model as is given by Eq. (15). We denote the experimental data set asx f that in principle is completely unrelated to our theoretical model x f (µ). Since we do not have an actual experimental data set for 100 TeV colliders for now, however, we takex f = x f (µ = 1) (for some fixed values of the EWIMP mass and charges) throughout our analysis, assuming that the EWIMP does exist. In particular, this choice tests the SM-only hypothesis if we take our theoretical model as x f (µ = 0).
If the expectation values of x f,i (µ) are precisely known, the sensitivity to EWIMPs can be studied only with statistical errors. In reality, however, the computation of x f,i (µ) suffers various sources of uncertainties, which results in systematic errors in our theoretical model. The sources include errors in the integrated luminosity, the beam energy, choices of the renormalization and the factorization scales, choices of PDF, the pile-up effect, higher order corrections to the cross section, and so on. In order to deal with these uncertainties, we introduce sets of free parameters θ f = {θ f,α } (i.e. nuisance parameters) which absorb (smooth) uncertainties of the number of events, and modify our theoretical model as where f sys,i (θ f ) is a function that satisfies f sys,i (0) = 1. We expect that, if the function f sys,i is properly chosen, the true distribution of the number of events in the SM is given byx f (θ f , 0) = {x f,i (0)f sys,i (θ f )} for some value of θ f . In our analysis, we adopt the five parameters fitting function given by [30] where p i = 2m i / √ s with m i being the central value of the lepton invariant mass (transverse mass) of the i-th bin for the NC (CC) processes. As we will see, the major effects of systematic errors can be absorbed into θ f with this fitting function.
In order to test the SM-only hypothesis, we define the following test statistic [31]: Here,θ f and {θ f ,μ} are determined so that f L(x f ; θ f , µ = 0) and f L(x f ; θ f , µ) are maximized, respectively. The likelihood function is defined as where The product in Eq. (20) runs over all the bins, while the product in Eq. (21) runs over all the free parameters we introduced. For each θ f,α , we define the "standard deviation" σ f,α , which parametrizes the possible size of θ f,α within the SM with the systematic errors. If the systematic errors are negligible compared with the statistical error, we can take σ f → 0. We identify (q 0 ) 1/2 = 5 (1.96) as the detection reach at the 5σ (95 % C.L.) level, since q 0 asymptotically obeys a chi-square distribution with the degree of freedom one. In order to determine σ f , we consider the following sources of the systematic errors: • Luminosity (±5 % uncertainty is assumed), • Renormalization scale (2Q and Q/2, instead of Q), • Factorization scale (2Q and Q/2, instead of Q), • PDF choice (We use 101 variants of NNPDF2.3QED with α s (M Z ) = 0.118 [27] provided by LHAPDF6 [32] with IDs ranging from 244600 to 244700).

Sources of systematic errors
σ ee,1 σ ee,2 σ ee,3 σ ee,4 σ ee,5 Luminosity: ±5 % (σ lumi.  The values of σ f are determined as follows. Let y f be the set of number of events in the SM for the final state f with the canonical choices of the parameters, and y f be that with one of the sources of the systematic errors being varied. We minimize the chi-square function defined as for each final state f , and determine the best-fit values of θ f for each set of y f . We repeat this process for different sets of y f , and σ f are determined from the distributions of the best-fit values of θ f . For example, let us denote the best-fit values for the fit associated with the luminosity errors ±5% as θ ± f . We estimate σ f associated with these errors, denoted here as σ lumi. f , as where N denotes the number of fitting procedures we have performed: N = 2 for this case.
We estimate σ f associated with the other sources of the errors, denoted as σ ren. f , σ fac. f , and σ PDF f , in a similar manner. Finally, the total values of σ f are obtained by combining all the sources together as In Table 2 and 3, we show the values of σ ee and σ eνe associated with each source of the systematic errors, respectively. These values can be interpreted as the possible size of the fit parameters within the SM, which is caused by the systematic uncertainties. As explained in Eq. (25), we combine these values in each column to obtain σ f . In Table 4, we summarize the result of the combination for all the final states. The values of σ f are independent of the final state lepton flavors since the energy scale of our concern is much higher than the lepton masses. However, we use different sets of fit parameters θ ee and θ µµ for the NC processes and θ eνe and θ µνµ for the CC processes because of the different detector response to electrons and muons.

Detection reach
Now we show the detection reach of EWIMPs at future 100 TeV colliders. In Fig. 3, we plot the value of √ q 0 as a function of the EWIMP mass, with the integrated luminosity L = 30 ab −1 . As representative scenarios, we show the cases for Higgsino (the red lines) and Wino (the blue lines). The dotted and dash-dotted lines are the result obtained only from the NC processes and the CC processes, respectively. We find that the CC processes are more sensitive to the effect of the EWIMPs than the NC processes because of the larger cross section. This result is consistent with Refs. [15,16]. The sensitivity of the CC processes is weakened for m 700 GeV because of the lepton p T cut we have applied. 5 The combined results of the NC and CC processes are shown by the solid lines. By combining the two types of processes, the 5σ discovery reaches (95 % C.L. bounds) for Higgsino and Wino are 850 GeV (1.7 TeV) and 1.3 TeV (2.3 TeV), respectively. We find that the combination of the NC and CC processes improves the sensitivity of the EWIMP mass. Furthermore, if we understand all the systematic uncertainties quite well and effectively take the σ f → 0 limit in the combined result, the detection reach will be pushed up significantly as shown by the dashed lines: 1.1 TeV Higgsino signal at well above 5σ level and a 4σ hint of the 2.9 TeV Wino. Therefore, it is essential to reduce the systematic uncertainties for the detection of EWIMPs through the NC and CC processes.

Determination of EWIMP properties
In this subsection, we show that it is possible to determine the properties of the EWIMPs from the NC and CC processes, thanks to the fact that we can study the m and m T distribution in great detail for these processes. Some information about the mass, charge, and spin of the EWIMPs can be extracted because the corrections to these distributions from the EWIMPs are completely determined by these EWIMP properties. Firstly, we can extract the EWIMP mass from the position of the dip-like structure in the correction since it corresponds to roughly twice the EWIMP mass as we have shown in Sec. 2. Secondly, the overall size of the correction gives us information about the SU (2) L and U (1) Y charges. The CC processes depend only on the SU (2) L charge, while the NC processes depend both on the SU (2) L and U (1) Y charges. Consequently, we can obtain information about the gauge charges of the EWIMPs from the NC and CC processes.
We now demonstrate the mass and charge determination of fermionic EWIMPs. This is equivalent to the determination of the parameter set (m, C 1 , C 2 ). We generate the data assuming the SM + EWIMP model (µ = 1) with some specific values of m, n, Y , and κ, with which we obtain (m, C 1 , C 2 ). We fix µ = 1 for our theoretical model as well, and hence the theoretical predictions of the number of events also depend on these three parameters, . We define the likelihood function L(x f ; θ f , m, C 1 , C 2 ) in the same form as Eqs. (16) and (19) with the theoretical prediction x f now understood as a function of Higgsino signal. The dotted and solid lines denote 1σ and 2σ contours, respectively, and the gray region corresponds to the parameter space that is in tension with the observation at more than 2σ level. The blue, green, and red lines correspond to the result from the NC processes, the CC processes, and the combined analysis, respectively. Each star marker annotated as "n Y " represents a point corresponding to a SU (2) L n-plet Dirac fermion with hypercharge Y , while that with "n Maj " corresponds to an SU (2) L n-plet Majorana fermion.
(m, C 1 , C 2 ), not of µ. 6 The test statistic is defined as where the parameters ({θ f },m,Ĉ 1 ,Ĉ 2 ) maximize f L(x f ; θ f , m, C 1 , C 2 ), whileθ f maximize L(x f ; θ f , m, C 1 , C 2 ) for fixed values of (m, C 1 , C 2 ). It follows the chi-squared distribution with three degrees of freedom in the limit of a large number of events [33]. The test statistic defined in this way examines the compatibility of a given EWIMP model (i.e. a parameter set (m, C 1 , C 2 )) with the observed signal. Once a deviation from the SM prediction is observed in a real experiment, we may determine (m, C 1 , C 2 ) using the above test statistic q. In the following, we show the expected accuracy of the determination of (m, C 1 , C 2 ) for the case where there exists 1.1 TeV Higgsino. 7 In Fig. 4, we show the contours of 1σ (dotted) and 2σ (solid) constraints, which correspond to the values √ q = 1.9 and √ q = 2.8, respectively, in the C 1 vs. C 2 plane for m = 1.1 TeV. The blue, green, and red lines denote the result obtained from the NC processes, the CC processes, and the combined analysis, respectively. The models in the gray region are in more than 2σ tension with the observation. We also show several star markers that correspond to the single SU (2) L multiplet contributions: the markers with "n Y " represent an SU (2) L n-plet Dirac fermion with hypercharge Y , while those with "n Maj " an SU (2) L n-plet Majorana fermion. Both the NC and CC constraints are represented as straight bands in the C 1 vs. C 2 plane since each process depends on a specific linear combination of C 1 and C 2 . In particular, the CC constraint is independent of C 1 , or Y . In this sense, the NC and CC processes are complementary to each other, and thus we can separately constrain C 1 and C 2 only after combining these two results. For instance, we can exclude a single fermionic SU (2) L multiplet with n = 2 at more than 2σ level, although each process by itself cannot exclude the possibility of 3 Maj . We can also constrain the hypercharge, yet it is not uniquely determined. In addition to the Higgsino, the EWIMP as an SU (2) L doublet Dirac fermion with |Y | 2 2 or an SU (2) L doublet Majorana fermion with |Y | 2 5 is still allowed. In Again, by combining the NC and CC results, we can significantly improve the determination of EWIMP properties, making 1σ and 2σ contours closed circles in the planes of our concern. In particular, as red lines show, the combined analysis allows us to determine the observed EWIMP mass at the level of O(10)%.
Finally, we comment on the possibility of discriminating between fermionic and scalar EWIMPs, whose difference comes from the loop function f (x) (see Eq. (2)). Here we repeat the same analysis explained above, assuming the 1.1 TeV Higgsino signal for example, but use the scalar loop function to evaluate the theoretical predictions x f (m, C 1 , C 2 ). In Figs. 6 and 7, we show the results in the C 1 vs. C 2 plane and the C 1 (or C 2 ) vs. m plane, respectively, where one of the three parameters is fixed to its best fit value. It is seen that, in the case of the 1.1 TeV Higgsino signal, it is hard to distinguish between the bosonic and fermionic EWIMPs only with our method. However, if a part of the EWIMP properties (in particular its mass) is determined from another approach, our method may allow us to determine its spin correctly. We also stress here that, with some favorable assumption about the observed signal, we may obtain some hint about its spin. For example, if we assume that the observed signal composes a fraction of the dark matter in our Universe, the choice of the EWIMP charges is significantly constrained. Note from Fig. 6 that the only choices of EWIMP charges that allow the EWIMP multiplet to contain an electrically neutral component are (n, |Y |) = (3, 0), (3, 1), (4, 1/2), (4, 3/2), and (5, 0) real . The last column of the table 5 shows proper choices of EWIMP masses in order for their thermal relic abundances become comparable with the dark matter abundance in the current Universe. All of those values are somewhat larger than the central value of the mass of the observed signal, which means that the scalar interpretation of the signal cannot explain the whole of the dark matter relic abundance without introducing some non-thermal production mechanism.

Conclusion
In this paper, we have discussed the indirect search of EWIMPs at future 100 TeV hadron colliders based on the precision measurement of the production processes of a charged lepton pair and that of a charged lepton and a neutrino. In particular, we have demonstrated that not only we can discover the EWIMPs, but also we can determine their properties such as their masses, SU (2) L and U (1) Y charges, and spins via the processes of our concern. It is based on two facts: the high energy lepton production channel enables us to study its momentum distribution in great detail, and the EWIMP correction shows characteristic features, including a dip-like structure as the final state invariant mass being twice the EWIMP mass. The latter feature also helps us to distinguish the EWIMP signals from backgrounds and systematic errors, as they are not expected to show a dip-like structure. In order to fully exploit the differences between the distributions the EWIMP signals and systematic errors, we have adopted the fitting based analysis as our statistical treatment. First, we have shown in Fig. 3 the detection reach of Higgsino and Wino from the neutral current (NC) processes (mediated by photon or Z-boson), the charged current (CC) processes (mediated by W -boson), and the combination of these two results. We have seen that the addition of the CC processes improves the detection reach from the previous analysis [14]. From the combined analysis, the bounds at the 5σ (95% C.L.) level for Higgsino and Wino are 850 GeV (1.7 TeV) and 1.3 TeV (2.3 TeV), respectively. This result, in particular that for short lifetime Higgsino, indicates the importance of our method for the EWIMP search.
Next, we have considered the determination of the mass and SU (2) L and U (1) Y charges of the observed EWIMP. By combining the NC and the CC events, the position and the (n, Y )   Fig. 6. The observed DM energy density is explained by the thermal relic of the EWIMP with m DM shown in the fourth column.
height of the dip in the EWIMP effect on the cross section gives us enough information for determining all the three parameters. In Figs. 4 and 5, we have shown the plots of the test statistics that test the validity of several choices of parameters. As a result, the SU (2) L charge of the observed signal is correctly identified under the assumption of a single EWIMP multiplet, and the U (1) Y charge and mass are also determined precisely. In order for the determination of the EWIMP spin, we have plotted the contours of the test statistics that test the validity of the scalar EWIMP models with some fixed values of masses and charges. The results are shown in Figs. 6 and 7, which reveals that the spin is not completely determined by solely using our method. Use of another approach to determine the EWIMP properties, or of some assumption like that the observed signal corresponds to the dark matter in our Universe, may help us to obtain further information regarding the EWIMP spin.