Measurement of three-jet production cross-sections in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$pp$$\end{document}pp collisions at 7 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {\ Te\,V}$$\end{document}TeV centre-of-mass energy using the ATLAS detector

*************Abstract Double-differential three-jet production cross-sections are measured in proton–proton collisions at a centre-of-mass energy of \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\mathrm \,TeV{}$$\end{document}s=7TeV using the ATLAS detector at the large hadron collider. The measurements are presented as a function of the three-jet mass \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(m_{jjj})$$\end{document}(mjjj), in bins of the sum of the absolute rapidity separations between the three leading jets \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\left| Y^{*}\right| )$$\end{document}(Y∗). Invariant masses extending up to 5 TeV are reached for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$8< \left| Y^{*}\right| < 10$$\end{document}8<Y∗<10. These measurements use a sample of data recorded using the ATLAS detector in 2011, which corresponds to an integrated luminosity of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4.51~\text{ fb }^{-1}$$\end{document}4.51fb-1. Jets are identified using the anti-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_{t}$$\end{document}kt algorithm with two different jet radius parameters, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R=0.4$$\end{document}R=0.4 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R=0.6$$\end{document}R=0.6. The dominant uncertainty in these measurements comes from the jet energy scale. Next-to-leading-order QCD calculations corrected to account for non-perturbative effects are compared to the measurements. Good agreement is found between the data and the theoretical predictions based on most of the available sets of parton distribution functions, over the full kinematic range, covering almost seven orders of magnitude in the measured cross-section values.


Introduction
Collimated jets of hadrons are a characteristic feature of highenergy particle interactions. In the theory of strong interactions, quantum chromodynamics (QCD), jets can be interpreted as the result of fragmentation of partons produced in a scattering process. In high-energy particle collisions two main phases can be distinguished. In the perturbative phase, partons with high-transverse momentum ( p T ) are produced in a hard-scattering process at a scale Q. This phase is described by a perturbative expansion in QCD. In the transition to the second (non-perturbative) phase, these partons e-mail: atlas.publications@cern.ch emit additional gluons and produce quark-antiquark pairs. The non-perturbative jet evolution is an interplay between the hadronisation process and the underlying event. The hadronisation process governs the transition from partons to hadrons and the underlying event represents initial-state radiation, multiple parton interactions and colour-reconnection effects [1]. In spite of these phenomena, the highly collimated sprays of particles, collectively identified as hadron jets, are observed in the final state. The effects of both hadronisation and the underlying event vary strongly with the jet radius parameter and are most pronounced at low p T . They are accounted for using phenomenological models that are tuned to the data.
The ATLAS Collaboration has measured the inclusive jet cross-sections at 7 TeV [2] and at 2.76 TeV [3] centre-ofmass energies in pp collisions for jets defined by the antik t algorithm [4] with two jet radius parameters, R = 0.4 and R = 0.6. Recent inclusive jet [5] and dijet [6] crosssection measurements at 7 TeV centre-of-mass energy in pp collisions have exploited improved jet energy calibration procedures [7] leading to smaller systematic uncertainties compared to those achieved in Refs. [2,3]. Similar measurements at 7 TeV centre-of-mass energy in pp collisions [8,9] have been carried out by the CMS Collaboration. These measurements test perturbative QCD (pQCD) at very short distances and have provided constraints on the gluon momentum distribution within protons at large momentum fraction. The impact of higher order effects on the inclusive jet crosssection ratios of anti-k t R = 0.5 and R = 0.7 jets has been studied in [10]. The inclusive three-jet to two-jet ratio [11] is used to determine the strong coupling constant. Theoretical predictions of the multi-jet cross-sections in pp collisions at 7 TeV centre-of-mass energy have been tested in Refs. [12,13].
Previous measurements of three-jet cross-sections in pp collisions were performed by the D∅ collaboration [14]. The measurements were compared to predictions, and agreement between data and theory was found within the uncertainties.
In this paper, measurements of double-differential threejet production cross-sections are presented as a function of the three-jet mass (m j j j ) and the sum of absolute rapidity separation between the three leading jets (|Y * |). The measurements are corrected for experimental effects and reported at the particle level. The three-jet mass distributions test the dynamics of the underlying 2 → 3 scattering process. The distributions are sensitive to both the transverse momentum ( p T ) spectra of the three leading jets and their angular correlations, since a massive three-jet system can be built either from highp T jets or from jets with large rapidity separation. Binning in |Y * | allows events with m j j j originating from these different regions of phase space to be separated.
The analysis presented in this paper tests the description of multi-jet events in next-to-leading-order (NLO) QCD and uses two different values of jet radius parameter, R = 0.4 and R = 0.6, since three-jet crosssections depend on the jet radius even at leading order (LO) in the perturbative expansion. The NLO QCD calculations corrected to account for non-perturbative effects are compared to the measured cross-sections. The measurements also provide constraints on the proton's parton distribution functions (PDFs) beyond those from inclusive and dijet crosssections, since they probe a different region of phase space in proton momentum fraction and squared momentum transfer (x, Q 2 ) and different combinations of initial-state partons.
The content of this paper is structured as follows. The ATLAS detector is briefly described in Sect. 2, followed by the definition of observables and description of Monte Carlo (MC) samples in Sects. 3 and 4, respectively. The trigger, data selection and jet calibration are presented in Sect. 5. Data unfolding and experimental uncertainties are described in Sects. 6 and 7. Section 8 describes the theoretical predictions for the measurements in this paper. The cross-section results are presented in Sect. 9 and the conclusions are given in Sect. 10.

