![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00040v1-figure-1dc2237f-6497-42ff-a813-5652cc92cc5c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00040v1-figure-1d01c216-3854-4180-aebc-2b768152157b.jpg)

## Probing the flavor dependence of proton's light-quark sea in the SeaQuest experiment

## Jason Dove ∗†

Department of Physics, University of Illinois at Urbana-Champaign,

Urbana, Illinois 61801

E-mail:

jdove2@illinois.edu

Surprisingly large flavor asymmetry of the light-quark sea in the proton was reported in deepinelastic scattering and Drell-Yan experiments. The Bjorkenx dependence of the ¯ d / ¯ ratio exu tracted from the Fermilab E866 experiment also revealed an intriguing drop at the highest values of x . The Fermilab E906/SeaQuest experiment was designed to measure ¯ d / ¯ with improved acu curacy at high Bjorkenx . By detecting high-mass dimuon events produced in the interaction of 120 GeV proton beam with liquid hydrogen (LH2) and deuterium (LD2) targets, the SeaQuest experiment probes the ¯ d / ¯ ratio up to u x ≈ 0 45. . Data collection has been completed and we report the status of data analysis for SeaQuest. Two different methods for extracting the LD2/LH2 Drell-Yan cross section ratios for 0 1 . ≤ x ≤ 0 45 are discussed. . These ratios are compared with calculations using various PDFs.

XXVII International Workshop on Deep-Inelastic Scattering and Related Subjects - DIS2019 8-12 April, 2019

Torino, Italy

∗ Speaker.

† Representing the Fermilab SeaQuest Collaboration.

https://pos.sissa.it/

## 1. Flavor asymmetry of the nucleon sea

First indications for a flavor asymmetric nucleon sea came from the measurements of the Gottfried sum [1]

$$S _ { G } = \int _ { 0 } ^ { 1 } ( F _ { 2 } ^ { p } ( x ) - F _ { 2 } ^ { n } ( x ) ) \frac { d x } { x } = \frac { 1 } { 3 } + \frac { 2 } { 3 } \int _ { 0 } ^ { 1 } [ \bar { u } ^ { p } ( x ) - \bar { d } ^ { p } ( x ) ] d x. \quad \.$$

For a flavor symmetric sea, ¯ u x ( ) = ¯ d x ( ) , the Gottfried Sum Rule (GSR), SG = 1 3 is obtained. First measurement of SG at SLAC [2] gave SG = 0 28, suggesting that the GSR is violated. This early re-. sult prompted Field and Feynman [3] to suggest that the Pauli-blocking effect from the uud valence quarks in the proton would inhibit the ¯ uu sea more than the ¯ dd sea, leading to a flavor asymmetric ¯ d &gt; ¯ sea. The most precise measurement of u SG was performed by the NMC Collaboration [4, 5], giving SG = 0 235 . ± 0 026, establishing the violation of the GSR at a 4 . s level.

While the violation of the GSR could be explained by a flavor asymmetric sea, it still left open the possibility of pathological behavior of F 2 in the unmeasured x region. A direct and independent experimental tool to probe the flavor asymmetry is to use the proton-induced Drell-Yan process [6]. By measuring the s ( p + ) d / s ( p + ) p Drell-Yan cross section ratio on hydrogen and deuterium targets, the ¯ d x ( ) / ¯ u x ( ) ratio can be extracted as a function of x . The NA51 Collaboration at CERN [7] found ¯ u d / ¯ = 0 51 . ± 0 04 . ± 0 05 at . 〈 x 〉 = 0 18, confirming a surprisingly large asymmetry between . the ¯ and u ¯ sea in the proton. d

The E866/NuSea experiment at Fermilab aimed at a high statistics measurement of ¯ d x ( ) / ¯ u x ( ) covering a broad range in x than the NA51 experiment. At the leading-order, the Drell-Yan cross section is given as

