Ratios of helicity amplitudes for exclusive $\rho^0$ electroproduction on transversely polarized protons

TL;DR

This study measures helicity amplitude ratios in exclusive $ ho^0$ electroproduction on transversely polarized protons, providing new insights into nucleon-helicity-flip processes and supporting models based on generalized parton distributions.

Contribution

It is the first to extract ratios of nucleon-helicity-flip amplitudes using transversely polarized targets in $ ho^0$ electroproduction, aligning experimental results with handbag approach predictions.

Findings

Ratios of nucleon-helicity-non-flip amplitudes agree with previous analyses.

First extraction of nucleon-helicity-flip amplitude ratios with transverse polarization.

Data support a positive $ ho- ext{pion}$ transition form factor.

Abstract

Exclusive $\rho^0$-meson electroproduction is studied by the HERMES experiment, using the 27.6 GeV longitudinally polarized electron/positron beam of HERA and a transversely polarized hydrogen target, in the kinematic region 1.0 GeV$^2$<Q$^2$<7.0 GeV$^2$, 3.0 GeV<W<6.3 GeV, and -t'<0.4 GeV$^2$. Using an unbinned maximum-likelihood method, 25 parameters are extracted. These determine the real and imaginary parts of the ratios of several helicity amplitudes describing $\rho^0$-meson production by a virtual photon. The denominator of those ratios is the dominant amplitude, the nucleon-helicity-non-flip amplitude $F_{0\frac{1}{2}0\frac{1}{2}}$, which describes the production of a longitudinal $\rho^{0}$-meson by a longitudinal virtual photon. The ratios of nucleon-helicity-non-flip amplitudes are found to be in good agreement with those from the previous HERMES analysis. The transverse target polarization allows for the first time the extraction of ratios of a number of nucleon-helicity-flip amplitudes to $F_{0\frac{1}{2}0\frac{1}{2}}$. Results obtained in a handbag approach based on generalized parton distributions taking into account the contribution from pion exchange are found to be in good agreement with these ratios. Within the model, the data favor a positive sign for the $\pi-\rho$ transition form factor. By also exploiting the longitudinal beam polarization, a total of 71 $\rho^0$ spin-density matrix elements is determined from the extracted 25 parameters, in contrast to only 53 elements as directly determined in earlier analyses.

Tables

Table 1. Table 1: Mean values for the kinematic variables W 𝑊 W , Q 2 superscript 𝑄 2 Q^{2} , − t ′ superscript 𝑡 ′ -t^{\prime} , and f b ​ g subscript 𝑓 𝑏 𝑔 f_{bg} under the exclusive peak in each of the ( Q 2 superscript 𝑄 2 Q^{2} , − t ′ superscript 𝑡 ′ -t^{\prime} ) cells.

Cell limits	$⟨W⟩$, GeV	$⟨Q^2⟩$, GeV2	$⟨-t^{\prime }⟩$, GeV2	$f_{bg}$
$1.0$ GeV GeV2	4.70	1.19	0.128	0.065
$1.4$ GeV GeV2	4.75	1.67	0.128	0.073
$2.0$ GeV GeV2	4.80	3.06	0.136	0.122
$0.00$ GeV GeV2	4.75	1.89	0.023	0.064
$0.05$ GeV GeV2	4.75	1.92	0.074	0.085
$0.10$ GeV GeV2	4.71	1.94	0.145	0.108
$0.20$ GeV GeV2	4.72	2.00	0.281	0.147

Table 2. Table 2: Parametrization of the helicity-amplitude ratios and parameter values extracted from the fit. The combinations of the helicity-amplitude ratios ξ 𝜉 \xi and ζ 𝜁 \zeta are defined in Eq. ( 64 ). An additional scale uncertainty of 8% originating from the uncertainty on the target polarization is present for the ratios t λ V ​ λ γ ( 2 ) subscript superscript 𝑡 2 subscript 𝜆 𝑉 subscript 𝜆 𝛾 t^{(2)}_{\lambda_{V}\lambda_{\gamma}} , u λ V ​ λ γ ( 2 ) subscript superscript 𝑢 2 subscript 𝜆 𝑉 subscript 𝜆 𝛾 u^{(2)}_{\lambda_{V}\lambda_{\gamma}} , ξ 𝜉 \xi and ζ 𝜁 \zeta , but not shown. An extra scale uncertainty of 2% originating from the uncertainty on the beam polarization is present for the ratios Im ​ { t λ V ​ λ γ ( 1 ) } Im subscript superscript 𝑡 1 subscript 𝜆 𝑉 subscript 𝜆 𝛾 {\rm{Im}}\{t^{(1)}_{\lambda_{V}\lambda_{\gamma}}\} , Re ​ { t λ V ​ λ γ ( 2 ) } Re subscript superscript 𝑡 2 subscript 𝜆 𝑉 subscript 𝜆 𝛾 {\rm{Re}}\{t^{(2)}_{\lambda_{V}\lambda_{\gamma}}\} , Re ​ { u λ V ​ λ γ ( 2 ) } Re subscript superscript 𝑢 2 subscript 𝜆 𝑉 subscript 𝜆 𝛾 {\rm{Re}}\{u^{(2)}_{\lambda_{V}\lambda_{\gamma}}\} , and ζ 𝜁 \zeta , but also not shown.The correlations between the 25 parameters are listed in Table  3 in Appendix  C .

Parametrization	value of parameter	statistical uncertainty	total uncertainty
$\mathrm{Re}\left\{t_{11}^{\left(1\right)}\right\}=b_1/Q$	$b_1=1.145$ GeV	0.033 GeV	0.081 GeV
$\left|u_{11}^{\left(1\right)}\right|=b_2$	$b_2=0.333$	0.016	0.088
$\mathrm{Re}\left\{u_{11}^{\left(2\right)}\right\}=b_3$	$b_3=-0.074$	0.036	0.054
$\mathrm{Im}\left\{u_{11}^{\left(2\right)}\right\}=b_4$	$b_4=0.080$	0.022	0.037
$\xi =b_5$	$b_5=-0.055$	0.027	0.029
$\zeta =b_6$	$b_6=-0.013$	0.033	0.044
$\mathrm{Im}\left\{t_{00}^{\left(2\right)}\right\}=b_7$	$b_7=0.040$	0.025	0.030
$\mathrm{Re}\left\{t_{01}^{\left(1\right)}\right\}=b_8\sqrt{-t^{\prime }}$	$b_8=0.471$ GeV-1	0.033 GeV-1	0.075 GeV-1
$\mathrm{Im}\left\{t_{01}^{\left(1\right)}\right\}=b_9\frac{\sqrt{-t^{\prime }}}{Q}$	$b_9=0.307$	0.148	0.354
$\mathrm{Re}\left\{t_{01}^{\left(2\right)}\right\}=b_{10}$	$b_{10}=-0.074$	0.060	0.080
$\mathrm{Im}\left\{t_{01}^{\left(2\right)}\right\}=b_{11}$	$b_{11}=-0.067$	0.026	0.036
$\mathrm{Re}\left\{u_{01}^{\left(2\right)}\right\}=b_{12}$	$b_{12}=0.032$	0.060	0.072
$\mathrm{Im}\left\{u_{01}^{\left(2\right)}\right\}=b_{13}$	$b_{13}=0.030$	0.026	0.033
$\mathrm{Re}\left\{t_{10}^{\left(1\right)}\right\}=b_{14}\sqrt{-t^{\prime }}$	$b_{14}=-0.025$ GeV-1	0.034 GeV-1	0.063 GeV-1
$\mathrm{Im}\left\{t_{10}^{\left(1\right)}\right\}=b_{15}\sqrt{-t^{\prime }}$	$b_{15}=0.080$ GeV-1	0.063 GeV-1	0.118 GeV-1
$\mathrm{Re}\left\{t_{10}^{\left(2\right)}\right\}=b_{16}$	$b_{16}=-0.038$	0.026	0.030
$\mathrm{Im}\left\{t_{10}^{\left(2\right)}\right\}=b_{17}$	$b_{17}=0.012$	0.018	0.019
$\mathrm{Re}\left\{u_{10}^{\left(2\right)}\right\}=b_{18}$	$b_{18}=-0.023$	0.030	0.039
$\mathrm{Im}\left\{u_{10}^{\left(2\right)}\right\}=b_{19}$	$b_{19}=-0.045$	0.018	0.026
$\mathrm{Re}\left\{t_{1-1}^{\left(1\right)}\right\}=b_{20}\frac{\left(-t^{\prime }\right)}{Q}$	$b_{20}=-0.008$ GeV-1	0.096 GeV-1	0.212 GeV-1
$\mathrm{Im}\left\{t_{1-1}^{\left(1\right)}\right\}=b_{21}\frac{\left(-t^{\prime }\right)}{Q}$	$b_{21}=-0.577$ GeV-1	0.196 GeV-1	0.428 GeV-1
$\mathrm{Re}\left\{t_{1-1}^{\left(2\right)}\right\}=b_{22}$	$b_{22}=0.059$	0.036	0.047
$\mathrm{Im}\left\{t_{1-1}^{\left(2\right)}\right\}=b_{23}$	$b_{23}=0.020$	0.022	0.026
$\mathrm{Re}\left\{u_{1-1}^{\left(2\right)}\right\}=b_{24}$	$b_{24}=-0.047$	0.035	0.039
$\mathrm{Im}\left\{u_{1-1}^{\left(2\right)}\right\}=b_{25}$	$b_{25}=0.007$	0.022	0.029

Table 3. Table 3: The correlation matrix for the 25 parameters defined in Table  2 .

	$b_1$	$b_2$	$b_3$	$b_4$	$b_5$	$b_6$	$b_7$	$b_8$	$b_9$	$b_{10}$	$b_{11}$	$b_{12}$	$b_{13}$	$b_{14}$	$b_{15}$	$b_{16}$	$b_{17}$	$b_{18}$	$b_{19}$	$b_{20}$	$b_{21}$	$b_{22}$	$b_{23}$	$b_{24}$	$b_{25}$
$b_1$	1.00
$b_2$	0.08	1.00
$b_3$	-0.01	-0.03	1.00
$b_4$	0.05	0.02	-0.20	1.00
$b_5$	0.07	0.00	-0.16	-0.19	1.00
$b_6$	-0.05	0.02	-0.19	-0.03	-0.20	1.00
$b_7$	-0.05	0.02	0.07	0.08	-0.39	-0.18	1.00
$b_8$	-0.01	0.10	0.01	0.01	0.01	0.03	0.01	1.00
$b_9$	0.23	-0.05	-0.03	-0.02	-0.02	-0.09	-0.06	-0.21	1.00
$b_{10}$	0.00	0.04	-0.03	0.02	0.04	-0.19	0.19	0.02	-0.01	1.00
$b_{11}$	0.00	0.02	-0.06	0.08	-0.02	-0.02	0.04	0.04	0.00	0.05	1.00
$b_{12}$	-0.05	-0.02	-0.01	0.12	0.04	-0.04	-0.02	0.01	-0.13	-0.17	0.08	1.00
$b_{13}$	-0.02	-0.01	0.10	-0.09	-0.07	-0.07	0.13	0.01	-0.05	0.08	-0.33	0.03	1.00
$b_{14}$	-0.22	0.13	0.05	0.03	-0.01	0.06	0.05	-0.11	-0.57	-0.02	-0.02	0.09	0.05	1.00
$b_{15}$	0.27	-0.11	-0.07	-0.03	0.00	-0.05	-0.06	-0.14	0.62	0.03	0.01	-0.10	-0.06	-0.71	1.00
$b_{16}$	0.05	-0.05	0.02	0.05	-0.08	-0.13	0.17	-0.02	0.12	0.46	0.15	-0.26	0.06	-0.16	0.21	1.00
$b_{17}$	0.06	-0.01	0.03	0.04	-0.07	-0.11	0.13	-0.02	0.08	0.49	-0.14	-0.41	-0.03	-0.07	0.12	0.42	1.00
$b_{18}$	0.00	0.04	0.23	0.04	0.09	-0.10	-0.02	-0.01	0.05	0.31	-0.06	-0.45	-0.11	0.00	-0.02	0.17	0.37	1.00
$b_{19}$	0.01	-0.03	0.10	0.07	0.02	-0.15	0.16	0.01	0.02	0.48	0.06	-0.44	0.15	-0.01	0.01	0.40	0.35	0.47	1.00
$b_{20}$	0.07	-0.10	-0.03	-0.03	0.02	-0.06	-0.03	-0.08	0.20	-0.02	-0.02	-0.02	-0.03	-0.23	0.28	0.07	0.02	-0.05	-0.04	1.00
$b_{21}$	-0.07	0.13	0.06	0.03	-0.04	0.10	0.00	0.15	-0.05	0.01	0.04	-0.01	0.01	0.14	-0.27	-0.07	-0.03	0.10	0.04	-0.53	1.00
$b_{22}$	-0.02	0.03	-0.01	-0.27	-0.23	0.34	-0.06	0.01	-0.03	-0.01	-0.09	-0.08	0.03	0.05	-0.07	-0.09	0.04	-0.06	-0.07	-0.07	0.13	1.00
$b_{23}$	0.02	-0.03	-0.25	0.40	-0.17	0.09	-0.11	-0.02	0.08	-0.07	0.01	0.04	-0.09	-0.03	0.03	0.03	-0.03	-0.01	-0.09	0.00	0.07	0.02	1.00
$b_{24}$	-0.06	-0.04	-0.09	-0.22	-0.16	0.31	-0.06	-0.01	-0.03	-0.03	0.04	-0.06	0.13	0.04	-0.02	0.02	-0.06	-0.19	-0.04	-0.02	0.03	0.08	-0.02	1.00
$b_{25}$	0.00	-0.02	-0.33	0.40	-0.21	0.14	-0.09	-0.02	0.04	0.00	0.01	0.10	-0.07	-0.03	0.04	0.05	0.01	-0.12	-0.13	0.00	0.00	0.00	0.37	0.14	1.00
	$b_1$	$b_2$	$b_3$	$b_4$	$b_5$	$b_6$	$b_7$	$b_8$	$b_9$	$b_{10}$	$b_{11}$	$b_{12}$	$b_{13}$	$b_{14}$	$b_{15}$	$b_{16}$	$b_{17}$	$b_{18}$	$b_{19}$	$b_{20}$	$b_{21}$	$b_{22}$	$b_{23}$	$b_{24}$	$b_{25}$

