Measurement of differential cross sections for the production of a pair of isolated photons in pp collisions at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}=7\,\text {TeV} $$\end{document}s=7TeV

A measurement of differential cross sections for the production of a pair of isolated photons in proton–proton collisions at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}=7\,\text {TeV} $$\end{document}s=7TeV is presented. The data sample corresponds to an integrated luminosity of 5.0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\,\text {fb}^{-1}$$\end{document}fb-1 collected with the CMS detector. A data-driven isolation template method is used to extract the prompt diphoton yield. The measured cross section for two isolated photons, with transverse energy above 40 and 25\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\,\text {GeV}$$\end{document}GeV respectively, in the pseudorapidity range \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|\eta |<2.5$$\end{document}|η|<2.5, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|\eta |\notin [1.44,1.57]$$\end{document}|η|∉[1.44,1.57] and with an angular separation \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta R > 0.45$$\end{document}ΔR>0.45, is \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$17.2 \pm 0.2\,\text {(stat)} \pm 1.9\,\text {(syst)} \pm 0.4\,\text {(lumi)} $$\end{document}17.2±0.2(stat)±1.9(syst)±0.4(lumi) \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {\,pb}$$\end{document}\,pb. Differential cross sections are measured as a function of the diphoton invariant mass, the diphoton transverse momentum, the azimuthal angle difference between the two photons, and the cosine of the polar angle in the Collins–Soper reference frame of the diphoton system. The results are compared to theoretical predictions at leading, next-to-leading, and next-to-next-to-leading order in quantum chromodynamics.


Introduction
The measurement of differential diphoton production cross sections offers an important test of both perturbative and non-perturbative quantum chromodynamics (QCD). At leading order (LO), diphotons are produced via quark-antiquark annihilation qq → γ γ . At next-to-leading order (NLO), diphoton production also includes the quark-gluon channel, while next-to-next-to-leading order (NNLO) adds the gluongluon channel, which includes a box diagram and represents a non-negligible fraction of the total cross section. Diphoton production is sensitive to the emission of soft gluons in the initial state and to the non-perturbative fragmentation of quarks and gluons to photons in the final state. Due to this rich phenomenology, theoretical predictions are challenging especially in restricted regions of phase space.
Diphoton production constitutes the major source of background in the diphoton decay channel of the newly discovered * e-mail: cms-publication-committee-chair@cern.ch Higgs boson [1][2][3], as well as to searches for physics beyond the standard model. New physics processes may also appear as non-resonant deviations from the predicted diphoton spectrum in events with large missing transverse energy, as in gauge-mediated SUSY breaking [4] or in models of universal extra dimensions [5]. Alternatively, some models predict narrow resonances, such as the graviton in the Randall-Sundrum model for warped extra dimensions [6,7].
The most recent diphoton measurements were published by the CDF and D0 Collaborations [8,9] at the Tevatron and by the ATLAS Collaboration [10] at the LHC. This paper presents an update of a previous CMS measurement at √ s = 7 TeV [11] and is based on the full 2011 data sample of 5.0 fb −1 . It probes a phase space defined by a highly asymmetric selection for the transverse energy (E T ) of the two photons. The leading-order kinematic configuration where photons are produced back-to-back in the transverse plane is suppressed, enhancing the sensitivity to higher-order diagrams. The ratio of the NNLO to the LO prediction is increased by 20 % with respect to the previous CMS measurement [11].
The main experimental challenge for the measurement of the diphoton cross section is distinguishing the "prompt" photon signal produced either directly or as a result of fragmentation from the background that arises mainly from energetic neutral mesons, predominantly π 0 and η mesons, inside jets. These mesons typically decay to two collimated photons that are reconstructed as a single photon candidate, which is referred to as "non-prompt" in this paper. The main features used to discriminate a prompt photon from a non-prompt one are the shape of the shower measured by the electromagnetic calorimeter (ECAL) [12] and the isolation energy in a cone around the photon direction [11,13]. This information can be used to trigger on diphoton candidate events and, at the analysis level, to statistically evaluate the fraction of prompt diphoton candidates.
The particle flow (PF) event reconstruction [14] consists in reconstructing and identifying each particle with an optimal combination of all sub-detector information. In this pro-cess, the identification of the particle type (photon, electron, muon, charged hadron, neutral hadron) plays an important role in the determination of the particle direction and energy.
In this analysis, the photon component of the PF isolation is used as the discriminating variable. The distributions of signal and background components are built from data and used in a maximum likelihood fit to estimate the signal fraction. An increased separation power with respect to previous results [11] is achieved by improving the identification and subtraction of the photon energy deposit in the isolation cone.
After a brief description of the CMS detector in Sect. 2, data and simulated samples are presented in Sect. 3, and the photon reconstruction in Sect. 4. The diphoton signal is estimated as described in Sect. 5. The number of signal events is then corrected for inefficiencies and unfolded as described in Sect. 6. Systematic uncertainties are assessed in Sect. 7, and the differential cross sections are presented and compared to theoretical predictions in Sect. 8.

