Distinguishing spins in decay chains with photons at the Large Hadron Collider
Abstract:
Several models for physics beyond the Standard Model predict new particles with a decay signature including hard photons and missing energy. Two wellmotivated examples are supersymmetry with gaugemediated breaking (GMSB) and the standard model with two universal extra dimensions. Both models lead to decay chains with similar collider signatures, including hard photon emission. The main discriminating feature are the spins of the new particles. In this paper we discuss how information about the spins of the particles can be extracted from leptonphoton or quarkphoton invariant mass distributions at the Large Hadron Collider. The characteristic shapes of the distributions are derived analytically and then studied in a realistic MonteCarlo simulation. We find that for a typical GMSB mass spectrum with particle masses below 1 TeV, already 10 fb integrated luminosity at 14 TeV centerofmass energy are sufficient to discriminate the two models with high significance.
1 Introduction
Several models for physics beyond the Standard Model (SM) introduce partner particles for all SM particles, the lightest of which is stable and could be the constituent of dark matter. The best known examples are supersymmetry (SUSY) with conserved Rparity and universal extra dimensions (UED) with conserved KaluzaKlein (KK) parity [1]. Since the observable signatures for these models at the Large Hadron Collider (LHC) look quite similar [2], it will be important to test the fundamental quantum numbers of the new particles in order to scrutinize the nature of the underlying physics.
For example, both SUSY and UED require that the couplings of the new particles are identical to the corresponding couplings of their SM partners; a prediction which can be tested at the LHC by measuring cross section ratios [3]. However, a crucial distinction between the two models is given by the spins of the new particles. While the SUSY partners differ from their SM counterparts by one halfunit of spin, the KK excitations in UED have the same spin as their SM partners. Recently, extensive work has been performed to determine the spins of SUSY or UED particles by exploiting angular correlations in the decay of those particles at the LHC^{1}^{1}1A few studies have explored more modeldependent discrimination methods based on total cross sections [4] and higher KK modes [6].. Many papers have focused on decay chains involving lepton pairs [5, 6, 7, 8, 9, 10, 11, 12]. A typical example of such a decay chain in SUSY is . Other studies have examined channels involving heavy gauge bosons [9, 13, 14, 15], sleptons [14, 16] and top quarks. However, in all of the existing studies it was assumed that the lightest new particle is either a neutralino or weak vector boson in SUSY or UED, respectively.
In this article the spin determination from angular correlations is extended to decay chains that involve hard photons in the final state. Such decay channels occur naturally in gaugemediated supersymmetry breaking (GMSB), where the lightest SUSY particle is the gravitino [18], as well as in the extension of the SM by two universal extra dimensions (UED6), where the lightest KK mode is typically a scalar component of a higherdimensional vector boson, called “scalar adjoint” [19, 20, 21]. It has been shown earlier that a highenergy collider could distinguish between GMSB and UED6 by studying angular correlations in pair production and decay of the selectron (KKelectrons) and neutralinos (KKgauge bosons) [22]. The purpose of the present paper is to study how such a distinction can be achieved at the LHC by analyzing decay chains involving leptons and photons. In particular, we are investigating how the spin of the lightest new particle, which escapes from the detector in form of missing momentum, can be inferred from invariant mass distributions of the leptons and photons.
In section 2 we describe analytical calculations of the relevant invariant mass distributions and compare the angular correlations predicted by GMSB and UED6. In order to evaluate the prospects for experimental measurements of these distributions, we present in section 3 results of a realistic MonteCarlo simulation, incorporating the spin correlation effects of the two models. We present our conclusions in section 4. The detailed analytical results are collected in the appendix.
2 Spin correlations in GMSB and UED6
In our notations and conventions we follow Ref. [23] for supersymmetry and Ref. [19, 20, 21] for the Standard Model in six dimensions (UED6).
GMSB is a promising candidate for a mechanism that generates TeVscale masses for the SUSY partners, most notably since it explains the absence of large flavorchanging neutral currents. In GMSB, the lightest SUSY particle is typically the gravitino , with a mass MeV. If the nexttolightest SUSY particle is a neutralino, it will mostly decay into the gravitino and a photon, . Depending on the neutralino mass, there could be a smaller branching fraction into a boson, which we will not investigate further.
In UED6, the lightest particle at KK level (1,0) is typically the scalar adjoint of the hypercharge boson, [20, 24]. The vector mode of the KKhypercharge boson, , can decay into the scalar adjoint via a loopinduced process, . This decay mode has a sizable branching fraction of about 34% [21]. Therefore the two models lead to very similar decay signatures. In particular, the decay chain
(1)  
which is typical in GMSB, is imitated by the equivalent process in UED6,  
(2) 
see Fig. 1. Both processes lead to a final state signature of a sameflavor, oppositesign lepton pair, one photon, and missing transverse momentum .
In the following we will study spin correlation effects in these decay chains in detail. We will also briefly analyze the following shorter decay modes of squarks:
(3)  
for GMSB, and  
(4) 
for UED6. Although these decay channels have considerably larger SM backgrounds at the LHC, they also have larger branching ratios compared to the decay channels with leptons. They could be very useful to obtain information on the spin of the quark partners.
While both GMSB and UED6 have similar mass hierarchies, which allow the decay channels in eqs. (1)–(4), the typical mass spectra are quantitatively rather different.
In minimal gauge mediation the gaugino mass parameter relation implies that the weak gauginos are rather light, while the gluino is much heavier. Furthermore, the squarks are also heavy, while the sleptons have masses close the the gauginos. For our numerical analysis in the next section we will use the reference scenario G1a from Ref. [25]. The masses of the particles appearing in our decay chain and their branching ratios are summarized in Tab. 1.
G1a  U1  

