Transversity distributions from difference asymmetries in semi-inclusive DIS

TL;DR

This paper presents a novel method to extract the ratio of u and d quark transversity distributions from SIDIS data alone, avoiding the need for $e^+e^-$ data, and applies it to COMPASS data.

Contribution

The paper introduces an alternative approach to determine transversity ratios using difference asymmetries in SIDIS without relying on $e^+e^-$ data.

Findings

Successfully extracted the $h_1^d/h_1^u$ ratio from COMPASS data.

Compared results with previous methods using SIDIS and $e^+e^-$ data.

Validated the new method's consistency with existing measurements.

Abstract

In recent years information on the transversity distribution $h_1$ has been obtained combining the Collins asymmetry results from semi-inclusive deep inelastic scattering (SIDIS) data on transversely polarized nucleon targets with the information on the fragmentation function of a transversely polarized quark from the asymmetries measured in $e^+e^-$ annihilation into hadrons. An alternative method was proposed long time ago, which does not require the $e^+e^-$ data, but allows one to get ratios of the $u$ and $d$ quark transversity distributions from the SIDIS data alone. The method utilizes the ratio of the difference of the Collins asymmetries of positively and negatively charged hadrons produced on transversely polarized proton and deuteron targets. We have applied this method to the COMPASS proton and deuteron data, and extracted the ratio $h_1^d/h_1^u$. The results are compared to those obtained in a previous point--by--point extraction based both on SIDIS and $e^+e^-$ data.

Equations

$${\sigma }_t^{\pm }({\Phi }_C)={\sigma }_{0,t}^{\pm }+f{P}_T{D}_{NN}{\sigma }_{C,t}^{\pm }\sin {\Phi }_C$$

$$A_{C,t}^{\pm }=\frac{{\sigma }_{C,t}^{\pm }}{{\sigma }_{0,t}^{\pm }}$$

$${\sigma }_{0,p}^{+}$$

$${\sigma }_{0,p}^{-}$$

$${\sigma }_{0,d}^{+}$$

$${\sigma }_{0,d}^{-}$$

$$H_{1}^{\perp(1/2)}(z,Q^{2})\equiv\int\mathrm{d}^{2}\mathchoice{\mbox{\boldmath$\displaystyle p$\unboldmath}}{\mbox{\boldmath$\textstyle p$\unboldmath}}{\mbox{\boldmath$\scriptstyle p$\unboldmath}}{\mbox{\boldmath$\scriptscriptstyle p$\unboldmath}}_{T}\,\frac{p_{T}}{zM_{h}}\,\,H_{1}^{\perp}(z,p_{T}^{2},Q^{2})\,.$$

$${\sigma }_{C,p}^{+}$$

$${\sigma }_{C,p}^{-}$$

$${\sigma }_{C,d}^{+}$$

$${\sigma }_{C,d}^{-}$$

$$A_{D,t}=\frac{{\sigma }_{C,t}^{+}-{\sigma }_{C,t}^{-}}{{\sigma }_{0,t}^{+}+{\sigma }_{0,t}^{-}}.$$

$$A_{D,t}^{\prime }=\frac{{\sigma }_{C,t}^{+}-{\sigma }_{C,t}^{-}}{{\sigma }_{0,t}^{+}-{\sigma }_{0,t}^{-}}.$$

$$A_{D,p}$$

$$A_{D,d}$$

$$\frac{A_{D,d}}{A_{D,p}}=3\left[\frac{(4f_1^u+4f_1^{\bar{u}}+f_1^d+f_1^{\bar{d}})(D_{1,\mathrm{fav}}+D_{1,\mathrm{unf}})+2(f_1^s+f_1^{\bar{s}})D_{1,s}}{5(f_1^u+f_1^d+f_1^{\bar{u}}+f_1^{\bar{d}})(D_{1,\mathrm{fav}}+D_{1,\mathrm{unf}})+4(f_1^s+f_1^{\bar{s}})D_{1,s}}\right]\frac{h_1^{u_v}+h_1^{d_v}}{4h_1^{u_v}-h_1^{d_v}},$$

$${\sigma }_t^D({\Phi }_C)=({\sigma }_{0,t}^{+}-{\sigma }_{0,t}^{-})+f{P}_T{D}_{NN}({\sigma }_{C,t}^{+}-{\sigma }_{C,t}^{-})\sin {\Phi }_C$$

