Spin density matrix elements in exclusive ρ0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho ^0$$\end{document} meson muoproduction

We report on a measurement of Spin Density Matrix Elements (SDMEs) in hard exclusive ρ0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho ^0$$\end{document} meson muoproduction at COMPASS using 160 GeV/c polarised μ+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mu ^{+}$$\end{document} and μ-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mu ^{-}$$\end{document} beams impinging on a liquid hydrogen target. The measurement covers the kinematic range 5.0 GeV/c2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c^2$$\end{document}<W<\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$< W<$$\end{document} 17.0 GeV/c2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c^2$$\end{document}, 1.0 (GeV/c)2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^2$$\end{document}<Q2<\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$< Q^2<$$\end{document} 10.0 (GeV/c)2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^2$$\end{document} and 0.01 (GeV/c)2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^2$$\end{document}<pT2<\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$< p_{\textrm{T}}^2<$$\end{document} 0.5 (GeV/c)2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^2$$\end{document}. Here, W denotes the mass of the final hadronic system, Q2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q^2$$\end{document} the virtuality of the exchanged photon, and pT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_{\textrm{T}}$$\end{document} the transverse momentum of the ρ0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho ^0$$\end{document} meson with respect to the virtual-photon direction. The measured non-zero SDMEs for the transitions of transversely polarised virtual photons to longitudinally polarised vector mesons (γT∗→VL\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma ^*_T \rightarrow V^{ }_L$$\end{document}) indicate a violation of s-channel helicity conservation. Additionally, we observe a dominant contribution of natural-parity-exchange transitions and a very small contribution of unnatural-parity-exchange transitions, which is compatible with zero within experimental uncertainties. The results provide important input for modelling Generalised Parton Distributions (GPDs). In particular, they may allow one to evaluate in a model-dependent way the role of parton helicity-flip GPDs in exclusive ρ0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho ^0$$\end{document} production.


Introduction
Exclusive vector meson production in lepton-nucleon scattering provides a convenient tool for studying the production mechanism and, in a model-dependent way, the structure of the nucleon.In this paper, exclusive  0 meson muoproduction on the proton is studied: In the one-photon-exchange approximation, this process is described by the interaction of a virtual photon  * with the target proton : At high virtuality  2 of the photon, this process is known as Hard Exclusive Meson Production (HEMP).A wealth of information is contained in the spin density matrix elements (SDMEs), which are the observables describing how the spin components of the virtual photon are transferred to those of the created vector meson [1,2].The comparison of the new  0 results presented in this paper to our previous  results [3] will provide insight into details of their respective production mechanism, because  0 and  vector mesons have different quantum numbers and hence different quark-flavour and gluon contributions to the cross section.
The colour dipole model describes HEMP as a fluctuation of the virtual photon into a quark-antiquark ( q) pair that scatters off the nucleon and then hadronises into the final vector meson [4].Regge phenomenology and perturbative QCD (pQCD) provide complementary approaches to describe the scattering of the  q pair off the nucleon.The interaction of the  q pair with the nucleon depends on the transverse separation between  and q.A pair with large transverse separation is thought to interact primarily softly, which is described in Regge phenomenology [5] by the exchange of a pomeron or a secondary reggeon.The interaction of a  q pair with small transverse separation is calculable in pQCD.In lowest order of the strong coupling constant   , this hard interaction is mediated by the exchange of a gluon-gluon or quark-antiquark system.In this approach, it is possible to calculate not only transitions without spin-flip induced by both longitudinally () and transversely () polarised virtual photons,  *  →  0  and  *  →  0  , but also to estimate single and double spin-flip transitions.In an alternative approach, the framework of General Parton Distributions (GPDs) [6][7][8][9][10] can be used to describe HEMP.These distribution functions contain a wealth of new information on the parton structure of the nucleon.For HEMP by longitudinally polarised virtual photons, the amplitude was proven to factorise into a hard-scattering part and a soft part [9,11].While the former is calculable in pQCD, the latter contains GPDs that describe the structure of the probed nucleon and a distribution amplitude that accounts for the structure of the produced meson.This factorisation is referred to as collinear because parton transverse momenta are neglected.For HEMP by transversely polarised virtual photons no corresponding proof of factorisation exists.Instead, phenomenological pQCD-inspired models [12][13][14][15] postulate the so-called  ⊥ factorisation, where  ⊥ denotes parton transverse momentum.In particular, the Goloskokov-Kroll (GK) model [13][14][15][16][17] allows for a simultaneous description of SDMEs as well as target and beam-spin asymmetries for HEMP induced by both longitudinally and transversely polarised virtual photons.
The chiral-even GPDs   and   are used to describe at leading twist the production of longitudinally polarised vector mesons by longitudinally polarised virtual photons.Here,  denotes a quark of a given flavour or a gluon.
In the GK model, the chiral-odd GPDs    and Ē   are used together with higher-twist effects in the threedimensional light-cone wave function to describe  *  →  0  transitions.These GPDs account for a helicity flip of the "active" quark and are hence related to the violation of s-channel helicity conservation (SCHC).The GPDs   and   , and also the pion-pole exchange mechanism treated in the GK model as one-boson exchange contribution, provide unnatural parity exchange (UPE) contributions to the transitions  *  →  0  and  *  →  0  .With this ansatz the GK model offers an explanation for the contrast between a substantial UPE contribution in exclusive  production and a small UPE contribution in exclusive  0 production.Spin density matrix elements are related to helicity amplitudes that describe transitions between specified spin states of virtual photon, target proton, produced vector meson, and recoil proton.In the case of an unpolarised nucleon target, SDMEs depend only on the helicities of virtual photon and produced meson, if the initial and final spin states of the proton are summed over.The interpretation of the measured SDME values is a rich field and in this paper we will address the following: the test of SCHC, the evaluation of UPE contributions, the determination of the phase difference between helicity amplitudes, and the calculation of the longitudinal-totransverse cross-section ratio.
However, among the quoted publications only three [23,26,29] are providing extensive sets of SDME values that were obtained through an analysis of three-dimensional angular distributions of  0 production and decay.Such an analysis allows for the determination of all 15 SDMEs that are not coupled to the beam polarisation ("unpolarised SDMEs"), as in Refs.[26,29].The complete set of 23 SDMEs, which includes the 8 SDMEs coupled to the beam polarisation ("polarised SDMEs"), was obtained for the first time by HERMES [23].The published results confirm the dominance of amplitudes for NPE transitions and the violation of the SCHC hypothesis that is observed for the transitions  *  →  0  .The present COMPASS results on SDMEs for exclusive  0 muoproduction have the potential to further constrain GPDs, in particular in conjunction with the published COMPASS results on SDMEs for exclusive  production [3].These additional constraints on GPD parameterisations are beyond those obtained from measurements of cross sections and spin asymmetries in HEMP.The COMPASS SDME results provide input to asses the role of chiral-odd, i.e., parton helicity-flip GPDs in exclusive vector-meson production, which are related to the mechanism of SCHC violation.

