## Jefferson Lab PAC 47

## Neutron DVCS Measurements with BONuS12 in CLAS12

M. Hattawy †‡ , M. Amaryan, S. B¨ltmann, G. Dodge, N. Dzbenski, C. Hyde, S. Kuhn , u † D. Payette, J. Poudel, L. Weinstein

Old Dominion University, Norfolk, VA 23529, USA

R. Dupr´ e , M. Guidal, D. Marchand, C. Mu˜oz, S. Niccolai, E. Voutier † n Institut de Physique Nucl´aire, CNRS-IN2P3, Univ. Paris-Sud, Universit´ Paris-Saclay, e e 91406 Orsay Cedex, France

K. Hafidi, Z. Yi

Argonne National Laboratory, Lemont, IL 60439, USA

T. Chetry, L. El-Fassi

Mississippi State University, Mississippi State, MS 39762, USA

N. Baltzell, G. Gavalian, F. X. Girod, S. Stepanyan Thomas Jefferson National Accelerator Facility, Newport News, VA 23606, USA

I. Albayrak, E. Christy, A. Nadeeshani Hampton University, Hampton, VA 23669, USA

## N. Kalantarians

Virginia Union University, 1500 N. Lombardy St., Richmond, VA 23220, USA

C. Ayerbe Gayoso College of William and Mary, Williamsburg, VA 23187, USA

## D. Jenkins

Virginia Tech, Blacksburg, VA 24061, USA

## A CLAS12 Run-Group Addition Proposal

‡ Contact person: Hattawy@jlab.org

†

Spokesperson

## Abstract

The three-dimensional picture of quarks and gluons in the nucleon is set to be revealed through deeply virtual Compton scattering (DVCS). With the absence of a free neutron target, the deuterium target represents the simplest nucleus to be used to probe the internal 3D partonic structure of the neutron. We propose here to measure the beam spin asymmetry (BSA) in incoherent neutron DVCS together with the approved E12-06-113 experiment (BONuS12) within the run group F, using the same beam time, simply with addition of beam polarization. The DVCS BSA on the quasi-free neutron will be measured in a wide range of kinematics by tagging the scattered electron and the real photon final state with the spectator proton. We will also measure BSA with all final state particles detected including the struck neutron. The proposed measurements is complementary to the approved CLAS12 experiment E12-11-003, which will also measure the quasi-free neutron DVCS by detecting the scattered neutron, but not the spectator proton. Indeed, besides providing more data for neutron DVCS, this experiment will allow a comparison of the measurement of the BSA of neutron DVCS from the approved E12-11-003 with the measurements using the two methods proposed herein. This comparison will help to understand the impact of nuclear effects, such as the final state interactions (FSI) and Fermi motion on the measurement of the neutron DVCS.

## Contents

| Abstract     | Abstract                                                                                      | Abstract                                                                                                                | Abstract                                                                                                                | 5     |
|--------------|-----------------------------------------------------------------------------------------------|-------------------------------------------------------------------------------------------------------------------------|-------------------------------------------------------------------------------------------------------------------------|-------|
| Introduction | Introduction                                                                                  | Introduction                                                                                                            | Introduction                                                                                                            | 7     |
| 1            |                                                                                               | Neutron                                                                                                                 | Partonic Structure                                                                                                      | 9     |
|              | 1.1                                                                                           | DVCS Formalism and Observables . . . . . . . . . .                                                                      | DVCS Formalism and Observables . . . . . . . . . .                                                                      | 9     |
|              | 1.2                                                                                           | Neutron GPDs . . . . . . . . . . . . . . . . . . . . . . . . Measuring the neutron DVCS . . . . . . . . . . . . . . . . | Neutron GPDs . . . . . . . . . . . . . . . . . . . . . . . . Measuring the neutron DVCS . . . . . . . . . . . . . . . . | 11 12 |
|              | 1.3                                                                                           | 1.3.1                                                                                                                   | Two new methods for two objectives . . . . .                                                                            | 12    |
|              |                                                                                               | 1.3.2                                                                                                                   | Proton-Tagged Neutron DVCS with BONuS12                                                                                 | 13    |
|              |                                                                                               | 1.3.3                                                                                                                   | Fully exclusive neutron DVCS with BONuS12                                                                               | 16    |
| 2            | Run conditions and Experimental setup . . . . . . . . . . . . . . . . . . . . . . . . . . . . | Run conditions and Experimental setup . . . . . . . . . . . . . . . . . . . . . . . . . . . .                           | Run conditions and Experimental setup . . . . . . . . . . . . . . . . . . . . . . . . . . . .                           | 17    |
|              | 2.1                                                                                           | Run Group F conditions .                                                                                                | Run Group F conditions .                                                                                                | 17    |
|              | 2.2                                                                                           | The CLAS12 Spectrometer . . . . . . . . . . . . . . . . . .                                                             | . . . . . . .                                                                                                           | 18    |
|              | 2.3                                                                                           | BONuS12 RTPC . . .                                                                                                      | . . . . . . .                                                                                                           | 18    |
|              | 2.4                                                                                           | Beam Polarization . .                                                                                                   | . . . . . . . . . . . . . . . . .                                                                                       | 21    |
| 3            |                                                                                               | Projections                                                                                                             | for the Proposed Measurements                                                                                           | 22    |
|              | 3.1                                                                                           | Monte-Carlo Simulation . . .                                                                                            | . . . . . . . . . . . . .                                                                                               | 22    |
|              | 3.2                                                                                           | Particle Identification . . . . . . . . . . . . . . . . . . .                                                           | . . . . . .                                                                                                             | 23    |
|              |                                                                                               | 3.2.1                                                                                                                   | Electron Identification . . . . . .                                                                                     | 23    |
|              |                                                                                               | 3.2.2                                                                                                                   | Proton Identification . . . . . . . . . . . . . .                                                                       | 23    |
|              |                                                                                               | 3.2.3                                                                                                                   | Neutrals Identification . . . . . . . . . . . . . . .                                                                   | 24    |
|              | 3.3                                                                                           | Beam-Spin Asymmetry . . .                                                                                               | . . . . . . . . . . . . .                                                                                               | 26    |
|              | 3.4                                                                                           | Projections . . . . . . .                                                                                               | . . . . . . . . . . . . . . .                                                                                           | 27    |
|              |                                                                                               | 3.4.1                                                                                                                   | Proton-Tagged neutron DVCS Projections . . . . .                                                                        | 27    |
|              |                                                                                               | 3.4.2                                                                                                                   | Fully exclusive n-DVCS projections . . .                                                                                | 28    |
| 4            | Additional Physics Measurements with Run Group F                                              | Additional Physics Measurements with Run Group F                                                                        | Additional Physics Measurements with Run Group F                                                                        | 36    |
| Summary      | Summary                                                                                       | Summary                                                                                                                 | Summary                                                                                                                 | 37    |

## Introduction

Inclusive deep inelastic scattering (DIS) experiments have been instrumental in advancing our understanding of the QCD structure of nucleons in the past. More recently, hard exclusive experiments such as Deeply Virtual Compton Scattering (DVCS) and Deeply Virtual Meson Production (DVMP) have provided important new probes that allows us to explore both the longitudinal motion of partons and their transverse spatial structure in nucleons through the generalized parton distribution (GPD) framework. The GPDs correspond to the coherence between quantum states of different (or same) helicity, longitudinal momentum, and transverse position. In an impact parameter space, they can be interpreted as a distribution in the transverse plane of partons carrying a certain longitudinal momentum [1-3]. A crucial feature of GPDs is the access to the transverse position of partons which, combined with their longitudinal momentum, leads to the total angular momentum of partons [4]. This information is not accessible to inclusive DIS which measures probability amplitudes in the longitudinal plane.