The ATLAS experiment
The ATLAS detector is described in detail in Ref. [15]. 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 pointing along the beam axis. The x-axis points from the IP to the centre of the LHC ring, and the yaxis points upward. Cylindrical coordinates (r , φ) are used in the transverse plane, φ being the azimuthal angle around the beam pipe. The pseudorapidity is defined in terms of the polar angle θ as η = − ln tan(θ/2). The rapidity is defined in terms of the energy E and longitudinal to the beam pipe momentum p z as y = 1/2 ln ((E + p z )/(E − p z )). The transverse momentum p T is defined as the component of the momentum transverse to the beam pipe.
The inner detector (ID) is used to measure the momenta and trajectories of charged particles. The ID has full coverage in the azimuthal angle φ and over the pseudorapidity range |η| < 2.5. The ID is immersed in a 2 T magnetic field provided by a superconducting solenoid magnet. The main detector system used for this analysis is the calorimeter. The electromagnetic calorimeters use liquid argon (LAr) as the active detector medium. They employ accordion-shaped electrodes and lead absorbers, and are divided into one barrel (|η| < 1.475) and two end-cap components (1.375 < |η| < 3.2). The technology used for the hadronic calorimeters depends on η. In the barrel region (|η| < 1.7), the detector is made of scintillator tiles with steel absorbers. In the end-cap region (1.5 < |η| < 3.2), the detector uses LAr and copper. A forward calorimeter consisting of LAr and tungsten/copper absorbers has both electromagnetic and hadronic sections, and extends the coverage to |η| = 4.9.
The muon spectrometer has one barrel and two end-cap air-core toroid magnets. Three layers of precision tracking stations provide muon momentum measurements over the range |η| < 2.7. The ATLAS trigger system consists of three levels of event selection: a first level implemented using custom-made electronics, which selects events at a design rate of at most 75 kHz, followed by two successive software-based levels. The level-2 trigger uses fast online algorithms, and the final trigger stage, event filter (EF), uses reconstruction software with algorithms similar to the offline versions.

Cross-section definition
Jets are defined using the anti-k t algorithm as implemented in the FastJet [16] package, with two different values of the radius parameter: R = 0.4 and R = 0.6.
Events containing at least three jets within the rapidity range |y| < 3.0 with p T > 50 GeV are considered. The leading, subleading and sub-subleading jets are required to have p T > 150 GeV, p T > 100 GeV and p T > 50 GeV, respectively.
Three-jet double-differential cross-sections are measured as a function of the three-jet mass m j j j = ( p 1 + p 2 + p 3 ) 2 and the summed absolute rapidity separation of the three leading jets Y * = |y 1 − y 2 | + |y 2 − y 3 | + |y 1 − y 3 | , where p i (y i ) are the four-momenta (rapidities) of the three leading jets. The measurements are made in five ranges of |Y * | < 10, in equal steps of two. In each range of |Y * |, a lower limit on the three-jet mass is imposed to avoid the region of phase space affected by the jet p T cuts. The measurement starts at m j j j = 380 GeV in the |Y * | < 2 bin, increasing to 1180 GeV for the 8 < |Y * | < 10 bin.
The three-jet mass distributions are corrected for detector effects, and the measured cross-sections are defined at the particle level. Here particle level refers to jets built using produced particles with a proper lifetime longer than 10 ps, including muons and neutrinos from decaying hadrons [17].

