# 1 Introduction

## 1 Introduction

A sizeable fraction of the final states produced in high energy collisions shows the characteristic feature of large amounts of hadronic energy in small angular regions. These collimated sprays of hadrons (called jets) are the observable signals of underlying short distance processes and are considered to be the footprints of the underlying partonic final states. Quantitative studies of jet production require a precise jet definition, which is given by a jet finding algorithm. 

Jets so defined exhibit an internal structure which is sensitive to the mechanism by which a complex aggregate of observable hadrons evolves from a hard process. The understanding of this mechanism involves higher orders of the strong coupling constant in perturbation theory as well as non-perturbative contributions. This is a challenging task for theory. Recently, for some specific hadronic final state quantities, encouraging results have been obtained by exploiting the characteristic power behaviour of non-perturbative effects and by analytical, approximate calculations of perturbative QCD parton evolution down to the semi-soft regime . Furthermore, since jet production rates are used to test the predictions of perturbative QCD, the understanding of their detailed properties and internal structure is an important prerequisite.

The internal structure of jets has been studied in $`e^+e^{}`$ and in hadron-hadron collisions . At the $`e^\pm p`$ collider HERA, these investigations can be performed in photoproduction ($`Q^20\mathrm{GeV}^2`$) and in deep-inelastic scattering (DIS) at large squared four momentum transfers $`Q^2`$. In a previous publication we have measured the $`E_T`$ dependence of the jet width  in photoproduction. Recently, the ZEUS collaboration has investigated jet shapes in photoproduction  and in DIS at $`Q^2>100`$ $`\mathrm{GeV}^2`$ . Both analyses are carried out in the laboratory frame. This means that for DIS at high $`Q^2`$ mostly events with only one jet enter the analysis.

The hadronization of the current jet in deep-inelastic scattering in the Breit frame has already been studied with event shape variables , charged particle multiplicities and fragmentation functions . In this paper we take the first steps towards a complete understanding of jet properties in DIS. We analyse the hadronization of jets in multijet production in the Breit frame. The Breit frame, where the virtual photon interacts head-on with the proton, has been chosen in this analysis because here the produced transverse<sup>1</sup><sup>1</sup>1transverse with respect to the $`z`$-axis which is given by the axis of the virtual photon and the proton. energy, $`E_{T,\mathrm{Breit}}`$, directly reflects the hardness of the underlying QCD process. We present measurements of internal jet structure in a sample of inclusive dijet events with transverse jet energies of $`E_{T,\mathrm{Breit}}>5`$ GeV, $`10<Q^2120\mathrm{GeV}^2`$ and $`210^4x_{\mathrm{Bj}}810^3`$. This is the $`E_{T,\mathrm{Breit}}`$ range where jet cross section measurements are currently performed at HERA and compared to perturbative QCD calculations (e.g. ). The analysis is based on data taken in 1994 with the H1 detector at HERA when $`27.5`$GeV positrons collided with $`820`$GeV protons. The data correspond to an integrated luminosity of $`_{\mathrm{int}}2\text{pb}^1`$.

Jets are defined in the Breit frame by $`k_{}`$ and cone jet algorithms. Two observables, jet shapes and, for the first time, subjet multiplicities, are studied. The jet shape measures the radial distribution of the transverse jet energy around the jet axis. For the $`k_{}`$ cluster algorithm we have also measured the multiplicity of subjets, resolved at a resolution scale which is a fraction of the jet’s transverse energy. Both observables are presented for different ranges of the transverse jet energy and the pseudo-rapidity<sup>2</sup><sup>2</sup>2The pseudo-rapidity $`\eta `$ is defined as $`\eta \mathrm{ln}(\mathrm{tan}\theta /2)`$ where $`\theta `$ is the polar angle with respect to the proton direction. This definition is chosen in both the laboratory frame and the Breit frame. of the jets in the Breit frame.

The paper is organized as follows. Section 2 gives a brief description of the H1 detector. In section 3 we introduce the jet algorithms used in the analysis and give the definition of the measured observables in section 4. In section 5 we give a short description of the QCD models which are used for the correction of the data and to which the results are later compared (in section 9). The data selection and the correction procedure are described in sections 6 and 7 and the results are discussed in section 8.

