# MLLA Parton Spectra Compared to ARIADNE

## 1 Introduction

The perturbative QCD approach to describing the inclusive energy spectra, via the modified leading log approximation (MLLA) in conjunction with local parton hadron duality (LPHD), has been very successful in both $`e^+e^{}`$ annihilation and deep inelastic scattering experiments . Using LPHD, the non-perturbative effects of such distributions are reduced to a simple factor of normalisation that relates the hadronic distributions to the partonic ones. Perturbative features of these distributions are calculated by MLLA which accounts for both the double and single logarithmic effects. The MLLA approach has two free parameters: a running strong coupling, governed by a QCD scale $`\mathrm{\Lambda },`$ and an energy cut-off, $`Q_0,`$ below which the parton evolution is truncated.

The MLLA evolution equation allow the parton spectra for the logarithmic scaled energy spectra, $`\xi ,`$ to be calculated . The variable $`\xi `$ is defined as $`\mathrm{ln}(E_0/E)\mathrm{ln}(1/x_p),`$ where $`E_0`$ is the original energy of the jet and $`E`$ is the parton’s energy. The cut-off, $`Q_0,`$ bounds the parton energy, $`Ek_TQ_0,`$ where $`k_T`$ is the transverse energy of the decay products in the jet evolution. In order to reconstruct the $`\xi `$ distributions one has to perform the inverse Mellin transformation:

$$\overline{D}(\xi ,Y,\lambda )=_{ϵı\mathrm{}}^{ϵ+ı\mathrm{}}\frac{d\omega }{2\pi ı}x_p^\omega D(\omega ,Y,\lambda )$$
(1)   
where the integral runs parallel to the imaginary axis on the right of all singularities in the complex $`\omega `$plane, $`Y=\mathrm{ln}(E_0/Q_0)`$ and $`\lambda =\mathrm{ln}(Q_0/\mathrm{\Lambda }).`$

The Mellin-transformed distributions, $`D(\omega ,Y,\lambda ),`$ can be expressed  in terms of confluent hypergeometric functions, $`\mathrm{\Phi }`$:

$$\begin{array}{cc}D(\omega ,Y,\lambda )\hfill & =\frac{t_1A}{B(B+1)}\mathrm{\Phi }(A+B+1,B+2;t_1)\mathrm{\Phi }(AB,1B;t_2)\\ & +\left(\frac{t_2}{t_1}\right)^B\mathrm{\Phi }(A,B;t_1)\mathrm{\Phi }(A,B+1;t_2),\end{array}$$
(2)   
where

$$\begin{array}{cc}t_1=\omega (Y+\lambda ),\hfill & \hfill t_2=\omega \lambda .\end{array}$$
(3)   
In addition $`A`$ and $`B`$ are defined as:

$$\begin{array}{cc}A=4N_c/b\omega ,\hfill & \hfill B=a/b,\end{array}$$
(4)   
where $`N_c`$ is the number of colours, $`a=11N_c/3+2n_f/3N_c^2`$ , $`n_f`$ is the number of flavours and $`b=11N_c/3+2n_f/3.`$

Equation 1 is then calculated using a numerical integration in the complex $`\omega `$plane.

The current region in the $`ep`$ Breit frame is analogous to a single hemisphere of $`e^+e^{}`$ annihilation. In $`e^+e^{}q\overline{q}`$ annihilation the two quarks are produced with equal and opposite momenta, $`\pm \sqrt{s}/2,`$ where $`\sqrt{s}`$ is the positron-electron centre of mass energy. The fragmentation of these quarks can be compared to that of the quark struck from the proton; this quark has an outgoing momentum $`Q/2`$ in the Breit frame, where $`Q^2`$ is the negative square of the four-momentum of the virtual exchanged boson in DIS. In the direction of this struck quark the scaled momentum spectra of the particles are expected by MLLA  to have a dependence on $`Q`$ similar to that observed in $`e^+e^{}`$ annihilation at energy $`\sqrt{s}=Q,`$ with no Bjorken-$`x`$ dependence.

