Measurement of the ratio of cross sections for inclusive isolated-photon production in pp collisions at s = 13 and 8 TeV with the ATLAS detector

: The ratio of the cross sections for inclusive isolated-photon production in pp collisions at centre-of-mass energies of 13 and 8 TeV is measured using the ATLAS detector at the LHC. The integrated luminosities of the 13 TeV and 8 TeV datasets are 3.2 fb (cid:0) 1 and 20.2 fb (cid:0) 1 , respectively. The ratio is measured as a function of the photon transverse energy in di(cid:11)erent regions of the photon pseudorapidity. The predictions from next-to-leading-order perturbative QCD calculations are compared with the measured ratio. The experimental systematic uncertainties as well as the uncertainties a(cid:11)ecting the predictions are evaluated taking into account the correlations between the two centre-of-mass energies, resulting in a reduction of up to a factor of 2 : 5 (5) in the experimental (theoretical) systematic uncertainties. The predictions based on several parameterisations of the proton parton distribution functions agree with the data within the reduced experimental and theoretical uncertainties. In addition, this ratio to that of the (cid:12)ducial cross sections for Z boson production at 13 and 8 TeV using the decay channels Z ! e + e (cid:0) and Z ! (cid:22) + (cid:22) (cid:0) is made and compared with the theoretical predictions. In this double ratio, a further reduction of the experimental uncertainty is obtained because the uncertainties arising from the luminosity measurement cancel out. The predictions describe the measurements of the double ratio within the theoretical and experimental uncertainties.

The production of prompt photons in proton-proton collisions, pp → γ + X, provides a means of testing perturbative QCD (pQCD) with a hard colourless probe.Since the dominant production mechanism in pp collisions at the LHC proceeds via the qg → qγ process, measurements of prompt-photon 1 production are sensitive to the gluon density in the proton [1,2].These measurements can also be used to tune Monte Carlo (MC) models to improve our understanding of prompt-photon production and aid those analyses for which events containing photons are an important background.
At leading order (LO) in pQCD, two processes contribute to prompt-photon production: the direct-photon process, in which the photon originates directly from the hard interaction, and the fragmentation-photon process, in which the photon is emitted in the fragmentation of a high transverse momentum (p T ) parton [3,4].
Measurements of prompt-photon production at a hadron collider necessitate an isolation requirement to reduce the large contribution of photons from neutral-hadron decays and the fragmentation component in which the emitted photon is close to a jet.The production of isolated photons in pp collisions has been measured previously by the ATLAS [5-9] and CMS [10,11] collaborations at centre-of-mass energies ( √ s) of 7, 8 and 13 TeV.
Comparisons of measurements of prompt-photon production and pQCD predictions are usually limited by the theoretical uncertainties associated with the missing higher-order terms in the perturbative expansion.
The measurements of inclusive isolated-photon cross sections performed by ATLAS at 13 TeV [9] and 8 TeV [8] were compared with the predictions of pQCD at next-to-leading order (NLO) [12,13].At both centre-of-mass energies, the uncertainties affecting the predictions are dominated by terms beyond NLO and are larger than those of experimental nature, preventing a more precise test of the theory.An avenue to reach a more stringent test is the inclusion of next-to-next-to-leading-order (NNLO) QCD corrections in the calculations [14].Another avenue is to make measurements of the ratio of cross sections for inclusive isolated-photon production at 13 and 8 TeV (R γ 13/8 ) and compare them with the predictions [15,16].The impact of the experimental systematic uncertainties and theoretical uncertainties on the ratio of the cross sections is reduced, allowing a more precise comparison between data and theory.This is achieved by accounting for inter-√ s correlations in the experimental systematic uncertainties affecting the measurements and in the uncertainties of the theory predictions.
A further reduction of the experimental uncertainty can be achieved by measuring a double ratio: the ratio of R γ 13/8 to the ratio of the fiducial cross sections for Z boson production at 13 TeV and 8 TeV (R Z  13/8 ≡ σ fid Z (13 TeV)/σ fid Z (8 TeV)) presented in Ref. [17].The measurements of the fiducial cross sections for Z boson production use the decay channels Z → e + e − and Z → µ + µ − .This observable, D γ/Z 13/8 ≡ R γ 13/8 /R Z 13/8 , can be viewed as the increase of the cross section for isolated-photon production as a function of √ s normalised to the increase for Z boson production as a function of √ s.Measuring D γ/Z 13/8 is beneficial because the uncertainties from the luminosity measurement cancel out, and D γ/Z 13/8 has only a slightly larger theory uncertainty than R γ 13/8 .This paper presents measurements of the ratio of cross sections for isolated-photon production in pp collisions at √ s = 13 TeV and 8 TeV with the ATLAS detector at the LHC.The phase-space region is given by the overlap of the ATLAS measurements at √ s = 13 and 8 TeV, defined by the photon transverse energy 2 (E γ T ) in the range E γ T > 125 GeV and the photon pseudorapidity (η γ ) in the region |η γ | < 2.37, excluding the region 1.37 < |η γ | < 1.56.The photon is isolated by requiring that the transverse energy inside a cone of size ∆R ≡ (∆η) 2 + (∆φ) 2 = 0.4 in the η-φ plane around the photon direction, E iso T , is smaller than E iso T,cut ≡ 4.2 • 10 −3 • E γ T + 4.8 GeV [8,9].Non-isolated prompt photons are not considered as signal.The measurements of the ratios are based on the ATLAS measurements at 13 TeV [9] and 8 TeV [8] and a detailed study of the correlations of the experimental systematic uncertainties between the two centre-of-mass energies is presented here.The measurement of the ratios is presented as a function of E γ T in different regions of η γ , namely |η γ | < 0.6, 0.6 < |η γ | < 1.37, 1.56 < |η γ | < 1.81 and 1.81 < |η γ | < 2.37.Next-to-leading-order pQCD predictions for the ratio are compared with the measurements.In addition, measurements of D γ/Z 13/8 are presented using the ATLAS results for R Z 13/8 [17]; the measurements are compared with available theory predictions.
The paper is organised as follows: the ATLAS detector is described in Section 2. The analysis strategy is summarised in Section 3. Fixed-order QCD predictions and their uncertainties are discussed in Section 4. Section 5 is devoted to the description of the experimental uncertainties.The results are reported in Section 6.A summary is given in Section 7. 2 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe.The x-axis points from the IP to the centre of the LHC ring, and the y-axis points upwards.Cylindrical coordinates (r, φ) are used in the transverse plane, φ being the azimuthal angle around the z-axis.The pseudorapidity is defined in terms of the polar angle θ as η = − ln tan(θ/2).