Theoretical formalism
Throughout this article, the theoretical formalism of SDMEs and helicity amplitudes introduced by K. Schilling and G. Wolf [1] is used, thereby following the notation from Refs.[3,31].

Definition of Spin Density Matrix Elements
In the hard exclusive process of vector-meson production on a nucleon  with helicity   ( ′  ) in the initial (final) state (Eq.2), the transition of a virtual photon  * with helicity   to a vector meson  with helicity   is described by helicity amplitudes     ′      , which depend on the three kinematic variables ,  2 , and Here  is the squared four-momentum transfer to the proton,  0 represents the smallest kinematically allowed value of || for given  2 and meson mass and  2 T the squared transverse momentum of the vector meson with respect to the virtual-photon direction.In the  * - centre-of-mass (CM) system, the vector-meson spin density matrix     ′  is related to the helicity amplitude     ′ amplitudes are related to helicity amplitudes as follows [1]: The asymptotic behaviour of amplitudes  at small  ′ was derived from angular-momentum conservation [32] and reads [2] Here and in the following  denotes the proton mass.Equations (6)(7)(8) show that double-helicity-flip amplitudes with |  −   | = 2 are suppressed at least by a factor of √  ′ /, and that their contributions to SDMEs are suppressed by  ′ / 2 .

Introducing the notation
and the symmetry properties [1,23] of the amplitudes , Eq. ( 9) becomes Note that the first product  * on the right-hand side represents contributions from NPE amplitudes without nucleon-helicity flip, while the second product of NPE amplitudes  * includes a nucleon-helicity flip.The relations for the UPE amplitudes can be written in an analogous way.For brevity, the nucleon-helicity indices will be omitted for amplitudes with   =  ′  in the following, i.e., The assumption that there exist only diagonal  * →  transitions (  =   ) is usually referred to as hypothesis of -channel helicity conservation.

Experimental access to SDMEs
Spin density matrix elements are extracted from COMPASS data in exclusive muoproduction of  0 mesons (Eq.( 1)).The SDMEs are fitted as parameters of the three-dimensional angular distribution W + (Φ, ,

Φ φ
lepton-scattering plane Fig. 1: Definition of angles in the process  →  ′  ′  0 with  0 →  +  − .Here, Φ is the angle between the  0 production plane and the lepton scattering plane in the centre-of-mass system of the virtual photon and the target nucleon, while  is the angle between the  0 production and decay planes.The variable Θ is the polar angle of the decay  + in the  0 meson rest frame.cos Θ), which is defined below, to the corresponding experimental distribution.The angles and reference frames for the production and decay process  →  ′  0 ( 0 →  +  − ) are shown in Fig. 1.
The right-handed "hadronic CM system" of virtual photon and target nucleon is identical to that used in Ref. [3].The directions of axes of the hadronic CM system and the  0 -meson rest frame coincide with the directions of axes of the helicity frame [1,23,33].The angles Φ, , and Θ involved in  0 -meson production and decay are defined as follows [23].The azimuthal angle Φ between  0 -meson production plane and lepton scattering plane in the hadronic CM system is given by: Here, k, k ′ , q = k − k ′ , and v are the three-momenta of the incoming and outgoing lepton, the virtual photon, and the  0 meson respectively.The azimuthal angle  between  0 -meson decay and production planes is defined by: where p  + is the three-momentum of the positively charged decay pion  + in the hadronic CM system.The polar angle Θ of the decay  + in the vector-meson rest frame, with the -axis aligned opposite to the outgoing nucleon three-momentum P ′ and the -axis directed along P ′ × q, is given by: where P  + is the three-momentum of the positively charged decay pion in the vector-meson rest frame.
The angular distribution W + is decomposed into contributions that are not coupled (W  -unpolarised) or coupled (W  -polarised) to the longitudinal beam polarisation  b : Since the data were collected with longitudinally polarised muon beams, both unpolarised and polarised SDMEs can be accessed, allowing the extraction of 15 unpolarised SDMEs from W  : and the extraction of 8 polarised SDMEs from W  : The virtual-photon polarisation parameter , which represents the ratio of fluxes of longitudinal and transverse virtual photons, is given by: where The symbols ,  and  denote the four-momenta of target proton, virtual photon and incident lepton, respectively.The energy of virtual photon and incident lepton in the target rest frame is denoted by  and , respectively.

