Search for high-mass diphoton resonances in $pp$ collisions at $\sqrt{s}=8$ TeV with the ATLAS detector

This article describes a search for high-mass resonances decaying to a pair of photons using a sample of $20.3$ fb$^{-1}$ of $pp$ collisions at $\sqrt{s}=8$ TeV recorded with the ATLAS detector at the Large Hadron Collider. The data are found to be in agreement with the Standard Model prediction, and limits are reported in the framework of the Randall-Sundrum model. This theory leads to the prediction of graviton states, the lightest of which could be observed at the Large Hadron Collider. A lower limit of $2.66$ ($1.41$) TeV at 95% confidence level is set on the mass of the lightest graviton for couplings of $k/\overline{M}_{\mathrm{Pl}} = 0.1$ ($0.01$).


Introduction
One of the primary goals of the experiments at the Large Hadron Collider (LHC) is the search for new phenomena which might manifest themselves in the high-energy regime made accessible for the first time by this machine.This article reports on a search for new high-mass resonances in the diphoton channel using the full data set of proton-proton (pp) collisions recorded with the ATLAS detector [1] at √ s = 8 TeV in 2012.
The diphoton channel is important for searches for physics beyond the Standard Model (SM) since new high-mass states decaying to two photons are predicted in many extensions of the SM, and since searches in the diphoton channel benefit from a clean experimental signature: excellent mass resolution and modest backgrounds.One popular example of theories predicting a new high-mass diphoton resonance is the Randall-Sundrum (RS) model [2].This model attempts to solve the so-called hierarchy problem of the SM, i.e. the fact that the electroweak scale is 16 orders of magnitude smaller than the Planck scale (M Pl ≈ 10 16 TeV).Despite significant efforts over the past few decades, this problem has not been solved so far.In the RS model, the observed hierarchy of scales is accommodated assuming an extra spatial dimension.In this paradigm, the fundamental strength of gravity is comparable to the strength of the electroweak interaction, but gravity appears much weaker than the other interactions because it is diluted by the presence of the extra spatial dimension.
The RS model postulates a "warped" five-dimensional spacetime configuration.In this configuration there are two 3-branes that are the boundaries of the five-dimensional spacetime.The 3-branes can support (3 + 1)-dimensional field theories.The SM fields are located on one of the two 3-branes, the so-called TeV brane, while gravity originates from the other brane, called the Planck brane.Gravitons propagate in the five-dimensional bulk.The fifth dimension is compactified with length r c .The solution to Einstein's equations for this configuration is given by a spacetime metric that couples the standard four dimensions to the extra dimension with an exponential "warp factor".In this configuration, a field with the fundamental mass parameter m 0 will appear to have the physical mass m = e −kπr c m 0 , where k is the curvature scale of the extra dimension.Scales of the order of 1 TeV are therefore generated from fundamental scales of the order of M Pl if kr c ≈ 12 [3].The compactification of the extra dimension gives rise to a Kaluza-Klein (KK) [4] tower of graviton excitations G * -a set of particles on the TeV brane with increasing mass.The phenomenology can be described in terms of the mass of the lightest KK graviton excitation (m G * ) and the dimensionless coupling to the SM fields, k/M Pl , where M Pl = M Pl / √ 8π is the reduced Planck scale.The masses and couplings of the individual graviton states are determined by the scale Λ π = e −kπr c M Pl ≈ O(TeV) [3].Values of k/M Pl in the range [0.01, 0.1] are preferred from theoretical arguments [3].The lightest graviton state is expected to be a fairly narrow resonance for k/M Pl < 0.3 [3].Its natural width varies like the square of k/M Pl .For k/M Pl = 0.1, the natural width increases from ∼ 8 GeV at m G * = 800 GeV to ∼ 30 GeV at m G * = 2200 GeV.In the first example, the natural width and the experimental mass resolution are similar, in the second example the natural width dominates the total width1 .
Searches for high-mass diphoton resonances at the LHC have been reported by both the ATLAS and CMS collaborations, and the latest results are summarized in Table 1.The most stringent experimental limits on RS gravitons to date have been obtained at the LHC, using a combination of the dielectron and dimuon channels.Limits from searches for narrow resonances in other channels are available in, for example, the dijet and WW channels, but these channels are not competitive for the specific case of RS gravitons.The limits quoted here are also included in Table 1, along with those from earlier searches at the Tevatron.2 The ATLAS detector The ATLAS detector is described in detail in Ref. [1].It consists of an inner detector (ID) surrounded by a solenoid that produces a 2 T magnetic field, electromagnetic (EM) and hadronic calorimeters, and a muon spectrometer which employs air-core toroidal magnets.In the pseudorapidity2 region |η| < 2.5, the ID provides measurements of the tracks originating from the primary pp collision as well as tracks from secondary vertices, permitting the efficient identification of photon conversions.The EM calorimeter (ECAL) is a lead/liquid-argon (LAr) sampling calorimeter.It is divided into a barrel section covering the region |η| < 1.475 and two endcap sections covering 1.375 < |η| < 3.2.The ECAL is segmented longitudinally in shower depth into three layers for |η| < 2.5 and two for 2.5 < |η| < 3.2.Up to |η| < 2.4, the first layer consists of highly granular strips segmented in the η direction for efficient candidate-bycandidate discrimination between single photon showers and two overlapping showers originating from π 0 decay.The second layer, with a typical granularity of ∆η × ∆φ = 0.025 × 0.025, collects most of the energy from photon showers.Significant energy deposited in the third layer is an indication for leakage beyond the ECAL from a high-energy shower, and the measurements from the third layer are used to correct for this effect.A thin presampler layer, installed in front of the ECAL and covering the pseudorapidity interval |η| < 1.8, is used to correct for energy loss before the ECAL.In the region |η| < 1.7, the hadronic sampling calorimeter uses iron absorbers and plastic scintillator tiles as the active material.In the region 1.5 < |η| < 4.9, LAr is used as the active material, with copper or/and tungsten absorbers.The muon spectrometer features three stations of precision chambers that allow for accurate measurements of the muon track curvature in the region |η| < 2.7.Fast detectors enable muon triggering in the range |η| < 2.4.A three-level trigger system is used to select the events to be recorded for subsequent offline analysis.

