Photons, photon jets, and dark photons at 750 GeV and beyond

In new physics searches involving photons at the LHC, one challenge is to distinguish scenarios with isolated photons from models leading to “photon jets”. For instance, in the context of the 750 GeV diphoton excess, it was pointed out that a true diphoton resonance S→γγ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S \rightarrow \gamma \gamma $$\end{document} can be mimicked by a process of the form pp→S→aa→4γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p p \rightarrow S \rightarrow a a \rightarrow 4\gamma $$\end{document}, where S is a new scalar with a mass of 750 GeV and a is a light pseudoscalar decaying to two collinear photons. Photon jets can be distinguished from isolated photons by exploiting the fact that a large fraction of photons convert to an e+e-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^+e^-$$\end{document} pair inside the inner detector. In this note, we quantify this discrimination power, and we study how the sensitivity of future searches differs for photon jets compared to isolated photons. We also investigate how our results depend on the lifetime of the particle(s) decaying to the photon jet. Finally, we discuss the extension to S→A′A′→e+e-e+e-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S\rightarrow A^\prime A^\prime \rightarrow e^+e^-e^+e^-$$\end{document}, where there are no photons at all but the dark photon A′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A^\prime $$\end{document} decays to e+e-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^+e^-$$\end{document} pairs. Our results will be useful in future studies of the putative 750 GeV signal, but also more generally in any new physics search involving hard photons.


Introduction
In their recent end-of-year jamboree, the ATLAS and CMS collaborations have reported an impressive cornucopia of LHC Run II results. One of them -a possible excess in the two photon final state at an invariant mass ∼750 GeV [1,2] -has caused a flurry of discussion in the community [3]. Most of these works introduce a new neutral scalar particle φ with a mass around 750 GeV and decaying to two photons. Both the production and the decay of this particle typically proceed through loop diagrams. Constraints from Run I data a e-mail: bdasgupta@theory.tifr.res.in b e-mail: jkopp@mpi-hd.mpg.de c e-mail: pedro.schwaller@desy.de imply that the production cross section of φ must be significantly larger at the Run II center-of-mass energy of 13 TeV than at the Run I energy of 8 TeV. Moreover, decay modes of φ other than φ → γ γ should not be too strong, but at the same time, φ should have a large total width ∼45 GeV to optimally fit the data.
A very appealing class of alternative models explaining the 750 GeV excess are those in which the final state is in fact not two body, but contains two "photon jets", i.e. groups of highly collinear photons [4][5][6][7][8][9][10]. If the photon jets are sufficiently collimated, they are indistinguishable from isolated photons using information from the electromagnetic calorimeter alone. Therefore, models of this type could explain the diphoton anomaly, as discussed in Refs. [11][12][13][14][15][16]. While the experiments have strong discriminating variables to reject e.g. photon pairs coming from neutral hadron decays, the studies [11][12][13][14][15][16] show that there are regions of parameter space where the photon jets are expected to pass the tight photon selection.
However, there is a catch: since photons have to travel through some amount of detector material before reaching the calorimeter, they have a high (e.g. ∼40 % at ATLAS [17]) probability of converting to an e + e − pair already in the inner detector, with nontrivial pseudorapidity dependence (see Fig. 1). Such conversions can occur in the strong electric field of an atomic nucleus through a process γ + Z → Z + e + + e − . "Converted photons" are routinely included in analyses involving photons.
For a highp T photon jet with ≥2 photons, it is clear that the probability that at least one of the photons inside the jet converts is higher than for isolated photons. Even if no further discrimination is performed, we will show below that the ratio of converted to unconverted photon events already provides a powerful discrimination between isolated photon and photon jet models. Furthermore this ratio can also be used to improve the sensitivity of searches for photon jet events. Going beyond conversion ratios, several other observables could be used to reveal the photon jet origin of signals involving photons, including non-resonant photons. This includes a mismatch of the track p T and the calorimeter E T if only one photon converts, a non-standard response of the signal to changes of the photon selection criteria, and converted photon candidates with more than two tracks.
In the rest of this note, we will first discuss the use of converted photon ratios to discriminate events with photon jets from isolated photons and to improve the sensitivity of searches for such models (Sects. 2 and 3). After that we will analyze the effects of finite lifetime of the intermediate states on this analysis (Sect. 4), and we will extend the discussion to models with dark photons decaying directly to displaced e + e − pairs (Sect. 5). Finally, we discuss the prospects of other observables in more detail (Sect. 6). While most of our numerical results are obtained using ATLAS 13 TeV data, we expect that at CMS similar results can be expected, since the conversion rate is similar in magnitude and rapidity dependence, as seen from Fig. 1.