Experimental setup and data selection
The fixed-target experiment COMPASS is located at CERN in the M2-beamline of the Super Proton Synchrotron (SPS).The experiment consists of a versatile setup that can use variety of targets, detectors and make use of different beams.It uses a two-stage spectrometer with a number of tracking and particle identification detectors placed over a length of approximately 60 m.Each stage of the spectrometer is built around one of the two dipole magnets (SM1 and SM2).The first stage covers large scattering angles up to 180 mrad, and the second stage smaller scattering angles down to 0.5 mrad.More detailed descriptions of the COMPASS experiment can be found in Refs.[34][35][36].
In this paper, we analyse data collected during 4 weeks in 2012 that were dedicated to the pilot run of a program designed to study GPDs and hadron tomography through Deeply Virtual Compton Scattering (DVCS) and HEMP processes.The experiment made use of the 160 GeV/ polarised muon beam and an unpolarised liquid-hydrogen target.The target was surrounded by a time-of-flight (TOF) system for the detection of recoil protons.The use of the recoil-proton detector (RPD) is important in the studies of the DVCS process, but for HEMP processes, like  0 -meson production, it restricts the kinematic coverage, so that the RPD information is not used in the present analysis.
The muon beam originates from in-flight decays of pions and kaons produced by SPS protons impinging on a primary target.Due to the weak nature of the decay, the muon beam is naturally polarised.The beam is negatively polarised for  + and positively polarised for  − , and the achieved absolute value of polarisation is (80 ± 4)%.The data has been taken using both  + and  − beams.The SM1 and SM2 spectrometer magnets polarities were changed accordingly with beam charge to ensure equal acceptance of the COMPASS spectrometer in both cases.
The analysis is focused on the process which defines the topology of the accepted events.The events are required to have two hadron tracks of opposite charge and one reconstructed vertex inside the target with incoming and outgoing muon associated.The outgoing muon is required to traverse more than 15 radiation lengths of material and to have the same charge as the incoming muon.Charged hadron tracks are identified by requiring to traverse less than 10 radiation lengths of material.

Kinematic selections
In order to select exclusively produced  0 mesons, events are required to meet the following kinematic constraints: -1.0 (GeV/) 2 <  2 < 10.0 (GeV/) 2 , which selects the region of perturbative QCD (lower limit) and suppresses background from hadrons produced in DIS, hereafter referred to as "SIDIS background" (upper limit); - < 0.9 to minimize the effect of radiative corrections; - > 5.0 GeV/ 2 to avoid significant fluctuations in the cross section that appear in the lower  region because of the production of resonances; - > 20 GeV, which is the energy of the virtual photon in the laboratory frame; -0.01 (GeV/) 2 <  2 T < 0.5 (GeV/) 2 to remove events with poorly determined azimuthal angle (lower limit) and to suppress SIDIS background (upper limit); -  0 > 15 GeV/ to reduce the SIDIS background contribution, where   0 is the  0 momentum in the laboratory frame; -0.5 GeV/ 2 <   +  − < 1.1 GeV/ 2 to select  0 mesons, where   +  − is the two-pion invariant mass.
In order to select exclusively produced  0 mesons, the missing energy is used.Here  2 X = (  +  −   + −   − ) 2 is the missing mass squared,   + (  − ) the pion four-momenta and  the mass of the proton.In order to account for experimental resolution the selection −2.5 GeV <  miss < 2.5 GeV is applied.The distribution of the missing energy is shown in Fig. 3, where the exclusive peak in the experimental data appears within the selection limits.
After having applied all the selection requirements, the data set for physics analysis consists of 23785 events taken with the  + beam and 28472 events with the  − beam.

Invariant mass distribution
The two-pion invariant mass distribution is shown in Fig. 2. A clear  0 signal is observed.Background coming from exclusive production of  and its decay  →  +  − , where the kaons are misidentified as pions, is expected and seen at   +  − < 0.4 GeV/ 2 .Therefore, applying the selection 0.5 GeV/ 2 <   +  − < 1.1 GeV/ 2 removes this background.The distribution for the experimental data is compared to that for the reconstructed events of Monte Carlo (MC) events obtained with the HEPGEN++  0 generator, in the following denoted by HEPGEN [37,38].As in HEPGEN only exclusive  0 production is generated, while the data contains exclusive production of both resonant and non-resonant pi+pi-pairs, as well as their interference, a difference in shapes between the experimental and simulated distributions is observed.The effect is in agreement with the expectation from the Söding model [39], which predicts a sizeable interference between the small amplitude for non-resonant  +  − pair production and the large one for resonant  +  − production.The characteristic prediction of the model is the change of the sign of the interference term at the maximum of the  0 resonance from being positive at smaller values of   +  − to negative at larger   +  − values.
In order to evaluate the contribution of non-resonant  +  − pair production, the following procedure is used.The invariant mass distribution from HEPGEN is normalised to the data in the region 0.75 GeV/ 2 <   +  − < 0.77 GeV/ 2 .The difference between the integrals of the distributions for the data and HEPGEN over the full range 0.5 GeV/ 2 <   +  − < 1.1 GeV/ 2 is approximately equal to the global contribution of nonresonant production.This contribution, which includes the interference between the amplitudes for resonant and non-resonant production as well as the squared amplitude for the latter one, is equal to about 3% and hence neglected.
In addition, the  →  +  − channel, with a branching fraction of 1.5%, gives an irreducible background to the  0 channel.As the branching fraction is small and the contribution of  0 −  interference was found to be very small [21], the contribution of  →  +  − channel is neglected in this analysis.

Backgrounds for exclusive 𝝆 0 production i) SIDIS background
The largest background contribution is due to the SIDIS background, i.e., events with hadrons produced in DIS.In order to determine the fraction of SIDIS background in the selected  0 events, the  miss distribution is used as shown in Fig. 3.The procedure is described in detail in Refs.[40,41].The SIDIS background simulation is performed using the LEPTO 6.5.1 generator with the COMPASS tuning of parameters [42] and processed with the simulation of the COMPASS setup [43].The simulated events are selected using the same criteria as for the experimental data.In order to improve the agreement between LEPTO events and the data, the simulated The COMPASS Collaboration events are reweighted.For this purpose, events with the same-sign hadron pairs are selected.The reweighting is applied on a bin-by-bin basis to the  miss distribution with the following weight: Here  sc D ( miss ) ( sc MC ( miss )) is the number of events with same-sign hadron pairs selected from experimental (D) or simulated (MC) data.
The distribution of the reweighted LEPTO events is normalised to the experimental data in the background dominated region of 7 GeV <  miss < 20 GeV.It is shown in Fig. 3 as the blue points.The procedure estimates the background fraction  bg for the selected  0 to be 0.17 in the signal region −2.5 GeV <  miss < 2.5 GeV.However, it was found that the fraction of SIDIS background changes within the kinematic coverage of this measurement, in particular it is increasing with increasing  2 and  2 T and with decreasing .Therefore, the background fraction is estimated in each kinematic bin separately, resulting in values of  bg from 0.10 to 0.32 for the determination of the SDME values as functions of kinematic variables.

