Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://www.atnf.csiro.au/people/Sebastian.Haan/thesis_haan.pdf
Äàòà èçìåíåíèÿ: Wed Feb 1 10:35:58 2012
Äàòà èíäåêñèðîâàíèÿ: Tue Oct 2 06:06:30 2012
Êîäèðîâêà:

Ïîèñêîâûå ñëîâà: hep-ex
Diploma Thesis

Deeply Virtual Compton Scattering on Neon
Sebastian Haan ¨ UNIVERSITAT LEIPZIG August 16, 2005


Sup ervisor at the University of Leipzig: Prof. Dr. T. Butz Sup ervisor at DESY: Dr. W.-D. Nowak

2


The HERMES Experiment

Deutsches Elektronen-Synchroton DESY

3


4


Contents

1 Intro duction 2 Generalized Parton Distributions and Deeply Virtual Compton 2.1 Prop erties of Generalized Parton Distributions . . . . . . . 2.2 The Energy-Momentum Tensor and the Spin of the Nucleon 2.3 Kinematics of Deeply Virtual Compton Scattering . . . . . 2.4 Beam-Spin and Beam-Charge Asymmetries . . . . . . . . . 2.5 Deeply Virtual Compton Scattering on Nuclei . . . . . . . . 3 The 3.1 3.2 3.3 HERMES Exp eriment The HERA Storage Ring Facility . . . The Internal Gas Target . . . . . . . . The HERMES Sp ectrometer . . . . . . 3.3.1 The Tracking System . . . . . 3.3.2 The Particle Identification . . . 3.3.3 The Trigger System . . . . . . 3.4 Data Acquisition . . . . . . . . . . . . 3.4.1 Data Taking . . . . . . . . . . 3.4.2 Data Acquisition and Structure . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Scattering ...... ...... ...... ...... ...... . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

9 11 12 15 17 19 22 28 28 30 31 33 33 35 36 36 37 38 38 38 38 39 40 41 44 47 49

. . . . . . . . . . . . . . . . . . . . . . .

. . . . . . . . . . . . . . . . . . . . . . .

4 Data Analysis 4.1 Event Selection . . . . . . . . . . . . . . . . 4.1.1 Data Sample . . . . . . . . . . . . . 4.1.2 Data Quality . . . . . . . . . . . . . 4.1.3 DIS Event Selection . . . . . . . . . 4.1.4 Single-photon Selection . . . . . . . 4.1.5 Exclusive Event Selection . . . . . . 4.2 Extraction of the Beam-Spin Asymmetry . 4.3 Hydrogen Data . . . . . . . . . . . . . . . . 4.4 Separation of Coherent and Incoherent Part

5 Systematic Studies 53 5.1 Monte Carlo Studies . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53 5.1.1 The HERMES Monte Carlo Pro duction . . . . . . . . . . . . . . . 53

5


5.1.2 Coherent and Incoherent Contribution . . . . . . 5.2 The Influence of the Separation Point on the BSA . . . 5.3 Comparison of the Calculation Metho d with Proton and 5.4 Normalization metho d . . . . . . . . . . . . . . . . . . . 5.5 Study of different Fit Metho ds . . . . . . . . . . . . . . 5.5.1 The Influence of Binning . . . . . . . . . . . . . . 5.5.2 The Fit Function . . . . . . . . . . . . . . . . . . 5.5.3 The Anti-symmetrization Fit Metho d . . . . . . 5.6 The Influence of the Calorimeter on the Results . . . . . 5.6.1 Study of the Photon Energy Reconstruction . . . 5.6.2 Study of the Photon Position Reconstruction . . 5.7 Smearing Effects . . . . . . . . . . . . . . . . . . . . . . 5.8 Determination of the Background Contribution . . . . . 5.9 Combined Systematic Uncertainties . . . . . . . . . . . . 5.10 Study of the Constant Term . . . . . . . . . . . . . . . . 5.10.1 The Time Dep endence of p0 . . . . . . . . . . . . 5.10.2 Study of Systematic Influences . . . . . . . . . .

... ... Neon ... ... ... ... ... ... ... ... ... ... ... ... ... ...

... ... Mass ... ... ... ... ... ... ... ... ... ... ... ... ... ...

... ... .. ... ... ... ... ... ... ... ... ... ... ... ... ... ...

. . . . . . . . . . . . . . . . .

53 54 56 57 58 58 59 61 63 63 63 65 65 66 67 68 69

6 Summary of the Results 72 6.1 Coherent Sample . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72 6.2 Incoherent Sample . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75 7 Interpretation in the Light of Theoretical Mo dels 77 7.1 Comparison with Hydrogen Data . . . . . . . . . . . . . . . . . . . . . . . 77 7.2 Mo del Calculations and Interpretation . . . . . . . . . . . . . . . . . . . . 81 8 Summary and Outlo ok 9 App endix 83 85

6


Abstract
The measurement of Deeply Virtual Compton Scattering (DVCS) provides access to the spin structure of nucleons. Its description within the theoretical framework of generalized parton distributions offers a p ossibility to determine the total angular momentum carried by the quarks in the nucleon. In this thesis a cross section asymmetry with resp ect to the b eam-helicity has b een measured for hard exclusive electropro duction of photons from a neon target. The resulting asymmetry app ears in the distribution of the pro duced photon in the azimuthal angle around the virtual-photon direction. Attributed to the interference b etween the BetheHeitler pro cess and the DVCS pro cess, the asymmetry gives access to the latter at the amplitude level. The DVCS reaction on nuclei pro ceeds through two different pro cesses the coherent pro cess that involves the nucleus as a whole and the incoherent pro cess as the reaction on a single nucleon. The b eam-helicity asymmetry and its dep endences are studied for the coherent and incoherent sample, separately. The data presented has b een accumulated by the HERMES exp eriment at DESY, scattering the HERA 27.6 GeV p olarized p ositron b eam off an unp olarized neon target.

Kurzfassung
Schnelle geladene Teilchen hab en sich schon seit l¨ angerer Zeit als nutzliches Instru¨ ment zur Untersuchung der Struktur von Nukleonen erwiesen. Bei dem Prozess der tiefvirtuellen Comptonstreuung (DVCS) werden durch die Streuung von ho chenergetischen Elektronen o der deren Antiteilchen (Positronen) an einem Parton (Quark o der Gluon) im Inneren des Nukleons Photonen erzeugt. Seit einigen Jahren existiert der Formalismus der sogenannten "`generalisierten Partonverteilungen"' (GPDs) zur Beschreibung derartiger Prozesse. Diese GPDs enthalten zus¨ atzlich zur Beschreibung der Partonstruktur des Nukleons, ausgedruckt durch die ¨ einfachen Partonverteilungen (PDFs), no ch zus¨ atzliche Information ub er Korrelationen ¨ der Quarks und Gluonen untereinander. Die Messung von GPDs erm¨ ht im Prinzip oglic die Bestimmung des Gesamtdrehimpulsanteils der Quarks und Gluonen innerhalb des Nukleons. In der vorliegenden Arb eit wurde die Asymmetrie in Bezug auf die Polarisationsrichtung des Positronstrahles fur die exklusive Erzeugung von Photonen an Neonkernen gemessen. ¨ Die Asymmetrie tritt in der Verteilung der Photonen b ezuglich des azimuthalen Winkels ¨ zwischen Positronenstreueb ene und Photonenerzeugungseb ene auf. Dab ei kommt es ¨ zu einer Ub erlagerung mit dem Bethe-Heitler (BH) Prozess, der auf der Abstrahlung von Photonen durch Bremsstrahlung b eruht. Aufgrund der Interferenz b eider Prozesse, oglic arteil der Interferenz-Amplitude zu b estimDVCS und BH, ist es m¨ h, den Imagin¨ men, der wiederum in Beziehung zu den generalisierten Partonenverteilungen steht. Im Fall von DVCS an Kernen ist zwischen zwei verschiedenen Reaktionen zu unterscheiaren den: Zum einen gibt es den koh¨ ten Prozess, b ei dem der Kern in seiner Gesamtheit

7


reagiert, und zum anderen den inkoh¨ ten Prozess, der durch die Reaktion an einem aren einzelnen Proton o der Neutron gekennzeichnet ist. In dieser Arb eit wurde fur b eide, den ¨ koh¨ ten und inkoh¨ ten Prozess, eine Strahl-Polarisations Asymmetrie gemessen, aren aren sowie deren kinematische Abh¨ angigkeiten untersucht. Die zu Grunde liegenden Daten stammen aus dem Jahr 2000 vom HERMES Exp eriment am DESY, wob ei der 27.6 GeV p olarisierte Positronen Strahl von HERA an unp olarisiertem Neongas gestreut wurde.

8


1 Introduction
Most of our knowledge ab out the structure of matter was gained by scattering exp eriments. After Rutherford discovered the existence of a nucleus in an atom more than 90 years ago, it was shown 50 years ago that they are comp osed of protons and neutrons [Hof57]. Over the past few decades, the research on the structure of matter has made enormous strides. In the late sixties, deep inelastic lepton nucleon scattering at SLAC (Stanford Linear Accelerator Center) [A + 76, B+ 83] showed that protons are comp osed of p oint-like ob jects called partons. In the quark-parton-mo del partons are identified with quarks carrying spin and a charge of multiples of one third of the elementary charge e. Further exp eriments showed that the nucleon consists not only of quarks but also of gluons, which are the mediating particles of the strong interactions. At HERA (Hadron Elektron Ring Anlage) of the DESY lab oratory in Hamburg, this research is continued to determine the parton structure of nucleons. In order to investigate in particular the spin structure of the nucleon, the HERMES exp eriment was built. The nucleon spin can b e decomp osed conceptually into the angular momentum contributions of its constituents according to the equation s
N z

=

1 1 = + Lq + J 2 2

g

where the three terms give the contributions to the nucleon spin from the quark spins, the quark orbital angular momentum, and the total angular momentum of the gluons, resp ectively. Since the results from the Europ ean Muon Collab oration (EMC) [Muon88, Muon89] led to the conclusion that 0.1 - 0.2, one of the most interesting questions is to understand how the remaining contribution is distributed b etween the orbital angular momentum of the quarks L q and the total angular momentum of the gluons Jg . The only known way nowadays to access the size of those contributions is to determine the total angular momentum of the quarks in the framework of generalized parton distributions (GPDs). This theoretical formalism describ es the quark-gluon structure of a nucleon and takes into account dynamical correlations b etween partons of different momenta in the nucleon. The ordinary parton distribution functions (PDFs) and form factors turn out to b e the limiting cases and moments of GPDs, resp ectively. The theoretically cleanest pro cess sensitive to GPDs and accessible nowadays is Deeply Virtual Compton Scattering (DVCS). DVCS denotes exclusive pro duction of high energy photons, where in contrast to Bremsstrahlung the photon is not emitted by the lepton but by one of the quarks inside the nucleon. In

9


the present analysis DVCS is studied through p ositron scattering on a nuclear target. This reaction receives contributions from b oth the DVCS pro cess, whose origin lies in the strong interaction, and the electromagnetic Bethe-Heitler (BH) pro cess, well known as Bremsstrahlung. Note that the BH cross section can b e precisely calculated in quantum electro dynamics (QED) using elastic form factors. At HERMES kinematics, the DVCS pro cess can not b e clearly separated from the interfering BH pro cess b ecause the final states are indistinguishable and hence their scattering amplitudes add coherently: | |2 = |
BH

|2 + |

DV C S

|2 +

B H D V C S

+

D V C S B H

.

I nter f er ence ter m I

The BH-DVCS interference term I offers the p ossibility to directly access the DVCS amplitudes. Thereby the imaginary part of the interference term can b e isolated by measuring the angular dep endence of the pro duced photon if p olarized lepton b eams are available. In order to access the real part in addition, lepton b eams of b oth charges are needed. A determination of the angular dep endences in b eam-spin and b eam-charge asymmetry provides the p ossibility to access DVCS amplitudes and through them certain combinations of GPDs. The measurement of DVCS on the proton has shown the p ossibility to provide a sensitive test of current mo dels of GPDs [Ell, Kra]. Such studies op en also access for questions to the prop erties of the quark and gluon matter inside nuclei. Consequently, nuclear DVCS allows the study of the mo difications of particle correlations enco ded in GPDs within the nuclear environment. Hence it may shed new light on the dynamical interplay of highly complex b ound hadronic systems [Pol03]. The aim of this thesis was to study and to extract b eam-spin asymmetries for an unp olarized neon target. The outline is as follows. The theoretical description of GPDs and their relation to the DVCS pro cess is explained in the next chapter. The third chapter will explain the relevant asp ects of the HERMES exp eriment for this study. In chapter 4 the selection criteria for the data and the extraction of the b eam-spin asymmetry is describ ed. The systematic studies are given in chapter 5, and the results are shown in chapter 6. In the end an interpretation of the results by a comparison to hydrogen results is p erformed in chapter 7. The last chapter of this thesis will give a short summary and outlo ok.

10


2 Generalized Parton Distributions and Deeply Virtual Compton Scattering
In the theory of Quantum Chromo dynamics (QCD), one of the main questions is the understanding of the substructure of hadrons. Therefore, most high energy exp eriments study inclusive reactions, such as deep inelastic scattering (DIS): e + p - e + X . These pro cesses are describ ed in QCD by parton distribution functions (PDFs), which enco de the one-dimensional distribution of longitudinal momentum and p olarization, carried by quarks, antiquarks and gluons. By definition PDFs do not contain information on the three-dimensional distribution of the hadronic substructure. In order to extract the additional information, particle correlation functions have to b e measured that dep end on additional variables, characterizing e.g. the momentum difference of the "in" and "out" state. Therefore, pro cesses are required in which the nucleon stays intact, otherwise the dynamical relationship b etween the different partons could not b e accessed anymore. In such pro cesses, either all pro duced particles have to b e detected and clearly separated from the intact final state nucleon, or the missing mass has to b e calculated in order to characterize the difference b etween the initial energy and the sum of the energy of all reconstructed particles in the detector. In particular Deeply Virtual Compton Scattering (DVCS) e + p e + p + provides such a pro cess, which constitutes a hard exclusive reaction. Note that hard means here a large momentum transfer from prob e to target. These exclusive hard reactions, e.g. DVCS, are describ ed by particle correlation functions, called Generalized Parton Distributions (GPDs). GPDs have b een extensively studied theoretically in recent years [M + 94, Ji97a, Rad97, Die03]. They came into the fo cus of interest, after Ji has found that GPDs offer the p ossibility to determine the total angular momentum carried by quarks and gluons [Ji97b]. Nowadays, GPDs start to play a central role in the studies of nuclear structure. Figure 2.1 presents a schematic overview ab out the exclusive pro cesses that are related to GPDs. In comparison to other exclusive pro cesses, DVCS is the theoretically cleanest pro cess accessible to day that allows access to GPDs in the nucleon. This can b e explained by the fact that the real photon in the final state is a p ointlike particle and not a b ound state like a meson, or an even more complicated state. Note that the non-p ointlike part of the real photon wave function is suppressed in DVCS [Rad97, JF99]. Exp erimentally very difficult to access is the pro cess of timelike Compton scattering describing the pro duction of a virtual photon that converts into a lepton-antilepton pair, as it is suppressed by

11


Figure 2.1: Reactions and nucleon prop erties related to GPDs.

2 compared to the DVCS pro cess. In the Wide-Angle Compton Scattering pro cess, em the initial and final photons are real and the involved momentum transfer is large.

2.1 Prop erties of Generalized Parton Distributions
GPDs are universal non-p erturbative ob jects, entering the description of hard exclusive electropro duction pro cesses [GPRV04]. First implicitly intro duced in [M + 94, Ji97a, Rad97], they are generally defined for each quark flavor (u,d,s) and gluon (g) as a parameterization of matrix elements of lightcone op erators [Rad96]. The matrix elements can b e describ ed as elements of the transition matrix b etween the initial and final hadron state. The GPDs can b e characterized by the following features: · They dep end on two longitudinal momentum fraction variables (x, ) and on Mandelstam variable t, defining the momentum transfer to the nucleon, which contains transverse comp onents. The involvement of GPDs in the DVCS pro and the relevant kinematic variables are shown in figure 2.2 which illustrates the also cess the

12


'

$

&

Figure 2.2: DVCS handbag diagram.

%

leading order pro cess in the generalized Bjorken limit. Note that GPDs are defined at a starting scale µ2 and their Q2 evolution is generated by p ertubative QCD with Q2 b eing the photon virtuality. In Off-forward parton-distributions (OFPDs) [Ji98] the variables x, and t are describ ed in longitudinal lightcone co ordinates (or the infinite momentum frame) by:

k



= xP + -+ = 2P + = P -P
2

+

x [-1, 1]

(2.1) (2.2) (2.3) (2.4)

= -t

with x (the average longitudinal momentum fraction of the involved parton with the momentum k ) and the skewedness variable (the longitudinal fraction of the momentum transfer ). The meaning of the longitudinal momentum fraction variables in the context of OFPD is that the parton with the longitudinal momentum fraction x + is removed from the proton and then put back with a longitudinal fraction x - . Note that these longitudinal momentum fractions are either p ositive or negative if the active parton is a particle or antiparticle. In the Bjorken limit of large photon virtuality Q2 and collision energy, the parameter is related to x B j , the momentum fraction of the proton carried by the struck quark, as xB j 2 - xB
j

