Measurement of the ratio of the inclusive 3-jet cross section to the inclusive 2-jet cross section in pp collisions at sqrt(s) = 7 TeV and first determination of the strong coupling constant in the TeV range

A measurement is presented of the ratio of the inclusive 3-jet cross section to the inclusive 2-jet cross section as a function of the average transverse momentum, 〈pT1,2〉, of the two leading jets in the event. The data sample was collected during 2011 at a proton-proton centre-of-mass energy of 7 TeV with the CMS detector at the LHC, corresponding to an integrated luminosity of 5.0 fb−1. The strong coupling constant at the scale of the Z boson mass is determined to be αS(MZ) = 0.1148± 0.0014 (exp.)± 0.0018 (PDF)+0.0050 −0.0000 (scale), by comparing the ratio in the range 0.42 < 〈pT1,2〉 < 1.39 TeV to the predictions of perturbative QCD at next-to-leading order. This is the first determination of αS(MZ) from measurements at momentum scales beyond 0.6 TeV. The predicted ratio depends only indirectly on the evolution of the parton distribution functions of the proton such that this measurement also serves as a test of the evolution of the strong coupling constant beyond 0.42 TeV. No deviation from the expected behaviour is observed. Submitted to the European Physical Journal C c © 2013 CERN for the benefit of the CMS Collaboration. CC-BY-3.0 license ∗See Appendix A for the list of collaboration members ar X iv :s ub m it/ 07 06 04 1 [ he pex ] 2 8 A pr 2 01 3


Introduction
As a consequence of the non-Abelian nature of quantum chromodynamics (QCD), the renormalization group equation (RGE) [1][2][3] predicts that the strong force becomes weaker at short distances corresponding to large momentum transfers, a property of QCD referred to as asymptotic freedom. The strength of the strong force, α S (Q), at a given distance or momentum scale Q is not predicted and has to be extracted from experiment. Measurements at different Q can then be compared for consistency with QCD via the RGE, which precisely describes the evolution of α S (µ r ), where the renormalization scale µ r is identified with Q. By convention, the consistency is tested by evolving all values of α S (Q) to the common scale µ r = Q = M Z , i.e. the precisely known mass of the Z boson. The current world average value is α S (M Z ) = 0.1184 ± 0.0007 [4].
Measurements of the running of α S (Q) provide a stringent test of QCD. Previous collider experiments at LEP and HERA have established the validity of the RGE up to momentum transfers Q of 208 GeV [4]. A recent publication by the D0 Collaboration extends this range up to 400 GeV [5]. The determination of α S (Q) from jet cross sections as in [6] or [7] depends directly on parton distribution functions (PDFs) that have been evolved from small to very high momentum scales via the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations [8][9][10], which assume the validity of the RGE. This dependence on the evolution of the PDFs can be reduced by investigating cross section ratios. The ratio R 32 of the inclusive 3-jet cross section to the inclusive 2-jet cross section is proportional to α S (Q) where Q is defined as the average transverse momentum of the two jets leading in p T , Many theoretical systematic uncertainties related to the choice of the renormalization and factorization scales, µ r and µ f , or to nonperturbative effects are reduced in the cross section ratio. In addition, experimental uncertainties such as those due to the jet energy scale largely cancel in the measurement of R 32 . The uncertainty on the integrated luminosity measurement cancels completely. The Compact Muon Solenoid (CMS) Collaboration has previously measured R 32 [11], and the predictions of various Monte Carlo (MC) event generators were found to be in general agreement with the measurement.
This measurement is performed using a sample of multijet events, collected during 2011 by the CMS experiment at the Large Hadron Collider (LHC), corresponding to an integrated luminosity of 5.0 fb −1 of pp collisions at a centre-of-mass energy √ s = 7 TeV. The transverse momentum p T and the rapidity y of a jet with energy E and momentum p = (p x , p y , p z ) (where p z is the momentum component along the direction of the anticlockwise proton beam) are defined as p T = , respectively. Jets are reconstructed using the infrared-and collinear-safe anti-k T clustering algorithm [12,13] with a size parameter of 0.7. This measurement uses jets with p T > 150 GeV and |y| < 2.5.