Event reconstruction and selection
The events used in this analysis were recorded using a diphoton trigger.At the highest, software-based level of the trigger system, this trigger requires at least one photon with transverse energy E γ T > 35 GeV and at least one other photon with E γ T > 25 GeV.Both photons need to satisfy requirements on the shape of the energy deposit in the calorimeter.This trigger is nearly fully efficient for diphoton pairs that pass the offline event selection.In the offline analysis, only events that pass the standard ATLAS data quality requirements are considered.Photons are reconstructed from fixed-size clusters of cells in the EM calorimeter.To ensure good diphoton mass resolution, photons are required to be located inside the precision region of the EM calorimeter (|η| < 1.37 or 1.52 < |η| < 2.37).They are further required to satisfy E γ T > 50 GeV as well as the loose identification requirements from Ref. [13] updated for 2012 run conditions.The two highest-E γ T photon candidates that pass these requirements are retained for further analysis.Events are retained for the main analysis if both candidates also satisfy the tight identification requirements [13] updated for the 2012 run conditions.These requirements include upper limits on the energy leakage into the hadronic calorimeter and on the width of the shower in the first and second layers of the EM calorimeter.Events where at least one candidate fails the tight photon identification are retained for the study of backgrounds with misidentified (fake) photons in data control samples.The invariant mass of the diphoton system is evaluated using the photon energies measured in the calorimeter as well as the direction of the photon momenta determined from the positions of the photons in the calorimeter and the position of the primary vertex of the hard interaction.Primary vertices are reconstructed from the tracks in the ID and are required to be formed from at least three tracks.The primary vertex of the hard interaction is identified using information on the directions of flight of the photons as determined using the longitudinal segmentation of the ECAL (calorimeter pointing), the parameters of the beam spot and the properties of the tracks associated with each reconstructed vertex [14].This procedure ensures good diphoton mass resolution despite the presence of pileup.Pileup refers to the additional pp interactions that occur in the same bunch crossing as the pp → γγ + X interaction or in the adjacent crossings.A calorimetric isolation observable is defined for photon candidates.This observable provides further discrimination between genuine photons and fake candidates from jets ( j).It is used for two purposes in this analysis: to quantify the composition of the data sample in terms of true and fake photon candidates, and to suppress the reducible background by requiring both photons to be isolated.The isolation energy E iso T is defined as the sum of the transverse energies of cells inside clusters of variable size (also called "topo-clusters" [15]) inside a cone of ∆R = (∆η) 2 + (∆φ) 2 = 0.4 around the photon candidate, excluding a 5 × 7 grid of EM cells in the center of the cone.The measurement of E iso T is corrected for the following two effects: i) soft energy contributions from the underlying event and pileup that fall inside the isolation cone and ii) leakage of photon energy outside the central core, causing the isolation energy to grow as a function of photon E γ T .The expected soft contributions are subtracted from E iso T on an event-by-event basis; the correction is based on a measurement of the properties of low-p T jets [16] in the event.The average expected contributions to E iso T from leakage are subtracted based on simulations.The fluctuations (from one photon shower to another) around the average expected leakage contribution to E iso T become more important as E γ T increases.In particular, the tail at large E iso T becomes more pronounced and would cause significant inefficiencies at large E γ T if a fixed upper limit on E iso T were imposed in the event selection.To avoid inefficiencies, a negative E γ T -dependent offset is added to E iso T , and the requirement E iso T < 8 GeV is imposed on both photon candidates after the addition of the offset.In the following, the symbol E iso T refers to the isolation energy with the offset applied.The E γ T -dependence of the offset is chosen such that the efficiency of the E iso T requirement is independent of E γ T in simulated samples, and the size of the offset varies from zero at low E γ T to −5 GeV at E γ T = 1 TeV.The corresponding efficiency of the isolation requirement is 96% per photon.For the purpose of the determination of the background composition, a less stringent requirement of E iso T < 14 GeV is used.Examples of the E iso T distributions will be shown later, after the discussion of the backgrounds.

