Radiative Corrections to Bremsstrahlung in Radiative Return

Scott A. Yost Baylor University

with C. Glosser, S. Jadach and B.F.L. Ward 1

Radiative Corrections to Bhabha Scattering

In the 1990s, S. Jadach, M. Melles, B.F.L. Ward and S.A. Yost calculated all two photon real and virtual corrections to the small angle Bhabha scattering process used in the luminosity monitor at LEP and SLC.

These corrections were used to bring the theoretical uncertainty in the luminosity measurement, as calculated by the BHLUMI Monte-Carlo program, to within a

0.06%

precision level.

The same amplitudes, with crossing, could be used to calculate the virtual photon corrections to another important Standard Model process, fermion pair production: We include single-photon hard Bremsstrahlung.

2

Radiative Corrections for Radiative Return

Radiating a photon from the initial state makes a convenient way to study low energy final states in a high-energy collider: a wide range of energies can be studied by varying the energy of the radiated photon.

The radiative corrections needed to calculate the cross-sections for radiative return are identical to those calculated in the fermion pair process. A common application is to pion pairs, for example.

Kuhn, Rodrigo, and collaborators have published a series of papers and a Monte Carlo program, PHOKHARA, to compute this process. The expressions they derived for the order a 2 corrections to single hard bremsstrahlung provide an important cross-check for the same radiative corrections computed earlier for ISR radiation in fermion pair production.

3

Radiative Corrections to Fermion Pair Production

Representative Feynman diagrams for the process calculated are shown here.

k

, s

p

1 , l 1

p

3 , l 3

p

4 , l 4

p

2 , l 2 4

Computational Method

The graphs shown were calculated by S. Jadach, M. Melles, B.F.L. Ward and S. Yost in Phys. Rev. D65 , 073030 (2002), based on earlier results for the corresponding

t

-channel graphs by the same authors,

Phys. Lett.

B377

, 168 (1996).

The results were obtained using • Helicity spinor methods • Vermaseren’s algebraic manipulation program FORM • Oldenborgh’s FF package of scalar one-loop Feynman integrals (later replaced by analytic expressions) • Mass corrections added via methods of Berends

et al

(CALCUL Collaboration), which were checked to show that all significant collinear mass corrections were included 5

Assumptions Made In The Calculation

A small angle approximation was made, which suppresses box diagrams of the type shown below. The addition of such diagrams is in progress, since they will be useful if the proposed NLC is built.

The lightness of the fermions permits masses to be neglected except when the multiply denominators which vanish in collinear limits. This was used to eliminate many very small mass terms from the calculation.

6

The Complete Result

The differential cross section for initial state radiation (ISR) can be expressed as

r i

= 2

p i .

k

and the differential cross section with virtual corrections can be expressed as 7

The Complete Result

The complete amplitude with virtual corrections can be expressed in terms of the amplitude without virtual corrections (pure hard Bremsstrahlung) as where

I

1 and

I

2 are spinor factors which vanish in the collinear limits

p i .

k

= 0 : Spinor product: |s l (p,k)| 2 k = 2p . with

p ij

=

p i

or

p j

when s = l

i

or l

j

.

8

For helicities s

Form Factors

= l 1 . Otherwise interchange

r

1 and

r

2 : (mass terms added later)

r i

= 2

p i .

k v

=

r

1 +

r

2

z

= 1 -

v

9

Special Functions

The form factors depend on the following special functions, designed to be stable in the collinear limits (small

r i

) which are most important in the cross section.

with Dilogarithm(Spence function) 10

Next to Leading Log Approximation

For Monte Carlo use, it is desirable to have the shortest possible equation with sufficient accuracy. For most of the range of hard photon energies, the leading log ( LL ) and next to leading log ( NLL ) contributions suffice. These include all terms important in the collinear limits (

r i

small).

• • To NLL order, the spinor terms

I

1 and

I

2 can be dropped, since most logarithmically divergent for small

r i . f

1 and

f

2 are at The spin-averaged collinear limit of

f

0 is ...

big log mass corrections 11