Event selection and reconstruction
is instrumented with various layers of particle detection systems. 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)], and θ is the polar angle with respect to the z axis. Surrounding the trackers are a lead tungstate crystal electromagnetic calorimeter (ECAL) with a preshower detector in the endcaps, and a brass and scintillator hadron calorimeter (HCAL), covering the region |η| < 3. In addition to the barrel and endcap detectors, CMS has extensive forward calorimetry which extends the coverage to |η| = 5. The steel flux return yoke outside the solenoid is instrumented with gas-ionization detectors used to identify and reconstruct muons. A more detailed description of the CMS detector can be found in [14].

Event selection and reconstruction
The CMS detector records events using a two-level trigger system consisting of a hardwarebased level-1 (L1) trigger and a software-based high level trigger (HLT). In this study, singlejet triggers that reconstruct jets from calorimeter energy deposits at L1 and HLT are used to select events based on three HLT p T thresholds, 190, 240, and 370 GeV. All except the highestthreshold trigger were prescaled during the 2011 run. The corresponding integrated luminosity L for each of the three samples is shown in Table 1. The efficiency of each of the triggers is estimated using lower-p T -threshold triggers. These three jet trigger thresholds ensure 100% trigger efficiency in the three jet samples for p T1,2 > 215, 269, and 409 GeV.  Each event is required to have at least one offline-reconstructed vertex [15] along the beam line that is within 24 cm of the nominal interaction point. The four-vectors of particle candidates reconstructed by the CMS global event reconstruction algorithm (also called particle-flow event reconstruction [16]) are used as input to the jet-clustering algorithm. The clustering is performed by the FASTJET package [13] using four-momentum summation. The global event reconstruction algorithm reconstructs and identifies each particle with an optimized combination of subdetector information. The energy of photons is obtained directly from the ECAL measurements after being corrected for zero-suppression effects. The energy of electrons is determined from a combination of the track momentum at the main interaction vertex, the corresponding ECAL cluster energy, and the energy sum of all bremsstrahlung photons attached to the track. The energy of muons is derived from the corresponding track momentum. The energy of charged hadrons is determined from a combination of the track momentum and the corresponding ECAL and HCAL energies, corrected for zero-suppression effects, and calibrated for the nonlinear response of the calorimeters. Finally, the energy of neutral hadrons is obtained from the corresponding calibrated ECAL and HCAL energies.
Jet energy corrections [17] are derived using simulated events, generated by PYTHIA 6.4.22 [18] and processed through the CMS detector simulation based on GEANT4 [19], and in situ measurements with dijet, photon+jet, and Z+jet events. An offset correction is applied to take into account the extra energy clustered into jets from additional proton-proton interactions within the same or neighbouring bunch crossings (in-time and out-of-time pileup) [17]. Pileup effects are important only for low-p T jets and become negligible for jets with p T > 200 GeV. The current measurement is therefore largely insensitive to pileup effects. The jet energy corrections, which depend on the η and p T of the jet, are applied to the jet four-momentum vector as a multiplicative factor. The multiplicative factor is in general smaller than 1.2, approximately uniform in η, with typical values of 1.1 for jets having p T = 100 GeV and decreasing to 1.0 for higher values of p T .
To suppress nonphysical jets, i.e. jets resulting from noise in the ECAL and/or HCAL calorimeters, tight identification criteria are applied: each jet should contain at least two particles, one of which is a charged hadron, and the jet energy fraction carried by neutral hadrons, photons, muons, and electrons should be less than 90%. These criteria have an efficiency greater than 99% for genuine jets.
The selection of multijet events requires two or more jets with transverse momentum greater than 150 GeV and |y| < 2.5. The final sample is extracted by rejecting events if either or both of the leading jets in p T have |y| > 2.5.