Particle  Mass [GeV]  Particle  Mass [GeV] 
747  696  
986  662  
989  
942  608  
939  606  
224  538  
119  487  
326  521  
164  508  
0  427  
Branching ratios  Branching ratios  
BR[]  16%  BR[]  50% 
BR[]  6.4%  
BR[]  26%  BR[]  1.5% 
BR[]  100%  BR[]  34% 
In universal extra dimensions, on the other hand, the masses of all particles of one KK level have the same value , where is the size of the extra dimensions. This degeneracy is lifted only by radiative corrections [20, 24], which shift the masses by up to 20%. As a result, a typical spectrum for the particles of KK level 1 in UED6 is much less hierarchical than in GMSB. In our numerical analysis, we will use the reference scenario U1, defined by GeV, with the masses and branching ratios [21] given in Tab. 1.
These qualitative features of the spectra could be used to distinguish GMSB and UED6 experimentally. However, there are several caveats to consider: The mass spectra of GMSB can vary substantially in nonminimal models, see e. g. Ref. [26]. Extra dimensional models are known to become strongly coupled at large energies and require some new physics to be present at the scale. The effects of this unknown highscale physics could generate mass contributions to the KK particles in UED [27]. However, the spins of the the new particles can serve as very robust discriminators between different models.
A nonzero spin of a particle can lead to angular correlations between its decay products. At the LHC, angular correlations are manifested in the invariant mass distributions of the visible decay products of a decay chain. The long decay chains (1) and (2) are of the general form
(5) 
with . Here we call the lepton that is emitted in the first decay step the “nearlepton” , while the lepton from the second decay step is named the “farlepton” . From this one can construct the invariant masses
(6)  
(7)  
(8)  
(9) 
which are related to , the angle between the nearlepton and the farlepton in the rest frame of , , the angle between the nearlepton and the photon in the rest frame, and , the angle between the farlepton and the photon in the rest frame, respectively. The maximum values for the invariant masses are given in eq. (30) in the appendix.
For a given decay matrix element, the distribution with respect to some invariant mass is then obtained by integrating over all remaining kinematical variables in a given reference frame, as described in detail in Refs. [7, 28].
In practice, the near and far leptons cannot be distinguished in a straightforward way. Instead the observable leptonphoton invariant mass distribution is the sum of . Additional information can be obtained from the distributions with respect to the minimum and maximum of the leptonphoton invariant masses,
(10) 
Since the total magnitude of the decay width does not carry any information about the spins of the particles involved, we will normalize the invariant mass distributions to unity,
(11) 
where is the integrated decay width of the given decay channel, and is defined as a differential probability density.
For the short decay chains (3) and (4) of the general form
(12) 
with , the only observable invariant mass distribution that can be constructed is
(13) 
Here .
2.1 Spin correlations in GMSB
In GMSB the final state is fed by the decay chain eq. (1), with and as intermediate particles. Since the is a scalar it does not transmit any angular correlations. The fermionic can lead to nontrivial spin correlations in the decay chain, but only if the couplings at both the production and decay vertex are chiral, i.e. left and righthanded components have different coupling strength [9]. While this condition is fulfilled for the  vertex, the  vertex has the form [29]
(14) 
Here the left and righthanded component of the neutralino couple with equal strength so that no observable angular correlations are generated. Thus the only features in the invariant mass distributions of the GMSB decay chain (1) are generated by the phase space. Analytical results for the leptonphoton and leptonlepton distributions are listed in the appendix.
Similarly, the short decay chain eq. (3) does not lead to visible angular correlation effects. The quarkphoton invariant mass distribution is given by
(15) 
By comparison with the formulas in the appendix one can see that eq. (15) is identical to the distribution of the long chain. This can be easily understood by the fact that the chirality of the slepton and squark couplings is identical.
2.2 Spin correlations in UED6
The typical mass hierarchy generated by radiative corrections in UED6, , enables the decay chain eq. (2). However, in general sizeable corrections to the KKparticle masses could be generated by the unknown physics that complete the theory at high energies [27]. Thus for completeness we will study all possible decay chains that, for suitable mass hierarchies, could lead to the final state , as illustrated in Fig. 2,
(VFVS)  (SFVS) 
(VFSV)  (SFSV) 
(16)  
(17)  
(18)  
(19) 
Here we have introduced shorthand notations for the four decay chains based on the KK particles at each decay stage being a scalar (S), fermion (F) or vector (V). In all cases we keep the couplings of the KK particles as they are predicted by the UED6 model, and mixing between gauge eigenstates is neglected. Analytical results for the invariant mass distributions for all four combinations are listed in the appendix.
In case of the short decay chain, eq. (4), there are two possible decay chains with the same final state, depending on the mass hierarchy,
(20)  
(21) 
The quarkphoton invariant mass distributions for the two cases read
(22)  
(23) 
As before, one can see that these are identical to the distributions of the long chain for the VFVS/SFVS and VFSV/SFSV combinations, respectively. In the FSV case, since the intermediate particle is a scalar, this decay process does not involve any spin correlations and is identical to the pure phase space distribution, and thus to the GMSB distribution.
2.3 Discussion of analytical results
Figs. 3 and 4 show the distributions for the four independent observable invariant mass combinations of the final state: the dilepton invariant mass , the “low” and “high” leptoninvariant masses and , respectively, and the leptonleptonphoton invariant mass . Each plot contains five curves corresponding to the five models (or spin assignments) GMSB, VFVS, SFVS, VFSV, and SFSV. In case of Fig. 3, for all five models the masses have been chosen from the G1a scenario, with , , , and . On the other hand, Fig. 4 shows the situation for the U1 spectrum with , , , and .
As evident from the plots, different distributions could discriminate between different spin assignments. The dilepton distribution is markedly different for the VFSV and SFSV models, where the particle is a scalar, compared to the other models. This can be understood from the fact that the chiral structure of the KKfermion couplings to KKscalars or KKvector bosons is different (see for example Ref. [10] for more details). On the other hand, the peak in the and spectra is relatively enhanced for the VFSV and SFVS models. As a result, the discriminative power between different spin assignments is maximized by including all four distributions in the analysis.
Fig. 5 shows the jetphoton invariant mass distribution for the short decay chain. For the distribution on the left hand side the masses have been chosen from the G1a scenario, , , and on the right hand side the distribution is shown for the U1 scenario, where , , . As mentioned before, the GMSB and FSV chains both do not generate any spin correlations and thus cannot be distinguished from each other. However, the plots show a discrimination potential between FVS and the other two cases, especially for the G1a mass scenario. The qualitative features of the jetphoton distribution are similar for the G1a and U1 scenarios, although the spin correlation effects are less pronounced for the more degenerate U1 mass spectrum than for the more hierarchical G1a spectrum.
3 MonteCarlo simulation and numerical analysis
In the previous section theoretical formulas for the invariant mass spectra of decay chains in GMSB and UED6 were discussed. However, in a realistic experimental setup one must take into account several effects that influence the measured distributions.