Equations

$${\gamma }^{\ast }({\lambda }_{\gamma })+N({\lambda }_N)\to V({\lambda }_V)+N({\lambda }_N^{\prime }),$$

$$e+N\to e+{\rho }^0+N,$$

$${\rho }^0\to {\pi }^{+}+{\pi }^{-}.$$

$$W^2=(p+q{)}^2=M^2+2M\nu -Q^2,$$

$$\nu =\frac{p\cdot q}{M}\stackrel{lab.}{=}E-E^{\prime },$$

$$t=(q-v{)}^2,$$

$$ϵ=\frac{1-y-\frac{Q^2}{4E^2}}{1-y+\frac{y^2}{2}+\frac{Q^2}{4E^2}},$$

$$\Delta E=\frac{M_X^2-M^2}{2M}\stackrel{lab.}{=}{E}_V-(E_{{\pi }^{+}}+E_{{\pi }^{-}}),$$

$$F_{{\lambda }_V{\lambda }_N^{\prime }{\lambda }_{\gamma }{\lambda }_N}=(-1{)}^{{\lambda }_{\gamma }}⟨v{\lambda }_V{p}^{\prime }{\lambda }_N^{\prime }|J^{\kappa }|p{\lambda }_N⟩e_{\kappa }^{({\lambda }_{\gamma })},$$

$$e^{(\pm 1)}$$

$$e^{(0)}$$

$$q_0=\frac{M\nu -Q^2}{W},q_Z=\frac{M\sqrt{{\nu }^2+Q^2}}{W}.$$

$$F_{{\lambda }_V{\lambda }_N^{\prime }{\lambda }_{\gamma }{\lambda }_N}=T_{{\lambda }_V{\lambda }_N^{\prime }{\lambda }_{\gamma }{\lambda }_N}+U_{{\lambda }_V{\lambda }_N^{\prime }{\lambda }_{\gamma }{\lambda }_N},$$

$$T_{{\lambda }_V{\lambda }_N^{\prime }{\lambda }_{\gamma }{\lambda }_N}$$

$$=\frac{1}{2}[F_{{\lambda }_V{\lambda }_N^{\prime }{\lambda }_{\gamma }{\lambda }_N}+(-1{)}^{{\lambda }_N-{\lambda }_N^{\prime }}{F}_{{\lambda }_V-{\lambda }_N^{\prime }{\lambda }_{\gamma }-{\lambda }_N}],$$

$$U_{{\lambda }_V{\lambda }_N^{\prime }{\lambda }_{\gamma }{\lambda }_N}$$

$$=\frac{1}{2}[F_{{\lambda }_V{\lambda }_N^{\prime }{\lambda }_{\gamma }{\lambda }_N}-(-1{)}^{{\lambda }_N-{\lambda }_N^{\prime }}{F}_{{\lambda }_V-{\lambda }_N^{\prime }{\lambda }_{\gamma }-{\lambda }_N}].$$

$$T_{{\lambda }_V{\lambda }_N^{\prime }{\lambda }_{\gamma }{\lambda }_N}$$

$$U_{{\lambda }_V{\lambda }_N^{\prime }{\lambda }_{\gamma }{\lambda }_N}$$

$$T_{{\lambda }_V{\lambda }_{\gamma }}^{(1)}$$

$$U_{{\lambda }_V{\lambda }_{\gamma }}^{(1)}$$

$$T_{{\lambda }_V{\lambda }_{\gamma }}^{(2)}$$

$$U_{{\lambda }_V{\lambda }_{\gamma }}^{(2)}$$

$$T_{{\lambda }_V{\lambda }_{\gamma }}^{(j)}$$

$$U_{{\lambda }_V{\lambda }_{\gamma }}^{(j)}$$

$$\displaystyle F_{\lambda_{V}\lambda^{\prime}_{N}\lambda_{\gamma}\lambda_{N}}=\sum_{k=0}^{\infty}c^{\lambda_{V}\lambda^{\prime}_{N}\lambda_{\gamma}\lambda_{N}}_{k}\Bigl{(}\frac{\sqrt{-t^{\prime}}}{M_{h}}\Bigr{)}^{s+2k},$$

$$t_{11}^{(1)},t_{01}^{(2)},t_{10}^{(2)},u_{11}^{(1)},u_{01}^{(2)},u_{10}^{(2)}$$

$$t_{{\lambda }_V{\lambda }_{\gamma }}^{(j)}$$

$$u_{{\lambda }_V{\lambda }_{\gamma }}^{(j)}$$

$$t_{00}^{(2)},t_{11}^{(2)},t_{01}^{(1)},t_{10}^{(1)},u_{11}^{(2)},u_{1-1}^{(2)},u_{01}^{(1)},u_{10}^{(1)}.$$

$$t_{1-1}^{(1)},u_{1-1}^{(1)}$$

$${\varrho }_{{\lambda }_{\gamma }{\lambda }_{\gamma }^{\prime }}^{U+L}={\varrho }_{{\lambda }_{\gamma }{\lambda }_{\gamma }^{\prime }}^U+P_B{\varrho }_{{\lambda }_{\gamma }{\lambda }_{\gamma }^{\prime }}^L,$$

$${\varrho }_{{\lambda }_{\gamma }{\lambda }_{\gamma }^{\prime }}^U(ϵ,\Phi )$$

$$\displaystyle=\frac{1}{2}\!\left(\!\!\begin{array}[]{ccc}1&\sqrt{\epsilon(1+\epsilon)}e^{-i\Phi}&-\epsilon e^{-2i\Phi}\\ \sqrt{\epsilon(1+\epsilon)}e^{i\Phi}&2\epsilon&-\sqrt{\epsilon(1+\epsilon)}e^{-i\Phi}\\ -\epsilon e^{2i\Phi}&-\sqrt{\epsilon(1+\epsilon)}e^{i\Phi}&1\\ \end{array}\!\right)\!,$$

$${\varrho }_{{\lambda }_{\gamma }{\lambda }_{\gamma }^{\prime }}^L(ϵ,\Phi )$$

$$\displaystyle=\frac{\sqrt{1-\epsilon}}{2}\left(\begin{array}[]{ccc}\sqrt{1+\epsilon}&\sqrt{\epsilon}e^{-i\Phi}&0\\ \sqrt{\epsilon}e^{i\Phi}&0&\sqrt{\epsilon}e^{-i\Phi}\\ 0&\sqrt{\epsilon}e^{i\Phi}&-\sqrt{1+\epsilon}\\ \end{array}\right),$$

Full text

1 Introduction

Exclusive electroproduction of vector mesons ($V$) on nucleons ($N$) has been investigated for many decades (see, for instance, Refs. Bauer ; INS ). Originally, the reaction mechanism was of primary interest, but now it has become apparent that this process also offers the possibility to study the structure of the nucleon and of the vector meson Bauer ; INS ; MD , especially at large virtuality $Q^{2}$ of the photon exchanged between electron and nucleon. In the one-photon-exchange approximation, all electroproduction observables can be expressed in terms of the virtual-photon spin-density matrix and the matrix elements of the electromagnetic current between quantum states of initial and final hadrons. The latter matrix elements are called helicity amplitudes $F_{\lambda_{V}\lambda^{\prime}_{N}\lambda_{\gamma}\lambda_{N}}$. They describe the process

[TABLE]

where $\gamma^{*}$ denotes the virtual photon and the helicities of the particles are given in parentheses. In the present paper, the helicity amplitudes are defined in the center-of-mass (CM) system of virtual photon and nucleon. The spin-density matrix of the virtual photon is well known from quantum electrodynamics and the spin-density matrix elements (SDMEs) of the produced vector meson, which describe its final spin states, are experimentally accessible. This opens in principle the possibility to extract the helicity amplitudes, as it is detailed below.

The formalism describing SDMEs of the produced vector meson was first presented in Ref. SW for unpolarized targets only, and expressions of SDMEs in terms of helicity amplitudes were also established. The formalism was then extended to the case of polarized targets in Ref. Fraas . An alternative, general formalism for the description of the process in Eq. (1) through SDMEs was presented in Ref. Diehl . In the latter formalism, which is used throughout this paper, the SDMEs describing the production on an unpolarized target are denoted by $u^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$, those describing the production on a longitudinally polarized target are denoted by $l^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$, and those describing the production on a transversely polarized target are denoted by $n^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$ and $s^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$. Here, longitudinal and transverse polarization are defined with respect to the momentum direction of the virtual photon in the CM system of the process in Eq. (1).

The exact expressions for SDMEs SW ; Fraas ; Diehl , which are dimensionless quantities, can be rewritten in terms of ratios of helicity amplitudes. When fitting the experimental angular distribution of the final-state particles, either the SDMEs or alternatively the amplitude ratios can be considered as independent free parameters. The first fit method is referred to as the “SDME method” in the rest of this paper, while the second one is referred to as the “amplitude method”.

Exclusive meson production in hard lepton-nucleon scattering was shown to offer the possibility of constraining generalized parton distributions (GPDs), which provide correlated information on transverse-spatial and fractional-longitudinal-momentum distributions of partons in the nucleon (see Refs. MD ; gpd1 ; VGG ; GPV ; GPRV ; BR and references therein). Vector-meson-production amplitudes contain various linear combinations of GPDs for quarks of various flavors and for gluons. In particular, exclusive $\rho^{0}$ production on an unpolarized target is sensitive to the nucleon-helicity-non-flip GPD $H$, while exclusive $\rho^{0}$ production on a transversely polarized target is sensitive to the nucleon-helicity-flip GPD $E$, as well. Through the Ji relation gpd2 , the sum of both GPDs $H$ and $E$ is related to the parton total angular momentum. Access to GPDs relies on the factorization property of the process amplitude, i.e., the amplitude can be written as a convolution of GPDs and vector-meson distribution amplitudes, which are both non-perturbative quantities, and amplitudes of hard partonic subprocesses, which are calculable within the frameworks of perturbative quantum chromodynamics (pQCD) and quantum electrodynamics.