Monte Carlo samples
The default MC generator used to simulate events is Pythia 6 [18] with the Perugia 2011 tune [19] and the CTEQ5L PDFs [20]. Usually, "tune" refers to a set of model parameters, which provide an optimal description of high-energy particle collisions. Data from previous colliders (LEP, TEVA-TRON, etc), as well as early LHC data are included in the process of tuning the model parameters [19,21,22]. The Pythia 6 is a generator with LO 2 → 2 matrix element calculations, supplemented by leading-logarithmic calculations of parton showers ordered in p T . A simulation of the underlying event, including multiple parton interactions, is also included. The Lund string model [23,24] is used to simulate the fragmentation process. The signal reconstruction is affected by multiple proton-proton interactions occurring during the same bunch crossing and by remnants of electronic signals from previous bunch crossings in the detectors (pileup). To simulate pileup, inelastic pp events are generated using Pythia 8 [25] with the 4C tune [26] and MRST LO * * proton PDF set [27]. The number of minimum-bias events overlaid on each signal event is chosen to reproduce the distribution of the average number of simultaneous pp collisions μ in an event. During the 2011 data-taking period μ changed from 5 to 18 with increasing instantaneous luminosity.
To estimate the uncertainties in the modelling of the hard scattering, hadronisation, the underlying event and of parton showers, events are also simulated using Alpgen [28], a multi-leg LO MC simulation, with up to six final-state partons in the matrix element calculations, interfaced to Herwig 6.5.10 [29-31] using the AUET2 tune [21] with the CTEQ6L1 PDF set [32] for parton showers and Jimmy 4. 31 [33] for the underlying event.
The outputs from these event generators are passed to the detector simulation [34], based on Geant4 [35]. Simulated events are digitised [36,37] to model the detector responses, and then reconstructed using the same software as used to process the data.

Data selection and jet calibration
This analysis is based on data collected with the ATLAS detector in the year 2011 during periods with stable pp collisions at √ s = 7 TeV in which all relevant detector components were operational. The resulting data sample corresponds to an integrated luminosity of 4.51 ± 0.08 fb −1 [38].
The presence of at least one primary vertex (compatible with the position of the beam spot), reconstructed using two or more tracks with p T > 500 MeV, is required to reject cosmic ray events and beam-related backgrounds. The primary vertex with the largest sum of squared transverse momenta of associated tracks is used as the interaction point for the analysis.
Due to the high instantaneous luminosity and a limited detector readout bandwidth, a set of single-jet triggers with increasing transverse energy (E T ) thresholds is used to collect data events with jets. Only a fraction of the events that fired the trigger are actually recorded. The reciprocal of this fraction is the prescale factor of the trigger considered. The triggers with lower E T thresholds were prescaled with higher factors and only the trigger with the highest E T threshold remained unprescaled during the whole data-taking period. The prescale factors are adjusted to keep the jet yield approximately constant as a function of E T .
An event must pass all three levels of the jet trigger system. The trigger is based on the E T of jet-like objects. Level-1 provides a fast hardware decision based on the summed E T of calorimeter towers using a sliding-window algorithm. Level-2 performs a simple jet reconstruction in a geometric region around the object that fired the Level-1 trigger. Finally, a full jet reconstruction using the anti-k t algorithm with R = 0.4 is performed over the entire detector by the third level trigger.
The trigger efficiencies are determined as a function of m j j j in each bin of |Y * | separately for R = 0.4 and R = 0.6 jet radius parameters. They are evaluated using an unbiased sample of events that fired the jet trigger with a p T = 30 GeV threshold at the EF level. This trigger is fully efficient in events with a leading jet passing the three-jet analysis requirements. For every |Y * | bin, the full range of three-jet mass is divided into subranges, each filled by only one of the several single-jet triggers. Triggers are used only where the trigger efficiency is above 99 %. Moreover, the lower m j j j bound for each trigger is shifted up by 15 % from the 99 % efficiency point to avoid any possible biases from the trigger strategy chosen for this measurement. This shift leads to a negligible increase in the statistical error on the measured cross-sections, compared to the total uncertainty.
Since the EF reconstructs jets with a radius parameter R = 0.4, the p T threshold at which the trigger for jets defined with R = 0.6 becomes fully efficient is significantly higher than for R = 0.4 jets. Using the same trigger subranges for both jet sizes would reduce the number of events with anti-k t R = 0.4 jets. To take advantage of the lower p T at which triggers are fully efficient for R = 0.4 jets, different assignments between triggers and m j j j ranges are considered for these jets and jets reconstructed with R = 0.6. After events are selected by the trigger system, they are fully reconstructed offline. The input objects to the jet algorithm are three-dimensional topo-clusters [39]. Each topocluster is constructed from a seed calorimeter cell with energy |E cell | > 4σ , where σ is the width of the total noise distribution of the cell from both the electronics and pileup sources. Neighbouring cells are added to the topo-cluster if they have |E cell | > 2σ . At the last step, all neighbouring cells are added. A local hadronic calibration (LC) that accounts for inactive material, out-of-cluster losses for pions, and calorimeter response is applied to clusters identified as hadronic by their energy density distribution [40]. The LC improves the topo-cluster energy resolution, and the jet clustering algorithm propagates this improvement to the jet level. The LC is validated using single pions in the combined testbeam [40].
Each topo-cluster is considered as a massless particle with an energy E = E cell , and a direction given by the energyweighted barycentre of the cells in the cluster with respect to the geometrical centre of the ATLAS detector. The fourmomentum of an uncalibrated jet is defined as the sum of the four-momenta of the clusters making up the jet. The jet is then calibrated in four steps: 1. An estimated mean additional energy due to pileup is subtracted using a correction derived from MC simulation and validated in situ using track-jets in dijet events and photons in γ -jet events as a function of the average number of pp collisions in the same bunch crossing, μ , the number of primary vertices, N PV , and jet η [41]. Here, track-jets are reconstructed from all tracks associated to the primary vertex using the anti-k t jet algorithm. 2. The direction of the jet is corrected such that the jet originates from the selected hard-scatter vertex of the event instead of the geometrical centre of ATLAS. 3. The energy and the position of the jet are corrected for instrumental effects (calorimeter non-compensation, additional inactive material, effects due to the magnetic field) using correction factors obtained from MC simulation. The jet energy scale is restored on average to that of the particle-level jet. For the calibration, the particle-level jet does not include muons and non-interacting particles. 4. An additional in situ calibration is applied to correct for residual differences between the MC simulation and data, derived by combining the results of dijet, γ -jet, Z -jet, and multi-jet momentum balance techniques.
The full calibration procedure is described in detail in Ref. [7].
Data-taking in the year 2011 was affected by a read-out problem in a region of the LAr calorimeter, causing jets in this region to be poorly reconstructed. In order to avoid a bias in the spectra, events with any of the three leading jets falling in the region −0.88 < φ < −0.5 were rejected. Approximately 15 % of events are removed by this requirement. This inefficiency is corrected for using MC simulation (cf. Sect. 6).
The three leading jets are required to satisfy the "medium" quality criteria as described in Ref. [42], designed to reject cosmic-rays, beam-halo particles, and detector noise. More than 5.3(2.5) × 10 6 three-jet events are selected with radius parameter R = 0.4(0.6).