Recent results from CLAS on incoherent DVCS from a 4 He target have shown that even in light nuclei, large nuclear effects can be observed [7]. As this recent measurement shows a deviation of ≈ 30% for the beam spin asymmetry (BSA) of DVCS of bound protons in helium, one can assume that a sizable correction will be necessary in the neutron case as well; although reduced in the much less tightly bound deuterium. This could be particularly dramatic as the BSA expected on neutrons are rather small, making the nuclear suppression observed in helium of a similar size as the signal expected for the neutron. We propose here to transpose the methods developed for BONuS12 for PDF measurements [8] to the GPDs. In particular, we are looking how the initial state and the break-up of the deuterium might affect the measured BSA. The upcoming Run group F (solely composed of the BONuS12

A high luminosity facility such as Jefferson Lab offers a unique opportunity to map out the three-dimensional quark and gluon structure of nucleons and nuclei. This is one of the flagship of the JLab 12 GeV scientific program, such that many approved proposals to the JLab Program Advisory Committee (PAC) have focused on studies of the 3D proton or neutron structure. We propose here to extend the DVCS measurements of neutrons to the effective quasi-free neutrons obtained using the BONuS12 experiment. Measurements on the neutron are critical in the Jefferson Lab GPD program for two reasons, first as neutrons offer the only access to flavor decomposition for GPDs, second because of the importance of the E GPD to measure the Ji sum rule [5,6].

experiment) scheduled for February-April 2020, is the perfect occasion to gather data on this topic, which can be done with the simple addition of beam polarization to the setting of the run group!

## Chapter 1

## Neutron Partonic Structure

## 1.1 DVCS Formalism and Observables

A wealth of information on the structure of hadrons lies in the correlations between the momentum and spatial degrees of freedom of the partons. These correlations can be revealed through deeply virtual Compton scattering (DVCS), i.e., the hard exclusive lepto-production of a real photon, which provides access to a three-dimensional (3-D) imaging of partons within the generalized parton distributions (GPDs) framework [5, 9-12]. The cross section for DVCS on a spin-1/2 target can be parametrized in terms of four helicity conserving GPDs: H q , E q , ˜ H q , and ˜ E q . The GPDs H E H , , ˜ and ˜ E are defined for each quark flavor (q = u, d, s, ... ). Analogous GPDs exist for the gluons, see references [12, 13] for further details. In this work, we are mostly concerned by the valence quark region, in which the sea quarks and the gluons contributions do not dominate the DVCS scattering amplitude. The GPDs H and ˜ H conserve the spin of the nucleon, while E and ˜ E flip it [14]. The H and E GPDs are called the unpolarized GPDs as they represent the sum over the different configurations of the quarks' helicities, whereas ˜ H and ˜ E are called the polarized GPDs because they are made up of the difference between the orientations of the quarks' helicities.

The differential cross section of leptoproduction of photons for a longitudinally-polarized electron beam and an unpolarized nucleon target can be written as:

$$\frac { d \sigma } { d x _ { B } \, d y \, d t \, d \phi \, d \varphi } = \frac { \alpha ^ { 3 } x _ { B } y } { 1 6 \pi ^ { 2 } Q ^ { 2 } \sqrt { 1 + \varepsilon ^ { 2 } } } \left | \frac { \mathcal { T } } { e ^ { 3 } } \right | ^ { 2 }$$

where ε ≡ 2 x B Q M n , x B = Q / 2 (2 p 1 · q 1 ) is the Bjorken variable, y = ( p 1 · q 1 ) / p ( 1 · k 1 ) is the photon energy fraction, φ is the angle between the leptonic and hadronic planes, ϕ is the scattered electron's azimuthal angle, Q 2 = -q 2 1 , and q 1 = k 1 -k 2 . The particle momentum definitions are shown in Figure 1.1. The momentum transfer where the nucleon is initially at rest, ∆ = p 1 -p 2 and t = ∆ . The Bjorken variable is related to another scaling variable 2

Figure 1.1: Leading-twist DVCS handbag diagram with the momentum definitions labeled.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-582c36cf-e156-4fba-ba34-ae2b1558467a.jpg)

called skewedness:

$$\xi = \frac { x _ { B } } { 2 - x _ { B } } + \mathcal { O } ( 1 / Q ^ { 2 } ).$$

The amplitude is the sum of the DVCS, the Bethe-Heitler (BH), and the interference amplitudes, and when squared has terms

$$\mathcal { T } ^ { 2 } = | \mathcal { T } _ { \text{BH} } | ^ { 2 } + | \mathcal { T } _ { \text{DVCS} } | ^ { 2 } + \mathcal { I }$$

where the first is the BH contribution, the second is the DVCS part, and the last term is the interference part,

$$\mathcal { I } = \mathcal { T } _ { \text{DVCS} } \mathcal { T } _ { \text{BH} } ^ { * } + \mathcal { T } _ { \text{DVCS} } ^ { * } \mathcal { T } _ { \text{BH} }.$$

The corresponding amplitudes are calculated with the diagrams shown in Figures 1.1 and 1.2. The details of contracting the DVCS tensor with various currents and tensors can be found in [15]. The resulting expressions for the amplitudes are

Figure 1.2: BH handbag diagrams.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-f4180ef8-9ae3-49e4-9b77-8d857fbda72a.jpg)

$$| \mathcal { T } _ { \text{BH} } | ^ { 2 } = \frac { e ^ { 6 } ( 1 + \varepsilon ^ { 2 } ) ^ { - 2 } } { x _ { B } ^ { 2 } y ^ { 2 } \, t \, \mathcal { P } _ { 1 } ( \phi ) \mathcal { P } _ { 2 } ( \phi ) } \left \{ c _ { 0 } ^ { \text{BH} } + \sum _ { n = 1 } ^ { 2 } \left [ c _ { n } ^ { \text{BH} } \cos ( n \phi ) + s _ { n } ^ { \text{BH} } \sin ( n \phi ) \right ] \right \} \quad \quad ( 1. 5 )$$

$$| \mathcal { T } _ { D V C S } | ^ { 2 } = \frac { e ^ { 6 } } { y ^ { 2 } Q ^ { 2 } } \left \{ c _ { 0 } ^ { D V C S } + \sum _ { n = 1 } ^ { 2 } \left [ c _ { n } ^ { D V C S } \cos ( n \phi ) + s _ { n } ^ { D V C S } \sin ( n \phi ) \right ] \right \} \, _ { \cdot } \,.$$

$$\mathcal { I } = \frac { \sigma \, \bar { \nu } } { e ^ { 6 } ( 1 + \varepsilon ^ { 2 } ) ^ { - 2 } } \left \{ c _ { 0 } ^ { \mathcal { I } } + \sum _ { n = 1 } ^ { 3 } \left [ c _ { n } ^ { \mathcal { I } } \cos ( n \phi ) + s _ { n } ^ { \mathcal { I } } \sin ( n \phi ) \right ] \right \}$$

The functions c 0 , c n , and s n are called Fourier coefficients . They depend on the kinematic variables and the operator decomposition of the DVCS tensor for a target with a given spin. At leading twist there is a straightforward form factor decomposition which relates the vector and axial-vector operators with the so-called Compton form factors (CFFs) [16]. The Compton form factors appearing in the DVCS amplitudes are integrals of the type

$$\mathcal { F } = \int _ { - 1 } ^ { 1 } d x F ( \mp x, \xi, t ) C ^ { \pm } ( x, \xi )$$

where the coefficient functions at leading order take the form

$$C ^ { \pm } ( x, \xi ) = \frac { 1 } { x - \xi + i \varepsilon } \pm \frac { 1 } { x + \xi - i \varepsilon }.$$

We plan on measuring the beam spin asymmetry as a function of φ

$$A _ { L U } ( \phi ) = \frac { d \sigma ^ { \uparrow } ( \phi ) - d \sigma ^ { \downarrow } ( \phi ) } { d \sigma ^ { \uparrow } ( \phi ) + d \sigma ^ { \downarrow } ( \phi ) }$$

where the arrows indicate the electron beam helicity.

## 1.2 Neutron GPDs

The measurement of free proton DVCS has been the focus of a worldwide effort [1725] involving several accelerator facilities such as Jefferson Lab, DESY and CERN. These measurements now enable the extractions of GPDs and a 3-D tomography of the free proton [26, 27]. The aim of this proposal is to enhance the neutron GPD measurements along the approved CLAS12 experiment E12-11-003 [28], which will also measure the quasi-free neutron DVCS by detecting the scattered neutron in deuterium.

In the fits of PDFs, for example [29] , neutrino and deuterium data allow to make a flavor dependent extraction, this option is not available for GPDs yet due to the lack of reliable data. Indeed, the observables of DVCS, such as cross sections and beam-spin asymmetries are much smaller on neutron targets, while the nuclear effects in deuterium increase uncertainties. These issues have lead to results [30] which are not precise enough to help in a flavor dependent GPD extraction. To achieve this performance, one will need a large quantity of