Background estimate
An important contribution to the background arises from prompt γγ production via Standard Model processes.This contribution is irreducible and, as quantified below, represents the dominant source of background.Another significant contribution is due to events in which one or both of the photon candidates arise from other objects, such as misidentified jets or electrons.This background component is dominated by γ + j and j + j events with one or two jets reconstructed as photons.Backgrounds with electrons misidentified as photons -e.g.electron-positron events from Drell-Yan production, in W/Z + γ or in tt events -were verified to be negligible after the event selection.The irreducible background is estimated using Monte Carlo (MC) simulations normalized to the data in a control region defined below, and the reducible background is estimated using control samples in data.A combination of the leading-order (LO) event generator Pythia [17] and the next-to-leading-order (NLO) parton-level simulation Diphox [18] is used to simulate the irreducible background.Monte Carlo samples are generated using Pythia 8.163 with the CTEQ6L1 [19] parton distribution functions (PDFs) and the AU2 Pythia parameter tune [20].The detector response is simulated using Geant 4 [21,22], including pileup conditions similar to those observed in data.The simulated events are reweighted such that the number of pp interactions per bunch crossing has the same distribution as in data.The NLO correction and the contributions of the fragmentation process are evaluated using Diphox 1.3.2 with the CTEQ6.6M[23] PDF set and lead to a correction of the Pythia prediction that depends on the generated mass of the diphoton system.The so-called box contribution gg → γγ through a quark loop is included in both the Pythia and Diphox predictions.From the point of view of power counting, this diagram is a next-to-next-to-leading-order (NNLO) contribution, but given the large gluon luminosity at the LHC compared to the quark-antiquark one, the contribution of the box diagram is comparable to that of the LO qq → γγ process.
To reduce the impact of the theory uncertainties in the absolute normalization of the production crosssection of the irreducible background, the background estimates are scaled to the data at low mass.In the mass region m γγ < 600 GeV, where m γγ denotes the measured invariant mass of the diphoton system, the most sensitive searches for narrow resonances in the diphoton channel rely on background estimates that make use of low-mass and high-mass sidebands [24].The search described in the present article focuses on higher diphoton masses where the sidebands, especially at high mass, provide insufficient event yields.This search therefore uses a different approach to obtain the background prediction at high m γγ .Specifically, the simulated irreducible background is normalized to the data in a low-mass control region.This control region is defined to coincide with the first 22 bins in Fig. 3, i.e. as 179 < m γγ < 409 GeV.To determine the composition of the data sample in the low-mass control region in terms of irreducible and reducible backgrounds, a template fit to the E iso T distributions of the leading and subleading photon candidates is performed.A similar technique was used in earlier analyses of low-mass diphoton events [25,26].The E iso T distributions, along with the result of the template fit, are shown in Fig. 1 (top and middle).The full event selection described in the text, with the relaxed requirement E iso T < 14 GeV, has been applied.The fit is performed simultaneously to the two distributions, and the normalizations for four background components are allowed to float in the fit, namely the irreducible γγ component plus three reducible components: γ + j (subleading photon candidate is fake), j + γ and j + j.Templates for the E iso T distribution of true photons and of fake photons are determined from data.The template for fake photons is obtained using a sample of photon candidates that satisfy the loose identification criteria and fail to meet a subset of the criteria in the tight selection.This subset has been chosen to include criteria that have been found to be only weakly correlated with E iso T .The template for true photons is obtained from the sample of tight photon candidates after subtraction of the fake-photon template normalized to the tight data sample at E iso T > 10 GeV where fake photons dominate.The fit is performed separately in five subsamples of the data, because the relative contributions of the four background components vary significantly from one subsample to another.The five subsamples correspond to different combinations of pseudorapidities of the leading and subleading photon candidates: both photons in the "central" region (|η| < 1.37), leading photon central and subleading photon in the "endcap" region (1.52 < |η| < 2.37), leading photon in endcap and subleading central, both photons in same endcap and both photons in endcaps of opposite-sign pseudorapidity.In Fig. 1, the data from the five subsamples, as well as the corresponding fit results are combined, and only the sum of the photon and jet components over the four background components is shown.To illustrate the broadening of the E iso T distribution at large E γ T described above, Fig. 1 also includes the E iso T distribution for photons with E γ T > 400 GeV.The shapes of the distribution of the isolation energy for genuine photons extracted from data (in bins of E γ T ) are compared to the shapes predicted using simulation.The backgrounds in the search region, i.e. at masses higher than the low-mass control region, are predicted using estimates of the m γγ shape of each background component and the normalizations obtained in the low-mass control region described above.As discussed above, the shape of the irreducible component is obtained using simulation.The shapes of the reducible components are obtained from data control samples.The control sample used to extract the shape of the γ + j background is selected in the same way as the signal sample, except that the subleading photon candidate has to satisfy the loose identification criteria and fail to meet the tight criteria.The control samples for the j + γ and j + j backgrounds are defined using the same approach.Since these data control samples contain relatively few events at high m γγ , a fit to a smooth function is used to extrapolate the shapes of the reducible backgrounds to high m γγ .Specifically, the functional form from Ref. [27] is used.It describes well the shapes of the data in the control samples of this analysis.This functional form was also successful in describing the shapes of the background control samples in earlier searches in the dijet [27,28], dilepton [29], photon-jet [30] and diphoton [5] channels.The impact of the small observed differences in the isolation distributions on the predicted m γγ shape of the irreducible background and on the signal efficiency is quantified and taken into account in the systematic uncertainties.
The uncertainties in the prediction of the total background are shown in Fig. 2 as a function of m γγ .The uncertainties in the shape of the m γγ distribution of the irreducible background are dominated by the uncertainties in the PDF set used for the Diphox simulation.The set of PDF variations provided with the CTEQ6.6MPDF set are used to propagate the PDF uncertainty to the shape of the irreducible background.The comparatively small difference between the irreducible background shapes obtained with the CTEQ6.6Mand the MSTW2008 [31] PDF sets is added in quadrature to the uncertainty.An additional uncertainty arises from higher-order contributions that are not included in the Diphox generator.This uncertainty is evaluated by varying in a correlated and anticorrelated way the renormalization scale and the factorization scales in the hard scattering process and in the photon fragmentation function [18] by a factor two around their nominal value, which is set to the mass of the diphoton system.A small uncer- tainty arises from residual imperfections in the simulation of the isolation energy E iso T .The uncertainties in the shape of the reducible background arise from the finite size of the corresponding control samples, and from the extrapolation of the background shapes from the control samples to the signal region.The latter are assessed by varying the definition of the loose selection requirement that is used to define the control samples.The uncertainties in the normalizations of the individual background components are dominated by the uncertainties in the template shapes that are used in the fit to the isolation energy.They are assessed by varying the definition of the non-tight selection requirements that are used to obtain the template shapes.