Photon jets
Before digging into the details, let us first review the type of models that can give rise to photon jets and which therefore can be probed by the methods we present below. Any particle that decays to two or more photons can produce a photon jet if it is sufficiently boosted. Consider a particle a with mass m a and Lorentz boost γ = E/m a decaying to two photons. The minimal opening angle between the two photons is Experimentally, photon pairs with opening angles below R ∼ 0.01 are difficult to distinguish from isolated photons in the calorimeter. Therefore if a is produced in the decay of a TeV scale resonance, one finds that for m a 2 GeV the photon pairs from each a decay can easily pass as isolated photon candidates. 1 Models of this type were considered before in the context of exotic Higgs decays [4,9,10] and more recently as alternative interpretations of the 750 GeV resonance [11][12][13][14][15][16].
Couplings of a light state a to photons are also constrained by low energy data [19][20][21][22][23][24][25]. This makes it impossible to choose m a arbitrarily small. Nevertheless, m a could be so small that its (laboratory frame) decay length becomes comparable to or even larger than the size of the ATLAS and CMS inner detectors (about a meter). If a decays to γ γ at a macroscopic distance from the beam pipe, but still within the inner detector, the two photons have a smaller conversion probability than for quasi-instantaneous a decay. We will consider this possibility in Sect. 4. If the decay length of a is so large that most decays occur outside the electromagnetic calorimeter, they can no longer mimic isolated photons. Such scenarios are, however, still of phenomenological interest in the context of displaced object searches, which look for objects decaying in the calorimeters or in the muon system [26,27].
To be as model independent as possible, we will consider scenarios where a resonance X is produced in proton-proton collisions and decays to two light particles a 1 , a 2 , each of which in turn decays to N i photons: As a concrete realization, consider the case of a scalar resonance S with loop induced couplings to gluons and tree level couplings to a light pseudoscalar a, which in turn couples to photons: Here, m S and m a are the masses of scalar and pseudoscalar, respectively, and 1/ , 1/ f , λ are coupling constants. An LHC process which is induced by these couplings is shown in Fig. 2. The five dimensionful parameters in Eq. (3) are a priori independent and can be extracted from the data. The position 1 To be more precise, the first layer of the EM calorimeter in ATLAS is very finely segmented with η ≈ 0.002 − 0.003, and shower shape variables are used to suppress backgrounds from π 0 decays. So the actual bound on m a could be as low as 500 MeV, as argued in [11].
The main point here is that there is a region of parameter space where collimated photon jets can pass as single isolated photons, so the precise value of the limit is not important.

Fig. 2
Feynman diagram illustrating the production of a scalar resonance S followed by the decay into collimated photon jets of the peak in the photon invariant mass peak determines M S , and the signal cross section together with the decay width of S determines and λ. m a and f have to be chosen such that the photon jets pass as regular photons, which is non-trivial since the coupling f of a light pseudoscalar to photons is strongly constrained [13].

