# A new analysis technique to measure the Production Charge Asymmetry at the Tevatron

###### Abstract

The charge asymmetry of bosons produced in collisions at is sensitive to the ratio of and quark distributions in the range of at . We propose an analysis technique to directly measure production charge asymmetry from events at the Tevatron and show the feasibility for this method using Monte Carlo simulations.

## I Introduction

The differential cross section for boson production in as a function of rapidity is

where () is the fraction of the proton (anti-proton) momentum carried by the struck quark, and are the quark and anti-quark parton distribution functions, and is the rapidity of the boson. The x values of the quark in the proton and antiquark in the antiproton are related to the rapidity, , of the W boson via the equation as shown in Fig 1. Here is the center of mass energy and is the mass of the W boson.

Since bosons are produced in collisions primarily by the annihilation of quarks in the proton and quarks in the anti-proton, and since and by CPT symmetry, the differential cross sections for are approximately

(2) | |||||

(3) |

Since the quark tends to carry a larger fraction of the proton’s momentum than the quark on average, the is boosted in the proton (anti-proton) direction as shown in Fig. 2. The production charge asymmetry, , in the leading-order parton model is therefore

(4) | |||||

where we use Eq. 2 and Eq. 3 and introduce the ratio . As we see in Eq. 4, there is a direct correlation between the production charge asymmetry and the ratio. A precise measurement of the production charge asymmetry serves as a valuable constraint on the and quark momentum distributions text1 ().

production at hadron colliders is identified through the process . Since the decay involves a neutrino whose longitudinal momentum is experimentally undetermined, the charge asymmetry previously has been constrained by the measured charge asymmetry of the decay leptons and as a function of the lepton pseudo-rapidity run1 (); run2 (); d0prelim (). However, as shown in Fig. 2, there is a “turn-over” in the lepton charge asymmetry due to a convolution of the production charge asymmetry and the decay. This convolution means leptons from a single pseudorapidity come from a range of rapdities and thus a range of parton values. Thus, the measured lepton asymmetry is more complicated to interpret in terms of quark distributions, and we expect the direct measurement of the asymmetry of the rapidity distribution to be a more sensitive probe of the differences between and quarks.

In this paper, we propose a new analysis technique which resolves the kinematic ambiguity of the longitudinal momentum of the neutrino to directly reconstruct the rapidity. We describe the details of our new analysis technique and outline the sources of systematic uncertainty of this measurement. Our studies are performed in the channel produced in collisions at the Tevatron. We use a realistic Monte Carlo simulation (MC@NLO) and include the effects of higher-order QCD corrections mcatnlo ().

## Ii Analysis Technique

The decay to leptons, in our case , involves a neutrino whose longitudinal momentum cannot be experimentally determined. However, we can determine the longitudinal momentum by constraining the mass in Eq. 5, which results in a two-fold ambiguity. This ambiguity can be partly resolved on a statistical basis from the known (vector-axial vector) decay distribution using the center-of-mass decay angle between the electron and the proton, , and from the and production cross-sections as a function of rapidity, . The mass constraint is

(5) |

where the mass, , is contrained to its experimentally measured value wmass-cdfrun2 (); pdg (). Events which cannot satisfy the mass constraint (and which get imaginary values of the neutrino -momentum) are due to a mis-reconstruction of the neutrino (missing) transverse energy, coord (). Therefore, in such cases, we re-scale the to the value which makes the imaginary part to be zero. This new is then used to correct the for the event.

