Advances in High Energy Physics

Volume 2013 (2013), Article ID 585809, 11 pages

http://dx.doi.org/10.1155/2013/585809

Research Article

Subjet Multiplicities at LHC Energies and the QCD Color Factor Ratio

M. Kaur and Anter P. Kaur

Department of Physics, Panjab University, Chandigarh 160014, India

Received 22 August 2013; Revised 3 November 2013; Accepted 5 November 2013

Academic Editor: Kingman Cheung

Copyright © 2013 M. Kaur and Anter P. Kaur. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Subjet multiplicity distributions of two jets with highest transverse momenta are simulated in proton-proton collisions at LHC energies of 7 TeV and 14 TeV. Jets are defined via the software package FastJet version 3.0.3 which includes the implementation of sequential recombination algorithm with a jet size of . Subjets are resolved by rerunning the algorithm on the jets and using the resolution cutoff of . The subjet multiplicity distributions and their average values are measured as a function of the jet momenta and upto rapidities . The results are used to calculate the ratio, .

1. Introduction

One of the main motivations of studying highly energetic proton-proton (p-p) collisions at Large Hadron Collider (LHC) at CERN has been the discovery of Higgs Boson(s). The aim has recently been achieved with the discovery of a Higgs Boson. However, many more studies are needed to establish its properties. In addition, searches for new particles will continue for many more years.

The proton-proton collisions at the LHC are dominated by jet production which constitute a large background to potential signals of new physics. Quantum chromodynamics (QCD) describes well the dynamics of jet production in terms of partons. In QCD, quarks and gluons have different color charges and hence different coupling strengths. The values for the color factors originate directly from the symmetry group SU of QCD, calculated to be for the gluon and for the quark [1].

In hadron-hadron collisions, quarks and gluons fragment to produce collimated streams of particles conventionally called “jets.” A jet, characterized by its energy and momentum, is reconstructed according to a clustering prescription called a “jet algorithm” from the four vectors of constituents that may represent tracks, energy deposits in the calorimeter, reconstructed particle candidates, particles in a Monte Carlo (MC) event generator, or partons of a theory calculation. In dijet events, the two jets leading in transverse momentum can be associated with the two partons at leading order (LO) in perturbative QCD. The internal structure of these jets is then expected to depend mainly on the type of the primary parton, that is, either (anti-) quark or gluon, from which they originated. QCD predicts that gluons, because of their larger color factor, fragment more than quarks. Consequently, gluon-initiated jets become broader and exhibit a larger constituent multiplicity than quark jets. The ratio of the constituent multiplicity for gluon jets versus quark jets at LO is asymptotically given by the ratio of their color factor as. Effects of higher orders in the strong coupling constant, however, could change this theoretical prediction.

In order to study the jet substructure, an experimental observable called subjet multiplicity, , is found to be very useful. It is defined as the number of subjets that can be resolved within a jet by reclustering the jet constituents with the same clustering algorithm. However, the spatial resolution is chosen to be finer than the one used in the jet reconstruction. In this way an infrared- and collinear-safe measure can be defined that is usable experimentally as well as in perturbative QCD. Various collinear and infrared safe jet algorithms exist [2].

Following a sequential recombination procedure, jets are defined byalgorithm which starts clustering first the softest objects, trying to undo the effects of parton showering [3–5] and is well suited to explore the jet substructure.algorithm is implemented in the FastJet package [6–8].

In the present work, samples used to measure the subjet multiplicity at two center of mass energies (c.m.), 7 TeV and 14 TeV, corresponding to the current and future LHC energies, respectively, are generated from the event generator PYTHIA8 [9]. A sample of 7 million events at each of the two c.m. energies is generated. The dijet events are selected in which the average subjet multiplicityis studied as a function of the jet   in the range of 100 to 1100 GeV and in the pseudorapidity interval, where pseudorapiditywithas the azimuthal angle of the jet.

2. Jet Algorithm and the Subjet Multiplicity

Thealgorithm [5] clusters all input objects,, according to the following iterative steps.

(i)For each pairof objects find

