{'fraction_non_alphanumeric': 0.03432500547885163, 'fraction_numerical': 0.01594345825115056, 'mean_word_length': 4.552091254752852, 'pattern_counts': {'":': 0, '<': 16, '<?xml version=': 0, '>': 20, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Search for contact interactions and large extra dimensions in dilepton events from pp collisions at sqrt(s) = 7 TeV with the ATLAS detector

A search for non-resonant new phenomena, originating from either contact interactions or large extra spatial dimensions, has been carried out using events with two isolated electrons or muons. These events, produced at the LHC in proton-proton collisions at sqrt(s) = 7 TeV, were recorded by the ATLAS detector. The data sample, collected throughout 2011, corresponds to an integrated luminosity of 4.9 and 5.0 fb-1 in the e+e- and mu+mu- channels, respectively. No significant deviations from the Standard Model expectation are observed. Using a Bayesian approach, 95% confidence limit lower limits ranging from 9.0 to 13.9 TeV are placed on the energy scale of llqq contact interactions in the left-left isoscalar model. Lower limits ranging from 2.4 to 3.9 TeV are also set on the string scale in large extra dimension models. After combination of these limits with results from a similar search in the diphoton channel, slightly more stringent limits are obtained.


I. INTRODUCTION
Extensions to the Standard Model (SM), such as quark/lepton compositeness and large extra dimensions, predict modifications to the SM dilepton invariant mass spectra. This paper presents a comparison of the number of expected and observed events at high mass in the dielectron and dimuon datasets collected by the ATLAS detector [1] in 2011. These events resulted from protonproton collisions produced at √ s = 7 TeV by the LHC [2]. The data are interpreted in the context of contact interactions (CI) and virtual graviton exchange in the Arkani-Hamed-Dimopoulos-Dvali (ADD) model [3].
In the SM, quarks and leptons are fundamental particles. However, if they are composite particles, with at least one common constituent, the interactions of these constituents would manifest themselves through an effective four-fermion contact interaction at energies well below the compositeness scale. This type of contact interaction could also describe a new force with a messenger too heavy for direct observation at the LHC, in analogy with Fermi's nuclear β decay theory [4].
The Lagrangian for a general contact interaction has the form [5]: where g is a coupling constant chosen so that g 2 /4π = 1; Λ is the contact interaction scale, which in the context of compositeness models, is the energy scale below which fermion constituents are bound; and ψ L,R are left-handed and right-handed fermion fields, respectively. The parameters η ij , where i and j are L or R (left or right), define the chiral structure of the new interaction. Specific models are constructed by assigning particular combinations of these parameters to be −1, 0 or +1. For example, the left-left isoscalar model (LLIM) is defined by setting η LL = ±1 and η RR = η LR = 0. The LLIM model, commonly used as a benchmark for contact interaction searches [6], is utilized in this analysis. The addition of the contact interaction Lagrangian to that of the SM modifies the Drell-Yan (DY) production cross section (qq → Z/γ * → ℓ + ℓ − ). The largest deviations in the dilepton invariant mass spectra, either constructive or destructive, are expected at high mass and are determined by the sign of the parameter η ij and the scale Λ. The differential cross section for the process qq → ℓ + ℓ − , including a contact interaction, can be separated into three components: a SM DY term, a pure contact interaction term (F C ) and a DY-CI interference (F I ) term: where m ℓℓ represents the final-state dilepton mass. The full form of this expression is given in Ref. [7]. Constructive (destructive) interference corresponds to η LL = −1 (+1). At the largest Λ values to which this analysis is sensitive, both interference and pure contact interaction terms play significant roles. For example, at dilepton masses greater than 400 GeV and Λ = 12 TeV, the magnitude of the middle term in Eq. (2), which depends on the interference, is about twice that of the last term. Nonresonant deviations in the high mass dilepton invariant mass spectra are also predicted in large extra dimension models. These models were introduced to address some of the major unresolved issues in particle physics such as the hierarchy problem. The latter deals with the question of why gravity appears weak in comparison to the other three SM interactions and why the electroweak scale (∼1 TeV) is 16 orders of magnitude smaller than the Planck scale (M Pl ≃ 10 16 TeV). Arkani-Hamed, Dimopoulos, and Dvali addressed these issues by postulating the existence of n flat additional spatial dimensions of common size R, compactified on an n-dimensional torus [3]. The fundamental Planck scale in (4+n)-dimensional spacetime, M D , is then related to the scale M Pl by Gauss's law: M 2 Pl = M n+2 D R n . Con-sequently, the hierarchy problem can be solved with a Planck scale resulting from a fundamental scale (M n+2 D ) near 1 TeV if the volume (∝ R n ) is large enough.
In the ADD model, the SM particles and their interactions are confined to a three-dimensional slice of the multidimensional world, but gravity permeates the additional dimensions of size R. This results in Kaluza-Klein (KK) modes of the graviton. The mass splitting of these KK modes is determined by the factor 1/R. Resolution of the hierarchy problem necessitates large extra dimensional volumes and consequently implies small values of 1/R. This results in an almost continuous spectrum of KK graviton states and hence a nonresonant increase in the expected rate of dilepton events at large invariant mass. Performing the sum over the KK modes in the virtual graviton exchange process leads to an integral which has to be regulated by an ultraviolet cutoff value (Λ UV ). The ADD model is a low-energy effective theory valid below the scale of the onset of quantum gravity, characterized by the scale M S . The convention used throughout this analysis is to equate the cutoff to the scale of the effective theory (Λ UV =M S ).
For virtual graviton exchange, it is standard practice to present limits on the size of the extra dimensions in terms of M S , taken to be the string scale, which is related to M D by the following expression [8]: The strength of gravity in the presence of extra dimensions is typically parametrized by η G = F /M 4 S , where F is a dimensionless parameter of order unity. The definition of F depends on the formalism chosen [8], with three popular conventions: Giudice-Rattazzi-Wells (GRW) [9], Hewett [10] and Han-Lykken-Zhang (HLZ) [11]. The different values are In the GRW and HLZ representations, gravitational effects interfere constructively with the SM processes, while in Hewett's convention there can be destructive or constructive interference. This is encapsulated in the parameter λ, which is equal to +1 (−1) for constructive (destructive) interference. The total cross section (σ tot ), including effects of qq-and gg-initiated virtual graviton exchange, may be parametrized as where σ SM is the SM cross section for the process being considered, and F int and F G are functions of the cross sections involving the interference and pure graviton effects, respectively. Note that the interference term has a linear dependence on η G , and therefore a quartic dependence on M S (i.e., η G ∝ 1/M 4 S ), whereas the pure graviton exchange term is quadratic in η G and therefore has a 1/M 8 S dependence. A study of signal yields in the kinematic range relevant to this analysis shows that the pure graviton term dominates those yields. This is in part due to the fact that the gg-initiated contribution to the graviton exchange process does not interfere with the qq-initiated DY process. Results are nevertheless presented for both 1/M 4 S and 1/M 8 S priors. Previous searches for contact interactions have been carried out in neutrino-nucleus and electron-electron scattering [12,13], as well as at electron-positron [14][15][16][17][18], electron-proton [19,20], and hadron colliders [21][22][23][24][25][26][27][28]. In the case of eeqq contact interactions, the limits in the LLIM for all quark flavors from e + e − experiments are Λ − > 7.2 TeV and Λ + > 12.9 TeV [14] at 95% confidence level (C.L.) for η LL = −1 and +1, respectively. These limits assume that contact interactions of electrons with all quark flavors are of the same strength. The best limits set in the specific case of first generation quarks are Λ − > 9.1 TeV and Λ + > 8.6 TeV [18] at 95% C.L. In the case of eeqq contact interactions, the best limit for constructive interference is Λ − > 10.1 TeV from the ATLAS analysis of the first 1 fb −1 of 2011 data [28]. The best limits in the case of µµqq contact interactions are from an analysis of the same data: Λ − > 8.0 TeV and Λ + > 7.0 TeV [28].
Previous searches for large extra dimensions in the ADD model via virtual graviton exchange have been performed at electron-positron [29][30][31][32][33][34], electron-proton [20,35], and hadron colliders [25,[36][37][38][39][40][41][42]. Presently, the most stringent mass limits in the dielectron and dimuon channels require M S > 2.8 TeV for each channel and 3.1 TeV when combined (in the GRW formalism with no K factor) [38]. The best limits to date arise from the combination of these dilepton results with those from a search in the diphoton final state, which increases the limit by ∼ 0.1 TeV [38]. The following sections describe the first virtual graviton exchange search performed by ATLAS using dilepton data and its combination with an ATLAS diphoton data search [42].

II. ATLAS DETECTOR
ATLAS is a multipurpose particle detector composed of three main subsystems: the inner tracking detector, the calorimeter system and the muon spectrometer. The inner detector is used to track charged particles within a pseudorapidity η 1 in the range |η| < 2.5. It comprises a silicon pixel detector, a silicon-strip tracker and a transition radiation tracker. An axial 2 T magnetic field is produced by a superconducting solenoid.
The calorimeter system, covering |η| < 4.9, surrounds the solenoid and provides three-dimensional reconstruction of electromagnetic and hadronic showers. The lead/liquid-argon electromagnetic sampling calorimeter covers |η| < 2.5 and is finely segmented with a readout granularity varying by layer and with cells as small as 0.025 × 0.025 in (η, φ) to provide precise energy and position resolution, as needed for electron and photon identification and energy measurement. Hadron calorimetry is provided by an iron/scintillator tile calorimeter in the central pseudorapidity range |η| < 1.7 and a lead/liquidargon calorimeter extending the pseudorapidity range up to |η| = 3.2. Both the electromagnetic and hadronic calorimeters have liquid-argon-based forward detectors, with copper or tungsten as an absorber, to extend coverage up to |η| = 4.9.
Outermost is the muon spectrometer, another key detector component for this analysis. Three layers of precision tracking chambers, comprising monitored drift tubes and cathode strip chambers, enable muon reconstruction up to |η| = 2.7. The magnetic field is provided by three large air-core toroidal magnet systems (one barrel and two end caps), each consisting of eight azimuthally symmetric superconducting coils. Triggering capability up to |η| = 2.4 is provided by fast resistive plate chambers in the barrel and thin-gap chambers in the end caps.

III. SIGNAL AND BACKGROUND MODELING
The dominant background contribution comes from the SM DY process with smaller contributions from tt and electroweak diboson (W W , W Z, and ZZ) production. In the dielectron channel, there is also a significant background from multijet and W +jets events in which jets are misidentified as electrons. Backgrounds are estimated using fully simulated Monte Carlo (MC) samples except for the combined multijet and W +jets background, which is determined from the data.
DY samples are generated with Pythia 6.421 [43] using MRST2007 LO * * parton distribution functions (PDFs) [44]. The diboson background is generated with Herwig 6.510 [45] using MRST2007 LO * * PDFs. For the tt background, event generation is performed with Mc@nlo 4.01 [46] and the CTEQ 6.6 PDFs [47], as well as Herwig to model the underlying event and parton showers. Production of diboson and tt events relies on Jimmy 4.31 [48] to describe multiple interactions.
For the contact interactions analysis, Pythia 6.421 and the MRST2007 LO * * PDFs are used to generate Cylindrical coordinates (r, φ) are used in the transverse plane, φ being the azimuthal angle around the beam pipe. The pseudorapidity is defined in terms of the polar angle θ as η = − ln tan(θ/2). both DY and CI as part of the same process in order to properly handle the interference between DY and CI contributions, as shown in Eq. (2). All quark flavors contribute to the DY+CI production.
Final-state radiation is simulated with Photos [49] for all the samples listed above. Higher-order QCD corrections are implemented via mass-dependent K factors defined as the ratio of the next-to-next-to-leading-order (NNLO) Z/γ * cross section calculated using Phozpr [50] and the MSTW2008 NNLO PDFs [51], to the LO Z/γ * cross section from Pythia. Higher-order electroweak corrections originating from virtual gauge-boson loops are computed using the Horace NLO event generator [52]. The mass-dependent QCD and electroweak K factors are applied to both DY and DY+CI samples.
For the large extra dimensions analysis, Sherpa 1.3.1 [53] and the CTEQ6L PDFs [54] are used to generate DY+ADD events at leading order. The generated dilepton mass is required to be less than the scale M S , which is chosen to be in the range between 2 and 3 TeV in this study, since the model is not valid at energies beyond that scale.
Finally, the generated samples are processed through a full simulation of the ATLAS detector [55] based on Geant 4 [56] and reconstructed with the same software as for the data. Several corrections derived from data control samples are applied to the simulated samples. Specifically, the energy scale and resolution for electrons are corrected so that the Z → ee mass distribution in simulation matches the data [57]. Similarly, the muon transverse momentum scale and resolution are adjusted to reproduce the muon tracking performance measured in Z → µµ data as well as several dedicated data samples [58]. The effect of pileup (multiple pp interactions in the same or nearby bunch crossings) is included by superimposing minimum bias events, in the same proportion as in data, on top of the hard scattering process generated. Small corrections are included in the analysis to properly describe the pileup conditions for the selected data sample.

IV. EVENT SELECTION AND COMPARISON BETWEEN EXPECTED AND OBSERVED YIELDS
This analysis follows the same event selection as the search for new heavy resonances [59] in the dielectron channel but uses a somewhat tighter selection in the dimuon channel.
The data sample was collected during LHC operation in 2011 and corresponds to a total integrated luminosity of 4.9 and 5.0 fb −1 in the e + e − and µ + µ − final states, respectively. The events recorded by the ATLAS detector were selected by requiring that they pass specific triggers. The trigger for the dielectron dataset required the presence of two electromagnetic clusters consistent with originating from electrons with transverse momentum p T above 20 GeV, whereas events in the muon dataset were required to pass at least one of two single-muon triggers with p T thresholds of 22 GeV and 40 GeV.
After passing the trigger selection, events are required to have a pair of either electrons or muons with p T greater than 25 GeV. Furthermore, events are required to be recorded during stable beam conditions and with detector components operational. To reject cosmic ray events and beam halo background, events are required to have a reconstructed vertex with at least three charged particle tracks with p T > 0.4 GeV. If more than one such vertex is found, the vertex with the largest Σp 2 T is selected as the primary vertex of the event, where the sum is over all charged particles associated with the given vertex. Electron candidates are confined to |η| < 2.47, with the calorimeter barrel-to-end-cap transition region 1.37 ≤ |η| ≤ 1.52 excluded due to the degraded energy resolution in this region. No explicit η requirement is placed on muon candidates, but the selection described below leads to negligible acceptance beyond |η| of approximately 2.5.
Electron candidates are formed from clusters of cells in the electromagnetic calorimeter where energy is deposited. Identification criteria based on the transverse shower shape, the leakage into the hadronic calorimeter, and the association to an inner detector track are applied to the cluster to satisfy the medium electron definition [57]. The electron energy is obtained from the calorimeter measurements and its direction from the associated inner detector track. A hit in the first layer of the pixel detector is required (if an active pixel module is traversed) to suppress background from photon conversions. Further jet background suppression is achieved by demanding that the highest-p T electron in the event be isolated. To this effect, the sum of the transverse energies, ΣE T , in calorimeter cells within a radius R = (∆η) 2 + (∆φ) 2 of 0.2 around the electron direction, is required to be less than 7 GeV. The core of the electron energy deposition is excluded and the sum is corrected for transverse shower leakage and pileup. The two electron candidates are not required to have opposite charge because of possible charge misidentification either due to bremsstrahlung or limited momentum resolution of the inner detector at high p T . If the event contains more than two selected electrons, the two electrons with the highest-p T sum are chosen. For these selection criteria, the overall event acceptance for DY events has a small dependence on the dielectron mass above 500 GeV, with a value of approximately 65% at 1 TeV.
Muon candidates are reconstructed independently in the inner detector and the muon spectrometer. The momentum is taken from a combined fit to the measurements from the two subsystems. To obtain optimal momentum resolution and accurate modeling by the simulation, muon candidates are required to have at least three hits in each of the inner, middle, and outer detector layers of the muon spectrometer, and to have at least one hit in each of two different layers in the nonbending xy plane. To suppress background from cosmic rays, requirements are imposed on the primary vertex (PV) position and the muon impact parameter relative to the PV: z coordinate of the PV |z PV | < 200 mm, muon transverse impact parameter |d 0 | < 0.2 mm and muon z coordinate |z 0 − z PV | < 1 mm. Furthermore, the muons are required to be isolated to reduce background from jets: Σp T (R < 0.3)/p T (µ) < 0.05, where the sum is over inner detector tracks within a radius of 0.3 around the muon direction. If more than one opposite-sign muon pair is found in an event, the pair with the highest-p T sum is chosen. The overall event acceptance for DY events has only a weak dependence on the dimuon mass, with a value of approximately 40% at 1 TeV. This is lower than the acceptance in the dielectron channel primarily due to the stringent requirements on the presence of hits in all three layers of the muon spectrometer and the extent of the three-layer geometrical coverage.
The W +jets background in the dimuon channel is estimated from simulated samples and is found to be negligible since the event must contain two well-measured highp T isolated muons. Likewise, the multijet background, estimated directly from the data by reversing the muon isolation criterion, is found to be negligible. The multijet and W +jets backgrounds are not negligible in the dielectron channel. They are estimated primarily from the data using several methods [59]. The first method determines the multijet background from the data and relies on the MC simulation for the W +jets contribution. The background is measured with a template built by reversing one of the electron identification criteria and normalized to data in the range 70 < m ee < 200 GeV. Another independent method that is sensitive to both multijet and W +jets backgrounds uses jet-enriched data samples either from jet triggers or from the same trigger used to select the events in this analysis. The method relies on jet misidentification rates, defined as the number of jets that pass the full electron selection divided by the number that pass a loose electron selection obtained by reversing one of the identification criteria. The background estimate is then constrained by a fit in the range 140 < m ee < 850 GeV. The final combined multijet and W +jets background is obtained with a simple average of the expected event yields from the different methods.
Extensive comparisons between data and MC simulation were performed at the level of single-lepton distributions to confirm that the simulation reproduces the selected data, especially at high momentum. Figure 1 shows good data-MC agreement in the lepton transverse momentum distributions for events passing all selection criteria. Figure 2 shows the dielectron and dimuon mass distributions for selected events. Also shown are the predicted contributions from SM and new phenomena (NP) for several choices of model parameters. The expected SM distribution is dominated by the DY process over the entire mass range and is found to describe the data well. The level of agreement with the SM expectation is also illustrated in Fig. 3, which shows the number of events above a minimum mass m min ℓℓ . A more quantitative comparison is provided in Tables I and II showing the numbers of observed and expected events in the dielectron and dimuon channels, respectively. The expected yields are normalized to the number of events observed in the Z peak control region (70 < m ℓℓ < 110 GeV). The mass region shown in these tables corresponds to the CI search region defined by m ℓℓ > 400 GeV. These tables also display the expected yields for the SM+CI signal for the two scenarios where the CI interferes either constructively or destructively with the SM.

V. SYSTEMATIC UNCERTAINTIES
Except for the multijet and W +jets background contributions to the dielectron channel, all signal and background event yield estimates are based on MC simulation. Because these yields are normalized in the Z peak control region, only mass-dependent systematic uncertainties affect the event yield estimates in the higher-mass signal region. The only exception is a 5% uncertainty applied to the signal yield to account for the uncertainty in the Z/γ * cross section which affects the signal normalization.
Experimental uncertainties arise from lepton energy/momentum scale and resolution, as well as trigger, reconstruction and identification efficiencies. In the dielectron channel, the largest experimental uncertainty comes from the combined multijet and W +jets background estimate. It is determined from the envelope of the three separate methods used, including the effect of varying the mass ranges in the background fits and the uncertainties in the η and p T dependence of the jet misidentification rates. Electron energy scale and resolution are determined from data via J/ψ → ee and Z → ee mass distributions, as well as studies of electron E/p in W → eν decays [57]. The uncertainty in the constant term that dominates the resolution at high energy has negligible impact on this analysis. A somewhat larger impact comes from the energy scale knowledge, resulting in a systematic error of 1.2% and 2.4% for dielectron masses of 1 and 2 TeV, respectively. A slight efficiency drop of 1.0% per TeV is predicted by the simulation due to the isolation requirement on the leading electron. To account for this effect, an uncertainty of the same magnitude is introduced.
In the dimuon channel, the largest contribution to the experimental systematic error comes from the muon reconstruction efficiency and muon resolution. A slight drop in reconstruction efficiency is predicted by the simulation at high p T due to the presence of additional hits in the muon spectrometer from muons undergoing large energy loss in the detector. An uncertainty of 3.0% (6.0%) at a dimuon mass of 1 (2) TeV is assessed, corresponding to the magnitude of this effect. The limited knowledge of the momentum scale determined from Z → µµ data has a negligible impact on the analysis. The momentum resolution in the simulation is adjusted based on Z → µµ and W → µν data, as well as dedicated straight muon track data collected with the toroids turned off and tracks passing through overlapping sectors in the muon spectrometer. The latter provide two independent momentum measurements for the same muon. The toroidoff and overlapping sector tracks are key to determining the muon reconstruction performance at high p T . The uncertainty in the muon resolution, taken as equal in magnitude to the correction applied to the simulation, results in a change in the event yield of 1.2% (12%) for m µµ = 1 (2) TeV.
The largest error contribution due to theory arises from limited knowledge of the PDFs, α S , and QCD K factors. Scale uncertainties are computed by taking the maximum deviations obtained by independently varying the renormalization (µ R ) and factorization (µ F ) scales by a factor of 2 but with the constraint that the ratio µ F /µ R does not change by more than a factor of 2. The α S and PDF uncertainties are determined with the MSTW2008NNLO eigenvector PDF sets and the different PDFs corresponding to variations of α S . The overall uncertainty is computed using 90% C.L. ranges and includes the envelope of the uncertainty bands for the following different PDF sets: MSTW2008, NNPDF2.1, CT10, and CT10W. PDFs are the largest source of uncertainty, with the envelope of all PDFs considered becoming the dominant contribution above a few hundred GeV. Uncertainties in the electroweak K factor [60] originate from the calculation of real boson radiation, O(αα S ) corrections, higher-order electroweak corrections, an assumed uncertainty of 10% in the contribution from photon-induced processes, and a difference in the electroweak renormalization scheme definition used in Pythia and in the calculation of the electroweak corrections with Horace. The latter source is the largest contribution to the electroweak uncertainty.
The systematic uncertainties are summarized in Table III. Although not explicitly listed in this table, the uncertainty due to limited MC statistics is also taken into account in the limit setting. For DY+CI MC samples, this uncertainty grows from about 4% at low m ℓℓ to about 30% at high m ℓℓ for Λ = 12 TeV.

VI. STATISTICAL ANALYSIS
The data analysis proceeds with a Bayesian method to compare the observed event yields with the expected yields for a range of different NP model parameters (where the NP corresponds to either contact interactions or large extra dimensions). Specifically, the number of expected events in a given search region is where n DY+NP (θ,ν) is the number of events predicted by the DY+NP simulation for a particular choice of NP model parameter θ, n non−DY bg (ν) is the number of non-DY background events, andν represents the set of Gaussian nuisance parameters that account for systematic uncertainties. The parameter θ corresponds to a choice of energy scale Λ and interference parameter η LL in the CI analysis or to a choice of string scale M S and formalism in the ADD analysis. In the case of the CI analysis, the input to evaluate the complete set of µ values is shown in Tables I and II for the dielectron and dimuon channels, respectively. For each mass bin, a second-order polynomial is used to model the dependence of µ on 1/Λ 2 . In the case of the ADD analysis, µ is also parameterized by a second-order polynomial but as a function of 1/M 4 S . The likelihood of observing a set ofn events in N invariant mass bins is given by a product of Poisson probabilities for each mass bin k: According to Bayes' theorem, the posterior probability for the parameter θ givenn observed events is where Z is a normalization constant and the marginalized likelihood L M corresponds to the likelihood after all nuisance parameters have been integrated out. This integration is performed assuming that the nuisance parameters are correlated across all mass bins; Table III shows which parameters are taken into account for either or both of the signal and background expectations. The prior probability P (θ) is chosen to be flat in either 1/Λ 2 or 1/Λ 4 for the CI analysis, and either 1/M 4 S or 1/M 8 S for the ADD analysis. These choices are motivated by the form of Eqs. (2) and (5). The 95% C.L. limit is then obtained by finding the value θ lim satisfying θ lim 0 P(θ |n) dθ = 0.95, where θ is chosen to be 1/Λ 2 , 1/Λ 4 , 1/M 4 S or 1/M 8 S . The above calculations have been performed with the Bayesian Analysis Toolkit (BAT) [61], which uses a Markov chain Monte Carlo technique to integrate over nuisance parameters.

VII. CONTACT INTERACTIONS ANALYSIS AND RESULTS
To test the consistency between the data and the SM in the CI search region (m ℓℓ > 400 GeV), a likelihood ratio test is performed by producing a set of SM-like pseudoexperiments and comparing the likelihood ratio between the signal+background and pure background hypotheses obtained in the data to the results of the pseudoexperiments. The signal+background likelihood is evaluated at the Λ value that maximizes it. The derived p-value, corresponding to the probability of observing a fluctuation in the pseudoexperiments that is at least as signal-like as that seen in the data (i.e., with a maximum likelihood ratio greater than or equal to that obtained in the data), is estimated to be 15% (76%) in the dielectron channel and 79% (59%) in the dimuon channel for constructive (destructive) interference. These values indicate that there is no significant evidence for contact interactions in the analyzed data and thus limits are set on the contact interaction scale Λ.
Limits are obtained with the Bayesian method described above. Electroweak corrections are applied to both DY and DY+CI samples for consistency, although part of the electroweak corrections cannot be computed reliably due to the unknown new phenomena represented by the contact interaction. This particular choice results in slightly more conservative limits.
The expected 95% C.L. lower limit values on the energy scale Λ are found to be 13.8 ± 1.7 TeV for constructive interference (Λ − ) and 10.4 ± 1.0 TeV for destructive interference (Λ + ) in the dielectron channel. The corresponding expected limits in the dimuon channel are 12.7 ± 1.5 TeV and 9.9 ± 1.1 TeV. The quoted uncertainties correspond to the 68% range of limits surrounding the median value (taken to be the expected limit) of all limits obtained with a set of pseudoexperiments. Limits are expected to be stronger in the dielectron channel than in the dimuon channel due to the significantly larger acceptance for the dielectron selection.
The observed limits (at 95% C.L.) are Λ − > 12.1 TeV and Λ + > 9.5 TeV in the dielectron channel for constructive and destructive interference, respectively. The corresponding limits in the dimuon channel are Λ − > 12.9 TeV and Λ + > 9.6 TeV. These limits are summarized in Table IV.
If instead of choosing the prior to be flat in 1/Λ 2 , it is selected to be flat in 1/Λ 4 to match the form of the pure CI term in Eq. (2), the observed limit in the dielectron channel becomes weaker by 0.7 TeV for constructive interference and 0.4 TeV for destructive interference. The corresponding respective shifts to lower values are 1.2 and 0.6 TeV in the dimuon channel, see Table IV.
Finally, a limit is set for the combination of the dielec- tron and dimuon channels, assuming lepton universality, by computing a combined posterior probability for the two channels. The following sources of systematic uncertainty are treated as fully correlated between the two channels: PDF and α S , QCD and electroweak K factors, and Z/γ * cross section for normalization. All other sources are treated as uncorrelated. The resulting combined limits are Λ − > 13.9 TeV and Λ + > 10.2 TeV for the 1/Λ 2 prior. Table IV summarizes all limits for the two priors considered in this analysis.

VIII. LARGE EXTRA DIMENSIONS ANALYSIS AND RESULTS
The search for large extra dimensions is carried out similarly to that for contact interactions. A difference from the CI analysis is that the DY component present in the Sherpa DY+ADD simulated samples is subtracted out to compute the net ADD contribution to the total event yield. The DY background is modeled with the same Pythia DY sample as is used for the CI analysis. Another difference is that the search is performed in only one mass bin with the minimum mass chosen at the value giving the strongest expected limit. This optimization results in a signal region with a minimum mass requirement of 1300 GeV as determined from a set of pseudoexperiments in each of the dielectron and dimuon channels. Table V presents the expected and observed event yields in the signal region, including the expectation for several M S values in the GRW formalism.
The consistency between the number of observed events in the data and the predicted SM contribution is assessed using a set of SM-like pseudoexperiments. Using the same likelihood ratio approach as for the CI analysis, p-values of 6% and 68% are obtained in the dielectron and dimuon channels, respectively. These values indicate that there is no significant evidence for large extra dimensions and thus limits are set on the scale M S . The observed limits are M S > 2.73 (2.62) TeV in the dielectron chan-  Table VI shows these observed limits along with the expected limits. Limits in the dielectron channel are slightly worse than expected due to the larger number of events observed in the data compared with the SM expectation. The dielectron and dimuon channels are combined taking correlated systematic uncertainties into account in a way identical to the CI analysis.
A search for large extra dimensions has also been carried out in the diphoton final state using the data sample collected by ATLAS in 2011 [42]. The results of that search are combined with the dilepton results presented here with the use of BAT. Correlated sources of systematic uncertainty are treated as follows. The PDF uncertainty in the SM diphoton and DY production is considered to be fully correlated between the ee, µµ and γγ channels, whereas the multijet background uncertainty is fully correlated between the ee and γγ final states.
It should be noted that the ee and γγ samples are statistically uncorrelated since ee candidates have been explicitly removed from the diphoton analysis at the event selection stage. The observed and expected combined limits are given in Table VI, Table VII. Limits are also obtained with a K factor applied to the ADD signal yield to account for next-to-leading-order QCD corrections. A constant K factor of 1.6 is applied in the dilepton channel [62] and 1.7 in the diphoton channel [63].  Table VIII.

IX. CONCLUSIONS
A search for contact interactions and large extra dimensions has been performed in dielectron and dimuon events produced in LHC proton-proton collisions at √ s = 7 TeV. The data sample corresponds to an integrated luminosity of 4.9 (5.0) fb −1 of pp collisions in the dielectron (dimuon) channel recorded with the ATLAS detector. No significant deviation from the Standard Model is observed in the dilepton mass distributions. Using a Bayesian approach with a prior flat in 1/Λ 2 , as was done in most previous searches at hadron colliders, the following 95% C.L. limits are set on the energy scale of contact interactions: Λ − > 12.1 TeV (Λ + > 9.5 TeV) in the di- The crucial computing support from all WLCG partners is acknowledged gratefully, in particular from CERN and the ATLAS Tier-1 facilities at TRIUMF (Canada), NDGF (Denmark, Norway, Sweden), CC-IN2P3 (France), KIT/GridKA (Germany), INFN-CNAF (Italy), NL-T1 (Netherlands), PIC (Spain), ASGC (Taiwan), RAL (UK) and BNL (USA) and in the Tier-2 facilities worldwide.