high precision data, a goal set by E-12-11-003. However, the impact of the uncertain initial state and the final state interactions on the integration of these data in global fits remain unclear. With this proposal, we propose to both provide more data on the neutron, which is always important, but most importantly these data will have completely independent systematic uncertainties from the E-12-11-003 data. That will make them complementary and indicate in what ways these methods are equivalent or potentially need corrections.

Neutron DVCS is also hoped to provide an important contribution to the extraction of the GPD E [31]. The reasoning behind this expectation is as follows: the GPD E never appears to be dominant in the usually measured DVCS observables, so as a sub-leading contribution it is always affected with large error bars. Actually, recent extractions [32, 33] show that we still barely have any constraint on E using all the world proton data. However, as form factors often appear in the expressions of DVCS observables because of the Bethe Heitler process (see below) the situation is very different for protons and neutrons. Indeed, the F 1 form factor of the neutron is very small, making the E GPD more prominent in some of the neutron DVCS observables. Most notably the sin component of the beam spin asymmetry is defined as:

$$A _ { L U } ^ { \sin \phi } = \frac { 1 } { \pi } \int _ { \pi } ^ { \pi } d \phi \sin \phi A _ { L U } ( \phi ) & & ( 1. 1 1 )$$

is proportional to the following combination of Compton form factors [26]

$$A _ { L U } ^ { \sin \phi } \, \infty \, \Im m \left [ F _ { 1 } \mathcal { H } + \xi ( F _ { 1 } + F _ { 2 } ) \tilde { \mathcal { H } } - \frac { t } { 4 M ^ { 2 } } F _ { 2 } \mathcal { E } \right ],$$

highlighting the effect of a suppressed F 1 on the main term. This lead to the idea that neutron DVCS will also help significantly to constrain the GPD E . This goal is of course driven by the long standing objective of GPD physics to measure the Ji sum rule in the nucleon:

$$J _ { q } \, = \, \frac { 1 } { 2 } \, \int _ { - 1 } ^ { + 1 } d x \, x \, \left [ H ^ { q } ( x, \xi, t = 0 ) + E ^ { q } ( x, \xi, t = 0 ) \right ],$$

which links the total angular momentum ( J q ) carried by each quark q to the sum of the second moments over x of the GPDs H and E , that will complete the decomposition of the nucleon spin in its various components [5,6].

## 1.3 Measuring the neutron DVCS

## 1.3.1 Two new methods for two objectives

As shown in Fig. 1.3, one can measure all final state particles of the DVCS reaction on deuterium using CLAS12 and Bonus12. This method, while perfect on paper, leads to an efficiency problem as both the measurement of the spectator proton and of the neutrons are challenging and have low efficiency (intrinsically for the neutron detection and by the limitation of phase space available for the proton). However, it is possible to ensure the

Figure 1.3: Fully exclusive neutron DVCS diagram in deuterium.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-cb0558dc-2a45-45c3-9461-9443e056a7f7.jpg)

exclusivity of a process when missing one of the final state particles by applying cuts on missing mass, momentum and energy. This strategy can be used, either to leave the spectator proton or the neutron undetected. The former approach is being used in the ongoing E1211-003. We propose here to add the other option, detecting the spectator proton but not the neutron, in order to increase the amount of available data, but most importantly, to confirm the equivalency of both methods. Moreover, we will be able to perform the fully exclusive measurement, but with rather limited statistics. We intend to use this over-constrained last measurement to study the systematic effects linked to Fermi motion and final state interaction effects on our DVCS observable of interest, i.e., beam-spin asymmetry.

## 1.3.2 Proton-Tagged Neutron DVCS with BONuS12

The proton-tagged neutron detection scheme is described in Fig. 1.4, where we see that a detector for low energy proton spectators is necessary. Here, we propose to use the Bonus radial TPC, which is designed to make a similar type of measurement for inclusive deeply inelastic scattering. The first goal of the present proposed experiment is simply to provide more data in the field of neutron GPD. Indeed, the measurement of neutron DVCS is very challenging and very little published data are available at this point.

In order to resolve the initial state issue, we evaluate the standard Lorentz invariant x

Second, we observe that the systematics from this measurement are going to be mostly independent of the ones from E12-11-003. Indeed, while this measurement will be missing one of the high energy product of the reaction leading to more uncertainty in exclusivity cuts, it will detect the spectator proton, significantly helping in the understanding of the nuclear effects. This has interestingly an impact on both initial state and final state effects. In the initial state, the neutron is in fact not at rest and carries some Fermi momentum, which can be directly inferred from the kinematic of the spectator proton, thanks to the simple two body deuterium. On the final state side, detecting the spectator proton in a certain range of momentum and angle allows to significantly reduce the probability of final state interactions to have occurred.

Figure 1.4: Proton-tagged neutron DVCS diagram in deuterium.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-012d01a6-5ae1-468b-996f-a89272e4849e.jpg)

and the γ n ∗ invariant mass ( W ) in terms of the spectator kinematics. Both invariants, x and W , acquire a star to indicate true invariants rather than the values calculated assuming a stationary, on-shell target:

$$x ^ { * } = \frac { Q ^ { 2 } } { 2 M _ { N } E y ( 2 - \alpha _ { s p } ) } = \frac { x _ { B } } { 2 - \alpha _ { s p } },$$

$$W ^ { * } = M ^ { * 2 } - Q ^ { 2 } + 2 M E y ( 2 - \alpha _ { s p } ),$$

where α sp = E s -p z s M N , with M N , E s , and p z s are the on-shell mass, energy, and z-momentum component of the spectator proton. The off-shell mass of the bound nucleon is given by M ∗ 2 = ( M D -E s ) 2 -p 2 s , where M D is the rest mass of the deuterium.

It is important to note that the detection of the neutron or the spectator proton are not equivalent. While it could appear to be so after applying the exclusivity cut, it is not the case because of the large uncertainty (a percent or more) in the energy measurement of photons. This uncertainty is larger than the momentum of the spectator proton and therefore completely hinders our capability to reconstruct it from missing momentum and energy methods.

To understand the regions where final state interactions are expected to be significant, we look at the spectator momentum and angle relative to the momentum transfer, θ s . In Figure 1.5 [34, 35], we see calculations for the inclusive case. At low recoil momentum and backwards spectator angle, the FSIs are negligible, whereas at high momenta perpendicular to the momentum transfer, FSIs are maximized.

One should mention that our recent measurement [7] of incoherent DVCS on helium conducted during the 6 GeV era (E08-024) have shown significant modification of the proton beam spin asymmetry in 4 He. These results are shown in Figure 1.6, where we observe a much smaller asymmetry for bound protons. However, to this day, it is not possible to say if this suppression should be associated to the parasitic nuclear effects presented above or to a more fundamental source linked to the EMC effect. While not directly comparable,

Figure 1.5: Ratio of cross sections for the FSI model from [34] to Plane Wave Impulse Approximation (PWIA) calculation as a function of the spectator momentum (left) and spectator angle (right). The caption on the left figure has been edited from the original paper to correct for a typo.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-918cb42f-64ce-420d-9bee-fdfb00be9f22.jpg)

Figure 1.6: The A LU ratio of the bound to the free proton at φ = 90 ◦ as a function of t . The black squares are from CLAS-eg6 experiment [7], the green circle is the HERMES measurement [36]. The error bars represent the statistical uncertainties, while the gray band represents the systematic uncertainties. The blue and red curves are results of off-shell calculations [37]. The solid and dashed black curves are from on-shell calculations [38].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-0d548d56-2787-4b2d-96eb-a5e892095649.jpg)

this unexpectedly large deviation seems to hint to large nuclear effects in incoherent nuclear DVCS that have not been explored at all on the theory side to this day.

In conclusion, the measurement of this channel will offer a slightly smaller amount of data as other CLAS12 measurements. However, the key lies in the different and independent systematic effects of this new measurement method. In particular, we will study and reduce the impact of nuclear effects, which will significantly help in the interpretation of E12-11-003 data [28].

## 1.3.3 Fully exclusive neutron DVCS with BONuS12