Data unfolding
The three-jet cross-sections as a function of m j j j are obtained by unfolding the data distributions, and correcting for detector resolutions and inefficiencies. This procedure includes a correction for the undetected presence of muons and neutrinos from hadron decays in jets. The unfolding procedure is based on the iterative, dynamically stabilised (IDS) unfolding method [43]. Further details can be found in Ref. [2]. To account for bin-to-bin migrations, a transfer matrix is built from the MC simulation, relating the particle-level and reconstruction-level three-jet masses. The reconstruction-level to particle-level event association is done in the m j j j -|Y * | plane, such that only a requirement on the presence of a threejet system is made. Since bin-to-bin migrations are usually due to jet energy smearing of the three-jet mass, and less often due to jet angular resolution, the migrations across |Y * | bins are negligible and the unfolding is performed separately in each |Y * | bin.
The data are unfolded to the particle level using a threestep procedure where i ( j) is the particle-level (reconstruction-level) bin index, and N P i (N R i ) is the number of particle-level (reconstruction-level) events in bin i. The quantities R i ( P i ) are the fractions of reconstruction-level (particle-level) events matching (associated with) particle-level (reconstructionlevel) events in each bin i. These efficiencies are used to correct for the matching inefficiency at the reconstruction and particle level, respectively. The element A i j of the transfer matrix is the probability for a reconstruction-level event in bin j to be associated with a particle-level event in bin i. It is used to unfold the reconstruction-level spectrum for detector effects.
A data-driven closure test is used to evaluate the bias in the unfolded data spectrum shape due to mis-modelling of the reconstruction-level spectrum shape in the MC simulation. The transfer matrix is improved through a series of iterations, where the particle-level distribution from simulation is reweighted such that the reconstruction-level distribution from simulation matches the data distribution. The modified reconstruction-level MC simulation is unfolded using the original transfer matrix, and the result is compared with the modified particle-level spectrum. The resulting bias is considered as a systematic uncertainty. For the analyses in this paper, one iteration is used, which leads to a bias in closure tests of less than one percent.
The statistical uncertainties in the unfolded results are estimated using pseudo-experiments. Each event in the data and in the MC simulation is counted n times, where n is sampled from a Poisson distribution with a mean of one. A fluctuated transfer matrix and efficiency corrections are calculated as the average over these pseudo-experiments in MC simulation. Then, each resulting pseudo-experiment of the data spectrum is unfolded using the fluctuated transfer matrix and efficiency corrections. Finally, the covariance matrix between bins of measured m j j j cross-section is calculated using the set of unfolded pseudo-experiments of the data. The random numbers for the pseudo-experiments are generated using unique seeds. The dijet [6] and inclusive jet [5] crosssection measurements use the same unique seeds to evaluate the statistical uncertainties. In this way, the statistical uncertainty and bin-to-bin correlations in both the data and the MC simulation are encoded in the covariance matrix and the statistical correlation between different measurements can be taken into account in combined fits.