Analysis strategy
The measurements of ratios of cross sections presented in this paper are based on the measurements presented in previous ATLAS publications [8,9,17], where details of the analyses are given.The strategies followed for the measurement of the ratios and for the theoretical predictions are described below.

Analysis strategy for
The measurements of dσ/dE Photon candidates are reconstructed from clusters of energy deposited in the EM calorimeter.Candidates without a matching track or reconstructed conversion vertex in the inner detector are classified as unconverted photons, while those with a matching reconstructed conversion vertex or a matching track consistent with originating from a photon conversion are classified as converted photons [22].The photon identification is based primarily on shower shapes in the calorimeter [22].It uses information from the hadronic calorimeter, the lateral shower shape in the second layer of the EM calorimeter and the shower shapes in the finely segmented first EM calorimeter layer to ensure the compatibility of the measured shower profile with that originating from a single photon impacting the calorimeter.The photon energy measurement is made using calorimeter and, when available, tracking information.An energy calibration [23] is applied to the candidates to account for upstream energy loss and both lateral and longitudinal leakage.Events with at least one photon candidate with calibrated E γ T > 125 GeV and |η γ | < 2.37 excluding the region 1.37 < |η γ | < 1.56 are selected.The isolation transverse energy E iso T is corrected for leakage of the photon energy into the isolation cone and the estimated contributions from the underlying event (UE) and additional inelastic pp interactions (pile-up).The latter two corrections are computed simultaneously on an event-by-event basis using the jet-area method [24,25].After these corrections, isolated photons are selected by requiring E iso T to be lower than E iso T,cut .A small background contribution still remains after imposing the photon identification and isolation requirements and is subtracted using a data-driven method based on background control regions [8,9].The selected samples of events are used to unfold the distribution in E γ T for each |η γ | region to a phase-space region close to that used for event selection.The phase-space region at particle level uses particles with a decay length cτ > 10 mm; these particles are referred to as "stable".The particle-level isolation requirement for the photon is built by summing the transverse energy of all stable particles, except for muons and neutrinos, in a cone of size ∆R = 0.4 around the photon direction after the contribution from the UE is subtracted; the same subtraction procedure and isolation requirement used on data are applied at the particle level.
An important part of this analysis is the evaluation of the experimental systematic uncertainties in the ratio of the cross sections at 13 and 8 TeV taking into account correlations.This study is described in Section 5. Given the dominance of the systematic uncertainty arising from the photon energy scale when measuring the cross sections, it is necessary to carefully study this source of uncertainty.This source of systematic uncertainty is decomposed into independent components [23] and the treatment of the correlations of these components between the measurements at 13 and 8 TeV results in a reduction of the systematic uncertainty of the ratio.
The measurements of the ratio of cross sections are compared with NLO pQCD predictions for which a proper evaluation of the theoretical uncertainties is also of importance.The theoretical uncertainties in the predictions for the cross sections are O(10-15%) for both centre-of-mass energies and are dominated by contributions from terms beyond NLO.These uncertainties are much larger than those of experimental nature and limit how precisely the predictions can be tested.The study of the theoretical uncertainties in the ratio is described in Section 4. As is the case for the experimental systematic uncertainties, it is imperative that for each source of theoretical uncertainty the degree of correlation between the two centre-of-mass energies is taken into account.As a result, the theoretical uncertainty is reduced in the ratio, thus allowing a more stringent test of the predictions.