ii) Background from proton-dissociation processes
This background is due to the processes  * +  →  0 +  * with a baryon  * decaying into a system of hadrons.As observed by HERA experiments, the cross section for such process, when integrated over all  * states, is at  Fig. 3: The missing-energy distribution from experimental data (red points) compared to the distribution of SIDIS events from a LEPTO MC simulation (blue points).Each LEPTO MC event is reweighted by a  miss -dependent weight that is calculated using both experimental and simulated data with same-charge hadron pairs.See text for a detailed explanation.
The reweighted MC distribution is normalised to the data in the region 7 GeV <  miss < 20 GeV.The SIDIS-backgroundcorrected distribution for the data is shown as shaded histogram.The vertical lines at | miss | = 2.5 GeV indicate the limits of the exclusive region.The shown distributions are obtained using all event selections except the selection on  miss .
the level of 20% of the cross section for exclusive  0 production [30].In the present analysis such processes are suppressed by the applied selections, in particular on  miss and  2 T , which reduce the contribution from proton-dissociation by a factor of 2. As it is reported in Refs.[27,29] that the angular distributions for  0 decay and production in the exclusive and proton-dissociation channels are compatible, no correction is applied for proton-dissociation events.

iii) Other backgrounds
In addition to the background contributions discussed above and in Sec.4.2, several others were considered in Ref. [27].They may originate from the processes  →  +  −  0 ,  → ,  →  +  −  0 and  ′ →  +  −  0  0 .These events are expected to be strongly suppressed after applying the selections on  miss and  2 T and their contributions are neglected in the present analysis.

Unbinned Maximum Likelihood method
The method to determine SDMEs was described in Sec. 3. Equations (17,18,19) relate the angular distribution W to the 23 SDMEs      ′

𝑉
. In order to extract SDME values in this measurement, the Unbinned Maximum Likelihood (UML) method is used to fit the experimental three-dimensional angular distribution of  0 production and decay to the function W (R; Φ, , cos Θ), where R is the set of the 23 SDMEs.In the fit, the negative log-likelihood function is minimised.Here,  represents the number of selected events and N (R) is the likelihood normalisation factor defined as where  MC is the number of simulated  0 events generated by the HEPGEN generator [37,38].In order to simulate exclusive  0 production, the option of an isotropic three-dimensional angular distribution of  0 production and decay was chosen.The generated events are further processed with the simulation of the COMPASS setup [43].Identical selection requirements are applied as for the experimental data.

Background-corrected SDMEs
The above described procedure of SDME extraction does not account for the SIDIS background contamination.
In order to determine the background-corrected SDMEs, a two-step approach is applied.
In the first step, the 23 "background SDMEs" are determined using a parameterisation of the background angular distributions.The SIDIS background events simulated by the LEPTO generator are treated by the same method as described above.The UML fit is performed in the signal region, −2.5 GeV<  miss < 2.5 GeV according to Eq. ( 23), resulting in the set R bg of background SDMEs.
In the second step, the set R bg and the background fraction  bg determined in Sec.4.3 are used to extract the set R sig of the background-corrected SDMEs by fitting the negative log-likelihood function Here, N is the likelihood normalisation factor defined as

Statistical uncertainties of the observables depending on SDMEs
The statistical uncertainties of the observables, which depend on SDMEs and are discussed in Sec. 7, were calculated by propagating the statistical uncertainties of SDMEs and using their covariance matrix as obtained from the fit.The corresponding correlation matrix for the data in the total kinematic range is presented in Table A.5.

Systematic uncertainties
The systematic uncertainties of the measured SDME values are considered to arise from the following sources: i) Difference between  + and  − beam In the measurement,  + and  − beams were used, which were not identical in terms of intensity.The intensity of the  + beam was higher by a factor of approximately 2.7 than that of the  − beam.In order to account for a possible impact of this difference on the measured SDMEs, the latter are extracted separately for  + and  − data, and half of the difference between the two results is assigned as systematic uncertainty.
ii) Position of the  miss peak As observed in Fig. 3, the signal peak in the  miss distribution is not centred at zero, but rather slightly shifted towards negative values.The reason for this shift is a small imbalance between the measured incoming muon energy and energies of the measured final state particles in the spectrometer.Some SDME values depend on the position of this peak [44], hence a systematic uncertainty is assigned based on the difference between the SDMEs extracted without and with a correction of +0.25 GeV/ to the beam momentum to centre the  miss peak at zero.It was checked that this method of the beam momentum correction yields a similar systematic uncertainty as in case of the method that rescales the measured momenta of the final-state particles to centre the  miss peak position at zero.

iii) Dependence on the background angular distribution
The method to evaluate the background-corrected SDMEs, described in Sec.5.2, uses the LEPTO generated events for the estimation of the SIDIS background SDMEs in the signal region −2.5 GeV <  miss < 2.5 GeV.Note that the LEPTO generator was not tuned to reproduce the experimental angular distributions in the specific phase space of this analysis.In order to account for a possible source of uncertainty, another procedure was applied to estimate the background SDMEs using the backgrounddominated region 7.0 GeV<  miss < 20.0 GeV in the experimental data.The systematic uncertainty is assigned based on the difference between the two methods of evaluating the background SDMEs.

iv) Uncertainty in the determination of the background fraction
Another contribution to the systematic uncertainty of SDME values is related to the uncertainty of the background fraction determination.It is estimated to be about 1% based on the comparison of background fraction values that were evaluated using two different methods to normalise LEPTO MC results with respect to experimental data.The difference between the respective SDME values is taken as systematic uncertainty.