Distinguishing photon jets from isolated photons
Consider a photon jet consisting of N collimated photons. A regular isolated photon corresponds to N = 1 in this notation.
A photon jet will be registered as a converted photon if at least one of the photons inside the jet converts and leaves a signal in the tracker. For a given conversion rate p conv for individual photons in a given jet, the probability that the photon jet appears as a converted photon is then given by Obviously, the probability that the photon jet appears as an unconverted event is For the moment, we neglect the possible issue arising from having more than two reconstructed tracks associated with the photon candidate, which could make the photon fail isolation criteria. We will come back to this point later. Now consider a diphoton event 2 with angular separation R > 0.4 so that they do not overlap. Microscopically, the event contains two photon jets, with the number of photons in them denoted by N 1 and N 2 . From the experimental point of view, we distinguish three event categories, namely events with i = 0, 1, 2 of the photon jets being reconstructed as converted photons. The probabilities p for an event to fall into each of these categories depend on N 1 and N 2 , and they thus offer a handle for distinguishing different theoretical 2 More precisely, an event with two reconstructed photon candidates which are well separated from each other. models underlying a diphoton signal. It is easy to see that In the following, we will in particular consider the prospects for distinguishing a real diphoton resonance, (N 1 N 2 ) = (11), from models with (N 1 N 2 ) = (12) or (22), which have been proposed in the literature as alternative explanations of the 750 GeV signal [11][12][13][14][15][16]. This discrimination is complicated by the fact that there is a significant number of SM background events in the signal region. Here we assume that all background events are of (11) type but we expect the results to remain similar for any other known background composition (see Appendix A for details). Perhaps the simplest statistical way of approximately quantifying the model discrimination power is a Pearson χ 2 test based on the following χ 2 function: .
Here, S and B are the numbers of signal and background events, respectively, (N true 1 N true 2 ) corresponds to the model we assume to be realized in nature, while (N test 1 N test 2 ) describes the model we wish to test against. In other words, the question we are asking here is how likely it is that the hypothesis (N test 1 N test 2 ) is accepted if the actual events are of type (N true 1 N true 2 ). Obviously, the right hand side of Eq.
. The two jets have different p T and pseudorapidities η, which lead to unequal p conv for the photons in different jets. To account for this, we take the p T -and η-dependent conversion probabilities p conv ( p T , η) given in Ref. [17]. As the p T -dependence of p conv ( p T , η) is weak for photons above 100 GeV, we neglect it in the following and work with p conv (η) depending only on the pseudorapidity. The value of p conv (η) in each η bin is listed in Table 1 in Appendix B. The χ 2 function in Eq. (9) is generalized to also include a sum over the events in different bins ( jk) labeled by the pseudorapidities (η j , η k ) of the two jets. The probabilities p for i conversions in an event in rapidity bin ( jk) are given by Eqs. (6)-(8) using the appropriate p conv N 1 (η 1 ) and p conv N 2 (η 2 ) for each jet. Additionally, both terms in the numerator of Eq. (9) as well as the first term in the denominator must now be multiplied by p jk S , the probability in the respective true/test model for signal events to fall into that rapidity bin. Similarly the second term in the denominator must now be multiplied by the analogous probability p jk B for background events. These probabilities p jk S and p jk B can be obtained by computing the differential cross sections for the signal and background (N 1 N 2 ). We do so using MadGraph 5 v2.3.3 [28,29], with a FeynRules / UFO [30] implementation of a simple (N 1 N 2 ) = (11) model that augments the Standard Model with a scalar S and the effective couplings Note that binning the data in pseudorapidity η introduces some model dependence since the differential rapidity distribution will be different from model to model. We assume in the following that the η distribution of the photon jets in models with N 1 , N 2 > 1 is identical to the η distribution of the isolated photons following from Eq. (10). In our numerical results, we will go somewhat beyond the χ 2 test based on Eq. (9), and instead employ a slightly more sensitive likelihood ratio test, as discussed in Appendix B.
For a given S/B, we can now ask how many expected events S + B in the signal region are needed to reject different hypotheses (N 1 N 2 ) at the 2σ and 5σ level. The results are shown in Fig. 3. For S/B of order one, we see that at most a hundred events are necessary to distinguish the different hypotheses at the 2σ level. Discrimination between models of type (11) and (22) model requires fewer events than discrimination between the (21) and (11) or between the (21) and (22) scenarios. The reason is simply that the conversion probabilities Eqs. (6)-(8) for the two alternative hypotheses are more different in the former case. For the particular case of the excess observed around 750 GeV, the present data could already be sufficient to discriminated between the (11) and (22) hypotheses at the 2σ level, while more data would be needed to tell the (21) hypothesis apart from either (11) or (22) scenarios.
We see that photon conversion rates are a promising tool to distinguish between different new physics models in diphoton events once a signal is observed. However, also without an observed event excess, the different conversion probabilities for isolated photons and N > 1 photon jets can be employed as an additional tool to discriminate photon jet signals from the background. In the following we illustrate this, again using the example of a search for a diphoton resonance in the mass range between 200 and 1500 GeV.
In Fig. 4, we show the expected and observed limits on such resonances in the ATLAS diphoton data with 3.2 fb −1 of 13 TeV data [1], and the expected future limits in 300 fb −1 of data. Note that the observed limits shown in Fig. 4 are based on the published ATLAS data, assuming that the η distributions and the conversion rates (which are not public) follow the predictions from simulations. Comparing the limits on (11) resonances to those on (22) resonances, we observe a mild improvement in the latter case.