For spin-1 particles, longitudinal (transverse) polarization is assigned by convention to the states with helicity $\lambda=0$ ($\lambda=\pm 1$). The helicity amplitudes $F_{0\frac{1}{2}0\pm\frac{1}{2}}$ describe the transition of a longitudinally ($L$) polarized virtual photon to a longitudinally polarized vector meson, $\gamma^{*}_{L}\to V_{L}$, and dominate at large photon virtuality $Q^{2}$. Although factorization was rigorously proven CFS only for these amplitudes, it was assumed in Refs. golos2 ; golos3 that factorization also holds for the amplitudes $F_{1\frac{1}{2}1\pm\frac{1}{2}}$ and $F_{0\frac{1}{2}1\pm\frac{1}{2}}$, which describe the transition from a transversely polarized virtual photon to a transversely polarized meson, $\gamma^{*}_{T}\to V_{T}$, and a longitudinally polarized meson, $\gamma^{*}_{T}\to V_{L}$, respectively. The agreement found between certain calculated SDMEs and those extracted from HERMES DC-24 , ZEUS ZEUS2 , and H1 H1-amp data supports this assumption. In general, the differential and total cross sections for $\rho^{0}$-meson production by virtual photons are reasonably well described in the GPD-based approach of Refs. golos2 ; golos3 , not only at the high energies of the HERA collider experiments ZEUS3 ; ZEUS4 ; H1-1 ; H1 ; H1-2 , but also at intermediate energies covered by the fixed-target experiments E665 E665 and HERMES rho-xsec .

The real parts of the amplitude ratios in $\rho^{0}$ and $\phi$ meson electroproduction on the proton were first studied by the H1 experiment H1-amp at the HERA collider. In the HERMES experiment DC-84 , $\rho^{0}$-meson production on unpolarized protons and deuterons was investigated. Both real and imaginary parts of the ratios of amplitudes without nucleon helicity flip were extracted at HERMES using a longitudinally polarized electron or positron beam. The results of the analysis of $\rho^{0}$-meson and $\omega$-meson production on the unpolarized targets at HERMES using the SDME method were published in Ref. DC-24 and Ref. DC-95 , respectively. The SDMEs for the electroproduction on transversely polarized protons were published in Ref. DC-71 . In this paper, the work of Ref. DC-71 is continued. Ratios of $\rho^{0}$ helicity amplitudes with respect to the amplitude $F_{0\frac{1}{2}0\frac{1}{2}}$ are extracted from HERMES data collected with longitudinally polarized electron and positron beams scattered off transversely polarized protons. The amplitude ratios that require measurements with a transversely polarized target are reported for the first time in this paper.

At fixed $Q^{2}$ and CM energy $W$ in the $\gamma^{*}N$ system the cross section $\text{d}\sigma/\text{d}t$, which is differential in the Mandelstam variable $t$, contains the linear combination of squares of all helicity amplitudes. Including $\text{d}\sigma/\text{d}t$ in an amplitude analysis of all beam and target-polarization states would allow the extraction of the moduli of all amplitudes and of the phase differences between them, while the common phase would remain undetermined.

The amplitude ratios measured at HERMES, as described in this paper, will also be compared to those evaluated within the GPD-based handbag approach by Goloskokov and Kroll golos2 ; golos3 , hereafter referred to as “GK model”.

The paper is organized as follows. In Section 2, the theoretical formalism is introduced. Section 3 briefly describes the experimental setup and specifies the applied data selection. The extraction procedure of the amplitude ratios is treated in Section 4. The obtained results are discussed in Section 5. Summary and conclusions are given in Section 6.

2 Formalism

2.1 Kinematics

The process under investigation is

[TABLE]

with

[TABLE]

In accordance with the notation of Ref. DC-24 , the kinematic variables of the process under study are defined as follows. The four-momenta of the incident and outgoing leptons are denoted by $k$ and $k^{\prime}$, respectively, the difference of which defines the four-momentum $q=k-k^{\prime}$ of the virtual photon $\gamma^{*}$. The photon virtuality $Q^{2}=-q^{2}$ is positive in leptoproduction. The squared invariant mass of the photon-nucleon system is given by

[TABLE]

with $M$ the nucleon mass, $p$ the four-momentum of the incident nucleon and

[TABLE]

the energy transfer from the incoming lepton to the virtual photon in the target rest frame (“lab.” frame). Here, $E$ ($E^{\prime}$) is the energy of the incident (scattered) lepton.

The Mandelstam variable $t$ is defined by the relation

[TABLE]

where $v$ is the four-momentum of the $\rho^{0}$ meson, equal to $p_{\pi^{+}}+p_{\pi^{-}}$, the sum of the $\pi^{+}$ and $\pi^{-}$ four-momenta. The variables $t$, $t_{0}$, and $t^{\prime}=t-t_{0}$ are always negative, where $-t_{0}$ is the minimal value of $-t$ for given values of $Q^{2}$, $W$, and the $\rho^{0}$-meson mass $M_{V}$. At small values of $-t^{\prime}$, the approximation $-t^{\prime}\approx v_{T}^{2}$ holds, where $v_{T}$ is the transverse momentum of the $\rho^{0}$ meson with respect to the direction of the virtual photon in the $\gamma^{*}N$ CM system.

The variable $\epsilon$ represents the ratio of fluxes of longitudinally and transversely polarized virtual photons:

[TABLE]

with $y=p\cdot q/p\cdot k\stackrel{{\scriptstyle lab.}}{{=}}\nu/E$.

The “exclusivity” of $\rho^{0}$ production in the process in Eq. (2) is characterized by the missing energy

[TABLE]

where $M_{X}=\sqrt{(k-k^{\prime}+p-p_{\pi^{+}}-p_{\pi^{-}})^{2}}$ is the reconstructed invariant mass of the undetected hadronic system (missing mass), $E_{V}=\nu+t/(2M)$ is the energy of the exclusively produced $\rho^{0}$ meson, and $(E_{\pi^{+}}+E_{\pi^{-}})$ is the sum of the energies of the two detected pions in the target rest frame.

2.2 Definition of angles and coordinate systems

The angles used for the description of the process are defined in the same way as in Ref. DC-24 , according to Ref. DESY2 , and are presented in Fig. 1. According to Ref. SW , the right-handed “hadronic CM system” of coordinates $XYZ$ of virtual photon and target nucleon is defined such that the $Z$-axis is aligned along the virtual-photon three-momentum $\vec{q}$ and the $Y$-axis is parallel to $\vec{q}\times\vec{v}$, where $\vec{v}$ is the $\rho^{0}$-meson three-momentum. The angle $\Phi$ is the angle between the $\rho^{0}$-meson production plane ($XZ$ plane, which coincides with the nucleon scattering plane) and the lepton scattering plane in the CM system. The angles $\theta$ and $\phi$ are defined in the right-handed $xyz$ system of coordinates (see Fig. 1) that represents the $\rho^{0}$-meson rest frame. The $y$ axis coincides with the $Y$ axis. The angle $\theta$ is the polar angle of the decay $\pi^{+}$-meson three-momentum with respect to the $z$ axis, where the latter is aligned opposite to the direction of the momentum of the outgoing nucleon. The azimuthal angle of the $\pi^{+}$ momentum with respect to the $\rho^{0}$-meson production plane in the CM system is denoted $\phi$. In the HERMES experiment, the vector $\vec{P}_{T}$ of the target polarization is orthogonal to the beam direction. The angle between the directions of the transverse part (with respect to the beam) of the scattered electron momentum and $\vec{P}_{T}$ is denoted by $\Psi$ and is defined in the target rest frame.

2.3 Natural and unnatural-parity-exchange helicity amplitudes

The helicity amplitudes $F_{\lambda_{V}\lambda^{\prime}_{N}\lambda_{\gamma}\lambda_{N}}$ describing exclusive $\rho^{0}$-meson production by the virtual photon are here defined in the hadronic CM system SW . These helicity amplitudes can be expressed as scalar products of the matrix element of the electromagnetic current vector $J^{\kappa}$ and the virtual-photon polarization vector $e_{\kappa}^{(\lambda_{\gamma})}$:

[TABLE]

where a summation over the Lorentz index $\kappa$ is performed. Here, $e_{\kappa}^{(\pm 1)}$ and $e_{\kappa}^{(0)}$ indicate transverse and longitudinal polarization of the virtual photon in the CM system, respectively:

[TABLE]

where $q_{0}$ and $q_{Z}$ are the energy and the $Z$ component of the three-momentum of the virtual photon in the CM system given by

[TABLE]

The ket vector $|p\lambda_{N}\rangle$ corresponds to the initial nucleon and the bra vector $\langle v\lambda_{V}p^{\prime}\lambda^{\prime}_{N}|$ represents the final state consisting of a $\rho^{0}$ meson and the scattered nucleon.

Any helicity amplitude $F_{\lambda_{V}\lambda^{\prime}_{N}\lambda_{\gamma}\lambda_{N}}$ can be decomposed into the sum of an amplitude $T_{\lambda_{V}\lambda^{\prime}_{N}\lambda_{\gamma}\lambda_{N}}$ for natural-parity exchange (NPE) and an amplitude $U_{\lambda_{V}\lambda^{\prime}_{N}\lambda_{\gamma}\lambda_{N}}$ for unnatural-parity exchange (UPE) SW ; Fraas ; Diehl :

[TABLE]

where the NPE and UPE amplitudes are defined as

[TABLE]

These amplitudes by their definition obey the symmetry relations

[TABLE]

Equations (15) and (16) permit the introduction of the following abbreviated notations for the amplitudes:

[TABLE]

which are diagonal with respect to the nucleon helicity ($\lambda_{N}=\lambda^{\prime}_{N}$), and

[TABLE]

hold for the amplitudes with nucleon-helicity flip. Due to parity conservation (see, e.g., Refs. SW ; Diehl ), the amplitudes $T^{(j)}_{\lambda_{V}\lambda_{\gamma}}$ and $U^{(j)}_{\lambda_{V}\lambda_{\gamma}}$ obey the symmetry relations for $j=1,\;2$:

[TABLE]

This implies that there is a linear dependence between certain amplitudes. Therefore, if some property of the amplitude $T^{(j)}_{\lambda_{V}\lambda_{\gamma}}$ ($U^{(j)}_{\lambda_{V}\lambda_{\gamma}}$) is established for some particular $\lambda_{V}$, $\lambda_{\gamma}$, and $j$, the amplitude $T^{(j)}_{-\lambda_{V}-\lambda_{\gamma}}$ ($U^{(j)}_{-\lambda_{V}-\lambda_{\gamma}}$) has the same property.

There are three important consequences of the symmetry relations (21) and (22) SW ; Diehl :

The number of linearly independent NPE (UPE) amplitudes is equal to 10 (8); 2. 2.

No UPE amplitude exists for the transition $\gamma_{L}\to\rho^{0}_{L}$, so that in particular $F_{0\frac{1}{2}0\frac{1}{2}}=T_{0\frac{1}{2}0\frac{1}{2}}\equiv T^{(1)}_{00}$; 3. 3.

For unpolarized targets there is no interference between NPE and UPE amplitudes SW ; Diehl .

At small values of $t$, in Regge phenomenology IW ; KS the exchange of a single natural-parity reggeon, i.e., with parity $P=(-1)^{J}$, such as a pomeron or secondary reggeons $\rho$, $f_{2}$, $a_{2}$, …, contributes only to the NPE amplitudes. The exchange of a single unnatural-parity reggeon, i.e., with parity $P=-(-1)^{J}$, such as $\pi$, $a_{1}$, $b_{1}$,…, contributes only to the UPE amplitudes CSM .

2.4 Asymptotic behavior of amplitudes at small $|t^{\prime}|$

Considering only the behavior of a helicity amplitude at $-t^{\prime}\to 0$, its magnitude relative to other amplitudes can be investigated by using the following parametrization

[TABLE]