Analysis strategy for
The measurement of the double ratio D γ/Z 13/8 is based on the measurement of R γ 13/8 described above as well as on the measurement of R Z 13/8 .It should be noted that R γ 13/8 is measured as a function of E γ T in different ranges of η γ , while R Z 13/8 is a single number.The measurement of R Z 13/8 used here is the one reported in Ref. [17].The fiducial cross section at a given √ s, σ fid Z ( √ s), is defined as the production cross section of a Z boson times the branching ratio of the decay into a lepton pair of flavour + − = e + e − or µ + µ − within the following phase space: the lepton transverse momentum p T > 25 GeV, the lepton pseudorapidity |η | < 2.5 and the dilepton invariant mass 66 < m < 116 GeV.The measurement at √ s = 13 TeV was performed in the aforementioned phase space while the measurement at √ s = 8 TeV was extrapolated to the same phase space as described in Ref. [17].Measurements of the fiducial cross sections were made using the decay channels Z → e + e − and Z → µ + µ − , and combined for the final result.The measured R Z 13/8 is 1.537 ± 0.001 (stat.)± 0.010 (syst.)± 0.044 (lumi.)[17], where "stat."denotes the statistical uncertainty, "syst."denotes the systematic uncertainty and "lumi."denotes the uncertainty due to the ratio of the integrated luminosities.The evaluation of the systematic uncertainty in the ratio takes into account correlations of systematic uncertainties across channels and √ s as described in Ref. [17].
The predictions for D γ/Z 13/8 are obtained from NLO pQCD calculations for R γ 13/8 [12,13] and NNLO pQCD calculations for R Z  13/8 [26,27].The evaluation of the uncertainties affecting the predictions for D γ/Z 13/8 requires considerations that account for the correlations arising from the parton distribution functions (PDFs) and the strong coupling constant, α s (m Z ).
4 Fixed-order QCD predictions The theoretical predictions for the ratios of cross sections are obtained using fixed-order QCD calculations.Details of the generators and of the estimations of the theoretical uncertainties are given below, especially emphasising the correlations between the two centre-of-mass energies.The theoretical predictions for R γ 13/8 presented here are based on NLO QCD calculations computed using the program Jetphox 1.3.1_2[12,13].This program includes a full NLO QCD treatment of both the directand fragmentation-photon contributions to the cross section for the pp → γ + X reaction.The number of quark flavours is set to five.The renormalisation (µ R ), factorisation (µ F ) and fragmentation (µ f ) scales are chosen to be µ R = µ F = µ f = E γ T .The calculations are performed using various parameterisations of the proton PDFs and the BFG set II of parton-to-photon fragmentation functions at NLO [28].The nominal calculation is based on the MMHT2014 PDF set [29].Predictions are also obtained with other PDFs, namely CT14 [30], HERAPDF2.0[31], NNPDF3.0 [32] and ABMP16 [33].For MMHT2014, CT14, HERAPDF2.0and NNPDF3.0parameterisations of the PDFs, the sets determined at NLO are used.For ABMP16, the set at NNLO is used.The strong coupling constant α s (m Z ) is set to the value assumed in the fit to determine the PDFs; as an example, in the case of MMHT2014 PDFs, α s (m Z ) is set to the value 0.120.
The calculations are performed using a parton-level isolation criterion for the photon, which requires a total transverse energy of the partons inside a cone of radius R = 0.4 around the photon direction below E iso T,cut .The predictions from Jetphox are at parton level, 3 while the measurements are at particle level.Corrections for the non-perturbative (NP) effects of hadronisation and the UE are estimated using samples from Pythia 8.186 [34] as described below.First, a correction factor (C NP √ s ) is derived for the isolated-photon cross section at each centre-of-mass energy as the ratio of the cross section at particle level for a Pythia sample with UE effects to the Pythia cross section at parton level without UE effects.Second, the ratio of the correction factor for √ s = 13 TeV to that for √ s = 8 TeV, C NP R = C NP 13 /C NP 8 , is evaluated.The ratio of correction factors is obtained using the ATLAS set of tuned parameters A14 [35] with the LO NNPDF2.3PDF set [36].The ratio of correction factors for non-perturbative effects applied to the ratio predictions from Jetphox is C NP R = 0.9964 ± 0.0020.The following sources of uncertainty in the theoretical predictions are considered: • The uncertainty in the NLO QCD predictions due to terms beyond NLO is estimated by repeating the calculations using values of µ R , µ F and µ f scaled by the factors 0.5 and 2. The three scales are either varied simultaneously or individually; in addition, configurations in which one scale is fixed and the other two are varied simultaneously are also considered.In all cases, the condition 0.5 ≤ µ A /µ B ≤ 2 is imposed, where A, B = R, F, f.The final uncertainty is taken as the largest deviation from the nominal value among the 14 possible variations.
• The uncertainty in the NLO QCD predictions related to the proton PDFs is estimated by repeating the calculations using the 50 additional sets from the MMHT2014 error analysis.
• The uncertainty in the NLO QCD predictions related to the value of α s (m Z ) is estimated by repeating the calculations using two additional sets of proton PDFs from the MMHT2014 analysis for which different values of α s (m Z ) were assumed in the fits, namely α s (m Z ) = 0.118 and 0.122 [37].
• The impact of the beam energy uncertainty is estimated by repeating the calculations with √ s varied by its uncertainty of 0.1% [38].
• The uncertainty in the corrections for non-perturbative effects is estimated by comparing the results of using variations of the A14 tune in which the parameter settings related to the modelling of the UE are changed [35].
For the individual differential cross sections and for both centre-of-mass energies, the dominant theoretical uncertainty arises from the estimate of contributions from terms beyond NLO [8, 9].
The predictions for R γ 13/8 are obtained by calculating the ratio of the individual differential cross sections at each centre-of-mass energy.To estimate the theoretical uncertainty in R γ 13/8 , the correlation between the two centre-of-mass energies for each source listed above needs to be considered.The uncertainties due to the PDFs, α s (m Z ), beam energy and non-perturbative effects are fully correlated between the two centre-of-mass energies.The relative uncertainties in R γ 13/8 due to the uncertainties in α s (m Z ), the PDFs and the beam energy exhibit a significant degree of cancellation with respect to the individual predictions.However, for the scale uncertainties, the correlation is a priori unknown.In the standard approach, varying the scales coherently or incoherently at both centre-of-mass energies leads to very different theoretical uncertainties: • In the coherent case, there are large cancellations in the uncertainties in the predictions for R A second approach is also investigated, which is free from ambiguity in the correlation.It consists of considering the difference between the LO and NLO predictions for R γ 13/8 .The LO predictions are obtained with Jetphox using the same parameter settings and PDF set as the baseline NLO predictions.The LO and NLO predictions for R γ 13/8 are compared and the differences are up to 3.5%, which are similar to the estimates based on the standard approach with coherent variations at the two centre-of-mass energies.Thus, the results of this second approach support the use of the standard approach with a coherent variation of the scales; an incoherent variation of the scales clearly leads to an overestimation of the theoretical uncertainty.The theoretical predictions for D γ/Z 13/8 presented here are based on NNLO QCD calculations for the predictions of R Z 13/8 computed using the program Dyturbo, which is an optimised version of the Dynnlo program [26,27], and NLO QCD calculations for the predictions of R γ 13/8 using Jetphox with the procedure described in Section 4.1.
The calculations using Dyturbo are based on sets of PDFs extracted using NNLO QCD fits, namely MMHT2014nnlo, CT14nnlo, HERAPDF2.0nnloand NNPDF3.0nnlo.The strong coupling constant α s (m Z ) is set to the value assumed in the fit to determine the PDFs.In the case of MMHT2014nnlo PDFs, α s (m Z ) is set to the value 0.118.
For consistency, and to properly take into account the correlations in the PDF uncertainties, the calculations of Jetphox for R γ 13/8 are repeated using the NNLO PDF sets mentioned above.It is consistent to use NLO matrix elements convolved with PDF sets determined at NNLO.The resulting predictions include partially NNLO corrections and, therefore, are understood to still have NLO accuracy.For these additional calculations, the same parameter settings for the number of flavours, scales and fragmentation functions mentioned in Section 4.1 are used.The change in the predictions for R T for different η γ regions arising from the scale variations (shaded area), the value of α s (dashed lines), the PDF (dotted lines) and the beam energy (dot-dashed lines).The total theoretical uncertainty is shown as the solid line.
The sources of uncertainty in the theoretical predictions based on MMHT2014nnlo are the same as those described in Section 4.1.The uncertainty related to the beam energy is neglected, due to the small size of its effect on R γ 13/8 .The uncertainties in the prediction of R Z 13/8 due to the scale variations, the PDFs and α s (m Z ) are +0.02−0.3 %, +0.9 −0.8 % and −0.03 −0.3 %, respectively.For the predictions of D γ/Z 13/8 , the uncertainties have been estimated as follows: • The scale variations are considered uncorrelated between Z boson production and isolated-photon production since they are different processes.
• The PDF uncertainties are considered fully correlated between Z boson production and isolatedphoton production.
• The α s (m Z ) uncertainties are considered fully correlated between Z boson production and isolatedphoton production.The uncertainty in the predictions due to that in α s (m Z ) is estimated by using PDF sets in which α s (m Z ) was fixed at 0.116 or 0.120.