v) Sensitivity to the shapes of the kinematic distributions generated by HEPGEN
The SDME values can be sensitive to the shapes of the kinematic distributions generated by the HEPGEN generator.In order to check for such an effect, the SDMEs were extracted using modified HEPGEN weights so that the reconstructed MC  2 and  distributions match those of the experimental data.
Although the effect on the measured SDMEs is fairly small, the difference between the extraction using the original simulated sample and the one with reweighting is assigned as systematic uncertainty.
The effect of different non-exclusive backgrounds on the extracted values of SDMEs was studied in the earlier COMPASS analysis of exclusive  production [3].Two event samples were used for the extraction of SDMEs.
The first one was obtained by applying selections similar to those described in the present  0 analysis.For the second one the more restrictive selections using the information from the RPD were added, which lead to a reduction of the non-exclusive background by a factor of about 10.As a limited  2 T -range is covered by the RPD, the same limited kinematic region was used to compare the SDMEs obtained with and without RPD.The SDME values extracted from the two data samples were found to be consistent within statistical uncertainties.This observation confirms the correctness of the method to extract the SIDIS-background-corrected SDMEs as explained in Sec.5.2.
The contributions from the aforementioned sources i)-v) of systematic uncertainties are shown individually in the Appendix in Table A.1.The largest source is typically from group i) (the difference between the  + and  − beams), followed by group ii) (shift in the  miss peak position) and group iii) (dependence on the background angular distribution).The systematic uncertainties arising from the above discussed sources are added in quadrature to obtain the total systematic uncertainty.The 23 SDMEs measured over the entire COMPASS kinematic region are given in Table 1 together with their statistical and total systematic uncertainties.For most SDMEs the total systematic uncertainty is larger than the statistical uncertainty.
The systematic uncertainties of the observables that are discussed in Sec.7 were estimated using a procedure analogous to that used for SDMEs.This means that for a given observable the contribution to its systematic uncertainty from each of the five sources indicated in points i) -v) is calculated directly as difference of the values of the concerned observable obtained when using the corresponding two sets of SDME values.For a given observable, the systematic uncertainties arising from the above discussed sources are added in quadrature to obtain the total systematic uncertainty.
The dominant contributions to the SDMEs in class A are related to the squared amplitudes for transitions from longitudinally polarised virtual photons to longitudinally polarised vector mesons,  *  →   , and from transversely polarised virtual photons to transversely polarised vector mesons,  *  →   .The former ones appear in the SDME  04 00 , and the latter ones in the SDMEs  1 1−1 and Im  2 1−1 , which approximately mirror each other value (see Fig. 4 and Table 1).The dominant terms in class B correspond to the interference between amplitudes for the two aforementioned transitions.The SDMEs of this class allow the determination of the phase difference between the amplitude  11 for

Test of the SCHC hypothesis
In the case of SCHC only three amplitudes,  00 ,  11 and  11 , may be different from zero.As a consequence all SDMEs of classes A and B may not vanish, while SDMEs from classes C, D, and E should be equal to zero.Six of the SDMEs in classes A and B have to fulfil the following relations [1] Using the extracted SDMEs one obtains: where total uncertainties are quoted.While the measurements of the first and the third relation in Eq. ( 27) are consistent with the expectation, a tension is observed for the second relation, which may indicate a contribution of single-helicity-flip amplitudes.In the case of the first relation only the contributions from squared small double-helicity-flip amplitudes violate SCHC.For the two other relations the contributions that violate SCHC are related to small terms corresponding to the interference of two single-helicity-flip amplitudes as well as the interference of the helicity-conserving amplitude  00 and the double-helicity-flip amplitude  1−1 .
However, for the transitions  *  →   of class C the non-zero values of five unpolarised SDMEs indicate a clear SCHC violation.In the GK model [15], these SDMEs are related to the chiral-odd GPDs  T and ĒT coupled to the higher-twist wave function of the meson.The kinematic dependences of these SDMEs, as presented in Section 6, may help to further constrain the model.

Contribution of the helicity-flip NPE amplitudes
The contributions of non-zero helicity-single-flip and helicity-double-flip amplitudes to the cross section can be quantified by the ratios    of the helicity-flip amplitudes    to the square root of the sum of all amplitudes squared Here, the normalisation factor N is given by N = N  + N  with The ratios    can be expressed in terms of SDMEs as shown in Ref. [23].
The quantity  10 , which is related to the amplitude  10 describing the transition  *  →  0  , is approximated by For the helicity-double-flip amplitude  1−1 describing the transition  * − →  0  the quantity  1−1 is given by In Fig. 8 the dependence of the quantities  01 ,  10 and  1−1 on  2 ,  2 T and  is presented.For  01 , values significantly different from zero are observed, while for  10 and  1−1 they are much smaller.This observation is consistent with the different degrees of SCHC violation seen for SDMEs in classes C, D and E. A squared ratio  2   represents the fractional contribution from amplitude    to the full cross section.Therefore the quantity represents the fractional contribution of helicity-flip NPE amplitudes to the cross section.The value of  2 NPE for the COMPASS entire kinematic range is small, equal to 0.023 ± 0.002 ± 0.004.

UPE contribution in exclusive 𝝆 0 meson production
By examining a linear combination of SDMEs such as the presence of a UPE contribution can be tested.The quantity  1 is expressed in terms of helicity amplitudes as  thus a positive value of  1 would indicate a non-zero contribution from UPE transitions.For the entire kinematic region of COMPASS  1 is equal to 0.047 ± 0.010 ± 0.029, which indicates a small UPE contribution.Additional information on UPE amplitudes can be obtained from the SDME combinations which in terms of helicity amplitudes can be combined into The value of  2 + iu 3 can vanish despite of the existence of UPE contributions.For COMPASS  2 = −0.008± 0.002 ± 0.013 and  3 = −0.010± 0.018 ± 0.037 are obtained, which are consistent with zero at the present accuracy of the data.In Fig. 9 the dependence of the quantities  1 ,  2 and  3 on  2 ,  2 T , and  is presented.The quantities  1 ,  2 and  3 are small and compatible with zero within experimental uncertainties.
The UPE fractional contribution to the cross section is given as where the contributions of the amplitudes  01 and  1−1 was neglected for the approximate relation to  1 .The value of Δ UPE ( 0 ) for the entire kinematic range is 0.024 ± 0.005 ± 0.014.
Altogether, for exclusive  0 production at COMPASS the contribution of UPE is very small.This is in sharp contrast to the significant UPE contribution observed by COMPASS [3] for exclusive  production in a similar kinematic range.There, this contribution is large over the entire kinematic range, Δ UPE () = 0.415 ± 0.037 ± 0.025.The UPE dominates  production at small  values and its contribution decreases with increasing  without vanishing towards large  values accessible in COMPASS.In the GK model, UPE is described by the GPDs   and   (non-pole), and by the pion-pole contribution treated as a one-boson exchange [17].The large difference in size of the UPE contributions for  and  0 production is mostly explained by the difference between  −  and  −  0 transition form factors, with the former one being about three times larger than the latter [17].