Signal response
Simulated event samples of RS graviton signals at different values of m G * are used to study the response of the detector for signal events.The samples are generated using Pythia 8.163 with the CTEQ6L1 PDF set, the AU2 tune and the same detector and pileup simulation as for the SM γγ samples described above.The sources of uncertainty in the expected signal yield are summarized in Table 2.The systematic uncertainty is dominated by the uncertainty in the efficiency due to residual imperfections in the simulation of the observables that are used in the tight photon identification criteria and in the simulation of the isolation energy.Smaller contributions arise due to the limited size of the simulated samples and the determination of the trigger efficiency.The uncertainty in the integrated luminosity is also propagated to the expected yield.Uncertainties in the RS resonance shape due to our current knowledge of the photon energy scale and resolution as well as imperfections in the simulation of the pileup conditions were verified to have a negligible effect on the result.

Results
The observed diphoton mass spectrum is shown in Fig. 3, together with the background expectation and the predicted signal for two examples of values of the RS model parameters.The acceptance (i.e. the fraction of simulated events in which both photons satisfy the requirements on |η| and E γ T described above) for the RS scenario with m G * = 2.0 TeV and k/M Pl = 0.1 is 89.3% and the selection efficiency for Table 3: Number of events expected from the reducible and irreducible background components as well as the total background and observed number of events in different mass bins.Each bin boundary in the table corresponds to a bin boundary in Fig. 3, but the binning of the table is more coarse than in the figure.Since the combined statistical and systematic uncertainties are strongly anticorrelated between the two components, the uncertainty in the total background prediction is smaller than the sum in quadrature of the uncertainties of the individual components.