$$A_{D,t}$$

$$\frac{A_{D,d}^{\prime }}{A_{D,p}^{\prime }}=\frac{4f_1^{u_v}-f_1^{d_v}}{f_1^{u_v}+f_1^{d_v}}\frac{h_1^{u_v}+h_1^{d_v}}{4h_1^{u_v}-h_1^{d_v}}$$

$$A_{D,t}^{\prime }=\frac{{\sigma }_{0,t}^{+}}{{\sigma }_{0,t}^{+}-{\sigma }_{0,t}^{-}}{A}_{C,t}^{+}-\frac{{\sigma }_{0,t}^{-}}{{\sigma }_{0,t}^{+}-{\sigma }_{0,t}^{-}}{A}_{C,t}^{-}.$$

$$A_{D,t}=\frac{\mathrm{var}(A_{C,t}^{-})}{\mathrm{var}(A_{C,t}^{+})+\mathrm{var}(A_{C,t}^{-})}{A}_{C,t}^{+}-\frac{\mathrm{var}(A_{C,t}^{+})}{\mathrm{var}(A_{C,t}^{+})+\mathrm{var}(A_{C,t}^{-})}{A}_{C,t}^{-}$$

Full text

I Introduction

Much interest has been dedicated in the past twenty years to the transversity distribution. Usually called $h_{1}$, it is a leading-twist parton distribution function (PDF) which describes the transverse polarization of quarks inside a transversely polarized nucleon (for reviews, see Barone:2010zz ; Aidala:2012mv ; Avakian:2016rst ).

Being chirally-odd, transversity cannot be measured in deep inelastic scattering (DIS). Over the last decade single-spin asymmetries clearly related to the transversity distribution function have been measured in semi-inclusive deep inelastic scattering (SIDIS) on transversely polarized nucleons, namely in DIS processes in which at least one hadron of the current jet is detected. In these processes the cross-section exhibits a spin-dependent azimuthal modulation which can be expressed in terms of a convolution of the transversity PDF and a fragmentation function (FF) which is also chiral-odd, thus guaranteeing the cross-section to be chirally-even. Two observables have been studied in so far. The first one is single hadron spin asymmetry, namely the amplitude of the target-spin dependent azimuthal modulation of each of the produced hadrons. The second one is the amplitude of the target spin-dependent azimuthal modulation of the plane defined by any two of the oppositely charged hadrons produced in the same SIDIS event. In the case of transversely polarized proton targets clear non zero azimuthal modulations have been measured for both observables by the HERMES Airapetian:2010ds ; Airapetian:2008sk and the COMPASS Adolph:2012sn ; Adolph:2014fjw Collaborations, assessing beyond any doubt that both the transversity PDFs and the single hadron and the dihadron FF are not zero. Corresponding measurements on a transversely polarized deuteron target by the COMPASS Collaboration Ageev:2006da ; Adolph:2012nw gave asymmetries compatible with zero, which have been interpreted as evidence of cancellation between $h_{1}^{u}$ and $h_{1}^{d}$.

The underlying physics of these processes Adolph:2015zwe ; Kerbizi:2018qpp is the left-right asymmetry in the hadronization of a transversely polarized quark, where left and right are relative to the plane defined by the quark direction of motion and its transverse spin. Such asymmetry is encoded, in the first case, in the so-called Collins FF $H_{1}^{\perp}$ Collins:1992kk , and, in the second case in the dihadron FF $H^{\sphericalangle}_{1}$  Collins:1994kq ; Radici:2001na . Independent evidence that both the Collins function and the dihadron FF are different from zero came from the measurements of azimuthal asymmetries in hadron inclusive production in $e^{+}e^{-}$ annihilation by the Belle Seidl:2008xc ; Vossen:2011fk , the BaBar TheBABAR:2013yha and the BESIII Ablikim:2015pta Collaborations.

Combining the SIDIS data and the $e^{+}e^{-}\rightarrow\mathrm{hadrons}$ measurements, first extractions of both the transversity functions and of the two transversely polarized quark FFs have been possible Anselmino:2015sxa ; Kang:2015msa . In all those works, in order to solve the convolution over the transverse momenta between the transversity PDF and the FF which appears in the cross-section, some parametrization for both $h_{1}$ and for the FFs had to be assumed. An exception is the recent extraction of transversity Martin:2014wua which has been performed point-by-point directly from the COMPASS SIDIS and the Belle $e^{+}e^{-}\rightarrow\mathrm{hadrons}$ data, without using any parametrization for the collinear variables.