The NPE-to-UPE asymmetry of the transverse cross section for the transition 𝜸 *
→   Another observable that is sensitive to the relative contributions of UPE and NPE amplitudes is the NPE-to-UPE asymmetry of the transverse differential cross section for the transition  *  →   .It is defined [17] as where the superscript  and  denotes the part of the cross section that is fed by NPE and UPE transitions, respectively.The value of  obtained in the entire kinematic region is 0.887 ± 0.016 ± 0.029, which indicates that the NPE contribution dominates when averaged over the whole kinematic range of COMPASS.The kinematic dependences of the asymmetry  are shown in Fig. 10.A small UPE contribution is observed only at small values of  and it becomes compatible with zero at larger .No significant  2 and  2 T dependences of the asymmetry are observed.
The COMPASS results for exclusive  production [3] on the asymmetry  and its kinematic dependences exhibit a different behaviour.For the whole kinematic region the value is compatible with zero, () = −0.007± 0.076 ± 0.125, which indicates that the UPE and NPE contributions averaged over the whole kinematic range are of similar size.The UPE contribution dominates at small values of  and decreases with increasing .At large values of  the NPE contribution for  production becomes dominant, while a non-negligible UPE contribution still remains.
A semi-quantitative explanation of the difference between the values of asymmetry  measured for  0 and  production is possible by considering only the dominant contributions to the UPE and NPE cross sections for each process.In the framework of GK model such a contribution to the UPE cross sections is due to pion-pole exchange.Due to the difference between  −  0 and  −  transition form factors, which was mentioned in Sec.7.3, the UPE cross section for  0 is about 9 times smaller than that for  production.For the NPE cross section, the dominant contribution is related to the gluon and sea-quark GPDs  and their relative contributions given in Ref. [45] imply that the cross section of exclusive  0 production is about 9 times larger than that for .Taken together, UPE contributions are close to zero for  0 and the  value approaches unity.

Longitudinal-to-transverse cross-section ratio
The longitudinal-to-transverse virtual-photon cross-section ratio T Fig. 11:  2 ,  2 T and  dependences of two estimates,  ′ and , of the longitudinal-to-transverse cross-section ratio .Inner error bars represent statistical uncertainties and outer ones statistical and systematic uncertainties added in quadrature.For better visibility the data points for  are presented with a small horizontal off-set.
is one of the most important observables in the study of light VM production since it is sensitive to the interaction dynamics.In order to evaluate  the quantity  ′ is commonly used: Using expressions defining  04 00 and 1 −  04 00 in terms of helicity amplitudes [1,23], the quantity  ′ may be interpreted [3] as longitudinal-to-transverse ratio of "effective" cross sections for the production of vector mesons that are polarised longitudinally or transversely irrespective of the virtual-photon polarisation.In case of SCHC,  ′ is equal to .In spite of the observed violation of SCHC, the approximate relation  ≈  ′ was used in most of the previous measurements.For the entire kinematic region, the ratio  ′ is found to be 0.980 ± 0.014 ± 0.089.The kinematic dependences of  ′ are shown in Fig. 11.
In order to evaluate the effect of helicity-changing amplitudes on the estimate of the longitudinal-to-transverse cross-section ratio one can use a quantity  that is defined by following relation: [23] where The quantity  can be approximately estimated as Here  01 and  10 are evaluated using Eqs.( 31) and ( 32), and the small contributions of the helicity-flip UPE amplitudes  01 and  10 in Eq. ( 45) are neglected.
For the entire kinematic region the ratio  is found to be 0.907 ± 0.014 ± 0.074.The values of  as functions of kinematic variables are shown in Fig. 11 for comparison with  ′ .The values of  and  ′ with their statistical and total uncertainties are also given in Table 3 for the kinematic bins and for the entire kinematic region.The Using  ′ as an estimate of  introduces an additional uni-directional systematic uncertainty on , which is due to the assumption of SCHC.It is estimated from  −  ′ and found to be about −0.07 on average.Thus for an estimate of the ratio  (defined by Eq. ( 42)) it is preferable to use the quantity , i.e.  ≈ , while  ′ values may be used when comparing to earlier measurements that assumed SCHC.
A strong increase of the   /  ratio with increasing  2 is observed, while the  2 T and  dependences are weaker (see Fig. 11).The  2 dependence of  indicates that   becomes dominant at  2 larger than about 2 (GeV/) 2 .

Phase difference between amplitudes 𝑻 11 and 𝑻 00
The phase difference between the amplitudes  11 and  00 can be evaluated as follows [23]: The result is || = 19.6 ± 0.9 ± 3.9 degrees.The sign of  can be obtained from sin [23] that depends on the polarised SDMEs: Using Eq. ( 48) one determines that the sign of  is positive and  = 12.9 ± 2.1 ± 0.8 degrees.Both results on  are compatible within large total uncertainties.For comparison to other experiments, which used cos  to estimate the phase difference, we use the positive value of  = 19.6 ± 0.9 ± 3.9 degrees.

