Measurement of the Collins and Sivers asymmetries on transversely polarised protons

The Collins and Sivers asymmetries for charged hadrons produced in deeply inelastic scattering on transversely polarised protons have been extracted from the data collected in 2007 with the CERN SPS muon beam tuned at 160 GeV/c. At large values of the Bjorken x variable non-zero Collins asymmetries are observed both for positive and negative hadrons while the Sivers asymmetry for positive hadrons is slightly positive over almost all the measured x range. These results nicely support the present theoretical interpretation of these asymmetries, in terms of leading-twist quark distribution and fragmentation functions.

After first indications of transverse spin effects in hadron physics in the 1970s [1,2] their importance was unambiguously established by the remarkably large single spin asymmetries (SSAs) found in pp collisions at Fermilab both for neutral and charged pions [3].Following the discovery by the EMC at CERN in 1988 that the quark spins contribute only little to the proton spin [4], the interest in the nucleon spin structure was revived and a more complete description including quark transverse spin and transverse momentum has been worked out.
The quark structure of the nucleon in the collinear approximation or after integration over the intrinsic quark transverse momentum k T is fully specified at the twist-two level by three parton distribution functions (PDFs) for each quark flavour [5]: the momentum distributions q(x), the helicity distributions ∆q(x) and the transverse spin distributions ∆ T q(x), where x is the Bjorken variable.The latter distribution-often referred to as transversity-is chiral-odd and thus not directly observable in deep inelastic scattering (DIS).In 1993 it was suggested [6] that transversity could be measured in semi-inclusive lepton-nucleon scattering (SIDIS) due to a mechanism involving another chiral-odd function in the hadronisation, known today as the Collins fragmentation function (FF).The mechanism leads to a left-right asymmetry in the distribution of the hadrons produced in the fragmentation of transversely polarised quarks.Thus a transverse spin dependence in the azimuthal distributions of the final state hadrons can be generated both in transversely polarised pp scattering and in SIDIS off transversely polarised nucleons.In the latter case the measurable Collins asymmetry, A Coll , is proportional to the convolution of the transversity PDF and the Collins FF.
Admitting a finite k T , in total eight PDFs are needed for a full description at leading twist and leading order in α S [7,8,9].All these functions lead to azimuthal asymmetries in the distribution of hadrons produced in SIDIS processes and can be disentangled measuring the different angular modulations.Amongst the transverse momentum dependent PDFs, the T -odd Sivers function [10] is of particular interest.This function arises from a correlation between the transverse momentum of an unpolarised quark in a transversely polarised nucleon and the nucleon spin.It can be different from zero because of final state interactions mediated by soft gluon exchange between the interacting quark and the target remnants [11].It is responsible for the Sivers asymmetry, A Siv , which is proportional to the convolution of the Sivers function and the unpolarised FF.The Sivers mechanism might also be the reason for the large asymmetries observed in pp collisions.
Transverse spin effects in SIDIS are investigated, at different beam energies, by the HERMES experiment at DESY and the COMPASS experiment at CERN.An experiment to measure transversity using a transversely polarised 3 He target has recently been performed at JLab [12].Transverse spin effects are also an important part of the scientific programme of the RHIC spin experiments at BNL.
Up to now, sizable Collins asymmetries for the proton were observed recently by HERMES using a proton target [13].This implies non-vanishing Collins fragmentation and transversity functions.Direct measurements at the KEK e + e − collider by the BELLE experiment established that this Collins FF is sizable [14,15].COMPASS measured vanishing asymmetries by scattering high energy muons off a deuteron target [16,17,18].All these data were well described by a global fit [19,20] which allowed for a first extraction of the u and d-quark transversity PDFs.
The Sivers asymmetry for the proton was measured by HERMES [13,21] to be different from zero for positive hadrons, while it was found to be compatible with zero for deuteron by COMPASS [16,17,18].These HERMES and COMPASS data could also be well described by theoretical calculations and fits, and allowed for extractions of the Sivers function [22], which turned out to be different from zero and opposite in sign for u and d-quarks.
In this Letter, we present the COMPASS results on the Collins and Sivers asymmetries for charged hadrons produced in SIDIS of high energy muons on transversely polarised protons.The data were collected in 2007 using NH 3 as target material and a 160 GeV/c beam with a momentum spread ∆p/p = ±5%.The beam was naturally polarised by the π-decay mechanism, with a longitudinal polarisation of about -80%.This measurement followed the measurements performed in 2002, 2003 and 2004 at the same energy with the transversely polarised 6 LiD target.
The COMPASS spectrometer [23] is in operation on the M2 beam line of CERN since 2002.Two magnetic stages are used to ensure large angular and momentum acceptance.A variety of tracking detectors is used to cope with the different requirements of position accuracy and rate capability at different angles.Particle identification is provided by a large acceptance RICH detector, calorimeters, and muon filters.Major upgrades in 2005 mainly concerned the polarised target, the tracking system, the RICH detector, and the electromagnetic calorimeters.The new target solenoid magnet provides a field of 2.5 T and has a polar angle acceptance of 180 mrad as seen from the upstream end of the target.In the earlier measurements with the 6 LiD target the polar angle acceptance was 70 mrad.The target material is cooled in a 3 He-4 He dilution refrigerator, and the protons in the H atoms are polarised to 0.80-0.90by dynamical nuclear polarisation.About 48 hours are necessary to reach 95% of the maximal polarisation.A pair of saddle-shaped coils can provide a 0.6 T vertical field which is used to rotate the target nucleon spin and to hold the polarisation vertical for the transversity measurements.In the frozen spin mode, and with the holding field at its operational value, the relaxation time of the polarisation exceeds 3000 hours.
The target consisted of three cylindrical cells with 4 cm diameter, one central cell of 60 cm length and two outer ones of 30 cm length, all separated by 5 cm.Neighbouring cells were polarised in opposite directions, so that data with both spin directions were recorded at the same time.In order to minimise the effects due to different spectrometer acceptance for different target cells, in each period of data taking a polarisation reversal was performed after 4-5 days by changing the microwave frequencies in the three cells.
The geometry of the polarised target and the data taking procedure were chosen such as to optimise the extraction of spin asymmetries.The principle of the measurement can be understood by considering the "ratio product" [17] where N ↑ inner and N ↓ outer are the number of hadrons produced in the first sub-period on oppositely polarised cells, and N ↓ inner and N ↑ outer are the corresponding numbers in the second sub-period, i.e. after polarisation reversal.The ratio product is constructed such that beam flux, spin-averaged cross-section, and the number of scattering centres cancel.As long as the ratios between the spectrometer acceptances of each cell are the same in the two sub-periods and the number of produced hadrons follows the generic azimuthal modulation N ↑↓ ∼ 1 ± ǫ sin Φ, one simply gets R = 1 + 4ǫ sin Φ, and the extraction of the amplitude ǫ of the azimuthal modulation is straightforward.
In 2007 data were taken at a mean beam intensity of about 5 × 10 7 µ + /s (typically 2.4 × 10 8 µ + /spill, for a spill length of 4.8 s every 16.8 s).Using up 4 × 10 13 muons, about 12 × 10 9 events were collected in six separate periods, corresponding to 440 TB of data.
In the data analysis, events were selected if they had at least one "primary vertex", defined as the intersection point of a beam track, the scattered muon track, and other possible outgoing tracks.The momenta of both incoming and outgoing charged particles were measured.The primary vertex was required to be inside a target cell.In order to guarantee the same muon flux along the target material, the extrapolated beam track had to traverse all the three target cells.For incoming and scattered muon tracks, as well as for the other reconstructed tracks, χ 2 cuts were applied to assure the quality of track reconstruction.Tracks from the primary vertex which traversed more than 30 radiation lengths were identified as scattered muons.The event was rejected if more than one of such tracks were found.
In order to be in the DIS regime, only events with a photon virtuality Q 2 > 1 (GeV/c) 2 , a fractional energy of the virtual photon 0.1 < y < 0.9, and a mass of the hadronic final state W > 5 GeV/c 2 were considered.The variable x covers the range from 0.004 to 0.7.
All particles emerging from the primary vertex were assumed to be hadrons if they traversed less than 10 radiation lengths of material.For tracks with an associated cluster in one of the hadronic calorimeters, a minimal amount of deposited energy was required to further reduce the electron and muon contamination.Finally, tracks reconstructed only in the fringe field of the first analysing magnet of the spectrometer were rejected.This roughly corresponds to a cut at 1.5 GeV/c in the hadron momenta.In order to reconstruct the hadron azimuthal angle with good precision, the hadron transverse momentum with respect to the virtual photon direction, p h T , was required to be above 0.1 GeV/c.A minimum value of 0.2 for z, the relative energy of the hadron with respect to the virtual photon energy, was chosen to avoid hadrons from the target fragmentation region.
As explained in detail in Ref. [17], the Collins effect shows up as a modulation [1 + ǫ C sin(φ h + φ S − π)] in the number of events, where φ h and φ S are the azimuthal angles of the hadron and of the target nucleon spin vector in a reference system in which the z-axis is the virtual photon direction and the x-z plane is the lepton plane according to Ref. [24]. .The amplitude of the modulation is ǫ C = D N N f P T A Coll , where D N N = (1 − y)/(1 − y + y 2 /2) is the transverse spin transfer coefficient from target quark to struck quark, f the dilution factor of the NH 3 material, and P T is the proton polarisation.Similarly, the Sivers effect results in a modulation [1+ǫ S sin(φ h −φ S )], where ǫ S = f P T A Siv .
The transverse spin asymmetries were obtained by comparing the azimuthal distributions of the detected hadrons as measured in the first sub-period of data taking with the corresponding distributions of the second half measured with opposite target polarisation.Since the two sets of data were taken typically one week apart, the stability of the apparatus is a central point in the measurement.As a first step in the data selection, the hit distributions of all trackers were scrutinised, as well as the number of reconstructed events, the number of vertices per events, and the number of tracks per event.In a second step, the stability of the average π + π − invariant mass in the K 0 region as well as the distribution of twelve kinematic quantities (x, y, W, z, ...) were investigated dividing the data in small time-ordered sub-samples.Each distribution of each sub-sample was compared with the corresponding ones of each other sub-sample within the same data taking period, and sub-samples were rejected when deviating more than 3.5 σ stat from the mean values.
As a final selection criterion, the data were tested for a possible dependence on either sin(φ h + φ S ) or sin(φ h − φ S ) of the acceptance ratio between two consecutive subperiods with opposite target polarisation.Combining the number of events reconstructed in the different target cells in two consecutive data taking sub-periods, one can construct two different estimators on the stability of the acceptance.The first estimator measures the mean modulation in the relevant azimuthal angle of the acceptance ratio between two sub-periods.The second one probes possible large differences in the acceptance ratios for the different target cells which could affect the physics asymmetry.These two pieces of information have been used to construct a χ 2 and the final selection of the data taking periods was done on the basis of its value.As a result of the quality control, all data collected in the six periods were used for the extraction of the Collins asymmetry, while only four periods were used for the Sivers asymmetry.This can be understood because the Sivers asymmetry is very sensitive to instabilities of the spectrometer since it is the amplitude of a modulation of the azimuthal angle of the hadron transverse momentum with respect to the target spin vector.On the contrary, the Collins asymmetry is an asymmetry in the azimuthal angle between the hadron transverse momentum and a direction which depends on the target spin direction and the lepton scattering plane, which is different for each event.The final sample contains 23.1×10 6 SIDIS events for the Collins asymmetry and 15.6×10 6 for the Sivers asymmetry.
The asymmetries were evaluated for positive and negative hadrons in bins of the three kinematic variables x, z and p h T .The binning is the same as used for the previous analyses of deuteron data and consists of 9 bins in x, 8 bins in z and 9 bins in p h T , integrating over the other two variables.For each period, the physics asymmetries were obtained by dividing the raw asymmetries by the target polarisation, the dilution factor, and, in the case of Collins analysis, by the D N N factor.The target polarisation was measured individually for each cell and each period.The dilution factor of the ammonia target was evaluated for each bin.It is 0.15 in average, and increases with x from 0.14 to 0.17.
The estimator used for the evaluation of the raw asymmetries is based on an extended unbinned maximum likelihood method [25].The likelihood function is built as the product of the probability densities p corresponding to each hadron i from each target cell.The likelihood for hadrons from a given target cell in one period is written as The + and − signs refer to the orientation of the target polarisation in the two subperiods and N ± is the corresponding total number of hadrons.The quantities I ± are the integrals of the probability densities over φ S and φ h .The probability densities p ± are the product of two parts, one corresponding to the acceptance description and the other to the SIDIS cross section of longitudinally polarised leptons on transversely polarised nucleons.Various parametrisations of the acceptance part were tested, resulting in a negligible dependence of the extracted asymmetries on the acceptance description.The cross section was parametrised taking into account both the unpolarised and polarised parts.The polarised part consists of all the expected eight modulations, namely sin(φ h + φ S − π), sin(φ h − φ S ), cos(φ h − φ S ), sin(2φ h − φ S ), cos(2φ h − φ S ), sin(φ S ), cos(φ S ), and sin(3φ h − φ S ), and all their amplitudes were extracted at the same time.The Collins and Sivers asymmetries are proportional to the amplitudes of the first two terms.Systematic studies for the other six amplitudes are still ongoing, and those results will be the subject The final Collins and Sivers asymmetries extracted with the likelihood method were compared with the asymmetries extracted using four other estimators, including those used in the previous publications which were based on the "ratio product" R of Eq. 1, finding an excellent agreement between all results.The correlation coefficient between the Collins and the Sivers asymmetries turned out to be small, less than 0.2 in absolute value over the whole x range.
Extensive studies were performed in order to assess the systematic uncertainty of the measured asymmetries.All the studies were done separately for positive and negative hadrons and for the Collins and the Sivers asymmetries.
The largest systematic error is due to residual acceptance variations within pairs of data taking sub-periods.To quantify these effects, two different types of false asymmetries were calculated, using the external cells and the internal cell divided in two parts, and assuming wrong sign polarisation for one of the two.Moreover, the physical asymmetries were also extracted using only the first and only the second half of the target.The difference between these two physical asymmetries, the false asymmetries, and the degree of compatibility of the results from different periods were all used to quantify the systematic uncertainty.
In the case of the Collins asymmetry, the systematic uncertainty is estimated to be 0.5 σ stat for positive and 0.6 σ stat for negative hadrons.In the case of the Sivers asymmetry, the systematic error is 0.8 σ stat for positive and 0.4 σ stat for negative hadrons.A further systematic uncertainty of ±0.01 is present in the absolute scale of the Sivers asymmetry for positive hadrons.It reflects a 0.02 difference in the mean value of the asymmetries extracted in the first two and in the second two periods of data taking used for this analysis.In spite of throughout studies, the origin of this difference, which affects only the Sivers asymmetry for positive hadrons, could not be identified and had therefore to be included in the systematic uncertainty.The results of this measurement of the Collins and Sivers asymmetries are shown in Fig. 1 and 2 as a function of x, z, and p h T , for positive and negative hadrons. Figure 3 displays the mean values of kinematic variables for positive hadrons in the x, z, and p h T bins.The corresponding quantities for negative  hadrons are very similar 1) .As it is clear from Fig. 1, the Collins asymmetry has a strong x dependence.It is compatible with zero at small x within the small statistical errors and increases in absolute value up to about 0.1 for x > 0.1.There, the values agree both in magnitude and in sign with the previous measurements of HERMES [13], which were performed at the considerably lower electron beam energy of 27.5 GeV.Also, the present results agree with the predictions of the global analysis of ref. [19,20] and thus strongly support the underlying interpretation of the Collins asymmetry in terms of a convolution of the twisttwo transversity PDF and the FF of a transversely polarised quark.An important issue is the Q 2 dependence of these functions.Our results at large x are compatible with the HERMES data in spite of the higher Q 2 values which exceed those of HERMES by a factor 2 to 3 with increasing x.This indicates that the possible Q 2 dependence should not 1) All numerical values have been put to HEPDATA.The results for the Sivers asymmetry for negative hadrons exhibit values compatible with zero within the statistical accuracy of the measurement.For positive hadrons, the data indicate small positive values, up to about 3% in the valence region.These values are somewhat smaller than but still compatible with the ones measured by HERMES at smaller Q 2 .Given the importance of the Sivers function in the present description of the transverse momentum structure of the nucleon, we looked at a possible kinematic dependence of our measurements.In particular, we evaluated the asymmetries as a function of W , separately for the "large-x" (x > 0.032) and "small-x" (x < 0.032) samples.The results are shown in Fig. 4. The mean values of Q 2 and x in all W bins are given in Fig. 5.As it is apparent from Fig. 4, no conclusion can be drawn about a possible W dependence of the Collins asymmetry.On the other hand, the signal of the Sivers asymmetry for positive hadrons seems to be concentrated at small W , in the region where HERMES measures, and goes to zero at large W , which for large x means large Q 2 .Thus our data give an indication for a possible W dependence of the Sivers asymmetry for positive hadrons.Definite conclusions will be possible only when new more precise data at high energy will become available.
In summary, for the first time the Collins and Sivers asymmetries for positive and negative hadron production in DIS off the proton have been measured at high energy.Our data extend the kinematic range to large Q 2 and large W values.The x range has been extended to considerably smaller values which are needed to evaluate the PDF first moments.For the Sivers asymmetry, a signal is seen for positive hadrons, which persists to rather small x values.The data give an indication for a possible W dependence of this asymmetry, but the present statistical and systematic uncertainties do not allow definite conclusions.The measured Collins asymmetry is sizable for both positive and negative hadrons also at high energies and Q 2 .Thus Collins asymmetries measured in SIDIS are an appropriate tool to investigate the transversity PDF.