Experimental uncertainties
The sources of systematic uncertainties that affect the measurements of the photon differential cross sections at √ s = 8 and 13 TeV are detailed in Refs.
[8] and [9], respectively.A proper estimation of the systematic uncertainties in this measurement of cross-section ratios requires taking into account inter-√ s correlations for each source of systematic uncertainty.Assuming no correlation provides a conservative estimate and full correlation is used only when justified.The estimation of the systematic uncertainties in the ratio has to take into account the changes in the data-taking conditions as well as changes in the detector conditions.The measurements at √ s = 8 (13) TeV are based on data taken when the LHC operated with a bunch spacing of 50 (25) ns.During the data-taking period at √ s = 8 (13) TeV there were on average 20.7 (13.5) proton-proton interactions per bunch crossing.Furthermore, the addition of the silicon-pixel insertable B-layer leads to extra material upstream of the calorimeters for data-taking at √ s = 13 TeV.The procedures used to account for the impact of each source of systematic uncertainty on the ratio R γ 13/8 are described below.

Photon energy scale
The systematic uncertainties associated with the photon energy scale and resolution represent the dominant experimental uncertainties in the measurements of the differential cross sections for inclusive isolatedphoton production at both centre-of-mass energies.The uncertainty arising from the photon energy scale (γES) in R γ 13/8 is estimated by decomposing it into uncorrelated sources for both the 8 TeV and 13 TeV measurements.A total of 22 individual components [23] influencing the energy scale and resolution of the photon are considered.Twenty of these components are common to both centre-of-mass energies.These components include the uncertainties in: the overall energy scale adjustment using Z → e + e − events; the non-linearity of the energy measurement at the cell level; the relative calibration of the different calorimeter layers; the amount of material in front of the calorimeter; the modelling of the reconstruction of photon conversions; the modelling of the lateral shower shape; the modelling of the sampling term; 4 and the measurement of the constant term in Z boson decays.The uncertainties depend on E γ T as well as on |η γ | and are larger in the region 1.56 < |η γ | < 1.81 due to the presence of more material than in other |η γ | regions.The remaining two components are specific to the 13 TeV measurement and take into account the differences in the configuration of the ATLAS detector between 2012 and 2015, namely changes in the LAr temperature, in the stability of the layer intercalibration and in the material in front of the calorimeters between Run 1 and Run 2 [39].
The procedure used to estimate the systematic uncertainty in R γ 13/8 is as follows: all the uncertainty components described above are taken as fully correlated except for the uncertainty in the overall energy scale adjustment using Z → e + e − events, which for 2015 includes the effects of the changes in the configuration of the ATLAS detector mentioned above, and the uncertainties specific to the 13 TeV measurement.Calibration differences due to a change of optimal filtering coefficients and LAr timing samples between Run 1 and Run 2 are considered as a source of uncertainty in R γ 13/8 .The uncertainties in the photon energy scale due to pile-up are small enough compared to other uncertainties that the specific treatment of the correlation does not impact the results.The uncertainties due the photon energy resolution are treated as uncorrelated between √ s = 13 TeV and 8 TeV since they include the effects of pile-up, which was different in the 2012 and 2015 data-taking periods.
The relative uncertainty due to the correlated components of the photon energy scale in R γ 13/8 as a function of E γ T is shown in Figure 2 for each region in η γ .For illustration purposes, the result of estimating that part of the systematic uncertainty assuming no correlation is also shown in this figure: the results obtained using the complete correlation model exhibit a large reduction in comparison with those in which the correlations are ignored.This demonstrates that a proper treatment of the inter-√ s correlations in the systematic uncertainties associated with the photon energy scale is important.In addition, the relative uncertainty due to the uncorrelated components of the photon energy scale and the components specific to 2015 is also shown in Figure 2.  13/8 as a function of E γ T for different η γ regions due to the γES correlated components (dashed lines).For comparison, the results of considering the components as uncorrelated are also shown (dotted lines) to illustrate the reduction in the size of the systematic uncertainty when the proper treatment is applied.The relative uncertainty due to the uncorrelated components of the photon energy scale and the components specific to 2015 is also shown (solid lines).