where $s=|(\lambda_{V}-\lambda^{\prime}_{N})-(\lambda_{\gamma}-\lambda_{N})|$ and $M_{h}$ represents the typical hadronic mass of the order of $1$ GeV. If $c^{\lambda_{V}\lambda^{\prime}_{N}\lambda_{\gamma}\lambda_{N}}_{0}$ vanishes, then the power series with respect to $t^{\prime}$ starts with a term $\propto(\sqrt{-t^{\prime}}/M_{h})^{s+2k}$ with $k\geq 1$. The asymptotic relations for $T_{\lambda_{V}\lambda^{\prime}_{N}\lambda_{\gamma}\lambda_{N}}$ and $U_{\lambda_{V}\lambda^{\prime}_{N}\lambda_{\gamma}\lambda_{N}}$, and hence for $T^{(1)}_{\lambda_{V}\lambda_{\gamma}}$, $T^{(2)}_{\lambda_{V}\lambda_{\gamma}}$, $U^{(1)}_{\lambda_{V}\lambda_{\gamma}}$, and $U^{(2)}_{\lambda_{V}\lambda_{\gamma}}$, follow from Eq. (23) and Eqs. (13) and (14).

From the asymptotic behavior, the NPE amplitudes $T^{(1)}_{00}$, $T^{(1)}_{11}$, $T^{(2)}_{01}$, $T^{(2)}_{10}$ and the UPE amplitudes $U^{(1)}_{11}$, $U^{(2)}_{01}$, $U^{(2)}_{10}$ are proportional to $(-t^{\prime})^{0}$ at $|t^{\prime}|\to 0$. Therefore, the amplitude ratios

[TABLE]

defined by the relations

[TABLE]

for $j=1,\;2$ can be non-zero for $|t^{\prime}|\rightarrow 0$. These ratios are expected to attain their largest values at small $-t^{\prime}$.

Similarly, one can conclude from Eq. (23) that the amplitude ratios proportional to $\sqrt{-t^{\prime}}/M_{h}$ at small $-t^{\prime}$ are

[TABLE]

However, if for a UPE amplitude appearing in Eq. (27) the pion exchange in the $t$ channel is significant, the typical scale for $t^{\prime}$ is about $m_{\pi}^{2}$. Therefore, this amplitude can be of the order of the dominant amplitude $T^{(1)}_{00}$ at $-t^{\prime}\sim m_{\pi}^{2}$. Hence some amplitude ratios from Eq. (27) can be of the same order of magnitude as those in Eq. (24) at $-t^{\prime}\sim m_{\pi}^{2}$.

The smallest amplitudes at $|t^{\prime}|\to 0$ are the double spin-flip amplitudes $T^{(1)}_{1-1}$ and $U^{(1)}_{1-1}$ and the amplitude ratios

[TABLE]

are proportional to $-t^{\prime}$.

2.5 Spin-density matrix of the virtual photon

The spin-density matrix of the virtual photon, normalized such that the flux of transversely polarized photons is equal to unity, embodies the unpolarized ($U$) and polarized ($L$) matrices SW :

[TABLE]

with

[TABLE]

and where $P_{B}$ is the longitudinal polarization of the beam and $\epsilon$ is defined by Eq. (7). In the above formulas, the spin-density matrices of the virtual photon are defined in the hadronic CM system.

2.6 Spin-density matrix of the initial nucleon

The angle between the three-momenta of the initial electron and the virtual photon in the target rest frame $\theta_{\gamma}$ can be calculated as

[TABLE]

The “hadronic rest system”, in which the initial proton is at rest, has the $\hat{X}$, $\hat{Y}$, and $\hat{Z}$ axes parallel to the $X$, $Y$, and $Z$ axes of the hadronic CM system, defined in Section 2.2 and shown in Fig. 1. The components of the target polarization vector $\vec{P}$ in the hadronic rest system are

[TABLE]

The spin-density matrix of the initial nucleon in the helicity representation can be written in this system as

[TABLE]

Since the nucleon spin is anti-parallel to the $\hat{Z}$ axis for $\lambda_{N}=\frac{1}{2}$ and parallel to the $\hat{Z}$ axis for $\lambda_{N}=-\frac{1}{2}$, Eq. (43) for the spin-density matrix of the initial nucleon does not coincide with the standard formula $(I+\vec{P}\vec{\sigma})/2$ for the quantization axis aligned along the $\hat{Z}$ axis. Here, $\vec{\sigma}=(\sigma_{x},\sigma_{y},\sigma_{z})$ are the Pauli matrices and $\vec{P}$ is the polarization vector in the target rest frame. The hadronic CM system can be obtained from the hadronic rest system by a boost along the virtual-photon three-momentum, which is antiparallel to the proton three-momentum. Since the value of the proton helicity is invariant under this boost the spin-density matrix is also boost invariant, hence it is given by Eq. (43) in the hadronic CM system.

2.7 Spin-density matrix of the $\rho^{0}$ meson

The spin-density matrix $\rho_{\lambda_{V}\lambda^{\prime}_{V}}$ of the produced $\rho^{0}$ meson is related through the von Neumann formula to those of the virtual photon, $\varrho^{U+L}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$, and the nucleon, $\tau_{\lambda_{N}\lambda^{\prime}_{N}}$:

[TABLE]

where the sum runs over $\lambda_{\gamma}$, $\lambda^{\prime}_{\gamma}$, $\lambda_{N}$, $\lambda^{\prime}_{N}$, and $\mu_{N}$. The normalization factor is given by

[TABLE]

with

[TABLE]

Equation (47) is obtained by using the symmetry relations (21) and (22).

2.8 SDMEs in the Diehl representation

The spin-density matrix elements calculated below are defined in accordance with Ref. Diehl as

[TABLE]

Here, the NPE and UPE helicity amplitudes are defined in Eqs. (13) and (14), and $\mathcal{N}$ denotes the normalization factor given by Eqs. (45-47).

2.9 Angular distribution of decay pions

The angular distribution $\mathcal{W}(\Phi,\Psi,\theta,\phi)$ of the pions from the $\rho^{0}$-meson decay in Eq. (3), which are produced in the process (2), is related to the spin-density matrix $\rho_{\lambda_{V}\lambda^{\prime}_{V}}$ through

[TABLE]

where $Y_{1\lambda_{V}}$ is the spherical function. The phases of $Y_{1\lambda_{V}}$ are chosen as in Ref. Diehl :

[TABLE]

These phases determine the phases of the extracted helicity amplitudes. Equation (52) shows explicitly the dependence of the angular distribution $\mathcal{W}$ on $\theta$ and $\phi$, while the dependences on $\Phi$ and $\Psi$ are hidden in the kinematic dependences of the spin-density matrix $\rho_{\lambda_{V}\lambda^{\prime}_{V}}$. The angular distribution depends on the kinematic variables $W$, $Q^{2}$, and $-t^{\prime}$ through the dependence of the helicity amplitudes on these variables in Eq. (44). For simplicity of notation, these dependences are omitted throughout the paper.

Since $\tau_{\lambda_{N}\lambda^{\prime}_{N}}$ depends linearly on the nucleon transverse polarization $P_{T}$, and $\varrho_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$ is a linear function of the beam polarization, the formula for the angular distribution contains four terms:

[TABLE]

Note that the angular dependent functions $\mathcal{W}_{m}$ for any $m$ in Eq. (53) are themselves independent of $P_{B}$ and $P_{T}$.

It can easily be shown that the angular distribution $\mathcal{W}(\Phi,\Psi,\theta,\phi)$ cannot be negative for any set of values of the complex amplitudes $F_{\lambda_{V}\mu_{N}\lambda_{\gamma}\lambda_{N}}$, even for unphysical ones. This property is of great importance for the fit procedure. It is worthwhile to note that using the SDME method one faces the problem of a possible negativity of $\mathcal{W}(\Phi,\Psi,\theta,\phi)$ for some angles when SDMEs assume unphysical values. As it is unknown in which region in the multi-dimensional space of SDMEs $\mathcal{W}(\Phi,\Psi,\theta,\phi)$ is not negative, serious problems may appear when applying the maximum-likelihood method. Hence the amplitude method is in that respect more reliable than the SDME method.

Altogether, Eqs. (52), (44-47), (29-37), and (38-43), with the substitutions $T_{00}^{(1)}\rightarrow 1$ and for all other amplitudes $T_{\lambda_{V}\lambda_{\gamma}}^{(j)}\rightarrow t_{\lambda_{V}\lambda_{\gamma}}^{(j)}$ and $U_{\lambda_{V}\lambda_{\gamma}}^{(j)}\rightarrow u_{\lambda_{V}\lambda_{\gamma}}^{(j)}$, constitute a basis for the amplitude method, in which the extracted quantities are the helicity-amplitude ratios.

3 Experiment and data selection

3.1 Experiment

A detailed description of the HERMES experiment can be found in Ref. herspec . The data analyzed in this paper were collected between the years 2002 and 2005. A longitudinally polarized positron or electron beam of $27.6$ GeV was scattered from a pure gaseous, transversely polarized hydrogen target internal to the HERA lepton storage ring. The helicity of the beam was typically reversed every two months. The beam polarization was continuously measured by two Compton polarimeters BAR ; BECK . The average value of the beam polarization for the events used in the analysis is about $\pm 0.3$, with a relative uncertainty of 2%. The target polarization was reversed every 60 s to 180 s HER-POL . The measured mean value of the target polarization is $\langle|P_{T}|\rangle=0.72\pm 0.06$ Ami-th ; Jeron-th .

The HERMES setup included a forward spectrometer herspec , in which the scattered lepton and the produced hadrons were detected within an angular acceptance of $\pm$170 mrad horizontally and $\pm(40-140)$ mrad vertically. The tracking system had a momentum resolution of about $1.5\%$ and an angular resolution of about $1$ mrad. Lepton identification was accomplished using a transition-radiation detector, a preshower scintillator counter, and an electromagnetic calorimeter. The particle-identification system included also a dual-radiator ring-imaging Cherenkov detector RICH to identify hadrons. Combining the responses of the detectors in a likelihood method leads to an average lepton-identification efficiency of 98%, with a hadron contamination of less than 1%.

3.2 Event selection

The event sample used in this analysis is almost the same as that used in Ref. DC-71 . The most important improvement is the application of a new tracking algorithm, which is based on a Kalman filter Kal . For the present analysis, the data are required to fulfill the following criteria:

The longitudinal beam polarization is restricted to the interval $15\%<|P_{B}|<80\%$. 2. 2.

Events with exactly two oppositely charged hadrons and one lepton with the same charge as the beam lepton are selected. All tracks are required to originate from the same vertex. 3. 3.

The scattered lepton has to have an energy larger than 3.5 GeV in order to not introduce effects from varying trigger thresholds. 4. 4.

The two-hadron invariant mass is required to lie around the $\rho^{0}$ mass, i.e., it is required to obey $0.6$ GeV $<M(\pi^{+}\pi^{-})<1.0$ GeV. 5. 5.

The photon virtuality is required to obey 1 GeV${}^{2}<Q^{2}<7$ GeV2. The lower limit is a minimum requirement for the application of pQCD, while the upper one delimits a well defined kinematic phase space. 6. 6.

The $t^{\prime}$ variable is restricted to $-t^{\prime}\leq 0.4$ GeV2 in order to reduce non-exclusive background of the reaction under study. 7. 7.

The invariant mass $W$ is required to obey 3 GeV $<W\leq 6.3$ GeV. The requirement $W>3$ GeV is imposed in order to be outside of the resonance region. The upper constraint delimits a well defined kinematic phase space. 8. 8.

For exclusive $\rho^{0}$-meson production, $\Delta E$ as defined in Eq. (8) must vanish. In the present analysis, taking into consideration the spectrometer resolution, the missing energy has to be in the region $-1.0$ GeV $<\Delta E<0.6$ GeV. This region is referred to as “exclusive region” in the following.

After application of all these constraints, the data sample contains 8741 events. These data are referred to in the following as data in the “entire kinematic region”. The applied requirements do not fully suppress contributions from background processes. The exclusive sample contains contributions from non-resonant $\pi^{+}\pi^{-}$ pair production, which is of the order of $1-2\%$ Ami-th , and from semi-inclusive deep-inelastic scattering (SIDIS) events. The presented results are not corrected for the former process, while a correction is applied for SIDIS background. The uncertainty of the correction for background from SIDIS events is considered to be one of the main contributions to the total systematic uncertainty.

4 Extraction of amplitude ratios

4.1 Fit of the angular distribution