The CMS detector
A detailed description of the CMS detector can be found elsewhere [15]. Its central feature is a superconducting solenoid, 13 m in length and 6 m in diameter, which provides an axial magnetic field of 3.8 T. The bore of the solenoid is instrumented with both the tracker (TRK) and the calorimeters. The steel flux-return yoke outside the solenoid is instrumented with gas-ionisation detectors used to reconstruct and identify muons. Charged-particle trajectories are measured by the silicon pixel and strip tracker, with full azimuthal (φ) coverage within |η| < 2.5, where the pseudorapidity η is defined as η = − ln[tan(θ/2)], with θ being the polar angle of the trajectory of the particle with respect to the counterclockwise beam direction. A lead tungstate crystal electromagnetic calorimeter (ECAL) and a brass/scintillator hadron calorimeter (HCAL) surround the tracking volume and cover the region |η| < 3. The ECAL barrel (EB) extends to |η| < 1.479 while the ECAL endcaps (EE) cover the region 1.479 < |η| < 3.0. A lead/silicon-strip preshower detector (ES) is located in front of the ECAL endcap in the region 1.653 < |η| < 2.6. The preshower detector includes two planes of silicon sensors measuring the x and y coordinates of the impinging particles. In the (η, φ) plane, and for |η| < 1.48, the HCAL cells map onto 5 × 5 ECAL crystal arrays to form calorimeter towers projecting radially outwards from points slightly offset from the nominal interaction point. In the endcap, the ECAL arrays matching the HCAL cells contain fewer crystals. A steel/quartz-fibre Cherenkov forward calorimeter extends the calorimetric coverage to |η| < 5.0.

Data sample
The data sample consists of proton-proton (pp) collision events collected at the LHC with the CMS detector in the year 2011, at a centre-of-mass energy ( √ s) of 7 TeV and corresponding to an integrated luminosity of 5.0 fb −1 .
Events are triggered [15] by requiring the presence of two photons with asymmetric transverse energy thresholds. The E T thresholds at trigger level are 26 (18) and 36 (22) GeV on the leading (sub-leading) photon, depending on the running period. Each candidate is required to satisfy either loose calorimetric identification requirements, based on the shape of the electromagnetic shower, or loose isolation conditions. The trigger efficiency is evaluated using a tag-and-probe technique on Z → e + e − events [16], with electrons treated as photons. The trigger efficiency for photons selected in this analysis is measured to be between 98.8 and 100 % depending on the pseudorapidity and the interaction with the material in front of the ECAL. The total trigger efficiency is found to be constant over the data taking period.
Several samples of simulated events are used in the analysis to model signal and background processes. Drell-Yan+jets and γ γ +jets signal events are generated with Mad-Graph 1.4.8 [17]. The gg→ γ γ box signal process, γ +jet, and QCD dijet background processes are generated with pythia 6.4.24 [18]. For all simulated samples the CTEQ6L1 [19] parton distribution functions (PDFs) are used. All generated events are then processed with pythia (Z2 tune) [20] for hadronization, showering of partons and the underlying event; a detailed simulation of the CMS detector based on Geant4 [21] is performed, and the simulated events are finally reconstructed using the same algorithms as used for the data.
The simulation includes the effects of in-time pileup (overlapping pp interactions within a bunch crossing) and out-oftime pileup (overlapping pp interactions from interactions happening in earlier and later bunch crossings) with a distribution matching that observed in data.