Long-lived intermediate states
So far, we have assumed that the photon jets in a model with N 1 > 1 or N 2 > 1 form instantaneously at the primary interaction vertex. We consider now a more general scenario, where the intermediate particle a has a non-negligible proper lifetime τ . In this scenario, a decays to photons only after traveling some macroscopic distance x in the inner detector. Since photon conversion cannot take place until the photons have been produced, the conversion probability for an individual photon is reduced. The reduction factor depends on  The expected (dashed) and observed (solid) 95 % CL s limits on a true diphoton signal (the (11) topology in our notation, black curves and Brazilian bands) and on a signal with two photon jets, each consisting of two photons (the (22) topology, red curves). To derive the observed limits, we have assumed that the η distributions of the data and the conversion ratios follow the predictions from Monte Carlo simulations. We show results for an integrated luminosity of 3.2 fb −1 , corresponding to the data published in [1], and for an integrated luminosity of 300 fb −1 . The Brazilian bands were obtained using the CL s method [31,32] as implemented in ROOT. The blue contours show the discrimination power between the (11) and (22) scenarios, defined here as the confidence level at which the (11) hypothesis can be rejected if the signal in the data consists of (22) photon jets many parameters, in particular on the distribution of material in the inner detector and on the efficiency for reconstructing tracks starting away from the beam axis. A full detector simulation is needed to determine this but a key ingredient is a knowledge of the radial dependence of the conversion probability p conv . In the following we outline two simplified approaches.
To obtain an intuitive understanding, it is useful to consider the highly simplistic assumption that the detector is homogeneous. The conversion probability then scales as 1 − x/L t (η), where L t (η) is the total distance from the primary vertex to the edge of the tracker. The probability that at least one photon in an N -photon jet converts to an e + e − pair inside the tracker is the probability that a decays between x − dx to x, and at least one of the N photons converts between x and L t (η), integrated over all x from 0 to L t (η). This is easy to compute and we find where γ is the Lorentz boost of a and p conv (η) on the right hand side is, as in Sect. 3, the probability for a photon to convert between the point of production at the origin and the edge of the tracker at a distance L t (η). We have in particular, for N = 1, 2: Analogously, the probability for an N -photon jet to be detected without any of the photons converting is Here, the first two terms give the probability that the photon jet is detected at all, i.e. that a decays before reaching the calorimeter at a distance L c (η) from the primary vertex. Note that, because of this factor, p conv N (η, τ ) + p no-conv N (η, τ ) < 1. For obtaining our numerical results, we model the conversion probability density as a function of the radial distance traveled using an approximate "two-zone" model based on Ref. [33]. In the central region (|η| < 0.6) there are 70 % conversions in the range (0 < r < 15 cm) and 30 % in (15 cm < r < 40 cm). In the forward region (1.3 < |η| < 1.7) there are 65 % conversions in (0 < r < 15 cm) and 35 % in (15 cm < r < 40 cm), where r is the radial distance. Thus, 1 − x/L t (η) is replaced by the above. The total conversion probability remains the same as before.
In an event with two photon jets, the boost factors γ 1 , γ 2 for the two jets are in general different. Therefore, in the following numerical analysis, we fold the conversion probabilities with the distribution of γ 1 , γ 2 in each (η 1 , η 2 ) bin, obtained from the same MadGraph simulation that determines the (η 1 , η 2 ) distribution (see Sect. 3). Afterwards, the analysis proceeds in the same way as in Sect. 3. In particular, the probability for zero, one or two of the photon jets in an event to convert are given by Eqs. (6)- (8), with the probabilities p conv N and p no-conv N on the right hand side replaced by the two-zone analog of Eqs. 11 and 14. The statistical analysis follows again the procedure described in Appendix B.
In Fig. 5, we show the number of expected signal events S + B required to discriminate between models of (11) and (22) type as a function of the lifetime of the intermediate particle a in the (22) model. We take the mass of the heavy resonance decaying to photon jets to be 750 GeV, the mass of a to be 1 GeV, and we assume the model predicts a signalto-background ratio S/B = 1. We emphasize that the vertical axis in Fig. 5 shows the expected number of detected For definiteness, we assume a signal-to-background ratio S/B = 1, and an a mass of 1 GeV. Note that (S + B) detected only counts those events for which a decays before reaching the EM calorimeter. For γ cτ larger than the inner e-cal radius, this is only a small fraction of the total number of required events, as indicated by the contours of constant total fiducial signal rate S fid diphoton events. Since for non-negligible τ , only those events where both a particles decay before entering the calorimeter are detected, we also show for comparison the total number of signal events in the (22) case (red contours in Fig. 5).
As is to be expected, the discrimination power is best when γ τ = 0, and worsens for longer lifetimes because a decay away from the beam axis (but still well within the tracker) leads to a decreased conversion probability in the (22) model. When γ τ L c , it is likely that a does not decay before reaching the calorimeter, so that events are no longer categorized as diphoton events. However, among those events which are detected, the fraction of converted events increases again. Since the vertical axis in Fig. 5 shows only the number (S + B) detected of detected diphoton events, the discrimination power based on (S + B) detected thus appears to improve again in this case. Note, however, that the condition S/B = 1 requires a significantly larger cross section σ fid BR γ γ when γ τ is large.