A maximum-likelihood method is used to fit the angular distribution of the $\rho^{0}$-meson decay pions and the scattered lepton. The probability to measure $\rho^{0}$ decay pions in the small angular region $\text{d}\Omega=\text{d}\Phi\text{d}\Psi\sin\theta\text{d}\theta\text{d}\phi$ is proportional to $\mathcal{W}(\Phi,\Psi,\theta,\phi)\mathcal{E}(\Phi,\Psi,\theta,\phi)\text{d}\Omega$. Here, the detector efficiency is denoted by $\mathcal{E}(\Phi,\Psi,\theta,\phi)$. It includes geometric detector acceptance, particle-detection and particle-identification efficiencies as well as track-reconstruction efficiency. In order to become a probability, this expression needs to be normalized to unity:

[TABLE]

Here, we have changed the notation of the angular distribution to $\mathcal{W}(\mathcal{R},P_{B},P_{T},\Phi,\Psi,\theta,\phi)$, as $\mathcal{W}$ depends on the set $\mathcal{R}$ of all amplitude ratios (due to Eqs. (52), (44) and (45-47)) as well as on the beam and target polarization. All factors that are independent of the set of amplitude ratios can be omitted in the expression of the likelihood function. The likelihood function is evaluated in $3\times 4$ ($Q^{2},-t^{\prime}$) cells. Within each of the 12 cells, the detector efficiency is considered to be independent of $Q^{2}$ and $-t^{\prime}$. The lower and upper boundaries of the cells and the mean values of the variables $W$, $Q^{2}$, and $-t^{\prime}$ are presented in Table 1. Within each of these cells, the logarithm of the likelihood function can be written as the sum over all experimental events in this cell ($i=1,2,...,I$):

[TABLE]

Here, $N_{i}$ is the normalization factor for the $i$-th event in the corresponding $(Q^{2},-t^{\prime})$ cell, with

[TABLE]

Using the expression for $\mathcal{W}$ from Eq. (53), the normalization factor can be represented as

[TABLE]

where the functions $\mathcal{K}_{m}$ ($m=1,\;2,\;3,\;4$) are given by the integrals

[TABLE]

In order to evaluate the integrals in Eqs. (58), a dedicated PYTHIA Monte Carlo (MC) pythia simulation was used to generate exclusively produced $\rho^{0}$ events with a uniform angular distribution, hereafter referred to as “uniform exclusive $\rho^{0}$ MC” Ami-th ; DC-24 . The total number of generated events is by a factor of one hundred larger than that of the experimental data. The same event selection requirements are applied to the reconstructed events from the simulation as to the events from experimental data. The integrals in Eqs. (58) can be approximated using the MC as follows:

[TABLE]

Here, $L$ is the total number of reconstructed MC events and $C$ at $L\to\infty$ is a constant equal to $2(2\pi)^{3}$, which is the total volume for the variables $\cos\theta$, $\Phi$, $\Psi$, and $\phi$. Since this constant does not depend on the free parameters, the knowledge of its value is unimportant for the maximum-likelihood method and it is set to unity below. In order to avoid computer calculations with very large values when building up the sum in Eqs. (59), the right-hand side is divided by $L$. The value of $\mathcal{K}_{m}(\mathcal{R})$ is evaluated separately in each of the $3\times 4$ ($Q^{2},-t^{\prime}$) cells. According to their experimental occurrence, the various settings of target and beam polarizations are assigned as $(P_{B},P_{T})=(+1,+1),$ $(+1,-1),$ $(-1,+1)$, and $(-1,-1)$ to the generated MC events. Thus four independent equations for $N_{i}$ can be obtained in each ($Q^{2},-t^{\prime}$) cell and the four functions $\mathcal{K}_{m}(\mathcal{R})$ can then be determined in each of the corresponding cells. The value of $\mathcal{K}_{m}(\mathcal{R})$ is also evaluated separately for electron and positron data. The finally maximized logarithm of the likelihood function is the sum of the logarithms of the likelihood function from Eq. (55) over all cells for both electron and positron data.

4.2 Background corrections

One of the main sources of background contamination to exclusive $\rho^{0}$-meson electroproduction in deep-inelastic scattering originates from SIDIS. A PYTHIA MC using GEANT3 geant to simulate the HERMES apparatus and tuned to the kinematics of the HERMES experiment th-pati , hereafter referred to as “SIDIS MC”, is used for the estimation of this background contribution. The same kinematic and geometrical requirements are imposed on both simulated and real data samples. The normalization of the MC data to the experimental data is performed in the region $2$ GeV $<\Delta E<20$ GeV (see Fig. 3 in Ref. DC-71 ), and the number of background events in the exclusive region is estimated. The fraction of SIDIS background $f_{bg}$ as estimated from the SIDIS MC is shown in the fifth column in Table 1.

It is assumed that the angular distribution of the SIDIS background events is reasonably well reproduced by the SIDIS MC simulation. The fit of the angular distribution of the SIDIS MC events under the exclusive peak for each ($Q^{2},-t^{\prime}$) cell is performed using Eq. (55) in which the substitutions $\mathcal{W}\to\mathcal{W}_{bg}$, $N_{i}\to N^{bg}_{i}$, and $\mathcal{R}\to\mathcal{S}$ must be performed, and the sum runs over all background MC events. The set of free parameters $\mathcal{S}$ represents the complete set of 15 “unpolarized” $u^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$ SDMEs describing the background. The normalization factor for the background, $N_{i}^{bg}$, is determined in an analogous way as done for the signal events, but in the present analysis the background angular distribution is considered to be independent of the beam and target polarizations, i.e., $P_{B}$ and $P_{T}$ are set to zero. After the fit of the SIDIS MC events, the angular distribution of the background is considered to be fixed, hence $\mathcal{W}_{bg}(\mathcal{S},\Phi_{i},\Psi_{i},\theta_{i},\phi_{i})$, $N^{bg}_{i}$, and $f_{bg}$ do not contain any free parameters.

The total probability $\text{d}w_{tot}$ to measure final-state particles from the reaction in Eqs. (2) and (3) or from SIDIS background in the small angular region $\text{d}\Omega$ is given by

[TABLE]

where $\text{d}w$ is given by Eq. (54) and $\text{d}w_{bg}$ is its analogue for the background process. In a similar manner as done in Eq. (55), the logarithm of the likelihood function $\mathcal{L}_{tot}$, which takes into account the background events, is given by

[TABLE]

The logarithms of the likelihood function are again calculated separately in each $(Q^{2},-t^{\prime})$ cell for both electron and positron data, and the finally maximized logarithm of the likelihood function is the sum of the logarithms of the likelihood function from Eq. (61) over all cells for both electron and positron data. As a result, the amplitude ratios in the entire kinematic region are obtained.

4.3 Choice of free parameters

As explained in Section 1, the angular distribution of the detected particles depends on the amplitude ratios. The total number of linearly independent helicity-amplitude ratios defined by Eqs. (25-26) is 17, which means that 34 real functions of $Q^{2}$, $-t^{\prime}$, and $W$ determine all SDMEs and angular distributions. As established in Ref. DC-84 , the large amplitudes at $-t^{\prime}\leq 0.4$ GeV2 and $Q^{2}\geq 1$ GeV2 are $T^{(1)}_{00}$, $T^{(1)}_{11}$, $U^{(1)}_{11}$, and $T^{(1)}_{01}$. For the ratios of large amplitudes, the parameterization of the $Q^{2}$ and $-t^{\prime}$ dependences is chosen as in Ref. DC-84 . For the small amplitude ratios, only the $-t^{\prime}$ dependence following from angular momentum conservation (see for instance Ref. Diehl ) is taken into account, while averaging over the kinematic range in $Q^{2}$. If all other amplitudes are expected to be significantly smaller than the large amplitudes, a possibility to extract the small amplitudes exists only if they are multiplied by large amplitudes. This means that they contribute linearly to the angular distribution.

The easiest way to interpret the extractability of the various helicity-amplitude ratios is through their contribution to the SDMEs, as detailed in Ref. Diehl . For the transversely polarized target, the SDMEs $n^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$ and $s^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$ contribute Diehl , while the contribution of the SDMEs $l^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$ is neglected in this analysis, since the latter are multiplied by the longitudinal component of the target polarization $|\hat{P}_{Z}|$. Indeed, $|\hat{P}_{Z}|$ is proportional to $\sin\theta_{\gamma}$ (see Eq. (42)), which in turn is proportional to $Q/\nu$ (according to Eq. (39)). At HERMES kinematics $Q$ is much smaller than $\nu$, with $Q/\nu$ of the order of 0.1. The helicity-amplitude ratios $t^{(2)}_{\lambda_{V}\lambda_{\gamma}}$ and $u^{(2)}_{\lambda_{V}\lambda_{\gamma}}$ can be extracted from data collected with a transversely polarized target, as they contribute linearly to the SDMEs $n^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$ and $s^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$, respectively. Contributions of squares of moduli of the helicity-amplitude ratios $u^{(1)}_{01}$, $u^{(1)}_{10}$, $u^{(1)}_{1-1}$ to $u^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$ are much smaller than the contribution of $|u^{(1)}_{11}|^{2}$ according to the hierarchy of amplitudes established in Refs. DC-24 and DC-84 . The small helicity amplitude ratios $u^{(1)}_{\lambda_{V}\lambda_{\gamma}}$ with $\lambda_{V}\neq\lambda_{\gamma}$ contribute linearly to the SDMEs $l^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$. The latter are multiplied by the small factor $\hat{P}_{Z}\sqrt{1-\epsilon}$ and cannot be extracted from the angular distributions of final-state pions. Therefore the helicity-amplitude ratios $u^{(1)}_{\lambda_{V}\lambda_{\gamma}}$, with the exception of $u^{(1)}_{11}$, are set equal to zero in the fit.

For an unpolarized target, only the SDMEs $u^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$ contribute to the angular distribution. As follows from the previous analysis at HERMES DC-84 the ratios $t^{(1)}_{\lambda_{V}\lambda_{\gamma}}$ and $|u^{(1)}_{11}|^{2}+|u^{(2)}_{11}|^{2}$ can be reliably extracted from data collected with an unpolarized target. The value of $|u^{(1)}_{11}|$ can be extracted from the unpolarized data, since the numerators of some SDMEs $u^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$ contain $|u^{(1)}_{11}|^{2}+|u^{(2)}_{11}|^{2}$. However, the phase $\delta_{u}$ of $u^{(1)}_{11}$ cannot be obtained reliably given the limited statistics in the present analysis. Another function that cannot be reliably extracted from the present data is ${\rm{Im}}\{t^{(1)}_{11}\}$. The reason is that it contributes mainly to the imaginary parts of the SDMEs $u^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$, which are multiplied by the small factor $P_{B}\sqrt{1-\epsilon}$. This factor is smaller than $0.15$, since $\epsilon$ is about $0.8$ and the mean value of $|P_{B}|$ is about $0.3$.

If ${\rm{Im}}\{t^{(1)}_{11}\}$ and the phase of $u^{(1)}_{11}$ are considered as free parameters, the fit to the angular distribution becomes unstable, leading to several local minima in the fit. In order to avoid such instabilities, the function ${\rm{Im}}\{t^{(1)}_{11}\}$ is taken from the previous analysis at HERMES DC-84 , where the number of events and the value of $|P_{B}|$ were larger than those in the present analysis, the fit was stable, and the minimum was unique. In order to take into account the $Q^{2}$ dependence of ${\rm{Im}}\{t^{(1)}_{11}\}$, the parametrization

[TABLE]

with $b$ taken from Ref. DC-84 is used. The phase $\delta_{u}$ is fixed from the data collected with a longitudinally polarized hydrogen target HERL1 ; HERL2 . A detailed discussion of this problem is given in Appendix A.

It is shown in Appendix B that

[TABLE]

does not contribute linearly to the angular distribution and hence is set to zero.

In the fit, 25 parameters $b_{i}$ (see Table 2) are extracted, which determine the following 25 real functions:

[TABLE]

Note that ${\rm{Im}}\{\frac{{t}^{(2)}_{11}}{{t}^{(1)}_{11}}\}$ is used rather than ${\rm{Im}}\{t^{(2)}_{11}\}$, since the latter is not independent of $\kappa$ and the inclusion of $\kappa$ in the fit leads to a divergence of the fit. Table 2 also shows the resulting parameters with their uncertainties. The correlations between the 25 parameters are listed in Table 3 in Appendix C.

4.4 Systematic uncertainties