Other sources of experimental uncertainty
The other sources of experimental uncertainty affecting the measurements are treated as listed below.For several of these sources, the uncertainties in the measurements at √ s = 13 TeV and 8 TeV are treated conservatively as uncorrelated since their impact is small.
• Statistical uncertainties.The statistical uncertainties in both the data and the Monte Carlo simulations at √ s = 13 and 8 TeV are treated as uncorrelated.
• Luminosity uncertainty.The luminosity uncertainties associated to the measurements of the photon cross sections at √ s = 8 TeV and 13 TeV are dominated by effects that are uncorrelated between different centre-of-mass energies and data-taking periods.The resulting relative uncertainty in R γ 13/8 amounts to ±2.8%.• Trigger uncertainty.The uncertainties in the trigger efficiency are treated as uncorrelated for data at different √ s.Different trigger requirements were used during 2012 and 2015.In addition, during 2012 a three-level trigger system was used to select events while in 2015 a two-level system was employed.
• Photon-identification uncertainty.In both measurements, the photon identification is based primarily on shower shapes in the EM calorimeter.These uncertainties are treated as uncorrelated since different methods are used at √ s = 13 TeV [40] and 8 TeV [22] to estimate the uncertainties; in addition, the photon identification criteria are re-optimised for data taken at 13 TeV.
• Modelling of the photon isolation in Monte Carlo.In both measurements, the photon candidate is required to be isolated.The in-time (out-of-time) pile-up, which is due to additional pp collisions in the same (earlier or later than) bunch crossing as the event of interest, was different in 2012 and 2015 due to the different LHC conditions, namely the instantaneous luminosity and the bunch spacing.
For simulated events, data-driven corrections to E iso T are applied such that the peak position in the E iso T distribution coincides in data and simulation.These uncertainties are treated as uncorrelated since different methods are used at √ s = 13 and 8 TeV for the corrections and uncertainties.
• Choice of background control regions.The background subtraction is performed using a datadriven two-dimensional sideband technique based on background control regions.A plane is formed by the variable E iso T and a binary variable that encapsulates the photon identification ("tight" vs. "nontight").A photon candidate is classified as "non-tight" if it fails at least one of four requirements on the shower-shape variables computed from the energy deposits in the first layer of the EM calorimeter, but satisfies the tight requirement on the total lateral shower width in the first layer and all the other tight identification criteria in other layers [22].The plane is divided into four regions: region A for tight isolated photons, region B for tight non-isolated photons, region C for non-tight isolated photons and region D for non-tight non-isolated photons.The background control regions B, C and D are specified by lower and upper limits on E iso T as well as by the definition of "non-tight" photon candidates.Variations of the limits and alternative definitions of the "non-tight" condition are used to estimate the uncertainties due to the choice of background control regions.These uncertainties are treated as uncorrelated since, as mentioned above, the photon-identification requirements are re-optimised for data-taking at 13 TeV.
• Photon identification and isolation correlation in the background.In the background subtraction method described above, the photon isolation and identification variables are assumed to be uncorrelated for background events.Uncertainties due to this assumption are estimated by using validation regions, which are dominated by background.These uncertainties are treated as uncorrelated since, as mentioned above, the photon-identification requirements are re-optimised for data-taking at 13 TeV.
• Signal modelling.MC simulations of signal processes are used to estimate the signal leakage fractions in the background control regions and to compute the unfolding corrections.For both measurements, at √ s = 13 and 8 TeV, the Pythia [34] generator is used for the nominal results and the Sherpa [41] generator for studies of systematic uncertainties related to the model dependence.The uncertainty due to the mixture of direct and fragmentation processes in the simulations is estimated using the MC simulations of Pythia.These uncertainties are treated as uncorrelated since different methods and versions of the generators are used at √ s = 13 TeV and 8 TeV to estimate the uncertainties.For √ s = 13 (8) TeV, Pythia 8.186 with the A14 tune (Pythia 8.165 with the AU2 tune) and Sherpa 2.1.1 (Sherpa 1.4.0) with the CT10 tune are used.For the 8 TeV analysis the results of using the default admixture of direct and fragmentation contributions in Pythia are compared with those using an optimal admixture obtained by fitting the two components to the data; for the 13 TeV analysis the results of enhancing the fragmentation contribution by a factor of two or removing it completely are compared with those using the default admixture.
• QCD-cascade and hadronisation model dependence.These uncertainties are treated as uncorrelated since different versions and tunes of the Monte Carlo generators Pythia and Sherpa are used at √ s = 13 and 8 TeV.
• Pile-up uncertainties.The in-time and out-of-time pile-up in the 2012 and 2015 data-taking periods were different.Conservatively and given the fact that the impact is rather small, these uncertainties are treated as uncorrelated.

