%Title: QCD Corrections to Jet Cross Sections in DIS
%Author: Erwin Mirkes and Dieter Zeppenfeld
%Published: Nucl. Phys. B (Proc. Suppl.) 51C (1996) 273-278
% hep-ph/9606332
%%%%%%%%%% espcrc2.tex %%%%%%%%%%
\documentstyle[twoside,fleqn,epscrc2]{article}
\setlength{\topmargin}{-1.8cm}
% put your own definitions here:
% \newcommand{\cZ}{\cal{Z}}
% \newtheorem{def}{Definition}[section]
\def\lsim{\mathrel{\raise.2ex\hbox{$<$}\hskip-.8em\lower.9ex\hbox{$\sim$}}}
\def\gsim{\mathrel{\raise.2ex\hbox{$>$}\hskip-.8em\lower.9ex\hbox{$\sim$}}}
\newcommand{\ttbs}{\char'134}
\newcommand{\AmS}{{\protect\the\textfont2
A\kern-.1667em\lower.5ex\hbox{M}\kern-.125emS}}
% add words to TeX's hyphenation exception list
\hyphenation{author another created financial paper re-commend-ed}
\vspace{-1cm}
\title{QCD Corrections to Jet Cross Sections in DIS\\[-1.8cm]
\thanks{
Invited talk presented by E. Mirkes.
To appear in proceedings of
"QCD and QED in Higher Orders" 1996
Zeuthen Workshop on Elementary Particle Theory, April 22-26, 1996.}
\hfill$\vtop{ \hbox{\normalsize TTP96-24}
\hbox{\normalsize MADPH-96-946}
\hbox{\normalsize June 1996}}$
}
\author{
Erwin Mirkes\address{Institut f\"ur Theoretische Teilchenphysik,
Universit\"at Karlsruhe, D-76128 Karlsruhe, Germany }
and
Dieter Zeppenfeld\address{Department of Physics,
University of Wisconsin, Madison, WI 53706, USA}}
\begin{document}
\begin{abstract}
Next-to-leading order corrections to
jet cross sections in deep inelastic scattering at HERA are
studied.
The predicted jet rates allow for a precise
determination of $\alpha_s(\mu_R)$ at HERA over a wide range of
$\mu_R$.
We argue, that the ``natural'' renormalization and
factorization scale is set by the average $k_T^B$
of the jets in the Breit frame and suggest to
divide the data in corresponding $$ intervals.
Some implications
for the determination of the gluon density and
the associated forward jet production in the low $x$ regime at HERA
are briefly discussed.
\end{abstract}
% typeset front matter (including abstract)
\maketitle
\section{Introduction}
Deep inelastic scattering (DIS) at HERA is a copious source of
multi-jet events. Typical two-jet cross sections
are in the 100~pb to few nb range and thus provide sufficiently high
statistics for precision QCD tests~\cite{exp_as}.
Clearly, next-to-leading order (NLO) QCD corrections
are mandatory on the theoretical side for such tests.
Full NLO corrections for one and two-jet production cross sections
and distributions are now available and implemented in the
$ep \rightarrow n$ jets event generator MEPJET,
which allows to analyze
arbitrary jet definition schemes and
general cuts in terms of parton 4-momenta\cite{mz1,mz2}.
A variety of topics can be studied with these tools.
1) The
determination of $\alpha_s(\mu_R)$ over a range of scales $\mu_R$
from dijet production at HERA:
The dijet cross section
is proportional to $\alpha_s(\mu_R)$ at leading order (LO), thus suggesting
a direct measurement of the strong coupling constant. However, the
LO calculation leaves the renormalization scale $\mu_R$ undetermined.
The NLO corrections substantially reduce the renormalization and
factorization scale dependencies which are present in the LO calculations
and thus reliable cross section predictions in terms of $\alpha_s(m_Z)$
(for a given set of parton distributions)
are made possible.
Clearly, a careful study of the choice
of scale in the dijet cross section
is needed in order to extract a reliable value for
$\alpha_s(M_z)$.
We will present some arguments
and studies of the scale dependence
of the NLO dijet cross section which suggest that the
average $k_T^B$ of the jets in the Breit frame provides
the ``natural'' scale for multi jet production in DIS.
Here, $(k_T^{B}(jet))^2$ is defined by
$2\,E_j^2(1-\cos\theta_{jP})$, where the subscripts $j$ and $P$
denote the jet and proton, respectively (all quantities are defined
in the Breit frame).
%We investigate the theoretical uncertainties for the
%dijet rate in different $$ bins and show,
%that these uncertainties are very small suggesting a precise
%determination of $\alpha_s()$ for different $$ bins.
2) The measurement of the
gluon density in the proton (via $\gamma g\to q\bar q$):
The gluon density can only be indirectly constrained
by an analysis of the structure function $F_2$ at HERA
\cite{exp_gluon}.
The boson gluon fusion subprocess
dominates the two jet cross section
at low $x$ and allows for a more direct measurement of the gluon density
in this regime. A first LO experimental analysis
has been presented in \cite{h1_gluon}.
NLO corrections reduce the factorization scale dependence
in the LO calculation (due to the initital state collinear factorization,
which introduces a mixture of the quark and gluon densities
according to the DGLAP (Altarelli-Parisi) evolution)
and thus reliable cross section predictions in terms of the
scale dependent parton distributions are made possible.
Some implications for the determination of the
gluon density have been discussed in \cite{mz2}.
3) Associated
forward jet production in the low $x$ regime as a signal of BFKL
dynamics:
BFKL evolution \cite{bfkl} leads to a larger cross section for
events with a measured forward jet (in the proton
direction) with transverse momentum $p_T^{lab}(j)$
close to $Q$ than the DGLAP \cite{dglap} evolution.
Clearly, next-to-leading order QCD corrections
for fixed order QCD, with Altarelli-Parisi (DGLAP) evolution,
are mandatory on the theoretical side in order to establish a
signal for BFKL evolution in the data.
We discuss these corrections in section 3.
The importance of higher order corrections
and recombination scheme dependencies
of the two jet cross sections
for four different jet algorithms
(cone, $k_T$ \cite{kt}, JADE, W) was already discussed
in \cite{mz1,mz2}.
While the higher order corrections and recombination scheme
dependencies in the cone
and $k_T$ schemes are small, very large corrections appear in the $W$-scheme.
We conclude from these studies that the
cone and $k_T$ schemes appear better suited for precision
QCD tests and will concentrate only on those schemes
in the following.
The goal of a versatile NLO calculation is to allow for an easy
implementation of an arbitrary jet algorithm and to impose any kinematical
resolution and acceptance cuts on the final state particles. This is best
achieved by performing all hard phase space integrals numerically,
with a Monte Carlo integration technique.
This approach also allows an investigation of the recombination scheme
dependence of the NLO jet cross sections. For dijet production at HERA such
a NLO Monte Carlo program is MEPJET\cite{mz1}.
The calculation is based on the phase space slicing method
and on the technique of universal ``crossing functions''
\cite{giele1,rgiele,rk}.
More details are described in \cite{mz1}
and we do not repeat them here.
This technique can also be extended to the case of massive
quark production \cite{laenen}.
An alternative technique for the calculation of NLO corrections in
jet physics is the ``subtraction'' method
\cite{catani}.
\section{Choice of the Renormalization and Factorization Scale in
Multijet Production}
Jet production in DIS is a multi-scale problem
and it is not a priori clear at which scale $\alpha_s$ is probed.
In the following we will discuss the choice of the renormalization and
factorization scale for dijet production at HERA.
As mentioned before, the NLO corrections
are expected to reduce the scale
dependencies in the LO predictions
provided the scale is of the order of the typical hardness
of the partonic process.
The following studies are done for the
cone algorithm (which is defined in the laboratory frame) and
the distance
$\Delta R=\sqrt{(\Delta\eta)^2+(\Delta\phi)^2}$ between two partons
decides whether they should be recombined to a single jet. Here the variables
are the pseudo-rapidity $\eta$ and the azimuthal angle $\phi$. We
recombine partons with $\Delta R<1$.
Furthermore, a cut on the jet transverse momenta of $p_T(j)>5$~GeV in the lab
is imposed. We employ the two loop formula for the strong coupling constant
with %${\Lambda_{\overline{MS}}}={\Lambda_{\overline{MS}}^{(4)}}=230$ MeV
a value for ${\Lambda_{\overline{MS}}^{(4)}}$
consistent with the value from the parton distribution functions.
The value of $\alpha_s$ is matched at the thresholds $\mu_R=m_q$ and the
number of flavors is fixed to $n_f=5$ throughout, {\it i.e.} gluons are
allowed to split into five flavors of massless quarks.
A running QED fine structure constant $\alpha(Q^2)$ is used.
In addition the following set of kinematical
cuts is imposed on the initial virtual photon
and on the final state electron and jets:
We require
40~GeV$^210$~GeV on the scattered
electron, and a cut on the pseudo-rapidity $\eta=-\ln\tan(\theta/2)$
of the scattered lepton and jets of $|\eta|<3.5$.
First studies of the scale dependence of the dijet cross section in the cone
scheme are presented in \cite{mz1}.
We have considered scales related to the
scalar sum of the parton transverse momenta in the Breit frame,
$\sum_i \,p_T^B(i)$, and the virtuality $Q^2$ of the incident photon.
In the following we will also consider scales related to
$\sum_i \,k_T^B(i)$.
Here, $(k_T^{B}(i))^2$ is defined by
$2\,E_i^2(1-\cos\theta_{iP})$, where the subscripts $i$ and $P$
denote the final parton (or jet) and proton,
respectively (all quantities are determined
in the Breit frame).
%The Breit frame is characterized by the property that
%the momentum of the exchanged photon $q_\mu$ is completely spacelike,
%having nonzero energy component
The Breit frame is characterized by the vanishing energy component of the
momentum of the exchanged photon. Both the photon momentum
\begin{equation}
q=(0,0,0,-2xE),\,\,\,-q^2=Q^2=4x^2E^2\\
\end{equation}
and the proton momentum
\begin{equation}
P=E(1,0,0,1)
\end{equation}
are chosen along the $z$-direction.
$x$ is the standard Bjorken scaling variable.
In the parton model, the incoming quark
with momentum $p=xE(1,0,0,1)$ collides elastically with the virtual
boson and is scattered
in the opposite direction (the 'current hemisphere')
with momentum $p'=xE(1,0,0,-1)$;
therefore, $(k_T^{B}(p'))^2= Q^2$ in the limit of the quark parton model,
whereas $p_T^B(p')=0$.
The kinematics for dijet production is more complex: the momentum fraction
$\eta$ of the incoming parton must be larger than $x$ since
$m_{jj}^2= Q^2(\eta/x-1)$ and, in general,
the jets have a nonvanishing transverse momentum
with respect to the $\gamma^*$-proton direction.
At LO, i.e. for massless jets, the
relation between $k_T^B(j)$ and $p_T^B(j)=E_j\sin\theta_{jP}$ reads:
\begin{equation}
k_T^B(j) = p_T^B(j)\sqrt{\frac{2}{1+\cos\theta_{{\,\!jP}}}}
\end{equation}
Obviously, $k_T^B(j)>p_T^B$, and one can also show that
$\sum_j \,k_T^B(j) > Q$. Thus, $\sum_j \,k_T^B(j)$ is approximately
given by the harder of the two scales $Q$ and $\sum_j \,p_T^B(j)$
\cite{mz3}.
For large dijet invariant masses (i.e. for ``true'' two jet kinematics)
one has $\eta >> x$, the dijet system will
be strongly boosted in the proton direction
and $k_T^B(j) \approx p_T^B(j)>>Q/2$,
as long as very small scattering angles in
the center of mass frame are avoided.
For $m_{jj}\approx Q$ or smaller, on the other hand,
$\sum_j \,k_T^B(j)\approx Q$ and, typically, both are
considerably larger than $\sum_j \,p_T^B(j)>m_{jj}$, which corresponds to the parton model limit).
Thus $\sum_j \,k_T^B(j)$ smoothly interpolates between the correct limiting
scale choices, it approaches $Q$ in the parton limit and it corresponds to
the jet
transverse momentum when the photon virtuality becomes negligible.
It appears to be the ``natural'' scale for multi jet production in DIS.
In \cite{mz1} we found that the scale dependence of the dijet cross section
does not markedly improve in NLO for $\mu^2=\xi Q^2$.
This is also shown in Fig.~1 (dotted curves)
where the dependence of the two-jet
cross section on the scale factor $\xi$ is shown.
We used the parton distribution functions set MRS D-$^\prime$
\protect\cite{mrs}
and employed the two loop formula for the strong coupling constant
both in the LO and NLO curves.
For scales related to $\sum_i \,p_T^B(i)$ the uncertainty from the variation
of the renormalization and factorization scale is markedly reduced
compared to the LO predictions (dashed curves in Fig.~1).
%For scales related to $\sum_i \,p_T^B(i)$,
Here $\xi$ is defined via
%
%
\begin{equation}
\mu_R^2 = \mu_F^2 = \xi\;(\sum_i \,p_T^B(i))^2\,.
\label{xidef}
\end{equation}
%
%
%and similar for $(\sum_i \,k_T^B(i))^2$.
The resulting $\xi$ dependence for
$ \mu_R^2 = \mu_F^2 = \xi\;(\sum_i \,k_T^B(i))^2$
is shown as the solid lines in Fig.~1.
In this case, the NLO two-jet cross section is essentially independent on
$\xi$ for $\xi<2$.
Hence, the theoretical uncertainties due to the scale variation
are very small suggesting a precise
determination of $\alpha_s()$ for different $$ bins,
where
\begin{equation}
=\frac{1}{2}\,\, (\sum_{j=1,2} \,k_T^B(j))
\label{akt}
\end{equation}
%
%
\begin{figure}[tp]
\vspace{5.5cm}
\begin{picture}(7,7)
%\special{psfile=scale.ps voffset=-150 hoffset=-10 hscale=40
\special{psfile=scale1.ps voffset=-170 hoffset=-80 hscale=60
vscale=60 angle=0}
\end{picture}
\caption{
Dependence of the two-jet exclusive cross
section in the cone scheme
on the scale factor $\xi$.
The dashed curves are for $\mu_R^2=\mu_F^2=\xi\;(\sum_i\;p_T^B(i))^2$.
Choosing $(\sum_i\;k_T^B(i))^2$ as the basic scale yields the solid curves.
Choosing $Q^2$ as the basic scale yields the dotted curves.
Results are shown
for the LO (lower curves) and NLO calculations.
}
\end{figure}
\begin{figure}[tp]
\vspace{5.5cm}
\begin{picture}(7,7)
%\special{psfile=scale.ps voffset=-150 hoffset=-10 hscale=40
\special{psfile=fig_kt.ps voffset=-170 hoffset=-80 hscale=60
vscale=60 angle=0}
\end{picture}
\caption{
$$ distribution for the two-jet exclusive cross section.
The parameters are explained in the text.
}
\end{figure}
\begin{figure}[tp]
\vspace{5.5cm}
\begin{picture}(7,7)
%\special{psfile=scale.ps voffset=-150 hoffset=-10 hscale=40
\special{psfile=fig_pt.ps voffset=-170 hoffset=-80 hscale=60
vscale=60 angle=0}
\end{picture}
\caption{
$$ distribution for the three $$ bins shown in Fig.~1.
}
\end{figure}
\begin{figure}[tp]
\vspace{5.5cm}
\begin{picture}(7,7)
%\special{psfile=scale.ps voffset=-150 hoffset=-10 hscale=40
\special{psfile=fig_qq.ps voffset=-170 hoffset=-80 hscale=60
vscale=60 angle=0}
\end{picture}
\caption{
$Q^2$ distribution for the three $$ bins shown in Fig.~1.
}
\end{figure}
Fig.~2 shows the $$ distribution for the NLO exclusive dijet
cross section.
We used the parton distribution functions set GRV \cite{grv}
and $\mu_R^2=\mu_F^2=1/4\;(\sum_i\;k_T^B(i))^2$.
In addition to the cuts imposed in Fig.~1,
we require $p_T^B>4$~GeV for each jet
%and $k_T^B>4$ GeV in the Breit frame,
and $\sum_j k_T^B(j)> 10$ GeV.
We have divided the NLO cross section in the following three
$$ bins.\\
bin 1: 5~GeV $<\,\,\,\,\,\, \,\,<$ 10~GeV \\
bin 2: 10~GeV $<\,\,\,\, \,\,\,\,<$ 15~GeV \\
bin 3: 15~GeV $<\,\,\,\,$.\\
Table~1 shows the corresponding NLO cross sections
for these three bins for two different sets of parton distributions
and for different values for the scale factor $\xi$.
The theoretical uncertainties
for the
NLO dijet cross section from these numbers
are very small, in particular for the first two bins.
\begin{table}[t]
\caption{Jet cross sections in pb for the three $$ bins.
}\label{table1}
\vspace{2mm}
\begin{tabular}{lccc}
\hspace{0.8cm}
& \mbox{bin 1 }
& \mbox{bin 2 }
& \mbox{bin 3 } \\
\hline\\[-3mm]
\mbox{GRV} $\xi=1/4 $ & 497 & 320 & 146 \\
\mbox{GRV} $\xi=1/16$ & 504 & 306 & 134 \\
\mbox{GRV} $\xi=1 $ & 488 & 322 & 151 \\
\mbox{MRSD-'} $\xi=1/4 $ & 487 & 322 & 148 \\
\end{tabular}
\end{table}
Fig.~3 and 4 show the $$ and $Q^2$
distribution for the NLO exclusive dijet
cross section for these three bins.
Whereas the $$ and $Q^2$ distributions are fairly similar
to the $$ distribution for the lowest bin, large differences
are found for the other bins.
This partly explains the rather different scale dependence observed
in Fig.1.
\section{Forward Jet Production in the Low $x$ Regime}
Deep inelastic scattering with a measured forward jet
with relatively large momentum fraction $x_{jet}$
(in the proton direction) and $p_T^{2\,lab}(j)\approx Q^2$
is expected to provide
sensitive information about the BFKL dynamics at low $x$
\cite{mueller,allen1}.
In this region there is not much phase space
for DGLAP evolution with transverse momentum ordering,
whereas large effects are expected for BFKL evolution in $x$.
In particular, BFKL evolution is expected to substantially enhance cross
sections in the region $x<8~$GeV$^2$ , $x<0.004$,
$0.1 < y < 1$, an energy cut of $E(e^\prime)>11$~GeV on the scattered
electron, and a cut on the pseudo-rapidity $\eta=-\ln\tan(\theta/2)$
of the scattered lepton of
$ -2.868 < \eta(e^\prime)< -1.735$
(corresponding to $160^o < \theta(l^\prime) < 173.5^o$).
Jets are defined in the cone scheme (in the laboratory frame) with
$\Delta R = 1$ and $|\eta(j)|<3.5$.
We require a forward jet with $x_{jet}=p_z(j)/E_{P} > 0.05$,
$E(j)>25$ GeV, $0.50.05$) and restricted transverse momentum ($0.50.05$ requirements. For $x<> Q^2\; ,
\end{eqnarray}
where $m_T$ and $y$ are the transverse mass and rapidity of the
partonic recoil system, respectively. Thus a recoil system with substantial
transverse momentum and/or invariant mass must be produced and this
condition favors recoil systems composed out of at least two additional
energetic partons.
As a result
one finds very large fixed order perturbative QCD corrections (compare
2 jet LO and NLO results with a forward jet in Table~\ref{table2}).
In addition, the LO $({\cal{O}}(\alpha_s^2))$ 3-jet cross section
is larger than the LO $({\cal{O}}(\alpha_s))$
2-jet cross section.
Thus, the forward jet cross sections in Table~\ref{table2} are dominated
by the $({\cal{O}}(\alpha_s^2))$ matrix elements.
The effects
of BFKL evolution must be seen and isolated on top of these fixed order QCD
effects. We will analyze these effects in a subsequent publication.
\section{Conclusions}
The calculation of NLO perturbative QCD corrections has received an enormous
boost with the advent of full NLO Monte Carlo programs
\cite{giele1,jim}. For dijet production at HERA the NLO Monte Carlo
program MEPJET \cite{mz1} allows to study jet cross sections for arbitrary
jet algorithms. Internal jet structure, parton/hadron recombination effects,
and the effects of arbitrary acceptance cuts can now be simulated at the full
${\cal O}(\alpha_s^2)$ level.
We found large NLO effects for some jet definition schemes (in particular the
$W$-scheme) and cone and $k_T$ schemes appear better suited for precision
QCD tests.
The extraction of gluon distribution functions is now supported by
a fully versatile NLO program. Preliminary studies show that large NLO
corrections are present in the Bjorken $x$ distribution for dijet events,
while these effects are mitigated in the reconstructed Feynman $x$ ($x_i$)
distribution, thus aiding the reliable extraction of $g(x_i,\mu_F^2)$.
For the study of BFKL evolution by considering events with a forward
``Mueller''-jet very large
QCD corrections are found at ${\cal O}(\alpha_s^2)$.
These fixed order effects form an important background to the observation
of BFKL evolution at HERA. They can now be studied systematically and for
arbitrary jet algorithms.
This research was supported by the University of Wisconsin
Research Committee with funds granted by the Wisconsin Alumni Research
Foundation and by the U.~S.~Department of Energy under Grant
No.~DE-FG02-95ER40896. The work of E.~M. was supported in part
by DFG Contract Ku 502/5-1.
\begin{thebibliography}{9}
\bibitem{exp_as}
H1 Collaboration, T.~Ahmed et. al.,
Phys. Lett. {\bf B346} (1995) 415;\\
ZEUS Collaboration, M. Derrick et al.,
Phys. Lett. {\bf B363} (1995) 201.
\bibitem{mz1}
E. Mirkes and D. Zeppenfeld,
TTP95-42, MADPH-95-916, (1995), hep-ph/9511448.
\bibitem{mz2}
E. Mirkes and D. Zeppenfeld,
TTP96-10, MADPH-96-935, (1996), hep-ph/9604281.
\bibitem{exp_gluon}
ZEUS Collaboration, M. Derrick et al.,
Z. Phys. {\bf C65} (1995) 379;\\
{\it ibid.} Phys. Lett. {\bf B345} (1995) 576;\\
H1 Collaboration, T.~Ahmed et. al.,
Nucl. Phys. {\bf B439} (1995) 471;\\
{\it ibid.} Phys. Lett. {\bf B354} (1995) 494.
\bibitem{h1_gluon}
H1 Collaboration, T.~Ahmed et. al.,
Nucl. Phys. {\bf B449} (1995) 3.
\bibitem{bfkl}
E.A. Kuraev, L.N. Lipatov and V.S. Fadin, {\bf JETP 45} (1972) 199;\\
Y.Y. Balitsky and L.N. Lipatov,
Sov. J. Nucl. Phys. {\bf 28} (1978) 282.
\bibitem{dglap}
G. Altarelli and G. Parisi, Nucl. Phys. {\bf 126} (1977) 297;\\
V.N. Gribov and L.N. Lipatov,
Sov. J. Nucl. Phys. {\bf 15} (1972) 438 and 675;\\
Yu. L. Dokshitzer, Sov. Phys. {\bf JETP 46} (1977) 641.
\bibitem{kt}
S.~Catani, Y.L.~Dokshitzer and B.R.~Webber,
Phys. Lett. {\bf B285} (1992) 291.
\bibitem{giele1}
W.~T.~ Giele and E.~W.~N.~Glover,
Phys. Rev. D {\bf 46} (1992) 1980.
W.~T.~ Giele, E.~W.~N.~Glover and D.A. Kosower,
Nucl. Phys. {\bf B403} (1993) 633.
\bibitem{rgiele} W. Giele, these proceedings.
\bibitem{rk} D. Kosower, these proceedings.
\bibitem{laenen} E. Laenen, these proceedings.
\bibitem{catani} Z. Kunszt, these proceedings;\\
S. Catani, these proceedings.
\bibitem{mz3}
E. Mirkes and D. Zeppenfeld,
in preparation.
\bibitem{mrs}
A.D. Martin, W.J. Stirling and R.G. Roberts,
Phys. Lett. {\bf B306} (1993) 145.
\bibitem{grv}
M. Gl\"uck, E. Reya and A. Vogt,
Z.~Phys.~{\bf C67} (1995) 433.
\bibitem{mueller}
A.H. Mueller, Nucl. Phys. {\bf B} (Proc. Suppl.) {\bf 18C} (1990) 125;
J. Phys. {\bf G17} (1991) 1443.
\bibitem{allen1}
J. Kwiecinski, A.D. Martin and P.J. Sutton, Phys. Rev. {\bf D46}
(1992) 921;\\
J. Bartels, A. De Roeck and M. Loewe,
Z. Phys. {\bf C54} (1992) 635;\\
W.K. Tang, Phys. Lett. {\bf B278} (1992) 363.
\bibitem{deroeck}
The cuts are similar to the cuts used in the present H1 forward
jet analyses. We thank A. De Roeck and E.M. Mroczko for this information.
\bibitem{jim}
H.~Baer, J.~Ohnemus, and J.~F.~Owens,
Phys. Rev. {\bf D40} (1989) 2844;
Phys. Lett. {\bf B234} (1990) 127;
Phys. Rev. {\bf D42} (1990) 61.
\end{thebibliography}
\end{document}