Dark photons
An interesting class of models that could in principle mimic a diphoton resonance signal are those where a new heavy particle S decays to two dark photons -the gauge bosons of a new U (1) gauge symmetry, often hypothesized in the context of dark matter models [22]. If the dark photon A is sufficiently light (< m μ /2), its dominant decay mode is A → e + e − , so that the process pp → S → (A → e + e − )+ (A → e + e − ) has the same final state as pp → S → γ γ , with both photons converting to e + e − pairs. 3 While prompt e + e − pairs will be vetoed in the photon reconstruction, if the dark photon lifetime is such that it mostly decays in the tracker, then these events could easily appear as a diphoton resonance.
At close inspections, the two topologies are of course different: first and foremost, A decays can only mimic converted photons. With good statistics, it should therefore be easy to tell an A A signal apart from true isolated diphoton signal. This is indeed the case, as illustrated in Fig. 6: at a signal-to-background ratio of O(1), even a handful of events is enough to discriminate between the A model, denoted here as (EE), and a (11) type diphoton model. Even for a 5σ test, 100 events are needed. However, the difference in the apparent photon conversion probabilities is not the end of the story. For short-lived A decaying quasi-instantaneously, the e + e − tracks in the (EE) model will come directly from the primary vertex, while the tracks from a converted γ can originate at any radius r inside the tracker. Therefore, taking the radial distribution of the secondary vertices into account, the discrimination power can be boosted further. We have refrained from doing so in Fig. 6 to be conservative and because dark photons could also have a macroscopic decay length. In the latter case, the r distribution predicted by the (EE) model is much more similar to that of a converted γ signal in the (11) scenario. Of course, small differences remain. For instance, an A decay can occur anywhere in the tracker, while photon conversion is only possible inside layers of detector material. Moreover, for A decay lengths comparable to the size of the tracker, the r distribution in the (EE) model is exponentially falling, while in the (11) model it is constant.
Note that, without inclusion of the r discrimination of the secondary vertices, the discrimination power depends on the laboratory frame A decay length γ τ only through the factor