Mass window
Background expectation (number of events) Observed 2.12 ± 0.61 0.41 ± 0.22 2.52 ± 0.30 3 events within acceptance is (55.9 ± 2.6)%.The bottom insert in Fig. 3 shows, bin-by-bin, the statistical significance of the difference between the data and the expected background.This significance is calculated using the prescription in Ref. [32] and is plotted as positive (negative) for bins with an excess (deficit) of data.The observed event yields and expected backgrounds are also summarized in Table 3.
No significant deviation of the data from the expected background is observed.To further quantify the level of agreement, the BumpHunter [33] algorithm is used to search for any "bump" on top of the expected background, i.e. any excess of events localized in a relatively narrow mass region.Specifically, the BumpHunter algorithm scans the mass distribution in windows of progressively increasing width in order to identify the window with the most significant excess of observed events over the background anywhere in a predefined search region.In this analysis, the binned m γγ spectrum from Fig. 3 is used, the search region is defined as 409 < m γγ < 3000 GeV, and the window size is allowed to vary from one bin to half the number of bins in the search region.The most significant excess is found in the region 494 < m γγ < 745 GeV.The probability to observe anywhere in the search region an excess that is at least as significant as the one observed in data and that arises from fluctuations in the background is 58%.This confirms the absence of any significant excess in the data.
As no evidence for a signal is found, upper limits at 95% confidence level (CL) are set on the production cross section times branching ratio, σ × BR(G * → γγ), for the lightest RS graviton.A Bayesian approach is used to calculate the limits, and the corresponding computer program is implemented using the toolkit from Refs.[34,35].A likelihood function is used to describe the probability to obtain the m γγ distribution observed in data, given the background and signal predictions.The likelihood is defined as the product of the per-bin probabilities over all bins in the search region (m γγ > 409 GeV) in Fig. 3.In each bin, the probability is computed from a Poisson distribution given the observed number of events and the model prediction, where the model prediction is the sum of the expected numbers of background and signal events.The expected signal yield depends on m G * and σ × BR(G * → γγ).Using Bayes' theorem, the likelihood is converted into a posterior probability.A uniform positive prior in σ × BR(G * → γγ) is assumed, and systematic uncertainties are incorporated using nuisance parameters with Gaussian priors.For a given value of m G * , the posterior probability is reduced to one parameter of interest, σ × BR(G * → γγ), and the 95% CL limit on σ × BR(G * → γγ) is obtained via integration of the posterior probability.
The resulting limits on σ × BR(G * → γγ), as a function of the mass hypothesis (value of m G * ), are shown in Fig. 4. The phenomenology of the RS model can be described in terms of the parameters m G * and k/M Pl , i.e. σ × BR(G * → γγ) can be calculated [36,37] for a given pair of values for these parameters.The results of this calculation are represented by the lines for given values of k/M Pl in Fig. 4.They are obtained using the Pythia generator, which implements the calculations from Ref. [37].A K-factor is applied to these predictions to account for QCD corrections at NLO.The value of the K-factor varies between 1.6 and 1.9, depending on m G * and k/M Pl .These results for the K-factor were provided by the authors of Refs.[38,39], using updated calculations at √ s = 8 TeV.Given these theory predictions for σ × BR(G * → γγ), the upper limits on σ × BR(G * → γγ) from Fig. 4 can be translated into upper limits on k/M Pl .The results are shown, as a function of m G * , in Fig. 5.Alternatively, the upper limits on σ × BR(G * → γγ) can be translated into lower limits on the graviton mass.The results are shown in Table 4, for different values of k/M Pl .The limits presented in Fig. 5 and in Table 4 do not take into account any theory uncertainties in the prediction of σ × BR(G * → γγ) for the RS signal.The size of one of the dominant theory uncertainties (the PDF uncertainty) is included in Fig. 4; uncertainties due to missing higher-order contributions in the calculation of the K-factor are estimated to be of comparable size.

