Systematic Analysis of Scaling Properties in Deep Inelastic Scattering
Abstract
Using the “Quality Factor” (QF) method, we analyse the scaling properties of deepinelastic processes at HERA and fixed target experiments for We look for scaling formulae of the form where is a scaling variable suggested by the asymptotic properties of QCD evolution equations with rapidity . We consider four cases: “Fixed Coupling”, corresponding to the original geometric scaling proposal and motivated by the asymptotic properties of the BalitskyKovchegov (BK) equation with fixed QCD coupling constant, two versions “Running Coupling I,II” of the scaling suggested by the BK equation with running coupling, and “Diffusive Scaling” suggested by the QCD evolution equation with Pomeron loops. The Quality Factors, quantifying the phenomenological validity of the candidate scaling variables, are fitted on the total and DVCS crosssection data from HERA and predictions are made for the elastic vectormeson and for the diffractive crosssections at fixed small or The first three scaling formulae have comparably good QF while the fourth one is disfavored. Adjusting initial conditions gives a significant improvement of the “Running Coupling II” scaling.
I Introduction
Geometric scaling Stasto:2000er () is a remarkable empirical property verified by data on high energy deep inelastic scattering (DIS) virtual photonproton crosssections. It was realized that one can represent with reasonable accuracy the cross section by the formula where is the virtuality of the photon, the total rapidity in the proton system and
(1) 
is the scaling variable. The value confirms the value found within the GolecBiernat and Wüsthoff model GolecBiernat:1998js () where geometric scaling was explicitely used for the parametrization.
The scaling using the variable defined in Formula (1) was intimately related to the concept of saturation saturation (), i.e. the behaviour of perturbative QCD amplitudes when the density of partons becomes high enough to exercise nonlinear effects ensuring the unitarity bound. Indeed, there were many theoretical arguments to infer that in a domain in and where saturation effects set in, geometric scaling is expected to occur. Within this framework, the function can be called the saturation scale, since it delineates the approximate lower bound of the saturation domain.
This type of geometric scaling is motivated by asymptotic properties of QCD evolution equations with rapidity. Using the nonlinear BalitskyKovchegov (BK) equation Balitsky:1995ub () which represents the “meanfield” approximation of high energy (or high density) QCD, geometric scaling could be derived from its asymptotic solutions Munier:2003vc (). This equation is supposed to capture some essential features of saturation effects. Considering the BK equation with coupling constant leads asymptotically to the original geometric scaling of Formula (1). Considering a coupling leads to the following scaling:
(2) 
Recently gb (), it was noticed that the scaling solution (2) of the BK equation with running coupling is only approximate and not unique. Another equivalent approximation leads to a different scaling variable, namely
(3) 
The key ingredient to theoretically prove geometric scaling of the asymptotic solutions of the nonlinear BK equation is the traveling wave method Munier:2003vc (). Indeed, the BK equation admits solutions in the form of traveling waves . has the interpretation of a space variable while , interpreted as time, is an increasing function of rapidity namely for the fixed coupling case and for the running one. is the critical velocity of the wave. These properties confirmed results previously obtained Mueller:2002zm (), for instance by replacing the nonlinear damping term by absorbing boundary conditions, and thus considering only the linear part of the BK equation, which is equivalent to the Balitsky Fadin Kuraev Lipatov (BFKL) equation Lipatov:1976zz ().
The subsequent elaboration on QCD evolution equations led to go beyond the meanfield approximation. The effect of fluctuations was examined in Ref. Hatta:2006hs () in the fixed coupling scheme and gives rise to a new “diffusive scaling”, the scaling variable being
(4) 
The aim of this paper from a theoretical pointofview is to test and compare the different scaling behaviors, arising from different versions of QCD evolution, using the data available from HERA. We shall study the phenomenological relevance of the four kinds of scaling and refer to them in the following as “Fixed Coupling” for the variable (1), “Running Coupling I” and “Running Coupling II” for the variables (2) and (3) respectively, and to “Diffusive Scaling” for (4). Our method can be easily applied to any new proposal of scaling.
From an experimental pointof view, one wants to include not only the newest published data on from HERA, but also from other processes where the small physics allow one to discuss scaling properties. We thus include the data on DVCS cross sections and, following the geometric scaling analysis performed in Ref. Marquet:2006jb (), we extend the analysis to elastic vectormeson production and to diffractive crosssections at fixed and small Moreover, and as a new aspect of scaling analysis, we are led to also include diffractive crosssections at fixed and small since scaling should also be considered for the cross section, following the IngelmanSchlein interpretation of diffractive structure functions Ingelman:1984ns ().
In order to define how good the scalings are, the authors of Ref. usgeom () introduce a (QF) as an estimator on the validity of scaling. The main property of the QFs is that it does not depend on a given parametrization of the scaling curve as in an ordinary approach. The aim of this paper is thus to extend on the same footing the QF studies to all scaling laws, including the new scaling Running Coupling II, and to compare them. We perform a systematic analysis using the QF techique investigating the scaling properties in deepinelastic scattering using all data available. It will include both the largest available set of published data where scaling is expected to occur and the new status on theoretical scaling properties using QCD evolution equations.
The plan of the paper is as follows. In section II we describe the theoretical motivation and the give the precise formulation of the various scaling hypotheses. In section III we present the Quality Factor method we shall use in the analysis and the set of data we consider for the fitting procedure. We provide and comment our results in section IV, and present our conclusions in the final section V.
Ii Scaling Variables
Let us sketch the theoretical motivation for the different forms of scaling (14) in deepinelastic scattering and their extension to DVCS crosssection data, elastic vectormeson and for diffractive crosssections at fixed and small or We shall focus on the general arguments leading to scaling which could be independent from too specific theoretical predictions. The different types of scaling we refer to originate from properties of nonlinear QCD evolution equations, and the existence of asymptotic traveling wave solutions. We shall first introduce the general hint for these traveling wave solutions in the four cases and subsequently study the effect of varying the initial values of the phase space variables and . This amounts to introduce some natural non asymptotic effects.
For the three first cases (1, 2, 3), let us start with a “meanfield” equation for the dipoletarget amplitude of the type
(5) 
where , is the gluon transverse momentum (Fourier conjugate to the dipole size), is an arbitrary constant and is the rapidity. In Eq. (5), the coupling constant will be considered in the following as fixed, or running such that where and are the numbers of flavours and colours respectively. As we will show, for the sake of a general derivation of scaling, the differential operator kernel can be considered to be general, provided the function admits a minimum value at some point Indeed, the kernel at the leading logarithm (LL) accuracy of the perturbative QCD expansion, the BFKL kernel Lipatov:1976zz () verifies this property but large corrections at nexttoleading logarithm (NLL) accuracy exist (see, for the initial references, e.g. Fadin:1998py (); Salam:1998tj (); Ciafaloni:1999yw (); us ()). We will thus consider a general argument, instead of using a specific form of the QCD kernel.
Let us first consider the “Fixed Coupling” scheme The solution to the linear part of equation (5) corresponds to a linear superposition of waves:
(6) 
where is the Mellin transform of the kernel. In particular, one can interpret (6) by stipulating that each partial wave of wavenumber has a velocity By contrast, the critical velocity is defined GLR (); Munier:2003vc () by the minimum of the exponential phase factor is thus the velocity of the waves. So, the balance between the initial velocity (for fast enough initial condition with see Munier:2003vc ()) and the damping due to the nonlinear term in Eq (5) leads to an asymptotic of the velocity (and of the wave front) on the minimal and critical value As a consequence, reporting this dominant value in formula (6), one finds
(7) 
which is the expression of scaling with the variable (1).
For the running coupling case of Eq. 5, the situation is more subtle, since different approximate traveling wave solutions may coexist. In the original derivation of Ref. Munier:2003vc (), one starts with the linear solution analogue to Eq. (6), namely
(8) 
where plays the role of the new effective kernel. The key point is that the equivalent of the variable is now instead of This comes from a saddlepoint integration over the conjugate variable to the rapidity Ciafaloni:1999yw (). Consequently, after the same “blocking” mechanism due to the nonlinear damping terms, the scaling variable is which is nothing else than the “Running Coupling I” case, see Eq. (2). However, the scaling inferred by Eq. (8) is neither an exact nor unique approximate scaling solution.
As an alternative to that Mellin transform method, it was proposed recently in Ref. gb () to search directly for the running coupling analog of partial wave solutions, i.e. scaling solutions of Eq. (5) with running coupling with some generic velocity . Imposing the scaling of the left hand and right hand side of the equation as a function of the scaling variable leads to two different constraints, namely the functions and should be both functions of only
(9)  
(10) 
which are conditions to have a scaling solution In fact conditions (9) and (10) cannot be fulfiled simultaneously (contrary to the ‘Fixed Coupling” case, see gb ()) and can only be considered in an approximate way. If one choses to fulfil condition (10) exactly, and (9) only approximately, one finds again the scaling variable which is the “Running Coupling I” case, see Eq. (2). By this way one recovers the scaling variable previously found in Refs. Mueller:2002zm (); Munier:2003vc () in a different way. If, conversely, one chooses to satisfy exactly (9) and only approximately (10), one finds a new form of scaling, with the variable which corresponds to the “Running Coupling II” case, see Eq.(3). As the “blocking” mechanism is more general than the approximations leading to the variables “Running Coupling I” or “Running Coupling II”, one expects that Eq. (5) with running coupling has no exact scaling solution, but admits at least two different approximate solutions, namely “Running coupling I” and “Running coupling II”. One aim of our study is to make a phenomenological comparison between both solutions.
For the case of diffusive scaling, one considers the modification of Eq. (5) due to Pomeronloop contributions, which appear when the system of QCD partons or dipoles is dilute. It may seem irrelevant to an equation describing the high density partonic effects on the amplitude, but the structure of the equation, with an exponential increase of the linear regime may lead to important modifications ^{1}^{1}1Recent results on a toy model with running coupling and fluctuations Dumitru:2007ew () seem not to give diffusive scaling. A more complete QCD study is still missing. w.r.t. the meanfield solutions of (5). Indeed, it was shown that, while travelingwave solutions of the type (7) are formed, the effect of eventbyevent fluctuations leads to a stochastic superposition of these waves around an average solution which is governed by some diffusion coefficient and thus to a different scaling variable which corresponds to the “Diffusive Scaling” case, see Eq.(4).
The following comments are in order. One would like to confront the QCD predictions on the scaling form of the amplitude ( the wave front) with data. However, the precise determination of the wave front seems at present to be less reliable than the scaling variable ( the wave structure). The form of the linear kernel and of the non linear damping at nexttoleading logarithmic order have just been derived BKNLL (). However, in the kinematical range where scaling is observed, higher order contributions may be relevant, leading to a modified kernel in Eq. (5). Hence we will stick to a rather general and model independent (for each case) prediction for the scaling variable. The phenomenological analysis may on contrary help the development of the theoretical investigation.
Another important comment concerns the nonasymptotic corrections. While some of them can be deduced from the mechanism itself (see Munier:2003vc ()) we shall consider the simplest and expected nonasymptotic corrections due to the definition of the typical momentum and rapidity scales and Hence we will consider, if necessary the modification of all scaling variables (14), by a shift and Note that the reference value (without shift) is always .
Let us now consider diffractive, elastic vectormeson and DVCS crosssections. As discussed in Ref. Marquet:2006jb (), geometric scaling was shown to be valid and the “Fixed Coupling” was verified within uncertainties for these reactions, using instead of the Bjorken in the scaling variable (1).
In our phenomenological discussion, we shall thus extend the scaling properties to all four cases including the one studied in Marquet:2006jb (), by considering the various scaling variables (14).
We shall thus investigate the following scaling properties. The scaling prediction is
(11) 
for the DVCS crosssection,
(12) 
for the elastic vector meson cross sections, where the hard scale is assumed to be from the known properties of the vector meson wave function, and
(13) 
for the diffractive crosssection at fixed and small
For completion, besides the cases studied in Marquet:2006jb (), we may also consider the diffractive cross section at fixed and small namely
(14) 
since in the hypothesis of factorisation of the Pomeron flux, it relates to ( the total “Pomeronphoton” crosssection) at small for which the scaling arguments should also hold.
The prediction for DVCS (see Eq. (11)) does not rely on any nonperturbative assumption. This interesting feature of DVCS crosssections will allow us to include the data in the fitting procedure. Note that we will not use the diffractive data nor the vector meson ones for the fit. Instead, we will consider the fit to total and DVCS crosssections as scaling predictions for diffractive and vector meson data.
Iii The Quality factor method
In order to fit the parameters in the scaling variables (14), it is useful to use a quantity called quality factor usgeom () (QF). The quality factor is used to find the best parameters for a given scaling and to compare quantitatively the different scalings.
Given a set of data points and a parametric scaling variable we want to know whether the crosssection can be parametrised as a function of the variable only. Since we do not know the function of the variable describing the data, we need to define the QF independently of the form of that function. Therefore the QF is defined in such a way that it quantitatively describes how close the data points are to the scaling law , for a given .
Let us consider a set of points , where ’s are ordered, and introduce the QF as follows usgeom ()
(15) 
where is a small constant that prevents the sum from being infinite in case of two points having the same value of . Using this definition, the contribution to the sum in (15) is large when two successive points are close in and far in . Thus, we expect a set of points lying close to a unique curve to have large QF (small sum in (15)) compared to a QF of a set of points that are more scattered.
Since the crosssection in data differs by orders of magnitude and is more or less linear in (see 13), we decided to take and . This ensures that the low data points contribute to the QF with a similar weight as the higher data points. To complete the definition of the QF, the set is also rescaled so that , and ordered in before entering the QF formula. All the QFs in this paper are calculated with .
In order to test a scaling law , we search for the parameter
that minimises the 1/QF variable.
To do that, we fit 1/QF with the MINUIT
package. Given the maximum value of the QF,
we are able to directly compare different scaling laws. In this paper, we test the scaling laws
(14). The different scaling variables are given in
Table 1.
scaling  formula  parameters  