The leading order boson production mechanism in collisions results in the boson being polarized in the direction by means of the structure of the weak interaction. The structure means that the weak current couples only to left-handed and quarks (or to right-handed and quarks). For ultra-relativistic quarks, where helicity and chirality are approximately equivalent, this results in full polarization of the produced bosons in the direction of the beam. The leptonic decay process also couples only to left-handed and right-handed (or right-handed and left-handed ). The conservation of angular momentum favors a decay with the final state lepton (neutrino or electron) at a small angle with respect to the initial state quark direction (and a similar small angle between the initial state anti-quark and final anti-lepton). The systematic shift in lepton pseudo-rapidity with respect to depending on the charge of the final state lepton is illustrated in Fig. 3 and 3, which shows the lepton pseudo-rapidity vs. rapidity for the different charges. This effect also explains the discrepancy at high rapidity between the lepton charge asymmetry and the charge asymmetry as illustrated in Fig. 2. The bias in the decay angle causes leptons at high rapidity to originate primarily from bosons produced in the opposite hemisphere.

bosons at the Tevatron are primarily produced from the valence quarks in the proton and the anti-quarks in the anti-proton and rarely from sea quarks simply because production requires at least one moderately high parton to be involved in the collision. At very large forward or backward rapidities where one very high parton must participate in the production, the production probability from the sea quarks nearly vanishes. Understanding of the sea quark contribution is important to exactly know the decay angle distributions from the structure because production by sea anti-quarks will result in the opposite polarization from valence quark production.

We use a Monte Carlo simulation with NLO QCD corrections mcatnlo () to determine the production probability with sea quarks by identifying initiating quarks as a function of . We verify the expected angular distribution of from production of with quarks in the proton and the opposite distribution with anti-quarks in the proton. For example, in Fig. 4, we show the cos distributions of in the rest frame for the case when a quark from the proton and an anti-quark from the anti-proton form the (labeled “quark”) and the case when an anti-quark from the proton and a quark from the anti-proton form the (labeled “anti-quark”). The ratio of quark (proton) and anti-quark (proton) induced production, therefore, determines the angular decay distribution. In the simulation, we measure the fraction of quark and anti-quark contributions, and parameterize the angular distributions for and the transverse momentum, . We find an empirical functional form that fits the data,

(6) |

(7) |

The parameters and are

(8) | |||||

where is the Landau distribution with most probable value and the RMS . The first term of Eq. 6 corresponds to the contribution from quarks in the proton and the second term from anti-quarks in the proton. The parameterization, , the ratio of the two angular distributions as a function of the rapidity and , is obtained from the fit to the distribution in Fig. 4.

A second relevant factor in the selection among the two rapidity solutions is the differential cross-section as a function of , . The boson production decreases sharply beyond because of the scarcity of high quarks. For instance, if one of the two possible solutions falls in the central region of rapidity and the other has , the former should receive more weight as the latter is very unlikely to be produced.

The information used to select among the two solutions can be represented by a weighting factor for each rapidity solution and charge, , can be represented as

(9) |

where the signs indicate the boson charge and indices of 1, 2 are for the two rapidity solutions.

In our analysis, we include kinematic cuts for detecting charged leptons. For event selection, we apply , , and . We also consider a multiplicative correction factor for the detector acceptance and event migration from smearing effects as shown in Figs. 5 and 5. In order to study smearing effects, we use the fact that the energy resolutions in the electromagnetic calorimeter of the Collider Dectector at Fermilab (CDF) are (central calorimeter) and (the end plug calorimeter) and in the hadronic calorimeter are (central) and (the end plug) TDR (). We randomly smear the electron and recoil hadronic energies in simulated events with a Gaussian distribution modeling their uncertainties prior to making the selection above. The correction factors are determined using a Monte Carlo program which includes both a model of the process under study as well as a simulation of the measuring apparatus. In Eq. 9, the weighting factor depends primarily on the and cross-sections, but does have some weak dependence on the assumed charge asymmetry, and thus the correction factors can be biased by computing the factors with different Monte Carlo models. Therefore, this method requires us to iterate the procedure to eliminate our measurement’s dependence on the input asymmetry. In order to confirm our analysis technique and take into account the bias from physics input variables (such as the charge asymmetry itself, the total differential cross-section and the angular distribution) we have studied the charge asymmetry measurement with different Monte Carlo models and evaluated systematic uncertainties, which are described in the next section.

## Iii Systematic Uncertainties