Conclusions
In summary, a search for high-mass diphoton resonances has been performed using the full 2012 dataset collected by the ATLAS detector at the LHC (20.3 fb −1 of pp collisions at √ s = 8 TeV).No significant excess over the expected background is observed, and upper limits on the production cross section times branching fraction σ × BR(G * → γγ) for narrow resonances are reported as a function of the resonance mass.High-mass diphoton resonances are predicted e.g. in models that postulate the existence of extra spatial dimensions in order to address the hierarchy problem, and could be observed at the LHC.Limits on the mass of the lightest graviton state in the framework of the Randall-Sundrum model of extra dimensions are derived from the limits on σ × BR(G * → γγ).The only free parameters of this model are the mass of the lightest graviton (m G * ) and the coupling to the SM fields (k/M Pl ), i.e. σ × BR(G * → γγ) can be calculated once the values of these two parameters are specified.The limits can therefore be expressed in terms of m G * for a given value of k/M Pl .A lower limit of 2.66 (1.41) TeV at 95% confidence level is obtained on the mass of the lightest graviton for a coupling k/M Pl = 0.1 (0.01).The results reported in this article constitute a significant improvement over the results in Ref. [5] obtained at √ s = 7 TeV.The upper limits on σ × BR(G * → γγ) are reduced by a factor of 2.1 (3.6) for a resonance with a mass of 0.5 TeV (2.75 TeV).The improvement in terms of limits on the parameters of any model is even larger, as the new limits are expressed in terms of the production cross section at √ s = 8 TeV, compared to √ s = 7 TeV for the old limits.The lower limits on the mass of the lightest RS graviton are also significantly improved compared to the previous analysis.d Also at TRIUMF, Vancouver BC, Canada e Also at Department of Physics, California State University, Fresno CA, United States of America f Also at Department of Physics, University of Fribourg, Fribourg, Switzerland g Also at Departamento de Fisica e Astronomia, Faculdade de Ciencias, Universidade do Porto, Portugal h Also at Tomsk State University, Tomsk, Russia i Also at CPPM, Aix-Marseille Université and CNRS/IN2P3, Marseille, France j Also at Università di Napoli Parthenope, Napoli, Italy k Also at Institute of Particle Physics (IPP), Canada l Also at Particle Physics Department, Rutherford Appleton Laboratory, Didcot, United Kingdom m Also at Department of Physics, St. Petersburg State Polytechnical University, St. Petersburg, Russia n Also at Louisiana Tech University, Ruston LA, United States of America o Also at Institucio Catalana de Recerca i Estudis Avancats, ICREA, Barcelona, Spain p Also at Department of Physics, National Tsing Hua University, Taiwan q Also at Department of Physics, The University of Texas at Austin, Austin TX, United States of America r Also at Institute of Theoretical Physics, Ilia State University, Tbilisi, Georgia s Also at CERN, Geneva, Switzerland t Also at Georgian Technical University (GTU),Tbilisi, Georgia u Also at Ochadai Academic Production, Ochanomizu University, Tokyo, Japan v Also at Manhattan College, New York NY, United States of America w Also at Institute of Physics, Academia Sinica, Taipei, Taiwan