## 2 The H1 Detector

A detailed description of the H1 detector can be found elsewhere . Here we briefly introduce the detector components relevant for this analysis: the liquid argon (LAr) calorimeter , the backward lead-scintillator calorimeter (BEMC) , and the tracking chamber system .

The hadronic energy flow is mainly measured by the LAr calorimeter extending over the polar angular range $`4.4^{}<\theta <154^{}`$ with full azimuthal coverage. The polar angle $`\theta `$ is defined with respect to the proton beam direction ($`+z`$ axis). The LAr calorimeter consists of an electromagnetic section ($`2030`$ radiation lengths) with lead absorbers and a hadronic section with steel absorbers. The total depth of both calorimeters varies between $`4.5`$ and $`8`$ interaction lengths. Test beam measurements of the LAr calorimeter modules show an energy resolution of $`\sigma _E/E0.50/\sqrt{E[\text{GeV}]}0.02`$ for charged pions . The absolute scale of the hadronic energy is known for the present data sample to $`4\%`$.

The scattered positron is detected by the BEMC with a depth of $`22.5`$ radiation lengths covering the backward region of the detector, $`155^{}<\theta <176^{}`$. The electromagnetic energy scale is known to an accuracy of $`1\%`$.

The calorimeters are surrounded by a superconducting solenoid providing a uniform magnetic field of $`1.15`$ T parallel to the beam axis in the tracking region.

Charged particle tracks are measured in two concentric jet drift chamber modules (CJC), covering the polar angular range $`15^{}<\theta <165^{}`$. The forward tracking detector covers $`7^{}<\theta <25^{}`$ and consists of drift chambers with alternating planes of parallel wires and others with wires in the radial direction. A backward proportional chamber (BPC) with an angular acceptance of $`151^{}<\theta <174.5^{}`$ improves the identification of the scattered positron. The spatial resolution for reconstructed BPC hits is about 1.5 mm in the plane perpendicular to the beam axis.

## 3 Jet Definitions

The jet algorithms used in this analysis are applied to the particles boosted into the Breit frame. Particle refers here either to an energy deposit in the detector (see section 6), to a stable hadron or a parton in a QCD model calculation. In all cases the scattered positron is excluded. The Breit frame is defined by $`\stackrel{}{q}+2x_{\mathrm{Bj}}\stackrel{}{P}=0`$, where $`\stackrel{}{q}`$ and $`\stackrel{}{P}`$ are the momenta of the exchanged boson and the incoming proton. The $`z`$-axis is defined as the direction of the incoming proton.

In the following analysis we use two different jet definitions: a cone algorithm and a $`k_{}`$ cluster algorithm. Both jet definitions are invariant under boosts along the $`z`$-direction. The recombination of particles is carried out in the $`E_T`$ recombination scheme, which is based on transverse energies $`E_T`$, pseudo-rapidities $`\eta `$ and azimuthal angles $`\varphi `$ of the particles. The transverse energy and the direction of a jet are defined by

$$E_{T,\mathrm{jet}}=\underset{i}{}E_{T,i},\eta _{\mathrm{jet}}=\frac{_iE_{T,i}\eta _i}{_iE_{T,i}},\varphi _{\mathrm{jet}}=\frac{_iE_{T,i}\varphi _i}{_iE_{T,i}},$$
(1)   
where the sums run over all particles $`i`$ assigned to the jet<sup>3</sup><sup>3</sup>3All particles are considered massless by setting $`E_i=|\stackrel{}{p_i}|`$..

### 3.1 Cone Algorithm

Based on the original proposal of Sterman and Weinberg  many different implementations of cone algorithms have been developed. While the basic idea of the cone algorithm is simple and very intuitive, an operational definition is non-trivial. The resulting jet cross sections depend on how the algorithm treats the choice of jet initiators and configurations of overlapping jet cones. It has repeatedly been pointed out that many definitions of cone algorithms are not infrared and/or collinear safe .