Total experimental uncertainties in
Using the prescription for the treatment of the correlations between the measurements described in the previous sections, the systematic uncertainties in R γ 13/8 are evaluated.Figure 3 shows the relative uncertainties in R γ 13/8 due to (i) the photon energy scale, which includes the correlated and uncorrelated contributions as well as the additional ones associated with 2015 data, (ii) the remaining sources of systematic uncertainty excluding that in the luminosity measurements and (iii) the sum in quadrature of the non-γES uncertainties and the uncertainty due to the luminosity determination.The uncertainty due to the photon energy scale increases as E γ T increases and is larger for the region 1.56 < |η γ | < 1.81 due to more material in front of the calorimeters than in the other regions.From Figure 3 it is concluded that the relative uncertainty in R γ 13/8 due to the photon energy scale is no longer the dominant uncertainty, except for E γ T > 300 GeV in the regions 0.6 < |η γ | < 1.37 and 1.56 < |η γ | < 1.81.The total relative experimental systematic uncertainty in R γ 13/8 is shown in Figure 4, as is its sum in quadrature with the relative statistical uncertainty.In all pseudo-rapidity regions, the systematic uncertainty is dominant compared to the statistical uncertainty up to E γ T ∼ 300 GeV, while the measurement becomes statistically limited for E γ T 600 GeV.

Total experimental uncertainties in
The total relative experimental uncertainty in D γ/Z 13/8 is obtained as follows: • The uncertainty in R 13/8 as a function of E γ T for different η γ regions due to different sources: γES uncertainties (solid lines), non-γES uncertainties excluding the luminosity uncertainty (dashed lines) and non-γES and luminosity uncertainties added in quadrature (dotted lines).
• The statistical (0.1%) and systematic (0.7%) uncertainties in R Z 13/8 are added in quadrature to the total uncertainty in R  uncertainty in R γ 13/8 , the effect of adding in quadrature such a contribution has a small impact.On the other hand, the luminosity uncertainty, which amounts to 2.8% for R 13/8 (solid lines) and in D γ/Z 13/8 (dashed lines) as functions of E γ T for different η γ regions.

Results
The measurements of the ratios of cross sections are presented and the main features exhibited by the data are described.The theoretical predictions are compared with the experimental results for both R  5.20 ± 0.75 ± 0.25 8.4 ± 1.3 ± 0.5 900-1100 9.9 ± 2.3 ± 0.5 7.9 ± 2.4 ± 0.5 1100-1500 13.9 ± 9.8 ± 0.8 The NLO QCD predictions based on the MMHT2014 PDFs are compared with the measured R γ 13/8 in Figures 6 and 7.Even though there is a tendency for the predictions to underestimate the data, the measurements and the theory are consistent within the uncertainties; in particular, the increase as E γ T increases and the dependence on η γ are reproduced by the predictions.To study in more detail the description of the measured R γ 13/8 by the NLO QCD predictions, the ratio of the predictions to the data is shown in Figures 6 and 7.In these figures, the predictions based on different PDFs, namely MMHT2014, CT14, NNPDF3.0,HERAPDF2.0and ABMP16 are included to ascertain the sensitivity of R are below 2% for most of the phase space considered and smaller than the experimental uncertainties.The experimental uncertainties in R γ 13/8 also benefit from a significant reduction since the systematic uncertainties partially cancel out, in particular those related to the photon energy scale, which is dominant in the measurement of dσ/dE γ T .The total systematic uncertainty in R γ 13/8 is below 5% for most of the phase space considered.Thus, the significant reduction of the experimental and theoretical uncertainties in R γ 13/8 allows a more stringent test of NLO QCD.The overall level of agreement between data and the NLO QCD predictions based on several parameterisations of the proton PDFs within these reduced uncertainties validates the description of the evolution of isolated-photon production in pp collisions with the centre-of-mass energy.13/8 by the theoretical predictions, the ratio of the predictions to the data is shown in Figures 8 and 9.In these figures, the predictions based on different PDFs, namely MMHT2014nnlo, CT14nnlo, NNPDF3.0nnlo and HERAPDF2.0nnloare included to estimate the sensitivity of D γ/Z 13/8 to the proton PDFs.The predictions generally agree with the measured D γ/Z 13/8 within the experimental and theoretical uncertainties for all PDFs considered within the measured range.13/8 (black lines).The ratios of the pQCD predictions based on different PDF sets to the measured D γ/Z 13/8 are also included.13/8 (black lines).The ratios of the pQCD predictions based on different PDF sets to the measured D γ/Z 13/8 are also included.