where= 0.6 is the jet size inspace

(ii)For each object calculate the beam distance

(iii)If the minimumof all possibleandis, then merge objectsandinto a single object by 4-momentum vector addition. If it isthen remove objectfrom the list and define it to be a final jet.(iv)Repeat the above steps until no clustering objects are left.

2.1. Subjet Multiplicities

Subjets are resolved within a jet by repeating the application ofcluster algorithm described above. The clustering is terminated when allare above the quantity (jet). All remaining objects are called subjets. While the stopping parameter defines the hard scale of the process, the parameteris known as the resolution parameter. The subjet structure depends upon the values chosen for the resolution parameter. The mean subjet multiplicityis defined as the average number of subjets in a jet at a given value of:

By definitionand. The mean subjet multiplicity is measured forand , whereis the subjet resolution parameter. A study to optimise the value ofat 7 TeV was done in an earlier work, documented in [10, 11] and the same used here. This optimum value ofmatches with the one used for the definition of subjets in [12]. In addition two subjets within a jet are resolved if they are well separated inspace.

3. Determination of Ratio

In the p-p collisions, the final state stable particles are arranged into jets by using a jet algorithm, as described above. These particles are produced in the hadronisation of hard partons. In a mixed sample of quark and gluon initiated particle jets,   is the subjet multiplicity which can be written as a linear combination of subjet multiplicity in gluon jets and quark jets :

whereis the fraction of gluon jets and () is the fraction of quark jets in the mixed sample. The gluon jet fractionis the number of outgoing gluons that pass the selection cuts divided by the total number of outgoing partons that pass the selection cuts. Considering the above equation for two similar samples of jets at  (is the center of mass energy (c.m.))TeV and 14 TeV and assuming and to be independent of, we get

The solutions are

where and are extracted from the total subjets in mixed samples of quark and gluon initiated jets in the data at and 7 TeV, respectively, andandare the gluon jet fractions at the two energies. This method relies on the extraction of quark and gluon fractions from the Monte Carlo information, obtained by tagging the particle jets with partons. This way the knowledge of the two gluon jet fractions is obtained. The subjet multiplicity distributions can be characterized by the mean values such that and gives the average number of subjets in a gluon and quark jet, respectively. The gluon jets are compared to the quark jets by having a ratio:

4. Event Samples and Jet Selection

To study the subjet multiplicity in highjets and to calculate the color factor ratio, events are generated using PYTHIA8 with CTEQ6L1 PDFs at TeV and 14 TeV. At each energy, 1 million events are generated in 4 bins ofin the pseudorapidity range. The bins inare 50–300, 300–550, 550–800, and 800–1050 GeV. One million events are generated in each pseudorapidity rangeand. One million events are also generated separately in each bin of jet pseudorapiditywhich are 0.0–0.5, 0.5–1.0, 1.0–1.5, and 1.5–2.0 for leading jets. Thus in total we have generated 7 million events for each energy for the present study.

Quark and gluon jets are unambiguously defined only in the leading order QCD where there are only two jets, that is, two partons in the final state. So a sample of dijet events in hard QCD scattering events is defined by selecting the two jets leading in in each event with the following selection cuts:

(i)100 GeV,(ii)(requiring the jets to be back-to-back in the azimuthal plane),(iii).

Foralgorithm, mean subjet multiplicity is measured for cone radius, and.

5. Discrimination between Quark and Gluons Jets

To determine the ratio of color factor, jet substructure of gluon and quark initiated jets is studied in terms of subjet multiplicities. For the Monte Carlo PYTHIA8 events used here, the quark and gluon jets are identified using parton level information. The direction of each primary parton is determined after its perturbative evolution has terminated. The jet closest to the direction of the primary parton (quark/antiquark or gluon) tagged to be the quark/gluon jet, respectively. The two highest jets are identified as quark or gluon initiated jets via the tagging in the spatial coordinates and the kinematical cuts. Thus we implement these cuts by allowing the maximum distance inspace as. The maximum relative deviation in transverse momentum between the jet and the associated parton with respect to the jet transverse momentum is(is abs (partonjet)/jet). Figure 1 shows the ratio for quark and gluon jets matched, respectively, to partons for TeV and TeV. At each energy the matching of particle jets with the partons is found to be good as all the jets have the ratio less than 3.