13


- x

-- x

x +

- x

x +

x -

-1

-

0



1 x

Figure 2.3: The parton interpretation of GPDs in the three x-intervals [-1,- ], [- , ], and [ ,1] [Die03 ]. · GPDs are defined in the interval x [-1, 1], which is sub divided into three regions, shown in figure 2.3: 1. for x [ , 1]: b oth momentum fractions x + and x - are p ositive, describing emission and reabsorption of a quark. 2. for x [- , ]: x + 0 and x - 0, interpreted as an antiquark with momentum fraction - x emitted from the initial proton. 3. for x [-1, - ]: b oth momentum fractions x + and x - are negative, b elonging to emission and reabsorption of antiquarks

The first and third case are commonly referred to as DGLAP (Dokshitzer-Grib ovLipatov-Altarelli-Parisi) regions and the second as ERBL (Efremov-RadyushkinBro dsky-Lepage) region [Die03]. · In the case of the spin-1/2 nucleon the most imp ortant GPDs are the chirally-even ~ ~ GPDs H q,g ,H q,g ,E q,g ,E q,g , defined for each quark flavor (q = u,d,s) and gluon (g). Chirally-even means here that they do not flip the parton helicity [HJ98]. They can b e divided into unp olarized GPDs (H q,g and E q,g ) and p olarized ones (H q,g and E q,g ). Of these, H q,g and H q,g conserve the nucleon helicity, while E q,g and E q,g can also flip it. The numb er of GPDs for spin-zero hadrons (pions and spin-zero nuclei, e.g. neon) is reduced to one GPD H q,g for each quark flavor and gluon [Die03]. · The first moments of the helicity-conserving GPDs are related to the elastic form factors of the nucleon [Ji97b]:
1 -1 1 -1 1 -1 1 -1

dxH q (x, , t) = F1 (t) dxE q (x, , t) = F2 (t) dxH q (x, , t) = GA (t) dxE q (x, , t) = GP (t)

Dirac form factor Pauli form factor Axial-vector form factor Pseudo-scalar form factor

Since the result of the integration do es not dep end on , one can cho ose = 0. In

14


the limit of vanishing momentum transfer t - 0, the helicity conserving GPDs H q (x, = 0, t) and H q (x, = 0, t) reduce to H q (x, = 0, t = 0) q (x) and H q (x, = 0, t = 0) q (x). For gluons they reduce to H g (x, = 0, t = 0) = xg (x) and H g (x, = 0, t = 0) = xg (x). The reduced distributions are the ordinary spin-indep endent density q (x) and spin-dep endent density q (x) for quarks and xg (x), xg (x) for gluons. The corresp onding relation for the distributions E q,g and E q,g , which would express the "magnetic" density, can not b e directly expressed in terms of any known parton distribution [GPRV04]. · The t-indep endent part of GPDs complies the p olynomiality condition following from Lorentz invariance [Ji97b]
1

dx xN H q (x, ) = h
-1 1

q (N ) 0

+h

q (N ) 2 2

+ ... + h

q (N ) N +1 N +1

,

(2.5)

dx xN E q (x, ) = e
-1

q (N ) 0

+e

q (N ) 2 2

+ ... + e

q (N ) N +1 N +1

,

(2.6)

where the coresp onding p olynomials contain only even p owers of the parameter due to time reversal invariance [MPW98, Ji98]. In order to satisfy the p olynomiality condition, a parameterization of GPDs can b e defined by the Double Distribution (DD) formalism, which is completed by adding a D-term intro duced by Polyakov and Weiss [PW99]. The present exp erimental knowledge on GPDs is presented in figure 2.4. The GPDs are placed in the middle of three concentric rings, the PDFs in the outermost and nucleon form factors in the innermost rings. To day's exp erimental knowledge of the different functions is illustrated in different colors from light (no data exist) to dark (well known). Note that for empty sectors no function exists or no strategy is known for its measurement. For more details see Ref. [Now02].

2.2 The Energy-Momentum Tensor and the Spin of the Nucleon
In order to understand the spin structure of the nucleon, the b est way is to start with the QCD angular momentum op erator in its gauge-invariant form . The sum of the quark and gluon contribution is given as [Ji03]: J where Jq = d3 x x â T
q QC D

= Jq + J

g

15


Figure 2.4: Visualization of the most relevant Generalized Parton Distributions and their limiting cases, forward Parton Distributions and Nucleon Form Factors [Now02].

and Jg = d3 x x â (E â B ).

Here Tq and (E â B ) are the quark and gluon momentum densities, resp ectively. The separation of the quark and gluon contributions to the nucleon spin can b e deduced from an analogy with the magnetic moment, if the form factors of the momentum density are

16


known at zero momentum transfer t: µ = GM (t = 0) = F1 (t = 0) + F2 (t = 0), where the spin-flip form factor GM (t) = q (F1 (t) + F2 (t)) yields the electric current distribution in the nucleon with the Dirac and Pauli form factors of the electromagentic current F1 (Q2 ) and F2 (Q2 ). In analogy, the second moments of the parton distributions yield the form factors of the energy-momentum tensor,
1

dx(H
-1

q ,g

(x, , t) + E

q ,g

(x, , t))x = A(t) + B (t),

(2.7)

where the dep endence drops out. Extrap olating this relation to t 0, the total quark and gluon contribution to the nucleon spin is obtained: J with
q ,g

1 = (A 2

q ,g

(0) + B

q ,g

(0))

(2.8)

1 = Jq + Jg . (2.9) 2 1 The relation Lq = Jq - 2 offers the p ossibility to determine L q through a measurement of Jq by using the information on the quarks spin contribution 1 , which is available 2 from inclusive and semi-inclusive p olarized DIS. In summary, the second moment of GPDs yield the total angular momentum
1

J

q ,g

= lim

t0 -1

x[H

q ,g

(x, , t) + E

q ,g

(x, , t)]dx,

(2.10)

well know as Ji's relation.

2.3 Kinematics of Deeply Virtual Compton Scattering
Deeply Virtual Compton Scattering (DVCS) is the hard exclusive pro duction of a real photon in lepton scattering e(k ) + P (p) e(k ) + P (p ) + (q ), with the four-momentum of the incoming (outgoing) lepton k (k ), the initial (final) hadron p(p ) and the real photon q . The exchanged virtual photon is describ ed by the four-momentum q = k - k . In analogy to ordinary Compton scattering, where a real photon is scattered on a charged particle, the DVCS pro cess refers to virtual -photon scattering with a real photon b eing pro duced.

17


Figure 2.5: Diagram for Deeply Virtual Compton Scattering with the azimuthal angle b etween lepton scattering plane and photon pro duction plane and the p olar angle b etween virtual and real photon.

The reaction of DVCS is defined in the generalized Bjorken limit of large photon virtuality and collision energy. DVCS is characterized through an reaction of the virtual photon with partons, where a real photon in the final state is pro duced, as illustrated in figure 2.5. The azimuthal angle is defined as angle b etween lepton scattering plane and photon pro duction plane, given as -- - - - - - q · q âk ( â k )( â q ) q q = arccos (2.11) -- - -- - - | â k || q â q | | q · q â k | q The angle b etween virtual and real photon is


, defined as . (2.12)

= arccos

- q - | q

- ·q - || q |

The following kinematic variables are defined: · The photon virtuality,
lab Q2 -q 2 = -(k - k )2 4E E sin2 ( ) = 2

(2.13)

with the lepton energies E and E . · The momentum fraction of the proton carried by the struck quark (Bjorken variable), Q2 Q2 = (2.14) xB j 2pq 2Mn with the mass of the nucleon Mn .

18


'

$

Figure 2.6: Diagrams for the DVCS pro cess (a) and the BH pro cess (b). In the latter a photon can b e radiated by the incoming or outgoing lepton.

&

%

· The energy of the virtual photon, pq lab = E-E . Mn (2.15)

· The fractional energy of the virtual photon to the b eam energy, y pq lab =. pk E (2.16)

· The invariant mass of the system of virtual photon and proton (squared energy in the photon-proton center-of-mass system),
2 W 2 = (q + p)2 = Mn + 2Mn - Q2 .

(2.17)

· The four-momentum transfer to the nucleon, t = (p - p )2 = (q - q )2 . (2.18)

2.4 Beam-Spin and Beam-Charge Asymmetries
In fact, not only the DVCS pro cess pro duces a real photon in the final state. Photons are also emitted by leptons through Bremsstrahlung, characterized by the so-called Bethe Heitler (BH) pro cess. In figure 2.6 b oth pro cesses are shown. The DVCS pro cess can not b e separated from the BH pro cess, b ecause the final states are indistinguishable; the scattering amplitudes of the two interfering pro cesses add coherently.

19


The square of the photon pro duction amplitude receives contributions from pure DVCS (DV C S ) and pure BH (B H ) and from their interference | |2 = |
BH

|2 + | I =

DV C S

|2 +

B H D V C S

+

D V C S B H

(2.19)

with the DVCS-BH interference term
B H D V C S

+

D V C S B H

.

(2.20)

The corresp onding five-fold differential cross section for a spin-0 target (e.g. neon) p er nucleon is given as [FS04] 1 A dx d = dy d|t|dd 16 2 Q
2 3x em B j

y
2 Bj
2 MN Q2

Bj

1 + 4x

A

3

|

| e3

2

(2.21)

with the fine-structure constant em , the charge of the lepton e, the atomic numb er A of the nucleus with the mass MN and the angle b etween the target p olarization vector and the hadron if the hadron is transversely p olarized. Note that the nuclear amplitudes are given p er nucleon, and in order to obtain the correct scaling of the cross section with A, the DVCS amplitude has to b e multiplied by A. In the case that the BH cross section is not dominating the DVCS cross section, it is p ossible to obtain the DVCS cross section through a measurement of the full cross section subtracting the BH cross section. Note that the BH cross section is exactly calculable in Quantum Electro dynamics (QED). The measurement of the DVCS cross section has b een done in the collider exp eriments H1 [H101] and Zeus [HZ] at DESY. Since HERMES is a fixed-target exp eriment, its kinematic region is disjunct from that of the collider exp eriments which measure at small x. Here the BH cross section mostly dominates the DVCS cross section and the subtraction of the calculated BH cross section might lead to large uncertainties. In contrast, the BH-DVCS interference term I offers the p ossibility to directly access the DVCS amplitudes. The three parts of the photon pro duction cross section 2.19, expanded in Fourier series, are for an unp olarized target [BMK02] | |
BH

| |

2

=- =

e x
2 y2 Bj

6

2

(1 +

2 4MN x Q2

2 Bj

c ) tP1 ()P2 () c
DV C S n 3 2

BH 0

+
n=1

c

BH n

cos n , (2.22)

DV C S

2

e6 y2 Q x

2 2

c

DV C S 0

+
n=1

cos n + s
I n

DV C S 1 2

sin , sI sin n , n

(2.23) (2.24)

I=

Bj

±e6 y 3 tP1 ()P2 ()

c+
n=1

I 0

c cos n +
n=1

where the +/- sign in the interference term I stands for a negatively/p ositively charged incident lepton. In the squared BH term | B H |2 and in the interference term I , an additional dep endence arises from the scaled lepton BH propagators P i () = Ai + Bi cos .

20


The cn and sn are the Fourier co efficients and dep end on the variables y , x an unp olarized target the helicity-dep endent co efficients are given by
I 1 DV C S 1 sI 2

Bj

, t, Q2 . For

s

s

= 8K y (2 - y ) m C = -8K/(2 - x
Bj

I DV C S I

)y

mC

= -16K 2 /(2 - x

Bj

)y

mC

with K (1 - xB j )(1 is the lepton helicity. (CFFs) F = {H, H, E , E order in s , the function

- y ) , where is the transverse momentum transfer, and Q The C 's are linear combinations of the Compton form factors }, which have b een worked out in Ref. [BMK02]. At leading C I for an unp olarized target is given by xB j 2 (F1 + F2 )H - F2 E 2 - xB j 4M 2

C I = F1 H +

where F1 and F2 are the nucleon Dirac and Pauli form factors, resp ectively. In order to get access to the helicity-dep endent Fourier co efficients, one has to measure the b eam-spin asymmetry, defined as [BMK02] A
LU

() =

d +d

d -d








.

(2.25)

This cross section asymmetry is built from two measurements with opp osite b eam helicity, denoted by arrows with opp osite orientation. The indices L and U of the asymmetry A denote the longitudinal p olarized b eam and the unp olarized target. By inserting the Fourier co efficients [BMK02] in the cross sections of equations 2.22 2.24, the b eam-spin-asymmetry ALU can b e expressed as d -d
e6 y2 Q
2

d +d



= |

s

DV C S 1

sin ±
DV C S 0

x

Bj

e6 y 3 tP1 ()P2 ()

(sI sin + sI sin 2) 1 2
e6 (cI + 0
3 P n=1

.

BH

|2 +

e6 y2 Q

2

2

(c

+
n=1

c

DV C S n

cI cos n) n

cos n) ±

x

Bj

y 3 tP1 ()P2 ()

To leading order in s and neglecting DVCS and interference part in the denominator it reduces to: xB j sI 1 sin (2.26) ALU () ± y cB H 0 where the +(-) sign stands for a negative (p ositive) charged lepton b eam. Note that the BH-co efficient cB H is exactly calculable in quantum-electro dynamics. In this ap0 proximation, the asymmetry ALU is prop ortional to sin and sI , the Fourier-co efficient 1 of the interference part, which is given as s I = 8K y (2 - y ) m C I (F ) . With 1 m {F } =
q

e2 (F q ( , , t) - F q (- , , t)) q

(2.27)

21


the imaginary part directly prob es the resp ective GPDs F q = {H, H , E , E } along the line x = ± . Another p ossibility to get access to the Fourier co efficients is provided by the b eamcharge asymmetry [BMK02] d + - d - (2.28) AC () = d + + d - where d + and d - denote the cross section measurement with a p ositron and an electron b eam, resp ectively. This cross section difference directly pro jects out the interference term I , since it is the only term dep ending on the charge, indicated by the sign ± in equation 2.24. Assuming unp olarized b eam ( = 0) and to leading p ower 1/Q, the b eam-charge asymmetry at leading order in s is approximately given as AC () = - xB j cI 1 cos y cB H 0 (2.29)

gives access to the integral over the resp ective GPDs, where P denotes the Cauchy's principal value. Since neon data at HERMES is presently only available for one lepton b eam charge (p ositrons), the b eam-charge asymmetry will not b e further discussed in this thesis.

with cI = -8K (2 - 2y + y 2 ) e C I (F ) . Consequently, the b eam-charge asymmetry 1 for an unp olarized b eam has a cos dep endence and is prop ortional to the real part of C I (F ). The real part of the Compton form factor 1 1 1 , (2.30) ± e {F } = - e2 P dx F q (x, , t) q x- x+ q
-1

2.5 Deeply Virtual Compton Scattering on Nuclei
The measurement of DVCS on the proton has shown the p ossibility to provide a sensitive test of current mo dels of GPDs [Ell, Kra]. Such studies op en also access to questions of the prop erties of quark and gluon matter inside nuclei. Consequently, nuclear DVCS allows the study of the mo dification of parton correlations enco ded in GPDs within a nuclear environment. Hence it may shed new light on the dynamical interplay of highly complex b ound hadronic systems. In the case of the simplest nucleus, the deuteron, which can b e seen as a weakly b ound state of a proton and a neutron, first results have already b een published [HERMES03]. In comparison to the deuteron, describ ed by two-b o dy wave functions p+n and their transition GPDs, it is even more complicated to work out a theoretical prediction for

22


heavier nuclei, as the numb er of wave functions increases with the atomic numb er of the hadron. The numb er of GPDs also dep ends on the spin of the hadron: for spin-1 nuclei (i.e. Deuteron) 9 GPDs, for spin- 1 nuclei 4 GPDs and for spin-0 nuclei (i.e. Neon) only one 2 GPD for each parton. In the case of Neon it reduces to the GPD H q,g [Die03]. For a nuclear target there exist two distinct pro cesses, b oth in DVCS and BH: 1. The scattering pro ceeds coherently, i.e. the target nucleus recoils as a whole while emitting a photon with momentum q , as illustrated in the left panels of figure 2.7 [LT05a]. 2. The scattering pro ceeds incoherently, i.e. the nucleus undergo es a breakup and the final photon is emitted from a quasi-elastically scattered nucleon, as shown in the right panels of figure 2.7. The extension to the case of heavier spin-0 nuclei has b een considered, among others, by Guzey and Strikman [GS03] and is explained in the following. In order to obtain an estimate of the change of ALU , the ratio of nuclear to proton asymmetries b een calculated in the limit t = 0:
Anucleus LU AP r oton LU

has

A H p ( , , 0) - H p (- , , 0) + ( Z - 1) (H n ( , , 0) - H n (- , , 0)) Anucleus() LU = . (2.31) H p ( , , 0) - H p (- , , 0) AP roton () LU

Note it was assumed that the pro cess is a coherent reaction. Since b oth coherent and incoherent contributions enter the total cross section, the exnucl pression for ALU should b e mo dified. The mo dified asymmetry A LU eus b ecomes [GS03]: A
nucleus LU

() = sin()

K 8y (2 - y )x y cB H 0

Bj

e.m. (Z (H p ( , , t) - H p (- , , t)) + Z (A - 1)FA (t ) (H p ( , , t ) - H p (- , , t ))) e.m. Z F1 (t) + Z (Z - 1)(FA (t ))2 (2.32) with t = t AA 1 . Figure 2.8 shows a schematic representation of the origin of the com- binatorial factors Z and Z (A - 1). The first term in nominator and denominator of equation 2.32 describ es the contribution from the "in" and "out" states to the same nucleon (incoherent term), which at small t is prop ortional to the numb er of protons, Z , times the GPD H of the free proton. This contribution has a slow t-dep endence due to the proton elastic form factor F1 (t). Note that the neutron contribution is suppressed by the smallness of the electromagnetic form factors and can therefore b e neglected. The contribution given by the second term in the nominator and denominator describ es the

â

23


q

q=q+

q

q=q+ P=P- P

P

P=P- PA-1

PA

P =PA- A

PA

(a)


(b)

P=P- P P P=P- P
A-1

PA

PA=PA-

PA

(c)

(d)

Figure 2.7: Diagrams for DVCS and BH pro cesses from a nuclear target at leading order in s . (a) DVCS, coherent pro cess; (b) DVCS, incoherent pro cess; (c) BH, coherent pro cess; (d) BH, incoherent pro cess [LT05a]. Z (A - 1) attachments of BH to the proton and DVCS to a different nucleon. It is mostly coherent and has a much steep er t-dep endence, forced by the nuclear charge form factor e.m. FA (t). Since the DVCS measurement do es not allow to extract only purely coherent DVCS events, the measured asymmetries present a sum of the coherent and incoherent contributions, as given by equation 2.32. A simulation of the ratio of asymmetries in the cases of neon (A = 20 and Z = 10) and krypton (A = 76 and Z = 36) is presented in figure 2.9. While
Anucleus LU Apr oton LU

is significantly larger than unity for coherent nuclear DVCS
Anucleus LU Apr oton LU

Consequently, the inclusion of the incoherent contribution should decrease the ratio of

(exp ected to b e close to the factor of 2 for t 0),

1 for the incoherent part.

24


J

T

J

T

A a)