In this analysis we use the definition implemented in the algorithm PXCONE  which does not suffer from the problems discussed in . This definition, which corresponds closely to the Snowmass proposal  and to the algorithm used in the CDF experiment , is also used by the OPAL collaboration .

Particles are assigned to jets based on their spatial distance $`R`$ in pseudo-rapidity and azimuth space ($`R^2=\mathrm{\Delta }\eta ^2+\mathrm{\Delta }\varphi ^2`$). The algorithm operates as follows:

1. Each particle is considered as a seed of a jet, for which steps 2-4 are performed.
2. The jet momentum is calculated from all particles within a cone of radius $`R_0`$ around the seed direction using eq. (1).
3. If the jet direction differs from the seed direction, the jet direction is taken as the new seed direction and step 2 is repeated.
4. When the jet direction is stable the jet is stored in the list of “protojets” (if it is not identical with a protojet already found).
5. The steps 2 to 4 are repeated for all midpoints of pairs of protojets as seed directions<sup>4</sup><sup>4</sup>4In practice it is sufficient to do this only for pairs of protojets with a distance between $`R_0`$ and $`2R_0`$.. This leads to the infrared safety of the procedure .
6. Protojets with transverse energies of $`E_{T,\mathrm{jet}}<ϵ`$ are removed from the list. The cut-off parameter $`ϵ`$ specifies below which transverse energies protojets are not considered in the overlap treatment (steps 7-8).
7. All remaining protojets that have more than a fraction $`f`$ of their transverse energy contained in a protojet of higher transverse energy are deleted.
8. All particles that are contained in more than one protojet are assigned to the protojet whose center is nearest in $`(\eta ,\varphi )`$.
9. The jet momenta are recalculated using eq. (1). All protojets with $`E_{T,\mathrm{jet}}<ϵ`$ are deleted and the remaining ones are called jets.

The jets with the highest transverse energies are considered in the analysis. Due to the reassignment of particles to jets and the recalculation of the jet axis (steps 7, 8) it may happen that single particles within a jet have a distance larger than $`R_0`$ to the jet axis. This analysis is made with the parameter settings $`ϵ=5\text{GeV}`$, $`f=0.75`$ and a cone radius of $`R_0=1.0`$.

### 3.2 Inclusive $`𝒌_{\mathbf{}}`$ Algorithm

The ambiguities that occur for cone jet definitions (choice of seeds, overlapping cones) are avoided in cluster algorithms which successively recombine particles to jets. One definition of such an algorithm (proposed in  and implemented in the KTCLUS algorithm ) has properties very similar to cone algorithms. As in the cone algorithm the clustering procedure is based on the longitudinally boost-invariant quantities $`E_T,\mathrm{\Delta }\eta ,\mathrm{\Delta }\varphi `$. The minimum of all distances between particles is determined and either the corresponding pairs of particles are merged into pseudo-particles or single (pseudo-) particles are declared as jets. This process is iterated until no particles are left:

1. We start with a list of all particles and an empty list of jets.
2. For each particle $`i`$ as well as for each pair of particles ($`i,j`$) the distances $`d_i`$ and $`d_{ij}`$ are calculated

$$d_i=E_{T,i}^2R_0^2\text{and}\text{ }d_{ij}=\mathrm{min}(E_{T,i}^2,E_{T,j}^2)R_{ij}^2\mathrm{with}R_{ij}^2=\mathrm{\Delta }\eta _{ij}^2+\mathrm{\Delta }\varphi _{ij}^2.$$
(2)
3. The smallest value of all the $`d_i`$ and $`d_{ij}`$ is labeled $`d_{\mathrm{min}}.`$
4. If $`d_{\mathrm{min}}`$ belongs to the set of $`d_{ij}`$, the particles $`i`$ and $`j`$ are merged into a new particle using the recombination prescription in eq. (1) and removed from the list of particles.
5. If $`d_{\mathrm{min}}`$ belongs to the set of $`d_i`$, the particle $`i`$ is removed from the list of particles and added to the list of jets.
6. When no particles are left (i.e. all particles are included in jets) the procedure is finished.