FC  “Fixed Coupling”  
RCI  “Running Coupling I”  
RCIb  “Running Coupling Ibis”  
RCII  “Running Coupling II”  
0.2 GeV  
RCIIb  “Running Coupling IIbis”  
DS  “Diffusive Scaling”  
1 GeV  
DSb  “Diffusive Scaling bis” 
Some comments are in order. We always define two versions of the scaling variables, the former depending only on a single parameter , and the later on two or three parameters . As an example, we define in this way the “Running Coupling II” and “Running Coupling IIbis” in this way respectively. The “Fixed Coupling” scaling does not depend on the scale and the rapidity since these additional variables disappear in the definition of the quality factor QF. The conclusion remains the same for the dependence of the “Running Coupling I” scaling. Concerning the “Running Coupling II” scaling, we fixed the parameter to a typical non perturbative scale inside the proton, GeV.
The QF defined above raise two additional technical issues which require further checks. First, it does not take into account the errors on data points. The errors can be introduced to the QF in an error weight function by modifying the definition (15) in the following way
(16) 
where is the error weight function of the data points . In the simplest form, we can take
(17) 
where is the relative statistical or uncorrelated error of point . This definition allows to take into account the statistical scatter of the different data points. This is particularly important if the statistical uncertainty differs from one data point to the other, and is significantly large. The data points we use in the following (mainly the structure function from the H1 and ZEUS experiments) show very little statistical uncertainties. We compared the results using either formula (15) or formula (16), and as expected, we get similar values of the parameters for both definitions. For sake of simplicity, we choose the QF (15) without any error weight function.
The second possible shortcoming is an ambiguity of the QF related to the point ordering if two or more ’s are equal. The order of these points in the data set is not unique; any permutation of these points satisfy the definition above. Obviously, we get diffrent QF for different permutations. The comparison of the fit results with different data point ordering was performed and the results were found similar (both in the value of the fitted parameters and the QF itself).
Iv Fits to and DVCS data
The first natural data set to test the different scaling laws are the data measured at low by the H1 H1 (), ZEUS ZEUS (), NMC NMC () and E665 E665 () experiments. These data are very precise and cover a wide range in and .
We choose to consider all available data from H1, ZEUS, NMC and E665 experiments with in the range GeV and . Fig. 1 displays the different kinematical domains reached by the H1, ZEUS, NMC and E665 experiments. While the H1 and ZEUS kinematic coverage is similar (with a tendency for H1 to reach lower and values), the data from NMC and E665 cover mostly the low , higher kinematic domain of our sample. We choose not to consider the data with since they are dominated by the valence and sea quark densities, and the formalism of saturation does not apply in this kinematical region. In the same way, the upper cut is introduced while the lower cut ensures that we stay away from the soft QCD domain. We will show in the following that the data points with GeV spoil the fit stability. Two kinds of fits to the scaling laws are performed, either in the full mentioned range, or in a tighter range GeV to ensure that we are in the domain where perturbative QCD applies. We will study the dependence of the fit parameters, as a possible indication of smooth scaling violations.
iv.1 Fits to data with GeV
We first give the results of the different scalings for the structure function with GeV. Fig. 26 show the normalised QF dependence on the value of the fitted parameters and the scaling curves in the GeV range for the “Fixed Coupling”, “Running Coupling I”, “Running Coupling II”, “Running Coupling II bis”, and “Diffusive Scaling” respectively. To see the impact of the different data sets, we choose to divide the data points in H1 alone, ZEUS alone, H1ZEUS, and all data. In the figures, we normalise the QF to 1. so that we are able to compare the values of the , and parameters for the different data sets. Table 2 give the fit parameters and absolute QF values — to distinguish the quality of the different scalings — for the various scaling laws (see Table I). On the different figures, we see that the parameter is well determined in the fit and the different data sets lead to values of which are quite close to one another. As an example, for fixed coupling, ranges between 0.33 to 0.38, for “Running Coupling I” between 1.74 and 1.84, for “Running Coupling II” between 3.23 and 3.44, and for “Diffusive Scaling” between 0.31 and 0.37. The scaling plots (always on the right of Fig.2, 3, 4, 5, and 6) show similar behaviour for the different scalings while it is clear that the scaling plot for diffusive scaling is slightly worse.
To improve the scaling quality, we checked the dependence of our results to additional parameters which can be introduced in the scaling. Introducing a shift in rapidity or a scale in the “Fixed Coupling” scaling does not change the results since these new parameters disappear in the definition of QF. This effect is the same if one introduces a scale in the “Running Coupling I” scaling. Introducing a new parameter to shift the rapidity in the “Running Coupling I” scaling does not lead to an improvement of the scaling. The fit leads to a value of compatible with 0 for all data sets considered. The results are similar when one introduces the parameter and the scale in “Diffusive Scaling” and the parameters are found to be close to 0 and 1 after the fit. We notice that “Diffusive Scaling” including the additional parameters and admits also other solutions leading to better QF, but these solutions are either instable from one data set to another or lead to unphysical values of the parameters (large values of and large negative values of ).
The QF plots in figure 5 show the dependence on one of the parameters while the remaining parameters are fixed to the fitted values. It is important to notice that in the case of the “Running Coupling IIbis” version, it is difficult to find the absolute minimum since the is not smooth and shows several minima. There is even a plateau of minimal values in , parameters and we picked the parameter values leading to the smallest values of QF. This is why the “Running Coupling IIbis” scaling law shows larger spread concerning the parameter values ( varies between 3.91 and 4.89). The typical values for the (respectively ) parameters vary between 1.2 and 2.5 (respectively 0.2 and 0.35). The values of are in a typical range of nonpertubative scale in the proton as expected.
We also notice in Table 2 that the best quality factor is obtained for “Running Coupling IIbis” for all data sets while the worst one is for the diffusive scaling. The comparison of the different QF is also shown in Fig. 9, left.
iv.2 Fits to data with GeV
Figure 7 shows the QF plots for different scalings in the range GeV. We notice that the spread on the parameter values is larger than after the cut on GeV. The values of the parameters are given in Table 2. We also notice a tendency to find lower values of for the “Fixed Coupling”, “Running Coupling I” and “Running Coupling II” scalings when one introduced more low data. The QF values are similar for the “Fixed Coupling”, “Running Coupling I”, and “Running Coupling IIbis” — with a tendency to be slightly better for “Running Coupling IIbis” — and is worse for diffusive scaling as before. The corresponding scaling curves are displayed in Fig. 8. They show similar behaviour with a tendency for “Diffusive Scaling” to show a larger spread of data. These figures also depict the GeV data points in grey, which were not used in the fit. These data points tend to show the same kind of scalings with a larger spread, and this is why they were not included in the fit.
Fig. 9, right, gives the curves for all data and different scaling laws in the regions GeV (on Fig. 9 left in the GeV range) so that the different scalings can be compared easily. As we mentionned already, in both ranges, the “Fixed Coupling”, the “Running Coupling I”, and the “Running Coupling II” are quite good, and the “Diffusive Scaling” is disfavoured. For GeV, the “Running Coupling IIbis” scaling shows the best QF.
iv.3 Dependence of the fitted parameters on
To study the dependence of the fitted parameters and to quantify the amplitude of scaling violations, the data points are divided into four separate samples: , , , and GeV. The number of points in each region is similar (see table 3). The lower and upper bounds of GeV and GeV are chosen for the same reason as before. The dependence of the parameter is depicted in figure 10, and the fitted parameters together with the QF values are shown in table 3. There is a slight increase of the parameter in the case of “Fixed Coupling” while “Running Coupling I” is quite flat. This can be easily understood since “Running Coupling I” shows a natural evolution. We notice a stronger increase in the case of “Running Coupling II”, showing the breaking of this scaling as a function of . Since this scaling gives already the best QF, it would be worth to study the breaking of scaling and introduce it in the model to improve further the description of the data. The parameter decreases strongly (especially in the last bin) for the diffusive scaling. The fact that depends strongly on confirms the fact that this scaling leads to the worst description of the data.
iv.4 Fits to DVCS data from H1 and ZEUS
After fitting all H1 and ZEUS data, it is worth studying whether the DVCS data measured by the same experiments DVCS () lead to the same results. The amount of data is smaller (34 points for H1 and ZEUS requiring as for data) and the precision on the parameter will be weaker. The kinematic coverage of the DVCS data covers a smaller region in and as : GeV and . The results of the fits can be found in table 4 and Fig.11 and 12. The scaling results are displayed only for “Fixed Coupling” and “Running Coupling IIbis” since all plots look similar. To facilitate the comparison between the results of the fits to and DVCS data, a star is put in Fig. 12 at the position of the value fitted to the H1+ZEUS data with in the range GeV. We note that the DVCS data lead to similar values to the data, showing the consistency of the scalings. The values of the QF show a tendency to favour fixed coupling, but all different scalings (even “Diffusive Scaling”) lead to reasonable values of QF.
V Implications for Diffraction and Vector Mesons
In this section, we check if the scalings found in the previous section can also describe diffractive and vector meson data. Since these data are much less precise than the or DVCS data and depend more on nonperturbative inputs (meson wave function, diffractive parton distribution inputs…), we choose to impose the same values of parameters found in the previous section and check if the scaling is also observed using this value. A fit to the diffractive data leads to large uncertainties on the parameter values since the number of data points is quite small and their uncertainties large.
Let us first describe the diffractive data. We use the latest integrated diffractive cross section measurements from H1 and ZEUS experiment F2d (), derived from the diffractive processes , where the proton stays either intact or turns into a low mass state . The diffractive cross section in terms of the integrated reduced cross section reads
(18)  
(19) 
The latter equation can be simplified to , since it is a very good approximation anywhere but at large . H1 large rapidity gap measurements are realised with GeV cut, whereas ZEUS measures with the GeV cut using the method. Both experiments do not measure exactly the same cross section. However, the difference is a known constant factor; the ZEUS data points can be converted to the same range as H1 by multiplying ZEUS values by factor of 0.85. In case of the processes where the proton is tagged in Roman pot detectors (both H1 and ZEUS), the data points have to be multiplied by 1.23.
In order to test the geometrical scaling properties we plot the dependence of , which can be expressed in terms of diffractive structure function as follows
(20) 
The definition of the variable is the same as in Table I replacing the variable by , remaining the same.
In figure 13, we show the H1 and ZEUS diffractive data for as a function of the fixed coupling scaling variable for six different fixed values of : 0.04, 0.1, 0.2, 0.4, 0.65, and 0.90. Only the data points in the region GeV and with are plotted. The restriction to the mentioned fixed values implies that the data set is dominated by H1 data. We just give in Fig. 13 the results for “Fixed Coupling” since the other scalings give similar results which are not distinguishable given the large uncertainties on . The parameters fitted to the H1 and ZEUS data are used to make the prediction on . Up to (the first five bins), the data points lie on a unique curve and thus confirm the geometrical scaling prediction. The data points at the highest values (low masses) are more sensitive to low mass resonances and are not expected to show a perfect scaling. The values of the QF are also given in Table 5. In many bins, “Running Coupling IIbis” gives the best QF.
Following the IngelmanSchlein model of the Pomeron, we also tested the scaling in , in the diffraction data. The definition of the variable is the same as in Table I replacing the variable by , remaining the same. Figure 14 presents the H1 and ZEUS diffractive data in the same region with at five different fixed values of : 0.0003, 0.001, 0.003, 0.01, and 0.03. The restriction to these fixed values results in favouring the H1 LRG data. Similarly as in the previous case, the parameter values from the fit to H1 and ZEUS data are used. The scaling is definitely not as good as the one for a fixed .
The fits to H1 and ZEUS data are also tested on vector meson data from both experiments vm (). Note the scale present in the scaling function in (12). The definition of the variable is the same as in Table I replacing the variable by being the mass of the vector meson, remaining the same. The QFs are given in table 6 and scaling curves for “Fixed Coupling” are given in Fig. 15. We only show the results for “Fixed Coupling” since all scalings give similar results. We see that the scaling is indeed verified for , , and and the data points show a tendency to be more dispersed for . Interestingly enough, “Diffusive Scaling” leads to the best QF.
Vi Conclusions
In this paper we analysed the scaling properties of deep inelastic observables such as the total cross section, DVCS, vector meson production, and the diffractive cross section, using all available low x data, i.e. from the H1 H1 (), ZEUS ZEUS (), NMC NMC () and E665 E665 () experiments. In order to be as model independant as possible, we used the quality factor method usgeom (), allowing to test the validity of scaling variables without any assumption concerning the shape of the scaling function. Four scaling variables (14), motivated by various approximations of QCD high energy evolution, were considered, including the new “Running Coupling II” one derived in gb (). Different versions of those scalings were tested, setting the rapidity shift to and the reference scale to either GeV or to GeV, or taking and as free parameters to be fitted.
We showed that the “Fixed Coupling”, “Running Coupling I” and “Running Coupling II” scaling variables lead to a good scaling behavior of the total cross section from the data, with similar quality factors. The “Diffusive scaling” is disfavoured compared to them. More precisely, the “Running Coupling IIbis”, with the free parameters and , gives the best scaling behavior in the range GeV GeV and , with reasonnable values of and extracted from the fit of the quality factor. The “Fixed Coupling” scaling is less deteriorated than the “Running Coupling II” scaling when one adds the points with GeV GeV to the fit, which corresponds to the boundary of the validity range of the perturbative scaling predictions. Dividing the data in four bins, we also studied the dependence of the optimal value for the parameter. “Running Coupling I” gives the most stable values and “Diffusive scaling” the less stable, and “Running Coupling II” gives stable values except at low .
The quality factors of the scaling variables were fitted independantly on the DVCS data, and give similar results as on the total cross section. In particular, for each scaling variable, the values of the parameter obtained from the fit in the two cases are very close.
We used the values of the parameters obtained from a fit to data to test the various scaling variables on the diffractive cross section and vector meson production data. For the former observable we tested both the fixed scaling behavior in and the fixed scaling behavior in . At fixed , we found a scaling behavior up to . For most bins the best results are given by “Running Coupling ” scaling variables, and the “Diffusive scaling” is often disfavoured.
At fixed , the scaling behavior of the diffractive cross section as a function of and is far less obvious. This is not a surprise, as not enough data is available in the genuine small region. A tendency of scaling is however observed for the bin.
Concerning , , and production, we found a
reasonable scaling behavior for all tested scaling variables, with the hard
scale , borrowed from vector mesons wave function studies.
Surprisingly, the best scaling is for all three vector mesons the “Diffusive
scaling”.
As a phenomenological outlook, it seems useful to work out models dipole amplitude wich could incorporate the successful scaling laws following the example of e.g. the Iancu Itakura Munier model IIM (). The quality factor method is a good tool to detect the scaling properties of data, and can guide the formulation of models with good . Our study is a good incentive for the formulation of appropriate models with running coupling scaling. There exists phenomenological models based on fixed coupling diffusive scaling Kozlov:2007wm (). The test of scaling and the formulation of model will be important to analyse future LHC data. In particular, the wider kinematical range open by the LHC could help disentangling the scaling laws in competition, or reveal a new scaling (such as the diffusive one) at higher energies.
On the theoretical ground, our phenomenological analysis can help to improve the theoretical analysis of scaling. Indeed, a mere comparison between the values of and the corresponding prediction based on the leading order BFKL kernel shows that one should go beyond this approximation to get an agreement. Thus, our determination of the basic scaling parameters can help the theoretical QCD analysis.
data  FC  RC I  RC II  RC II bis  DS  