As highlighted above, the fully exclusive measurement of neutron DVCS is difficult since two hard-to-detect particles have to be measured, a slow proton and a fast neutron. However, it can still provide relevant information in a limited kinematic space. As this measurement is fully exclusive, it provides over-constraints by using missing mass, momentum and energy cuts. This will give us a particularly clean sample of data with a minimal amount of corrections to be applied and the best control over systematic errors. We propose to use this sample, that can only be obtained using a recoil detector, to study the effects of Fermi motion, final state interactions and incomplete detection of the final state described above. This will confirm assumptions made in other neutron DVCS measurements and help understand better their systematic errors.

## Chapter 2

## Run conditions and Experimental setup

## 2.1 Run Group F conditions

Run Group F (E12-06-113) has been approved to collect 35 PAC days (100% efficiency) of data on deuterium with 11 GeV electron beam and another six days on hydrogen. One of the days of hydrogen data taking will be carried out with a low energy electron beam of about 2.2 GeV. The 40 cm long target filled with 7 atm deuterium gas at room temperature and the 200 nA electron beam will yield a combined nucleon luminosity of about 2 × 10 34 cm -2 s -1 , about a factor of five below the standard CLAS12 nominal luminosity.

The Run Group F approved measurement (BONuS12) [8], as well as the proposed measurements here, will use the same trigger as in RG-B [28], that is an electron in Forward Tracker with road matching.

For the detection of the recoil protons, Run Group F is going to install a new and enlarged radial time projection chamber (RTPC) and target gas cell assembly, very similar to the ones used by the BONuS (E03-012) and eg6 (E08-024) experiments. The RTPC can detect proton recoil momenta down to a lower limit of 70 MeV/c while being insensitive to minimum ionizing particles. The RTPC will be replacing the central detector's silicon tracker and barrel micromegas, while keeping an updated design of the forward micromegas (FMT). In the updated version of the FMT, only three layers of micromegas will be kept to improve the electron's reconstructed vertex resolution while reducing the material in the path of the electrons. In the following sections we briefly introduce the experimental setup of Run Group F.

## 2.2 The CLAS12 Spectrometer

The CLAS12 spectrometer is designed to operate with 11 GeV beam at an electronnucleon luminosity of L = 1 × 10 35 cm -2 s -1 . The baseline configuration of the CLAS12 detector consists of the forward detector and the central detector packages [39] (see Figure 2.1). The CLAS12 Central Detector [39] is designed to detect various charged particles over a wide momentum and angular range. The main detector package includes:

- · Solenoid Magnet: provides a central longitudinal magnetic field up to 5 Tesla, which serves to curl emitted low energy Møller electrons and determine particle momenta through tracking in the central detector.
- · Central Tracker: consists of 3 double layers of silicon strips and 6 layers of Micromegas. The thickness of a single silicon layer is 320 µm.
- · Central Time-of-Flight (CTOF): an array of scintillator paddles with a cylindrical geometry of radius 26 cm and length 50 cm; the thickness of the detector is 2 cm with designed timing resolution of σ t = 50 ps, used to separate pions and protons up to 1.2 GeV/ . c

We will use the forward detector for electron, photon, and neutron detection. The central detector's silicon tracker and barrel micromegas will be removed to leave room for the BONuS12 RTPC for the detection of the slow recoiling protons in deuterium. In addition to the main CTOF of CLAS12 Central detector, we will be using the Central Neutron Detector (CND) for the detection of the final state recoiling neutrons.

The scattered electron, photon, and some neutrons will be detected in the forward detector which consists of the High Threshold Cherenkov Counters (HTCC), Drift Chambers (DC), the Low Threshold Cherenkov Counters (LTCC), the Time-of-Flight scintillators (TOF), the Forward Calorimeter and the Preshower Calorimeter. The charged particle identification in the forward detector is achieved by utilizing the combination of the HTCC and TOF arrays with the tracking information from the Drift Chambers. The HTCC together with the Forward Calorimeter and the Preshower Calorimeter will provide a pion rejection factor of more than 2000 up to a momentum of 4.9 GeV/c, and a rejection factor of 100 above 4.9 GeV/c. The photons and the neutrons are detected using the calorimeters. As will be showing later on, the majority of the final state recoiling neutrons will be detected using the sub-detectors of CLAS12 Central Detector, in particular CTOF and CND.

## 2.3 BONuS12 RTPC

The new CLAS12 RTPC (BONuS12) is 400 mm long cylinder of 160 mm diameter. The electric field is directed perpendicularly to the beam direction, such that drifting electrons are pushed away from the beam line. These electrons are amplified by three layers of cylindrical gas electron multipliers (GEM) and detected by the readout system on the external shell of

Figure 2.1: (Top) The schematic layout of the CLAS12 baseline design with BONuS12 RTPC replacing the silicon tracker and the barrel micromegas. (Bottom) Schematic layout showing BONuS12 RTPC with the modified design of the forward micromegas.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-090dc8ae-2a0a-46e8-af16-ac9c67174bdb.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-ae9e7a16-0129-4461-98f2-8921aff63376.jpg)

Figure 2.2: (Left) Schematic layout showing BONuS12 RTPC showing the readout padboard and few adaptor boards in addition to the gas lines. (Right) Schematic drawing of the CLAS RTPC in a plane perpendicular to the beam direction. See text for description of the elements.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-2ef49d9d-e074-4c63-821f-ed6fc2c4f0cd.jpg)

the detector as illustrated in Figure 2.2. The BONuS12 RTPC covers almost 100% of the azimuthal angle range.

We detail here the different regions shown in Figure 2.2 starting from the beam line towards larger radius:

- · The 7 atm Deuterium gas target extends along the beamline forming the detector central axis. It is a 6 mm diameter Kapton straw with a 50 µ m wall of 492 mm length such that its entrance and exit 15 µ m aluminum windows are placed outside of the detector volume. The detector and the target are placed in the center of the solenoid.
- · The first gas gap covers the radial range from 3 mm to 20 mm. It is filled with 4 He gas at 1 atm to minimize secondary interactions from Møller electrons scattered by the beam. This region is surrounded by a 4 µ mthick window made of grounded aluminized Mylar.
- · The second gas gap region extends between 20 mm and 30 mm and is filled with the gas mixture of 80% 4 He and 20% CO . 2 This region is surrounded by a 4 µ m thick window made of aluminized Mylar set at -4260 V to serve as the cathode.
- · The drift region is filled with the same 4 He-CO 2 gas mixture and extends from the cathode to the first GEM, 70 mm away from the beam axis. The electric field in this region is perpendicular to the beam and averages around 550 V/cm.

- · The electron amplification system is composed of three GEMs located at radii of 70, 73 and 76 mm. The first GEM layer is set to ∆ V = -1620 V relative to the ground and then each subsequent layer is set to a lower voltage relative to the previous to obtain a strong ( ∼ 1600 V/cm) electric field between the GEM foils. A 275 V bias is applied across each GEM for amplification.
- · The readout board has an internal radius of 79 mm and collects charges after they have been multiplied by the GEMs. Adaptor circuit boards are plugged directly on its outer side and transmit the signal to the Hitachi cables connected to the BMT DAQ electronics.

## 2.4 Beam Polarization

For our proposed measurements we ask for the electron beam to be highly polarized, i.e., 85% longitudinally polarized beam, which is the average achieved polarization of the beam during Run Group A and Run Group B using the upgraded 12 GeV experimental setup. Regarding measuring the beam polarization, no additional commissioning time is required. The spokespersons of Run Group F [40] has accepted to schedule time to perform Møller runs once or twice a week, which is enough for the proposed measurements.

## Chapter 3

## Projections for the Proposed Measurements

We propose to measure the neutron's DVCS beam-spin asymmetry in the following two channels:

- 1. Proton-tagged neutron DVCS: D + γ ∗ -→ γ + ( n ) + p
- 2. Fully exclusive n-DVCS: D + γ ∗ -→ γ + n + p

To demonstrate the experimental feasibility and to extract projections for our proposed measurements, we present here a simulation study using a realistic n-DVCS event generator and the official simulation-reconstruction chain of the CLAS12 spectrometer augmented with BONuS12 RTPC.

## 3.1 Monte-Carlo Simulation