The last jets that entered the list are the ones with highest transverse energies. These jets are considered in the analysis. This jet definition implies that particles with $`R_{ij}<R_0`$ are subsequently merged, so that all final jets are separated by distances $`R_{ij}>R_0`$. It is still possible that particles inside a jet have a distance $`R_{ij}>R_0`$ to the jet axis and that particles with $`R_{ij}<R_0`$ are not part of the jet. The parameter $`R_0`$ is set to $`R_0=1.0`$.

## 4 The Observables

Two observables of internal jet structure are investigated in this analysis. They are sensitive to different aspects of jet broadening.

The jet shapes are studied for the cone and the $`k_{}`$ algorithm. This observable measures the radial distribution of the transverse jet energy only and is affected by hard and by soft processes over the whole radial range.

A natural choice for studying the internal structure of jets with the $`k_{}`$ cluster algorithm is the multiplicity of subjets, resolved at a resolution scale which is a fraction of the jet’s transverse energy. These subjet multiplicities are sensitive to more local structures of relative transverse momentum within a jet. Here the perturbative and the non-perturbative contributions are better separated. While at larger values of the resolution parameter perturbative contributions dominate, at smaller values non-perturbative contributions become increasingly important.

### 4.1 The Jet Shape

The jet shape $`\mathrm{\Psi }(r)`$ is defined as the fractional transverse jet energy contained in a subcone of radius $`r`$ concentric with the jet axis, averaged over all considered jets in the event sample

$$\mathrm{\Psi }(r)\frac{1}{N_{\mathrm{jets}}}\underset{\mathrm{jets}}{}\frac{E_T(r)}{E_{T,\mathrm{jet}}},$$
(3)   
where $`N_{\mathrm{jets}}`$ is the total number of these jets. As proposed in , only particles assigned by the jet algorithm to the jet are considered.

Usually the denominator in the definition of $`\mathrm{\Psi }`$ is given by the summed $`E_T`$ of all particles within a radius $`R_0`$ to the jet axis. This means that $`\mathrm{\Psi }(r/R_0=1)=1`$. In our definition (3) of $`\mathrm{\Psi }`$ the denominator is given by the transverse energy of the jet. Since neither for the cone nor for the $`k_{}`$ definition are all particles necessarily assigned to a jet within a radius of $`r/R_0<1`$ to the jet axis, $`\mathrm{\Psi }(r/R_0=1)`$ is not constrained to have the value of one. With this choice of our observable we are also sensitive to the amount of transverse jet energy outside the radius $`R_0`$.

### 4.2 Subjet Multiplicities

For each jet in the sample the clustering procedure is repeated for all particles assigned to the jet. The clustering is stopped when the distances $`y_{ij}`$ between all particles $`i,j`$ are above some cut-off $`y_{\mathrm{cut}}`$

$$y_{ij}=\frac{\mathrm{min}(E_{T,i}^2,E_{T,j}^2)}{E_{T,\mathrm{jet}}^2}\frac{\left(\mathrm{\Delta }\eta _{ij}^2+\mathrm{\Delta }\varphi _{ij}^2\right)}{R_0^2}>y_{\mathrm{cut}}$$
(4)   
and the remaining (pseudo-)particles are called subjets. The parameter $`y_{\mathrm{cut}}`$ defines the minimal relative transverse energy between subjets inside the jet and thus determines the extent to which the internal jet structure is resolved. From this definition it follows that for $`y_{\mathrm{cut}}>0.25`$ no subjet is resolved (therefore the number of subjets is one), while for $`y_{\mathrm{cut}}0`$ every particle in the jet is a subjet. The observable that is studied in this analysis is the average number of subjets for a given value of the resolution parameter, for values $`y_{\mathrm{cut}}10^3`$.

## 5 QCD Models

A simulation of the detailed properties of the hadronic final state is available in the form of Monte Carlo event generators. They include the matrix element of the hard subprocess in first order of the strong coupling constant $`\alpha _s`$, approximations of higher order QCD radiation effects, and a model to describe the non-perturbative transition from partons to hadrons.