In this subsection, the sources of systematic uncertainties and their effect on the extracted amplitude ratios are discussed. All systematic uncertainties except the one due to the uncertainty on the target and beam polarization measurements are added in quadrature to calculate the total systematic uncertainty. The statistical uncertainty and the total systematic uncertainty are added in quadrature to form the total uncertainty.

4.4.1 Systematic uncertainties due to beam and target polarization uncertainties

The measured mean value of the target polarization is $\langle|P_{T}|\rangle=0.72\pm 0.06$ Ami-th ; Jeron-th , i.e., the fractional uncertainty of the target polarization amounts to $0.08$. The ratios $t^{(2)}_{\lambda_{V}\lambda_{\gamma}}$ and $u^{(2)}_{\lambda_{V}\lambda_{\gamma}}$ have a corresponding scale uncertainty of 8%, since through their linear contribution to the “transverse” SDMEs $n^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$ and $s^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$, they are multiplied by $\langle|P_{T}|\rangle$. It was checked that the amplitude ratios $t^{(1)}_{11}$, $t^{(1)}_{10}$, $t^{(1)}_{1-1}$, $t^{(1)}_{01}$, and $|u^{(1)}_{11}|$, which can be extracted from data taken with an unpolarized target (see Ref. DC-84 ), are effectively insensitive to the uncertainty on the target polarization.

The fractional uncertainty on the beam polarization amounts to $2\%$ Sob . This results in an additional scale uncertainty on ${\rm{Im}}\{u^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}\}$, ${\rm{Re}}\{n^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}\}$, and ${\rm{Re}}\{s^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}\}$ of $2\%$, since these SDMEs enter the expression of the angular distribution of final-state particles multiplied by the beam polarization Diehl . From the expression of SDMEs in terms of helicity-amplitude ratios, it follows that there is an additional scale uncertainty of $2\%$ for ${\rm{Im}}\{t^{(1)}_{\lambda_{V}\lambda_{\gamma}}\}$, ${\rm{Re}}\{t^{(2)}_{\lambda_{V}\lambda_{\gamma}}\}$, and ${\rm{Re}}\{u^{(2)}_{\lambda_{V}\lambda_{\gamma}}\}$, while the influence of the uncertainty on the beam polarization can be neglected for ${\rm{Re}}\{t^{(1)}_{\lambda_{V}\lambda_{\gamma}}\}$, ${\rm{Im}}\{t^{(2)}_{\lambda_{V}\lambda_{\gamma}}\}$, and ${\rm{Im}}\{u^{(2)}_{\lambda_{V}\lambda_{\gamma}}\}$. The scale uncertainty arising from the uncertainty on the beam and target polarizations is not shown in the figures but quoted separately.

4.4.2 Systematic uncertainty due to the extraction method

In order to estimate the uncertainty due to the extraction method, yet another MC data sample was produced using a uniformly distributed angular distribution for exclusive $\rho^{0}$ production, which in contrast to the unweighted one used above for the normalization procedure, is weighted in order to mimic experimental data. It is obtained with the accept/reject method based on the experimental angular distribution and making use of the relevant parameters extracted from data. Using values of $\pm 0.30$ for the beam polarization and $\pm 0.72$ for the target polarization, it is analyzed in a way similar to the experimental data. The MC sample is divided into 20 independent sets such that each set contains the same number of exclusive events as the experimental data. In order to evaluate the systematic uncertainty, the difference between the output value of the $j$-th amplitude ratio in the $k$-th set $r_{j,k}^{(out)}$ and the input value of the same amplitude ratio $r_{j}^{(in)}$ is compared with the statistical uncertainty $\delta r_{j,k}^{(out)}$ of the output amplitude ratio in the $k$-th set. Averaging over the twenty sets, the relation

[TABLE]

with $K_{0}=20$ is used to calculate the systematic uncertainty due to the extraction method. If the sum in Eq. (65) is positive, then the obtained value of $\Delta r^{meth}_{j}$ is set as systematic uncertainty; otherwise the systematic uncertainty is set to zero.

4.4.3 Systematic uncertainty due to the background contribution

The helicity-amplitude ratios are extracted from the experimental data once taking into account the background contribution (see Eq. (61)) and once neglecting this contribution (see Eq. (55)). The systematic uncertainty from the background contribution of each amplitude ratio is computed as the modulus of the difference of the amplitude ratios obtained for these two cases. This conservative approach is used, since the background correction is estimated from MC data instead of experimental data. The SDMEs $n^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$ and $s^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$ are, as shown in Ref. DC-71 , much less sensitive to the background contribution than the SDME $u^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$, since they enter the formula for the angular distribution multiplied by the target polarization. As the amplitude ratios $t^{(2)}_{\lambda^{\prime}_{N}\lambda_{N}}$ and $u^{(2)}_{\lambda^{\prime}_{N}\lambda_{N}}$ contribute linearly to these SDMEs, they are expected to be less sensitive to the background contribution than the amplitudes relevant for scattering off an unpolarized target. It was checked that this is indeed the case. The small influence of the background correction to the nucleon-helicity-flip amplitude ratios $t^{(2)}_{\lambda^{\prime}_{N}\lambda_{N}}$ and $u^{(2)}_{\lambda^{\prime}_{N}\lambda_{N}}$ can be explained by the statistical correlations between these amplitude ratios and $t^{(1)}_{\lambda^{\prime}_{N}\lambda_{N}}$, $|u^{(1)}_{11}|$.

4.4.4 Systematic uncertainty due to the omission of inaccessible amplitude ratios

Another source of systematic uncertainty originates from setting $u^{(1)}_{10},\;u^{(1)}_{01},\;u^{(1)}_{1-1}$ and $\kappa$, given by Eq. (63), equal to zero in the fit. In order to estimate this systematic uncertainty, the following procedure is applied. Since $|T^{(2)}_{00}/T^{(1)}_{00}|$ and $|T^{(2)}_{11}/T^{(1)}_{11}|$ are proportional to $\sqrt{-t^{\prime}}/(2M)$, the unmeasured parameter $\kappa$ is estimated as

[TABLE]

The calculations are performed for both signs in Eq. (66) and the corresponding systematic uncertainties are averaged in quadrature (see Eq. (70)).

For one-pion exchange, the ratios of the UPE amplitudes $U^{(1)}_{10}/U^{(1)}_{11}$, $U^{(1)}_{01}/U^{(1)}_{11}$, and $U^{(1)}_{1-1}/U^{(1)}_{11}$ are known. Supposing that these ratios have approximately the same value for the full amplitudes, one obtains

[TABLE]

For the calculation of the systematic uncertainty, the value $|u^{(1)}_{11}|$ is obtained from the fit in the present analysis, while the phase shift $\delta_{u}=-39.2$ degrees, which is taken from the results of the measurement of the longitudinal double-spin asymmetry in exclusive $\rho^{0}$-meson electroproduction from Refs. DC-24 ; HERL1 ; HERL2 , corresponds to the value $A_{1}^{\rho}=0.24$ (see appendix A). The systematic uncertainty due to the omission of all inaccessible amplitude ratios is estimated from the relation

[TABLE]

where the values of $r_{j}$ are those obtained in the 25-parameter fit when all the amplitudes not extracted are set to zero. The ratios $r_{j}^{(in+)}$ and $r_{j}^{(in-)}$ denote the extracted $j$-th amplitude ratio obtained in the fit in which the amplitude ratios $u^{(1)}_{10}$, $u^{(1)}_{01}$, and $u^{(1)}_{1-1}$ are calculated using Eqs. (67-69), while $\kappa$ is taken according to Eq. (66), once with positive sign and once with negative sign, respectively.

4.4.5 Systematic uncertainty due to the experimental uncertainty of $A_{1}^{\rho}$

The uncertainty on the asymmetry $A_{1}^{\rho}$, as explained in appendix A, leads to a range of $\delta_{u}$ between $-26.2$ degrees and $-51.1$ degrees, corresponding to $A_{1}^{\rho}=0.24+0.14$ and $A_{1}^{\rho}=0.24-0.14$.

In order to estimate the systematic uncertainty of the obtained amplitude ratios due to the experimental uncertainty of $A_{1}^{\rho}$, the 25-parameter fits with $\delta_{u}=-26.2$ degrees and $\delta_{u}=-51.1$ degrees are performed in addition to the fit with $\delta_{u}=-39.2$ degrees. Two sets of the amplitude ratios $r_{j}^{(+)}$ and $r_{j}^{(-)}$ are obtained, correspondingly. The systematic uncertainty is calculated as

[TABLE]

The systematic uncertainty due to the experimental uncertainty on $A_{1}^{\rho}$ is less than the statistical uncertainty. It was found that the most important ratios $t^{(1)}_{\lambda_{V}\lambda_{\gamma}}$ and $|u^{(1)}_{11}|$ are almost insensitive to the value of $A_{1}^{\rho}$.

4.4.6 Systematic uncertainty due to the experimental uncertainty of ${\rm{Im}}\{t^{(1)}_{11}\}$

The parameter $b$, which enters the parametrization of ${\rm{Im}}\{t^{(1)}_{11}\}$ as given in Eq. (62), is equal to $(0.340\pm 0.025)$ GeV*-1* DC-84 . Its uncertainty is another source of systematic uncertainty. The basic fit performed with $b=0.340$ GeV*-1* gives the amplitude ratios $r_{j}$, while the values $(b=0.340\pm 0.025)$ GeV*-1* correspond to the ratios $r_{j}^{(im+)}$ and $r_{j}^{(im-)}$. The systematic uncertainty of the helicity-amplitude ratio $r_{j}$ is calculated as

[TABLE]

5 Results on the amplitude ratios

5.1 Discussion of results on the amplitude ratios

The result obtained from the 25-parameter fit is presented in Fig. 2. The results for the large amplitudes are calculated at $-t^{\prime}=0.132$ GeV2 and $Q^{2}=1.93$ GeV2, while integrating $W$ over the entire kinematic region. The results for the small amplitudes are calculated at $-t^{\prime}=0.132$ GeV2, while integrating $Q^{2}$ and $W$ over the entire kinematic region. Here, the values $-t^{\prime}=0.132$ GeV2 and $Q^{2}=1.93$ GeV2 are the mean values of the kinematic variables over the entire kinematic region, $0.0$ GeV${}^{2}\leq-t^{\prime}\leq 0.40$ GeV2, $1.0$ GeV${}^{2}\leq Q^{2}\leq 7.0$ GeV2, and $3.0$ GeV $\leq W\leq 6.3$ GeV. The mean value of $W$ over the entire kinematic region is $4.73$ GeV. The NPE amplitude ratio without nucleon-helicity flip, $t_{11}^{(1)}$, is the dominant amplitude ratio. Its real and imaginary parts differ from zero by more than five standard deviations. As also already known from the previous analysis DC-84 , ${\rm{Re}}\{t_{01}^{(1)}\}$ is significantly non-zero. In this analysis, the UPE amplitude ratios without nucleon-helicity flip, ${\rm{Re}}\{u^{(1)}_{11}\}$ and ${\rm{Im}}\{u^{(1)}_{11}\}$, are individually extracted and found to be nonzero with a significance of about four standard deviations of the total uncertainty. The values of $|u^{(1)}_{11}|$ and $|u^{(2)}_{11}|$, with $\sqrt{|u^{(1)}_{11}|^{2}+|u^{(2)}_{11}|^{2}}=0.35\pm 0.06$, agree with the result $\sqrt{|u^{(1)}_{11}|^{2}+|u^{(2)}_{11}|^{2}}\approx 0.40\pm 0.02$ obtained in the previous HERMES analysis DC-84 . The extracted values of the amplitude ratios show that the main contribution to the term $\sqrt{|u^{(1)}_{11}|^{2}+|u^{(2)}_{11}|^{2}}$ comes from the amplitude $U^{(1)}_{11}$ without nucleon-helicity flip, and in particular they show that $|U^{(1)}_{11}|^{2}\gg|U^{(2)}_{11}|^{2}$. The amplitude ratios ${\rm{Im}}\{t^{(2)}_{01}\}$, ${\rm{Im}}\{u^{(2)}_{11}\}$, and ${\rm{Im}}\{u^{(2)}_{10}\}$ deviate from zero by about two standard deviations, while the other extracted amplitude ratios with nucleon-helicity flip are consistent with zero within two standard deviations. The amplitude ratios ${\rm{Im}}\{t^{(2)}_{01}\}$ and ${\rm{Im}}\{u^{(2)}_{10}\}$ are part of those ratios in Eq. (24), which can be nonzero at $-t^{\prime}=0$. Among the amplitude ratios that can be zero at $-t^{\prime}=0$, only the amplitude ratio ${\rm{Im}}\{u^{(2)}_{11}\}$, which is proportional to $\sqrt{-t^{\prime}}$ at $-t^{\prime}\to 0$, differs from zero by about two standard deviations of the total uncertainty.