A

A b)

A

Figure 2.8: Schematic representation of the interference b etween the BH (J ) and DVCS (T ) amplitudes in nuclei: There are Z attachments of b oth J and T to the same proton (a) and Z (A - 1) attachments of J to the proton and T to a different nucleon (b) [GS03]. the asymmetries. A quantitative evaluation of the coherent and incoherent contributions to the ratio of the nuclear BSA over the proton BSA has b een worked out in a different mo del [LT05a, LT05b]. In this mo del off-shell effects are quite noticeable and an enhancement nucl of the ratio ALU ear /Aproton for higher t is predicted. In figure 2.10 the ratio of nuclearLU to-proton asymmetries AA /Ap for 4 He is shown, calculated including only coherent LU LU scattering terms in b oth the DVCS and BH contributions to the asymmetry (dashed line), and including only the incoherent terms (dot-dashed line). In these calculations the nuclear mo del included off-shell effects. An additional dep endence on the atomic numb er A has b een predicted in a study [Pol03] using a mo del for a large nucleus. A characteristic A-dep endence of the D-term was obtained in connection with the spatial comp onents of the energy-momentum tensor. In fact, a conclusion of this study implies, that the contribution of the D-term to the real part of the DVCS amplitude grows with an increase of the atomic numb er as A 4/3 . It is also mentioned as a general p ossibility that detailed information ab out deviations of energy, pressure, and shear forces distributions inside nuclei can b e accessed.

25


Figure 2.9: The ratio of nuclear to proton asymmetries A nuclear /Aproton for Neon (thick) LU LU and Krypton (thin), solid curves for b oth coherent and incoherent, dashed curves include only the coherent part [GS03].

26


A(A)/A(p)

LU

LU

1.3

1.25 1.2 1.15 1.1 1.05 1 0.95 0.9 0.85 0.8 0 coherent only incoherent only 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5

t (GeV2)
Figure 2.10: The ratio of nuclear to proton asymmetries A A /Ap for 4 He, calculated LU LU including only coherent scattering terms in b oth the DVCS and BH contributions to the asymmetry (dashed line), and including only the incoherent terms (dot-dashed line). The nuclear mo del including off-shell effects was used in these calculations [LT05a].

27


3 The HERMES Experiment
The aim of the HERMES (HERa MEasurement of Spin) exp eriment is the investigation of the spin structure of the nucleon. The HERMES exp eriment is lo cated in the East Hall of the HERA (Hadron Elektron Ring Anlage) storage ring complex at DESY, using the p olarized lepton b eam for scattering off a fixed gas target. The physics program for HERMES is very broad. The exp eriment was inspired by the flavour decomp osition of the nucleon spin. In addition it contributes inclusive data with qualitatively different systematic uncertainties to improve the world data set for the x dep endence and the integral of the spin structure function g 1 (x). A novel asp ect of the exp eriment is its capability to p erform measurements of exclusive pro cesses, even though it was originally conceived as an exp eriment of inclusive and semi-inclusive physics.

3.1 The HERA Storage Ring Facility
The underground storage ring facility HERA has a circumference of 6.3 km and consists of two b eam lines for the 27.57 GeV/c electron (p ositron) and the 920 GeV/c proton b eam resp ectively. Both b eams can b e used simultaneously by four different exp eriments. Two of these exp eriments, ZEUS and H1, are colliding-b eam exp eriments. HERA-B uses only the proton b eam, whereas HERMES accumulates data only with the p olarized lepton b eam. After the start of HERMES in 1995 the lepton ring was op erating with p ositrons, except for a short p erio d in 1998. As recently as the b eginning of 2005, data taking with electrons has again started. The lepton b eam consists of 189 individual bunches, which have a 96 ns (29 m) ration. Each lepton bunch is 27 ps (8 mm) long and contains 2 · 10 10 particles. average lifetime of the lepton b eam is around 10 hours with a starting intensity mA which decreases to 10 mA. After that, the b eam is usually dump ed or a one sp ecial HERMES run with high gas density follows. sepaThe of 50 hour

The lepton b eam is initially unp olarized after injection, but due to the Sokolov-Ternov effect [ST64], a parallel p olarization to the magnetic dip ole fields is naturally build up. That means the b eam is p olarized in the transverse direction, i.e. the b eam spin orientation is p erp endicular to the momentum. The Sokolov-Ternov effect can b e describ ed as a small asymmetry in the spin-flip probability in the emission of synchroton radiation. Hence the b eam p olarization increases exp onentially in time, with a rise-time of ab out 25

28


Figure 3.1: Schematic view of the HERA collider

minutes. In fact, dep olarization effects limit the achievable p olarization. The main contribution arises from non-vertical magnetic field comp onents with resp ect to the b eam orbit. These non-vertical magnetic fields are caused by small magnet misalignments and not p erfectly homogenous magnetic fields. Furthermore, the interaction of the lepton and proton b eam in the interaction regions of HERA may reduce the b eam p olarization. In order to reduce the influence of dep olarization effects, the so-called harmonic-bump metho d is used [B+ 94]. This scheme intro duces additional vertical closed-orbit corrections at strategic lo cations to comp ensate the effect of the spin-orbit distortions. At HERA there are eight harmonic bumps available and help ed to achieve up to 60% p olarization in the year 2000. Note that the p olarization is defined as N - N N + N


whereby N (N ) refers to the numb er of leptons with their spin aligned parallel (antiparallel) to the magnetic dip ole field. For the study of the helicity structure of the nucleon, a longitudinally p olarized lepton b eam is required, i.e., the lepton spins have to b e aligned parallel to their propagation direction. Therefore, b efore and after the interaction region, a pair of 90 o spin rotators have b een installed, to turn the spin into the b eam direction and then back into the vertical again. In figure 3.1 the spin rotators can b e seen, one upstream and the other one downstream of the HERMES exp eriment. They consist of a sequence of horizontal and vertical b eam normal conductor magnets. The b eam p olarization is measured by two indep endent p olarimeters, one for the transverse p olarization (TPOL) and the other one for measuring the longitudinal b eam p olarization (LPOL). For the determination of the p olarization, b oth use asymmetries in

29


the Compton scattering of p olarized laser light off the lepton b eam. The transverse p olarimeter [B+ 93, B+ 94], measures the p olarization of the lepton b eam at a p oint where it is p olarized in transverse direction. The p olarization measurement is based on a spatial asymmetry effect in the back-scattering of laser light off the p olarized lepton b eam. These back-scattered photons are measured in a split lead-scintillator sampling calorimeter, where the change in the p osition of the photons with initial circular p olarization determines the p olarization of the lepton b eam. The second p olarimeter, the LPOL [B + 02], measures the longitudinal p olarization of the b eam when crossing the exp eriment. It is also based on Compton back-scattering of laser light, but instead of the spatial distribution the asymmetry in the total crosssection is used. The larger asymmetry in this case allows a more precise measurement of the p olarization. In addition, the p olarization is measured for each individual lepton bunch in HERA. Furthermore the existence of two p olarimeters provides the p ossibility for a cross-check of the p olarization measurement.

3.2 The Internal Gas Target
Since a storage ring can not b e op erated with a liquid or solid target, the use of a gas target was mandatory for the HERMES exp eriment. The target system comprises the Atomic Beam Source (ABS), the internal storage cell and the target monitor system with the Target Gas Analyzer (TGA) and the Breit Rabi Polarimeter (BRP). The ABS-system consists of a disso ciator to form atoms, a skimmer and a collimator for the formation of the atomic b eam, a sextuple magnet system to sample out unwanted electron p olarization, and high frequency units to transfer the p olarization of electrons to nucleons. In figure 3.2 a schematic view of the HERMES ABS-system is presented. It is also p ossible to inject unp olarized gases into the storage cell by the unp olarized gas feed system (UGFS). The density of unp olarized gas can b e up to two orders of magnitude larger in comparison to the density of p olarized gas, which is at the maximum achievable value for present technology. Since the start of HERMES, a variety of pure unp olarized gas targets was used like H, D, He, N, Ne, Kr, Xe. The target region is shown schematically in figure 3.3. The gas enters an op en-ended tub e that confines the gas atoms in a region around the p ositron b eam. The storage cell increases the areal target density by ab out two orders of magnitude compared to a free atomic b eam. The cryogenically co oled storage cell is an elliptical tub e, 9.8 mm high by 29 mm wide and 400 mm in length made of ultra-pure aluminum with a uniform wall thickness of 75 µm [HERMES98]. At the op en ends of the target cell the gas atoms, which leak out, are pump ed away. Particles, emerging from the interaction, which are scattered into the sp ectrometer acceptance, exit the target chamb er through a thin (0.3 mm) stainless steel foil.

30


r (mm)

60 40 20 0 -20 -40

discharge tube

1st sexp. magn. syst. coll.

SFT

MFT

2nd sexp. magn. syst.

SFT / WFT

nozzle

injection tube

cell

-1400

-1200

-1000

-800

-600

-400

-200

0 z (mm)

Figure 3.2: Schematic view (seen against the electron b eam) of the HERMES ABS with disso ciator and collimator for b eam formation. Two sets of sextuple magnets are lo cated along the axis of the gas jet. The axis of the ABS is tilted by 30 downwards with resp ect to the horizontal plane.

Figure 3.3: Schematic view of the target cell.

In order to decrease the amount of synchroton radiation, two collimator are installed near the target. In addition, the synchroton radiation is reduced by two weak dip oles downstream of the last b ending magnet in the arc of the accelerator.

3.3 The HERMES Sp ectrometer
The HERMES sp ectrometer [HERMES98] is an op en-ap erture forward sp ectrometer, consisting of two symmetric halves ab ove and b elow a central shielding plate in the b eam plane. A schematic view of the sp ectrometer is shown in figure 3.4. The co ordinate

31


Figure 3.4: Schematic side view of the HERMES exp eriment. The incoming lepton b eam enters from the left side of the picture.

system of HERMES is defined by the z axis p ointing along the b eam direction, the y axis oriented vertical upwards, and the x axis horizontal, p ointing towards the outside of the ring. The p olar and azimuthal scattering angles, as well as the initial tra jectory for the determination of the momentum of the particle, are measured by the front tracking system, which consists of the Vertex Chamb er (VC) and the Drift Chamb ers (DVC, FCs). Two sets of drift chamb ers b ehind the magnet (BCs) are used for the momentum measurement of leptons. In addition, three chamb ers in the magnet (MCs) help to match front and back tracks as well as to track low momentum particles. The particle identification (PID) is provided by a combination of signals from the electromagnetic calorimeter, the preshower detector (H2), the Transition Radiation Detector (TRD), a threshold Cerenkov detector (C) [b efore 1998] or a Ring-Imaging Cerenkov counter (RICH) [after 1998]. Since the electromagnetic calorimeter is the only detector able to measure the p osition and energy of the photon, it is the most imp ortant detector for this analysis. The acceptance is limited at small angles by an iron plate in the b eam plane, which shields the lepton and proton b eams from the magnetic field of the sp ectrometer magnet. Therefore, the range of scattering angles is 40 mrad to 220 mrad. The range for the variable x, which defines the momentum fraction of the proton carried by the struck quark, is 0.02-1.0 in the HERMES kinematics.

32


3.3.1 The Tracking System
The kinematics of charged particles traversing the sp ectrometer are reconstructed by the tracking system, which comprises several detectors, b efore, inside and b ehind the sp ectrometer magnet. The following tasks have to b e p erformed by the tracking system: · Measurement of the scattering angles for kinematic reconstruction. · Determination of the event vertex in the target region. · Measurement of the particle momentum from the track deflection in the sp ectrometer dip ole magnet · Identification of the hits in the PID detectors asso ciated with each track. The tracking system consists of 51 planes of wire chamb ers. The vertex reconstruction and the determination of the scattering angles is provided by the front tracking system, which consists of the drift vertex chamb ers (DVC) and the front drift chamb ers FC1,2. In order to determine the track momentum of the charged particles, a combination of the front tracking information with the signals of the back drift chamb ers BCs is p erformed. These BCs measure the charged tracks b ehind the sp ectrometer magnet. The latter has an integrated field strength of 1.3 Tm and contains three prop ortional chamb ers. These so-called Magnet Chamb ers (MCs) helps resolving multi-tracks and provide the analysis of low energy tracks, which do not reach the BCs. The momentum resolution of the HERMES sp ectrometer ranges from 1.5 to 2.5% [HERMES05]. The uncertainty in the scattering angle is b etter than 0.6 mrad.

3.3.2 The Particle Identification
In order to achieve a high efficiency lepton-hadron separation, the resp onses of various dedicated detectors are combined. The particle identification system consists of the RICH counter, the TRD, two sets of ho doscop es (H1,2), and the electromagnetic calorimeter. At HERMES, particles are identified using a probabilistic algorithm that utilizes the resp onse of these four detectors. The resulting logarithmic ratio of the hadron and lepton probability is commonly called PID. The resp onses of the four PID detectors are combined into probabilities using a Bayesian algorithm, which maximizes the lepton-hadron separation. The PID co de asso ciates each track with the resp onses of each PID detector mo dule along the track. With these resp onses, calculations are p erformed and the result are parameters called PID2 through PID5, which corresp ond to the probability that a particular track was a lepton rather than a hadron. For a track with a particular momentum p and the

33


NDIS

700 600 500 400 300 200 100 0 4 2 0 -2 -4 -6 -10 -8 -6 -4 -2 0 2 4 6 8

PID2

PID5

Figure 3.5: Distribution of DIS events versus the PID2 and PID5 values, based on unp olarized neon data.

resp onses of the PID detectors, denoted by R, the quantity PID is then calculated as P I D = log
10

Pe (p, R) Ph (p, R)

whereby the probability distribution that an electron (hadron) of momentum p caused a resp onse R, is presented by Pe (p, R) (Ph (p, R)). A large p ositive PID value means that the track was very likely an electron, and a large negative value indicates a hadron. For the case PID = 0, the probability for b oth is the same. At HERMES the following combinations of PID values are commonly defined, P I D2 P I D
cal cal tr d

+ P ID + P ID

pr e pr e

(3.1) + P ID
cer

P I D3 P I D

(3.2) (3.3)

P I D5 P I D

The DIS-event distribution in simultaneous dep endence of PID2 and PID5 is presented in figure 3.5. For this analysis the condition P I D 2 + P I D 5 2 is used, indicating that only those tracks are selected which were at least 10 times more likely to b e a p ositron than a hadron. The distribution for the sum of PID2 and PID5 is shown in figure 3.6. In fact, only the preshower-detector (ho doscop e H2) and the calorimeter are able to detect photons and have therefore a particular imp ortance for this analysis. The preshower detector in each sp ectrometer half consists of 42 vertical 1 cm thick scintillator paddles.

34


NDIS

1200 1000 800 600 400 200 0

-10

-5

0

5

10 PID2 + PID5

Figure 3.6: Distribution of DIS events versus the sum PID2+PID5, based on unp olarized neon data.