The LEPTO Monte Carlo  incorporates the $`𝒪(\alpha _s)`$ QCD matrix element and takes higher order parton emissions to all orders in $`\alpha _s`$ approximately into account using the concept of parton showers  based on the leading logarithm DGLAP equations . QCD radiation can occur before and after the hard subprocess. The formation of hadrons is performed using the LUND string model  implemented in JETSET .

The HERWIG Monte Carlo  also includes the $`𝒪`$$`(\alpha _s)`$ QCD matrix element, but uses another implementation of the parton shower cascade which takes coherence effects fully into account. The hadronization is simulated with the cluster fragmentation model .

In ARIADNE  gluon emissions are treated by the colour dipole model  assuming a chain of independently radiating dipoles spanned by colour connected partons. The first emission in the cascade is corrected to reproduce the matrix element to first order in $`\alpha _s`$ .

DJANGO  provides an interface between the event generators LEPTO or ARIADNE and HERACLES  which makes it possible to include $`𝒪(\alpha )`$ QED corrections at the lepton line.

## 6 Data Selection

The analysis is based on H1 data taken in $`1994`$ corresponding to an integrated luminosity of $`_{int}2\text{pb}^1`$. The event selection closely follows that described in a previous publication . DIS events are selected where the scattered positron is measured in the acceptance region of the BEMC at energies where trigger efficiencies are approximately 100 %. To ensure a good identification of the scattered positron and to suppress background from misidentified photoproduction events the following cuts are applied:

* The cluster of the positron candidate must have an energy-weighted mean transverse radius below $`5\text{cm}`$.
* A reconstructed BPC hit within $`5\text{cm}`$ of the straight line connecting the shower center with the event vertex is required.
* The $`z`$ position of the reconstructed event vertex must be within $`\pm 30\text{cm}`$ of the nominal position.
* A cut on $`35\text{GeV}<(Ep_z)<70\text{GeV}`$ is applied, where the sum runs over all energy deposits in the calorimeter. In neutral current DIS events without undetected photon radiation the quantity $`(Ep_z)`$ is expected to be equal to twice the energy of the initial state positron. This cut reduces the contribution from photoproduction events as well as events where hard photons are radiated collinear to the incoming positron.

The event kinematics are calculated from the polar angle $`\theta _{el}`$ and the energy $`E_{el}^{}`$ of the scattered positron via $`Q_{el}^2=2E_0E_{el}^{}(1+\mathrm{cos}\theta _{el})`$, $`y_{el}=1E_{el}^{}/(2E_0)(1\mathrm{cos}\theta _{el})`$ and $`x_{\mathrm{Bj}}=Q^2/(sy)`$. $`E_0`$ denotes the energy of the incoming positron and $`s`$ the $`ep`$ centre-of-mass energy squared. Events are only accepted, if $`E_{el}^{}>11`$ GeV, $`156^{}<\theta _{el}<173^{}`$, $`Q^2>10\mathrm{GeV}^2`$ and $`y>0.15`$. The resulting kinematic range is $`10<Q^2120\mathrm{GeV}^2`$ and $`210^4x_{\mathrm{Bj}}810^3`$.

Jets are defined by the algorithms described in section 3. The input for the jet algorithms consists of a combination of energy clusters from the calorimeter and track momenta measured in the central and forward trackers (as described in ). While all energy clusters are considered, the four momentum of each single track is only allowed to contribute up to a momentum of $`350\text{MeV}`$. This procedure partly compensates for energy losses in the calorimeter due to dead material and noise thresholds. It reduces the dependence of the jet finding efficiency on the pseudo-rapidity of the jet and improves the reconstruction of the transverse jet energy .

The objects from tracking and calorimeter information are boosted to the Breit frame where the jet algorithms are applied. We select events with at least two identified jets with transverse energies of $`E_{T,\mathrm{Breit}}>5\text{GeV}`$ in $`1<\eta _{\text{jet,lab}}<2`$. The two jets with the highest $`E_{T,\mathrm{Breit}}`$ are considered in the analysis. The event sample for the inclusive $`k_{}`$ algorithm (the cone algorithm) consists of 2045 (2657) dijet events.