5.2 Comparison of calculated SDMEs with directly extracted SDMEs

A comparison of the SDMEs obtained from the SDME method in Refs. DC-24 and DC-71 to those calculated from the amplitude ratios extracted in the present analysis is presented in Figs. 3-5. The SDMEs are calculated in each individual bin using the average kinematics in the parameterizations obtained for the amplitude ratios. Furthermore, their mean value is then determined by weighting the SDME value calculated in a given bin by the number of events in this bin. The correlation matrix for the 25 parameters is taken into account for the calculation of the statistical uncertainties of the SDMEs $u^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$, $n^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$, and $s^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$ obtained in the amplitude method. As already mentioned in Section 1, the SDMEs $n^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$ and $s^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma\ }}$, presented in Figs. 4 and 5, can only be extracted from measurements with a transversely polarized target so that the helicity-flip amplitude ratios $t^{(2)}_{\lambda_{V}\lambda_{\gamma}}$ and $u^{(2)}_{\lambda_{V}\lambda_{\gamma}}$ are extracted in this paper for the first time. The systematic uncertainties of the SDMEs from the amplitude method are determined in an analogous way as for the amplitude ratios by varying the relevant parameters, as explained in Sections 4.4.2 to 4.4.6, and recalculating the corresponding SDMEs. The total uncertainty is the sum in quadrature of the statistical and the total systematic uncertainties.

The SDMEs in Figs. 3-5 are reordered according to the SDME classes proposed in Refs. DC-24 ; Ami-th . In these figures also class-F SDMEs are shown. Although the double-helicity-flip contribution was a priori not fitted, non-zero values are obtained for this class of SDMEs because these SDMEs also receive contributions from other helicity transitions.

Those SDMEs that can be extracted only from data taken with a longitudinally polarized lepton beam are shown in shaded areas. Figure 3 shows that for each SDME $u^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$ determined from our present results, there exists an SDME $u^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$ published in Ref. DC-24 . However, Figs. 4 and 5 show that for some of the SDMEs $n^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$ and $s^{\lambda_{V}\lambda^{\prime}_{V}}_{\lambda_{\gamma}\lambda^{\prime}_{\gamma}}$ determined in this analysis no published results from Ref. DC-71 exist, because the beam polarization was not exploited in the analyses presented in Ref. DC-71 . While in Refs. DC-24 and DC-71 a total of 53 SDMEs could be extracted, the amplitude method presented here allows for the calculation of 71 SDMEs based on the extraction of 25 parameters.

As seen from the figures, there is reasonable agreement between SDMEs obtained with the SDME method and those from the amplitude method. It is possible that the values of the SDMEs obtained in these two methods do not coincide, because the parameter space for SDMEs in the SDME method is different from that in the amplitude method. Indeed, the SDMEs should belong to a special region in the 71-dimensional real space to give a non-negative angular distribution. However, at present the equations determining the boundaries of this region are unknown. The physical SDMEs can be represented in terms of 17 helicity-amplitude ratios. This restricts the region in the 71-dimensional space. This requirement is not taken into account in the SDME method, but it suppresses statistical fluctuations especially when a SDME value is close to the boundary of the allowed region. Note that the positivity requirement on the angular distribution is inherent to the amplitude method, while it is not to the SDME method, where it is usually imposed artificially.

5.3 Comparison to amplitudes calculated in a GPD-based handbag model

Within the handbag approach (see, e.g., Refs. golos3 ; golos5 ), the amplitudes for $\gamma^{*}_{L}\to V_{L}$ and $\gamma^{*}_{T}\to V_{T}$ transitions are given by convolutions of appropriate subprocess amplitudes, ${\cal H}_{\lambda_{V}\lambda^{\prime}_{q}\lambda_{\gamma}\lambda_{q}}$, with the GPDs $H^{a}$, $E^{a}$, $\tilde{H}^{a}$, $\tilde{E}^{a}$,

[TABLE]

for the case of proton helicity non-flip and

[TABLE]

for the case of proton helicity flip.

Here, $\lambda_{q}$ and $\lambda^{\prime}_{q}$ are the helicities of the emitted and reabsorbed quarks from the proton, respectively, with $\lambda_{q}=\lambda^{\prime}_{q}$ for quark helicity non-flip GPDs, $\mu=\lambda_{V}=\lambda_{\gamma}$, $e_{a}$ are the quark charges in units of the positron charge $e_{0}$, and $x$ is an internal integration variable. The skewness $\xi$ is related to $x_{B}$ by $\xi=x_{B}/(2-x_{B})$ up to corrections of order $1/Q^{2}$, where $x_{B}$ is the Bjorken scaling variable defined as

[TABLE]

The coefficients ${\cal C}_{\rho}^{a}$ are appropriate flavor factors (${\cal C}_{\rho^{0}}^{u}=-{\cal C}_{\rho^{0}}^{d}=1/\sqrt{2})$. Because of parity invariance, the first of the each two terms in square brackets in Eqs. (LABEL:eq:non-flip) and (74) of both the quark and gluon parts behave like natural-parity-exchanges (see Eq. (13)), while the second terms are of the unnatural-parity type (see Eq. (14)). For the $\gamma^{*}_{L}\to V_{L}$ transition, there is a rigorous proof of factorization in hard subprocesses and GPD amplitudes CFS in the generalized Bjorken regime of large $Q^{2}$, large $W$ but fixed $x_{B}$. Contributions to longitudinal amplitudes come from GPDs $H$ and $E$, only. In contrast, the $\gamma^{*}_{T}\to V_{T}$ amplitudes are infrared singular in collinear approximation. In order to regularize this singularity, the so-called modified perturbative approach has been used in Refs. golos2 ; golos3 in which quark transverse momenta are retained in the subprocess, while the emission and reabsorption of the partons from the proton are still treated collinear to the incoming and outgoing proton momenta. The quark transverse momenta in the subprocess imply a separation of color charges, which results in gluon radiation, as it was calculated in Ref. BS in next-to-leading-log approximation and resummed to all orders of perturbative QCD. This gluon radiation is also taken into account in Refs. golos2 ; golos3 .

Measurements of the spin asymmetry with a transversely polarized target in exclusive $\pi^{+}$ leptoproduction Hermes-pi+ revealed that the transversity or helicity-flip GPDs play an important role in the $\gamma^{*}_{T}\to V_{L}$ transitions golos5 . The corresponding amplitudes read

[TABLE]

According to the discussion presented in Section 2.3, the amplitude $F_{0\frac{1}{2}\pm 1\frac{1}{2}}$ is of natural-parity type, while the proton helicity-flip amplitude $F_{0-\frac{1}{2}1\frac{1}{2}}$ has no specific parity. In the subprocess amplitude ${\cal H}_{0-\frac{1}{2}1\frac{1}{2}}$, quark and antiquark forming the longitudinally polarized $\rho^{0}$ meson have the same helicity. This fact necessitates the use of a twist-3 meson wave function.

An important role is played by the pion-pole contribution. As it was discussed in Ref. golos4 , the pion pole is to be treated as a one-pion-exchange term, since its evaluation through the GPD $\widetilde{E}$ underestimates the contribution grossly. An important element of the pion-pole contribution is the $\pi-\rho$ transition form factor $g_{\pi\rho}(Q^{2})$. It is estimated to be a third of the $\pi-\omega$ form factor golos4 that was extracted from the HERMES measurement of the $\omega$ SDMEs DC-95 . The factor 1/3 arises from the different quark content of the $\omega$ and $\rho^{0}$ mesons. The pion pole represents an unnatural-parity contribution and, as can be shown, it contributes dominantly to the amplitudes for transversely polarized vector mesons:

[TABLE]

with $\mu=0,1$ and $m_{\pi}$ the mass of the pion. The explicit expressions for $F^{pole}$ are given in Ref. golos4 . They are to be added to the amplitudes from Eqs. (LABEL:eq:non-flip) and  (74). The $\gamma_{L}^{*}\rightarrow V_{T}$ amplitudes receive contributions from only the pion pole. With regard to the large $Q^{2}$-behavior of the $\pi-V$ transition form factor, these transition amplitudes are suppressed by $1/Q$ and $1/Q^{2}$, respectively, as compared to the asymptotically leading $\gamma^{*}_{L}\to V_{L}$ amplitudes. The pion-pole contributions to the $\gamma^{*}_{T}\to V_{-T}$ and $\gamma^{*}_{T}\to V_{L}$ transition amplitudes are suppressed more strongly and are therefore neglected. As already discussed in Section 2.3, there is no UPE contribution for the $\gamma^{*}_{L}\to V_{L}$ amplitudes.

Details of the calculations of the amplitudes as well as the parametrization of the GPDs, the meson wave-function and the $\pi$-$\rho$ transition form factor can be found in the original papers golos3 ; golos4 . The evaluation of the amplitudes represent an intermediate step of the calculation of the observables discussed in these papers. These amplitudes are divided by $F_{0\frac{1}{2}0\frac{1}{2}}=T^{(1)}_{00}$ in order to obtain the amplitude ratios that can be compared to the ones discussed above. The phase convention from Eq. (9) is taken into account.

Figure 6 shows the comparison of all amplitude ratios determined by fitting the HERMES data to those calculated using the GK model. The HERMES data are obtained as explained in Section 5.1, while the values obtained in the GK model are calculated at $W=4.73$ GeV, $Q^{2}=1.93$ GeV2, and $-t^{\prime}=0.132$ GeV2. As shown by the following detailed comparison, good overall agreement is found.

• $t_{11}^{(1)}$:

Contributions come from GPDs $H-\xi^{2}/(1-\xi^{2})E\simeq H$. Good agreement is observed for the real part, which is by far the largest amplitude ratio. The calculated imaginary part appears to be too small. Note that a part of this difference is due to the known underestimation of the relative phase between the $\gamma^{*}_{T}\to V_{T}$ and $\gamma^{*}_{L}\to V_{L}$ amplitudes in the GK model DC-24 .

• $u_{11}^{(1)}$:

Contributions come from GPDs $\widetilde{H}$ and the pion pole. The GK calculations underestimate the unnatural-parity contribution to the $\gamma_{T}^{*}\to V_{T}$ amplitude, which is related to the small unnatural-parity cross section used in the GK model golos4 . It may be traced back to the neglect of the non-pole contribution of the GPD $\widetilde{E}$ or to a too small value for the $\pi-\rho$ transition form factor in the GK model.

• $t_{11}^{(2)}$:

Contributions come from GPDs $E$. The calculated imaginary part agrees with the measurement.

• $u_{11}^{(2)}$:

In GK calculations only the pion pole contributes since $\widetilde{E}$ is neglected, so that the GK result is mirror symmetric upon sign change of the $\pi-\rho$ transition form factor. Good agreement with the data is seen for the positive sign.

• $t_{00}^{(2)}$:

Contributions come from GPDs E. Agreement is observed with the measurement.

• $t_{01}^{(1)}$:

Contributions come from GPDs $\bar{E}_{T}$. Agreement is observed with the measurement.

• $t_{01}^{(2)}$:

Contributions from GPDs $H_{T}$. There is no pion-pole contribution to this ratio, hence data cannot decide on the sign of the form factor. The measured imaginary part seems to be lower than the GK calculation.

• $u_{01}^{(2)}$:

Contributions from GPDs $H_{T}$. Since these GPDs have no specific parity, $u_{01}^{(2)}$ is equal to $-t_{01}^{(2)}$ in the GK calculation.

• $u_{10}^{(2)}$:

Contributions come from the pion pole only, so that the GK result is mirror symmetric upon sign change of the $\pi-\rho$ transition form factor. The positive sign is favored by the data.