Experimental uncertainties
The uncertainty in the jet energy scale (JES) calibration is the dominant uncertainty in this measurement. The uncertainties in the central region are determined using a combination of the transverse momentum balance techniques, such as Z -jet, γ -jet and multi-jet balance measurements performed in situ. In each of the methods, the uncertainties in the energy of the well-measured objects, e.g. Z /photon or system of lowp T jets, are propagated to the energy of the balancing jet. The JES uncertainty in the central region is propagated to the forward region using transverse momentum balance between a central and a forward jet in events with two jets. The difference in the balance observed between MC simulation samples generated with Pythia and Herwig is treated as an additional uncertainty in the forward region. The JES uncertainty in the highp T range is evaluated using the in situ measurement of the single isolated hadron response [44]. The total JES uncertainty is described by the set of fully correlated in p T independent uncertainty sources. Complete details of the JES derivation and its uncertainties can be found in Ref. [7].
The uncertainty in the p T of each individual jet due to the JES calibration is between 1 and 4 % in the central region (|η| < 1.8), and increases to 5 % in the forward region (1.8 < |η| < 4.5).
The uncertainties due to the JES calibration are propagated to the measured cross-sections using the MC simulation. The energy and p T of each jet in the three-jet sample are scaled up or down by one standard deviation of a given uncertainty component, after which the luminosity-normalised three-jet event yield is measured from the resulting sample. The yields from the nominal sample and the samples where all jets were scaled up and down are unfolded, and the difference between each of these variations and the nominal result is taken as the uncertainty due to that JES uncertainty component. For example, the uncertainty in the three-jet cross-section in the 8 < |Y * | < 10(|Y * | < 2) bin due to the LAr electromagnetic energy scale uncertainty increases from 2(3) to 10(8) % with the m j j j increasing from 1(0.4) TeV to 4(3) TeV. In the same |Y * | bins, the uncertainty in the three-jet cross-section due to the uncertainty in the jet energy measurements in the forward region varies from 15(4) to 30(0.5) %, as a function of m j j j . Since the sources of JES calibration uncertainty are uncorrelated with each other by construction, the corresponding uncertainty components in the cross-section are also taken as uncorrelated.
Each jet is affected by the additional energy deposited in the calorimeters due to pileup effects. Additional energy due to pileup is subtracted during the jet energy calibration procedure [7]. To check for any residual pileup effects in the measured cross-sections, the luminosity-normalised threejet yields in all three-jet mass and rapidity-separation bins are split into bins of different pileup conditions under which the data were collected. No statistically significant deviation from the nominal result is observed.
The jet energy resolution (JER) is measured in the data using the bisector method in dijet events [45], where good agreement with the MC simulation is observed. The uncertainty in the JER is affected by selection parameters for jets, such as the amount of nearby jet activity, and depends on both jet p T and jet η.
Jet angular resolution (JAR) is studied by matching particle-level jets to reconstruction-level jets in simulation. Jets are matched by requiring that the angular distance R = ( φ) 2 + ( y) 2 between the particle-level and reconstruction-level jet is less than the jet radius parameter. The angular resolution is obtained from a Gaussian fit to the distribution of the difference of reconstruction-level and particle-level jet rapidity.
The difference between the JAR determined from the nominal MC simulation and that from the Alpgen sample is taken as a systematic uncertainty. The resolution varies between 0.005 radians and 0.03 radians depending on the jet η and p T values. The JAR uncertainty is about 10-15 % for p T < 150 GeV and decreases to ∼ 1 % for p T > 400 GeV. The jet angular bias is found to be negligible.
The JER and JAR uncertainties are propagated to the measured cross-section through the unfolding transfer matrix. The energy and direction of each jet in the MC sample are smeared according to their uncertainties. To avoid being limited by statistical fluctuations this procedure is repeated 1000 times in each event. The average transfer matrix derived from these pseudo-experiments is used to unfold the threejet yields, and the deviation from the three-jet yield unfolded using the nominal transfer matrix is taken as a symmetrised systematic uncertainty.
The uncertainty due to the jet reconstruction inefficiency as a function of jet p T is estimated by comparing the efficiency for reconstructing a calorimeter jet, given the presence of an independently measured track-jet of the same radius, in data and in MC simulation [7,46]. Since this method relies on tracking, its application is restricted to jets with |η| < 1.9 to ensure that both the R = 0.4 and R = 0.6 jets are fully within the tracker acceptance. For jets with p T > 50 GeV, relevant for this analysis, the reconstruction efficiency in both the data and the MC simulation is found to be 100 % for this rapidity region, leading to no additional uncertainty. The same efficiency is assumed for the forward region, where jets of a given p T are more energetic and, therefore, their reconstruction efficiency is expected to be at least as good as that of jets in the central region.
The efficiencies for single-jet selection using the "medium" criteria agree within 0.25 % in data and MC simulation [42]. Because three jets are considered for each event selected for the analysis, a 0.75 % systematic uncertainty in the cross-section is assigned.
The impact of a possible mis-modelling of the shape of m j j j spectra in MC simulation, introduced through the unfolding as described in Sect. 6, is also included. The luminosity uncertainty is 1.8 % [38] and is fully correlated between all data points.
The total experimental uncertainty in the three-jet crosssection is summarised in Fig. 1. The total uncertainty ranges from 8-10 % at low three-jet mass to 28 % at high three-jet mass for the range |Y * | < 6 (see Appendix), and increases slightly for larger |Y * | bins. In the 8 < |Y * | < 10 bin the total uncertainty ranges from 18 to 38 %, where it is dominated by the jet energy scale uncertainty component for forward jets.