Figure 3 :
Figure3: Observed invariant mass distribution of the selected diphoton events (black dots; the vertical and horizontal axes on logarithmic scales).To compensate for the rapid decrease of the spectrum, the bins are chosen to have constant logarithmic width.Specifically, the ratio of the upper to lower bin boundary is equal to 1.038 for all bins, and the first bin starts at 179 GeV.Superimposed are the SM background prediction including irreducible and reducible components and two examples of signal predictions.The low-mass control region is indicated by an arrow.(Bottom insert) Bin-by-bin significance of the difference between data and background expectation.

Figure 4 :
Figure 4:  Expected and observed upper limits on σ × BR(G * → γγ) expressed at 95% CL, as a function of the assumed value of the graviton mass.At large m G * , the −1σ and −2σ variations of the expected limit tend to be particularly close to the expected limit.This is expected, since signals with high m G * would appear in regions of m γγ where the background expectation is small and the Poissonian fluctuations around the mean expected background are highly asymmetric.The curves show the RS model prediction for given values of k/M Pl as a function of m G * .They are obtained using the Pythia generator plus a K-factor to account for NLO corrections (see text).The thickness of the theory curve for k/M Pl = 0.1 illustrates the PDF uncertainties expressed at 90% CL.

Figure 5 :
Figure 5: Expected and observed upper limits on k/M Pl expressed at 95% CL, as a function of the assumed value of the graviton mass.

x
Also at LAL, Université Paris-Sud and CNRS/IN2P3, Orsay, France y Also at Academia Sinica Grid Computing, Institute of Physics, Academia Sinica, Taipei, Taiwan z Also at Moscow Institute of Physics and Technology State University, Dolgoprudny, Russia aa Also at Section de Physique, Université de Genève, Geneva, Switzerland ab Also at International School for Advanced Studies (SISSA), Trieste, Italy ac Also at Department of Physics and Astronomy, University of South Carolina, Columbia SC, United States of America ad Also at School of Physics and Engineering, Sun Yat-sen University, Guangzhou, China ae Also at Faculty of Physics, M.V.Lomonosov Moscow State University, Moscow, Russia a f Also at National Research Nuclear University MEPhI, Moscow, Russia ag Also at Department of Physics, Stanford University, Stanford CA, United States of America ah Also at Institute for Particle and Nuclear Physics, Wigner Research Centre for Physics, Budapest, Hungary ai Also at Department of Physics, Oxford University, Oxford, United Kingdom a j Also at Department of Physics, The University of Michigan, Ann Arbor MI, United States of America ak Also at Discipline of Physics, University of KwaZulu-Natal, Durban, South Africa al Also at University of Malaya, Department of Physics, Kuala Lumpur, Malaysia * Deceased

Table 1 :
Overview of previous limits on the mass of the lightest RS graviton state, expressed at 95% confidence level.The notation + − refers to a combination of the dielectron and dimuon channels.The CDF and D0 experiments operated at the Tevatron pp collider.

Table 2 :
Summary of systematic uncertainties on the expected signal yield, excluding those on the production cross section and the acceptance.The total systematic uncertainty is obtained by summing the individual contributions in quadrature.

Table 4 :
Expected and observed lower limits at 95% CL on m G * for different values of k/M Pl .