Each scintillator is optically combined with a photomultiplier at the outside of the detector. In front of the scintillators a lead shield with a thickness of 11 mm causes the ma jority of particles to pro duce electromagnetic showers. The lepton-hadron separation uses the fact that the probability for pro ducing electromagnetic showers is considerably larger for leptons than for hadrons. Hence leptons dep osit more energy in the scintillators and are thus energy distinguishable from hadrons. In contrast to the preshower detector, the electromagnetic calorimeter is capable to provide a measurement of the energy and the p osition of the photon. The calorimeter consists of 840 radiation-resistant lead-glass blo cks, divided in two parts, ab ove and b elow the b eam pip e, as indicated in figure 3.7. Each blo ck is viewed from the rear by a photomultiplier tub e. The blo cks have an area of 9â9 cm 2 and a length of 50 cm. In order to prevent the lead-glass blo cks to suffer from radiation damage, b oth calorimeter walls can b e moved away vertically from the b eam pip e for b eam injection. The hit p osition of the photons has to b e determined from the energy distribution inside a cluster, which refers to a 3â3 array of lead glass blo cks. Note that the energy sum in a cluster is almost indep endent of the hit p osition.

3.3.3 The Trigger System
The Trigger System selects events that are p otentially interesting for physics analysis. At HERMES, different trigger schemes are implemented, which select Deeply-Inelastic Scattering (DIS) events and photopro duction pro cesses (without detection of the scat-

35


Figure 3.7: Schematic view of the preshower detector and the calorimeter.

tered lepton). There exist additional triggers for detector monitoring and calibration. The main physics trigger is trigger-21 for the DIS events, requiring hits in the three scintillator ho doscop es and a signal from the calorimeter corresp onding to a dep osition of at least 1.4 GeV in two adjacent calorimeter rows in coincidence with the HERA bunch signal (HERA clock ). The data used in this analysis was collected with a threshold of 3.5 GeV in the calorimeter, corresp onding to the fact that two thirds of the trigger had tracks, where 95% of reconstructed tracks came from the target, and one third had accompanying leptons.

3.4 Data Acquisition
3.4.1 Data Taking
At HERMES exists three mo des of data taking. The most common mo de is standard p olarized-target data-taking where the luminosity is limited by the attainable target densities. To collect data from unp olarized targets there are two ways. Unp olarizedtarget data can b e collected either by using target densities that are limited by the

36


b eam lifetime, or alternatively in high-density running, which is done regularly at the end of each fill when the collider exp eriments have switched off. Under these conditions the luminosity is limited by the trigger rates or the dead time. For this analysis only unp olarized-target data exists for the neon gas target, from which 97% of the data was accumulated in high density runs.

3.4.2 Data Acquisition and Structure
The backb one of the data acquisition system is constructed in Fastbus. It consists of 10 front-end crates, the event collector crate, and the event receiver crate connected to the online workstation cluster, whereby the CERN Host Interfaces (CHI) act as Fastbus master. Their readout p erformance is enhanced in most places with Struck Fastbus Readout Engines (FRE). HERMES data sets are divided into fills, runs and bursts. A fill is the time p erio d in which the HERA lepton ring accelerates and stores the same fill of leptons until the lepton b eam is dump ed. The usual fill time is 8-14 hours. Data is taken after a trigger o ccurs during the measurement. Then the acquisition system translates the resp onses from all detectors into a digital form and stores it in a file on a hard disk. A run is determined by a file size of approximately 450 MB, corresp onding to 10-30 min data taking. In addition, all data not asso ciated to a single event but to the hardware, target and b eam p erformance is stored in the slow-control data. The slow-control refers to the reading and recording of hardware information that changes on a slow time scale, which is read out every 10 sec. This chosen time p erio d defines a burst. The run-file is written in the Exp erimental Physics Input Output (EPIO) format on a hard disc and copied to the DESY tap e rob ot. After that, several programs convert this raw data into useful physics quantities: · At first the edge of the ture of the are defined HERMES DECo der (HDC) deco des the EPIO files by using the knowldesign of the sub detectors. The output file is stored in a tabular strucAleph DAta MOdel (ADAMO). Before writing a program, the tables using a sp ecific Data Definition Language (DDL).

· The ADAMO database is sent to the HERMES ReConstruction (HRC) program. This program combines individual detector resp onses with alignment and calibration information, in order to reconstruct the tra jectories in the sp ectrometer. · After that, the HRC output and the slow-control data are combined into a single file for every run. The output is written as a micro-Data Summary Tap e file (µDST), which is then normally used for the physics analysis.

37


4 Data Analysis
In order to get access to the Beam Spin Asymmetry (BSA), discussed in section 2.4, several data analysis steps have to b e p erformed. An exclusive event has to fulfill sp ecific conditions to b e considered as a DVCS-BH event. In a first step only analyzable data are selected. In a second step events within certain kinematic b oundaries according to a DIS event are accepted. After that, only those DIS events are selected that have exactly one charged track (the scattered p ositron), and one photon within certain constraints. These events are called single-photon events. Since the HERMES sp ectrometer can not detect the recoiling nucleon, the missing mass is calculated by using the kinematics of the scattered p ositron and the real photon. In order to ensure exclusivity, only single-photon events are selected which fulfill a certain missing mass constraint. In the last step the selected exclusive event sample has to b e separated in a coherent and incoherent part. The several steps of event selections and the according data treatment are discussed in the following sections.

4.1 Event Selection
4.1.1 Data Sample
The data for this analysis were accumulated during the 2000 running p erio d of HERA using a p olarized p ositron b eam with an average p olarization of 54%. The b eam p olarization in b oth helicity states was approximately the same, for parallel p olarization 52% and for antiparallel p olarization 55%. The p ositrons were scattered off an unp olarized neon gas target. The measurements are based on an integrated luminosity of ab out 82 pb-1 . The µDST-pro duction 00d0 is used which amounts to 845 runs.

4.1.2 Data Quality
In order to select data with sufficient quality for the physics analysis two different lists have to b e checked. The first quality check is done on the run level. Only runs are taken into account that are marked with the following criteria in the electronic logb o ok:

38


· The target was normal or high density neon. · The run was marked as analyzable. The second quality check is based on the slow-control data, written for every burst. For that the burstlist-file with the condition badbit &0x503E13DC is used, containing the following data-quality cuts: · The run is marked as analyzable in the logrun file. · There are no HV trips in the tracking chamb er. · There are no dead blo cks in the calorimeter or luminosity monitor. · The preshower detector is working correctly. · The TRD is fine during the burst. · The DAQ live time is greater than 0.8. · The burst is not the first one in a run. · The b eam current is reasonable (5 mA < I < 50 mA). Three additional cuts are applied: · The live time for trigger 21 is b etween 0.8 and 1 · The raw luminosity is b etween 5 and 3000. · The fitted b eam p olarization is b etween 30% and 80%.

4.1.3 DIS Event Selection
The DIS event selection is based on several detector acceptance cuts and kinematic constraints for p ositrons. In addition the numb er of the recorded DIS events is imp ortant for normalization purp oses. The following cuts for identifying a DIS event are required: · The combined resp onse from preshower and calorimeter (P I D 2) and the resp onse of the TRD (P I D 5) satisfy the requirement P I D 2 + P I D 5 > 2. This corresp onds to a 102 times larger probability, that the track is a lepton and not a hadron. · Each event contains exactly one p ositron. · Each event is triggered by trigger 21.

39


· The track passes the fiducial volume cuts. · The track traverses the full length of the sp ectrometer. · The distance of closest approach of the track to the b eam should b e in the longitudinal direction -18 < zvtx < 18 cm and in the transverse direction |t vtx | < 0.75 cm.
l · The (x, y ) p osition of the lepton in the calorimeter has to ob ey |x calo | < 175 cm l and 30 < |ycalo | < 108 cm. This corresp onds to the track not b eing incident in the outermost 2/3 of the outer row or column of calorimeter blo cks, in order to avoid that the energy might leak out of the sides of the calorimeter.

· The kinematical cuts on the detected lepton are: ­ The photon virtuality Q2 should b e large in order to b e in the hard scattering regime. Therefore the data is constrained to Q 2 > 1 GeV2 . ­ The invariant mass of the virtual-photon-proton system is required to b e W >9GeV2 . ­ The energy of the virtual photon should b e < 22 GeV. The distribution of DIS events for different variables is shown in figure 4.1. Note that the variables are calculated by using the proton mass for the reacting part of the neon nucleus. The fact that the distributions for b oth cases, proton and neon mass, are very similar, will not b e further discussed here. Hence the variables in this chapter are all calculated on the basis of the proton mass. A comparison b etween the two metho ds of using proton or neon mass is discussed in more detail in section 5.3.
2

4.1.4 Single-photon Selection
Since the single photon event is characterized by only one DIS lepton and one photon in a certain kinematic range, the following cuts are required: · Only DIS events with exactly one charged track are taken. · The photon requires exactly one cluster in the calorimeter with no track assigned to it. · A signal in the preshower detector with E
pr esh

> 0.001 GeV.

· The photon energy in the calorimeter is E > 5 GeV. · The fiducial volume of the calorimeter for photons is given by |x p 33 < |ycalo | < 105 cm.
p calo

| < 125 cm and

40


600 500 400 300 200 100 0 4

NDIS

NDIS 6 8 10 12 14 16 18 20 22 24 P (GeV) NDIS

2500 2000 1500 1000 500 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 xBj

3500 3000 2500 2000 1500 1000 500 0 0

NDIS

2

4

6

8

10

12

14 16 2 2 Q (GeV ) NDIS

600 500 400 300 200 100 0

10

15

20

25

30

35 40 2 W2 (GeV )

600 500 400 300 200 100 0

NDIS

2500 2000 1500 1000 500 0

4

6

8

10 12

14 16 18

20 22 (GeV)

-60

-40

-20

0

20 40 60 Polarization (%)

Figure 4.1: Distribution of DIS events scattered off a neon target in dep endence on different kinematical variables.

· The angle b etween virtual and real photon,



, has to b e b etween 3 and 45 mrad.

The distributions of single-photon events are presented in figure 4.2 for different kinematic variables.

4.1.5 Exclusive Event Selection
The selection of exclusive events has to b e based on a missing mass calculation, since there is no p ossibility to detect the recoiling nucleon with the present HERMES detector. 2 The missing mass squared Mx is given as
2 Mx = (q + p - q ) 2

(4.1)

with the four momenta of the virtual photon q , of the target nucleon (nucleus) p and the real photon q . Since the energy resolution of the HERMES sp ectrometer is not p erfect,

41


0.25 0.2 0.15 0.1 0.05 0 -5 0 5 10 15 20 25 30 2 M2 (GeV) x

0.16 0.14 0.12 0.1 0.08 0.06 0.04 0.02 0

N/(1000*NDIS)

N/(1000*NDIS)

5

10

15

20

25 E

0.45 0.4 0.35 0.3 0.25 0.2 0.15 0.1 0.05 0 0

N/(1000*NDIS)

N/(1000*NDIS)

0.5 0.4 0.3 0.2 0.1 0 0

0.1

0.2

0.3

0.4

0.5

0.6 X
Bj

2

4

6

8

10

12 14 2 2 Q (GeV)

N/(1000*NDIS)

N/(1000*NDIS)

0.08 0.07 0.06 0.05 0.04 0.03 0.02 0.01 0 0 0.0050.01 0.0150.02 0.025 0.030.035 0.040.045 * (rad)

0.1 0.08 0.06 0.04 0.02 0 -3 -2 -1 0 1 2 3 (rad)

Figure 4.2: Distributions of the single-photon event yields in dep endence on different kinematic variables, normalized p er 1000 DIS events.

2 2 events can b e reconstructed with negative values of M x . The Mx -distribution for neon is shown in figure 4.3. The exclusive p eak for this distribution is smeared out (negative 2 values for Mx are defined as Mx = - -Mx ).

In order to ensure exclusive events, the missing mass constraint -1.5 GeV < M x < 1.7 GeV is used, as indicated by the two lines in figure 4.3. This sp ecific missing mass interval is found by optimizing the signal-to-background ratio using semi-inclusive background [Ell]. In fact, the missing mass distribution is not the only kinematic distribution which shows an unphysical extension to negative values. Also the -t distribution is smeared, caused by the energy resolution of the calorimeter. Since the extraction of the BSA as a function of -t is imp ortant for theoretical predictions, the following reconstruction of -t is p erformed [Ell]. The variable Mx and -t are related to each other via
2 2 Mx = Mp + 2Mp ( - E ) + t.

(4.2)

By assuming that the missing mass squared is equal to the proton mass M x = Mp , it

42


N/(1000*NDIS)

0.25

0.2

0.15

0.1

0.05

0

-5

0

5

10

15

20

25

30 2 M2 (GeV ) x

2 Figure 4.3: Distribution of single photon events versus the missing mass squared M x .

is p ossible to calculate -t without the knowledge of the photon energy E . Hence the photon energy E is calculable via E = t + 2Mp 2 + Q2 cos ) (4.3)