Theoretical predictions and uncertainties
The NLO QCD predictions by the parton-level MC crosssection calculator NLOJET++ [47], corrected for hadronisation effects and underlying-event activity using Monte Carlo simulation with Perugia 2011 tune [19] of Pythia 6, are compared to the measured three-jet cross-sections. Total systematic uncertainty in the three-jet cross-section for anti-k t R = 0.6 jets as a function of m j j j (a) in |Y * | < 2 and (b) 8 < |Y * | < 10 bins. The bands shows the uncertainties due to jet energy scale, jet angular resolution, jet energy resolution and the combined uncertainty due to jet quality selection and unfolding. The outer band represents the total experimental uncertainty 8.1 Fixed-order predictions The fixed-order QCD calculations are performed with the NLOJET++ program interfaced to APPLgrid [48] for fast convolution with various PDF sets. The renormalisation (Q R ) and factorisation (Q F ) scales are set to the mass of the three- To estimate the uncertainty due to missing higher-order terms in the fixed-order perturbative expansion, the renormalisation scale is varied up and down by a factor of two. The uncertainty due to the dependence of the theoretical predictions on the factorisation scale, which specifies the separation between the short-distance hard scattering and long-distance non-perturbative dynamics, is estimated by varying the factorisation scale up and down by a factor of two. All permu-tations of these two scale choices are considered, except the cases where the scales are shifted in opposite directions. The maximum deviations from the nominal prediction are taken as the scale uncertainty. The scale uncertainty is generally 10−20 % depending on the m j j j .
The multiple uncorrelated uncertainty components of each PDF set, as provided by the various PDF analyses, are also propagated through the theoretical calculations. The PDF groups generally derive these from the experimental uncertainties in the data used in the fits. For the results shown in Sect. 9, the standard Hessian sum in quadrature [55] of the various independent components is calculated taking into account asymmetries of the uncertainty components. The NNPDF 2.3 PDF set is an exception, where uncertainties are expressed in terms of replicas instead of independent components. These replicas represent a collection of equally likely PDF sets, where the data used in the PDF fit were fluctuated within their experimental uncertainties. For the plots shown in Sect. 9, the uncertainties in the NNPDF 2.3 PDF set are evaluated as the RMS of the replicas in each bin of m j j j , producing equivalent PDF uncertainties in the theoretical predictions. These uncertainties are symmetric by construction. Where needed, the uncertainties of PDF sets are rescaled to the 68 % confidence level (CL). HERAPDF provides three types of uncertainties: experimental, model and parameterisation. The three uncertainty sources are added in quadrature to get a total PDF uncertainty.
The uncertainties in the cross-sections due to the strong coupling, α s , are estimated using two additional proton PDF sets, for which different values of α s are assumed in the fits, such that the effect of the strong coupling value on the PDFs is included. This follows Ref. [56]. The resulting uncertainty is approximately 3 % across all three-jet mass and |Y * | ranges considered.
The scale uncertainties are dominant in low and intermediate three-jet mass regions, while the PDF uncertainties become dominant at high m j j j . The uncertainties in the theoretical predictions due to those on the PDFs range from 5 % at low m j j j to 30 % at high three-jet mass for the range of |Y * | values up to four. For the values of |Y * | between four and ten, the PDF uncertainties reach 40-80 % at high three-jet mass, depending on the PDF set and the |Y * | value.