Figure 1: Distribution offor the matched quark and gluon jets with GeV at TeV compared to TeV foralgorithm.

Figure 2 shows the fraction of gluon initiated jets via matching of the two particle jets leading in to the hard partons of the LO process as a function of jet for the two different pseudorapidity bins,and. As expected the fraction of gluon jets decreases with increasing jet in both the cases. Also, quark jets make up a dominant fraction of the jets only at high jet transverse momenta.

Figure 2: Fraction of gluon initiated jets differentiated via matching of the two particle leading jets to the hard partons of the leading-order process as a function of jet fromalgorithm at TeV compared to TeV in the pseudorapidity bins ofand.

6. Systematic Uncertainties

The uncertainties in the subjet multiplicities are of two categories:

(i)statistical uncertainties,(ii)systematic uncertainties.

For the present analysis using Monte Carlo samples, the uncertainties are mainly on account of the statistics which are found to be negligibly small. The systematic uncertainties arise due to the change in gluon fraction due to minimum imposed. This is measured as follows.

(i) Gluon Jet fraction: The uncertainty in gluon jet fraction introduces a systematic error on the subjet multiplicity. To calculate this, the gluon jet fraction is obtained as a function of minimum jet (GeV) by varying minimum jet from 90 to 110 GeV. This introduces an error of   on gluon fraction at TeV and on gluon fraction at TeV. If the gluon fraction increases or decreases at both the energies, the change in ratio is not much effective. But the increase in gluon jet fraction at one energy and decrease at other or vice versa introduces a largest change in ratio , giving minimum and maximum value of . The difference of the calculated value given in Table 5 from the minimum value and maximum value of gives the systematic uncertainty due to gluon fraction as.

The same analysis for measuring , when performed on the data from a detector such as CMS at LHC, would involve the following systematic errors also. In one such measurement with 36data used for measuring the subjet multiplicity at TeV energy [10], we estimated the following uncertainties.

(ii) Jet-Energy calibration: The maximum systematic uncertainty in Jet-Energy calibration was estimated to be 5% depending upon jet and jet. This introduces an error of about 2% on the subjet multiplicity average value.

(iii) Jet-Energy resolution (JER): The uncertainty on the jet resolution (JER) was estimated to be less than 10% and this results in less than 0.25% uncertainty in the subjet multiplicity distributions.

(iv) Uncertainties due to the different physics models employed for the simulations: The uncertainty due to different physics models implemented in the Monte Carlo generators is estimated to be less than 3%.

Compounding these systematic errors (ii)–(iv), the maximum uncertainty in the subjet multiplicity of 3.28 at 7 TeV was estimated to be 0.10. The details can be found in [10]. We use the same percent value of each of the above mentioned systematic uncertainties in subjet multiplicities at the two energies which will be reflected in ratio . The systematic errors on ratio are given in Table 6 and added in quadrature to have the total systematic uncertainty in the calculated value of the ratio .

7. Results and Discussions

Figure 3 shows jet constituent multiplicity distributions foralgorithm at TeV and at TeV in the two pseudorapidity bins,and. It may be observed that the average number of jet constituents is larger at 14 TeV. The subjet fractions for the two pseudorapidity regions and for the two energies are shown in Figure 4. Inalgorithm, the procedure to resolve subjets is dependent on of jet and the of subjets in a particular jet is distributed in such a way that there is no sudden decrease. Thus as shown in Figure 4, we get monotonically decreasing values for subjet fraction.

Figure 3: Distributions of the jet constituent multiplicity in particle jets having100 GeV fromalgorithm at TeV compared to TeV in the pseudorapidity bins ofand.

Figure 4: Distributions of the subjet fraction in particle jets having GeV fromalgorithm at TeV compared to TeV in the pseudorapidity bins ofand.