Photon reconstruction
Photon candidates are reconstructed from the energy deposits in the ECAL by grouping its channels into superclusters [22]. About half of the photons convert into an e + e − pair in the material in front of the ECAL. Conversion-track pairs are reconstructed from a combination of Gaussian-sum filter (GSF) electron tracks [23] and ECAL-seeded tracks fit to a common vertex and then matched to the photon candidate. The superclustering algorithms achieve an almost complete collection of the energy of such converted photons. In the barrel region, superclusters are formed from five-crystalwide strips in η, centred on the locally most energetic crystal (seed), and have a variable extension in the azimuthal direction (φ). In the endcaps, where the crystals are arranged according to an x-y rather than an η-φ geometry, matrices of 5 × 5 crystals (which may partially overlap) around the most energetic crystals are merged if they lie within a narrow φ road. The photon candidates are reconstructed within the ECAL fiducial region |η| < 2.5 but excluding the barrelendcap transition regions 1.44 < |η| < 1.57. This exclusion of the barrel-endcap transition regions ensures containment of the shower of the selected photon candidate in either the ECAL barrel or one of the ECAL endcaps. The fiducial region requirement is applied to the supercluster position (defined as the log-weighted barycentre of the supercluster's active channels) in the ECAL.
The photon energy is computed starting from the raw crystal energies measured in the ECAL. In the region covered by the preshower detector the energy recorded in that subdetector is added. The variation of the crystal transparency during the run is continuously monitored and corrected using a factor based on the change in response to light from a laser and light-emitting-diode based monitoring system. The single-channel response of the ECAL is equalised by exploiting the φ symmetry of the energy flow, the mass constraint on the energy of the two photons in decays of π 0 and η mesons, and the momentum constraint on the energy of isolated electrons from W and Z decays. A correction factor compensates for the imperfect containment of the shower in the cluster crystals. The absolute energy scale and the residual long term drifts in the response are further corrected using Z → e + e − decays [22].
Interaction vertices are reconstructed from charged tracks and the vertex of the diphoton event is taken as the one with the largest sum of squared transverse momenta ( p 2 T ) of the associated tracks. The photon four-momentum is recalculated with respect to this vertex.

Photon selection
The photon candidates are first required to pass a sequence of filters that aim to remove beam backgrounds or identified detector issues and to satisfy more stringent criteria than the trigger requirements. The preselection is based on the shape of the electromagnetic shower in the ECAL and on the degree of isolation of the photon (i.e. the amount of energy deposited in the vicinity of the photon). The variables used are: -Photon supercluster raw energy E raw SC : the sum of the calibrated crystal energies; -Preshower energy E ES SC : the sum of the energy deposits reconstructed in the preshower detector (ES) and associated with the supercluster; -R 9 : the energy sum of 3 × 3 crystals centred on the most energetic crystal in the supercluster divided by the raw energy of the supercluster; -H/E: the ratio of the energy deposited in HCAL that is inside a cone of size R = ( η) 2 + ( φ) 2 = 0.15 centred on the photon direction, to the supercluster energy; -σ ηη : the shower transverse extension along η that is defined as: where the sum runs over all elements of the 5 × 5 matrix around the most energetic crystal in the supercluster, and η i = 0.0174η i in EB, η i = 0.0447η i in EE withη i denoting the index of the ith crystal along the η direction. The individual weights w i are given by where E i is the energy of the ith crystal andη = η i E i / E i is the weighted average pseudorapidity; -Iso 0.3 ECAL (ECAL isolation): the scalar sum of the E T of the deposits in the electromagnetic calorimeter lying inside a cone of size R = 0.3, centred on the direction of the supercluster but excluding an inner cone of size 3.5 crystals and an η-slice region of 2.5 crystals; -Iso 0.3 HCAL (hadronic calorimeter isolation): the scalar sum of the E T of the deposits in the hadron calorimeter that lie inside a hollow cone of outer radius of size R = 0.3 and inner radius of size R = 0.15 in the η-φ plane, centred on the direction of the supercluster; -Iso 0.3 TRK (tracker isolation): the scalar sum of the p T of the tracks that are consistent with originating from the primary vertex in the event, and lie inside a hollow cone of outer radius of size R = 0.3 and inner radius of size R = 0.04 in the η-φ plane, centred around a line connecting the primary vertex with the supercluster but excluding an η-slice region ( η = 0.015).
The isolation requirements are kept loose because the isolation is used as the discriminating variable in the signal extraction procedure. The selection criteria are defined to be slightly tighter than the trigger selection. The shower shape variables in the simulation are corrected to compensate for their imperfect modeling, mainly connected with (a) the simulation of effective readout noise in ECAL channels, (b) the effect of overlapping energy deposits from collisions in adjacent bunch crossings, and (c) the description of the material budget in the detector geometry. The correction factors are extracted from a sample of photons in Z → μ + μ − γ events, and validated as a function of E T and η in a sample of electrons from Z boson decays. The list of preselection criteria is presented in Table 1. The preselected photons must satisfy additional requirements to be considered as photon candidates. These consist of the absence of reconstructed electron track seeds in the pixel detector which match the candidate's direction, and a tighter selection on the hadronic leakage of the shower and the σ ηη shower shape variable. The list of additional selection criteria is shown in Table 2.
In the simulation, prompt photons are defined as candidates satisfying the analysis selection requirements and geometrically matched to an isolated generator-level photon, either directly produced or originating from a fragmentation process. The generator-level isolation is defined as the p T sum of stable particles in a cone of size R = 0.4, and is required to be less than 5 GeV.

Signal yield determination
The diphoton signal is extracted from events containing two photon candidates with transverse energy greater than 40 (25) GeV for the leading (sub-leading) photon, and with a separation of R > 0.45. If more than two photon candidates are selected, the two with highest E T are retained. The minimum separation requirement ensures that the energy deposit of one photon does not enter the isolation cone centered on the other one. The signal fraction is statistically separated from jets misidentified as photons by means of a binned maximum likelihood fit that uses the photon component of the PF isolation as the discriminating variable.
The diphoton signal is then studied as a function of the diphoton invariant mass m γ γ , the diphoton transverse momentum p γ γ T , the azimuthal angle difference φ γ γ between the two photons, and the cosine of the polar angle θ * in the Collins-Soper frame of the diphoton system [24]. A maximum likelihood fit is performed for each bin of the distributions in the above variables.

Particle flow isolation
The photon component of the PF isolation (Iso) is used to discriminate signal from background. The choice of the isolation variable is optimized to obtain the smallest total uncertainty of the measured cross section. This variable is computed, in a cone of size R = 0.4 around each selected photon candidate, as the E T sum of photons reconstructed with the PF algorithm [14]. The PF isolation deals more effectively with cases of overlapping particles than the calorimetry-based isolation.
When calculating the isolation, the energy deposited by the selected photon candidate is subtracted by removing from the cone the area where the photon is expected to have deposited its energy ("footprint"), since photon energy leaking into the cone could bias the isolation sum. This is done on an event-by-event basis relying on simple geometrical considerations. The directions of the momenta of reconstructed photon candidates around the selected photon are extrapolated from the interaction vertex to the inner surface of the ECAL, and whenever they overlap with a crystal belonging to the supercluster these photon candidates are removed from the isolation sum. For the matching between the propagated trajectory and the crystal front width, a tolerance of 25 % of the face size is applied.
This procedure does not use any generator-level information and can therefore be applied in both data and simulated events.
The pile-up introduces a spurious correlation between the two candidate photons' isolation sums. For this reason the PF isolation sums for both photons are corrected, event by event, for the presence of pile-up with a factor proportional to the average pile-up energy density (ρ) calculated with FastJet [25].

Template construction
The diphoton signal is extracted through a two-dimensional binned maximum likelihood fit that uses the isolation of the two selected photon candidates as discriminating variables. Different templates are built for the prompt-prompt ( f pp ), prompt-non-prompt ( f pn ), non-prompt-prompt ( f np ), and non-prompt-non-prompt ( f nn ) components in the (Iso 1 , Iso 2 ) plane, where Iso 1 and Iso 2 represent the isolation variables for the two selected photon candidates in the event. The probability distribution function has the following form: where T kk (Iso 1 , Iso 2 ) is the function describing the isolation distribution (template) for the component f kk . Techniques have been developed to extract the templates from data to avoid possible biases coming from an imperfect modeling of the events in the simulation. Samples of events where at least one photon passes the photon selection are used to create prompt-prompt, prompt-non-prompt, non-prompt-prompt and non-prompt-non-prompt templates with high statistical precision, as described in the following. The "random cone" technique is used to extract the prompt photon template with high statistical accuracy. In this procedure we compute the isolation energy in a region separated from the candidate photon. Starting from the photon (η, φ) axis, a new axis is defined at the same pseudorapidity η but with a random separation in azimuthal angle φ RC between 0.8 and 2π − 0.8 radians from the photon φ. This new axis is used to define the random cone provided that no jet with p T > 20 GeV or photon or electron with p T > 10 GeV is reconstructed within R < 0.8 and no muon is reconstructed within R < 0.4 from this axis. In the case where the new axis does not meet these requirements, a new azimuthal angle is generated. The isolation energy, which is defined as the energy collected in a cone of size R < 0.4 about the new axis once the fraction corresponding to the area of the photon supercluster has been removed, is then used to populate the prompt photon template.
The distribution of the template variable has been studied in Z → e + e − events and found to be in agreement with the template built with the random cone technique.
The background (non-prompt) template cannot be defined by simply inverting the photon preselection, because the candidates entering the analysis, i.e. fulfilling the preselection requirements, have "photon-like" characteristics, while the set of candidates not fulfilling the photon preselection criteria includes a large number of genuine jets. To avoid this bias, the candidates selected to populate the non-prompt photon template are chosen from those that fulfil all the photon selection criteria, except the σ ηη shower shape, which is not strongly correlated with the isolation variable as a result of the footprint removal technique described in the previous section. The events in a "sideband" close to the photon selection criterium are used to populate the non-prompt photon template. The sideband is defined as 0.011 < σ ηη < 0.014 for candidates reconstructed in the ECAL barrel and 0.030 < σ ηη < 0.034 for candidates reconstructed in the ECAL endcaps.  Fig. 1 Comparison of prompt photon templates in data and simulation: prompt photons in the simulation (squares), prompt photon templates extracted with the random cone technique from simulation (triangles) and from data (dots); (top) candidates in the ECAL barrel, (bottom) candidates in the ECAL endcaps. All histograms are normalized to unit area The same procedure (Sect. 5.1) is used for subtracting the pile-up energy from the photon isolation sums. The templates obtained using the random cone and the sideband techniques in the simulation are compared with the one-dimensional PF isolation distribution for prompt and non-prompt photons in simulated events and with the templates obtained from data (Figs. 1, 2).
The residual differences in the simulation between the isolation distribution and the templates defined with the random cone and the sideband techniques are accounted for as systematic uncertainties on the template shapes.
The two-dimensional templates are built selecting candidate photons from data with the same kinematics as the  Fig. 2 Comparison of non-prompt photon templates in data and simulation: non-prompt photons in the simulation (squares), non-prompt photon templates extracted with the sideband technique from simulation (triangles) and from data (dots); (top) candidates in the ECAL barrel, (bottom) candidates in the ECAL endcaps. All histograms are normalized to unit area diphoton events to be fitted. The procedure presented below correctly models the isolation distribution even in the case of overlap between the isolation cones of the two photon candidates.
The prompt-prompt template is built from events where the pileup energy density matches that of the event to be fitted, and where the two random cone directions are found having the same pseudorapidity and the same azimuthal angular separation as the selected photons.
The prompt-non-prompt template is built from events where a sideband photon is selected. The isolation sum around the sideband photon is used for the candidate to be fitted under the non-prompt-hypothesis. A direction satisfying the random cone criteria is then searched for in the same tem-plate event (oriented as the second candidate in the selected diphoton event) and used to calculate the isolation sum for the candidate to be fitted under the prompt hypothesis.
The non-prompt-non-prompt template is built selecting two events, each of which contains one sideband photon and such that their orientation matches the orientation of the candidate photons in the event to be fitted. Then, depending on the fraction of photon candidates with R γ γ < 1.0 present in the bin of the observable under analysis, a choice between two different strategies is made. If the fraction is below 10 %, the effect of the overlapping isolation cones can be neglected. The two-dimensional non-prompt-non-prompt template is then built by calculating each of the two isolation sums in the separate events. If the fraction is above 10 %, an additional requirement is imposed: the sum of the FastJet ρ of the two selected template events has to match the one of the diphoton event to be fitted. Then, the sets of reconstructed particles in the two template events are merged, and the isolation sums are calculated from this merged set of reconstructed particles along the direction of each sideband photon.
In this procedure, the pileup energy density of the template events is used to model the pileup energy density of the event to be fitted, and this allows us to describe the correlation between the isolation sums. The effect of the residual correlation mis-modeling is added to the template shape systematic uncertainty in the final result.

Fitting technique
The fit is performed separately for the cases where both candidates are reconstructed in the ECAL barrel, one in the ECAL barrel and one in the ECAL endcaps, or both in the ECAL endcaps. If both candidates are in the same detector region (EB-EB and EE-EE categories), the leading selected photon is assigned randomly to axis 1 or 2 of the two-dimensional plane, and the prompt-non-prompt ( f pn ) and non-promptprompt ( f np ) fractions are constrained to have the same value.
The fit, performed in each bin of the differential variables, is restricted to the region where the isolation of the photons is smaller than 9 GeV. To guarantee its stability even in the less populated bins, the fit is performed in steps. First the size of the bins in the two-dimensional plane (Iso 1 , Iso 2 ) is optimised to reduce statistical fluctuations of template shape in the tails; then a first fit is performed on the projections of the isolation distributions on the two axes of the plane using the one-dimensional templates described above. In a subsequent step, the fractions of prompt-prompt, prompt-nonprompt, non-prompt-prompt, and non-prompt-non-prompt, which are constrained to sum up to unity, are fit in the two-dimensional plane using as a constraint the results of the previous fit. The final likelihood maximisation is then performed after removing all constraints, and using as ini- An example of the first step of the procedure is obtained by fitting the one-dimensional projections of the isolation distributions as shown in Fig. 3. An example of the results of the final two-dimensional fit (projected on the axes for the sake of clarity) is shown in Fig. 4. The fractions of promptprompt, prompt-non-prompt, and non-prompt-non-prompt components are shown in Fig. 5 for the observables of the differential analysis. We fit about 69000 prompt diphoton events in the whole acceptance of the analysis.
The reported purity suffers from a contamination of electrons coming predominantly from Drell-Yan e + e − and incorrectly reconstructed as photons. The contamination is most significant in the Z peak region, where it reaches about 25 % of the raw diphoton yield. The fraction of electron pairs passing the analysis selection and contributing to the promptprompt fitted fraction is estimated from simulation, where correction factors are applied to obtain the electron to photon mis-identification probability measured in data, and used to subtract the contamination.

Efficiencies and unfolding
Following the methodology presented in the previous sections, a "raw" diphoton production cross section is extracted. To obtain the final result, this cross section is corrected for inefficiencies and finally unfolded from the measured quantities to the corresponding particle-level quantities.
The total diphoton efficiency can be separated into the trigger efficiency and reconstruction/selection efficiency, and can be written as: where trig is the trigger efficiency and sel is the diphoton reconstruction/selection efficiency from simulation. The factors C Z→e + e − The diphoton reconstruction/selection efficiency sel is about 85 % when both photons are in the barrel, 75 % when one photon is in the barrel and the other in one endcap, and 64 % when both photons are in the endcaps. All these correction factors are estimated from data and range from 0.99 to 1.02, depending on the photon E T and η.
The detector effects are unfolded from the measured yields for a direct comparison of experimental measurements with theoretical predictions. The number of unfolded diphoton events in each bin of the differential observables is The unfolding matrix is calculated using the iterative Bayesian technique [26,27]. The diphoton simulated sample from Mad-Graph hadronized with pythia is used. The distributions of diphoton candidates in the simulation are reweighted to the distributions of the raw diphoton yields from data as obtained from the fit procedure, for all the observables. The difference between the weighted and unweighted results is taken into account as a systematic uncertainty, and amounts to about 1 %. The unfolding correction amounts to 7 % of the raw yield at maximum, for the bins where the slope of the kinematic distributions is the steepest. Table 3 summarises the main sources of systematic uncertainty in the measurement of the integrated cross section.

Systematic uncertainties
The dominant uncertainty in the template shapes arises from the difference in shape between the templates built with the techniques described in Sect. 5 and the distributions of the isolation variable for prompt or non-prompt isolated photons for simulated events. The latter are used to generate data samples for each bin of the differential variables, with the fractions measured in data. Then, each of these datasets is fitted with templates built in the simulation with the same techniques used on data, and the average difference between the fitted fractions and those used for the generation is quoted as a systematic uncertainty. It amounts to 3 % (barrel template) and 5 % (endcap template) for the prompt component, and between 5 % (barrel template) and 10 % (endcap template) for the non-prompt component. The uncertainty in the template shape for fragmentation photons is evaluated in the simulation by doubling the probability of the fragmentation process, and that yields an additional 1.5 % uncertainty in the measured cross section. In the case of the non-prompt-non- prompt template, and only for the bins where a significant fraction of the diphoton candidates are close in R γ γ , an additional uncertainty ranging from 3 to 5 % is introduced to account for the imperfections on the template shape description due to the effect of ECAL noise and PF thresholds on the combination of two different events to build the template. The systematic uncertainty arising from the statistical uncertainty in the shape of the templates is evaluated generating modified templates, where the content of each bin is represented by a Gaussian distribution centred on the nominal bin value and with standard deviation equal to the statistical uncertainty of the bin. The root mean square of the distribution of the fitted purity values, divided by the purity measured with the original template, is used as systematic uncertainty in the purity measurement and amounts to about 3 %.
A possible bias associated with the fitting procedure is evaluated using pseudo-experiments. Pseudo-data samples are generated with given fractions of prompt-prompt, prompt-non-prompt, and non-prompt-non-prompt contributions, using the templates from simulation as generator probability density functions. Each data sample is then fitted with the same templates used for the generation. The average bias is negligible in all bins.
The systematic uncertainty associated with the subtraction of Drell-Yan e + e − events is evaluated by propagating the uncertainty in the electron to photon misidentification probability to the subtracted yield. The uncertainty in the fraction of such events that is fitted as prompt-prompt is also taken into account. This contribution is maximal for m γ γ close to the Z-boson mass. The relative contribution to the total systematic uncertainty is below 0.5 %.
The systematic uncertainty in the trigger efficiency is found to be below 0.5 %. The systematic uncertainty in the reconstruction and selection efficiencies is dominated by the uncertainty in the data-to-simulation corrections from the Z → e + e − and Z → μ + μ − γ control samples, and it ranges from 2 % in the barrel to 4 % in the endcap.
The systematic uncertainty in the integrated luminosity that corresponds to our data sample is 2.2 % [28].
The total systematic uncertainty in the measurement amounts to approximately 8 % when both candidates are reconstructed within the ECAL barrel, and to 11 % for the full acceptance of the analysis.  Fig. 6 The comparisons of the differential cross section between data and the sherpa, diphox + gamma2mc, resbos, and 2γ nnlo predictions for m γ γ . Black dots correspond to data with error bars including all statistical and systematic uncertainties. Only the scale uncertainty is included for the sherpa prediction. Scale, PDF and α S uncertainties are included for diphox + gamma2mc and resbos. Only statistical and scale uncertainties are included for the 2γ nnlo prediction (three extra jets) and with the box contribution at the matrix element level. The diphox NLO generator includes the direct and fragmentation contributions and uses a full fragmentation function for one or two partons into a photon at NLO. The direct box contribution, which is formally part of the NNLO corrections since it is initiated by gluon fusion through a quark loop, is computed at NLO with gamma2mc. The resbos NLO generator features resummation for Born and box contributions, and effectively includes fragmentation of one quark/gluon to a single photon at LO. The latter process is regulated to avoid divergences and does not include the full fragmentation function. The resbos p γ γ T spectrum benefits from a soft and collinear gluon resummation at nextto-next-to-leading-log accuracy. 2γ nnlo predicts the direct γ γ +X processes at NNLO. The sherpa sample is used after hadronization while diphox + gamma2mc, resbos, and 2γ nnlo are parton-level generators only and cannot be interfaced with parton shower generators.
The predictions have been computed for the phase space E γ 1 T > 40 GeV, E γ 2 T > 25 GeV, |η γ | < 1.44 or 1.57 < |η γ | < 2.5, R(γ 1 , γ 2 ) > 0.45. An isolation requirement is applied at the generator level. In sherpa, the E T sum of stable particles in a cone of size R = 0.4 has to be less than 5 GeV (after hadronization). In diphox, gamma2mc, and resbos the E T sum of partons in a cone of size R = 0.4 is required to be less than 5 GeV. In 2γ nnlo, the smooth Frixione isolation [35] is applied to the photons to suppress the fragmentation component: where E Iso T is the E T sum of partons in a cone of size R, R 0 = 0.4, = 5 GeV, and n = 0.05. This criterion, tested with diphox, is found to have the same efficiency as that used for the other generators within a few percent. A nonperturbative correction is applied to diphox, gamma2mc,  Fig. 7 The comparisons of the differential cross section between data and the sherpa, diphox + gamma2mc, resbos, and 2γ nnlo predictions for p γ γ T . Black dots correspond to data with error bars including all statistical and systematic uncertainties. Only the scale uncertainty is included for the sherpa prediction. Scale, PDF and α S uncertainties are included for diphox + gamma2mc and resbos. Only statistical and scale uncertainties are included for the 2γ nnlo prediction and 2γ nnlo predictions to correct for the fact that those generators do not include parton shower or underlying event contributions to the isolation cone. The fraction of diphoton events not selected due to underlying hadronic activity falling inside the isolation cone is estimated using the pythia 6.4.22 [18] event generator with tunes Z2, D6T, P0, and DWT [20]. A factor of 0.95 ± 0.04 is applied to the parton-level cross section to correct for this effect.
Theoretical predictions are performed using the CT10 [36] NLO PDF set for sherpa, diphox + gamma2mc, and resbos, and the MSTW2008 [37] NNLO PDF set for 2γ nnlo. The diphox and gamma2mc theoretical uncertainties are computed in the following way: the factorization and renormalization scales in gamma2mc are varied independently up and down by a factor of two around m γ γ (configurations where one scale has a factor of four with respect to the other one are forbidden). In diphox, the factorization, renormalization and fragmentation scales are varied in the same way. In resbos, the factorization and renormalization scales are varied simultaneously by a factor of two. The maximum and minimum values in each bin are used to define the uncertainty. In diphox, gamma2mc, and resbos, the 52 CT10 eigenvector sets of PDFs are used to build the PDF uncertainty envelope, also considering the uncertainty in the strong coupling constant α S , determined according to the CT10 α S PDF set. In 2γ nnlo, a simplified and less computationally intensive estimate of the renormalization and factorization scale uncertainties is performed by varying these scales simultaneously by a factor of two up and down around m γ γ ; no PDF uncertainty is computed. The same procedure is used in sherpa, using the internal METS scale, where scales are defined as the lowest invariant mass or negative virtuality in the core 2→2 configuration clustered using a k T -type algorithm.
The total cross section measured in data for the phase space defined above is: σ = 17.2 ± 0.2(stat.) ± 1.9(syst.) ± 0.4(lum.) pb, Fig. 8 The comparisons of the differential cross section between data and the sherpa, diphox + gamma2mc, resbos, and 2γ nnlo predictions for φ γ γ . Black dots correspond to data with error bars including all statistical and systematic uncertainties. Only the scale uncertainty is included for the sherpa prediction. Scale, PDF and α S uncertainties are included for diphox + gamma2mc and resbos. Only statistical and scale uncertainties are included for the 2γ nnlo prediction  Figures 6, 7, 8 and 9 show the comparisons of the differential cross section between data and the sherpa, diphox + gamma2mc, resbos, and 2γ nnlo predictions for the four observables.
The NLO predictions of diphox + gamma2mc are known to underestimate the data [11], because of the missing higherorder contributions. Apart from an overall normalization factor, the phase space regions where the disagreement is the largest are at low m γ γ , low φ γ γ . The resbos generator shows a similar trend, with a cross section closer to the data than diphox + gamma2mc; its prediction is improved at high φ γ γ due to soft gluon resummation. With higher-order diagrams included, 2γ nnlo shows an improvement for the overall normalization. It also shows a better shape description, especially at low φ γ γ , but it still underestimates the data in the same region. sherpa generally reproduces rather well the shape of the data, to a similar level as 2γ nnlo. One can note that 2γ nnlo and sherpa predict the p γ γ T ∼ 65 GeV observed in the data. This is expected since sherpa includes up to three extra jets at the matrix element level.

Summary
A measurement of differential cross sections for the production of a pair of isolated photons in pp collisions at √ s = 7 TeV has been presented. The data sample corresponds to an integrated luminosity of 5.0 fb −1 recorded in 2011 with the CMS detector. To enhance the sensitivity to higher-order diagrams, this measurement covers a phase space defined by an asymmetric E T selection by requiring two isolated photons with E T above 40 and 25 GeV respectively, in the pseudorapidity range |η| < 2.5, |η| / ∈ [1.44, 1.57] and with an angular separation R > 0.45.  Fig. 9 The comparisons of the differential cross section between data and the sherpa, diphox + gamma2mc, resbos, and 2γ nnlo predictions for |cos θ * |. Black dots correspond to data with error bars including all statistical and systematic uncertainties. Only the scale uncertainty is included for the sherpa prediction. Scale, PDF and α S uncertainties are included for diphox + gamma2mc and resbos. Only statistical and scale uncertainties are included for the 2γ nnlo prediction A data-driven method based on the photon component of the particle flow isolation has been used to extract the prompt diphoton yield. The isolation is calculated so that the energy leakage from the photon deposit inside the isolation cone is effectively subtracted. The measured total cross section is in agreement with the 2γ nnlo prediction. The sherpa and resbos predictions are compatible with the measurement within the uncertainties, while diphox + gamma2mc underestimates the total cross section. Differential cross sections for prompt diphoton production have been measured as a function of the diphoton invariant mass m γ γ , the diphoton transverse momentum p γ γ T , the azimuthal angular separation φ γ γ between the two photons, and the cosine of the polar angle θ * in the Collins-Soper frame of the diphoton system. The 2γ nnlo and sherpa pre-dictions show an improved agreement in shape with the data for the kinematic distributions with respect to the diphox + gamma2mc and resbos predictions, especially in the low m γ γ , low φ γ γ regions, which are the most sensitive to higher-order corrections. Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited. Funded by SCOAP 3 / License Version CC BY 4.0.

Appendix: Cross section tables
The numerical values of the cross sections for each bin of the diphoton invariant mass m γ γ , the diphoton transverse momentum p γ γ T , the azimuthal angle difference φ γ γ between the two photons, and the cosine of the polar angle θ * in the Collins-Soper frame of the diphoton pair are presented in Tables 4, 5, 6, and 7.