Non-perturbative effects
Non-perturbative corrections (NPC) are evaluated using leading-logarithmic parton-shower generators, separately for each value of the jet radius parameter. The corrections are calculated as bin-by-bin ratios of the three-jet differential cross-section at the particle level, including hadronisation and underlying-event effects, to that at parton-level after the parton shower (before the hadronisation process starts) with the underlying-event simulation switched off. The nominal    Fig. 3 The three-jet double-differential cross-section as a function of m j j j in bins |Y * |, as denoted in the legend. The jets are identified using the anti-k t algorithm with R = 0.4. For convenience, the cross-sections are multiplied by the factors indicated in the legend. Also shown is the comparison with the NLOJET++ prediction with the CT 10 PDF set corrected for non-perturbative effects. The statistical uncertainties are smaller than the size of the symbols. Where visible, the sum in quadrature of the statistical and experimental systematic uncertainties is plotted  Fig. 4 The three-jet double-differential cross-section as a function of m j j j in bins |Y * |, as denoted in the legend. The jets are identified using the anti-k t algorithm with R = 0.6. For convenience, the cross-sections are multiplied by the factors indicated in the legend. Also shown is the comparison with the NLOJET++ prediction with the CT 10 PDF set corrected for non-perturbative effects. The statistical uncertainties are smaller than the size of the symbols. Where visible, the sum in quadrature of the statistical and experimental systematic uncertainties is plotted corrections are calculated using Pythia 6 with the Perugia 2011 tune. The non-perturbative corrections as a function of three-jet mass are shown in Fig. 2 for the range |Y * | < 2 for R = 0.4 and R = 0.6 jets. The NPC are smaller than 10 % in all m j j j and |Y * | bins. The uncertainties in the non-perturbative corrections, arising from the modelling of the hadronisation process and the underlying event, are estimated as the maximum deviations of the corrections from the nominal ones, using the  57,58] with the UE-EE-3 tune [59] using the CTEQ6L1 set. The uncertainty in the non-perturbative corrections ranges up to ∼ 10 % depending on the three-jet mass in all |Y * | bins.
The total theoretical uncertainty is calculated as a sum in quadrature of PDF, scale, α s and NPC uncertainties.

Cross-section results
Measurements of the double-differential three-jet crosssections as a function of the three-jet mass in various ranges of |Y * | are shown in Figs. 3 and 4 for anti-k t jets with values of the radius parameter R = 0.4 and R = 0.6, respectively. The cross-section decreases rapidly as a function of the three-jet mass. The NLO QCD calculations using NLO-JET++ with the CT 10 PDF set corrected for non-perturbative effects are compared to the measured cross-sections. Good agreement between the data and the theoretical predictions is found over the full kinematic range, covering almost seven orders of magnitude in the measured cross-section values.
The ratios of the theoretical predictions calculated with various PDF sets to the measured cross-sections are presented in Figs The three-jet cross-sections are well described by the calculations that use CT 10, NNPDF 2.3, GJR 08, MSTW 2008 and HERAPDF 1.5 PDFs. Disagreement between data and the predictions using ABM 11 PDFs is observed for most of the cross-sections measured with both jet radius parameters.
For all PDF sets, the predictions for anti-k t R = 0.4 jets agree well with measured cross-sections, while the calculations that use the ABM 11 PDF set are systematically below all other theory curves. Theory predictions for anti-k t R = 0.6 jets underestimate the data across the full m j j j -|Y * | plane. This shift is within the experimental and theoretical uncertainties. The jet radius dependence of theory-to-data ratios is similar for all PDF sets considered, demonstrating that this tendency is independent of the assumptions made in different PDF determinations.