Summary and conclusions
The ratio of cross sections for inclusive isolated-photon production in pp collisions at √ s = 13 and 8 TeV (R γ 13/8 ) is measured using the ATLAS detector at the LHC.The integrated luminosities of the 13 TeV and 8 TeV datasets are 3.2 fb −1 and 20.2 fb −1 , respectively.The ratio of differential cross sections as a function of E γ T is measured in different regions of |η γ | for photons with 125 < E γ T < 1500 GeV and |η γ | < 2.37, excluding the region 1.37 < |η γ | < 1.56.In the estimation of the experimental systematic uncertainties for R γ 13/8 , the correlations between the measurements at the two centre-of-mass energies are taken into account.The systematic uncertainty arising from the photon energy scale, which is dominant for the individual cross sections, is reduced significantly in R γ 13/8 and no longer the dominant uncertainty.The total systematic uncertainty for R γ 13/8 is below 5% in most of the phase space of the measurement.The measurements can be useful for tuning models of prompt-photon production in pp collisions.
The predictions from NLO QCD calculations are compared with the measured R γ 13/8 .The theoretical uncertainties affecting these predictions are also evaluated taking into account the correlations between the two centre-of-mass energies, resulting in a significant reduction in the uncertainty of the predicted R γ 13/8 .The theoretical uncertainties in R γ 13/8 are below 2% for most of the phase space of the measurement, in contrast with those in the individual cross-section predictions, which have approximately 10-15% uncertainties.Thus, the comparison of the predictions with the measured R γ 13/8 represents a stringent test of the pQCD calculations.Within these reduced experimental and theoretical uncertainties, the NLO QCD predictions based on several parameterisations of the proton PDFs agree with the data.Even though there is a tendency of the predictions to underestimate the data, the measurements and the theory are consistent within the uncertainties.The level of agreement achieved validates the description of the evolution of isolated-photon production in pp collisions from √ s = 8 to 13 TeV.
A double ratio of cross sections is also measured: the ratio of R γ 13/8 to the ratio of the fiducial cross sections for Z boson production at 13 and 8 TeV (R Z 13/8 ).In D γ/Z 13/8 ≡ R γ 13/8 /R Z 13/8 , the uncertainty due to the luminosity cancels out at the expense of a small increase in the systematic uncertainty from all other sources, leading to a more precise measurement of the evolution of the inclusive-photon cross section with the centre-of-mass energy normalised to the evolution of the Z boson cross section.The theoretical prediction, based on NNLO (NLO) QCD calculations for Z boson (inclusive-photon) production, describes the measurements within the theoretical uncertainties and the reduced experimental uncertainties.ARRS and MIZŠ, Slovenia; DST/NRF, South Africa; MINECO, Spain; SRC and Wallenberg Foundation, Sweden; SERI, SNSF and Cantons of Bern and Geneva, Switzerland; MOST, Taiwan; TAEK, Turkey; STFC, United Kingdom; DOE and NSF, United States of America.In addition, individual groups and members have received support from BCKDF, CANARIE, CRC and Compute Canada, Canada; COST, ERC, ERDF, Horizon 2020, and Marie Skłodowska-Curie Actions, European Union; Investissements d' Avenir Labex and Idex, ANR, France; DFG and AvH Foundation, Germany; Herakleitos, Thales and Aristeia programmes co-financed by EU-ESF and the Greek NSRF, Greece; BSF-NSF and GIF, Israel; CERCA Programme Generalitat de Catalunya, Spain; The Royal Society and Leverhulme Trust, United Kingdom.
(13 TeV)  used in the measurement of R γ 13/8 are based on an integrated luminosity of 20.2 ± 0.4 fb −1 (3.16 ± 0.07 fb −1 ).The measurement of the ratio covers the range E γ T > 125 GeV and is performed separately in the four regions of η γ defined in Section 1.A summary of the analyses leading to the measurements of the differential cross sections for inclusive isolated-photon production at √ s = 13 and 8 TeV is given below.

4. 1
Theoretical predictions for R γ 13/8 particularly in the variation of µ R , which is O(10%) for the individual predictions and below 1% for R γ 13/8 .The envelope of the scale variations for R γ 13/8 shrinks in comparison with the envelopes for the individual predictions: from O(10%) for the individual predictions to below 2% for R γ 13/8 across most of the range in E γ T .• In the incoherent case, the envelope of the scale variations for R γ 13/8 is O(14%) in all regions of phase space.

Figure 1
Figure 1 shows an overview of the theoretical uncertainties in R γ 13/8 .The total relative uncertainty is below 2% (4%) at low (high) E γ T in all regions of |η γ |.The uncertainty due to the variation of the scales is dominant everywhere.At high E γ T for |η γ | < 0.6 and 0.6 < |η γ | < 1.37, the uncertainty due to the PDFs can be as large as the contribution from the scale variations.The NLO pQCD predictions of Jetphox for R γ 13/8 based on the MMHT2014 parameterisations of the proton PDFs are about 2 at E γ T = 125 GeV and increase as E γ T increases, to about 10 for E γ T = 1300 (1000) GeV for |η γ | < 0.6 (0.6 < |η γ | < 1.37).For 1.56 < |η γ | < 1.81 (1.81 < |η γ | < 2.37), the predicted R γ 13/8 increases from about 2 at E γ T = 125 GeV to around 10 (25) at E γ T = 600 GeV.The increase is greater for the forward regions than for the central regions.Predictions based on different parameterisations of the proton PDFs are compared.Those based on MMHT2014, NNPDF3.0 and CT14 are found to be similar in all η γ and E γ T regions.The predictions of R γ 13/8 based on HERAPDF2.0and ABMP16 show some differences from the predictions based on the other PDFs in some regions of phase space, especially at high E γ T (see Section 6).

4. 2
Theoretical predictions for D

γ 13 /Figure 1 :
Figure 1: Relative theoretical uncertainty in R γ13/8 as a function of E γ T for different η γ regions arising from the scale variations (shaded area), the value of α s (dashed lines), the PDF (dotted lines) and the beam energy (dot-dashed lines).The total theoretical uncertainty is shown as the solid line.