all data  217  =0.330  =1.841  =3.436  =3.905  =0.362  
=1.200  
=1.63  =1.62  =1.69  =1.82  =1.44  
H1  87  =0.331  =1.736  =3.228  =4.891  =0.311  
=2.516  
=5.72  =5.88  =5.79  =6.31  =5.28  
ZEUS  127  =0.368  =1.809  =3.395  =4.327  =0.366  
=1.917  
=3.03  =3.08  =3.00  =3.28  =2.32  
H1ZEUS  214  =0.379  =1.839  =3.436  =4.147  =0.321  
=1.182  
=1.99  =1.76  =1.83  =2.02  =1.54  
all data  308  =0.321  =1.700  =2.932  =3.154  =0.369  
=0.199  
=1.30  =1.20  =1.07  =1.27  =1.02  
H1  135  =0.314  =1.710  =3.073  =3.159  =0.353  
=0.367  
=3.43  =3.56  =3.51  =3.52  =2.79  
ZEUS  147  =0.358  =1.809  =3.331  =3.747  =0.313  
=1.290  
=3.28  =3.20  =3.05  =3.29  =2.22  
H1ZEUS  282  =0.368  =1.797  =3.226  =3.918  =0.367  
=1.201  
=1.69  =1.71  =1.61  =1.74  =1.32 
FC  RC I  RC II  DS  
91  =0.279  =1.603  =1.627  =0.461  
=0.594  =0.575  =0.600  =0.571  
98  =0.301  =1.800  =3.219  =0.357  
=0.584  =0.544  =0.547  =0.526  
86  =0.367  =1.794  =3.521  =0.340  
=3.53  =3.06  =3.22  =2.67  
53  =0.397  =1.877  =4.135  =0.108  
=8.33  =8.00  =8.26  =6.37  