$$\frac { d ^ { 2 } \sigma _ { D v } } { d x _ { 1 } d x _ { 2 } } = \frac { 4 \pi \alpha ^ { 2 } } { 9 s x _ { 1 } x _ { 2 } } \sum _ { i } e _ { i } ^ { 2 } [ q _ { i } ( x _ { 1 } ) \bar { q } _ { i } ( x _ { 2 } ) + \bar { q } _ { i } ( x _ { 1 } ) q _ { i } ( x _ { 2 } ) ],$$

where qi ( x ) or ¯ qi ( x ) refers to quark or antiquark parton distributions of flavor i evaluated at momentum fraction x of the hadron ( x 1 and x 2 correspond to beam and target hadrons). At the forward rapidity region ( x 1 &gt; x 2 ), the Drell-Yan cross section ratio is approximately given by

$$\frac { \sigma _ { _ { D Y } } ( p + d ) } { 2 \sigma _ { _ { D Y } } ( p + p ) } \approx \frac { 1 } { 2 } \left ( 1 + \frac { \bar { d } ( x _ { 2 } ) } { \bar { u } ( x _ { 2 } ) } \right )$$

Equation 1.3 shows that the x -dependence of ¯ d / ¯ u can be obtained from a measurement of the s DY ( p + ) d / 2 s DY ( p + ) p ratio as a function of x 2 .

Figure 1 shows the ¯ d / ¯ ratios extracted from the NA51 [7] and E866 [8] experiments over a u broad range of x . At the low x region, ¯ d / ¯ u approaches 1 as x → 0, consistent with a symmetric light-quark sea. As x increases, ¯ d / ¯ rises linearly, reaching a maximum at u x ≈ 0 18, then starting . to drop as x further increases. The striking flavor asymmetry between ¯ and d ¯ has inspired many u theoretical models, reviewed in several articles [9, 10, 11, 12, 13]. While these models are able to explain the ¯ d &gt; ¯ asymmetry at low u x , none of the models could reproduce ¯ d &lt; ¯ at the largest u x . The E906 experiment was proposed to obtain better precision in the region x &gt; . 0 2.

## 2. E906/SeaQuest Experiment

The E906 experiment measures high-mass dimuons produced with a 120 GeV proton beam incident on liquid hydrogen (LH2) and liquid deuterium (LD2) targets. From Eq. (1.2), there is an enhancement of the Drell-Yan cross section at lower beam energy compared with E866 due to the 1 / s factor. The background due to the open-charm leptonic decay is also expected to be reduced at lower beam energy. In addition, E906 has a beam intensity of 10 12 protons per second, about 10 times that of E866.

\_

\_

The SeaQuest spectrometer includes two large dipole magnets and four hodoscope/tracking stations. The first dipole is a solid iron magnet which also serves as a beam dump. The second magnet has an open aperture allowing the momentum of a charged track to be measured. Hit information from the hodoscopes was utilized for selecting dimuons originating from the target in an FPGA-based trigger system. Two 50-cm long

Figure 1: NA51 and E866 measurement of ¯ d / ¯ u [8]. Systematic uncertainty is included for the E866 data. Parametrizations of various PDFs are compared with the data.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00040v1-figure-66699aa1-f751-4e91-8c98-e0804e8dfd91.jpg)

LH2 and LD2 targets, together with three solid nuclear targets and an evacuated flask (emptytarget), were placed on a movable table so that the targets can be changed between two beam spills. More details on the E906 spectrometer can be found in [14].

## 3. E906 Data Analysis

Data collection for E906 ran from April 2014 to July 2017. The analysis presented here is performed on the data collected until August 2015, about half of the entire data set. After applying various analysis cuts to select candidate dimuon events originating from the LD 2 target, the dimuon invariant mass spectrum is shown in Fig. 2.

The J / y peak and the y ′ shoulder are clearly visible. By selecting events with dimuon mass greater than 4.2 GeV the J / y and y ′ events can be effectively removed leaving only Drell-Yan and some background events.

In order to help remove the background from interactions that did not originate from the target, data were also collected on an empty flask. The empty flask data properly normalized can be used to subtract events which pass the analysis cuts but are in fact originating from either the dump, the flask windows, or the air upstream or downstream of the target.