An event generator for DVCS/BH and exclusive π 0 electroproduction on the neutron inside a deuterium target has been developed [41]. The DVCS amplitude is calculated according to the BKM formalism [15], while the GPDs have been taken from the VGG model [42,43]. The Fermi-motion distribution is calculated with the Paris potential [44].

The output of the event generator was fed through the CLAS12 official simulation (GEMC 4.3.0) [45] and reconstruction (COATJAVA 5b.7.4) [46] chain, to simulate the detectors' acceptance and resolutions for the following final state particles, electrons, photons, and neutrons, within the proposed experimental setup of Run Group F. Spectator protons are reconstructed with a fastMC that was developed based on the performance of the eg6 CLAS run.

## 3.2 Particle Identification

The final state of proton-tagged neutron DVCS event consists of three particles: an electron, a proton, and a real photon, while in addition to these particles a neutron is required to be detected in the fully exclusive neutron's DVCS events. To identify the DVCS events, we first identify the different particles of interest. Then, events with three and four detected final-state particles will be further filtered by imposing the energy-momentum conservation laws.

For the identification of the different final state particle we use the official CLAS12 EventBuilder [47]. In the following subsections we detail the main requirement for the different particles.

## 3.2.1 Electron Identification

For charged particles in the forward detectors, the Event Builder first assigns e -(PID= 11) or e + (PID= -11) (depending on the bending direction of the reconstructed track in the DC) if a particle satisfies all corresponding HTCC and ECAL requirements, and has an associated FTOF hit:

- · 2.0 photoelectrons in HTCC.
- · 60 MeV in PCAL.
- · 5-sigma cuts on a parameterized momentum-dependent sampling fraction where 'sampling fraction' is ECAL visible energy deposition (PCAL+Inner+Outer) divided by DC momentum.

In addition to these initial selections cuts, some regions of CLAS12 have to be excluded from the analysis to ensure an accurate detection of the different particles. For instance, an electron that hits the edge of the EC will have only part of its electromagnetic shower contained within the detector. Also, the structure of the torus magnet divides CLAS12 into six separate sectors, which makes edge effects non-negligible. Figure 3.1 shows the kinematic distributions of the DVCS electrons being detected and reconstructed in the forward detector of CLAS12 spectrometer in terms of the energy as a function of the polar angle ( θ ) and the azimuthal angle ( φ ) as a function of θ .

## 3.2.2 Proton Identification

The slow recoiling final state protons will be detected within the BONuS12 RTPC. As the RTPC is physically within the CLAS12 simulation, but the track reconstruction is not finalized yet within CLAS12 reconstruction, we refer to a realistic fastmc that reproduces the expected performance of this detector. This fastmc has be developed from the BONuS6 experiment and tuned very precisely after EG6 experiment [48], which have used very similar

Figure 3.1: Electron's energy as a function of it's polar angle (left) and the azimuthal angle as a function of the polar angle (right), for n-DVCS events. Forward-CLAS12 acceptance and physics cuts are included.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-85613d4a-10b0-4046-877a-70fdd1774e89.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-c9049334-f6bc-4f47-9208-10cadd5952d8.jpg)

RTPCs. The geometry parameters in this fastmc have been updated with BONuS12's geometry. This fastmc smears the proton's kinematics and applies acceptance functions. Regarding the smearing, the momentum, polar angle, azimuthal angle, and z-vertex of the protons are smeared with Gaussians using the observed tracking resolutions of the RTPC (see chapter 2, section 3 in [48]). For the acceptance, the fastmc:

- · ensures that the proton's track intersects with the cathode.
- · rejects the track if it goes to the upstream or the downstream detector windows.
- · applies the RTPC's thresholds on the momentum and the polar angle.

We do not specifically apply energy loss and multiple scattering in our fastmc, but we do apply resolution effects based on experimental data that include both effects. Figure 3.2 presents the kinematic distributions of the recoiling final state protons within the active volume of the BONUS12 RTPC.

## 3.2.3 Neutrals Identification

For neutrals, only photon (PID= 22) and neutron (PID= 2112) are considered, and particle identification is assigned based on simple beta cut at 0.9. Currently only one timing response is used for this, and for ECAL the prioritization is: PCAL, EC Inner, EC Outer. The momentum direction is assigned based on the neutral's ECAL cluster position and the vertex of the charged particle used to determine the start time. For photons, the energy (magnitude of the momentum) is calculated from ECAL visible energy and momentumdependent sampling fraction, while for neutrons energy is calculated from beta. For the central detector, CND is treated similarly as ECAL, except only neutron PID is assigned based on beta and the cut is at 0.8. Figure 3.3 presents the kinematic distributions of the neutron-DVCS photons being detected in the CLAS12 forward detector. Figure 3.4 shows

Figure 3.2: Recoiling proton's momentum as a function of it's polar angle (left) and the azimuthal angle as a function of the polar angle (right), from n-DVCS events. BONuS12 RTPC acceptance and physics cuts are included.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-74fec16b-cc56-4f29-b222-91c8e6da794b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-dc0f58cd-4b72-44af-ad2c-76a06c73e397.jpg)

similar distributions of the detected final state neutrons being detected in both the forward and the central detectors of CLAS12 spectrometer. As can be seen from figures 3.1, 3.2, 3.3, 3.4, the DVCS electrons and photons are produced very forward, which is in a total agreement with our DVCS measurements during the 6 GeV era of CLAS, while slow recoiling protons are emitted more evenly in the polar angle range within the acceptance of the BONuS12 RTPC. The final neutrons are mostly emitted above θ = 40 , which is outside the acceptance of the ◦ forward detector of CLAS12 and was the main reason to upgrade CLAS12 with the Central neutron Detector during Run Group B to measure the neutron DVCS semi-exclusively (E1211-003) [28].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-5baa9002-08a7-42b5-a40d-4f844ec58c6b.jpg)

Figure 3.3: Photon's energy as a function of it's polar angle (left) and the azimuthal angle as a function of the polar angle (right), for n-DVCS events. Forward-CLAS12 acceptance and physics cuts are included.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-51e6f273-7657-4eab-b740-9fce0a6a2b85.jpg)

Figure 3.4: Recoiling neutron's momentum as a function of it's polar angle (left) and the azimuthal angle as a function of the polar angle (right), from n-DVCS events. BONuS12 forward and central acceptance and physics cuts are included.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-f6b138c0-4e48-4990-94fd-ad9e4ff48032.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-118e8c61-1e11-40ee-a5c8-9900069a1f89.jpg)

## 3.3 Beam-Spin Asymmetry

The beam spin asymmetry of a longitudinally polarized electron beam on an unpolarized target ( A LU ) is defined as:

$$A _ { L U } = \frac { 1 } { P _ { B } } \frac { N ^ { + } - N ^ { - } } { N ^ { + } + N ^ { - } }.$$

where P B is the beam polarization, and N + and N -are the number of events detected with positive and negative electron helicity, respectively, normalized to the life-time gated Faraday cup charge for each helicity. The statistical uncertainty of A LU is

$$\sigma _ { A _ { L U } } = \frac { 1 } { P _ { B } } \sqrt { \frac { 1 - ( P _ { B } A _ { L U } ) ^ { 2 } } { N } }$$

where N (= N + + N -) is the total number of measured DVCS events in each bin.

In order to evaluate the systematic uncertainties stemming from the DVCS selection cuts, the eg6-analysis was repeated with changing the width of the exclusive cuts. The resulting systematic uncertainty to the A LU asymmetry was around 6% for the incoherent DVCS channel. Because of the important improvement we expect with BONuS12 RTPC in terms of resolutions, we expect this uncertainty to be reduced to 4%.

It is particularly convenient to use the BSA A LU as a DVCS observable, because most of the experimental systematic uncertainties, such as normalization and efficiencies that appear in the cross sections cancel out in the asymmetry ratio. However, some systematic uncertainties remain and they still contribute to the measured A LU . The main known sources of systematic uncertainties are: the DVCS selection cuts, the fitting sensitivity to our binning, the beam polarization and the background (exclusive π 0 ) acceptance ratio. In the following, we present estimates of the contribution from each source based on our prior knowledge during CLAS-eg6 DVCS analysis [48].