The calculation of the four-momentum transfer squared is given as t = (q - q )2 = -Q2 - 2E ( -


(4.4)

By inserting equation 4.3 in equation 4.4, the momentum transfer squared can b e calculated without any dep endence on E via tc = -Q2 - 2 ( - 1+
1 Mp

2 + Q2 cos 2 + Q2 cos



)

( -



)

.

(4.5)

The result is called the constrained momentum transfer t c . The effect of this calculation is shown in figure 4.4, comparing the momentum transfer -t to the constrained momentum transfer -tc . Note that this metho d is completely correct for the imp ortant "3-particle final states" (e p ), while it is not exactly correct for the (less imp ortant) background events ( 0 , , K 0 ). Summarizing, the exclusive sample is constrained by the following kinematic b oundaries: 0.03 2.25 GeV < < xB j 2 Mx -tc < < < 0.35 2.89 GeV 0.7 GeV2

2

2

43


N/(1000*NDIS)

0.16 0.14 0.12 0.1 0.08 0.06

N/(1000*NDIS)

0.3

0.25

0.2

0.15

0.1

0.04 0.02 0
0.05

-0.5

0

0.5

1 1.5 2 -t (GeV )

0 0

0.1

0.2

0.3

0.4

0.5

0.6 0.7 2 -t c (GeV )

Figure 4.4: Comparison b etween the momentum transfer -t and the constrained momentum transfer -tc . Both distributions are for the exclusive sample, normalized p er 1000 DIS events.

The kinematic distributions of exclusive events are shown in figure 4.5 for different variables.

4.2 Extraction of the Beam-Spin Asymmetry
The azimuthal dep endence of the BSA can b e extracted by fitting the cross-section asymmetry as a function of the azimuthal angle . The b eam spin asymmetry A LU is defined as
2

A

LU

() =

d -d






d =
0 2

- d d - d d

- +

- d d

d +d

,
- d d

(4.6)

d
0

where L denotes the longitudinal p olarized b eam and U the unp olarized target. In order to extract the b eam spin asymmetry A LU from the data, the cross-section ratio is calculated for every -bin by counting the resp ective exclusive events. A LU is then

44


N/(1000*NDIS)

N/(1000*NDIS) 0.05 0.1 0.15 0.2 0.25 0.3 0.35

0.18 0.16 0.14 0.12 0.1 0.08

0.1 0.08 0.06 0.04 0.02 0 0

0.06 0.04 0.02 0 2 4 6 8
2

10

x
N/(1000*NDIS)

Bj

Q (GeV )
N/(1000*NDIS)

2

0.045 0.04 0.03 0.02 0.015 0.01 0.005 0 4 6 8 10 12 14 16 18 20 22 24 0.035 0.025

0.06 0.05 0.04 0.03 0.02 0.01 0 0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045

E (GeV)

* (rad)

Figure 4.5: Kinematic distributions of exclusive events with normalization p er 1000 DIS events

given as A
LU

() =

1 P

- N () -- -- N or m - N () -- -- N or m

- +

- N () -- -- N or m - N () -- -- N or m

(4.7)

-- with N ( N ) representing the exclusive yield in the helicity state parallel (antiparallel) to the b eam direction. The average b eam p olarization P for the neon data in 2000 -- -- -- -- was 54%. For normalization, denoted by N or m and N or m,the numb ers of DIS events in each helicity state were used. Another p ossibility for a normalization is based on the luminosity, calculated as the pro duct of the rate measured in the luminosity monitor, the live time of trigger 21 and of the burst length. Both metho ds for normalization are compared in section 5.4. for the neon data. The b eam spin asymmetries on neon for each -bin are plotted in figure 4.6 with the resp ective statistical error bars. The harmonics sin and sin 2 of the BSA are obtained by fitting the asymmetries over . As already discussed in section 2.4. the sin amsin plitude, denoted by ALU , is correlated to the amplitude of the BH-DVCS-interference

45


0.3 0.2 0.1 0 -0.1 -0.2 -0.3 -0.4 -0.5 -3 -2 -1 0

2 / ndf 9.777 / 7 -0.004471 ± 0.03008 p0 -0.2084 ± 0.04387 p1 0.09999 ± 0.0426 p2

A

LU

Fit: p0 + p1sin() + p2sin(2)

1

2

3 (rad)

Figure 4.6: The calculated b eam-spin asymmetry A LU on neon for each -bin with the resp ective statictical error bars. The fit-function is given as p 0 + p1 sin + p2 sin 2.

term. Therefore the asymmetry fit for the BSA is usually done with the fit function f () = p0 + p1 sin + p2 sin 2. (4.8)

Other fit functions are discussed in section 5.5. The sin 2 term is allowed to app ear, due to higher order in s , as explained in section 2.4. A p ossible constant term p 0 can only arise from a helicity-dep endent, but -indep endent term in the cross-sections. The constant term p0 is studied in section 5.10. The calculation of the statistical error for A dA
stat LU LU

is based on the Poisson distribution: (4.9)

=

2 P

-- -- -- -- N or m N or m N N (N + N ) - - - - . - -- -- ( N N or m+ N N or m)2

Note that the Poisson distribution b ecomes a Gaussian distribution for higher values of N. In figure 4.7 the amplitude of the sin harmonic, A sin , is plotted vs. missing mass Mx LU for neon. The upp er cut on tc and the upp er cut on have b een removed here in order to reduce the error bars in the non-exclusive region. The asymmetry is largest for the first two bins in Mx where the contribution of the background pro cesses is smallest. The third bin shows a slightly smaller asymmetry since background contributions start to dilute the asymmetry. The fourth bin yields a mixture of the exclusive and the

46


sin( ) LU

0.1

A -0 -0.1 -0.2 -0.3 -0.4 -2 0 2 4 6 Mx (GeV)

Figure 4.7: The asymmetry amplitude A sin vs. missing mass Mx for neon. The upp er LU cut on tc and the upp er cut on have b een removed in order to reduce the error bars in the non-exclusive region.

non-exclusive sample and consequently shows a decreased asymmetry value. The last three bins should not yield any contribution from the BH-DVCS interference anymore and there is indeed no sign for a significant negative amplitude anymore. The slightly p ositive asymmetry in the non-exclusive region shows the known b eam spin asymmetry for 0 pro duction, which is carried through the 0 decay to the one detected photon.

4.3 Hydrogen Data
Since hydrogen data is required for the interpretation in chapter 7, the currently available hydrogen results will b e presented and compared to the results obtained in the present study. In figure 4.3 the extracted asymmetry A LU for hydrogen is compared b etween previous released results [Now05] and non-released results obtained from this analysis. Note that the constraints are changed to 0.002 rad < < 0.7 rad for this comparison. The only difference b etween b oth data sets is due to different pro ductions of 2000 data: the results of [Now05] are extracted from the pro duction 00c0, whereas the present analysis is based on the newest pro duction 00d0. The asymmetry amplitudes are

47


ALU

0.8 0.6 0.4 0.2 0

+

e pe X
+

(Mx< 1.7 GeV)
(refined)

HERMES PREL. 2000

P1 + P2 sin + P3 sin 2

-0.2 -0.4 -0.6 -0.8 -1 -3 -2 -1 0 1 2 3
P1 = -0.04 ± 0.02 (stat) P2 = -0.18 ± 0.03 (stat) P3 = 0.00 ± 0.03 (stat) <-t > = 0.18 GeV , = 0.12, = 2.5 GeV
2 2 2

(rad)
LU

0.4 0.3 0.2 0.1 0 -0.1

A

f()=p0+p1sin+p2sin2

-0.2 -0.3 -0.4 -0.5

2 / ndf p0 p1 p2
-3 -2

7.95 / 7 -0.03 ± 0.02 -0.19 ± 0.03 -0.00 ± 0.03
-1 0 1 2 3 (rad)

Figure 4.8: Comparison of ALU for hydrogen b etween previous released results [Now05] and non-released results obtained from this analysis.

const p sin ALU Asin 2 LU

0

Previous Study Released [Now05] -0.04±0.02 -0.18±0.03 0.002±0.03

Present Study Not Released -0.03±0.02 -0.19±0.03 0.00±0.03,

and are found in go o d agreement. A comparison of A sin in dep endence on tc is presented LU in figure 4.9. The only available tc dep endence for Asin is presented in [Ell] and is LU

48


sin LU

0.3

New 3 Bins
0.2 0.1 -0

A

Ellingh 4 Bins

-0.1 -0.2 -0.3 -0.4 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 2 tc (GeV )

Figure 4.9: Comparison of the asymmetry amplitude A sin in dep endence of tc for hydroLU gen b etween previous results [Ell] and non-released results obtained from this analysis.

taken for this comparison. Note that there exists a different t c binning. In summary, the hydrogen results obtained in the present analysis are compatible to the released ones and therefore they will b e taken for further comparison with neon in chapter 7.

4.4 Separation of Coherent and Incoherent Part
The DVCS reaction on neon pro ceeds through two different pro cesses. On the one hand, there is the coherent pro cess that involves the nucleus as a whole, on the other hand, the incoherent pro cess as the reaction on a single nucleon. Since coherent and incoherent pro cesses contribute to the photon pro duction cross section, b oth parts have to b e separated. As has b een explained in section 2.5., the distribution of the momentum transfer -t is a reasonable to ol for a separation b etween coherent and incoherent part. In figure 4.10 the tc distribution is shown for the exclusive sample of neon. Note that the yield is on a logarithmic scale. In fact, the t c distribution shows a different slop e at low and high tc , corresp onding to coherent and incoherent part. In order to estimate the separation p oint psep , a double exp onential function is fitted to b oth parts simultaneously: (4.10) f () = ep0 +p1 x + ep2 +p3 x
exp
coh

exp

incoh

For the calculation of the transition p oint x 0 , b oth exp onential functions have to b e set

49


log N

10

3

Fit-Function: exp[0]+[1]*x+exp
10
2

[2]+[3]*x

10

1 0

p0 6.74 ± 0.06 p1 -55.0 ± 3.0 p3 4.68 ± 0.07 p4 -8.13 ± 0.29
0.1 0.2 0.3 0.4 0.5 -t
c

Figure 4.10: Distribution for the exclusive sample of neon. The y-axis is set logarithmically. The distribution is fitted by e p0 +p1 x + ep2 +p3 x , corresp onding to the coherent and incoherent part. The vertical dashed line indicates the p oint where b oth exp onential functions are equal.

equal exp and from that x0 is obtained as x0 =
coh

= exp

incoh 2 1

(4.11) (4.12)

By using the fit result, given in 4.10, x 0 was found to b e at tc = 0.044 ± 0.006 GeV2 . A comparison of the -tc distribution b etween neon and proton is presented in figure 4.11. The parameter of the fit results are given as Neon 6.74 ± 0.07 -54.99 ± 2.95 4.68 ± 0.07 -8.13 ± 0.30 Hydrogen 5.4 ± 0.2 -19.1 ± 3.7 5.6 ± 0.2 -7.9 ± 0.5 ,

p0 - p p3 - p

coh. constant coh. slop e incoh. constant incoh. slop e

p0 p1 p2 p3

indicating a clear difference in the slop e b ehavior b etween b oth pro cesses, as exp ected. Another p ossibility for the separation of coherent and incoherent part is offered by Monte-Carlo (MC) studies. Such studies allow in addition to estimate the contribution

50


N/(1000*NDIS)

10

-1

Neon Hydrogen

10

-2

10

-3

0

0.05

0.1

0.15

0.2

0.25

0.3 2 -t c (GeV )

Figure 4.11: Comparison of the tc distribution b etween neon and proton, normalized p er 1000 NDI S .

of coherent and incoherent to the full exclusive sample and will b e discussed in section 5.1.2. In the following the separation p oint, p sep , was set at tc = 0.045 GeV2 . It is imp ortant to note that psep do es not make a clear cut b etween coherent and incoherent part, b ecause each part contaminates the other. Further studies on the influence of the separation p oint on the extracted BSA will b e discussed in section 5.2. A comparison b etween coherent and incoherent part for the kinematic distributions is shown in figure 4.12. The separated BSA results for neon are presented in figure 4.13

51


N/(1000*NDIS)

0.04

N/(1000*NDIS)

0.08 0.07 0.06 0.05

0.035 0.03

0.025 0.02 0.015 0.01 0.005 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 x 0.2
Bj

0.04 0.03 0.02 0.01 0 1 2 3 4 5 6 2 2 Q (GeV)

N/(1000*NDIS)

0.025 0.02

N/(1000*NDIS) 6 8 10 12 14 16 18 20 22 E (GeV)

0.03

0.035 0.03

0.025 0.02 0.015

0.015 0.01

0.01 0.005 0 0.005 0 0.005 0.01 0.015 0.02 0.025 * (rad)

Figure 4.12: Kinematic distributions of exclusive events for coherent data (dashed) and incoherent data (solid), normalized p er 1000 DIS events.

Coherent
LU

2 / ndf p0 p1

11.53 / 7 0.06476 ± 0.03473 -0.2031 ± 0.05072 0.06222 ± 0.04883

Incoherent
LU

A

0.4

p2

A

0.4

6.589 / 7 2 / ndf -0.06155 ± 0.04148 p0 -0.2444 ± 0.05998 p1 0.04812 ± 0.05912 p2

0.2 0.2 0 0 -0.2 -0.2 -0.4 -0.4 -3 -2 -1 0 1 2 3 (rad) -0.6 -3 -2 -1 0 1 2 3 (rad)

Figure 4.13: BSA for coherent (left) and incoherent part (right) on neon.

52


5 Systematic Studies
In order to evaluate the systematic uncertainties of the BSA results, a numb er of different studies were p erformed and are describ ed in the following sections.

5.1 Monte Carlo Studies
5.1.1 The HERMES Monte Carlo Pro duction
The HERMES Monte Carlo (HMC) program simulates a variety of physics pro cesses using different programs on the event generation level. At first, a Generator Monte Carlo (GMC) program is used to generate the final state particles in each event. This information is then pro cessed by HMC, whereby each particle is tracked through the HERMES sp ectrometer by the detector description co de GEANT, and the resp onse of the detector is digitized. The output of HMC is very similar to the output of the data deco der, HDC, except that it also includes all the information generated by GMC. The Monte Carlo co de is sub divided into units, called packages. For example, the JETSET package [Sjo94] mo dels the fragmentation pro cess, based on the LUND string hadronization mo del [AGIS83]. In order to simulate the BH pro cess, the RADGEN [ABR98] program is used, which takes into account the lowest order radiative pro cesses in QED for p olarized and unp olarized DIS. Note that the unp olarized case can b e also calculated directly using the Mo and Tsai formalism [MT69].

5.1.2 Coherent and Incoherent Contribution
In comparison to data, where no simple distinction b etween coherent and incoherent contribution is p ossible, the MC simulation offers the advantage to provide a clear separation of coherent and incoherent pro cesses. For this study, DIS events on neon are generated via HMC and exclusive events are extracted by using the kinematic constraints discussed in section 4.1. A comparison b etween the exclusive MC event and data event distribution for tc is presented in figure 5.1. Both distributions show the same t c b ehavior and are found to b e in go o d agreement.

53


log N/(1000*NDIS)

10

-1

Neon data Neon MC

10

-2

10

-3

0

0.05

0.1

0.15

0.2

0.25

0.3

0.35

0.4

0.45 2 -t c (GeV )

Figure 5.1: Comparison of the -tc distribution b etween MC (dashed line) and data (solid line) for DVCS on neon, normalized p er 1000 DIS events.

In order to estimate the coherent and incoherent contributions to the full cross-section, the MC events can b e clearly separated in coherent and incoherent pro cesses. In the left panel of figure 5.2 the tc distribution is shown for coherent, incoherent, semi-inclusive and asso ciated BH pro cesses. The right panel of the figure shows the fractional contribution of these pro cesses to the full cross section. As exp ected, the low t c region is dominated by coherent events, whereas the higher t c region is dominated by incoherent events. The transition p oint, where b oth pro cesses have ab out the same contribution to the full cross-section, is found to b e at -t c = 0.045 GeV2 . This is in go o d agreement to the separation p oint tsep as calculated from data and discussed in 4.3. By setting tsep at 0.045, the incoherent (coherent) admixture in the coherent (incoherent) part is calculated to b e 16% (8%). Note that in the coherent (incoherent) part exists a 4% (27%) contribution from asso ciated BH and semi-inclusive events.

5.2 The Influence of the Separation Point on the BSA
In the following study a variation of the coherent-incoherent separation p oint t sep is p erformed. In figure 5.3 the BSAs Asin and Asin 2 are shown in dep endence of tsep . LU LU The differences of Asin in dep endence on tsep are found to b e less then 0.02 in the region LU of ± 0.015 from the separation p oint at -t c = 0.045 GeV2 for the coherent as well as the incoherent part. In the case of Asin 2 the differences are less then 0.03 for the coherent LU and less then 0.05 for the incoherent part.

54


N/(1000*NDIS)

0.12 0.1

Full sample Coherent Incoherent
Assoc. BH and semi-incl.

Fraction

Coherent
1

Incoherent
Assoc. BH und semi-incl

0.8

0.08 0.6 0.06 0.04 0.02 0 0 0.4

0.2

0.05

0.1

0.15

0.2

0.25 0.32 tc (GeV )

0 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 2 0.1 -tc (GeV )

Figure 5.2: MC study of the -tc distribution, separated in coherent, incoherent, semiinclusive and asso ciated BH pro cesses (left panel). The right panel shows the fractional contribution to the full cross-section.
sin(2 ) ) LU sin(2 ) ) LU

0.2

