# Subjet Multiplicity in Quark and Gluon Jets at DØ

## I Introduction

The Tevatron proton-antiproton collider is a rich environment for studying high energy physics. The dominant process is jet production, described in Quantum Chromodynamics (QCD) by scattering of the elementary quark and gluon constituents of the incoming hadron beams. In leading order (LO) QCD, there are two partons in the initial and final states of the elementary process. A jet is associated with the energy and momentum of each final state parton. Experimentally, however, a jet is a cluster of energy in the calorimeter. Understanding jet structure is the motivation for the present analysis. QCD predicts that gluons radiate more than quarks. Asymptotically, the ratio of objects within gluon jets to quark jets is expected to be in the ratio of their color charges $`C_A/C_F=9/4`$.

## II The $`k_T`$ Jet Algorithm

We define jets in the DØ detector  with the $`k_T`$ algorithm . The jet algorithm starts with a list of energy preclusters, formed from calorimeter cells or from particles in a Monte Carlo event generator. The preclusters are separated by $`\mathrm{\Delta }=\sqrt{\mathrm{\Delta }\eta ^2+\mathrm{\Delta }\varphi ^2}>0.2`$, where $`\eta `$ and $`\varphi `$ are the pseudorapidity and azimuthal angle of the preclusters. The steps of the jet algorithm are:

1. For each object $`i`$ in the list, define $`d_{ii}=E_{T,i}^2`$, where $`E_T`$ is the energy transverse to the beam. For each pair $`(i,j)`$ of objects, also define $`d_{ij}=min(E_{T,i}^2,E_{T,j}^2)\frac{\mathrm{\Delta }_{ij}^2}{D^2}`$, where $`D`$ is a parameter of the jet algorithm.

2. If the minimum of all possible $`d_{ii}`$ and $`d_{ij}`$ is a $`d_{ij}`$, then replace objects $`i`$ and $`j`$ by their 4-vector sum and go to step 1. Else, the minimum is a $`d_{ii}`$ so remove object $`i`$ from the list and define it to be a jet.

3. If any objects are left in the list, go to step 1.

The algorithm produces a list of jets, each separated by $`\mathrm{\Delta }>D`$. For this analysis, $`D=0.5`$.

The subjet multiplicity is a natural observable of a $`k_T`$ jet. Subjets are defined by rerunning the $`k_T`$ algorithm starting with a list of preclusters in a jet. Pairs of objects with the smallest $`d_{ij}`$ are merged successively until all remaining $`d_{ij}>y_{cut}E_T^2(jet)`$. The resolved objects are called subjets, and the number of subjets within the jet is the subjet multiplicity $`M`$. The analysis in this article uses a single resolution parameter $`y_{cut}=10^3`$.

## III Jet Selection

In LO QCD, the fraction of final state jets which are gluons decreases with $`xE_T/\sqrt{s}`$, the momentum fraction of initial state partons within the proton. For fixed $`E_T`$, the gluon jet fraction decreases when $`\sqrt{s}`$ is decreased from 1800 GeV to 630 GeV. We define gluon and quark enriched jet samples with identical cuts in events at $`\sqrt{s}=1800`$ and 630 GeV to reduce experimental biases and systematic effects. Of the two highest $`E_T`$ jets in the event, we select jets with $`55<E_T<100`$ GeV and $`|\eta |<0.5`$.

## IV Quark and Gluon Subjet Multiplicity

There is a simple method to extract a measurement of quark and gluon jets on a statistical basis, using the tools described in the previous sections. $`M`$ is the subjet multiplicity in a mixed sample of quark and gluon jets. It may be written as a linear combination of subjet multiplicity in gluon and quark jets:

$$M=fM_g+(1f)M_q$$
(1)   
The coefficients are the fractions of gluon and quark jets in the sample, $`f`$ and $`(1f)`$, respectively. Consider Eq. (1) for two similar samples of jets at $`\sqrt{s}=1800`$ and 630 GeV, assuming $`M_g`$ and $`M_q`$ are independent of $`\sqrt{s}`$. The solutions are