The $\gamma^{*}_{T}\to V_{-T}$ amplitudes, corresponding to the amplitude ratios $t_{1-1}^{(1)}$, $t_{1-1}^{(2)}$, and $u_{1-1}^{(2)}$, are neglected in the GK model. This is seen to be in reasonable agreement with the data. Only gluon transversity GPDs could contribute and the contribution from the pion pole is suppressed by $1/Q^{3}$ as compared to the longitudinal amplitudes. Both are neglected in the GK model.

As discussed in Section 4.3, the ratios $u_{01}^{(1)}$, $u_{10}^{(1)}$ and $u_{1-1}^{(1)}$ cannot be determined experimentally in the present analysis and are hence put equal to zero. In the GK model, $u_{01}^{(1)}$ and $u_{1-1}^{(1)}$ are also set equal to zero, while $u_{10}^{(1)}$ is non-zero due to a contribution from the pion pole, but small. Apart from the $\gamma^{*}_{T}\to V_{-T}$ amplitudes, $u_{01}^{(1)}$ and $u_{1-1}^{(1)}$, also $t_{10}^{(1)}$ and $t_{10}^{(2)}$ are set equal to zero. This is consistent with what is extracted from the data.

As the unnatural-parity amplitudes depend on the sign of the $\pi-\rho$ transition form factor, a conclusion on the sign of the latter can be drawn when comparing the calculated GK amplitude ratios to the data. Only the amplitude ratios $u_{11}^{(2)}$ and $u_{10}^{(2)}$ appear sensitive to the sign of the form factor and are hence used to calculate the $\chi^{2}$ per degree of freedom, i.e., $ndf=4$. For the positive sign $\chi^{2}/ndf=1.8/4$ is obtained and for the negative sign $\chi^{2}/ndf=30.3/4$. Hence the positive sign of this form factor is clearly favored.

6 Summary and conclusions

Exclusive electroproduction of $\rho^{0}$ mesons is studied in the HERMES experiment, using data collected with a $27.6$ GeV longitudinally polarized electron/positron beam and a transversely polarized hydrogen target in the kinematic region $1.0$ GeV${}^{2}<Q^{2}<7.0$ GeV2, $3.0$ GeV $<W<6.3$ GeV, and $-t^{\prime}<0.4$ GeV2. The fit to these data using an unbinned maximum-likelihood method with 25 free parameters permits the extraction of ratios of natural-parity-exchange amplitudes $T_{\lambda_{V}\lambda^{\prime}_{N}\lambda_{\gamma}\lambda_{N}}$ without nucleon-helicity flip ($\lambda^{\prime}_{N}=\lambda_{N}$) and, for the first time, both the natural-parity-exchange and unnatural-parity-exchange

amplitudes ($T_{\lambda_{V}\lambda^{\prime}_{N}\lambda_{\gamma}\lambda_{N}}$ and $U_{\lambda_{V}\lambda^{\prime}_{N}\lambda_{\gamma}\lambda_{N}}$) with nucleon-helicity flip ($\lambda^{\prime}_{N}\neq\lambda_{N}$), all obtained relative to the amplitude $T_{0\frac{1}{2}0\frac{1}{2}}$, which is the largest amplitude in the kinematic region of $Q^{2}>1$ GeV here considered. In particular, the modulus of the amplitude ratio $U_{1\frac{1}{2}1\frac{1}{2}}/T_{0\frac{1}{2}0\frac{1}{2}}$, the real part of $T_{1\frac{1}{2}1\frac{1}{2}}/T_{0\frac{1}{2}0\frac{1}{2}}$ as well as the real and imaginary parts of the amplitude ratios $T_{1\frac{1}{2}0\frac{1}{2}}/T_{0\frac{1}{2}0\frac{1}{2}}$, $T_{1\frac{1}{2}-1\frac{1}{2}}/T_{0\frac{1}{2}0\frac{1}{2}}$ and $T_{0\frac{1}{2}1\frac{1}{2}}/T_{0\frac{1}{2}0\frac{1}{2}}$ are extracted. They were also obtained in the previous HERMES analysis DC-84 and the amplitude ratios are in a good agreement with one another. The values of ${\rm{Im}}\{T_{1\frac{1}{2}1\frac{1}{2}}/T_{0\frac{1}{2}0\frac{1}{2}}\}$ and the phase of the ratio $U_{1\frac{1}{2}1\frac{1}{2}}/T_{0\frac{1}{2}0\frac{1}{2}}$ are taken from the HERMES results DC-24 ; DC-84 ; HERL1 ; HERL2 . By performing the fit, the ratios of small nucleon-helicity-flip natural-parity-exchange amplitudes

[TABLE]

and unnatural-parity-exchange amplitudes

[TABLE]

to $T_{0\frac{1}{2}0\frac{1}{2}}$ as well as

[TABLE]

are obtained for the first time, as the data presented here were taken using a transversely polarized hydrogen target and a longitudinally polarized lepton beam.

Within the total experimental uncertainty, all determined amplitude ratios with nucleon-helicity flip are consistent with zero. The extracted values of the amplitude ratios show that the main contribution to the quantity $\frac{\sqrt{|U_{1\frac{1}{2}1\frac{1}{2}}|^{2}+|U_{1\frac{1}{2}1-\frac{1}{2}}|^{2}}}{|T_{0\frac{1}{2}0\frac{1}{2}}|}$ obtained in Ref. DC-84 originates from the unnatural-parity-exchange amplitude $U_{1\frac{1}{2}1\frac{1}{2}}$ and that $|U_{1\frac{1}{2}1\frac{1}{2}}|^{2}\gg|U_{1\frac{1}{2}1-\frac{1}{2}}|^{2}$. Furthermore, it is shown that the 53 SDMEs extracted in Refs. DC-24 ; DC-71 can be described with good accuracy using the 25 amplitude ratios obtained in the present analysis. By also exploiting the longitudinal beam polarization, 18 additional $\rho^{0}$ SDMEs are determined from the extracted 25 parameters for the first time.

The unnatural-parity amplitudes depend on the sign of the $\pi-\rho$ transition form factor, so that the comparison of certain amplitude ratios to calculations within a GPD-based handbag model taking into account the contribution from pion exchange allows the conclusion that the positive sign of this form factor is favored.

Together with precise data on the unpolarized differential cross section $\text{d}\sigma/\text{d}t$ of exclusive $\rho^{0}$ production in deep-inelastic scattering, the extracted amplitude ratios could be used to obtain the amplitude $T_{0\frac{1}{2}0\frac{1}{2}}$, for which the factorization property is proven.

Acknowledgements.

**Acknowledgments ** We gratefully acknowledge the DESY management for its support and the staff at DESY aas well as the collaborating institutions for their significant effort. This work was supported by the State Committee of Science of the Republic of Armenia, Grant No. 15T-1C401; the FWO-Flanders and IWT, Belgium; the Natural Sciences and Engineering Research Council of Canada; the Alexander von Humboldt Stiftung, the German Bundesministerium für Bildung und Forschung (BMBF), and the Deutsche Forschungsgemeinschaft (DFG); the Italian Istituto Nazionale di Fisica Nucleare (INFN); the MEXT, JSPS, and G-COE of Japan; the Dutch Foundation for Fundamenteel Onderzoek der Materie (FOM); the Russian Academy of Science and the Russian Federal Agency for Science and Innovations; the Basque Foundation for Science (IKERBASQUE) and MINECO (Juan de la Cierva), Spain; the U.K. Engineering and Physical Sciences Research Council, the Science and Technology Facilities Council, and the Scottish Universities Physics Alliance; as well as the U.S. Department of Energy (DOE) and the National Science Foundation (NSF).

Appendix A The longitudinal double-spin asymmetry

The phase $\delta_{u}$ can be determined using the HERMES data HERL1 ; HERL2 taken with a longitudinally polarized hydrogen target with better accuracy than using measurements with a transversely polarized hydrogen target. The longitudinal double-spin asymmetry in exclusive $\rho^{0}$-meson electroproduction is defined as HERL1

[TABLE]

Here, $\text{d}\sigma_{\frac{1}{2}}/\text{d}t$ and $\text{d}\sigma_{\frac{3}{2}}/\text{d}t$ denote the differential cross section for $\rho^{0}$-meson production with a transverse virtual photon, where $1/2$ and $3/2$ are the total projections of the spins of $\gamma^{*}$ and $p$ onto the photon momentum in the $\gamma^{*}p$ CM system, respectively. For the transformations in Eq. (78), Eqs. (12), (15-16) and (25-26) are used. Equation (78) can be rewritten in terms of the phase $\delta_{u}$ and the phase $\delta$ of the amplitude ratio $t^{(1)}_{11}$ (which is nothing else than the phase difference between the amplitudes $T^{(1)}_{11}$ and $T^{(1)}_{00}$) as

[TABLE]

Since the moduli of the amplitude ratios $|t^{(1)}_{11}|$ and $|u^{(1)}_{11}|$ and the phase $\delta$ were extracted from HERMES data DC-24 ; DC-84 , the measured value of the asymmetry $A_{1}^{\rho}=0.24$ HERL1 ; HERL2 gives a value of $\delta_{u}-\delta\approx\pm 69.8$ degrees. The sign of the latter is obtained from

[TABLE]

The UPE amplitude $U^{(1)}_{11}$ is mainly the amplitude of pion exchange, which is a real positive function at small $-t^{\prime}$. At high energies, the imaginary part of the NPE amplitude $T^{(1)}_{11}$ is positive and much larger than its real part, so that the ratio $U^{(1)}_{11}/T^{(1)}_{11}$ has to have a negative imaginary part leading to $\sin(\delta_{u}-\delta)<0$. This gives $\delta_{u}-\delta\approx-69.8$ degrees and hence $\delta_{u}=-39.2$ degrees if the value $\delta=30.6$ degrees is taken from Ref. DC-24 . The range of $\delta_{u}$ considering the values 0.38 and 0.1 for the asymmetry $A_{1}^{\rho}$ (corresponding to $\pm$ one standard deviation) is $\delta_{u}=-26.2$ and $\delta_{u}=-51.1$ degrees.

Appendix B Linear contribution of amplitudes $T^{(2)}_{00}$ and $T^{(2)}_{11}$ to the angular distribution

Let us use the following parameterization for the ratios of the small amplitudes $T^{(2)}_{00}$ and $T^{(2)}_{11}$ to the big amplitudes $T^{(1)}_{00}$ and $T^{(1)}_{11}$:

[TABLE]

Hence the terms in the angular distribution that contain linear contributions from the small amplitudes $T^{(2)}_{00}$ and $T^{(2)}_{11}$ can be written as Diehl

[TABLE]

where the normalization factor $\mathcal{N}$ is defined by Eqs. (45-47). Here, the azimuthal angle $\phi_{s}$ between the transverse component of the target polarization with respect to the virtual-photon momentum and the lepton scattering plane is defined as in Ref. Diehl . It is related to the angle $\Psi$ through the following equations:

[TABLE]

where the angle $\theta_{\gamma}$ is defined by Eqs. (38) and (39). Note that the angle $\Phi$ used in Eq. (83) is related to $\phi_{[6]}$ used in Ref. Diehl by the equation $\Phi=2\pi-\phi_{[6]}$. Also, the angle $\phi$ of the present paper is denoted in Ref. Diehl by ${\varphi}$.

As seen from Eq. (83) there is a contribution from the combination $(\beta_{11}-\beta_{00})$ of the parameters $\beta_{11}$ and $\beta_{00}$, but there is no contribution from $(\beta_{11}+\beta_{00})$ to the angular distribution. Therefore, the three parameters $\alpha_{11}$, $\alpha_{00}$, and $(\beta_{11}-\beta_{00})/2$ can be extracted from the angular distribution, while

[TABLE]

cannot be reliably obtained from data on a transversely polarized target. Note that the second-order contributions proportional to $|\beta_{11}|^{2}$ and $|\beta_{00}|^{2}$ exist, since at least the normalization factor $\mathcal{N}$ in formula (44) defined by Eqs. (45-47) contains squared moduli of all amplitudes. However, these second-order contributions are negligibly small.

Appendix C Correlation matrix for fit parameters

Table 3 lists the correlations between the 25 parameters of the fit of helicity-amplitude ratios defined in Table 2.