(A

(A A
sin( ) LU

0.1

A

sin( ) LU

0.05 0

0.1

0

-0.05 -0.1 -0.15

sin( ALU )

Asin(2) LU

-0.1

A

sin() LU

Asin(2) LU

-0.2 -0.2 -0.25 0.03 0.035 0.04 0.045 0.05 0.055 0.06 2 cutting point -tc (GeV ) -0.3 0.03 0.035 0.04 0.045 0.05 0.055 0.06 2 Cutting point -tc (GeV )

Figure 5.3: The BSA results Asin and Asin 2 in dep endence of different separation LU LU p oints tsep for coherent (left) and incoherent part (right). In order to avoid a sharp cut in tc , the asymmetries can also b e calculated by weighting the corresp onding events. The weighting dep ends on the fraction of the coherent or incoherent contribution to the full cross-section, shown in the right figure 5.2. The result of the asymmetries, extracted by weighting, is given for the coherent data as A A
sin LU sin 2 LU

= -0.206 ± 0.057 (t = 0.057 ± 0.055 (t

sep

sep

= -0.045 : 0.062 ± 0.049),

= -0.045 : -0.203 ± 0.051)

55


shown in comparison to the results calculated by a sharp cut t sep . Note that the is still constrained at -tc < 0.07, b ecause the contribution of semi-inclusive and ated BH overwhelm the coherent contribution b eyond this p oint. In comparison separation p oint metho d, which uses t sep = 0.045, the differences are very small. the separation p oint metho d is correct and will b e used in the further analysis.

sample asso cito the Hence

5.3 Comparison of the Calculation Metho d with Proton and Neon Mass
The calculation of the kinematic variables can b e done either by using the proton mass (0.938 GeV) or the neon mass (18.79 GeV for Ne-20). This calculation must in principle b e done in accordance with the actual pro cess, which could b e either coherent or incoherent. As the coherent pro cess describ es a reaction on the nucleus as a whole and the incoherent pro cess is characterized by a reaction on a single nucleon, the corresp onding masses should b e the neon or the proton/neutron mass. It is imp ortant to note that this systematic study refers only to the mass definition and not to a comparison b etween proton and neon data. The calculation based on the neon mass changes some of the kinematic values, hence several kinematical cuts have to b e adapted: · Invariant mass of the system of virtual photon and proton: W 2 > 520 GeV2 . · Bjorken variable: 0.016 < x
Bj

< 0.0186

· The missing mass cut Mx is changed to 289 < Mx < 397, derived from a comparison of the Mx distribution for b oth calculation metho ds. The comparison b etween the neon mass and proton mass calculation for the M x distribution is shown in figure 5.4. Note that the distribution calculated with the neon mass is scaled to that with the proton mass as M
P r oton x

=

M

N eon x

21

- 16,

(5.1)

from where the changed b oundaries in M x were calculated. A p ossible effect on the tc distribution by the chosen mass can not b e seen from figure 5.4. Concerning to the BSA results, the differences are negligible, presented in the following table for the coherent part:

56


N

log N

800 700 600 500 400

Proton mass Neon mass

10

3

Proton mass Neon mass

10 300 200 100 0 -5 0 5 10 15 2 20 25 30 M2 (GeV ) (scaled on proton mass) x

2

0

0.1

0.2

0.3

0.4

0.5 2 -t c (GeV )

Figure 5.4: Comparison of Mx distribution (left) and tc distribution (right) when calculated based on neon mass or proton mass, resp ectively. Note that the M x distribution calculated with the neon mass is scaled to that with the proton mass. Proton mass 0.065±0.035 -0.203±0.051 0.062±0.049 Neon mass 0.055±0.034 -0.20±0.05 0.042±0.048

const p Asin LU Asin 2 LU

0

In order to avoid p ossible mismatches, the calculation with the proton mass is applied b oth for the coherent and the incoherent part in this analysis.

5.4 Normalization metho d
The asymmetry, given as 1


- N () -- -- N or m - N () -- -- N or m

A

LU

() =

- +

- N () -- -- N or m - N () -- -- N or m

,

-- -- -- -- has to b e normalized for different helicity states to the factors N or m and N or m. The normalization can b e done either by the corresp onding numb er of DIS events or by the measured luminosity. Since the normalization could effect a constant-term p 0 in the asymmetry, it is useful to check the normalization on consistency. In the following, b oth metho ds are compared. The normalization p er numb er of DIS events is the standard metho d and was applied in the previous sections (e.g. figure 4.13). The kinematic b oundaries and acceptance cuts for the DIS event selection were given in 3.1.2. In order to apply the luminosity metho d,

57


LU

A

A

LU

0.6

11.55 / 7 2 / ndf 0.09515 ± 0.03468 p0 -0.2029 ± 0.05067 p1 0.06173 ± 0.04876 p2

0.4

0.4 0.2 0.2 0 0 -0.2 -0.2 -0.4 -0.4 -3 -2 -1 0 1 2 3 (rad) -3 -2 -1 0

2 / ndf 6.418 / 7 -0.03452 ± 0.04167 p0 -0.2445 ± 0.06021 p1 0.05209 ± 0.05945 p2

1

2

3 (rad)

Figure 5.5: The BSA for neon, p erformed by using the luminosity metho d, for the coherent (left) and incoherent (right) part. the pro duct of the rate in the luminosity monitor, the live time of trigger 21 and of the burst length has to b e calculated. The BSA extracted by using the luminosity metho d instead is shown in figure 5.5 for the coherent and incoherent part. A comparison b etween b oth normalization metho ds is given in the following table: Norm p er DIS 0.065±0.035 -0.203±0.051 0.062±0.049 -0.062±0.042 -0.244±0.060 0.048±0.059 Norm p er Lumi 0.095±0.035 -0.203±0.051 0.062±0.049 0.035±0.042 -0.244±0.060 0.052±0.059

Coherent

Incoherent

const p Asin LU Asin 2 LU const p Asin LU Asin 2 LU

0

0

Summarizing, it is shown that the metho d of normalization do es not effect the BSAs Asin and Asin 2 , as exp ected. The constant term p0 is shifted up by ab out one sigma LU LU in the luminosity metho d.

5.5 Study of different Fit Metho ds
5.5.1 The Influence of Binning
In order to estimate the effect of binning for the azimuthal angle for the BSA, a sin calculation with different numb er of -bins is p erformed. In figure 5.6 A LU and Asin 2 LU are plotted in dep endence on the numb er of -bins. For the coherent part the average

58


sin(2) ) LU

sin( ALU ) (A

0.05 0

A

sin( ) LU

(A

0.1

sin(2 ) ) LU

0.15


sin(2) >=0.05± LU

0.02

0.1

sin(2) > LU

= 0.053 ± 0.024

-0

-0.05 -0.1 -0.15 -0.2 -0.25


sin() LU >=-0.206±0.021

-0.1


sin() > LU

= -0.241 ± 0.027

-0.2

-0.3
-0.3 8 9 10 11 12 13 Number of -bins

8

9

10

11

12 13 Number of -bins

Figure 5.6: Asin and Asin 2 in dep endence on the numb er of -bins for the coherent LU LU (left) and incoherent (right) part. Asin (Asin 2 ) is -0.206 (0.05) with an uncertainty of ±0.021 (±0.02). The numb er of LU LU -bins varied from 8 to 13. For the incoherent part the average A sin (Asin 2 ) is -0.241 LU LU (0.053) with an uncertainty of ±0.027 (±0.024). In comparison to the usual fit metho d with 10 -bins the difference is less then 0.01 for A sin and less than 0.02 for Asin 2 , i.e. LU LU can b e neglected at the present level of statistics.

5.5.2 The Fit Function
The standard fit-function is given as f () = p0 + p1 sin + p2 sin 2 (5.2)

In order to study the influence of additional harmonics, the following fit functions are applied: f () = p1 sin f () = p0 + p1 sin + p2 sin 2 + p3 sin 3 + p4 sin 4 f () = p0 + p1 cos + p2 cos 2 (5.3) (5.4) (5.5)

The corresp onding BSA results are shown in figure 5.7 for the coherent part. Note that the amplitudes can not change by using additional parameters in the fit-function, e.g. shown in the top part in figure 5.7 for A sin : LU p1 sin A
sin LU sin LU

p0 + p1 sin + p2 sin 2 A

= -0.203 ± 0.051

= -0.208 ± 0.05

The results of this study can b e summarized in the following:

59


2 / ndf 16.75 / 8 0 ± 1.414 p0 -0.2084 ± 0.05048 p1
0.4

0.4

2 / ndf 11.53 / 7 0.06476 ± 0.03473 p0 -0.2031 ± 0.05072 p1 0.06222 ± 0.04883 p2

LU

A

0.2

A 0.2 0 -0.2

0

-0.2

LU

-0.4 -3

Fit: p1sin()
-2 -1 0 1 2 3 (rad)

p0 + p1sin() + p2sin(2)
-0.4 -3 -2 -1 0 1 2 3 (rad)

0.4

2 / ndf p0 p1 p2 p3 p4

0.06601 -0.1989 0.05529 -0.0754 0.0513

± ± ± ± ±

7.921 / 5 0.03473 0.05086 0.04897 0.04996 0.04793

0.4

2 / ndf 26.33 / 7 0.06596 ± 0.035 p0 0.07158 ± 0.04813 p1 0.07052 ± 0.04982 p2

LU

A

0.2

A 0.2 0 -0.2

0

-0.2 p0 + p1sin() + p2sin(2) + p3sin(3) + p4sin(4 ) -3 -2 -1 0 1 2 3 (rad)

LU

p0 + p1cos() + p2cos(2)
-0.4 -3 -2 -1 0 1 2 3 (rad)

-0.4

Figure 5.7: The BSA results with different fit functions for the coherent part. ·A
sin 3 LU

and A

sin 4 LU

are given as A A
sin 3 LU sin 4 LU

= -0.051 ± 0.048

= -0.075 ± 0.05

A theoretical interpretation of these higher harmonics is only known for proton and not for nuclei. In addition they have to o small statistics and therefore they will not b e further discussed here. · The cos -dep endence is negligible within the statistical errors. Note that a nonzero cos amplitude is forbidden from theory. The same conclusions hold for the incoherent sample.

60


5.5.3 The Anti-symmetrization Fit Metho d
The anti-symmetrization fit metho d is characterized by mapping the region [- , 0] to the region [0, ], which pro jects out the odd part in in the BSA. Hence the even harmonics disapp ear, esp ecially acceptance effects. This means, anti-symmetrized BSAs, by ansatz, are not sensitive (anymore) to a p ossible constant term p 0 , as this is also an even contribution. For the b eam-spin asymmetry on an unp olarized target, two different formulae exist for the anti-symmetrization: 1st way: A
LU

can b e calculated for each -bin as A = 1 < |P | >
- - N ()- N (-) -- -- N or m - - N (-) N () -- + -- -- -- N or m N or m

LU

- +

- - N ()- N (-) -- -- N or m - - N () N (-) -- + -- -- -- N or m N or m

(5.6)

In terms of the cross-section d = d dU d
ev en in

+

dL P d

L

(5.7)

odd in

the ab ove defined asymmetry turns out to b e theoretically correct: A
LU

1P = PL

dL L d - dU d

(-PL ) +
dU d

dL d

=

dL d dU d

(5.8)

2nd way: Alternatively, ALU can b e calculated -bin: - - - 1 N () ALU = - - < |P | > N () - - - 1 N () ALU = - - < |P | > N () In terms of cross sections - P - - ALU = - P
dL d dU d dL d dU d

separately for each helicity state and - - N (-) - + N (-) - - N (-) - + N (-) - - = -ALU (5.9)

(5.10)

=

(5.11)

also the 2nd way to calculate the anti-symmetrized BSA is correct. In comparison with the standard (non-symmetrized) fit metho d b oth ways share the fact that lo cal acceptance effects cancel out. In addition, the second way provides the

61


0.1

6.866 / 3 2 / ndf -0.2176 ± 0.05076 p0 0.06373 ± 0.04893 p1

3.183 / 3 2 / ndf -0.2742 ± 0.05976 p0 0.03898 ± 0.05887 p1
-0

LU

A

-0 -0.1 -0.1 -0.2 -0.2

A -0.3 -0.4

-0.3

0.5

1

1.5

2

2.5

3 (rad)

LU

0.5

1

1.5

2

2.5

3 (rad)

Figure 5.8: The BSA obtained via the anti-symmetrization metho d for the coherent (left) and incoherent (right) part. advantage that any helicity-dep endent acceptance effects (as e.g. b eam envelop es and/or vertex p ositions dep ending on helicity) drop out, as long as they are even in . The final BSA, in the second way of antisymmetrization, can b e obtained as the weighted average of b oth measurements: -- -i -i - - w ALU - wk Ak i,k LU ALU = , (5.12) - -i k w +w i,k where i(k ) runs over helicity p erio ds (). The weighting can b e done either by 1 normalization w = N or m or by the statistical error w = (ALU )2 , which in fact should b e equivalent. Note that although the statistical uncertainty p er bin decreases by ab out 2, the resulting uncertainties of the fit stays the same, as the numb er of bins was halved. The result of the anti-symmetrization metho d via the 2nd way is plotted in figure 5.8. The asymmetries ALU are fitted with the function p0 sin + p1 sin 2.
1 Note that the weighting is done by the statistical error w = (ALU )2 for each helicity state. The result of this anti-symmetrization metho d in comparison with the standard (non-symmetrized) fit metho d is given in the following table:

non-symmetrized anti-symmetrized

Coherent Asin = -0.20 ± 0.05 LU Asin = -0.22 ± 0.05 LU

Incoherent Asin = -0.24 ± 0.06 LU Asin = -0.27 ± 0.06 LU

Summarizing, the BSA results obtained via the anti-symmetrization fit metho d are in go o d agreement with the results of the standard fit metho d.

62


N/(1000*NDIS)

0.22 0.2

0.18 0.16 0.14 0.12 0.1 0.08 0.06 0.04 0.02 0 0.8 0.85 0.9 0.95 1 1.05

Parallel Antiparallel

1.1

1.15

1.2 E/P

Figure 5.9: Distribution of the ratio b etween the lepton energy dep osited in the calorimeter E and the momentum of the lepton P , derived from the x-slop e difference b etween the front and the back track, for the parallel and antiparallel helicity state.

5.6 The Influence of the Calorimeter on the Results
5.6.1 Study of the Photon Energy Reconstruction
Since the calorimeter is the only detector able to measure the energy and p osition of the photon, it is very imp ortant for this analysis. In order to check the reconstruction of the photon energy from the calorimeter measurement, the ratio b etween the lepton energy dep osited in the calorimeter E and the momentum of the lepton P is used. The latter is derived from the x-slop e difference b etween the front and the back track. The distribution for the value E /P for the lepton sample is shown in figure 5.9 for the parallel and antiparallel helicity state. The mean value is extracted by a gaussian fit function and the result is given as E = 0.991 ± 0.001, P M ean b eing the same for parallel as well as antiparallel helicity. Hence there is no difference in the energy reconstruction from the calorimeter when comparing b oth helicity states.

5.6.2 Study of the Photon Position Reconstruction
Another imp ortant quantity for this analysis is the resolution of the p osition reconstruction of the calorimeter. In order to cross check the p osition measurement of the

63


LU

LU

A

A

0.4

0.4

No Calo shift 5mm Calo shift

No Calo shift 5mm Calo shift

0.2 0.2 0 0 -0.2 -0.2 -0.4 -0.4 -0.6 -3 -2 -1 0 1 2 3 (rad) -3 -2 -1 0 1 2 3 (rad)

Figure 5.10: Comparison b etween the results achieved with and without -5 mm correction of the photon p osition in the top part of the calorimeter for the coherent (left) and incoherent sample (right). calorimeter, the reconstructed cluster p osition for charged tracks can b e used. Note that the p osition reconstruction in the calorimeter b ehaves similarly for leptons and photons, b ecause the photons used in this analysis convert into lepton pairs in the preshower detector, which is lo cated in front of the calorimeter. Investigating the difference y track - ycalo separately for the top and the b ottom part of the calorimeter reveals that the mean of this distribution in the top part is ab out zero while it is ab out -5 mm in the b ottom part for the year 2000 [Ell]. The influence on the asymmetry ALU of this p osition shift is calculated by using a -5 mm correction of the photon p osition in the top part of the calorimeter. A comparison b etween the results achieved with and without the correction is shown in figure 5.10 for the coherent and incoherent sample. The results for the BSAs are given in the following table. No Correction 0.065±0.035 -0.203±0.051 0.062±0.049 -0.062±0.042 -0.244±0.06 0.048±0.059 5mm Correction 0.055±0.035 -0.199±0.051 0.072±0.049 0.045±0.042 -0.242±0.06 0.046±0.06

Coherent

Incoherent

const p Asin LU Asin 2 LU const p Asin LU Asin 2 LU
sin 2 LU

0

0

The influence on A

sin LU

(A

) is less than 0.005 (0.011).

64


rec. - gen. (rad)

1.4 1.2 1 0.8 0.6 0.4 0 0.002 0.004 0.006 0.008 0.01 0.012 * (rad)

Figure 5.11: Average difference b etween generated and reconstructed in dep endence on .

5.7 Smearing Effects
Smearing effects are mainly of interest as they change the relative direction of the real photon with resp ect to that of the virtual photon. Therefore they change the value of the azimuthal angle as well as the p olar angle and the momentum transfer tc . The smearing of the photon can b e studied by using MC and comparing the generated photon direction with the reconstructed photon direction. In figure 5.11 the average difference b etween the generated and the reconstructed is shown in dep endence of . The lower cut for was set to 3 mrad, in order to reduce smearing effects in . The remaining smearing for an asymmetry has b een obtained in a MC simulation in 4 bins in t. The generated asymmetry amplitude is 0.5. It is shown in figure 5.12 that smearing effects reduce the generated asymmetry by ab out 0.02 (0.04) in the coherent (incoherent) part. Hence an relative error of 4% (8%) of the asymmetry amplitude in the coherent (incoherent) sample is estimated.

5.8 Determination of the Background Contribution
The background for DVCS originates from the pro duction of photons pro duced in a semi-inclusive or exclusive background (= exclusive reaction other than DVCS). It was derived from MC studies that the main part stems from 0 meson pro duction ( 80%) and the rest is dominated by photons from the decay of the meson. The general

65


A 0.5 0.4 0.3 0.2 0.1 0 0

sin LU

0.6

0.05

0.1

0.15

0.2

0.25

0.3

0.35

0.4 0.45 2 tc (GeV )

Figure 5.12: Reconstructed Asin from a MC simulation with a generated asymmetry LU amplitude of 0.5, for different tc bins.

problem in the data analysis is not only the exclusive 0 -or pro duction, which are almost indistinguishable from BH-DVCS under certain kinematic constraints, but also semi-inclusive pro cesses can leak into the exclusive sample. Since there are statistical fluctuations in the measured energy dep osition in the calorimeter, the energy of the single photon cluster can b e overestimated such that the missing mass gets closer to the exclusive range. The fraction of the background in dep endence of t is studied in [Kra]. The upp er limit for the coherent sample is 7% and it increases to 14% for the incoherent sample. In order to determine the uncertainty of the b eam-spin asymmetry results caused by the background contribution, the asymmetry of the background has to b e estimated. Since the largest contribution stems from the 0 meson, it is imp ortant to determine the b eam-spin asymmetry in semi-inclusive 0 pro duction, which is found to b e less than 3% [HERMES04a]. In summary the uncertainty caused by background contributions is approximately given as: Abkg = 0.03 (t) where denotes the fraction of the dataset that originates from background: =
Nbkg Ntot

5.9 Combined Systematic Uncertainties
In addition to the statistical uncertainty, as given in equation 4.9, the estimation of the systematic uncertainties is also imp ortant. Therefore the results describ ed in the

66


previous sections are used. The total systematic uncertainty arises from the following uncertainties: · Photons from the semi-inclusive background and exclusive reactions other than DVCS: It was derived from MC studies that the main part stems from 0 meson ( 80%) and the rest is dominated by photons from the decay of the meson. For the coherent (incoherent) part the upp er limit of the uncertainty is 0.003 (0.005), as discussed in section 5.8. · Detector smearing: The relative error is estimated to b e 4% (8%) of the asymmetry amplitude in the coherent (incoherent) sample, as explained in section 5.7. These values include binning effects as discussed in section 5.5.1. · Uncertainty of the p olarization measurement: The measurement of the b eamp olarization is done by two p olarimeters (see section 3.1.) in parallel. Usually, the p olarization value is taken from the longitudinal p olarimeter (LPOL) since it has a smaller systematic uncertainty (1.6%). Only if there are some problems in the LPOL during op eration, the values of the transverse p olarimeter (TPOL) with a systematic uncertainty of 3.4% are taken. The resulting systematic uncertainty is 1.9% of the b eam p olarization. Hence the uncertainty of the asymmetry is 0.004 for the coherent part and 0.005 for the incoherent part. · The systematic uncertainty of a vertical calorimeter shift (alignment) was found to b e 0.005 (0.003) for the coherent (incoherent) part, as describ ed in section 5.6.2. The combined systematic uncertainty is obtained by adding in quadrature the uncertainties listed ab ove. The systematic uncertainty for the coherent part is dA and for the incoherent part dA
sin LU

(sy st) = 0.011

(5.13) (5.14)

sin LU

(sy st) = 0.021.

5.10 Study of the Constant Term
In the fits to extract the BSAs, as shown in figure 4.13, a non-negligible constant term p0 is found for the coherent part (p0 = 0.065 ± 0.035), as well as for the incoherent part (p0 = -0.062 ± 0.041). Note that the sign of p 0 is different for each part and that consequently the sum of b oth parts in figure 4.6 shows no constant term p 0 . For the coherent part, p0 is two times the statistical error (2 ), which corresp onds to a probability of 95% not to b e caused by statistical effects. Since there is no known theoretical reason for a constant term, which could explain a -indep endent different

67


Ratio

Ratio

2 / ndf 49.72 / 45 8.365 ± 0.1642 p0 + + -

16 14 12 10 8 6

12000 14000 16000 18000 20000 22000 24000 26000 28000 30000 Runnumber

Figure 5.13: Ratio of coherent DVCS to DIS events p er fill. The vertical lines indicate different helicity running p erio ds.

cross-section b etween parallel and antiparallel helicity states, the only remaining reason can b e due to systematic differences. Because of the -indep endence, p 0 is defined by the difference b etween the ratio of DVIC S for b oth helicity states, given as DS 1 p0 =


-- - - - ND V C S -- - ND I S -- - - - ND V C S -- - ND I S

- +

- - - -- ND V C S - -- ND I S - - - -- ND V C S - -- ND I S

(5.15)

In order to find the systematic reason for p 0 , several systematic studies are p erformed in the following.

5.10.1 The Time Dep endence of p

0

Since there is no physical reasons for a different ratio of DVIC S events, a strong evidence DS for a time-dep endent change of some detector prop erties exists. For example, a different ratio could b e caused by a change in the detector efficiency b etween two run p erio ds, which corresp ond to different b eam helicity states. In order to estimate the time dep endence, the ratio DVIC S is calculated p er fill for the coherent sample, as shown in figure DS 5.13. In particular, the region from run 14000 to run 16000 shows a significantly low ratio.

68


5.10.2 Study of Systematic Influences
In order to identify a systematic origin for the non-zero value of p 0 , a numb er of studies has b een done: · The standard normalization is done p er DIS events. In addition, the normalization p er luminosity is p erformed as describ ed in section 5.4. The result by using the luminosity is shown in figure 5.5 with p0 = 0.095 ± 0.035. Since p 0 do es not disapp ear, p0 can not b e caused by wrong normalization. · A different photon energy reconstruction of the calorimeter b etween b oth helicity states was also considered. Since the energy reconstruction E /P is the same for b oth helicity states, as describ ed in section 5.6.1, it could not cause a constant term. In addition there is also no difference in E /P , dep ending on the p osition in the calorimeter. · The photon p osition reconstruction was reviewed and describ ed in 5.6.2. The result, calculated with a correction of -5 mm in the top part of the calorimeter, gives a constant term p0 reduced by 0.01, as describ ed in section 5.6.2. · Another p ossibility could arise from a variation of the efficiency in the ho doscop e H0. In order to calculate the H0 inefficiency, the numb er of DIS events was counted by using trigger-18 in comparison to trigger-21. Since the only difference b etween trigger-18 and trigger-21 is the required signal in H0, the efficiency is given by the ratio of the numb er of trigger-18 events with and without trigger-21. The H0 efficiency is calculated for a grid of 2â2cm 2 across the H0 scintillator plane using the lepton tracks that passed trigger-18. The correction is done by weighting the resp ective events with the inverse efficiency for each run p erio d. The results show a negligible influence (less than 2%) for p 0 as well as for Asin and Asin 2 . LU LU · Usually the threshold of the calorimeter is set at 3.5 GeV for high density runs and 1.4 GeV for standard-density runs. In fact, the threshold was also set for several high-density runs at 1.4, which could have an influence on the acceptance. In order to estimate this influence, the BSA is calculated without these runs and compared to the standard result. The differences for A sin and Asin 2 are negligible. For p0 LU LU the difference is ab out 0.006 and can b e included in its systematic uncertainty. · Since high density runs with neon are usually following standard ABS density hydrogen data taking, a p ossible systematic deviation should b e also b e seen in the hydrogen data. Therefore, the hydrogen data of the according time p erio ds is analyzed by using the same cuts as in the neon analysis. In the hydrogen results no constant term app ears. · In addition to the studies discussed ab ove, the following kinematic variables are reviewed for differences in the distributions b etween b oth helicity states: Q 2 , W 2 ,

69


2 , xB j , yB j , tc , , lepton , P , , E , Mx , zvtx , tvtx and the b eam energy. Summarizing, a significant variation can b e only seen in t c and as well as in the Mx distribution, illustrated in figure 5.14

A combined systematic uncertainty for p 0 is obtained by adding in quadrature the uncertainties of the photon p osition and the H0 inefficiency. From that the systematic uncertainty for p0 is given as dAp0 (sy st) = 0.012. (5.16) LU This is clearly not covering the observed 2 deviations of p 0 . As a non-zero value for p0 can only b e due to exp erimental systematic problems, it must at present b e assigned to unknown sources of systematic uncertainties.

70


p0

0.4 0.3 0.2 0.1 -0 -0.1

-0.2 -0.3 0 0.05 0.1 0.15 0.2 0.25 0.3 2 -t c (GeV )

p0 0.2 0.1 0 -0.1 -0.2 -0.3 0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 (rad) p0 0.4 0.2 0 -0.2 -0.4 -0.6 -2 -1 0 1 2 M2 (GeV ) x
2

Figure 5.14: The constant term p0 in dep endence of the kinematic variables t c , 2 Mx for full sample. 71



and


6 Summary of the Results
In this chapter the b eam-spin asymmetries extracted in Deeply Virtual Compton Scattering using a neon target and their kinematic dep endences are presented . The interpretation and a comparison with hydrogen data will b e discussed in the next chapter. All data for this study was accumulated during the 2000 running p erio d of HERA using a p olarized p ositron b eam with an average p olarization of 54%. The measurements amount to an integrated luminosity of ab out 82 pb -1 . For the mostly coherent (incoherent) part the exclusive single photon sample contains 1394 (909) events in the parallel helicity state and 1470 (1102) in the antiparallel helicity state. The asymmetries have b een extracted in dep endence on the variables x their kinematic constraints, given as: 0.03 < xB j -tc Q2 < < > 0.35 0.7 GeV 1 GeV2
Bj

, -tc , Q2 within

2

A separation b etween coherent and incoherent exclusive sample is done at -t c = 0.045, whereby the coherent (incoherent) sample is constrained as -t c < 0.045 (-tc >= 0.045). The average values of the kinematic variables for the coherent and incoherent exclusive sample are: Coherent 0.068 0.019 GeV 1.75 GeV2 Incoherent 0.107 0.163 GeV2 2.76 GeV2

xB j -t c Q2

2

6.1 Coherent Sample
The result for A
LU

is plotted in figure 6.1 in dep endence on using the fit function f () = p0 + p1 sin + p2 sin 2. (6.1)

72


A

LU

0.6 0.4 0.2 0 -0.2 -0.4 -3

f()=p0+p1sin()+p2sin(2)

11.535 / 7 2 / ndf 0.065 ± 0.035 p0 -0.203 ± 0.051 p1 0.062 ± 0.049 p2
-2 -1 0 1 2 3 (rad)

Figure 6.1: Beam spin asymmetry on unp olarized neon (2000) in 10 bins of for the coherent sample. The three-parameter fit function f () and its co efficients are also shown.

Note that the parameters p0 , p1 and p2 corresp ond to the constant term, to the amplitude Asin and to the amplitude Asin 2 . A clear left/right asymmetry is visible. This is related LU LU to the negative amplitude of Asin with a statistical significance of 4 . The amplitude LU Asin 2 is p ositive with a very small statistical significance of 1.3 and therefore no LU kinematic dep endences are shown for it. In addition there is a p ositive constant term p 0 with a statistical significance of 1.9 . The asymmetry A sin is clearly non-zero: LU A
sin LU

= -0.203 ± 0.051(stat) ± 0.011(sy st)

with a significance of 3.9 when quadratically combining statistical and systematic uncertainties. The extracted amplitudes Asin for the coherent sample are shown in figure 6.2 in deLU p endence on xB j , -tc , Q2 , whereby the resp ective other two variables were integrated over in each case. The numerical results and the average kinematics of the resp ective other two variables are presented in the app endix (table 9.1). The dep endences on x B j and Q2 , shown in the top and b ottom panel, app ear to b e not very strong within the statistical uncertainty. Lo oking at the dep endence on -t c in the middle panel, the very low tc region (-tc 0.03) seems to have a significant non-zero asymmetry.

73


sin LU

A

0.1

-0

-0.1

-0.2

-0.3

0.02

0.04

0.06

0.08

0.1

0.12

0.14

0.16

0.18

0.2

x

Bj

A

sin LU

-0

-0.1

-0.2

-0.3

-0.4 0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045

-tc (GeV )
sin LU

2

A

-0

-0.1

-0.2

-0.3

-0.4

1

2

3

4

5
2

6

Q (GeV )

2

Figure 6.2: A xB j , -tc , Q2 .

sin LU

for the coherent sample in dep endence on the kinematical variables

74


A

LU

0.4 0.2 0 -0.2 -0.4 -0.6 -3

f()=p0+p1sin()+p2sin(2)

2 / ndf p0 p1 p2
-2

6.589 / 7 -0.062 ± 0.041 -0.244 ± 0.060 0.048 ± 0.059
-1 0 1 2 3 (rad)

Figure 6.3: Beam spin asymmetry on unp olarized neon (2000) in ten bins of for the incoherent sample. The three-parameter fit function f () and its co efficients are also shown.

6.2 Incoherent Sample
For the incoherent sample the extracted asymmetry A LU is shown in figure 6.3, fitted with the same function as in the case of the coherent sample (equation 6.1). A significant negative amplitude of Asin is visible: LU A
sin LU

= -0.244 ± 0.06(stat) ± 0.021(sy st)

with a significance of 3.8 when adding quadratically statistical and systematic uncertainties. A2 sin is found to b e compatible with zero and therefore no kinematic LU dep endences are shown. The constant term p 0 is small and negative with a significance of 1.5 . The extracted amplitudes Asin for the incoherent sample are shown in figure 6.4 in LU dep endence on xB j , -tc , Q2 . The numerical results are presented in the app endix (table 9.2) in conjunction with the corresp onding average kinematic values. Because of the small statistics no significant dep endence on x B j , -tc and Q2 can b e seen.

75


A
0 -0.1 -0.2 -0.3 -0.4 -0.5 0

sin LU

0.1

0.05

0.1

0.15

0.2

0.25

0.3

0.35

x

Bj

A

sin LU

-0

-0.1

-0.2

-0.3

-0.4 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7

-tc (GeV )
sin LU

2

A

0 -0.1 -0.2 -0.3 -0.4 -0.5 2 4 6 8
2

10

Q (GeV )

2

Figure 6.4: A

sin LU

for the incoherent sample in dep endence on x

Bj

, -tc , Q2 .

76


7 Interpretation in the Light of Theoretical Models
In the following, the results obtained on neon in the last chapter are compared with hydrogen results and nuclear effects predicted from theoretical mo dels are discussed. Since hydrogen consists only of one proton in comparison to a neon nucleus (Ne-20) with 10 protons and additional 10 neutrons, some nuclear effects can b e exp ected, as explained in section 2.5.

7.1 Comparison with Hydrogen Data
The hydrogen data was accumulated in the running p erio d 2000 and analyzed using the same kinematic cuts as for neon, describ ed in section 4.1. Note that the hydrogen results are compatible with released results [Now05, Ell] as shown in section 4.3. According to the fact that the neon data is separated in a coherent and incoherent sample, also the hydrogen data is separated into the same t c regions: -tc < 0.045 GeV2 for the coherent sample and -tc > 0.045 GeV2 for the incoherent sample. In the following table a comparison b etween the results for neon and hydrogen is presented. Neon 0.065±0.035 -0.203±0.051 0.062±0.049 -0.062±0.042 -0.244±0.060 0.048±0.059 Hydrogen -0.011±0.036 -0.152±0.051 -0.021±0.050 -0.048±0.026 -0.213±0.037 0.010±0.037
sin LU

Coherent

Incoherent

const p Asin LU Asin 2 LU const p Asin LU Asin 2 LU

0

0

For b oth coherent and incoherent part the amplitudes A within one sigma of the statistical uncertainty.

on neon and hydrogen agree

sin In figure 7.1 the comparison of ALU and Asin 2 in dep endence of tc is shown. The LU binning in tc has b een chosen, since it is the b est variable to discriminate b etween different pro cesses. In the b ottom part the corresp onding ratios of neon over proton

77


sin LU

sin 2 LU

0.1

A

0.3 0.25 0.2

A

Neon
Hydrogen

-0

-0.1

Neon
Hydrogen

0.15 0.1 0.05

-0.2

-0.3

0 -0.05

-0.4 0 0.1 0.2 0.3 0.4 0.5 -0.1 0
2

0.05

0.1

0.15

0.2

0.25

0.3

0.35

0.4
2

-tc (GeV )
sin sin /A LU,hyd LU,neon sin 2 sin 2 /A LU,hyd LU,neon

tc (GeV )
25 20 15 10 5

4 3.5 3 2.5 2 1.5

A

A
1 0.5 0 -0.5 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4

0 -5 -10 -15 0

tc (GeV )

2

0.05

0.1

0.15

0.2

0.25

0.3

0.35

0.4
2

tc (GeV )

Figure 7.1: Comparison of the asymmetry amplitude A sin (left panel) and Asin 2 (right LU LU panel) vs. tc for neon (filled squares) and hydrogen (empty squares). In the b ottom part the ratios of neon over proton asymmetry are shown as A sin ,neon /Asin ,hyd and LU LU
2 2 Asin ,neon /Asin ,hyd . Note that error bars with a larger value than the y-axis are cut LU LU off and hence not shown in their full length.

asymmetry are shown. Note that error bars with a larger value than the y-axis are cut off and therefore not shown in their full length. Lo oking at the left panel, the asymmetry Asin is systematically slightly larger for neon than for hydrogen in all t c bins although LU fully compatible within error bars. Also the asymmetries A sin 2 seems to b e larger LU for neon than for hydrogen in a similar way. Because the latter are compatible with zero, the ratio
A
sin 2 LU,neon sin 2 ALU,hy d

will not b e further discussed here. In the app endix (table 9.3)
sin LU,neon

the numerical results of A

and A

sin LU,hy d

with the average values of the kinematic

variables are given for each tc bin. A comparison of Asin b etween neon and hydrogen in LU dep endence on the variables xB j and Q2 is presented in figure 7.2. No significant effect Aneon is found within statistical uncertainties. A study of the ratio LU d in dep endence on hy
A

the azimuthal angle , presented in figure 7.3, shows no visible harmonics for all three

LU

78


sin LU

A

A
0 -0.1 -0.2 -0.3 -0.4 -0.5 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35

0.05 0

-0.05 -0.1 -0.15 -0.2 -0.25 -0.3 -0.35

x

sin LU

0.1

0.1

Bj

x

Bj

sin LU

0.3 0.2 0.1

A

A
0 -0.1 -0.2 -0.3 -0.4 -0.5 1 2 3 4 5
2

0 -0.1 -0.2 -0.3 -0.4 -0.5 6

sin LU

2
2

4

6

8
2

10

Q (GeV )

Q (GeV^2)

Figure 7.2: Comparison of Asin for the kinematic variables xB j and Q2 for neon (filled LU squares) and hydrogen (op en squares). The asymmetries are shown for the coherent sample (-tc < 0.045 GeV2 ) in the left panels and for the incoherent sample (-t c > 0.045 GeV2 ) in the right panels. tc bins. Hence no dep endent nuclear effects are found.

79


0.00<-tc<0.045
hyd neon /A LU LU

10

A
5 0 -5 -10

-3

-2

-1

0

1

2

3

(rad)

0.045<-tc<0.18
hyd neon LU /A LU

10

A

5

0

-5

-10

-3

-2

-1

0

1

2

3

(rad)

0.18<-tc<0.7
hyd neon LU /A LU

10

A

5

0

-5

-10

-3

-2

-1

0

1

2

3

(rad)

Figure 7.3: The ratio

A

neon LU Ahy d LU

in dep endence on for all three tc bins.

80


7.2 Mo del Calculations and Interpretation
Exploratory studies using GPDs were p erformed on the generalized EMC effect, i.e. the mo difications of the nuclear GPDs with resp ect to the free-nucleon (more precisely, deuteron) GPDs, normalized to their resp ective form factors [GS03, LT05a]. In particular it was shown in [LT05a] that the role of partonic transverse degrees of freedom, accounted by a consideration of nucleon off-shellness, is enhanced in the generalized Anucleus EMC effect. In the mo del of [LT05a] an enhancement of the ratio LUoton for higher t pr
A

is predicted. An influence on A is also predicted by Guzey and Strikman [GS03], and describ ed in section 2.5. It is shown, apart from the combinatorial enhancement of the asymmetries b ecause of the neutron contribution, that there are two additional effects: while Fermi motion of the nucleons enhances the asymmetries, the presence of the incoherent scattering for larger t reduces the asymmetry. Following the study of Anucleus Guzey and Strikman, the ratio of neon to proton asymmetry ( ALU oton ) is significantly Pr larger than unity for coherent nuclear DVCS (exp ected to b e close to the factor of two for t 0). In the case of higher t the inclusion of the incoherent contribution should decrease the ratio of the asymmetries, and part.
Anucleus LU Apr oton LU
LU

sin LU

LU

1 is exp ected for the incoherent

From the results of Asin given ab ove a difference of 0.051 (0.031) b etween neon and LU proton is obtained for the coherent (incoherent) part. The extracted ratio is 1.34±0.77 (1.15±0.48) for the coherent (incoherent) part. In figure 7.4 a comparison of the extracted ratio with the theoretical mo dels of Guzey/Strikman [GS03] and Liuti/Taneja [LT05a], as explained in section 2.5, is shown in dep endence on t c . Note that in the theoretical mo dels b oth coherent and incoherent contribution are used for the calculations. In fact, the statistical error is yet to o large to distinguish b etween the theoretical mo dels. Summarizing, the asymmetry is in all t c bins slightly higher for neon than for proton. This seems to b e in go o d agreement to all known theoretical results [GS03, LT05a]. Because of limited statistics no significant conclusions on the b ehavior of asymmetry ratios vs. tc can b e extracted.

81


sin LU,hyd

Ratio ALU,neon/A

4

sin

Model GS total Model LT total

3

2

1

0

-1 0

0.1

0.2

0.3

0.4

0.5 2 tc (GeV )

Figure 7.4:

Comparison of the ratio

A

sin LU,neon Asin ,hy d LU

vs.

t

c

to theoretical mo dels of

Guzey/Strikman (GS) and Liuti/Taneja (LT).

82


8 Summary and Outlook
The aim of this thesis was to study and to extract the b eam-spin asymmetry in hard electropro duction of real photons off neon. The data presented has b een accumulated by the HERMES exp eriment at DESY, scattering the HERA 27.6 GeV p olarized p ositron b eam on an unp olarized neon gas target. Attributed to the interference b etween the Bethe-Heitler pro cess and the deeply-virtual Compton scattering (DVCS) pro cess, the asymmetry gives access to the latter at the amplitude level. Its description is expressed in the theoretical framework of generalized parton distributions (GPDs), which offers a p ossibility to determine the total angular momentum carried by the quarks in the nucleon. The measurement of DVCS on the proton has shown the p ossibility to provide a sensitive test of current mo dels of GPDs. Similarly, studying DVCS on a nucleus op ens access to the prop erties of quark and gluon matter inside nuclei by measuring the mo dification of particle correlations enco ded in GPDs due to the nuclear environment. The DVCS reaction on neon pro ceeds through two different pro cesses, the coherent pro cess that involves the nucleus as a whole and the incoherent pro cess as the reaction on a single nucleon. Since coherent and incoherent pro cesses contribute to the photon pro duction cross section, b oth parts have to b e separated. In the present analysis the leading amplitude Asin of the b eam-spin asymmetry is obtained for the coherent and LU the incoherent sample. The extraction of these asymmetries has b een describ ed in detail and the results are presented as functions of one of the variables x B j , tc and Q2 . In addition, several systematical studies are p erformed in order to determine consistency and uncertainties of the obtained results. An explicit extraction of non-leading order effects is not p ossible, since the amplitudes of the higher harmonics are not significant within the present exp erimental uncertainty. Finally, although the amplitudes of the leading harmonic A sin have significant values, LU their interpretation is still difficult. With the present Hermes detector it is not p ossible to make a separation b etween the following pro cesses: · coherent BH-DVCS on neon, · incoherent BH-DVCS on neutron, · incoherent BH-DVCS on proton, · incoherent BH-DVCS on the neutron with excitation of a nuclear resonance.

83


· incoherent BH-DVCS on the neutron with excitation of a nuclear resonance. Hence there is presently no way to clearly disentangle the different contributions. In order to solve this problem and to improve the exclusivity of the DVCS pro cess, the installation of a recoil detector is currently under preparation [HERMES04b, Kra]. In summary, the amplitudes of the leading harmonics for the coherent and the incoherent sample are extracted and a clear evidence for a b eam-spin asymmetry in each part has b een found. In addition, a comparison b etween proton and neon asymmetries has shown a slightly higher asymmetry for neon than for proton, which is in agreement with theoretical predictions for nuclear DVCS. In the near future, the analysis of BSAs will b e extended to other nuclei. Hence, p ossible nuclear effects could b ecome more visible, including dep endences of the atomic numb er. In fact, with more statistics a distinction b etween different theoretical mo dels for nuclear mo dels seems p ossible. In addition, also the theoretical development in the field of nuclear DVCS probably needs more time. A b etter theoretical interpretation may eventually allow, through comparison to data of higher precision and also more complete kinematics information, to extract new information ab out the dynamical interplay of complex hadronic systems.

84


9 Appendix
Tables of Results for the Coherent Sample
xB j 0.03 0.06 0.09 bin - 0.06 - 0.09 - 0.20 xB j 0.049 0.073 0.110 xB j 0.058 0.071 0.081 xB j 0.055 0.075 0.106 -t (GeV 2 ) 0.016 0.019 0.027 Q
2

c

(GeV 2 ) 1.29 1.83 2.79 (GeV 2 ) 1.47 1.84 2.07 (GeV 2 ) 1.26 1.95 3.20

Asin LU -0.196 -0.236 -0.166 Asin LU -0.196 -0.305 -0.109

± ± ± ± ± ± ± ±

stat. 0.078 0.083 0.119 stat. 0.078 0.086 0.104

± sy st. ± 0.011 ± 0.012 ± 0.010 ± sy st. ± 0.011 ± 0.014 ± 0.008 ± ± ± ± sy st. 0.011 0.010 0.014

-tc bin 0 - 0.014 0.014 - 0.026 0.026 - 0.045 Q2 1.0 1.6 2.5 bin - 1.6 - 2.5 - 6.0

-t

c

(GeV 2 ) 0.009 0.019 0.034 (GeV 2 ) 0.017 0.020 0.026

Q

2

-t

c

Q

2

Asin LU -0.214 -0.158 -0.298

± stat. ±0.069 ±0.090 ±0.131

Table 9.1: The asymmetry Asin of the coherent sample p er kinematic bin at the reLU sp ective average kinematic value.

85


Tables of Results for the Incoherent Sample
xB j 0.03 0.08 0.13 bin - 0.08 - 0.13 - 0.35 xB j 0.059 0.102 0.182 xB j 0.099 0.111 0.121 xB j 0.067 0.108 0.185 -t (GeV 2 ) 0.157 0.148 0.191 (GeV 2 ) 0.075 0.174 0.374 (GeV 2 ) 0.135 0.162 0.222 Q
2

c

(GeV 2 ) 1.62 2.60 4.58 (GeV 2 ) 2.40 2.85 3.50 (GeV 2 ) 1.48 2.80 5.23

Asin LU -0.255 -0.203 -0.355 Asin LU -0.305 -0.155 -0.273 Asin LU -0.138 -0.379 -0.215

± ± ± ± ± ± ± ± ± ± ± ±

stat. 0.098 0.104 0.119 stat. 0.089 0.103 0.130 stat. 0.096 0.095 0.131

± sy st. ± 0.022 ± 0.018 ± 0.030 ± sy st. ± 0.026 ± 0.015 ± 0.023 ± sy st. ± 0.014 ± 0.031 ± 0.019

-tc bin 0.045 - 0.12 0.12 - 0.25 0.25 - 0.7 Q2 1.0 2.0 3.8 bin - 2.0 - 3.8 - 10

-t

c

Q

2

-t

c

Q

2

Table 9.2: The asymmetry Asin of the incoherent sample p er kinematic bin at the LU resp ective average kinematic value.

Table of Results for the Comparison of Neon and Hydrogen
-tc bin 0 - 0.045 0.045 - 0.18 0.18 - 0.7 xB j 0.07 0.10 0.12 -t (GeV ) 0.02 0.10 0.31
sin LU 2

c

Q

2

(GeV ) 1.7 2.5 3.3

2

Neon A ± stat. ± sy st. -0.203±0.051±0.011 -0.295±0.075±0.025 -0.157±0.099±0.015
sin LU

Hydrogen A ± stat. ± sy st. -0.152±0.051±0.010 -0.257±0.045±0.022 -0.124±0.065±0.012
sin LU

Table 9.3: Comparison of A

b etween neon and hydrogen p er tc bin.

86


Bibliography

[A+ 76] [ABR98]

M. J. Alguard et al., Deep inelastic scattering of p olarized electrons by p olarized protons, Phys. Rev. Lett. 37, 1261 (1976). I. Akushevich, H. Bottcher and D. Ryckb osch, RADGEN 1.0: Monte Carlo generator for radiative events in DIS on p olarized and unp olarized targets, (1998), hep-ph/9906408. B. Andersson, G. Gustafson, G. Ingelman and T. Sjostrand, PARTON FRAGMENTATION AND STRING DYNAMICS, Phys. Rept. 97, 31 (1983). G. Baum et al., A new measurement of deep inelastic e p asymmetries, Phys. Rev. Lett. 51, 1135 (1983). D. P. Barb er et al., The HERA p olarimeter and the first observation of electron spin p olarization at HERA, Nucl. Instrum. Meth. A329, 79­111 (1993). D. P. Barb er et al., High spin p olarization at the HERA Electron Storage Ring, Nucl. Instrum. Meth. A338, 166­184 (1994). M. Beckmann et al., The longitudinal p olarimeter at HERA, Nucl. Instrum. Meth. A479, 334­348 (2002), physics/0009047. A. V. Belitsky, D. Muller and A. Kirchner, Theory of deeply virtual Compton scattering on the nucleon, Nucl. Phys. B629, 323­392 (2002), hep-ph/0112108. M. Diehl, Generalized parton distributions, Phys. Rept. 388, 41­277 (2003), hep-ph/0307382. F. Ellinghaus, Beam charge and b eam spin azimuthal asymmetries in deeply- virtual Compton scattering, DESY-THESIS-2004-005. A. Freund and M. Strikman, DVCS on nuclei: Observability and consequences, Eur. Phys. J. C33, 53­61 (2004), hep-ph/0309065.

[AGIS83]

[B+ 83] [B+ 93]

[B+ 94] [B+ 02] [BMK02]

[Die03] [Ell] [FS04]

87


[GPRV04]

M. Guidal, M. V. Polyakov, A. V. Radyushkin and M. Vanderhaeghen, Nucleon form factors from generalized parton distributions, (2004), hepph/0410251. V. Guzey and M. Strikman, DVCS on spinless nuclear targets in impulse approximation, Phys. Rev. C68, 015204 (2003), hep-ph/0301216. C. Adloff et al. (H1 Collab oration), Measurement of deeply virtual Compton scattering at HERA, Phys. Lett. B517, 47­58 (2001), hepex/0107005. K. Ackerstaff et al. (HERMES Collab oration), HERMES sp ectrometer, Nucl. Instrum. Meth. A417, 230­265 (1998), hep-ex/9806008. F. Ellinghaus, R. Shanidze and J. Volmer (HERMES Collab oration), Deeply-virtual Compton scattering on deuterium and neon at HERMES, AIP Conf. Pro c. 675, 303­307 (2003), hep-ex/0212019.

[GS03] [H101]

[HERMES98] [HERMES03]

[HERMES04a] E. Avetisyan, A. Rostomyan and A. Ivanilov (HERMES Collab oration), Beam-spin azimuthal asymmetries in pion electropro duction at HERMES, (2004), hep-ex/0408002. [HERMES04b] B. Seitz (HERMES Collab oration), The HERMES recoil detector: A combined silicon strip and scintillating fibre detector for tracking and particle identification, Nucl. Instrum. Meth. A535, 538­541 (2004). [HERMES05] A. Airap etian et al. (HERMES Collab oration), Quark helicity distributions in the nucleon for up, down, and strange quarks from semiinclusive deep-inelastic scattering, Phys. Rev. D71, 012003 (2005), hep-ex/0407032. P. Ho o dbhoy and X.-D. Ji, Helicity-flip off-forward parton distributions of the nucleon, Phys. Rev. D58, 054006 (1998), hep-ph/9801369. R. Hofstadter, Nuclear and nucleon scattering of high-energy electrons, Ann. Rev. Nucl. 7, 231­316 (1957). P. R. B. Saull (H1 and ZEUS Collab orations), Deeply virtual Compton scattering, Prepared for International Europhysics Conference on HighEnergy Physics (HEP 2001), Budap est, Hungary, 12-18 Jul 2001. J.C.Collins and A. Freund, Phys. Rev. D 59, 074009 (1999). X.-D. Ji, Deeply-virtual Compton scattering, Phys. Rev. D55, 7114­ 7125 (1997), hep-ph/9609381. X.-D. Ji, Gauge invariant decomp osition of nucleon spin, Phys. Rev. Lett. 78, 610­613 (1997), hep-ph/9603249.

[HJ98] [Hof57] [HZ]

[JF99] [Ji97a] [Ji97b]

88


[Ji98] [Ji03] [Kra] [LT05a] [LT05b] [M+ 94]

X.-D. Ji, Off-forward parton distributions, J. Phys. G24, 1181­1205 (1998), hep-ph/9807358. X.-d. Ji, Generalized parton distributions and the spin structure of the nucleon, Nucl. Phys. Pro c. Suppl. 119, 41­49 (2003), hep-lat/0211016. B. Krauss, Deeply virtual Compton scattering and the HERMES recoildetector, DESY-THESIS-2005-008. S. Liuti and S. K. Taneja, Microscopic description of deeply virtual Compton scattering off spin-0 nuclei, (2005), hep-ph/0505123. S. Liuti and S. K. Taneja, Nuclear medium mo difications of hadrons from generalized parton distributions, (2005), hep-ph/0504027. D. Muller et al., Wave functions, evolution equations and evolution kernels from light-ray op erators of QCD, Fortschr. Phys. 42, 101 (1994), hep-ph/9812448. L. Mankiewicz, G. Piller and T. Weigl, Hard exclusive meson pro duction and nonforward parton distributions, Eur. Phys. J. C5, 119­128 (1998), hep-ph/9711227. L. W. Mo and Y.-S. Tsai, RADIATIVE CORRECTIONS TO ELASTIC AND INELASTIC e p AND mu p SCATTERING, Rev. Mo d. Phys. 41, 205­235 (1969). J. Ashman et al. (Europ ean Muon Collab oration), MEASUREMENT OF THE RATIOS OF DEEP INELASTIC MUON - NUCLEUS CROSSSECTIONS ON VARIOUS NUCLEI COMPARED TO DEUTERIUM, Phys. Lett. B202, 603 (1988). J. Ashman et al. (Europ ean Muon Collab oration), An investigation of the spin structure of the proton in deep inelastic scattering of p olarized muons on p olarized protons, Nucl. Phys. B328, 1 (1989). W.-D. Nowak, Physics ob jectives for future studies of the spin structure of the nucleon, (2002), hep-ph/0210409. W.-D. Nowak, Deeply virtual Compton scattering: Results and future, (2005), hep-ex/0503010. M. V. Polyakov, Generalized parton distributions and strong forces inside nucleons and nuclei, Phys. Lett. B555, 57­62 (2003), hepph/0210165. M. V. Polyakov and C. Weiss, Skewed and double distributions in pion and nucleon, Phys. Rev. D60, 114017 (1999), hep-ph/9902451.

[MPW98]

[MT69]

[Muon88]

[Muon89]

[Now02] [Now05] [Pol03]

[PW99]

89


[Rad96] [Rad97] [Sjo94] [ST64]

A. V. Radyushkin, Asymmetric gluon distributions and hard diffractive electropro duction, Phys. Lett. B385, 333­342 (1996), hep-ph/9605431. A. V. Radyushkin, Nonforward parton distributions, Phys. Rev. D56, 5524­5557 (1997), hep-ph/9704207. T. Sjostrand, High-energy physics event generation with PYTHIA 5.7 and JETSET 7.4, Comput. Phys. Commun. 82, 74­90 (1994). A. A. Sokolov and I. M. Ternov, On Polarization and spin effects in the theory of synchrotron radiation, Phys. Dokl. 8, 1203­1205 (1964).

90


Acknowledgements
This thesis could not have b een written without the help of many p eople and I would like to thank all of them for their supp ort. First and foremost, I would like to thank my advisor Dr. Wolf-Dieter Nowak who encouraged me to b egin this work at HERMES and always found the time and patience to discuss all questions and to listen to my presentations. I am also grateful for his careful reading of this manuscript and his detailed comments and suggestions on it. My sp ecial thanks go to Prof. Dr. Butz who agreed on sup ervising this dissertation as my official advisor at the University of Leipzig. I want to express my gratitude to all Hermes memb ers, esp ecially to Elke Aschenauer and Delia Hasch, for their supp ort throughout the work on the diploma thesis. In particular I would like to thank Frank Ellinghaus, who intro duced me in the DVCS analysis, Hayg Guler and Frank Zhenyu Ye for their helps on my data analysis and for all the inspiring discussions on physics as well as on unrelated topics. Outside the DVCS group I b enefited from discussions with Achim Hillenbrand and Markus Diefenthaler, who help ed me with their patient advices and with software problems of all kind. I would like to thank Jim Stewart and Eduard Avetisyan for their knowledge and explanation of the Hermes detectors. My thanks also go to all Hermes memb ers who are not mentioned ab ove and I very much appreciated the friendly athmosphere in the entire group, who made work and life in Hamburg so enjoyable. Finally, I would like to thank my friends and my family for their supp ort and understanding.

91


Declaration
This do cument is the result of my own work. Material from the published or unpublished work of otheres, which is referred to in the do cument, is credited to the author in the text.

Selbst¨ andigkeitserkl¨ arung
Hiermit erkl¨ ich, die vorliegende Diplomarb eit selbst¨ are andig und nur unter Verwendung der angegeb enen Hilfsmittel und Quellen angefertigt zu hab en.

Leipzig, den 16.08.2005

92