## 7 Correction of the Data

The data are corrected for detector effects and QED radiation from the lepton. The detector response is determined using events from Monte Carlo event generators that were subjected to a detailed simulation of the H1 detector. The following event generators are used: ARIADNE interfaced in DJANGO (with and without the inclusion of QED corrections) and LEPTO. Both generators give a good description of the kinematic variables of the inclusive DIS data sample as well as of the angular and transverse energy distributions of the jets . We also observe a reasonable description of the observables introduced in section 4 (see section 9).

The measured data points are corrected bin-by-bin for detector effects. Using the generated event samples, the correction factor for each bin is determined as the ratio of the generated value of the observable and the value that is reconstructed after detector simulation. These correction factors are independent of the inclusion of QED radiation effects as included in DJANGO. Their dependence on details of the modeling of the hadronic final state is taken into account by considering the difference between the correction factors from ARIADNE and LEPTO as systematic uncertainty.

For the $`k_{}`$ (cone) algorithm the corrections for $`\mathrm{\Psi }(r)`$ are below $`10\%`$ ($`13\%`$) for subcone radii $`r>0.3`$ and always below $`27\%`$ ($`23\%`$). The corrections for $`N_{\mathrm{subjet}}(y_{\mathrm{cut}})`$ are always below $`7\%`$. The correction factors from both QCD models are in good agreement (they differ typically by not more than $`2\%`$) for the jet shapes as well as for the subjet multiplicities . The final correction factors are taken to be the mean values of the two models, taking the spread as the error. In addition we have varied the calibration of the hadronic energy scale in the data sample in the range of $`\pm 4\%`$ around the nominal value. The error is estimated as the maximal deviation from the results at the nominal value. For all observables it is at most $`5\%`$. The overall systematic error is calculated by adding the errors from the model dependence and from the uncertainty of the hadronic energy scale in quadrature. In all figures the statistical and systematic errors are added in quadrature. Since each jet enters in all bins of a distribution, all errors are correlated.

The background from misidentified photoproduction events is estimated with a sample of photoproduction events generated with PHOJET  and is found to be negligible.

## 8 Results

The jet shape and the subjet multiplicity are presented as functions of quantities directly related to the single jets, namely the transverse jet energy ($`E_{T,\mathrm{Breit}}`$) and the pseudo-rapidity ($`\eta _{\text{Breit}}`$) in the Breit frame. We also investigated whether the observables depend on the event kinematics. The jet shapes and subjet multiplicities were compared for two bins of $`Q^2`$ ($`Q^2<20\mathrm{GeV}^2`$ and $`Q^2>20\mathrm{GeV}^2`$) and $`x_{\mathrm{Bj}}`$ ($`x_{\mathrm{Bj}}<810^4`$ and $`x_{\mathrm{Bj}}>810^4`$) respectively. No dependence on $`Q^2`$ and $`x_{\mathrm{Bj}}`$ has been observed.

### 8.1 Jet Shapes

The radial dependence of the jet shape $`\mathrm{\Psi }(r)`$ for the $`k_{}`$ algorithm is shown in Fig. 1 in different ranges of the pseudo-rapidity in the Breit frame. The results for jets of transverse energies $`5<E_{T,\mathrm{Breit}}<8\mathrm{GeV}`$ and $`E_{T,\mathrm{Breit}}>8\mathrm{GeV}`$ are superimposed. The jet shape $`\mathrm{\Psi }(r)`$ increases faster with $`r`$ for jets at larger transverse energies, indicating that these jets are more collimated. The same tendency is seen for the jets defined by the cone algorithm which are compared to the jets found by the $`k_{}`$ algorithm in Fig. 2. For both jet definitions we also observe a dependence of the jet shape on the pseudo-rapidity of the jets. Jets towards the proton direction (at larger values of $`\eta _{\text{Breit}}`$) are broader than jets towards the photon direction (smaller $`\eta _{\text{Breit}}`$). In the region where the jets are most collimated ($`E_{T,\mathrm{Breit}}>8\mathrm{GeV}`$ and $`\eta _{\text{Breit}}<2.2`$), very similar jet shapes are observed for the $`k_{}`$ and cone algorithms. The broadening of the jets for smaller $`E_{T,\mathrm{Breit}}`$ and larger $`\eta _{\text{Breit}}`$ is more pronounced for the cone jet definition.