Measurement of R 32 and comparison with theoretical predictions
The measured ratio R 32 as a function of p T1,2 is the ratio of the number of selected inclusive 3-jet events to the number of selected inclusive 2-jet events in each p T1,2 bin. The ratio R 32 is corrected for detector smearing effects and unfolded to stable-particle level. The unfolding method is the iterative Bayesian method [20], as implemented in the ROOUNFOLD software package [21]. Unfolding uses a response matrix that maps the true distribution onto the measured one. The response matrix is derived from a simulation, which uses as input the true R 32 distribution from PYTHIA6 tune Z2 and introduces the smearing effects by taking into account the p T1,2 resolution [17]. After unfolding R 32 to stable-particle level, the final statistical uncertainties include the correlation among the various p T1,2 bins.
Two main sources of systematic uncertainties on R 32 are considered. The first is due to the jet energy scale (JES) and the second due to the unfolding.
The JES uncertainty has been estimated to be 2.0-2.5% for particle-flow jets [22], depending on the jet p T and η. All mutually uncorrelated JES uncertainty sources are considered following the procedure described in Ref. [23]. The total systematic uncertainty on R 32 due to the JES uncertainty is 1.2%.
The unfolding method takes into account three different mutually uncorrelated uncertainty sources. The first arises from insufficient knowledge of the simulated inclusive 3-jet and 2-jet p T1,2 spectra, which are used to construct the simulated ratio R 32 used in the unfolding. The uncertainty is estimated by varying the 3-jet and 2-jet spectra slopes by ±10%. This is a conservative estimate and is motivated by the observed difference in the 3-jet and 2-jet spectra slopes between simulations using the event generators PYTHIA6 tune Z2 and HERWIG++ [24] with its default tune in version 2.3. Simulations of 3-jet and 2-jet spectra using the MADGRAPH [25,26] event generator are in agreement with those of PYTHIA6 . The second uncertainty arises from the insufficient knowledge of the p T1,2 resolution and is estimated by varying it by ±10%. This variation is motivated by the observed difference between data and simulation in the jet energy resolution [17]. Finally, the third uncertainty arises from non-Gaussian components in the p T1,2 resolution and is estimated by adding non-Gaussian tails to the simulation. The overall systematic uncertainty on R 32 due to unfolding is less than 1%. A potential bias originat-ing from the unfolding technique is studied by comparing the unfolding result of the Bayesian method with that of the singular-value decomposition (SVD) method [27]. The bias is found to be negligible.
The theoretical calculation of the ratio R 32 is based on next-to-leading-order (NLO) perturbative QCD (pQCD) calculations multiplied by a nonperturbative factor, which corrects for multiparton interactions (MPI) and hadronization effects. The NLO calculations assume N f = 5 massless quark flavours and are based on the parton-level generator NLOJET++ [28,29]. The computations with NLOJET++ are performed within the FASTNLO framework [30,31] using the following four PDF sets: NNPDF2.1 [32,33], ABM11 [34], MSTW2008 [35,36], and CT10 [37,38]. In each case the PDF version employing next-to-next-to-leading order (NNLO) evolution code is chosen, and for comparisons the respective default values of α S (M Z ), which are 0.119, 0.1134, 0.1171, and 0.118, are used. NNPDF2.1, MSTW2008, and CT10 utilize a variableflavour-number scheme with the maximal number of flavours N f ,max equal to 6, 5, and 5 respectively, while the ABM11 PDF set was developed in a fixed-flavour-number scheme with N f = 5. The renormalization and factorization scales are set to the average transverse momentum p T1,2 .
The nonperturbative effects are estimated using the PYTHIA6 tune Z2 and HERWIG++ tune 2.3 event generators. The chosen MC models feature different descriptions of the phenomena and are representative of the possible values of the nonperturbative corrections. The nonperturbative correction (NPC) factor is defined as the ratio of R 32 predicted with the nominal generator settings to that obtained with the MPI and hadronization switched off. This factor is calculated considering the average of the two MC generators and has typical values of ≈1.02 for p T1,2 = 250 GeV, decreasing to 1.0 for higher p T1,2 . The uncertainty is considered to be half of the difference between the NPC values obtained using the two MC generators and amounts to ≈0.1%, leading to a negligible influence on the final result.
Scale uncertainties have a very similar behaviour for all PDF sets and dominate the region up to p T1,2 ≈ 400 GeV. A comparison to jets with a size parameter of R = 0.5 reveals consistent results but with larger scale uncertainties.
The PDF uncertainties are different for each individual PDF set. The CT10 set exhibits the largest PDF uncertainties, which are of the order of 2% at p T1,2 = 400 GeV, increasing to 2.5% in the 1 TeV region. For the NNPDF2.1 set PDF uncertainties are of the order of 1.5% at 400 GeV, increasing to 2.3% at 1 TeV. Finally, for the MSTW2008 and ABM11 sets PDF uncertainties are of the order of 1% throughout the range of this measurement.
The comparison of data with the predictions of pQCD in Fig. 1 demonstrates that the NLO calculations using the NNPDF2.1, MSTW2008, and CT10 PDF sets are in agreement with the measured ratio R 32 throughout the range of this measurement. The NLO result employing the ABM11 PDF set underestimates R 32 , especially for p T1,2 < 600 GeV.