Detector effects, such as limited resolution, and reconstruction effects from particle identification can smear out the visible invariant mass distributions. Further, detector acceptance can distort the shape of the distribution.

For the separation of the new physics signal from the SM backgrounds suitable selection cuts have to be implemented, which can distort the shape of the distributions.

Since both SUSY and KK particles are produced in pairs, often there will be two hard photons in one event, leading to a twofold ambiguity in the reconstruction of the correct decay chain.
Therefore we have performed a realistic experimental simulation of the new physics signal from the long decay chain for GMSB and UED6, with the goal of reconstructing the spinsensitive invariant mass distributions and discriminating between models from the simulated data.
The decay chains in eqs. (1), (16)–(19) have been computed with CompHEP 4.4 [30] for 14 TeV protonproton collisions, using the GMSB model file from Ref. [31] and our own implementation of the UED6 model. Partonlevel events generated with CompHEP were then passed on to Pythia 6.4.12 [32]. Finally, the ATLAS detector[33] is simulated using the parametrized fast detector simulation ATLFAST[34, 35], which includes detector acceptance, resolution and some basic particle identification.
For the G1a and U1 mass spectra, the aforementioned decay chains are initiated mostly by squark (KKquark) and gluino (KKgluon) production processes, respectively. Since our analysis does not depend crucially on the details of the hadronic decay products of the squarks (KKquarks) and gluinos (KKgluons) we have only generated events for squark (KKquark) pair production as the primary hard process, and then normalized the total event count according to the total production cross section including gluino (KKgluon) processes.
As a further simplification, we have only generated matrix elements with the decay chain of one of the two squarks (KKquarks) within CompHEP, which correctly implements all spin correlations. The decay of the second squark (KKquark) was simulated in Pythia, without spin correlations. Nevertheless, this procedure provides a good approximation to the complete matrix elements since the branching ratio for the decay chain leading to the final state is relatively small and thus only a very small fraction of the events contains two decay cascades of this type. Moreover, Pythia was used for simulating initial and final state radiation and hadronization.
The event selection was performed according to Ref. [25]. First the effective mass
(24) 
is defined, where denotes the missing transverse energy and the transverse momenta of the 4 hardest jets. The selected events need to fulfill the following conditions:

4 jets with transverse momenta ,

,

,

2 photons with transverse momenta ,

2 electrons or muons with transverse momenta .
After application of these cuts the SM background is reduced to a negligible level [25], while about 20% of the signal is retained.
Since the selected events include two photons, one of them has to be selected to compute the invariant mass distributions. Good results are obtained, when choosing the one photon, which gives the smaller invariant mass.
It would be interesting to also analyze the short decay chain in a complete simulation, since the expected rates are large (about 5 pb both for G1a and U1). However, the signature of this final state, two hard jets, two photons, and missing energy, is very sensitive to issues related to jets faking photons. Therefore, it would require a more careful analysis of QCD backgrounds, which will be left for a future publication.
3.1 Cross sections
For a realistic analysis, the modeldependent crosssections and event numbers have to be calculated. At the G1a point, the total squark and gluino production cross section is [25]. Since the gluino is lighter than the squarks, the squarks decay through cascades involving a gluino. With the branching ratios from Tab. 1 a good approximation for the cross section for the decay chain is
(25) 
The factors of result from the two generations of leptons in the dileptonic decay chain and the fact that squarks or gluinos are produced in pairs and both of them can decay through this channel.
In the U1 model, the KKquarks and KKgluons are produced with the following cross sections [21]
(26) 
In this scenario, KKquarks are lighter than the KKgluon, such that the latter will mostly decay into KKquarks. Therewith, the total KKquark production cross section is
(27)  
As above, factors of 2 account for the two sides of the pair production process. Then the total cross section for the decay chain adds up to
(28) 
Using the signal efficiency of the cuts in Ref. [25] and assuming an integrated luminosity of 10 fb we obtain selected events for the G1a point and selected events for the U1 model.
It should be noted that the event rates depend strongly on the underlying model and its parameters, as well as the choice of selection cuts. Furthermore, the signal efficiency after cuts depends also on the decay chain of the second squark or KKquark, whose branching ratios vary between different models. For the purpose of this study, we do not vary the choice of cuts, cross sections and signal efficiency between different models when we compare them for one given mass spectrum (i. e. for all spin assignments we assume 2500 selected events for the G1a mass spectrum and 250 selected events for the U1 mass spectrum). Rather, our numerical analysis of the two scenarios should only serve as concrete examples for a spin determination of a new physics signal, in particular since our method does not rely on information about total rates.
3.2 analysis
In order to discriminate the histograms for the different spin configurations we used the test implemented in ROOT [36]. It returns the probability, i. e. the probability that two histograms with identical underlying distribution functions have a bigger value than the two compared ones. Since these values depend on the number of bins, the discrimination was performed with 5 bins, which showed the best discriminative power at . In Tab. 2 the minimal probabilities for each pair of spin configurations are listed.
GMSB  VFVS  VFSV  SFVS  SFSV  