An alternative way to measure transversity from the Collins asymmetries alone is via the so-called “difference asymmetries”, which allow extracting combinations of the $u$ and $d$ quark transversity without knowing the Collins FF. This method was proposed a long time ago Frankfurt:1989wq ; Christova:2000nz ; Sissakian:2006vz to access the helicity PDFs, and has been used by the SMC Collaboration Adeva:1995yi . It was also included in the COMPASS proposal Baum:1996yv , to measure both longitudinal and transverse spin asymmetries. At that time it looked particularly interesting, since the Collins FF was completely unknown. Later on it has been used to measure the helicity PDFs in COMPASS Alekseev:2007vi , and recently it has been proposed again in the context of the Sivers, Boer-Mulders and transversity distributions Christova:2015waz . In the present work the difference asymmetries are used for the first time to access transversity using the COMPASS measurements of the Collins asymmetries on $p$ Adolph:2012sn and $d$ targets Ageev:2006da .

The paper is organized as follows. In Section II we describe the formalism and the procedure. Section III is dedicated to the Monte Carlo studies. The results are discussed in Section IV.

II Cross sections and difference asymmetries

In this paper we will extract the asymmetries of differences from the Collins asymmetries measured by the COMPASS Collaboration impinging a 160 GeV/c momentum muon beam either on a transversely polarized deuteron (6LiD) target or a transversely polarized proton (NH3) target. The data we have considered were taken in the years 2002–2004 Ageev:2006da and 2010 Adolph:2012sn .

In order to ensure the DIS regime, only events with photon virtuality $Q^{2}>1$ (GeV/c)2, fractional energy of the virtual photon $0.1<y<0.9$, and mass of the hadronic final state system $W>5$ GeV/c2 were considered in the data analysis. The charged hadrons were required to have at least 0.1 GeV/c transverse momentum $P_{T}$ with respect to the virtual photon direction and a fraction of the available energy $z>0.2$. All the details of the event selection and of the analysis can be found in Ageev:2006da ; Adolph:2012sn . The published data are binned in $x$, the target nucleon momentum fraction carried by the struck quark, in $z$ and in $P_{T}$. In our analysis we will only consider the asymmetry data binned in $x$, in 9 bins, from 0.003 to 0.7.

In the following, for simplicity we will write explicitly only the Collins part of the SIDIS transverse spin dependent cross-section, and consider charged pions, even if, at the end, we will use the results for charged hadrons assuming they are all pions, as it was done, for instance, in Martin:2014wua . The SIDIS cross section can be written as

[TABLE]

where $\Phi_{C}$ is the Collins angle, $f$ is the target dilution factor, $P_{T}$ is the nucleon polarization, and $D_{NN}$ is the mean transverse-spin-transfer coefficient not included in $\sigma_{C}$ to simplify the expressions used in the following. Only the deuteron (or hydrogen) nuclei in the targets were polarized, and the target dilution factor $f$ is given by the ratio of the absorption cross-sections on the deuteron (or proton) to that of all nuclei in the target. The signs $\pm$ refer to the pion charge and $t=p,\ d$ is the target type. The Collins angle $\Phi_{C}=\phi_{h}+\phi_{S}-\pi$ is the sum of the azimuthal angles $\phi_{h}$ of the hadron transverse momentum and of the spin direction $\phi_{S}$ of the target nucleon with respect to the lepton scattering plane, in a reference system in which the $z$ axis is the virtual photon direction.

The Collins asymmetry is defined as

[TABLE]

In terms of the ordinary PDFs and FFs the unpolarized part of the cross-sections in eq. (1) can be written as (omitting a kinematic factor that cancels out when taking the ratios of cross sections):

[TABLE]

where $D_{1,{\rm fav}}$ ($D_{1,{\rm unf}}$) is the favored (unfavored) unpolarized FF, $D_{1,s}$ is the strange sea unpolarized FF, and $f_{1}^{q}$ are the unpolarized PDFs.

Following Martin:2014wua , the corresponding spin–dependent cross sections are obtained by replacing $f_{1}^{q}$ with the transversity PDFs $h_{1}^{q}$ and the FFs $D_{1}$ with the “half moments” of the Collins function, $H_{1}^{\perp(1/2)}$, defined as