Further observables
Photon reconstruction in the LHC detectors offers additional handles that could be used to further discriminate photon jets from isolated photons, and possibly pin down underlying structures like the multiplicity of photons inside the jet. In the following we briefly discuss the most promising ideas: • The photon pairs coming from a boosted decay a → γ γ carry roughly equal energy. If only one of them converts, the ratio of track p T to calorimeter energy E cal should differ substantially from one, the value expected for single photons. This is a very powerful variable that is not currently being used. An accurate measurement of the electron track p T is complicated by their relatively large p T ∼ 100 GeV and the fact that they only traverse part of the tracker, depending on where they convert. Therefore a sufficient number of events is needed such that the measurement can be made on those events where both electrons from the conversion are well reconstructed. • When more than one photon inside a photon jets converts to e + e − , up to 2N tracks could be reconstructed for an N -photon jet. Such multiple conversions might be rejected by the standard photon reconstruction algorithms, for example in ATLAS [34] a cut is placed on the p T sum of tracks within R = 0.3 of the photon candidate which are not associated with the photon candidate itself. Therefore we expect that the reconstruction efficiency for N > 1 photon jets is reduced. However, once a resonance is found, events with larger track mul-tiplicities can be explicitly searched for in loose photon samples to get additional information on the signal. • Photon jet events will react differently from single photons to changes in the isolation criteria. While variables which mainly cut on nearby hadronic activity are insensitive to the photon multiplicity, those using electromagnetic calorimeter shower shapes could be very sensitive.
In the ATLAS search for Higgs decays to pairs of photon jets [35], it was already shown that the variable F side , which considers the ratio of energies deposited in 3 vs. 7 bins centered on the highest bin, is very sensitive to the mass of the intermediate pseudoscalar m a . A large change in efficiency of F side in the signal region compared to the side bands would therefore be a strong indication of photon jets, and even give direct access to m a . In this context, it is worth commenting on photon jets with N > 2 constituent photons. At first glance it seems unlikely that such final states could successfully mimic an isolated photon signal, given how difficult it already is to sufficiently collimate two photons. However, one should also note that for N > 2, the energy does not have to be distributed evenly between the photons. For instance, if one is significantly harder than the others, the energy deposit in the electromagnetic calorimeter would have a single peak structure such that rejection methods based on the shape of the calorimeter cluster would fail.
Each of these strategies can provide additional insight into the nature of a diphoton signal which might be discovered in the future, or more general into any new physics signals involving photons. Comparison with control regions and side bands can be used to verify that abnormal behavior of the photon candidates in the signal region is indeed due to photon jets and not just from e.g. QCD backgrounds.

Conclusions
To summarize, we have discussed from a phenomenologist's point of view how the conversion of photons to e + e − pairs inside the LHC detectors can be exploited to discriminate between final states involving isolated photons and events containing jets of multiple highly collimated photons. Such photon jets arise, for instance, when a light new particle is produced on-shell and decays to two photons. We have illustrated that, even with modest statistics, a resonance decaying to two isolated photons can be distinguished from a new particle decaying to two photon jets. For instance, in the context of the possible 750 GeV resonance observed in ATLAS and CMS data, ∼30 events, are sufficient to make this distinction at the 2σ level, while O(100) events are required for a 5σ discrimination.
We have also illustrated how the sensitivity to photon jet signals mimicking a diphoton resonance depends mildly on the multiplicity of the photon jets. Finally, we have studied scenarios in which photon jets emerge at a macroscopic distance from the beam pipe in the decay of a long-lived intermediate particle.
We conclude that photon candidates in the LHC detectors offer an extremely rich substructure which can be exploited for highly efficient model discrimination. This substructure is theoretically well modeled and seems readily accessible experimentally. We hope that the results presented in this note will be useful in this endeavor.
where L(data|M m ) denotes the likelihood of the data D if the model hypothesis M m is true. A model hypothesis is characterized here by the multiplicities (N 1 N 2 ) of the two photon jets in each signal event, and by the associated prediction P i M m for the event rate in the ith bin. The likelihood is given by In the simplest case, the bin index i = 0, 1, 2 denotes the number of photon jets in the event that are reconstructed as converted photons. However, since the conversion probability depends on the rapidities η 1 , η 2 of the photon jets and since the rapidity dependence is different for signal and background events, we bin the data also in |η 1 |, |η 2 |. This turns i into a multi-index (i jk), with i = 0, 1, 2 the number of converted photon jets, and j, k denoting the rapidity bins. We use 4 of the latter for each photon jet, as given in Table 1.
Note that this binning excludes the transition region between the barrel and the endcap. In the notation of Sect. 3, Table 1 Rapidity bins used in our analysis, together with the associated probabilities p conv for a single photon in a given bin to convert into an e + e − pair bin boundaries in η p conv (η)  While this Monte Carlo-based method for evaluating confidence intervals is very general and, by the Neyman-Pearson lemma, offers optimal discrimination power, it could be replaced by a much simpler χ 2 test. Namely, note that −2 log L(D|M 2 ) follows a χ 2 distribution if the total number of events predicted by M 2 is not too small, The number of degrees of freedom of the χ 2 distribution is given by the number of bins. The expected confidence level at which model M 2 is disfavored if M 1 is true is thus given by the CDF of the χ 2 distribution, evaluated at −2 log L(M 1 |M 2 ). If the number of events in each bin is 10, so that the Poissonian likelihood Eq. (B2) is well approximated by the Gaussian likelihood we recover the χ 2 from Eq. (9).