We consider potentially significant sources of systematic uncertainty on the charge asymmetry measurement from the assumed parton distributions, the detector resolutions and misidentifications and backgrounds. Input PDFs are used to determine the parameters of the weighting factor, and may affect the final result. The detector resolutions affect the rapidity reconstruction due to uncertainties in the calorimeter energy scale and its energy resolution, and the missing transverse energy scale also has a significant uncertainty from the boson recoil energy scale. Finally, the detector may misidentify the charge, especially from leptons at high , and there are backgrounds to at the Tevatron.

The uncertainties on the weighting factor (Eq. 9) arise from uncertainties on the momentum distribution of quarks and gluons in the proton modeled with the PDF sets used. The choice of PDF set has an effect on the shape of the distribution as well as on the ratio of quark and anti-quark in the angular decay distribution. We use the CTEQ6 error PDF sets cteq6m () and re-determine the production cross section and the angular distribution of for each error PDF set. We evaluate the uncertainty on the charge asymmetry by checking the deviation of the asymmetry values based on each calculation from the central value obtained using the best-fitted PDF set.

We also consider several experimental sources of systematic uncertainty. The scale and resolution of the electromagnetic calorimeter energy and the missing transverse energy () can change the rapidity and thus the asymmetry measurement. We use the energy uncertainties measured in wzxsec (), where the electromagnetic (EM) calorimeter energy scale and resolution was tuned in the simulation so as to fit to data. The uncertainty on the energy scale and resolution was measured to be (scale), (resolution) for central electron and (scale), (resolution) for plug electron. These values correspond to a 3 variation. The asymmetry uncertainties are estimated as the changes in the measured asymmetry when the energy scale and resolution are changed between its default and the value. The missing transverse energy () in our sample is determined by the assumption that the vector sum of all transverse energy should be zero. Since hadronic transverse energy is due to the boson recoil energy, we consider the transverse recoil energy, which is affected by multiple interactions in the event. The uncertainty on the transverse recoil energy scale is (3) wzxsec (). The charge misidentification rate and background estimates are crucial for the charge asymmetry measurement since both can directly change the measurement. We estimate these uncertainties using the charge fake rates (CFR) and background fractions (BKG) from the previous lepton charge asymmetry result from CDF run2 (). The charge fake rate is about 0.01 for and 0.04 for . The upper bound on the background fraction is for and increases to about for . We also investigate sources of any charge bias and dependence in the kinematic and geometrical acceptance of the event. An uncorrected acceptance shift of 3 central and forward electrons and 5 far forward electrons () based on measurements of data dzdy () are taken to address the effects of systematic on charge asymmetry measurement.

Table 1 summarizes the systematic uncertainties on the boson production charge asymmetry for rapidities . We compare the expected statistical uncertainty obtained by assuming an analysis using an integrated luminosity of 1 , where we also extrapolate the expected statistical uncertainty from the number of events from the previous lepton charge asymmetry result of CDF with 0.2 run2 ().

() | ||||||||
---|---|---|---|---|---|---|---|---|

rapidity | EM | Recoil | CFR | BKG | Accep | Syst. | Stat. | |

0.00.2 | 0.1 | 0.1 | 0.1 | 0.1 | 0.2 | 0.0 | 0.3 | 6.3 |

0.20.4 | 0.4 | 0.1 | 0.1 | 0.1 | 0.3 | 0.3 | 0.6 | 6.4 |

0.40.6 | 0.5 | 0.2 | 0.2 | 0.2 | 0.5 | 0.9 | 1.3 | 6.4 |

0.60.8 | 0.6 | 0.3 | 0.3 | 0.2 | 0.5 | 1.6 | 1.9 | 6.5 |

0.81.0 | 0.6 | 0.2 | 0.3 | 0.2 | 0.6 | 2.2 | 2.4 | 6.6 |

1.01.2 | 0.6 | 0.2 | 0.2 | 0.3 | 0.6 | 2.8 | 2.9 | 6.7 |