Regarding the sensitivity of the fit results to our binning, the eg6 data were binned into two different bin sizes in φ and the reconstructed asymmetries were compared. The associated systematic uncertainty for A LU at φ = 90 ◦ was found to be 7.1%. For the proposed measurements, we expect to achieve higher statistics and smaller bin sizes, and therefore we reduced the expected systematics to 3%.

The total systematic uncertainty is estimated to be around 11%. To be conservative, in particular because expected asymmetries on neutrons are much smaller than on protons, we used an increased total of 20% systematic uncertainty for our projections. This is added quadratically to the statistical error bars in each bin of the reconstructed asymmetry.

The beam polarization will be measured during the experiment by the Hall B Møller polarimeter. This polarimeter measures the helicity dependence of the Møller electron yield to obtain the beam polarization. The precision of the Hall B Møller polarimeter was measured to be around 3.5% [49]. We assume a 3.5% systematic uncertainty on the measured asymmetries similar to what was achieved during the 6 GeV run.

## 3.4 Projections

Once the final state particles have been reconstructed and identified after passing all the physics and the geometry cuts. To consider an event as a DVCS, it has to pass sets of requirements: DVCS characteristic cuts and exclusivity cuts.

## DVCS characteristic cuts

- · Q 2 &gt; 1 GeV : 2 to ensure that the interaction occurs at the partonic level and the applicability of the factorization in the DVCS handbag diagram.
- · a cut on the invariant mass ( W ∗ ) of the virtual photon and the target neutron system to be greater than 2 GeV/c 2 . This cut avoids the region of excitation of the neutron to resonances.

We reconstruct two types of neutron DVCS events as listed at the beginning of this chapter. Based on the approved Run Group F electron-nucleon running luminosity, a total of 9 million tagged and 850K fully exclusive neutron DVCS events will be collected during the 35 PAC days. In the following two subsections, we present the DVCS exclusivity cuts and our proposed binning for each set of data selection, and finally the projections of the proposed measured beam-spin asymmetries.

## 3.4.1 Proton-Tagged neutron DVCS Projections

After identifying the proton-tagged neutron DVCS events, i.e., having only one electron, one photon, and one proton in the final state, we further filter them by imposing the energymomentum conservation laws. Figure 3.5 shows the missing mass squared distribution of the

identified proton-tagged neutron DVCS events in addition to the missing energy distribution of these events. We apply an additional cut on the reconstructed missing mass squared to further clean the selected events as shown by the vertical red-dashed lines.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-e9a80345-c3ed-48cb-a4ea-eaeec30c6897.jpg)

s

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-49a72227-3242-45ec-b961-c1bcecfec716.jpg)

s

Figure 3.5: The distributions of the missing mass squared (left) and the missing energy (right) of the identified proton-tagged neutron DVCS events. The DVCS exclusivity cuts are represented by the vertical red-dashed lines. The black distributions represent the incoherent DVCS event candidates before the exclusivity cut. The shaded distributions represent the DVCS events that passed the cut on the missing mass squared.

The spectator approximation assumes that the recoil proton is on its mass shell when the electron strikes the neutron and it gains neither energy nor momentum during the interaction. In order to study the effect of Fermi motion on our measured DVCS observable, we specifically use the kinematics of the photons to determine the transferred momentum squared t . Figure 3.6 shows the distribution of Q 2 as a function of x ∗ and x ∗ as a function of t of the identified proton-tagged neutron DVCS events. On the right side plot of figure 3.6 we show the proposed binning in x ∗ versus -t space. The data will be binned three-dimensionally into 108 bins. That is, 12 bins in x ∗ vs. -t , and then nine bins in the azimuthal angle ( φ ) for each of the 12 bins. The data are integrated over the full measured range of Q 2 .

Figure 3.7 presents the reconstructed proton-tagged neutron DVCS A LU as a function of φ in bins of x ∗ vs -t .

## 3.4.2 Fully exclusive n-DVCS projections

Our aim is to investigate the Fermi motion and the final state interaction (FSI) effects on the measured neutron DVCS beam-spin asymmetry. This can be achieved by comparing

Figure 3.6: The distributions of the proton-tagged neutron DVCS events in terms of Q 2 versus x ∗ (left) and x ∗ versus -t (right). On the right we show the binning we propose in x ∗ versus -t space.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-e224616c-df6e-4661-8d57-fe5668d79f2b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-1b29c9eb-6832-442f-9e3e-ef256918c732.jpg)

the measured proton-tagged neutron DVCS beam-spin asymmetry, that is d e, e p γ ( ′ s ) X , to the measured full exclusive neutron DVCS channel, that is d e, e np γ ( ′ s ). The selection of the proton-tagged neutron DVCS events have been presented in the previous subsection, although a different way of binning the data will be used here, as will be shown in this section.

As mentioned previously, a total of 850K fully exclusive events will be collected within the approved running luminosity and experimental setup of Run Group F. Similarly to the protontagged neutron DVCS events, we use the kinematics of the detected spectator proton to define the modified Lorentz invariant, x ∗ , and we define the transferred momentum squared, t , using the photons to investigate the initial Fermi motion effect on our DVCS observable of interest, A LU . Figure 3.9 shows the distributions of Q 2 as a function of x ∗ and x ∗ as a function of -t for the identified fully exclusive events. Before binning our data in the space of the momentum and the polar angle of the detected spectator momentum, we apply an initial cut on x ∗ vs. -t , as can be seen in Figure 3.9 which stands for the fully exclusive detected n-DVCS events. Similar cut is applied on the x ∗ vs. -t space of the proton-tagged neutron DVCS events. After that, both data sets are binned into 6 bins in the momentum of the spectator proton and its polar angle as shown in Figure 3.10. Finally, the data of each bin in p s versus θ s is binned into 9 φ bins for the fully exclusive events and 12 φ bins for the proton-tagged neutron DVCS events. The reconstructed A LU is presented in Figure 3.11 as a function of the hadronic angle φ for the proton-tagged neutron DVCS events in black

Regarding the fully exclusive neutron DVCS events selection, events with the four final state particles will be identified after applying all the geometry and physics cuts on the individual final state particle. After identifying these events, the exclusivity of the DVCS events is ensured by imposing a set of constraints based on the four-momentum conservation in the reaction ed → e p nγ ′ s . The distributions for the exclusive variables are shown in figure 3.8.

points and for the fully exclusive n-DVCS events in blue points. The error bars in these projections include both statistical and systematic uncertainties. We are considering 20% systematic uncertainties in our projections here as well. As one can see from our projections, the systematic uncertainties have the major contribution in the precision of our asymmetries and even with the extremely conservative assumption of 20% consideration, we will observe very precise beam-spin asymmetries.

Figure 3.12 presents the predicted precision on the reconstructed A LU ratio between the proton-tagged neutron DVCS events and the fully exclusive neutron DVCS events as a function of the spectator proton polar angle ( θ s ) in the two momentum bins shown in figure 3.10. For each data selection (tagged or exclusive) in each bin in p s versus θ s phasespace, each reconstructed A LU signal versus φ has been fitted with the form a 0 sin( φ ) 1 + a 1 cos( φ ) . Then, A LU ratio is defined as A LU at φ = 90 ◦ of the reconstructed fitted tagged DVCS events to the equivalent A LU at φ = 90 ◦ of the reconstructed fitted fully exclusive DVCS events. In other words, A LU ratio is the a 0 fit parameter of the tagged DVCS to a 0 fit parameter of the fully exclusive DVCS beam-spin asymmetry. Such a measurement would give us insights about the size effect of the FSI and the Fermi motion on deep exclusive reactions. Since only models predicting the FSI effect on the PDFs are available, such valuable measurement will trigger theorists to come up with new models explaining our results.

Figure 3.7: Projected beam-spin asymmetries as a function of the hadronic angle φ h in the 12 bins in x ∗ vs -t space. The error bars include both the statistical and the systematic uncertainties added quadratically.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-b2f77671-5910-4753-94d2-fc0c27ecee90.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-0bc70238-fe49-4567-8648-f7a389b4c67c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-8b376ef0-33b9-4069-bde5-299fea489174.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-9e9ff302-a81b-4464-b6d0-8c317aba2338.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-5d00e0e7-8d25-4515-8cca-8e867f948fad.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-37b8593a-07ee-4456-97cc-f2cce68b18d3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-2287e8c9-d246-4654-9b32-c4ef99b6fecd.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-ad16d1af-514b-42f9-9a24-1eaa26e6d508.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-199bf00d-2947-4529-859e-73fd7d712e95.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-dfaeeb10-e64f-4190-847e-55fc6a7adab2.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-2c63d10c-109f-4ae0-bbc8-b7ed4fb69324.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-237eef24-50ef-4038-87d7-0535d6c0c396.jpg)