Recently jet shapes have been measured in dijet production in photon-photon collisions  for jets defined by a cone algorithm at transverse energies comparable to those presented here. The jet shapes in photon-photon collisions (where no $`\eta `$ dependence is observed) are very similar to those measured in DIS in the Breit frame at $`\eta _{\text{Breit}}<1.5`$.

### 8.2 Subjet Multiplicities

The subjet multiplicities for the $`k_{}`$ algorithm are displayed in Fig. 3. The average number of subjets $`N_{\mathrm{subjet}}(y_{\mathrm{cut}})`$ as a function of the subjet resolution parameter at $`y_{\mathrm{cut}}10^3`$ is plotted. Towards smaller values of $`y_{\mathrm{cut}}`$, an increasing number of jet fragments with smaller relative transverse momenta is resolved. The number of subjets at a given value of $`y_{\mathrm{cut}}`$ reflects the amount of relative transverse momentum with respect to the jet axis. The subjet multiplicity is therefore a measure of the broadness of the jet.

At $`y_{\mathrm{cut}}=10^3`$ a jet is on average resolved into $`4.1`$$`4.6`$ subjets, depending on $`E_{T,\mathrm{Breit}}`$ and $`\eta _{\text{Breit}}`$ of the jet<sup>5</sup><sup>5</sup>5On average the jets in the data (as in the simulated events) consist of eleven calorimetric energy clusters. For the LEPTO generator this is also approximately the average multiplicity of stable particles inside the jets.. For almost all values of $`y_{\mathrm{cut}}`$ the subjet multiplicity is larger for jets at smaller $`E_{T,\mathrm{Breit}}`$ and larger $`\eta _{\text{Breit}}`$, indicating broader jets.

A summary of the results for both observables is given in Fig. 4. Here the $`E_{T,\mathrm{Breit}}`$ and the $`\eta _{\text{Breit}}`$ dependence of the jet shape and the average number of subjets are shown at an intermediate value of the resolution parameter (jet shape: $`r=0.5`$ and subjet multiplicity: $`y_{\mathrm{cut}}=10^2`$).

Although the subjet multiplicities are sensitive to the jet broadening in a different way than the jet shapes, consistent conclusions can be drawn for both measurements. The jet broadening depends on both the transverse jet energy as well as the pseudo-rapidity in the Breit frame. While the pseudo-rapidity dependence is most pronounced at smaller transverse jet energy, the transverse energy dependence is stronger in the forward region (at larger pseudo-rapidities).

## 9 Comparison with QCD Model Predictions

The predictions of different QCD models are compared in Fig. 5 to the jet shapes measured for the $`k_{}`$ algorithm. The models LEPTO, ARIADNE and HERWIG all show $`E_{T,\mathrm{Breit}}`$ and $`\eta _{\text{Breit}}`$ dependences similar to that seen in the data. LEPTO gives the best description of the measured shapes for $`\eta _{\text{Breit}}<2.2`$ while at $`\eta _{\text{Breit}}>2.2`$ the predicted jet shapes are too broad. A reasonable description is also obtained by the ARIADNE model except for jets at smaller pseudo-rapidities where the jet shapes have the tendency to be too narrow. For the HERWIG model the jet shapes are narrower than those in the data in all $`E_{T,\mathrm{Breit}}`$ and $`\eta _{\text{Breit}}`$ regions. The same observations as above are made when comparing these QCD models with the subjet multiplicities and with the jet shapes for the cone algorithm (not shown here).