$$M_q=\frac{f^{1800}M^{630}f^{630}M^{1800}}{f^{1800}f^{630}}$$
(2)   
$$M_g=\frac{\left(1f^{630}\right)M^{1800}\left(1f^{1800}\right)M^{630}}{f^{1800}f^{630}}$$
(3)   
where $`M^{1800}`$ and $`M^{630}`$ are the experimental measurements in the mixed jet samples at $`\sqrt{s}=1800`$ and 630 GeV, and $`f^{1800}`$ and $`f^{630}`$ are the gluon jet fractions in the two samples. The method relies on knowledge of the two gluon jet fractions.

## V Results

The HERWIG 5.9 Monte Carlo event generator provides an estimate of the gluon jet fractions. The method is tested using the detector simulation and CTEQ4M PDF. We tag every selected jet in the detector as either quark or gluon by the identity of the nearer (in $`\eta \times \varphi `$ space) final state parton in the QCD 2-to-2 hard scatter. Fig. 1 shows that gluon jets in the detector simulation have more subjets than quark jets. The tagged subjet multiplicity distributions are similar at the two center of mass energies, verifying the assumptions in § IV.

We count tagged gluon jets and find $`f^{1800}=0.59\pm 0.02`$ and $`f^{630}=0.33\pm 0.03`$, where the uncertainties are estimated from different gluon PDF’s. The nominal gluon jet fractions and the Monte Carlo measurements at $`\sqrt{s}=1800`$ and 630 GeV are used in Eqs. (2-3). The extracted quark and gluon jet distributions in Fig. 1 agree with the tagged distributions and demonstrate closure of the method.

Figure 2 shows the raw subjet multiplicity in DØ data at $`\sqrt{s}=1800`$ GeV is higher than at $`\sqrt{s}=630`$ GeV. This is consistent with the prediction that there are more gluon jets at $`\sqrt{s}=1800`$ GeV compared to $`\sqrt{s}=630`$ GeV, and gluons radiate more than quarks. The combination of the distributions in Fig. 2 and the gluon jet fractions gives the raw subjet multiplicity distributions in quark and gluon jets, according to Eqs. (2-3).

The quark and gluon raw subjet multiplicity distributions need separate corrections for various detector-dependent effects. These are derived from Monte Carlo, which describes the raw DØ data well. Each Monte Carlo jet in the detector simulation is matched (within $`\mathrm{\Delta }<0.5`$) to a jet reconstructed from particles without the detector simulation. We tag detector jets as either quark or gluon, and study the subjet multiplicity in particle jets $`M^{ptcl}`$ vs. that in detector jets $`M^{det}`$. The correction unsmears $`M^{det}`$ to give $`M^{ptcl}`$, in bins of $`M^{det}`$. Figure 3 shows the corrected subjet multiplicity is clearly larger for gluon jets compared to quark jets.

The gluon jet fractions are the largest source of systematic error. We vary the gluon jet fractions by the uncertainties in an anti-correlated fashion at the two values of $`\sqrt{s}`$ to measure the effect on $`R`$. The systematic errors listed in Table I are added in quadrature to obtain the total uncertainty in the corrected ratio $`R=\frac{M_g1}{M_q1}=1.91\pm 0.04(\mathrm{stat})_{0.19}^{+0.23}(\mathrm{sys})`$.

## VI Conclusion

We extract the $`y_{cut}=10^3`$ subjet multiplicity in quark and gluon jets from measurements of mixed jet samples at $`\sqrt{s}=1800`$ and 630 GeV. On a statistical level, gluon jets have more subjets than quark jets. We measure the ratio of additional subjets in gluon jets to quark jets $`R1.9\pm 0.2`$. The ratio is well described by the HERWIG parton shower Monte Carlo, and is only slightly smaller than the naive QCD prediction 9/4.

## VII Acknowledgements

We thank the Fermilab and collaborating institution staffs for contributions to this work and acknowledge support from the Department of Energy and National Science Foundation (USA), Commissariat à L’Energie Atomique (France), Ministry for Science and Technology and Ministry for Atomic Energy (Russia), CAPES and CNPq (Brazil), Departments of Atomic Energy and Science and Education (India), Colciencias (Colombia), CONACyT (Mexico), Ministry of Education and KOSEF (Korea), and CONICET and UBACyT (Argentina).