GMSB  0.000 ()  0.000 ()  0.006 ()  0.000 ()  
VFVS  0.056 ()  0.000 ()  0.000 ()  0.000 ()  
VFSV  0.577 ()  0.155 ()  0.000 ()  0.000 ()  
SFVS  0.025 ()  0.065 ()  0.084 ()  0.000 ()  
SFSV  0.000 ()  0.000 ()  0.000 ()  0.000 () 
The results of the probabilities reflect the general features that can be seen in the histograms of Fig. 6 and 7. As expected, the discrimination between different spin combinations is far more effective for the G1a scenario than for the U1 scenario, as a result of the more degenerate mass spectrum and lower cross section in the latter case. For the G1a scenario, even with 10 fb almost all models can be distinguished with a confidence level better than 99.9%. In the case of the U1 scenario only the SFSV spin assignment can be distinguished at this confidence level with 10 fb luminosity, owing to the distinctly different shape of the distribution of the dilepton invariant mass . In other cases however, in particular for the GMSB spin assignment and the extradimensional chain VFSV, it is not possible to make a distinction even at the 95% confidence level. Here a much larger integrated luminosity would be required for a significant discrimination.
As an example Fig. 8 shows the reconstructed invariant mass distributions corresponding to G1a masses for an integrated luminosity of 30 fb. The distributions are divided in 10 bins, which gives the best discriminative power for this luminosity. More shape information is available due to more bins and events allowing for a better separation between the different models. For the U1 mass scenario the GMSB spin assignment can be distinguished from the extradimensional chain VFSV at almost 90% confidence level.
Nevertheless, in summary the analysis of invariant mass distributions proves to be a powerful tool to identify the spins of particles in decay chains with hard photons. In this analysis the form the of the heavy particle interactions was fixed to the GMSB and UED6 predicitions. In principle, the shapes of the invariant mass distributions are also sensitive to the couplings of those particles [12], but a more general analysis, where these parameters are taken as free variables, will be left for future work.
4 Conclusion
In this article we have analyzed the prospects for determining the spins of new particles in decay chains with photons and missing energy at the LHC. As concrete model incarnations of such signatures we considered supersymmetry with gaugemediated breaking (GMSB) and the standard model with two universal extra dimensions (UED6). Each of these models predicts new partners of the SM particles, the lightest of which is stable on grounds of a conserved parity. The supersymmetric or KK partner of the photon can decay into this stable particle and a hard photon. At the LHC the partners of the colored particles (squarks/gluinos or KKquarks/KKgluons) are produced with large cross sections and subsequently could decay in several steps until the final photon emission step. As became well known in recent years, the distributions of the invariant masses of two or more of the decay products is sensitive to the spins of the decaying particles. The measurement of the spins of the new particles is of central interest since this is the crucial difference between supersymmetry and extra dimensions.
Following this approach, we first derived analytical expressions for the invariant mass distributions for characteristic decay chains of those two models. In greater detail we analyzed a class of decay chains that lead to two leptons, one photon, and missing energy. This signature can stem from five different decay processes (one in GMSB, four in UED6), which differ by the spins of the intermediate particles. It was found that different invariant mass distributions show different distinctive features between those five cases, so that the discriminative power is maximized by combining information from all distributions. In a second step, we performed a realistic phenomenological MonteCarlo analysis including a fast detector simulation for those five processes. Due to detector effects and cuts for background rejection, the reconstructed invariant mass distributions are distorted compared to the analytical partonlevel results. Nevertheless, the essential characteristic features of the five different spin assignments are preserved, so that the spins of the intermediate particles can be determined.
For a typical GMSB mass spectrum with superpartner masses below 1 TeV, we found that with 10 fb integrated luminosity almost all of the five different spin assignments, and in particular the GMSB and UED6 models, can be distinguished with 99.9% confidence level. A typical UED6 spectrum, on the other hand, is more degenerate, leading to smaller branching ratios to leptons and to suppressed spin correlation effects in the invariant mass distributions. As a result, for such a mass spectrum, GMSB and UED6 can only be distinguished with a confidence level of less than 50% with 10 fb. A much higher luminosity would be needed for better discrimination in this scenario.
Additional information could be obtained by looking at other decay chains. We briefly investigated decay chains leading to a hard jet, a photon, and missing energy, but no leptons. Such a decay chain has a large branching fraction for UED6 scenarios and thus might be useful for discriminating models in this case. A more conclusive answer would require a detailed analysis of backgrounds to this process, which is beyond the scope of this work.
In summary, it appears feasible to distinguish GMSB and sixdimensional UED models with LHC data alone if the mass spectrum is not very degenerate.
Acknowledgments.
This work was supported in part by the Schweizer Nationalfonds. A. F. is grateful for warm hospitality at the Universität Zürich, Argonne National Laboratory and the University of Chicago, where part of this work was performed. W. E. is grateful for warm hospitality at the Universität Zürich, where part of this work was performed. We thank members of the ATLAS Collaboration for helpful discussions. We have made use of the ATLAS physics analysis framework and tools which are the result of collaborationwide efforts.Appendix A Analytical results for invariant mass distributions
In the appendix we list analytical results for the invariant mass distributions of the various GMSB and UED6 decay chains discussed in section 2. All those decay chains are of the form
(29) 
with . Owing to kinematical constraints, the invariant mass distributions are divided into sections, which are bounded by the kinematical edges
(30) 
Therefore the distributions take on the following structures:
The farlepton–photon mass distribution :
(31) 
The nearlepton–photon mass distribution :
(32) 
The high lepton–photon mass distribution :
Hierarchy A11:
(33) 
Hierarchy A12:
(34) 
Hierarchy A2:
(35) 
Hierarchy B1:
(36) 
Hierarchy B2:
(37) 
The low leptonphoton mass distribution :
Hierarchy A1:
(38) 
Hierarchy A2:
(39) 
Hierarchy B:
(40) 
Many of these coefficients are related:
(41) 
Below the results for the independent coefficients are given:
GMSB = phase space:
(42)  
(43)  
(44)  
(45)  
(46)  
(47) 