FC  RC I  RC II  RC II bis  DS  

34  =0.361  =1.829  =3.481  =5.717  =0.335 
=1.89  
=3.75  =3.62  =3.24  =3.52  =3.38 
FC  RC I  RC II  RC II bis  DS  

0.04  14  2.54  2.80  2.91  3.02  1.64 
0.1  30  0.610  0.579  0.600  0.605  0.660 
0.2  40  0.951  1.13  1.20  1.46  1.14 
0.4  64  1.05  0.952  0.984  0.998  1.00 
0.65  60  1.34  1.47  1.93  1.35  1.20 
0.9  59  0.380  0.510  0.572  0.492  0.372 
VM  FC  RC I  RC II  RC II bis  DS  

28  2.16  1.83  1.88  1.93  2.50  
62  1.02  0.814  0.803  0.964  1.28  
28  1.95  2.44  2.41  3.17  3.18 
References
 (1) A. M. Staśto, K. GolecBiernat, and J. Kwiecinski, Phys. Rev. Lett. 86, 596 (2001).
 (2) K. GolecBiernat and M. Wusthoff, Phys. Rev. D 59, 014017 (1999).
 (3) A.H. Mueller, “Parton saturation: An overview”, [arXiv:hepph/0111244].
 (4) I. Balitsky, Nucl. Phys. B 463, 99 (1996); Y. V. Kovchegov, Phys. Rev. D 60, 034008 (1999), D 61, 074018 (2000).
 (5) S. Munier and R. Peschanski, Phys. Rev. Lett. 91, 232001 (2003), Phys. Rev. D69, 034008 (2004), D 70, 077503 (2004).
 (6) G. Beuf, arXiv:0803.2167 [hepph].
 (7) E. Iancu, K. Itakura and L. McLerran, Nucl. Phys. A 708, 327 (2002); A. H. Mueller and D. N. Triantafyllopoulos, Nucl. Phys. B 640, 331 (2002); D. N. Triantafyllopoulos, Nucl. Phys. B 648, 293 (2003).
 (8) L. N. Lipatov, Sov. J. Nucl. Phys. 23, 338 (1976); E. A. Kuraev, L. N. Lipatov, and V. S. Fadin, Sov. Phys. JETP 45, 199 (1977), I. I. Balitsky and L. N. Lipatov, Sov. J. Nucl. Phys. 28, 822 (1978).
 (9) E. Iancu, A. H. Mueller and S. Munier, Phys. Lett. B 606, 342 (2005) [arXiv:hepph/0410018]; Y. Hatta, E. Iancu, C. Marquet, G. Soyez and D. N. Triantafyllopoulos, Nucl. Phys. A 773, 95 (2006) [arXiv:hepph/0601150].
 (10) C. Marquet and L. Schoeffel, Phys. Lett. B 639, 471 (2006) [arXiv:hepph/0606079].
 (11) G. Ingelman and P. E. Schlein, Phys. Lett. B 152, 256 (1985).
 (12) F. Gelis, R. Peschanski, G. Soyez and L. Schoeffel, Phys. Lett. B 647, 376 (2007) [arXiv:hepph/0610435].
 (13) V. S. Fadin and L. N. Lipatov, Phys. Lett. B 429, 127 (1998); M. Ciafaloni and G. Camici, Phys. Lett. B 430, 349 (1998).
 (14) G. P. Salam, JHEP 07, 019 (1998).
 (15) M. Ciafaloni, D. Colferai, and G. P. Salam, Phys. Rev. D 60, 114036 (1999).

