An approach to NLO QCD analysis of the semi-inclusive DIS data with modified Jacobi polynomial expansion method

It is proposed the modification of the Jacobi polynomial expansion method (MJEM) which is based on the application of the truncated moments instead of the full ones. This allows to reconstruct with a high precision the local quark helicity distributions even for the narrow accessible for measurement Bjorken x region using as an input only four first moments extracted from the data in NLO QCD. It is also proposed the variational (extrapolation) procedure allowing to reconstruct the distributions outside the accessible Bjorken x region using the distributions obtained with MJEM in the accessible region. The numerical calculations encourage one that the proposed variational (extrapolation) procedure could be applied to estimate the full first (especially important) quark moments. The extraction of the quark helicity distributions is one of the main tasks of the semi-inclusive deep inelastic scattering (SIDIS) experiments (HERMES [1], COMPASS [2]) with the polarized beam and target. At the same time it was argued [3] that to obtain the reliable distributions at relatively low average Q available to the modern SIDIS experiments, the leading order (LO) analysis is not sufficient and next to leading order analysis (NLO) is necessary. In ref. [4] the procedure allowing the direct extraction from the SIDIS data of the first moments of the quark helicity distributions in NLO QCD was proposed. However, in spite of the special importance of the first moments, it is certainly very desirable to have the procedure of reconstruction in NLO QCD of the polarized densities themselves. However, it is extremely difficult to extract the local in xB distributions directly, because of the double convolution product entering the NLO QCD expressions for semi-inclusive asymmetries (see [4] and references therein). On the other hand, operating just as in ref. [4], one can directly extract not only the first moments, but the Mellin moments of any required order. The simple extension of the procedure proposed in ref. [4] gives for the n-th moments ∆nq ≡ ∫ 1 0 dx x q(x) of the valence distributions the equations ∆nuV = 1 5 A p +A (n) d L(n)1 − L(n)2 ; ∆ndV = 1 5 4A (n) d −A (n) p L(n)1 − L(n)2 , (1) where all quantities in the right-hand side are the same as in ref. [4] (see Eqs. (18-23)) with the replacement of ∫ 1 0 dx by ∫ 1 0 dx x . E-mail address: [email protected] E-mail address: [email protected] E-mail address: [email protected] For example, HERMES data [1] on semi-inclusive asymmetries is obtained at Qaverage = 2.5GeV .