Determination of α S (M Z )
The measurement of the ratio R 32 is used for the determination of the strong coupling constant α S (M Z ). The value of α S (M Z ) is determined by minimizing the χ 2 between the experimental measurement and the theoretical predictions. The χ 2 is defined as where M is the vector of the differences between the data (R i 32 ) and the theoretical values (T i 32 ) in each bin i, and C is the covariance matrix including all experimental (statistical, JES, and unfolding) uncertainties. C is defined as where Cov Stat is the statistical covariance matrix that accounts for the correlations due to unfolding, and Cov JES Sources , Cov Unfolding Sources are the covariance matrices that account for the JES and unfolding systematic uncertainty sources, respectively. Each systematic uncertainty source for the JES and unfolding is treated as 100% correlated across the p T1,2 bins.
To avoid the region with large scale uncertainties close to the minimal jet p T requirements, visible in Fig. 1, α S (M Z ) is extracted only for p T1,2 > 420 GeV. The central result is obtained by minimizing the χ 2 (Eq. (2)) with respect to α S (M Z ) for the NNPDF2.1 PDF set, which is the only one that permits the propagation of PDF uncertainties to the fits for each value of α S (M Z ). The contribution of PDFs to the uncertainty of the measurement is evaluated by repeating the fit for each of the 100 PDF replicas of the NNPDF2.1 set at the relevant value for α S (M Z ). In this way 100 determinations of α S (M Z ) are obtained, whose distribution corresponds to the propagation of the underlying probability density from the PDFs to the fitted strong coupling. The PDF uncertainty of the measurement is then computed as the standard deviation of this distribution. A more detailed description of the method can be found in Ref. [39].
The uncertainties due to the renormalization and factorization scales are treated separately by varying the default choice of µ r = µ f = p T1,2 between p T1,2 /2 and 2 p T1,2 in six combinations as explained in Section 4. The χ 2 minimization with respect to α S (M Z ) is repeated for these six combinations. The contribution from the µ r , µ f scale variations to the uncertainty in the measurement is evaluated by considering the differences between the NNPDF2.1 α S (M Z ) central value and the highest and lowest values found in these six scale combinations. Out of all scale combinations the lowest α S (M Z ) value corresponds to the default scale choice of µ r = µ f = p T1,2 and the highest to the scale choice of µ r = µ f = p T1,2 /2. The frequent observation of asymmetric scale uncertainties with larger downward uncertainties in the case of NLO cross sections is transformed into a purely upward uncertainty for the ratio, as can be seen in Table 2. Table 2: The values of α S (M Z ) at the central scale and for the six scale factor combinations. A cross check on the impact of the top quark by imposing N f = 6 massless flavours in the NLO matrix elements revealed an increase by +0.0009 in the fitted value of α S (M Z ). Further effects, for example from the evolution of α S and the PDFs with five or six flavours, multijet production via fully hadronic decays in the reaction pp → tt + X, or an incomplete cancellation of electroweak corrections between numerator and denominator, are estimated to contribute each at a ±1% level to the theoretical uncertainty. These residual effects are taken into account by symmetrizing the scale uncertainty such that the largest deviation is adopted as the total symmetric theory uncertainty.
The determination of α S (M Z ), which is based on the NNPDF2.1 PDF set, is also in agreement with the results obtained using the MSTW2008 or CT10 PDF sets with χ 2 /N dof = 20.6/20 and 21.1/20, respectively. If PDF sets with NLO evolution are used instead the impact on the results of the fits to the ratio observable R 32 is negligible. This is in contrast to fits to cross sections, where NNLO PDF sets usually lead to smaller values of α S (M Z ) than NLO ones, and confirms the reduced dependence of R 32 on details of the PDF evolution.
In the case of the ABM11 PDFs the series in values of α S (M Z ) ends at 0.120, which is insufficient for a derivation of the complete shape of the χ 2 curve at minimum such that a fit value for It is observed that with ABM11 PDFs a higher value of α S (M Z ) is preferred. This is in accord with the fact that the ABM11 gluon density in the phase space relevant for this analysis is significantly smaller than that of all other PDF sets. Thus, the fit favors a larger α S (M Z ) value to compensate for this effect. In summary, the ABM11 PDF set does not describe the data as well as the alternative PDF sets, as shown in Fig. 1, which leads to an inferior fit quality and a less consistent result for the strong coupling.
To investigate the running of the strong coupling constant in more detail, the fitted region of 420-1390 GeV is split into three bins of p T1,2 and the fitting procedure is repeated in each of these bins. The three separate extractions of α S (M Z ) are reported in Table 3. The experimental uncertainties in the three obtained values are correlated. These α S (M Z ) determinations are then evolved back to the corresponding values α S (Q) using the 3-loop solution to the RGE from the NNPDF2.1 set. For each fit region the cross-section-weighted average of p T1,2 from the inclusive dijet calculation at NLO with NLOJET++ is chosen as the momentum scale Q and is computed to be Q = 474, 664, and 896 GeV, respectively. These values, derived again with the FASTNLO framework, are identical within about 1 GeV for different PDFs and vary at most by a few GeV when using inclusive 3-jet events.
To emphasize that theoretical uncertainties limit the achievable precision, Table 4 presents the decomposition of the total uncertainty for the three bins in p T1,2 into the experimental, PDF, and theory components. Figure 3 presents the strong coupling α S (Q) (solid line) and its total uncertainty (band) as evolved from the CMS determination, α S (M Z ) = 0.1148 ± 0.0055, using the 3-loop solution to the RGE from NNPDF2.1, as before. The extractions of α S (Q) in three separate ranges of Q as presented in Table 3 Figure 3: The strong coupling α S (Q) (solid line) and its total uncertainty (band) evolved from the CMS determination α S (M Z ) = 0.1148 ± 0.0055 using a 3-loop solution to the RGE as a function of the momentum transfer Q = p T1,2 . The extractions of α S (Q) in three separate ranges of Q as presented in Table 3