[TABLE]

Thus we have:

[TABLE]

where we have assumed $H_{1,s}^{\perp(1/2)}=0$.

We now define the difference asymmetries as

[TABLE]

In Christova:2000nz an alternative definition was proposed, namely

[TABLE]

As we will see, the two definitions turn out to give the same results. For the sake of simplicity, our discussion in the following we be centered on the definition (12), but we will also briefly summarize the results obtained with eq. (13).

Writing explicitly the asymmetries one gets:

[TABLE]

When taking the ratios of the asymmetries on deuteron and proton, the Collins FFs cancel out:

[TABLE]

and the only unknowns are the transversity PDFs. Thus, by measuring $A_{D}$ on $p$ and $d$, one obtains the ratio $h_{1}^{d_{v}}/h_{1}^{u_{v}}$ in terms of known quantities.

In order to determine $A_{D,t}$, one should in principle fit the quantity

[TABLE]

and extract the amplitude of the $\sin\Phi_{C}$ modulation. Since usually the acceptances for positively and negatively charged particles are not the same, one should correct the number of events for the acceptance before taking the differences, and treat carefully the statistical errors.

The measurements are much simpler if the $\Phi_{C}$ acceptance for positively charged particles is equal to that for negatively charged ones. In this case it is not necessary to evaluate the difference asymmetries from the amplitude of the modulation, as it is possible to get them from the measured Collins asymmetries. One has in fact

[TABLE]

where the ratios of the cross sections are known. In order to apply this procedure, extensive Monte Carlo studies have been performed. They are described in the next Section.

Notice that if one uses instead the definition (13), the ratio of the difference asymmetries has the form

[TABLE]

and the equivalent of eq. (18) is

[TABLE]

III Monte Carlo studies

The acceptance of the COMPASS spectrometer for positively charged and negatively charged hadrons have been investigated with Monte Carlo simulations. In the case of the deuteron data, collected in the years 2002-2004, this work was a prerequisite to the extraction of the $\sin\phi_{h}$, $\cos\phi_{h}$ and $\cos 2\phi_{h}$ modulations  Adolph:2014pwc which are expected in the unpolarized SIDIS cross-section. Within the statistical errors, the acceptance turned out to be essentially the same for positive and negative hadrons. In 2010, when the proton data were collected, the spectrometer was substantially different from the one utilized for the deuteron data taking were taken, thus the whole work had to be repeated. To this end we have used a full Monte Carlo chain using LEPTO Ingelman:1996mq as event generator and TGEANT, a GEANT4 Asai:2015xno based program, for the simulation of the particle interaction with the COMPASS apparatus and the detector response. The Monte Carlo events have been reconstructed with the COMPASS package CORAL Abbon:2007pq , and analyzed to extract the acceptances and the acceptance ratios. The same kinematic selections used for the analysis of the real data have been applied on the generated variables and on the reconstructed ones. While integrating over the other kinematical variables, the acceptances ${\mathcal{A}}(x)={N^{\rm rec}(x)}/{N^{\rm gen}(x^{\rm gen})}$ have been obtained by taking the ratio of the reconstructed and generated events counted in every bin using respectively the generated $x^{\rm gen}$ and the reconstructed $x^{\rm rec}$ values. In this way also the smearing due to the experimental resolution is accounted for.

The acceptances, which include both the geometrical acceptance of the apparatus and the reconstruction efficiency, are shown in Fig. 1 (left). The acceptances for positively (red points) and negatively (black points) charged hadrons are in good agreement and the small differences are compatible with the statistical fluctuations. Their ratios is constant over the full $x$ range of the measurement, with an average value of $1.003\pm 0.006$.

A possible Collins modulation in the acceptance was also studied, separately for positively and negatively charged hadrons, by fitting in each $x$ bin the $\Phi_{C}$ distribution with a function $c\,(1+a_{C}\sin(\Phi_{C}))$. The results for $a_{C}$ are shown in Fig. 1 (right). The amplitudes of the modulation are compatible with zero over the full $x$ range for both positive and negative hadrons. This result stays true also when repeating the procedure for the ratio of the acceptances.

IV Results