Figure 5 shows the subjet multiplicity distributions for jets matched to gluons in the two pseudorapidity bins and jets matched to quarks in the two pseudorapidity binsand for. The mean values are given in Table 1. Average number of subjets in gluon jets is found to be larger than in the quark jets indicating a higher degree of collimation in quark jets.

Table 1: Average subjet multiplicities for jets matched to gluons and quarks in the bins with  GeV for algorithm at = 7 TeV and = 14 TeV.

7 TeV
14 TeV

Gluon jets
Quark jets
Gluon jets
Quark jets

Stat.

Stat.

Stat.

Stat.

3.39
0.01
2.27
0.01
3.44
0.01
2.27
0.01

3.46
0.01
2.26
0.01
3.45
0.01
2.32
0.01

Figure 5: Subjet multiplicity distributions for jets matched to gluons (top) and quarks (bottom), having GeV, in the pseudorapidity bins of(left) and(right) at TeV and TeV fromalgorithm.

Average subjet multiplicity,dependence on absolute jet pseudorapidity for TeV and at TeV for the wholerange, is shown in Figure 6. Figure 7 shows the subjet multiplicity distributions of jets as a function of the pseudorapidity bins (0.0–0.5, 0.5–1.0, 1.0–1.5, and 1.5–2.0) for fourintervals 50–300, 300–550, 550–800, and 800–1050 GeV. The values are given in Table 2. We observe that from inner to outer rapidities windows, with same rapidity window size, the subjet multiplicity tends to get smaller very slightly.

Table 2: Average subjet multiplicities as a function of in different bins, for jets with 100 GeV for algorithm at = 7 TeV and = 14 TeV.

(GeV)

7 TeV
14 TeV

Stat.

Stat.

50–300

2.99
0.05
3.25
0.06

2.90
0.07
3.06
0.07

2.89
0.09
3.14
0.08

2.71
0.11
3.02
0.09

300–550

2.11
0.01
2.21
0.01

2.06
0.01
2.18
0.01

2.01
0.01
2.17
0.01

1.95
0.02
2.10
0.01

550–800

1.82
0.01
1.94
0.01

1.77
0.01
1.90
0.01

1.70
0.01
1.86
0.01

1.65
0.04
1.80
0.01

800–1050

1.68
0.01
1.81
0.01

1.63
0.01
1.77
0.01

1.51
0.01
1.72
0.01

1.49
0.08
1.64
0.02

Figure 6: Average subjet multiplicityas a function of absolute jet pseudorapidityfor the whole range inat TeV and TeV foralgorithm.

Figure 7: Average subjet multiplicityas a function of absolute jet pseudorapidityin the fourbins at TeV and TeV foralgorithm.

Figure 8 shows the jet versus for pseudorapidityfor quarks, gluons, and matched jets for the two energies under study. The mean values are given in Table 3.

Table 3: Average subjet multiplicities as a function of (GeV) for mixed jets and jets matched to gluons and quarks, for algorithm at = 7 TeV and = 14 TeV.

(GeV)

Stat.

Stat.

Stat.

7 TeV

100–300
2.63
0.01
3.06
0.01
2.05
0.01

300–500
2.03
0.01
2.45
0.01
1.72
0.01

500–700
1.82
0.01
2.25
0.01
1.58
0.01

700–900
1.71
0.01
2.17
0.01
1.53
0.01

900–1100
1.64
0.01
2.11
0.01
1.49
0.01

14 TeV

100–300
2.76
0.01
3.07
0.01
2.13
0.01

300–500
2.11
0.01
2.43
0.01
1.70
0.01

500–700
1.93
0.01
2.28
0.01
1.60
0.01

700–900
1.82
0.01
2.17
0.01
1.56
0.01

900–1100
1.75
0.01
2.11
0.01
1.52
0.01

Figure 8: Average subjet multiplicityas a function of jet in the pseudorapidity range ofat TeV and TeV for  algorithm.