Comparison to other experiments
We compare the COMPASS results on SDMEs and related observables to those from the experiments that measured extensive sets of SDMEs for exclusive  0 electroproduction: the HERMES [23], H1 [26,27] and ZEUS [29,30] experiments.Only the HERMES experiment measured the complete set of 15 unpolarised and 8 polarised SDMEs, while H1 and ZEUS measured 15 unpolarised SDMEs.Compilations of selected SDMEs-related observables, including other experiments, can be found e.g. in Refs.[23,27].
The complete sets of SDMEs from HERMES and COMPASS for their entire kinematic ranges are compared in Fig. 12.The kinematic range for HERMES is 1.0 (GeV/) 2 <  2 < 7.0 (GeV/) 2 , 3.0 GeV/ 2 <  < 6.3 GeV/ 2 , | ′ | < 0.4 (GeV/) 2 , while that for COMPASS 1.0 (GeV/) 2 <  2 < 10.0 (GeV/) 2 , 5.0 GeV/ 2 <  < 17.0 GeV/ 2 , 0.01 (GeV/) 2 <  2 T < 0.5 (GeV/) 2 .The ranges of  2 and the momentum transferred to the recoil proton are similar, but the  ranges overlap only marginally and the COMPASS range extends significantly towards large  values.In consequence, the contribution of gluons and sea quarks involved in exclusive meson production is higher by a factor of about 2.5 in COMPASS, while the remaining contribution (from valence quarks and interference terms) is a little smaller than in HERMES (see e.g.Ref. [46]).Significant differences are observed for the SDME  04 00 , which is proportional to the square of the leading helicity-conserving amplitude  00 , and for the unpolarised SDMEs from class C.
A more detailed comparison of selected observables in similar kinematic ranges for both experiments is presented in Table 4.The published HERMES results for the entire kinematic range (⟨⟩ = 4.8 GeV/ 2 ) are compared to the COMPASS results for the lowest  bin (⟨⟩ = 7.0 GeV/ 2 ).The quoted uncertainties are the total ones.For most of the observables the results from both experiments are compatible within one standard deviation, except  01 and  1−1 , for which the agreement is within two standard deviations.
The comparison of the  2 dependence of  between COMPASS and HERMES is not straightforward because  11), the estimates of  from both experiments are compatible within experimental uncertainties as shown in Fig. 13.
Deviations from zero are observed for five unpolarised SDMEs from classes A and B, which depend on the helicity-conserving amplitudes  00 and  11 .All other SDMEs are compatible with zero except  5 00 , which indicates the violation of SCHC for  *  →  0  transitions.In order to quantify the size of SCHC violation, the ratios were evaluated.The approximate expressions for τ  are given in Ref. [29].In contrast to Eqs. (31,32,33) those expressions rely on the assumptions of zero phase difference between the considered amplitude    and the corresponding dominant amplitude ( 00 or  11 ), neglecting UPE contributions and assuming  = 1.For the kinematic ranges of the HERA experiments the values of τ01 are equal to 0.079 ± 0.026 for ZEUS and 0.08 ± 0.03 for H1, while the values of τ10 and τ1−1 are compatible with zero within experimental uncertainties.These results indicate that the helicity-flip amplitude  01 does not vanish even at the highest available energies.
results with  2 > 1.0 (GeV/) 2 and with moderate to large values of .The HERMES and COMPASS results are corrected for contributions of the spin-flip amplitudes  01 and  10 .For those from H1 the contribution of  01 is taken into account, whereas the SCHC approximation is used for the other data.Despite small differences due to different treatments of small contributions of spin-flip amplitudes, and also due to a possible weak  dependence, all the results consistently show a main characteristic feature, i.e. the fast increase of  as a function of  2 within the wide energy range, from the fixed target experiments to the HERA collider measurements.
In leading-order pQCD and for  = 0 the ratio  is predicted to be  =  2 / 2  [47], where   is the mass of the produced vector meson.The experimental data on  for exclusive  0 ,  and / production confirm the scaling with   , but they lie systematically below this prediction (see, e.g., Fig. 38 from Ref. [27]).Deviations from this dependence, which become more pronounced as  2 increases, are due to effects of QCD evolution and quark transverse momentum [45,48].
In the framework of the colour dipole model, different transverse sizes are predicted for virtual  q pair fluctuations originating from longitudinally and transversely polarised virtual photons, which leads to different kinematic dependences of   ,   and .The transverse size of these colour dipoles is on average smaller for longitudinal photons than for transverse ones.This results in a more shallow  (or  2 T ) dependence of the cross section for the longitudinal photons.In the unseparated cross section this effect leads to a decrease of the value of the -slope parameter with increasing  2 (see e.g.Refs.[27,30]).

Summary
Using exclusive  0 meson muoproduction on the proton, we have measured 23 SDMEs at the average COMPASS kinematics, ⟨ 2 ⟩ = 2.4 (GeV/) 2 , ⟨⟩ = 9.9 GeV/ 2 and ⟨ 2 T ⟩ = 0.18 (GeV/) 2 .The SDMEs are extracted in the kinematic region 1.0 (GeV/) 2 <  2 < 10.0 (GeV/) 2 , 5.0 GeV/ 2 <  < 17.0 GeV/ 2 and 0.01 (GeV/) 2 , <  2 T < 0.5 (GeV/) 2 , which allows us to study their  2 ,  2 T and  dependences.Several SDMEs that are depending on amplitudes describing  *  →  0  transitions indicate a considerable violation of the SCHC hypothesis.These SDMEs are expected to be sensitive to the chiral-odd GPDs  T and ĒT , which are coupled to the higher-twist wave function of the meson.A particularly prominent effect is observed for the SDME  5 00 , which strongly increases with increasing  2 and  2 T .Using specific observables that are constructed to be sensitive to the relative contributions from transitions with unnatural-parity exchanges, such as  1 and the NPE-to-UPE asymmetry for the transverse cross section, we observe a dominance of NPE exchanges.The UPE contribution is very small and compatible with zero within experimental uncertainties.
The COMPASS results presented in this paper are obtained in the kinematic range that partially overlaps with the kinematic range of HERMES experimental data on SDMEs for exclusive  0 electroproduction, but extends considerably towards higher  values.In the overlap region the results from both experiments are compatible.Important data on the behaviour of helicity amplitudes at very large values of  2 and  are provided by ZEUS and H1.Characteristic features are the dominant contribution of the amplitude describing the transition  *  →  0  , which increases with increasing  2 , as well as negligible contributions of spin-flip and UPE amplitudes.These results allow one to better constrain extrapolations of trends observed at fixed target experiments.
The present results provide important input for modelling GPDs, in particular they may help to better constrain the chiral-odd GPDs and the amplitudes for UPE transitions in exclusive  0 leptoproduction.