1.21.4 | 0.7 | 0.4 | 0.1 | 0.7 | 0.5 | 3.0 | 3.2 | 7.0 |

1.41.6 | 0.7 | 0.3 | 0.2 | 1.4 | 0.2 | 2.9 | 3.3 | 7.5 |

1.61.8 | 0.7 | 0.5 | 0.5 | 2.9 | 0.2 | 2.8 | 4.1 | 8.1 |

1.82.05 | 0.6 | 0.5 | 0.5 | 4.4 | 0.8 | 2.9 | 5.4 | 8.0 |

2.052.3 | 0.7 | 0.4 | 0.4 | 6.3 | 1.5 | 3.4 | 7.4 | 9.6 |

2.32.6 | 0.7 | 0.7 | 0.3 | 6.1 | 1.3 | 4.9 | 8.0 | 11.7 |

2.63.0 | 2.9 | 1.0 | 0.9 | 1.1 | 0.3 | 5.0 | 6.0 | 17.0 |

## Iv Results

We compare the expected statistical uncertainties in 1.0 fb of data at the Tevatron with the uncertainties coming from parton distribution functions (PDFs) using CTEQ6M in Fig. 6. In particular, we notice that at high rapidities ( 1.4) there is a large difference in the precision with which the as yet unmeasured production asymmetry and the previously measured asymmetry from the decay leptons scaled to 1.0 fb of integrated luminosity are known. The total systematic and statistical uncertainties on the production charge asymmetry measurement is shown in Fig. 7 with the uncertainties coming from parton distribution functions (PDFs) using CTEQ6M. Since the systematic uncertainty estimates, as summarized in Table 1, are lower than the statistical error, a direct measurement of the charge asymmetry with this method should significantly improve parameterizations of the PDFs.

In this paper, we present a study of the boson production charge asymmetry with the decaying leptonically to an electron and neutrino in collisions at the Tevatron. We propose a new analysis technique which resolves the ambiguity in the neutrino longitudinal momentum, using a realistic Monte Carlo simulation. We show that the charge asymmetry can be directly measured at the Tevatron. We conclude that by measuring the production charge asymmetry with reconstructed rapidity, the result should be one of the best determinations of the proton momentum ratio, and play an important role in global PDF fits.

## References

- (1) A. D. Martin, R. G. Roberts, and W.J. Stirling, Phys. Rev. D 50, 6734 (1994).
- (2) F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 81, 5754 (1998).
- (3) D. Acosta et al. [CDF Collaboration], Phys. Rev. D 71, 051104 (2005).
- (4) V. M. Abazov et al. [D0 Collaboration], arXiv:0709.4254 [hep-ex].
- (5) S. Frixione and B.R. Webber, J. High Energy Phys. 0206, 029 (2002), 0308, 007 (2003).
- (6) T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 99, 151801 (2007).
- (7) W.-M. Yao et al., Journal of Physics G 33, 1 (2006).
- (8) We assume a (, , ) coordinate system with the -axis in the direction of the proton beam; and are the azimuthal and polar angle respectively. The pseudorapidity is defined as , and the transverse momentum and energy as and , respectively. Missing transverse energy () is defined as where is a unit vector in the transverse plane that points from the beam-line to the calorimeter tower.
- (9) R. Blair et al.The CDF-II detector: Techinical Design Report, FERMILAB-PUB-96-390-E.
- (10) J. Pumplin, D.R. Stump, J. Huston, H.L. Lai, P. Nadolsky and W.K. Tung, J. High Energy Phys. 0207, 012 (2002).
- (11) D. Acosta et al. (CDF collaboration), Phys. Rev. Lett. 94, 091803 (2005); A. Abulencia et al. [CDF Collaboration], J. Phys. G: Nucl. Part. Phys. 34, 2457 (2007)
- (12) Jiyeon Han et al., in Proceedings of the 2007 Europhysics Conference on High Energy Physics, Manchester, 2007.