The selected dijet sample consists of mixture of gluon and quark initiated jets as in the LO hard QCD processes. Following the recipe from Section 3, we use (7) to separate these two types and obtain the average subjet multiplicitiesfor gluon jets and quark jets in the five bins of jet. The results are given in Table 5 (the weighted meanand the uncertainties are calculated by using and . These values are then used to calculate the ratio from (8) as given in Table 5.

8. Summary

The subjet multiplicities in proton-proton collisions at TeV and at TeV are estimated in 14 million proton-proton collisions by selecting dijet samples simulated using PYTHIA8.

The jets and subjets have been resolved with thejet algorithm for a jet size of and a subjet resolution cutoff. The subjet multiplicity distribution and its averagehave been determined in the jet range of 100–1100 GeV for jet pseudorapidity. In the whole pseudorapidity range, the average subjet multiplicityfor particle jets is found to decrease from   (stat.) down to   (stat.) for 7 TeV and from   (stat.) down to   (stat.) at 14 TeV with increasing jet . The average subjet multiplicity for jets matched to gluons, , is also found to decrease from   (stat.) down to   (stat.) for 7 TeV and from   (stat.) down to   (stat.) at 14 TeV with increase in jet from 100 GeV to 1100 GeV. The average subjet multiplicity for jets matched to quarks,, is also found to decrease from   (stat.) down to (stat.) for 7 TeV and from   (stat.) down to   (stat.) at 14 TeV with increase in jet .

This is a clear indication of a higher degree of collimation with rising jet . In Table 2, we can observe that from inner to outer rapidities the subjet multiplicity tends to get slightly smaller as well. Evidently, a direct conclusion on the gluon-quark fraction cannot be drawn since the subjet multiplicity not only depends on this fraction but also on details of the Multiple Parton Interactions (MPI), parton shower, and hadronization models. Exploiting data at two different center-of-mass energies, the model dependence can be potentially eliminated. This technique was applied for data at two center-of-mass energies at Tevatron, as described in [12] to measure the ratio . This method has been exploited for present and future LHC energies to predict the value of ratio . The values of average subjet multiplicities in quark and gluon jets in different jet bins are calculated from (7) using the values given in Tables 3 and 4 and the ratio from (8). The results are given in Table 5. The mean value of ratio is calculated by using maximum likelihood estimator of the subjet distributions as   (stat.). The total systematic uncertainty from different sources comes out to be, given in Table 6. So the extracted value of ratio of subjet multiplicities of gluon jets and quark jets is   (stat.)(sys.). This value converges to the ratio of the color factors,. In [17], it is shown that the ratio of the multiplicities from quark and gluon jets should converge to the color factor ratio. The extracted value is in good agreement with the value predicted from the symmetry group SU of QCD and also with the values measured at previous experiments listed in Table 7.

Table 4: Fraction of gluon and quark jets in phase space in different (GeV) bins for algorithm.

(GeV)
7 TeV
14 TeV

Gluon Jets
Quark Jets
Gluon Jets
Quark Jets

Stat.

Stat.

Stat.

Stat.

100–300
0.57
0.01
0.43
0.01
0.68
0.01
0.32
0.01

300–500
0.44
0.01
0.56
0.01
0.56
0.01
0.44
0.01

500–700
0.33
0.01
0.67
0.01
0.48
0.01
0.52
0.01

700–900
0.27
0.01
0.73
0.01
0.42
0.01
0.58
0.01

900–1100
0.22
0.01
0.78
0.01
0.37
0.01
0.63
0.01

Table 5: Average subjet multiplicities (calculated using (7)) and ratio of subjet multiplicities (from (8)) in the phase space in different (GeV) bins for algorithm.

(GeV)
Gluon jets
Quark jets
=
Stat.

Stat.

Stat.

100–300
3.17
0.18
1.90
0.17
2.42
0.25

300–500
2.39
0.06
1.75
0.05
1.87
0.07

500–700
2.33
0.04
1.56
0.03
2.36
0.06

700–900
2.25
0.04
1.52
0.03
2.42
0.06

900–1100
2.19
0.04
1.49
0.02
2.43
0.06

Table 6: Systematic uncertainties on ratio of subjet multiplicities from different sources.

Source
r

Gluon jet fraction
(+0.08)

(−0.05)

Jet-energy calibration
0.15

Jet-energy resolution
0.02

Different physical model
0.22

Total
(+0.28)

(−0.27)

Table 7: Color factor ratio for different experiments.

Experiment
Ratio
Stat.
Sys.

D (2 jet events) [12]
1.91
0.04
0.21

DELPHI (3 jet events) [13]
1.74
0.03
0.10

OPAL (4 jet events) [14]
2.11
0.16
0.28

DELPHI (4 jet events) [15]
2.12
0.35
0.11

ALEPH (4 jet events) [16]
2.20
0.09
0.13

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The authors would like to acknowledge Klaus Rabbertz for his contributions in steering the work reported in [14, 15] from where the estimation of some errors has been adopted in the present paper.

References

R. K. Ellis, W. J. Stirling, and B. R. Webber, QCD and Collider Physics, Cambridge University Press, Cambridge, UK, 1996.

G. P. Salam, “Towards Jetography,” The European Physical Journal C, vol. 67, no. 3-4, pp. 637–686, 2010.

S. Catani, Y. L. Dokshitzerc, M. Olssond, G. Turnocka, and B. R. Webbera, “New clustering algorithm for multijet cross sections in e⁺ e⁻ annihilation,” Physics Letters B, vol. 269, no. 3-4, pp. 432–438, 1991.

S. Catani, Y. L. Dokshitzer, and B. R. Webber, “The -clustering algorithm for jets in deep inelastic scattering and hadron collisions,” Physics Letters B, vol. 285, no. 3, pp. 291–299, 1992.

S. D. Ellis and D. E. Soper, “Successive combination jet algorithm for hadron collisions,” Physical Review D, vol. 48, no. 7, pp. 3160–3166, 1993.

M. Cacciari, G. P. Salam, and G. Soyez, “FastJet user manual,” The European Physical Journal C, vol. 72, 2012.

M. Cacciari, G. P. Salam, and G. Soyez, “fastJet homepage,” http://www.lpthe.jussieu.fr/salam/fastjet/

M. Cacciari and G. P. Salam, “Dispelling the myth for the jet-finder,” Physics Letters B, vol. 641, no. 1, pp. 57–61, 2006.

T. Sjostrand, S. Mrenna, and P. Z. Skands, “A brief introduction to PYTHIA 8.1,” Computer Physics Communications, vol. 178, no. 11, pp. 852–867, 2008.

The CMS Collaboration, “Measurement of the Subjet Multiplicity in Dijet Events from proton-proton collisions at  TeV,” CMS PAS QCD-10-041, 2012

M. Zubin Mehta, M. Kaur, K. Rabbertz, et al., “Measurement of the Subjet Multiplicity in Dijet Events from proton-proton collisions at  TeV,” CMS AN-2010-252, 2011

DØ Collaboration, “Subjet multiplicity of gluon and quark jets reconstructed with the algorithm in fip collisions,” Physical Review D, vol. 65, no. 5, Article ID 052008, 20 pages, 2002.

J. William Gary, “Determination of the QCD color factor ratio from the scale dependence of multiplicity in three jet events,” Physical Review D, vol. 61, no. 11, Article ID 114007, 8 pages, 2000.

R. Akers, G. Alexander, J. Allison, et al., “A measurement of the QCD colour factor ratios and from angular correlations in four-jet events,” Zeitschrift für Physik C, vol. 65, no. 3, pp. 367–377, 1995.

P. Abreu, W. Adam, T. Adye, et al., “Measurement of the triple-gluon vertex from 4-jet events at LEP,” Zeitschrift für Physik C, vol. 59, no. 3, pp. 357–368, 1993.

The ALEPH Collaboration, “A measurement of the QCD colour factors and a limit on the light gluino,” Zeitschrift für Physik C, vol. 76, pp. 1–14, 1997.

E. D. Malaza, “Multiplicity distributions in quark and gluon jets to ,” Zeitschrift für Physik C, vol. 31, no. 1, pp. 143–150, 1986.