Research Article  Open Access
Comparison of HORACE and PHOTOS Algorithms for Multiphoton Emission in the Context of Boson Mass Measurement
Abstract
boson mass measurement is sensitive to QED radiative corrections due to virtual photon loops and real photon emission. The largest shift in the measured mass, which depends on the transverse momentum spectrum of the charged lepton from the boson decay, is caused by the emission of real photons from the finalstate lepton. There are a number of calculations and codes available to model the finalstate photon emission. We perform a detailed study, comparing the results from HORACE and PHOTOS implementations of the finalstate multiphoton emission in the context of a direct measurement of boson mass at Tevatron. Mass fits are performed using a simulation of the CDF II detector.
1. Introduction
The measurement of boson mass () is one of the most interesting precision electroweak observables. In the standard model (SM), the mass of boson can be calculated with higher precision [1] than the existing measurement uncertainty [2–6], thus providing motivation for improving the statistical and systematic uncertainties on the measurement. The comparison between the theoretical prediction and the measurement provides a stringent test of the SM and constrains beyondstandard model (BSM) theories.
At hadron colliders, the mass of boson is extracted from inclusively produced bosons decaying to electrons or muons and the associated neutrinos. At the Tevatron, almost pure samples of such candidate events have been identified with backgrounds typically smaller than 1%. The momenta of the decay electrons and muons have been measured with a precision of ~0.01%, allowing a boson mass measurement with precision of 0.02% [6].
The calibration of the electron and muon momenta is the single most important aspect of the measurement. In the approximation that boson undergoes a twobody decay, the distribution of the transverse momentum (, defined as the component of the momentum perpendicular to the beam axis) of the charged lepton has the characteristic Jacobian edge at half the mass of boson. In practice, electroweak radiative corrections modify the lepton spectrum, mainly due to the emission of photons from the decay lepton. If no correction was applied for this radiative process, the measurement would be biased by ≈200 MeV [7]. Throughout this paper, we use the convention .
In the first Run II measurement [7] of boson mass, WGRAD [8] and ZGRAD [9] programs were used to calculate the QED radiative correction. WGRAD and ZGRAD are exact nexttoleading order (NLO) electroweak matrix element calculations of the and processes, respectively. The effect of higherorder radiative corrections has been estimated to be about 10% [7, 10] of boson mass shift estimated from these NLO calculations. In order to increase the precision of the QED radiative correction for boson mass measurement, higherorder calculations were used, as implemented in HORACE [10–14] and PHOTOS [15, 16] programs. These programs calculate the emission of multiple photons with the appropriate rates, energy, and angular distributions.
PHOTOS uses the exact firstorder matrix element of and boson decay for the photon emission kernel. For multiphoton radiation, PHOTOS uses an iterative solution for this kernel, developed on the basis of an exact and complete phase space parametrization. This ensures not only resummation of leadinglogarithm contributions of higher orders but also the infrared region of the phase space being accurately simulated [16] (Z. Was, private communication).
HORACE is a partonlevel electroweak Monte Carlo generator for precision simulations of chargedcurrent and neutralcurrent DrellYan processes. HORACE uses the full phase space for each radiated photon and there is no ordering of the photons (i.e., in angle or transverse momentum) in multiphoton emission (C. M. Carloni Calame, G. Montagna, and A. Vicini, private communication). Two versions of the HORACE program are available. The OLD version [10, 11] implements a multiphoton emission QED parton shower algorithm for the simulation of finalstate radiation (FSR) in the leadinglogarithmic approximation, without initialstate radiation (ISR) and without interference between ISR and FSR. In this sense the OLD HORACE program is similar to the PHOTOS program, which also implements multiphoton FSR. OLD HORACE does not include full oneloop electroweak corrections, but it mimics the real radiation matrix element for the description of the photon radiation in and boson decays, in the leadinglogarithmic approximation (C. M. Carloni Calame, G. Montagna, and A. Vicini, private communication). There is also a NEW HORACE program [12, 13], which implements multiphoton ISR and FSR with interference and also matches each photon to the exact matrix element calculation of oneloop electroweak corrections and singlephoton emission (C. M. Carloni Calame, G. Montagna, and A. Vicini, private communication).
The PHOTOS program provides a generic interface to any other event generator such that all charged leptons produced by the latter can be passed through the PHOTOS FSR algorithm. We use this feature as follows. We generate and boson events for Tevatron collisions at TeV, including higherorder QCD matrix elements and QCD resummation effects, but without loops or emission of electroweak bosons. We interface these events to PHOTOS such that the events from the chain contain the QEDFSR photons added by PHOTOS. We save photons with MeV and the events are processed with a detector simulation [6, 7] to make the pseudodata and the massfitting templates. Lowering the photon threshold further has negligible effect on the results presented here, within the uncertainties quoted.
In this paper we present comparisons between the distributions and the massfitting results obtained from the OLD HORACE and PHOTOS programs.
2. Electron Channel Comparisons
To make direct comparisons between quantities sensitive to QED physics, we need to ensure that the underlying boson and lepton distributions are identical between OLD HORACE and PHOTOS. For this purpose we use the “Born” mode of OLD HORACE to generate Bornlevel and events, which are then processed through PHOTOS. The Born mode generates these purely parton processes with no radiative photons. These events are compared with events from OLD HORACE run in the QED multiphoton emission FSR mode. Both OLD HORACE and PHOTOS are run in the “exponentiation” mode, which exercises their full physics content. All of the events used in these comparisons have unit weights. For all generated events we make a generatorlevel cut on the partonic centerofmass energy GeV to remove the contribution of the photon pole for neutralcurrent events. For consistency, we also apply this cut on the chargedcurrent events.
In Figure 1 we compare the distributions for photon emission rates as well as the energy and angular distributions for the process. For these comparisons we consider photons with energy MeV; photons with lower energy than this threshold are not counted and ignored in the distributions. In addition to the number of photons emitted, we find that the distributions of the following quantities are useful to compare: , the fractional photon energy (where is the energy of the finalstate lepton), and (the angular separation between a photon and the finalstate lepton, where the pseudorapidity , being the polar angle with respect to the beam axis, and is the azimuthal angle about the beam axis).
Figure 2 shows the photon distributions separately for lowenergy and highenergy photons. The photon energy cut value of 400 MeV is chosen because roughly half of the photons are above and below this cut. These distributions show that the angular distribution is almost independent of the photon energy, allowing us to draw conclusions from the inclusive photon distributions.
In Figure 3 we compare the photon distributions for the process. These distributions are identical for the charge conjugate processes due to the collisions at the Tevatron. The comparisons between OLD HORACE and PHOTOS show good agreement in the photon emission rates and the photon energy distributions. For the process, the photon angular distribution shows about 10% difference at small angles, a difference that is not correlated with photon energy. As we show in Section 4, this difference does not cause a relative shift in the fitted mass between the two algorithms.
3. Muon Channel Comparisons
We repeat the above comparisons for the muon channel. In Figures 4 and 5 we compare the distributions for photon emission rates as well as the energy and angular distributions for the and processes, respectively. Figure 6 shows the photon distributions separately for lowenergy and highenergy photons. The rates and distributions are in good agreement between OLD HORACE and PHOTOS for the muon channel.
The photon emission at wide angle to the lepton is similar between the muon and electron channels, for both the energy and angular distributions. In the collinear region the photon emission off muons is highly suppressed due to the larger muon mass. This also reduces the number of photons emitted per event for the muon channel, relative to the electron channel.
4. Mass Fits
We quantify the impact of the small differences in photon rates and distributions between OLD HORACE and PHOTOS in terms of shifts in the fitted boson masses. For this purpose, we propagate highstatistics OLD HORACE and PHOTOS samples through the parameterized CDF detector simulation used in boson mass measurement [6]. We generate templates and pseudodata from both samples and perform the mass fits to these pseudodata using these templates, in the same manner that templates are used to fit the collider data [3–6].
We perform fits to the distributions of transverse quantities in boson events: the charged lepton , neutrino , and transverse mass . The transverse mass is defined as , where is the azimuthal angle between the charged lepton and neutrino momenta in the transverse plane. In the electron channel, we also perform a fit to the ratio of electron calorimeter energy to track momentum (), which is used by the CDF experiment to obtain the calorimeter calibration using the electron track. In boson events, we fit the distributions of boson invariantmass obtained from electron calorimeter deposition measurements (cluster mass) and from track momentum measurements of electrons and muons (track mass).
We obtain the difference between the mass fits to the HORACE pseudodata and the PHOTOS pseudodata which quantifies the relevant differences between the two QED codes. Tables 1 and 2 show the differences along with their statistical uncertainties. Table 1 uses PHOTOS templates and is essentially identical to Table 2 which uses HORACE templates. The two tables provide validation that the comparison of the two pseudodata samples does not depend on the template choice, as long as the same templates are used for the fits being compared.