(16)
R. Peschanski, C. Royon and L. Schoeffel,
Nucl. Phys. B 716, 401 (2005);
O. Kepka, C. Marquet, R. Peschanski and C. Royon, Phys. Lett. B 655, 236 (2007); Eur. Phys. J. C 55, 259 (2008).  (17) L. V. Gribov, E. M. Levin and M. G. Ryskin, Phys. Rept. 100 (1983) 1.
 (18) A. Dumitru, E. Iancu, L. Portugal, G. Soyez and D. N. Triantafyllopoulos, JHEP 0708, 062 (2007) [arXiv:0706.2540 [hepph]].
 (19) Y. V. Kovchegov and H. Weigert, Nucl. Phys. A 784 (2007) 188 [arXiv:hepph/0609090]; I. Balitsky, Phys. Rev. D 75, 014001 (2007) [arXiv:hepph/0609105]; I. Balitsky and G. A. Chirilli, Phys. Rev. D 77, 014019 (2008) [arXiv:0710.4330 [hepph]].
 (20) C. Adloff et al. [H1 Collaboration], Eur. Phys. J. C 21 (2001) 33 [arXiv:hepex/0012053]; Eur. Phys. J. C 30 (2003) 1 [arXiv:hepex/0304003].

(21)
J. Breitweg et al. [ZEUS Collaboration],
Eur. Phys. J. C 7 (1999) 609 [arXiv:hepex/9809005];
Phys. Lett. B 487 (2000) 273 [arXiv:hepex/0006013];
S. Chekanov et al. [ZEUS Collaboration],
Eur. Phys. J. C 21 (2001) 443 [arXiv:hepex/0105090].
Phys. Rev. D 70 (2004) 052001 [arXiv:hepex/0401003].
The data in the last ZEUS paper include contributions for and but those can be neglected within the kinematical domain we consider.  (22) M. Arneodo et al. [New Muon Collaboration], Nucl. Phys. B 483 (1997) 3 [arXiv:hepph/9610231].
 (23) M. R. Adams et al. [E665 Collaboration], Phys. Rev. D 54 (1996) 3006.
 (24) A. Aktas et al. [H1 Collaboration], Phys. Lett. B659 (2008) 796806 A. Aktas et al. [H1 Collaboration], Eur. Phys. J. C 44 (2005) 111 [arXiv:hepex/0505061]
 (25) A. Aktas et al. [H1 Collaboration], Eur. Phys. J. C 48 (2006) 715748 [arXiv:hepex/0606004] S. Chekanov et al. [ZEUS Collaboration], Nucl. Phys. B 713 (2005) 380 S. Chekanov et al. [ZEUS Collaboration], Eur. Phys. J. C 38 (2004) 4367
 (26) S. Chekanov et al. [ZEUS Collaboration], Nucl. Phys. B 718 (2005) 331 A. Aktas et al. [H1 Collaboration], Eur. Phys. J. C 46 (2006) 585603 [arXiv:hepex/0510016] A. Aktas et al. [H1 Collaboration], Eur.Phys.J. C 13 (2000) 371396 [arXiv:hepex/9902019]
 (27) E. Iancu, K. Itakura and S. Munier, Phys. Lett. B 590 (2004) 199 [arXiv:hepph/0310338].
 (28) M. Kozlov, A. Shoshi and W. Xiang, JHEP 0710 (2007) 020 [arXiv:0707.4142 [hepph]].