A major challenge for the E906 experiment was to properly account for the accidental background originating from the coincidence of two independent single-muon tracks. Due to the high proton beam intensity in E906, the optimal removal of the accidental background is essential.

Two different methods have been used to account for the accidental background in E906. The first method examines the beam-intensity dependence of the LD2 and LH2 candidate Drell-Yan events. Because accidental events originate from two independent interactions, the rate for accidental events should have a strong intensity dependence. By plotting the ratio of deuterium yield

over hydrogen yield as a function of incident proton intensity and extrapolating to zero intensity, one can obtain the Drell-Yan cross section ratio without the contamination of accidental events.

Figure 3: a) Second-order polynomial fit to the beam-intensity dependence of the ratios of deuterium yield to hydrogen yield (LD2/LH2). The data included all Drell-Yan candidate events. b) The LD2/LH2 ratios for an x 2 bin compared with the curve from a correlated fit to all x 2 bins using Eq. (3.1)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00040v1-figure-2aebda0d-d3c2-4b13-ab30-7d22971b792f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00040v1-figure-7c625e60-b4a9-4174-b069-b8dfe205fc7e.jpg)

Fig. 3a shows the LD2/LH2 ratios for DrellYan candidate events as a function of the proton beam intensity I . From a second-order polynomial fit to the data, shown as the solid curve in Fig. 3a, the value of the s pd / 2 s pp is given as the intercept of the fit. In order to extract the s pd / 2 s pp cross section ratios as a function of x 2 , the Drell-Yan candidate events are split into seven x 2 bins and the intensity dependencies of the LD2/LH2 ratios for all x 2 bins are fitted with correlated second-order polynomials as follows:

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00040v1-figure-af935ee4-3c28-4ec3-a7cb-e23dce61e4f7.jpg)

$$R _ { n } ( I ) = P _ { 0 n } + ( P _ { 1 } + P _ { 1 } ^ { \prime } x ) I + ( P _ { 2 } + P _ { 2 } ^ { \prime } x ) I ^ { 2 } \ \ ( 3. 1 ) \quad \text{and}$$

Figure 2: Dimuon invariant mass spectrum after applying analysis cuts to select dimuon events from the target [14]. By applying an additional mass cut mass &gt; . 4 2 GeV, the J/ y and y ′ can be effectively removed.

where n ranges from 1 to 7 for the seven x 2 bins. The P ′ 1 and P ′ 2 fitting parameters accommodate the expected smooth variation of the I and I 2 dependencies among various x 2 bins. Fig. 3b shows the result of this correlated fit to one of the seven x 2 bins ( 0 16 . &lt; x 2 &lt; . 0 195 . )

The s pd / 2 s pp cross section ratios as a function of x 2 , obtained from the intercepts P n 0 in Eq. (3.1) are shown in Fig. 4. Calculations of the Drell-Yan cross section ratios in the Nextto-leading order framework, weighted by the E906 spectrometer acceptance, are shown for three sets of nucleon parton distribution functions (PDFs). Fig. 4 shows that the E906 result is in good agreement with the calculations for all PDFs at x &lt; . 0 24. The high statistics of the E866 data at this low x region put stringent constraints on the various PDFs. At large x , the curves begin to disagree with both E906 data and each other. The large discrepancy between PDFs is indicative

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00040v1-figure-62b7dbce-c312-4b6f-88f0-1c260a826b9b.jpg)

2

x

Figure 4: Comparison of cross section ratio from E906 data and various PDF fits. All PDFs and data are in good agreement at low x where the data from previous E866 was well determined. At higher x &gt; . 0 2 the uncertainty in the E866 data leads to underdetermined ¯ d / ¯ and more variety in u the PDFs. The addition of E906 data at this high x range will provide additional constraints on the PDFs.

of a larger uncertainty in the constraints from E866. The new data from E906 will clearly provide new constraints for future PDFs at the large x region