The lepton track momentum is more sensitive to photon emission as compared to the calorimeter cluster energy. Radiated FSR photons directly subtract from the track momentum. The electron calorimeter cluster absorbs most of the photon radiation and is therefore less sensitive to the emission pattern. The comparison of the mass shift obtained from the calorimeter cluster and the track measurements provides an additional, sensitive test of the two programs.
5. Conclusions
We find that the QED generators OLD HORACE and PHOTOS agree with each other in the photon rates and distributions. The only noticeable difference is in the photon angular distribution for the process, at small angular separation from the nearest lepton. We quantify the comparison by computing relative and boson mass shifts and find them to be consistent with ≈0.7 MeV within statistical uncertainties. We conclude that a systematic uncertainty of 0.7 MeV would account for any differences in the FSR multiphoton emission between the HORACE and PHOTOS algorithms.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
The authors wish to thank Ilija Bizjak for his assistance with the HORACE program and Zbigniew Was for providing the interface to the PHOTOS program. The authors wish to thank William Ashmanskas, Franco Bedeschi, Daniel Beecher, Ilija Bizjak, Kenichi Hatakeyama, Christopher Hays, Mark Lancaster, Sarah Malik, Larry Nodulman, Peter Renton, Tom Riddick, Ravi Shekhar, Melvyn Shochet, Oliver StelzerChilton, Siyuan Sun, David Waters, Yu Zeng, and other colleagues in the CDF Collaboration for helpful discussions. They also thank Carlo Carloni Calame, Guido Montagna, Alessandro Vicini, Doreen Wackeroth, and Zbigniew Was for discussions regarding electroweak radiative corrections. They acknowledge the support of the U.S. Department of Energy, Office of High Energy Physics, and the Fermi National Accelerator Laboratory (Fermilab). The computational resources used in this study were provided by Fermilab. Fermilab is operated by Fermi Research Alliance, LLC, under Contract no. DEAC0207CH11359 with the United States Department of Energy.
References
 M. Baak, J. Cúth, J. Haller et al., “The global electroweak fit at NNLO and prospects for the LHC and ILC,” The European Physical Journal C, vol. 74, article 3046, 2014. View at: Publisher Site  Google Scholar
 T. Aaltonen, V. M. Abazov, B. Abbott et al., “Combination of CDF and D0 $W$Boson mass measurements,” Physical Review D, vol. 88, no. 5, Article ID 052018, 11 pages, 2013. View at: Publisher Site  Google Scholar
 T. Aaltonen, B. A. González, S. Amerio et al., “Precise measurement of the Wboson mass with the CDF II detector,” Physical Review Letters, vol. 108, no. 15, Article ID 151803, 2012. View at: Publisher Site  Google Scholar
 V. M. Abazov, B. Abbott, B. S. Acharya et al., “Measurement of the W Boson Mass with the D0 Detector,” Physical Review Letters, vol. 108, Article ID 151804, 2012. View at: Publisher Site  Google Scholar
 V. M. Abazov, B. Abbott, B. S. Acharya et al., “Measurement of the W boson mass with the D0 detector,” Physical Review D, vol. 89, no. 1, Article ID 012005, 2014. View at: Publisher Site  Google Scholar
 T. Aaltonen, S. Amerio, D. Amidei et al., “Precise measurement of the Wboson mass with the collider detector at Fermilab,” Physical Review D, vol. 89, no. 7, Article ID 072003, 2014. View at: Publisher Site  Google Scholar
 T. Aaltonen, A. Abulencia, J. Adelman et al., “First run II measurement of the W boson mass at the Fermilab Tevatron,” Physical Review D, vol. 77, Article ID 112001, 2008. View at: Publisher Site  Google Scholar
 U. Baur, S. Keller, and D. Wackeroth, “Electroweak radiative corrections to $W$ boson production in hadronic collisions,” Physical Review D, vol. 59, no. 1, Article ID 013002, 18 pages, 1999. View at: Publisher Site  Google Scholar
 U. Baur, O. Brein, W. Hollik, C. Schappacher, and D. Wackeroth, “Electroweak radiative corrections to neutralcurrent DrellYan processes at hadron colliders,” Physical Review D, vol. 65, no. 3, Article ID 033007, 2002. View at: Publisher Site  Google Scholar
 C. M. Carloni Calame, G. Montagna, O. Nicrosini, and M. Treccani, “Higherorder QED corrections to Wboson mass determination at hadron colliders,” Physical Review D, vol. 69, no. 3, Article ID 037301, 2004. View at: Publisher Site  Google Scholar
 C. M. Carloni Calame, G. Montagna, O. Nicrosini, and M. Treccani, “Multiple photon corrections to the neutralcurrent DrellYan process,” Journal of High Energy Physics, vol. 2005, article 019, 2005. View at: Publisher Site  Google Scholar
 C. M. Carloni Calame, G. Montagna, O. Nicrosini, and A. Vicini, “Precision electroweak calculation of the charged current DrellYan process,” Journal of High Energy Physics, vol. 2006, no. 12, article 16, 2006. View at: Publisher Site  Google Scholar
 C. M. Carloni Calame, G. Montagna, O. Nicrosini, and A. Vicini, “Precision electroweak calculation of the production of a high transversemomentum lepton pair at hadron colliders,” Journal of High Energy Physics, vol. 2007, no. 10, article 109, 2007. View at: Publisher Site  Google Scholar
 HORACE version 3.1, http://www2.pv.infn.it/~hepcomplex/horace.html.
 E. Barberio and Z. Wąs, “PHOTOS—a universal Monte Carlo for QED radiative corrections: version 2.0,” Computer Physics Communications, vol. 79, no. 2, pp. 291–308, 1994. View at: Publisher Site  Google Scholar
 A. B. Arbuzov, R. R. Sadykov, and Z. Was, “QED bremsstrahlung in decays of electroweak bosons,” The European Physical Journal C, vol. 73, no. 11, article 2625, 18 pages, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Ashutosh V. Kotwal and Bodhitha Jayatilaka. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP^{3}.