The ARIADNE Monte Carlo generator  is based on the colour dipole model, CDM . In the CDM, all gluon emissions constituting the QCD cascade start as radiation from the colour dipole formed between the quark and the anti-quark in the case of $`e^+e^{}`$ annihilation or the struck quark and the proton remnant in the case of DIS. All subsequent radiation arises from independent colour dipoles formed either from $`q\overline{q}`$ pairs or softer gluons radiated by the previously produced gluons. In the DIS scenario, the proton remnant is treated as an extended object which results in a suppression of radiation , generally in the proton direction. In addition the struck quark is treated as extended, as the photon only probes it to a distance inversely proportional to the transferred momentum. Treating the remnant and the struck quark as extended objects, rather than point like, results in a reduction in the available phase space for gluon radiation in DIS.

The QCD cascade in ARIADNE is governed by a number of parameters in the Monte Carlo models. Two of the most important are the QCD scale, $`\mathrm{\Lambda },`$ (PARA(1)) and the parameter that determines the $`k_T`$ cut off for the shower (PARA(3)). An additional parameter, PARA(28), also allows the user to bound the lower energy of the emitted parton as well. For this study ARIADNE version 4.10 has been used.

## 2 Comparisons with MLLA

Before investigating the evolution of the shower in DIS, the evolution in the simpler case of $`q\overline{q}`$ pair production in $`e^+e^{}`$ annihilation was studied. The spectra for both MLLA and ARIADNE were generated with a $`\mathrm{\Lambda }=150`$ MeV and a cut-off $`Q_0=2\mathrm{\Lambda }=\mathrm{PARA}(3)=\mathrm{PARA}(28).`$ Below $`\xi 1`$ there are instabilities in the numerical integration of equation 1 so all subsequent comparisons are for $`\xi >1.`$ Except for an overall normalisation discrepancy (a factor of 1.4 greater parton multiplicity in ARIADNE) the $`\xi `$ spectra are in very good agreement as illustrated in Figure 1. This normalisation discrepancy is constant, independent of the $`\sqrt{s}`$ at which the events were generated at. There is a slight tendency for the MLLA calculation to fall off quicker at large values of $`\xi `$ than the ARIADNE predictions.

Using LEPTO  to generate the electroweak cross section and colour flow configuration for DIS, ARIADNE was then used to generate the subsequent QCD cascade. The event was boosted to the Breit frame and those partons in the current fragmentation region selected. The DIS events were generated with fixed kinematics that are accessible in the HERA regime. The corresponding value of Bjorken-$`x`$ with $`Q`$ are shown in Table 1. Using the same values of $`\mathrm{\Lambda }`$ and $`Q_0`$ that was used for the $`e^+e^{}`$ annihilation study, the MLLA prediction is again compared to the ARIADNE generated spectra.

Figure 2 shows the default version of ARIADNE for DIS compared to the MLLA predictions. As $`Q^2`$ increases the discrepancy between ARIADNE and the MLLA calculations becomes more pronounced. The $`\xi `$ distribution of ARIADNE peaks at higher values than the MLLA calculation. In addition, the MLLA calculations are narrower than the ARIADNE predictions. Again the parton multiplicity of the two distributions are different. Unlike the $`e^+e^{}`$ situatuion this normalisation factor seems to exhibit a $`Q`$ dependence. At low $`Q`$ the height of the peak for ARIADNE compared to MLLA is a factor of 1.1 higher whilst in the highest $`Q`$ bin it is 1.1 lower.

In the default ARIADNE, the mechanism for soft suppression of radiation due to the extended source of the proton remnant results in a suppression of radiation in the current region of the Breit frame at high $`Q^2.`$ Figure 3 shows the high $`Q^2`$ modified version of ARIADNE  for DIS, where this suppression in the current region is removed, compared to the MLLA predictions. As expected, this modification to ARIADNE leads to a much better agreement between the MLLA calculations and ARIADNE. The situation with the parton multiplcity is similar. The $`Q`$dependence of the ratio of the peak heights is less than the default ARIADNE, with ARIADNE being a factor $`1.21.3`$ higher.

In both options of the ARIADNE program there are discrepancies evident in the lower $`(x,Q)`$ bins compared to the MLLA predictions. One possible explaination of this discrepancy is given in Ref. , where it shown that high $`p_T`$ emissions in DIS can lead to the situation where the current region of the Breit frame is depopulated.

## 3 Conclusions

The QCD cascade as implemented in the ARIADNE Monte Carlo program is in good agreement with the shape of the MLLA prediction in the simple scenario of $`q\overline{q}`$ production in $`e^+e^{}`$ annihilation. In the more complex situation of DIS the agreement is not as good, unless account is taken of the additional suppression introduced into the model in the current fragmentation region caused by the suppression of phase space due to the extended nature of the proton remnant.