Another method we adopted to account for the accidental background is by fitting the dimuon mass spectrum including various sources of dimouns. As shown in Fig. 5, the dimuon mass spectrum is well fitted by a combination of different processes. The empty flask mass spectrum is taken from data, properly normalized. The shapes of the J / y y , ′ , and Drell-Yan mass distributions are taken from Monte-Carlo simulations.

The remaining task is to produce an appropriate shape for the accidental events. Since the accidental events primarily arise from coincidence of single muon tracks, the accidental background is produced by mixing single muon tracks. Events which contain only a single muon are mixed to ensure the kinematics of the individual tracks is consistent with a single muon, and because the events are separated in time, there is no chance that a true dimuon will be erroneously constructed. Once the relative normalization between the various components has been fixed, the fit can be projected onto other variables (for example x 2 ) and the backgrounds can be appropriately subtracted and the Drell-Yan yields can be extracted. Preliminary results from the mass-fit to LD2 and LH2 data show that the LD2/LH2 cross section ratios are consistent with those extracted from the intensity extrapolation method discussed earlier.

Figure 5: Dimuon mass spectrum after analysis cuts designed to suppress dump events and background [14]. Monte Carlos are run to produce expected shapes for J / y y , ′ and Drell-Yan. Accidental background is produced by mixing uncorrelated single tracks from data. Events with the empty flask are normalized by the number of protons incident on targets and are subtracted. The remaining data is fit to the Monte Carlo and mixed background. By applying a mass cut &gt; . 4 2 GeV, the J/ y and y ′ can be removed.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00040v1-figure-931ec21d-eabe-413b-bfd6-af9c966b1cbb.jpg)

## References

- [1] K. Gottfried. Physical Review Letters , 18:1174, 1967.
- [2] E. D. Bloom. Proc. 6th Int. Sym. on Electron and Photon Interactions at High Energies , page 227, 1974.
- [3] R. D. Field and R. P. Feynman. Physical Review D , page 2590, 1977.
- [4] P. Amaudruz et al. Physical Review Letters , 66:2712, 1991.
- [5] M. Arneodo et al. Physical Reivew D , 50:121, 1994.
- [6] S. D. Ellis and W. J. Stirling. Physical Review Letters B , 256:258, 1991.
- [7] A. Baldit et al. Physical Review Letters B , 332:244, 1994.
- [8] R. S. Towell et al. Physical Review D , 67:052002, 2001.
- [9] J. Speth and A. W. Thomas. Advanced Nuclear Physics , 24:83, 1997.
- [10] S. Kumano. Phys. Rep , 303:183, 1998.
- [11] G. T. Garvey and J. C. Peng. Prog. Part. Nucl. Phys. , 47:203, 2001.
- [12] W. C. Chang and J. C. Peng. Prog. Part. Nucl. Phys. , 79:95, 2014.
- [13] D. F. Geesaman and P. E. Reimer. Rep. Prog. Phys. , 82:042301, 2019.
- [14] Aidala et al. Nucl. Instrum. Meth. , A930:49, 2019.
Jason Dove
Probing the flavor dependence of proton's light-quark sea in the SeaQuest experiment
Surprisingly large flavor asymmetry of the light-quark sea in the proton was
reported in deep-inelastic scattering and Drell-Yan experiments. The
Bjorken-$x$ dependence of the $\bar{d}/\bar{u}$ ratio extracted from the
Fermilab E866 experiment also revealed an intriguing drop at the highest values
of $x$. The Fermilab E906/SeaQuest experiment was designed to measure
$\bar{d}/\bar{u}$ with improved accuracy at high Bjorken-$x$. By detecting
high-mass dimuon events produced in the interaction of 120 GeV proton beam with
liquid hydrogen (LH$_2$) and deuterium (LD$_2$) targets, the SeaQuest
experiment probes the $\bar{d}/\bar{u}$ ratio up to $x \approx 0.45$. Data
collection has been completed and we report the status of data analysis for
SeaQuest. Two different methods for extracting the LD$_2$/LH$_2$ Drell-Yan
cross section ratios for $0.1\leq x\leq 0.45$ are discussed. These ratios are
compared with calculations using various PDFs.