Conclusions
Cross-section measurements of three-jet production in pp collisions at 7 TeV centre-of-mass energy as a function of  Fig. 8 The ratio of NLO QCD predictions, obtained by using NLO-JET++ with different PDF sets (NNPDF 2.3, ABM 11, HERAPDF 1.5) and corrected for non-perturbative effects, to data as a function of m j j j in bins of |Y * |, as denoted in the legend. The ratios are for jets identified using the anti-k t algorithm with R = 0.6. The experimental error bands are centered at one and designate the relative statistical (thin dashed line) and total (statistical and systematic uncertainties added in quadrature) experimental uncertainties (thick solid line). The theoretical predictions are represented by thick lines with the hatched or filled band around it. The line show the central values and the band represent the total theory uncertainty the three-jet mass, in bins of the sum of the absolute rapidity separations between the three leading jets are presented. Jets are reconstructed with the anti-k t algorithm using two values of the radius parameter, R = 0.4 and R = 0.6. The measurements are based on the full data set collected with the ATLAS detector during 2011 data-taking at the LHC, corresponding to an integrated luminosity of 4.51 fb −1 . The measurements are corrected for detector effects and reported at the particle level. The total experimental uncertainty in these measurements is dominated by the jet energy scale calibration uncertainty. The measurement uncertainties are smaller than, or similar to, those in the theoretical predictions.
The measurements probe three-jet masses up to ∼ 5 TeV and are well described by perturbative QCD at NLO accuracy across the full m j j j -|Y * | plane. The comparison of NLO QCD predictions corrected for non-perturbative effects to the measured cross-sections is performed using several modern PDF sets. The data are well described by the theoretical predictions when using CT 10, NNPDF 2.3, HERAPDF 1.5, GJR 08 and MSTW 2008 PDFs. The theoretical calculations based on the ABM 11 PDFs are systematically below all the other predictions. Comparison of measured cross-sections to theoretical predictions for two different jet radius parameters shows good agreement for R = 0.4 jets but shifted theory-to-data ratios for R = 0.6 jets. This shift is covered by the experimental and theoretical uncertainty bands and it has only a minor dependence on the PDF set used. are the statistical uncertainties in the data (MC simulation). The γ components are the uncertainty in the jet energy calibration from the in situ, the pileup, the close-by jet, and flavour components.
The u components show the uncertainty for the jet energy and angular resolution, the unfolding, the quality selection, and the luminosity. While all columns are uncorrelated with each other, the in situ, pileup, and flavour uncertainties shown here are the sum in quadrature of multiple uncorrelated components m j j j (bin #) m j j j -range (TeV)   are the statistical uncertainties in the data (MC simulation). The γ components are the uncertainty in the jet energy calibration from the in situ, the pileup, the close-by jet, and flavour components.
The u components show the uncertainty for the jet energy and angular resolution, the unfolding, the quality selection, and the luminosity. While all columns are uncorrelated with each other, the in situ, pileup, and flavour uncertainties shown here are the sum in quadrature of multiple uncorrelated components m j j j (bin #) m j j j -range (TeV)  Table 3 Measured double-differential three-jet cross-section,   Table 4 Measured double-differential three-jet cross-section, σ , for R = 0.6 jets and 2 ≤ |Y * | < 4, along with uncertainties in the measurement. All uncertainties are given in %, where δ data stat (δ MC stat ) are the statistical uncertainties in the data (MC simulation). The γ components are the uncertainty in the jet energy calibration from the in situ, the pileup, the close-by jet, and flavour components. The u components show the uncertainty for the jet energy and angular resolution, the unfolding, the quality selection, and the luminosity. While all columns are uncorrelated with each other, the in situ, pileup, and flavour uncertainties shown here are the sum in quadrature of multiple uncorrelated components m j j j (bin #) m j j j -range (TeV)     Table 6 Measured double-differential three-jet cross-section, σ , for R = 0.6 jets and 4 ≤ |Y * | < 6, along with uncertainties in the measurement. All uncertainties are given in %, where δ data stat (δ MC stat ) are the statistical uncertainties in the data (MC simulation). The γ components are the uncertainty in the jet energy calibration from the in situ, the pileup, the close-by jet, and flavour components. The u components show the uncertainty for the jet energy and angular resolution, the unfolding, the quality selection, and the luminosity. While all columns are uncorrelated with each other, the in situ, pileup, and flavour uncertainties shown here are the sum in quadrature of multiple uncorrelated components m j j j (bin #) m j j j -range (TeV)  Table 7 Measured double-differential three-jet cross-section,