Summary
The ratio R 32 of the inclusive 3-jet cross section to the inclusive 2-jet cross section, for jets with p T > 150 GeV and |y| < 2.5, has been measured in the range 250 < p T1,2 < 1390 GeV for proton-proton collisions at a centre-of-mass energy of 7 TeV. The results have been compared with predictions of QCD at NLO obtained with various PDF sets. The NLO calculations using the NNPDF2.1, MSTW2008, and CT10 NNLO PDF sets are in agreement with the measured ratio R 32 throughout the range of p T1,2 studied. However, calculations using the ABM11 PDF sets underestimate R 32 for p T1,2 < 600 GeV.
Measurements of R 32 over the range 420 < p T1,2 < 1390 GeV have been used to determine the strong coupling constant α S at the scale of the Z boson mass. The final result is α S (M Z ) = 0.1148 ± 0.0014 (exp.) ± 0.0018 (PDF) ± 0.0050 (theory) = 0.1148 ± 0.0055, where experimental, PDF, and theory uncertainties have been added quadratically to give the total uncertainty. The result is in agreement with the world average value of α S (M Z ) = 0.1184 ± 0.0007 [4] and represents the first determination of the strong coupling constant from jet measurements with momenta of the order of 1 TeV. The dominant uncertainties are of theoretical origin and limit the currently achievable precision. The predicted ratio depends only indirectly on the evolution of the parton distribution functions of the proton and consequently this measurement also serves as a test of the evolution of α S (Q). No deviation from the expected behaviour is observed. and personnel of the Worldwide LHC Computing Grid for delivering so effectively the computing infrastructure essential to our analyses. Finally, we acknowledge the enduring support for the construction and operation of the LHC and the CMS detector provided by the following funding agencies: [41] H1 Collaboration, "Jet production in ep collisions at high Q 2 and determination of α S ", Eur. Phys. J. C 65 (2010) 363, doi:10.1140/epjc/s10052-009-1208-7, arXiv:0904.3870.