Fig. 2 :
Fig. 2: Distribution of the  +  − invariant mass for experimental data (shaded histogram) and HEPGEN MC (open histogram).The distributions are obtained applying all event selections except the selection on the invariant mass.The invariant mass distribution from HEPGEN is normalised to the data in the region 0.75 GeV/ 2 <   +  − < 0.77 GeV/ 2 .The vertical lines indicate the applied limits.

Table 1 :
The 23 unpolarised and polarised SDMEs for the entire COMPASS kinematic region, shown in the same order as in Fig.4for classes A to E. The first uncertainties are statistical, the second systematic.± 0.0085 ± 0.0043

)Fig. 5 :
Fig. 5:  2 dependence of the measured 23 SDMEs.The capital letters A to E denote the class, to which the SDME belongs.Inner error bars represent statistical uncertainties and outer ones statistical and systematic uncertainties added in quadrature.

Fig. 6 :
Fig. 6:  2 T dependence of the measured 23 SDMEs.The capital letters A to E denote the class, to which the SDME belongs.Inner error bars represent statistical uncertainties and outer ones statistical and systematic uncertainties added in quadrature.

Fig. 7 :
Fig. 7:  dependence of the measured 23 SDMEs.The capital letters A to E denote the class, to which the SDME belongs.Inner error bars represent statistical uncertainties and outer ones statistical and systematic uncertainties added in quadrature.

Fig. 10 :
Fig. 10:  2 ,  2 T and  dependences of the NPE-to-UPE asymmetry of the transverse cross section for the transition  *  →   .Inner error bars represent statistical uncertainties and outer ones statistical and systematic uncertainties added in quadrature.

Table 2 .
transitions and the amplitude  00 for  *  →   transition (cf Sec.7.6).In class C, the main terms in most of the SDMEs are proportional to the interference between the helicity-flip amplitude  01 , describing  *  →   transitions, and the large helicity-conserving amplitudes, either  11 (for Re  04 10 , Re  1 10 , Im  2 10 , Im  3 10 ) or  00 (for  5 00 ,  8 00 ).The dominant terms in the SDMEs of classes D and E are proportional to the interference between the amplitude  11 and small amplitudes describing  *  →   and  *  →  − transitions, respectively.The experimental uncertainties of the polarised SDMEs are in most of the cases larger than those of the unpolarised ones because the lepton-beam polarisation is smaller than unity (| b | ≈ 80%), and in the expressions for the angular distributions (see Eq. (19)) they are multiplied by the small factor | b | √ 1 − , where  ≈ 0.90.The SDMEs values extracted in four kinematic bins of  2 ,  2 T , or  are shown in Figs. 5, 6 and 7.The limits of the kinematic bins and the mean values of kinematic variables in each bin are given in The value of the SDME  04 00 , which corresponds to the fractional contribution of | 00 | 2 from longitudinally polarised virtual photons to the cross section, increases with  2 and  2 T , while the opposite trend is observed for the absolute values of the SDMEs  1 1−1 and Im  2 1−1 , which represent the fractional contribution of | 11 | 2 from transversely polarised virtual photons.In class C a sizeable increase of  5 00 with  2 is observed.As a consequence of angular-momentum conservation the helicity single and double-flip amplitudes should vanish as  2 T → 0, which is consistent with the measured  2 T -dependence of SDMEs in classes C, D and E. No clear -dependence is observed for any of 23 SDMEs.

Table 2 :
Kinematic binning and mean values for kinematic variables.

Table 3 :
Values of the ratios  and  ′ in kinematic bins and for the entire kinematic region.The first uncertainty values correspond to the statistical errors, and the second ones to the statistical and systematic uncertainties added in quadrature.both observables are dominated by the systematic ones, which are similar for both cases.

Table 4 :
[23]arison of selected observables measured by HERMES and COMPASS in similar kinematic regions.The HERMES results for the proton target[23]are integrated over the entire kinematic region.The COMPASS results are given for 5.0 <  < 7.3 GeV/ 2 .The total uncertainties are given.
the results are integrated over different  ranges for each experiment.Despite the moderate  dependence of  observed by COMPASS (cf Fig.

Table A . 5 :
The correlation matrix for the measured 23  0 SDMEs for data in the entire kinematic range.
.de) (Germany) e Also at ORIGINS Excellence Cluster, 85748 Garching, Germany f Also at Institut für Theoretische Physik, Universität Tübingen, 72076 Tübingen, Germany g Present address: NISER, Centre for Medical and Radiation Physics, Bubaneswar, India h Also at University of Eastern Piedmont, 15100 Alessandria, Italy h1 Supported by the Funds for Research 2019-22 of the University of Eastern Piedmont i Also at Chubu University, Kasugai, Aichi 487-8501, Japan j Also at KEK, 1-1 Oho, Tsukuba, Ibaraki 305-0801, Japan k Also at Dept. of Physics, Pusan National University, Busan 609-735, Republic of Korea k1 Also at Physics Dept., Brookhaven National Laboratory, Upton, NY 11973, USA l Also at Dept. of Physics, National Central University, 300 Jhongda Road, Jhongli 32001, Taiwan m Also at Dept. of Physics, National Kaohsiung Normal University, Kaohsiung County 824, Taiwan † Deceased A Supported by MEYS, Grants LM2015058, LM2018104 and LTT17018 and Charles University Grant PRIMUS/22/SCI/017 (Czech Republic) B Supported by BMBF -Bundesministerium für Bildung und Forschung (Germany) C Supported by B. Sen fund (India) D Supported by the Israel Academy of Sciences and Humanities (Israel) E Supported by MEXT and JSPS, Grants 18002006, 20540299, 18540281 and 26247032, the Daiko and Yamada Foundations (Japan) F Supported by NCN, Grant 2020/37/B/ST2/01547 (Poland) G Supported by FCT, Grants CERN/FIS-PAR/0022/2019 and CERN/FIS-PAR/0016/2021 (Portugal) H Supported by the Ministry of Science and Technology (Taiwan) I Supported by the National Science Foundation, Grant no.PHY-1506416 (USA)