In QCD models the evolution of a jet is described by perturbative contributions (radiation of partons) and non-perturbative contributions (hadronization). Studies based on the LEPTO and HERWIG parton shower models show that all observables studied in this analysis are strongly influenced by hadronization. This process has the largest impact on the jet broadening in our kinematic region (Fig. 6). Basic characteristics of the perturbative contributions are however still visible after hadronization. The model prediction suggests that the large difference between quark and gluon-initiated jets before hadronization survives the hadronization process. This especially applies to jets with large transverse energies .

Fig. 6 shows the jet shapes and the subjet multiplicities as predicted by the LEPTO parton shower model for the $`k_{}`$ algorithm, separately for quark and gluon jets at $`E_{T,\mathrm{Breit}}>8\mathrm{GeV}`$ and $`\eta _{\text{Breit}}<1.5`$. Gluon jets are broader than quark jets. The same prediction is obtained by the HERWIG parton shower model. Although the jets in HERWIG are slightly narrower, the differences between gluon and quark jets are equally large. In the phase space considered here, LEPTO and HERWIG (in agreement with next-to-leading order calculations) predict a fraction of approximately $`80\%`$ photon-gluon fusion events with two quarks in the partonic final state. The jet samples of these models are therefore dominated by quark jets. Both model predictions for the jet shapes and the subjet multiplicities therefore mainly reflect the properties of the quark jets as can be seen in Fig. 6. These predictions give a reasonable description of the data. Thus, we conclude, that the jets we observe are consistent with being mainly initiated by quarks.

## 10 Summary

Measurements of internal jet structure in dijet events in deep-inelastic scattering in the kinematic domain $`10<Q^2120\mathrm{GeV}^2`$ and $`210^4x_{\mathrm{Bj}}810^3`$ have been presented. Jet shapes and subjet multiplicities have been studied for jets of transverse energies $`E_{T,\mathrm{Breit}}>5\mathrm{GeV}`$ defined by $`k_{}`$ and cone jet algorithms in the Breit frame.

The radial dependence of the jet shape and the dependence of the average number of subjets on the subjet resolution parameter $`y_{\mathrm{cut}}`$ are both sensitive to different aspects of jet broadening. For both observables a dependence of the jet broadness on the transverse energy $`E_{T,\mathrm{Breit}}`$ and on the pseudo-rapidity in the Breit frame $`\eta _{\text{Breit}}`$ is seen. With increasing $`E_{T,\mathrm{Breit}}`$ jets are narrower. Jets of the same $`E_{T,\mathrm{Breit}}`$ become broader towards the proton direction. This effect is more pronounced at lower $`E_{T,\mathrm{Breit}}`$.

At lower $`E_{T,\mathrm{Breit}}`$ jets defined by the $`k_{}`$ algorithm are more collimated than jets defined by the cone algorithm, while at higher $`E_{T,\mathrm{Breit}}`$ both algorithms produce very similar jets.

The QCD models LEPTO, ARIADNE and HERWIG roughly reproduce the dependence of the jet shape and the subjet multiplicities on $`E_{T,\mathrm{Breit}}`$ and $`\eta _{\text{Breit}}`$ as seen in the data. LEPTO has a tendency to produce broader jets in the proton direction than measured. HERWIG and ARIADNE produce jets which are too collimated especially at higher transverse energies. We have reported earlier that in the same kinematic domain the predicted jet rates from LEPTO and HERWIG are about a factor of two below the data . Since these models are able to reproduce the internal jet structure, this failure must be largely connected to an inadequate modeling of the underlying hard partonic subprocess.

According to the parton shower models LEPTO and HERWIG, quark and gluon initiated jets differ both at the parton and at the hadron level. Both models predict that the jet sample is dominated by quark initiated jets. Since these models describe our data, we conclude that the observed jet structures are compatible with those of quark initiated jets.

## 11 Acknowledgments

We are grateful to the HERA machine group whose outstanding efforts have made and continue to make this experiment possible. We thank the engineers and technicians for their work in constructing and now maintaining the H1 detector, our funding agencies for financial support, the DESY technical stuff for continual assistance, and the DESY directorate for the hospitality which they extend to the non-DESY members of the collaboration.