0.1

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-63efed25-b921-4372-8aaf-1f920d6d8cfa.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-3264e597-db49-41dc-861c-c998c948adb8.jpg)

s

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-f3ebc457-2255-466e-8982-b65a9c687108.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-836d8060-6386-417a-b519-83f67c875e64.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-3c39b7fa-3330-4242-b65e-1c32e60b07b3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-449d7d93-7d56-4cc2-bec2-f80bd3210b89.jpg)

0.5

X

Figure 3.8: The distributions from left to right and from top to bottom are: missing mass squared, missing energy, missing total momentum, the x-component of the missing momentum, the y-component, and the transverse missing momentum in the e p nγ ′ s finalstate system. The DVCS exclusivity cuts are represented by the vertical red-dashed lines. The black distributions represent the DVCS event candidates before the exclusivity cuts. The shaded distributions represent the DVCS events that passed all of these cuts.

Figure 3.9: The distributions of the fully exclusive neutron DVCS events in terms of Q 2 versus x ∗ (left) and x ∗ versus -t (right). On the right we show the bin of interest in x ∗ versus -t space.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-1a0929cd-d77a-4230-9e7b-2ec4f8a5d8cf.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-84d8de6a-eb48-4fb5-bcad-a6417165459e.jpg)

Protons: P vs.

:: epngamma

q

Figure 3.10: Data binning in p s versus θ s .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-9d7490e6-7044-4a7d-a3a3-eb6d1508a07b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-49f5db4a-2d43-42a9-b928-097e3ce2a2f8.jpg)

h

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-8b428bc9-c28c-4b9d-b83d-6674dd9f846d.jpg)

h

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-7e77314c-955b-45eb-a3bf-0c137235e97e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-2f768c03-1ae5-41d6-ad85-8940c4482891.jpg)

h

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-0361c22e-5539-4beb-aa54-4fd70395556b.jpg)

h

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-bbf9e66d-0c47-4c78-a3ba-d66dd9ab0583.jpg)

h

Figure 3.11: Projected beam-spin asymmetries as a function of the hadronic angle φ h in the binning of p s vs θ s space for the proton-tagged neutron DVCS events in black points and for the fully exclusive n-DVCS events in blue points. The error bars include both the statistical and the systematic uncertainties added quadratically.

Figure 3.12: The projected precision on the desired measurable A LU ratio at φ = 90 , from ◦ fitting A LU signals in figure 3.11, of the proton-tagged neutron DVCS events to the fully exclusive neutron DVCS events as a function of the spectator proton θ s in the two bins of the spectator proton momentum ( p s ).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00949v1-figure-57b3cd8e-9341-4164-acd1-3f84db168cea.jpg)

## Chapter 4

## Additional Physics Measurements with Run Group F

The combination of the high luminosity available at Jefferson Lab, the large acceptance of CLAS12 detector and the BONuS12 RTPC offers an amazing opportunity to advance our understanding of long standing questions in QCD. This new development in detection capabilities will allow the study of medium modification with a handle on Fermi motion uncertainties and FSI effects. It is therefore clear that the focus of the neutron DVCS proposal is only a fraction of the physics that can be achieved by successfully analyzing the Run Group F data with highly polarized beam. This data will be a gold mine, which will allow us to investigate in a unique way several important physics questions and conquer new territories in the nuclear QCD land scope. Some of the topics of interest to increase the physics outcome of BONuS12 polarized beam data are:

- · Coherent DVCS and deeply virtual mesons production (DVMP) off deuteron. For DVMP, we can study for example π 0 , φ , ω and ρ mesons.
- · Incoherent proton DVCS and DVMP off deuterium.
- · Deep virtual π 0 production off neutron, which is interesting by itself but also the background of DVCS measurements.
- · Semi-inclusive reaction p(e,e'p)X to study the π 0 cloud of the proton and D e, e pp ( ′ S ) X to study the π -cloud of the neutron, at very low proton momenta.
- · Transverse momentum distributions (TMDs) on the neutron (twist-3).
- · The medium modification of the transverse momentum dependent parton distributions.
- · Final state interactions through the 5 th structure function in D e, e p ( ′ s ) n .

## Summary

In summary, polarizing the electron beam during the approved E12-06-113 experiment (BONuS12) will allow us to investigate in a unique way many aspects of QCD within the GPD framework. The approved E12-11-003 experiment, 'Deeply Virtual Compton Scattering on the Neutron with CLAS12 at 11 GeV' is set to measure the n-DVCS beam-spin asymmetry by directly detecting the struck neutron in the reaction γ ∗ + d → n + γ +( ). p In contrast, we intend to measure the neutron DVCS beam-spin asymmetry by tagging the spectator slow-recoiling proton in addition to measuring the fully exclusive neutron DVCS channel. The first channel will enrich our knowledge about the partonic structure of the quasi-free neutrons, while the fully exclusive neutron DVCS measurement will be a golden data set to understand the Fermi motion and final state interaction effects on the measured DVCS beamspin asymmetries. While this proposal is focusing only on the neutron DVCS measurements, highly polarizing the beam during Run Group F will be giving us a golden chance to measure additional physics topics, increasing the physics outcome of the approved beam time and advancing our understanding on many aspects of QCD.

## Bibliography

- [1] M. Burkardt, 'Impact parameter dependent parton distributions and off-forward parton distributions,' Phys. Rev. , vol. D62, p. 071503, 2000.
- [2] M. Diehl, 'Generalized parton distributions in impact parameter space,' Eur.Phys.J. , vol. C25, pp. 223-232, 2002.
- [3] A. V. Belitsky and D. Mueller, 'Nucleon hologram with exclusive leptoproduction,' Nucl. Phys. , vol. A711, pp. 118-126' 2002.
- [4] M. Burkardt, 'Transverse deformation of parton distributions and transversity decomposition of angular momentum,' Phys. Rev. , vol. D72, p. 094020, 2005.
- [5] X.-D. Ji, 'Gauge-Invariant Decomposition of Nucleon Spin,' Phys. Rev. Lett. , vol. 78, pp. 610-613, 1997.
- [6] E. Leader and C. Lorc´, 'The angular momentum controversy: What's it all about and e does it matter?,' Phys. Rept. , vol. 541, no. 3, pp. 163-248, 2014.
- [7] M. Hattawy et al. , 'Deeply Virtual Compton Scattering measurement off bound protons,' 2018.
- [8] M. Amaryan et al. , 'The Structure of the Free Neutron at Large x-Bjorken (E12-06113),' A proposal to PAC 30 , 2006.
- [9] D. M¨ uller, D. Robaschik, B. Geyer, F. M. Dittes, and J. Horejsi, 'Wave functions, evolution equations and evolution kernels from light ray operators of QCD,' Fortsch. Phys. , vol. 42, pp. 101-141, 1994.
- [10] X. Ji, 'Deeply virtual compton scattering,' Phys. Rev. D , vol. 55, pp. 7114-7125, Jun 1997.
- [11] A. V. Radyushkin, 'Scaling limit of deeply virtual Compton scattering,' Phys. Lett. , vol. B380, pp. 417-425, 1996.
- [12] A. V. Radyushkin, 'Nonforward parton distributions,' Phys. Rev. D , vol. 56, pp. 55245557, Nov 1997.

- [13] K. Goeke, M. V. Polyakov, and M. Vanderhaeghen, 'Hard exclusive reactions and the structure of hadrons,' Prog. Part. Nucl. Phys. , vol. 47, pp. 401-515, 2001.
- [14] M. Diehl, 'Generalized parton distributions with helicity flip,' Eur. Phys. J. , vol. C19, pp. 485-492, 2001.
- [15] A. V. Belitsky, D. Mueller, and A. Kirchner, 'Theory of deeply virtual Compton scattering on the nucleon,' Nucl. Phys. , vol. B629, pp. 323-392, 2002.
- [16] A. V. Belitsky, D. Mueller, L. Niedermeier, and A. Schafer, 'Leading twist asymmetries in deeply virtual Compton scattering,' Nucl. Phys. , vol. B593, pp. 289-310, 2001.
- [17] S. Stepanyan et al. , 'Observation of exclusive deeply virtual compton scattering in polarized electron beam asymmetry measurements,' Phys. Rev. Lett. , vol. 87, p. 182002, Oct 2001.
- [18] A. Airapetian et al. , 'Measurement of the beam-spin azimuthal asymmetry associated with deeply-virtual compton scattering,' Phys. Rev. Lett. , vol. 87, p. 182001, Oct 2001.
- [19] A. Airapetian et al. , 'Beam-charge azimuthal asymmetry and deeply virtual compton scattering,' Phys. Rev. D , vol. 75, p. 011103, Jan 2007.
- [20] F. X. Girod et al. , 'Measurement of Deeply virtual Compton scattering beam-spin asymmetries,' Phys. Rev. Lett. , vol. 100, p. 162002, 2008.
- [21] M. Defurne et al. , 'E00-110 experiment at jefferson lab hall a: Deeply virtual compton scattering off the proton at 6 gev,' Phys. Rev. C , vol. 92, p. 055202, Nov 2015.
- [22] M. Mazouz et al. , 'Deeply virtual compton scattering off the neutron,' Phys. Rev. Lett. , vol. 99, p. 242501, Dec 2007.
- [23] G. Gavalian et al. , 'Beam spin asymmetries in deeply virtual compton scattering (dvcs) with clas at 4.8 gev,' Phys. Rev. C , vol. 80, p. 035206, Sep 2009.
- [24] E. Seder et al. , 'Longitudinal target-spin asymmetries for deeply virtual compton scattering,' Phys. Rev. Lett. , vol. 114, p. 032001, Jan 2015.
- [25] H. S. Jo et al. , 'Cross sections for the exclusive photon electroproduction on the proton and Generalized Parton Distributions,' Phys. Rev. Lett. , vol. 115, no. 21, p. 212003, 2015.
- [26] M. Guidal, H. Moutarde, and M. Vanderhaeghen, 'Generalized Parton Distributions in the valence region from Deeply Virtual Compton Scattering,' Rept. Prog. Phys. , vol. 76, p. 066202, 2013.
- [27] R. Dupr´ e, M. Guidal, and M. Vanderhaeghen, 'Tomographic image of the proton,' Phys. Rev. D , vol. 95, p. 011501, Jan 2017.

- [28] A. Fradi et al. , 'Deeply Virtual Compton Scattering on the Neutron with CLAS12 at 11 GeV (E12-11-003),' A proposal to PAC 37 , 2011.
- [29] R. D. Ball et al. , 'Parton distributions for the LHC Run II,' JHEP , vol. 04, p. 040, 2015.
- [30] M. Mazouz et al. , 'Deeply virtual compton scattering off the neutron,' Phys. Rev. Lett. , vol. 99, p. 242501, 2007.
- [31] N. d'Hose, S. Niccolai, and A. Rostomyan, 'Experimental overview of Deeply Virtual Compton Scattering,' Eur. Phys. J. , vol. A52, no. 6, p. 151, 2016.
- [32] R. Dupr´ e, M. Guidal, S. Niccolai, and M. Vanderhaeghen, 'Analysis of Deeply Virtual Compton Scattering Data at Jefferson Lab and Proton Tomography,' Eur. Phys. J. , vol. A53, no. 8, p. 171, 2017.
- [33] H. Moutarde, P. Sznajder, and J. Wagner, 'Border and skewness functions from a leading order fit to DVCS data,' Eur. Phys. J. , vol. C78, no. 11, p. 890, 2018.
- [34] C. Ciofi degli Atti, L. P. Kaptari, and B. Z. Kopeliovich, 'Final state interaction effects in semiinclusive DIS off the deuteron,' Eur. Phys. J. , vol. A19, pp. 145-151, 2004.
- [35] C. Ciofi degli Atti and B. Z. Kopeliovich, 'Final state interaction in semiinclusive DIS off nuclei,' Eur. Phys. J. , vol. A17, pp. 133-144, 2003.
- [36] A. Airapetian et al. , 'Nuclear-mass dependence of azimuthal beam-helicity and beamcharge asymmetries in deeply virtual Compton scattering,' Phys. Rev. , vol. C81, p. 035202, 2010.
- [37] S. Liuti and S. K. Taneja, 'Microscopic description of deeply virtual Compton scattering off spin-0 nuclei,' Phys. Rev. , vol. C72, p. 032201, 2005.
- [38] V. Guzey, A. W. Thomas, and K. Tsushima, 'Medium modifications of the bound nucleon GPDs and incoherent DVCS on nuclear targets,' Phys. Lett. , vol. B673, pp. 914, 2009.
- [39] 'CLAS12 Technical Design Report,' 2008.
- [40] S. Kuhn , Private communication .
- [41] A. E. Alaoui and E. Voutier CLAS Note 2009-024 .
- [42] M. Vanderhaeghen, P. A. M. Guichon, and M. Guidal, 'Deeply virtual electroproduction of photons and mesons on the nucleon: Leading order amplitudes and power corrections,' Phys. Rev. D , vol. 60, p. 094017, Oct 1999.

- [43] M. Guidal, M. V. Polyakov, A. V. Radyushkin, and M. Vanderhaeghen, 'Nucleon formfactors from generalized parton distributions,' Phys. Rev. , vol. D72, p. 054013, 2005.
- [44] M. Lacombe, B. Loiseau, J. M. Richard, R. V. Mau, J. Cˆt´, P. Pir`s, and R. de Tourreil, o e e 'Parametrization of the paris n -n potential,' Phys. Rev. C , vol. 21, pp. 861-873, Mar 1980.
- [45] CLAS12 GEMC simulation, https: // gemc. jlab. org/ gemc/ html/ index. html .
- [46] CLAS12 CLARA, https: // claraweb. jlab. org/ clara/ .
- [47] CLAS12 Event-Builder, https: // clasweb. jlab. org/ wiki/ index. php/ CLAS12\_ EventBuilder# Particle\_ Identification .
- [48] M. Hattawy and others (EG6 Working Group), 'Deeply Virtual Compton Scattering off 4 He,' CLAS internal analysis note , 2016.
- [49] J. M. Grames et al. , 'Unique electron polarimeter analyzing power comparison and precision spin-based energy measurement,' Phys. Rev. ST Accel. Beams , vol. 7, p. 042802, Apr 2004.
M. Hattawy M. Amaryan S. Bültmann G. Dodge N. Dzbenski C. Hyde S. Kuhn D. Payette J. Poudel L. Weinstein R. Dupré M. Guidal D. Marchand C. Muñoz S. Niccolai E. Voutier K. Hafidi Z. Yi T. Chetry L. El-Fassi N. Baltzell G. Gavalian F. X. Girod S. Stepanyan I. Albayrak E. Christy A. Nadeeshani N. Kalantarians C. Ayerbe Gayoso D. Jenkins
Neutron DVCS Measurements with BONuS12 in CLAS12
The three-dimensional picture of quarks and gluons in the nucleon is set to
be revealed through deeply virtual Compton scattering (DVCS). With the absence
of a free neutron target, the deuterium target represents the simplest nucleus
to be used to probe the internal 3D partonic structure of the neutron. We
propose here to measure the beam spin asymmetry (BSA) in incoherent neutron
DVCS together with the approved E12-06-113 experiment (BONuS12) within the run
group F, using the same beam time, simply with addition of beam polarization.
The DVCS BSA on the quasi-free neutron will be measured in a wide range of
kinematics by tagging the scattered electron and the real photon final state
with the spectator proton. We will also measure BSA with all final state
particles detected including the struck neutron. The proposed measurements is
complementary to the approved CLAS12 experiment E12-11-003, which will also
measure the quasi-free neutron DVCS by detecting the scattered neutron, but not
the spectator proton. Indeed, besides providing more data for neutron DVCS,
this experiment will allow a comparison of the measurement of the BSA of
neutron DVCS from the approved E12-11-003 with the measurements using the two
methods proposed herein. This comparison will help to understand the impact of
nuclear effects, such as the final state interactions (FSI) and Fermi motion on
the measurement of the neutron DVCS.