Mass Corrections

Mass corrections were added following Berends,

et al

(CALCUL collaboration). We checked that all significant mass corrections are obtained in this manner.

The most important corrections for a photon with momentum

k

radiated collinearly with each incoming fermion line

p

1 and

p

2 are added via the prescription At the cross-section level, the net effect is that the spin-averaged form factor

f

0 receives an additional mass term 12

YFS Residuals

The Monte Carlo program will calculate YFS residuals, which are obtained by subtracting the YFS factors containing the infrared singularities. This amounts to subtracting a term from the form factor

f

0 .

At NLL order, we would have In our comparisons we will actually subtract this NLL term and look at the NNLL behavior of each expression.

13

IN

Comparisons

Igarashi and Nakazawa, Nucl. Phys. B288 (1987) 301

-

spin-averaged cross section, fully differential in

r

1 and

r

2 , no mass corrections

BVNB

Berends, Van Neerven and Burgers, Nucl. Phys. B297 (1988) 429 - spin-averaged cross section, differential only in

v

=

r

1 +

r

2 , includes mass corrections

KR

Kuhn and Rodrigo, Eur. Phys. J. C25 (2002) 215 - spin-averaged Leptonic tensor, fully differential in

r

1 and

r

2 , includes mass corrections The

KR

comparison is new, and closest to our calculation in its assumptions.

14

The New Comparison

The new comparison is to the leptonic tensor of Kuhn and Rodrigo, which was constructed for radiative return in hadron production, but can be adapted to fermion pairs by changing the final state tensor. The ISR result is expressed as a “ leptonic tensor ” with Leptonic tensor which can be contracted with a final state tensor to get the squared matrix element: 15

The New Comparison

The coefficients in the leptonic tensor can be separated into a tree-level term and a one-loop correction, An infrared term is subtracted to get a finite result: with cutoff on Calculating the coefficients requires considerable attention to numerical stability!

16

Analytical Comparison

The expression for the leptonic tensor is very different from our earlier exact result, but it is possible to show that in the massless NLL limit, they agree analytically. The virtual correction to the YFS residual can be expressed in the NLL limit as with coefficients We find that 17

Results of a KK Monte Carlo run with 10 8 events at

E

CMS = 200 GeV.

This figure shows the complete real + virtual photon radiative correctionto muon pair production.

The standard YFS infrared term 4 p B YFS has been subtracted to create a finite result.

Monte Carlo Results

Real + Virtual Photon Correction

Maximum hard photon energy fraction 18

Results of a KK Monte Carlo run with 10 8 events at

E

CMS = 200 GeV.

This figure shows only the pure virtual photon correction to single hard bremmstrahlung.

The standard YFS infrared term 4 p B YFS has been subtracted to create a finite result.

Monte Carlo Results

Pure Virtual Correction

19

Results of a KK Monte Carlo run with 10 8 events at

E

CMS = 200 GeV.

This figure shows the next to leading log (NLL) contribution to the real + virtual photon cross section.

The leading log (LL) contribution has been subtracted from each expression.

Monte Carlo Results

NLL

20

Results of a KK Monte Carlo run with 10 8 events at

E

CMS = 200 GeV.

This figure shows the sub-NLL contribution to the real + virtual photon correction to muon pair production The NLL expression of JMWY has been subtracted in each case to reveal the NNLL contributions.

Monte Carlo Results

NNLL

21

Summary

• The size of the NNLL corrections for all of the compared “exact” expressions is less than 2 x 10 -6 in units of the Born cross section for photon energy cut

v

max < 0.75.

• For

v

max < 0.95 (5 Gev photon), all the results except BVNB agree to within 2.5

x 10 -6 of the Born cross ection.

• For the final data point,

v

max = 0.975, the KR and JMWY results differ by 3 x 10 -5 of the Born cross ection.

• These comparisons show that we have a firm understanding of the precision tag for an important part of the order fermion pair production in precision studies of the final LEP2 data analysis, radiative return at F a 2 corrections to and B-factories, and future NLC physics.

22