On the basis of the Monte Carlo results, the difference asymmetries have been calculated using eq. (18) with the Collins asymmetries from the 2010 COMPASS data. Actually, since $\sigma^{\pm}_{0,t}\sim N^{\pm}_{t}$ and $\mathrm{var}(A^{\pm}_{C,t})\sim 1/N^{\pm}_{t}$, where $N^{\pm}_{t}$ is the total number of hadrons which has been used to extract the Collins asymmetries, in a given $x$ bin, eq. (18) can be rewritten as:

[TABLE]

The calculation of the difference asymmetries can thus be performed using the published COMPASS data for $A^{\pm}_{C,t}$ and their statistical uncertainties Adolph:2012sn . An interesting remark is that $A_{D,t}$ is equal to the weighted mean of the Collins asymmetries for positive and negative hadrons, after changing sign to $A^{-}_{C,t}$. The results for proton and deuteron are shown in Fig. 2.

The ratio $A_{D,d}/A_{D,p}$ is shown in Fig. 3. Only the four points at larger $x$ are plotted in the figure. The points at smaller $x$ have much too large uncertainties since the proton asymmetries in that region are compatible with zero.

From the ratios $A_{D,d}/A_{D,p}$ the quantities $(h_{1}^{u_{v}}+h_{1}^{d_{v}})/(4h_{1}^{u_{v}}-h_{1}^{d_{v}})$ have been extracted using eq. (16) and standard parametrizations and tables for the unpolarized PDFs Lai:1999wy and FFs deFlorian:2007aj .

Finally, from the quantities $(h_{1}^{u_{v}}+h_{1}^{d_{v}})/(4h_{1}^{u_{v}}-h_{1}^{d_{v}})$ the ratios $h_{1}^{d_{v}}/h_{1}^{u_{v}}$ are determined. They are shown as closed circles in Fig. 4. Again, the values in the first five $x$ bins have very large uncertainties, are compatible with zero and are not plotted in the figure. At larger $x$ the values are negative, in agreement with previous extractions. The same procedure has been carried through starting from the difference asymmetries $A^{\prime}_{D,t}$ and using eq. (19), getting essentially the same values and the same statistical uncertainties, which are shown as closed squares in Fig. 4. In the same figure we also compare our results with the values of $h_{1}^{d_{v}}/h_{1}^{u_{v}}$ calculated from the transversity values obtained in Martin:2014wua (open circles). In the evaluation of the uncertainty of the ratio $h_{1}^{d_{v}}/h_{1}^{u_{v}}$ from Martin:2014wua , proper account has been taken of the correlations between the extracted values of $h_{1}^{u_{v}}$ and $h_{1}^{d_{v}}$, and use has been made of the correlation coefficients as evaluated in addendum3 . The results of the three determinations are in very good agreement, but some reduction (up to $\sim 20\%$) of the uncertainties can be observed in the ratios obtained in the present work from the difference asymmetries.

V Conclusion

We have determined for the first time the transverse-spin difference asymmetries of positively and negatively charged hadrons using the SIDIS $p$ and $d$ COMPASS data. Thanks to the good COMPASS spectrometer acceptance it could be easily obtained from the measured Collins asymmetries. From the ratio between the difference asymmetries on deuterons and on protons we have extracted the quantity $h_{1}^{d_{v}}/h_{1}^{u_{v}}$, the ratio between the valence $d$-quark and $u$-quark transversity PDF.

At small $x$ the difference asymmetries on the protons are compatible with zero, thus the statistical uncertainty on the ratio $h_{1}^{d_{v}}/h_{1}^{u_{v}}$ is too large and no useful information is provided by the present analysis. On the other hand, for larger $x$ ($x\leq 0.05$) the extracted ratio $h_{1}^{d_{v}}/h_{1}^{u_{v}}$ has negative sign and is in very good agreement with the results of a previous point by point extraction.

The method we applied is interesting and simple, and does not require any knowledge of the Collins fragmentation functions. Hence it strengthens the validity of the methods utilized so far to extract the transversity distributions, based on a combined analysis of SIDIS and $e^{+}e^{-}$ data, and can be used as a useful cross-check for more elaborated extractions.

Acknowledgements.

This work has been possible thanks to the project FRA2015, supported by the Università degli Studi di Trieste. V.B. has been partially supported by “Fondi di Ricerca Locale ex–60%” of the University of Piemonte Orientale. We would like to thank E. Christova and E. Leader for interesting discussions. We are grateful to the COMPASS Collaboration for the use of the Monte Carlo chain.