γT
≤ 650 (650) GeV and it rises to ≈ 4.5% (3.3%) for E γ T = 1300 (1000) GeV.In both η γ regions, the total uncertainty is mostly dominated by the variation of the scales.For |η γ | < 0.6 and E γ T 300 GeV, the uncertainties in the PDFs are dominant, and for 0.6 < |η γ | < 1.37 and E γ T 750 GeV, the contributions from the scale variations and the PDFs are equally large.In the region 1.56 < |η γ | < 1.81 (1.81 < |η γ | < 2.37), the total relative uncertainty is below 2% (3%) for 125 ≤ E γ T ≤ 350 (470) GeV and it rises to ≈ 3% (3.6%) for E γ T = 600 GeV.For 1.56 < |η γ | < 1.81, the uncertainty due to the variation of the scales is dominant, but for 1.81 < |η γ | < 2.37 and E γ T 550 GeV, the contributions from the scale variations and the PDFs are equally important.The theoretical predictions based on the MMHT2014nnlo parameterisations of the proton PDFs for D γ/Z 13/8 are about 1.4 at E γ T = 125 GeV and increase as E γ T increases, to 6-17 at the high end of the spectrum, depending on the η γ region.The increase is larger for the forward regions than for the central regions.Predictions based on different parameterisations of the proton PDFs are compared; those based on MMHT2014nnlo, NNPDF3.0nnlo and CT14nnlo are found to be similar in all η γ and E γ T regions.The predictions of D γ/Z 13/8 based on HERAPDF2.0nnloshow some differences from the predictions based on the other PDFs in some regions of phase space, especially at high E γ T (see Section 6).

Figure 2 :
Figure 2: Relative systematic uncertainty in R γ13/8 as a function of E γ T for different η γ regions due to the γES correlated components (dashed lines).For comparison, the results of considering the components as uncorrelated are also shown (dotted lines) to illustrate the reduction in the size of the systematic uncertainty when the proper treatment is applied.The relative uncertainty due to the uncorrelated components of the photon energy scale and the components specific to 2015 is also shown (solid lines).

γ 13 /Figure 3 :
Figure 3: Relative systematic uncertainty in R γ13/8 as a function of E γ T for different η γ regions due to different sources: γES uncertainties (solid lines), non-γES uncertainties excluding the luminosity uncertainty (dashed lines) and non-γES and luminosity uncertainties added in quadrature (dotted lines).

γ 13 /
8 (see Section 5.3).The systematic uncertainty in R Z 13/8 is dominated by the uncertainty in the lepton reconstruction and efficiency; the correlation between the small contribution due to the electron energy scale and the photon energy scale in D γ/Z 13/8 can be safely neglected.The relative total systematic uncertainty in the measured D γ/Z 13/8 as a function of E γ T is shown for each η γ region in Figure5.For comparison, the relative total systematic uncertainty in R γ 13/8 is also shown.Since the total systematic uncertainty in R Z 13/8 is at least a factor of three smaller than the total systematic

Figure 4 :
Figure 4: Total relative systematic uncertainty in R γ 13/8 as a function of E γ T for different η γ regions (turquoise coloured band) and the sum in quadrature of the total relative systematic and statistical uncertainties (solid line).
a function of E γ T in different regions of |η γ | is shown in Figures 6 and 7 and Table 1.The measured R γ 13/8 increases with E γ T from approximately 2 at Eγ T = 125 GeV to approximately 8-29 at the high end of the spectrum.In the forward regions the increase of R γ 13/8 with E γ T is larger than in the central regions.At a fixed value of E γ T , the measured ratio increases as |η γ | increases.

γ 13 /
8 to the proton PDFs.The predictions generally agree with the measured R γ 13/8 within the experimental and theoretical uncertainties for all PDFs considered within the measured range.The comparison of the NLO QCD predictions for dσ/dE γ T and the measured differential cross sections in the ATLAS analyses at 8 and 13 TeV is limited by the theoretical uncertainties, which are larger than those of experimental nature and dominated by the uncertainties due to the terms beyond NLO.The theoretical uncertainties in dσ/dE γ T are 10-15%; in contrast, the theoretical uncertainties for R γ 13/8

Figure 6 :Figure 7 :
Figure 6: The measured R γ 13/8 (dots) as a function of E γ T in different regions of |η γ |.The NLO QCD predictions based on the MMHT2014 PDFs (black lines) are also shown.The inner (outer) error bars represent the statistical (total)uncertainties.The shaded band represents the theoretical uncertainty in the predictions.For most of the points, the error bars are smaller than the marker size and, thus, not visible.The lower part of the figures shows the ratio of the NLO QCD predictions based on the MMHT2014 PDFs to the measured R γ 13/8 (black lines).The ratios of the NLO QCD predictions based on different PDF sets to the measured R γ 13/8 are also included.

Figure 8 :
Figure 8: The measured D γ/Z 13/8 (dots) as a function of E γ T in different regions of |η γ |.The pQCD predictions based on the MMHT2014nnlo PDFs (black lines) are also shown.The inner (outer) error bars represent the statistical (total) uncertainties.The shaded band represents the theoretical uncertainty in the predictions.For most of the points, the error bars are smaller than the marker size and, thus, not visible.The lower part of the figures shows the ratio of the pQCD predictions based on the MMHT2014nnlo PDFs to the measured D γ/Z13/8 (black lines).The ratios of the pQCD predictions based on different PDF sets to the measured D γ/Z 13/8 are also included.

Figure 9 :
Figure 9: The measured D γ/Z 13/8 (dots) as a function of E γ T in different regions of |η γ |.The pQCD predictions based on the MMHT2014nnlo PDFs (black lines) are also shown.The inner (outer) error bars represent the statistical (total) uncertainties.The shaded band represents the theoretical uncertainty in the predictions.For most of the points, the error bars are smaller than the marker size and, thus, not visible.The lower part of the figures shows the ratio of the pQCD predictions based on the MMHT2014nnlo PDFs to the measured D γ/Z13/8 (black lines).The ratios of the pQCD predictions based on different PDF sets to the measured D γ/Z 13/8 are also included.

Table 2 :
The measured D γ/Z 13/8 as a function of E γ T together with the statistical and total systematic uncertainty in different regions of |η γ |.