Figure 1 :
Figure 1: Collins asymmetry as a function of x, z, and p h T , for positive (closed points) and negative (open points) hadrons.The bars show the statistical errors.The point to point systematic uncertainties have been estimated to be 0.5 σ stat for positive and 0.6 σ stat for negative hadrons and are given by the bands.

Figure 2 :
Figure 2: Sivers asymmetry as a function of x, z, and p h T , for positive (closed points) and negative (open points) hadrons.The bars show the statistical errors.The point to point systematic uncertainties have been estimated to be 0.8 σ stat for positive and 0.4 σ stat for negative hadrons and are given by the bands.For positive hadrons only, an absolute scale uncertainty of ±0.01 has also to be taken into account.

Figure 3 :
Figure 3: Mean values of some kinematic variables in the final data sample.From left to right: mean values of p h T , z and Q 2 as functions of x; mean values of p h T , x and Q 2 as functions of z; mean values of x, z and Q 2 as functions of p h T .

Figure 4 :
Figure 4: Collins (upper row) and Sivers (lower row) asymmetry as a function of W , for positive (left) and negative (right) hadrons.The closed and open points give the values for the "large x" and the "small x" samples respectively.The errors are statistical only.

Figure 5 :
Figure 5: Mean values of Q 2 (left) and x (right) as functions of W .The closed and open points give the values for the "large x" and the "small x" samples respectively.