\documentstyle[12pt,epsf]{article}
%\documentstyle[aps,epsfig,preprint]{revtex}
\setlength{\topmargin}{.2in}
\setlength{\oddsidemargin}{0in} 
\setlength{\textwidth}{6.5in}
\setlength{\textheight}{9.0in}
%\textwidth              16.4cm
%\oddsidemargin           2.5cm
%\advance\oddsidemargin  by -0.9in
%\evensidemargin          0.0cm
%\advance\evensidemargin by -0.9in
%\marginparwidth          1.9cm
%\marginparsep            0.6cm
%\marginparpush           0.6cm
%\topmargin              0cm
%\advance\topmargin      by -0.20in
%\textheight             24.00cm


\newcommand{\etal}{{\it et al.}}

\def\Journal#1#2#3#4{{#1} {\bf #2}, #3 (#4)}

%% Some useful journal names

\def\NCA{Nuovo Cimento}
\def\NIM{Nucl. Instrum. Methods}
\def\NIMA{{Nucl. Instrum. Methods} A}
\def\NPB{{Nucl. Phys.} B}
\def\PLB{{Phys. Lett.}  B}
\def\PRL{Phys. Rev. Lett.}
\def\PRD{{Phys. Rev.} D}
\def\ZPC{{Z. Phys.} C}
\def\NPV{{Nucl. Phys.} V}


\begin{document}

\title{\boldmath Prospects for Searching for Excited
Leptons during Run~II of the Fermilab Tevatron}

\maketitle
%%author list commands

%\author


\begin{centering}
{
  E. Boos$^1$,  A. Vologdin$^1$, D. Toback$^2$, and J. Gaspard$^2$\\
  \hfil\\
  $^1$Institute of Nuclear Physics, Moscow State University,
      119899, Moscow, Russia \\
  $^2$Texas A\&M University, College Station, TX 77843  USA}

\end{centering}


\begin{abstract}
This letter presents a study of prospects of searching for excited
leptons during Run~II of the Fermilab Tevatron.  We concentrate on
single and pair production of excited electrons in
the photonic decay channel in one CDF/D$\O$ detector equivalent for 
2 fb$^{-1}$.  By the end of Run~IIa, the limits should be easily extended 
beyond those set by LEP and HERA for excited leptons with mass above about 
190 GeV.

\end{abstract}


%\begin{flushright}
%\pacs{PACS numbers: xx.xx.Ly, yy.yy.Jv, zz.zz.Rm, ww.ww.Pb}
%\end{flushright}

\clearpage


%\section{Introduction}

The Standard Model (SM) of particle physics is known to give results that 
match the current experimental data with a high precision.
However, because of well known theoretical problems and disadvantages,
it is widely believed it can not be a complete theory of elementary 
particles, but rather, likely a kind of effective theory at energies 
below some scale of the order of a TeV.  While many models of ``new physics" 
beyond the SM have been suggested over the years, one of the most 
straight-forward ideas proposes that quarks and leptons are composite 
particles.
Such a scenario is a recurring theme in nature;  molecules are
composed of  atoms, which are in turn composed of nuclei, which are in
turn composed of  nucleons and so on.  Furthermore, composite models 
can explain, in principle, family replication, mixing in the quark 
and lepton sectors, and even make the fermion masses and weak mixing 
angles calculable.

In most composite models fermions possess an
underlying substructure which is characterized by a
scale $\Lambda$~\cite{Other} with  $\Lambda$ about 1 TeV or higher~\cite{g-2}.
While there is no unified model of compositeness, a model-independent
effective Lagrangian for excited leptons, originally proposed
in~\cite{Hagiwara}, can be used to model single and pair production in experiments.  The Lagrangian:

\begin{eqnarray}\label{EffLagrangian}
L_{l^{*}l^{*}}  &=&  \overline{l}^{*}\gamma^{\mu}(\frac{\tau}{2}W_{\mu}+g^{'}\frac{Y}{2}B_
{\mu})l^{*}\\
L_{l^{*}l}  &=&  \frac{1}{2\Lambda} \overline{l}^{*}_{R} \sigma^{\mu\nu}
(fg\frac{\tau}{2}W_{\mu\nu} + f^{'}g^{'}\frac{Y}{2}B_{\mu\nu} ) l_{L}
+h.c.
\end{eqnarray}


\noindent
where $f$ and $f'$ are coupling constants, has been used extensively in a number of phenomenological papers which 
presents ideas on searching for excited fermion production and decay to 
final state gauge bosons in $e^+e^-$, $p\bar{p}$, $ep$, and $e \gamma$ 
collisions~\cite{Boudjema2,Baur,Boos}.  Direct searches for lepton 
compositeness has been done extensively at  LEP~\cite{LEP} and 
HERA~\cite{HERA}, each with no discovery, but with ever more sensitive limits. 
Unfortunately, only direct searches for quark compositeness have 
been done at the Fermilab Tevatron~\cite{Tevatron Quark}.  

%Each imposes
%constraints on $\frac{f}{\Lambda}$ as a function of
%the excited lepton mass, $m^*$. Searches for single production have produced
%limits
%****
%in $\frac{f}{\Lambda}$ and $m^*$ parameter plane.
%****
% with a comparison of LEP and HERA results is shown on Fig~\ref{LEPHERA:Fig1}.

%???? I do not see this Fig. Alexander, could you put it in????

%Pair production at LEP give a more coupling-independent
%limits which exclude excited leptons with mass $m^*$
%lower than 90 GeV, see
%****
%%Fig~\ref{LEPHERA:Fig1}.
%****


%As far as we know searches for
%excited leptons have never been carried out
%at Tevatron.
%***

In this Letter we present a study for searching for excited lepton production
and decay for the upcoming Run~II of CDF and D$\O$ detectors.
We begin with an updated simulation of excited lepton production, and optimize
the sensitivity by studying the kinematics of excited lepton production and 
decay.  With these results we estimate the mass reach and compare to the 
recent results from LEP and HERA. 
%of the prospects an updated simulation of excited lepton production and 
%propose a new search strategy, as well as a study of the kinematics of excited
%lepton production and decay and estimate the mass reach for the upcoming 
%Run~II of CDF and D$\O$ detector.
  %As we will show, searching for excited leptons at the Tevatron has many advantages; background signatures are much cleaner than those of previous mentioned colliders, and the Tevatron has the potential to greatly extend the mass search range.


%for compositeness at the Fermilab
%Tevatron by looking for excited leptons,
%which in contrast to the compositeness sequences like quark contact interactions
%and excited quarks was not carried out at the Tevatron before.

%****


%\section{Simulations and tools}

To study the mass reach of the Tevatron, we simulated single and pair
production and decay of excited leptons using the upgraded Fermilab accelerator
(1.8 $\rightarrow$ 2.0 TeV), and CDF and D$\O$ detectors for Run~II.  The 
Feynman rules from the effective Lagrangian (Eqn.~\ref{EffLagrangian}) are
implemented into {\sc CompHEP}~\cite{CompHEP} using the {\sc LanHEP}~\cite{LanHEP} 
software package.
We have included into this simulation a complete tree-level calculation
which takes into account all the spin correlations between excited
states production, subsequent decays, and incorporate the known NLO 
corrections.  All the partial width and known 2$ \rightarrow$ 2 
cross-section have been cross-checked at the symbolic level.  Events 
on parton level generated by means of {\sc CompHEP} have been used as an 
external process for {\sc pythia} with the help of the {\sc CompHEP}-{\sc pythia} 
interface \cite{interface}.  The underlying event, jet fragmentation, 
ISR and FSR are modeled using the {\sc pythia}~\cite{Pythia} Monte Carlo 
with the {\sc cteq4l}~\cite{CTEQ4L} structure functions.
Since Drell-Yan production of excited leptons production 
is similar to that of Supersymmetric leptons, we take the K-factors 
 given in \cite{K-signal}, which only depend on  the  masses of 
final particles, and vary between 1.23 and 1.24 in the mass range 
of the search.

While excited leptons can come in three flavors: $e^{*}$, $\mu^{*}$ and $\tau^{*}$, 
we chose to concentrate on the electron since the result of 
$\mu^{*}$ expected to be similar to $e^{*}$, and $\tau$'s at the Tevatron 
are still difficult to trigger on and identify.  
Excited electrons can decay via:
 $e^{*} \rightarrow e \gamma $, $e^{*}
\rightarrow e Z $, and $ e^{*}\rightarrow W \nu$ channels with the branching 
ratios shown in Fig.~\ref{BranchingRatios}.  The backgrounds to the $W$ and $Z$ 
decay channels in the hadronic channels are fairly large and the leptonic  branching fractions needed to
identify the $W$ and $Z$ channels are small; however, the backgrounds to the 
photonic final states are relatively small in comparison making them a 
gold-plated signature.

In the case of the single excited electron production mode, there are two
possible signal 
$p \overline{p} \rightarrow  e^*e \rightarrow e^+ e^- \gamma $ and
$p \overline{p} \rightarrow e^*\nu \rightarrow e \nu \gamma$.
Similarly, pair production gives
$p \overline{p} \rightarrow  e^{*+}e^{*-} \rightarrow e^+ e^- \gamma\gamma $ 
and
 $p \overline{p} \rightarrow e^*\nu^* \rightarrow e \nu \gamma\gamma$.
For simplicity we concentrate on inclusive $ee\gamma$ final state for 
both single and pair production as it covers most of the production in both 
single and pair production cases.  Before taking into account the branching 
ratio to photons, the total production cross section is shown in 
Fig.~\ref{CrossSections} for the case $f/\Lambda = 10^{-2}~{\rm GeV}^{-1}$ and 
$f=f^{'}$.

There are a number of backgrounds to the $ee\gamma$ channel.  The dominant 
backgrounds are $W\gamma$+jets and $Z\gamma$+jets production.  Others include 
$W$+jets and $Z$+jets, and multijets, where jets can fake leptons and/or 
photons.  Studies have shown~\cite{DZero Diboson Studies} that the fake 
backgrounds can be modeled by the kinematics of the irreducible backgrounds.
The backgrounds are modeled using the same {\sc CompHEP} simulation structure 
described above.
Since there is an infrared singularity at $p_T^{\gamma} = 0$,
we require  $ P_{T}^{\gamma}> 10$ GeV and $\Delta R_{ij}>0.1 $, where 
$\Delta R_{ij}\equiv \sqrt{(\Delta\eta)^{2}+(\Delta\varphi)^{2}}$, and  
at the generator level, $i$ and $j$ are any lepton-photon combination.  
We take the K-factor, which has a value of about 1.36 depending on 
$p_{T}^{\gamma}$, for the background processes from the literature 
\cite{K-background}.


For both signal and backgrounds, we use a parametric simulation to model
the detector response. The SHW detector simulation~\cite{SHW} has been 
shown to be an effective averaging between the CDF~\cite{CDF} and 
D\O~\cite{DZero} detectors for Run~II.  After detector simulation the kinematic
distributions for both the signal and estimated backgrounds are shown 
in Figs~\ref{Kine1} and~\ref{Kine2}.  We find acceptances at about 
the 0.3 level for signal.

%\section{ Limits and Exclusion plot}

To maximize our sensitivity we assume, for simplicity, that taking a set 
of cuts which minimizes the expected cross section limit also
maximizes our sensitivity to calculate these limits. To calculate these limits, we use the signal 
acceptance and background estimates from the simulations above. We add in a 
single detector with 2 fb$^{-1}$ of luminosity and we assume the experimental 
parameters of Table~\ref{Parameters} for systematic uncertainties.  We use a 
frequentist method~\cite{Frequentest} to incorporate the errors and assume 
that all errors on  acceptance, background and luminosity are uncorrelated.  
To determine the expected limit as a function of a given cut using Run~II data,
under the assumption of no observed signal, the 95\% confidence level (C.L.) cross 
section upper limit, for a given cut, is uniquely determined by number of events 
observed in the data.  With no signal, the number of events observed 
in the data, $N$, fluctuates around the mean number of expected events, $M$, 
according to Poisson statistics.  Thus, an average expected cross section 
limit is given by:

\begin{equation}
 \label{Bruce's Relation}
   \sigma_{95}^{cut}=\sum_{N=0}^{\infty} \sigma_{95}(N,cut)\frac{ e^{-M}M^{N}}{N!}
\end{equation}


\begin{table}[tb]
\centering
 \begin{tabular}{|c|c|}
  \hline
       Luminosity (fb$^{-1}$)&  2.0 $\pm$ 0.1\\
  \hline
        Background error  $\frac{\delta N_{back}}{N_{back}}$  &  0.1\\

  \hline
        Acceptance error  $\frac{\delta A}{A}$  &  0.1\\

  \hline
 \end{tabular}

 \caption{Parameters used to calculate the 95\% C.L.  cross section upper limits.}
 \label{Parameters}
\end{table}


 
%where $M$ is the number of background events expected from SM backgrounds and $N$ is the number of observed events.  
\noindent In this way we get an average expected cross section limits as a function of each cut. We can then optimize for each mass
point as a function of the cuts.

After a study of multi-dimensional distributions for background and signal, we determined that the $e^+e^-\gamma$ final state has the best signal to background ratio and kinematical distributions. Furthermore, we find that applying a single cut on the kinematical variable $M_{e\gamma}$
 (where $e$ is the electron with highest $p_{T}$) gives the best relation
 between acceptance and $N$ to find minimal value of $\sigma_{95}$
 for a given mass of the excited electron. This is readily apparent
using Figs.~\ref{Kine1} and \ref{Kine2} which shows that the best separation
between signal and background is the $M_{e\gamma}$ variable.


\begin{table}[tb]
 \centering
 \begin{tabular}{|c|c|c|c|c|c|c|}
 \hline
 $M_{e^*}$ & $\sigma^{prod}$ & $\cal B$ ($e^{*} \rightarrow e\gamma$) &  Cut & N$_{back}$ & $\sigma_{\rm 95\%~C.L}^{expected}$ & $\sigma_{\rm 95\%~C.L}^{expected}\cdot\cal B$ \\
 (GeV) & (fb) &    & (GeV) & & (fb) & (fb)\\ \hline
  150 & 712 & 0.400 & 140.0  & 42.0 & 41.3 & 16.5\\
 \hline
  200 & 312  & 0.334  & 185.0 & 18.0 &  30.9 & 10.3\\
 \hline
  250 & 160  & 0.309 & 230.0 & 7.64 & 23.7 & 7.3\\
 \hline
  300 & 87 & 0.297  & 275.0 & 4.36 & 20.0 & 5.9\\
 \hline
\end{tabular}

\caption{Values for $\frac {f}{\Lambda}$ = $\frac {1}{100}$ GeV$^{-1}$ optimized points.}
\label{optimized}
\end{table}


Figure~\ref{sigma95} shows the expected 95\% C.L. cross
section upper limit, $\sigma_{95}$, as a function of the $M_{e\gamma}$ cut for different masses of $e^{*}$ for $\frac{f}{\Lambda}$ = $\frac{1}{100}$ GeV$^{-1}$. Placing our cut at the minimum of the curve gives our optimization, and final expected cross section limit. These results are shown in Fig.~\ref{cross95}. Using the
 Feynman rules for the Lagrangian in Eqn.~\ref{EffLagrangian} leads to a simple relation  between cross-section and $\frac{f}{\Lambda}$ for signal:

\begin{eqnarray}
(\Lambda_{1})^{2}\sigma(\frac{f}{\Lambda_{1}},M_{e^{*}})=(\Lambda_{2})^2\sigma(\frac{f}{\Lambda_{2}},M_{e^{*}}) \\
\frac{f}{\Lambda} = \frac{ \sigma_{95} }{10^{4}\sigma(10^{-2},M_{  e^{*} })}
\end{eqnarray}\\

Taking into account this equation, numerical values of signal
cross-sections (where $\frac{f}{\Lambda} = 10^{-2}$ GeV$^{-1}$)
 and results for $\sigma_{95}$ we can make an  exclusion plot in the 
$\frac{f}{\Lambda}$ vs. $M_{e^{*}}$ plane. This result is shown in 
Fig.~\ref{ExclPlot}.  We compare our results to those of LEP and HERA
in Fig.~\ref{Result} which would give the most stringent limits, to date, for 
masses above 190 GeV.

The prospects for searching for excited electrons
at the upgraded Fermilab Tevatron are excellent. We expect that with a single detector and 2 fb$^{-1}$ of data should significantly extend the mass
reach, especially in the low $\frac{f}{\Lambda}$ region and large masses.
We also point out that similar results are obtained for excited muons
which significantly improve on the current limits which are not producible at HERA.


The authors would like to thank John Hobbs for the use of his limit calculator, Bruce Knuteson for his relation in Eqn.~\ref{Bruce's Relation}, and Bhaskar Dutta for helpful discussions.  We would also like to thank the D$\O$
collaboration for the use of their computers to do the simulation work.
We also would like to thanks the Institute of Nuclear Physics of Moscow State University, the University of Maryland and Texas A\&M University for their support.  Finally we would like to thank the US DOE and Russian Ministry of Industry, Science and Technology for their support during this project. The work of E.B. and A.V. was partly supported by the RFBR-DFG 99-02-04011, RFBR 00-01-00704, CERN-INTAS 99-377, and Universities of Russia
990588 grants.




\clearpage


%\mydrawfig{branch_ratios.eps}{Branching ratio for excited leptons to decay as a function of the mass of the excited lepton $M_{e^*}$.}{BranchingRatios}


%\mydrawfig{cross_section_e_e_ph.eps}{Total cross sections production for single and pair production excited electrons to produce the $ee\gamma$ final state.}{CrossSections}

%\mydrawfig{comp_shw1.eps}{A comparison of the kinematical distributions for excited leptons and background processes in the $ee\gamma$ final state. }{Kine1}

%\mydrawfig{comp_shw2.eps}{A comparison of the kinematical distributions for excited leptons and background processes  in the $ee\gamma$ final state. Note that the $M_{e\gamma}$ variable gives the best separation between signal and background.}{Kine2}

%\mydrawfig{cross_95_cut_mey.eps}{The expected 95\% confidence level cross section upper limit, $\sigma_{95}$, as a function of the $M_{e\gamma}$ cut for different masses of $e^{*}$.}{Sigma95}

%%\mydrawfig{ex_plot_150_300.eps}{Exclusion Plot}{ExclPlot}

%\mydrawfig{result.eps}{Comparison of exclusion plots for HERA, LEP and Tevatron}{Result}


%******************************************************************************


\clearpage
 
  \begin{figure}[htc]
  \vspace*{-.75in}
  %\hspace*{.5in}
  \epsfbox{branch.ps}
  \vspace*{-1in}
  \caption{The branching ratios for excited leptons as a function of the mass of the excited lepton.}
  \label{BranchingRatios}
  \end{figure}
    
  \begin{figure}[ht]
  \vspace*{-.75in}
  %\hspace*{1cm}  
  \epsfbox{cross_section_eey.ps}
  \vspace*{-1.25in}
  \caption{Total cross section for the sum of both single and pair 
  production and decay of excited electrons in the $ee\gamma$ final state.}
  \label{CrossSections}
  \end{figure}
  
  \begin{figure}[ht]
  %\vspace*{-.75in}
  %\hspace*{1cm}  
  \epsfbox{comp_shw1.ps}
  \vspace*{-1.25in}
  \caption{A comparison of the kinematic distributions of the final state 
  electrons and photons for excited electron and standard model background 
  processes.  The figures show the results for single and pair production of 
  ${e^*}e \rightarrow ee\gamma$ and ${e^*}{e^*} \rightarrow ee\gamma$.}
  \label{Kine1}
  \end{figure}

  \begin{figure}[ht]
  \vspace*{-.75in}
  %\hspace*{1cm}  
  \epsfbox{comp_shw2.ps}
  \vspace*{-1.25in}
  \caption{A comparison of some of the variables studied which gave good 
  separation between excited leptons and standard model background processes 
  in the $ee\gamma$ final state. Note that the $M_{e\gamma}$ variable gives 
  the best separation between signal and background for $M_{e^{*}}=200$ 
  GeV at $M_{e\gamma}>185$ GeV.}
  \label{Kine2}
  \end{figure}
 
  \begin{figure}[ht]
  \vspace*{-.75in}
  %\hspace*{1cm}
  \epsfbox{cross_95_cut_mey.ps}
  \vspace*{-1.25in}
  \caption{The expected 95\% C.L. cross section upper limit as a function 
  of the $M_{e\gamma}$ cut for different masses of $e^{*}$.  There the cross 
  section is defined as $\sigma (p\bar{p} \rightarrow ee^{*} \rightarrow 
  ee\gamma)$ + $\sigma(p\bar{p} \rightarrow e^{*}e^{*} \rightarrow 
  ee\gamma)$ and $\mathcal {L}$ = 2 fb$^{-1}$.}
  \label{sigma95}
  \end{figure}
 
  \begin{figure}[ht]
  \vspace*{-1.75in}
  %\hspace*{1cm}
  \epsfbox{cross95_150_300.ps}
  \vspace*{-1.25in}
  \caption{Optimized expected 95\% C.L. for $e^* \rightarrow ee\gamma$ 
  production using the $ee\gamma$ final state and 2 fb$^{-1}$ of data with a 
single detector.}
  \label{cross95}
  \end{figure}

  \begin{figure}[ht]
  \vspace*{-.75in}
  %\hspace*{1cm}
  \epsfbox{ex_plot_150_300.ps}
  \vspace*{-1.25in}
  \caption{Exclusion plot for $\frac{f}{\Lambda}$ as a function of $e^*$ mass
  using 2 fb$^{-1}$ worth of data and a single detector.}
  \label{ExclPlot}
  \end{figure}

  \begin{figure}[ht]
  \vspace*{-1.75in}
%  \hspace*{-1.65in}
  \epsfbox{estar_final.ps}
  \vspace*{-1.25in}
  \caption{For masses above about 190 GeV, we expect that the Tevatron should, with 2 fb$^{-1}$ of data, and one detector, produce the most stringent limits.}
  %Comparing the different results for various experiments.  
  %The expected minimum mass reach of the Tevatron starts where the maximum 
  %mass reach of LEP and HERA end. }
  \label{Result}
  \end{figure}

\clearpage

%******************************************************************************


\begin{thebibliography}{20}

\bibitem{Other} On the other hand an existence of the excited states may be
a results of a presence of extra dimensions.
This may occur if the quark, lepton and boson fields can propagate
in extra dimensions. However in that case one would expect an appearance
of excited states in pairs as discussed in  T.~Appelquist, H.~Cheng and B.~A.~Dobrescu, \Journal{\PRD}{64}{035002}{2001}.
  %%CITATION = ;%%
  %``Bounds on universal extra dimensions,''
  %T. Appelquist, \etal, .

\bibitem{g-2} Current experimental constraints come from the analysis of the $g-2$ properties of muon and electron, and from the estimate of quark and
lepton radii. See A.~Czarnecki, W.~J.~Marciano, , talk given at the Fifth International Symposium on Radiative Corrections (RADCOR 2000), Carmel, CA, September 2000; K.~Lane, 1].
  %A.~Czarnecki and W.~J.~Marciano, .
  %%CITATION = ;%%
  %``The muon anomalous magnetic moment: Standard model theory and beyond,''

\bibitem{Hagiwara} K.~Hagiwara \etal, \Journal{\ZPC}{29}{115}{1985}.

\bibitem{Boudjema2} F.~Boudjema \etal, \Journal{\ZPC}{57}{425}{1993}.

\bibitem{Baur} U.~Baur \etal, \Journal{\PRD}{42}{815}{1990}.

\bibitem{Boos} E.~Boos \etal, \Journal{\NPV}{56}{5}{1993}.

\bibitem{LEP}  L3 Collaboration, M.~Acciarri \etal, \Journal{\PLB}{473}{177}{2000};
  %%CITATION = ;%%
  %``Search for excited leptons at s**(1/2) = 189-GeV,''
ALEPH Collaboration, D.~Decamp \etal,  \Journal{\PLB}{236}{501}{1990};
  %%CITATION = PHLTA,B236,501;%%
  %``Search For Excited Leptons In Z0 Decay,''
ALEPH Collaboration, D.~Buskulic \etal, \Journal{\PLB}{385}{445}{1996}.
  %%CITATION = PHLTA,B385,445;%%
  %``Search for excited leptons at 130-GeV - 140-GeV,''

\bibitem{HERA} 
ZEUS Collaboration, M.~Derrick \etal, \Journal{\PLB}{316}{207}{1993};
ZEUS Collaboration, M.~Derrick \etal, \Journal{\ZPC}{65}{627}{1995};
ZEUS Collaboration, J.~Breitweg \etal, \Journal{\ZPC}{76}{631}{1997};
H1 Collaboration, C.~Adloff \etal, Eur. Phys. J. C17, 567 (2000); 
H1 Collaboration, S.~Aid \etal, \Journal{\NPB}{483}{44}{1997}; 
H1 Collaboration, T.~Ahmed \etal, \Journal{\PLB}{340}{205}{1994}

\bibitem{Tevatron Quark} D$\O$ Collaboration,  B.~Abbott \etal, \Journal{\PRL}{82}{2457}{1999};
  %%CITATION = ;%%
  %``The dijet mass spectrum and a search for quark compositeness in anti-p  p collisions at s**(1/2) = 1.8-TeV,''
D$\O$ Collaboration, B. Abbott \etal, \Journal{\PRL}{80}{666}{1998};
  %%CITATION = ;%%
  %``Measurement of dijet angular distributions and search for quark  compositeness,''
D$\O$ Collaboration, B.~Abbott \etal, \Journal{\PRD}{62}{031101}{2000};
  %%CITATION = ;%%
  %``Limits on quark compositeness from high-energy jets in anti-p p  collisions at 1.8-TeV,''
CDF Collaboration, F.~Abe \etal, \Journal{\PRL}{71}{2542}{1993}.
  %%CITATION = PRLTA,71,2542;%%
  %``Search for quark compositeness, axigluons and heavy particles using the dijet invariant mass spectrum observed in p anti-p collisions,''

\bibitem{CompHEP} E.~E.~Boos \etal,   %preprint INP MSU 94-36/358 and SNUTP-94-116; 
  %E.~E.~Boos, M.~N.~Dubinin, V.~A.~Ilin, A.~E.~Pukhov and V.~I.~Savrin,
  %``CompHEP: Specialized package for automatic calculations of elementary particle decays and collisions,''
  %.
  %%CITATION = ;%%
P.~Baikov \etal, Workshop on High Energy Physics and Quantum Field Theory, ed. by B.~Levtchenko and V.~Savrin, Moscow, 1996, 101;
  %B. B. Levchenko and V. I. Savrin, Moscow, Russia: State U. (1996) 401;
  %``High-energy physics and quantum field theory. Proceedings, 10th International Workshop, Zvenigorod, Russia, September 20-26, 1995,''
A.~Pukhov \etal, . 
  %Preprint INP MSU 98-41/542.
  %%CITATION = ;%%
  % ``CompHEP: A package for evaluation of Feynman diagrams and integration  over multi-particle phase space. User's manual for version 33,''

\bibitem{LanHEP} A.~Semenov, \Journal{\NIMA}{393}{293}{1997}; .
  %%CITATION = INP MSU Preprint 98-2/503;%%
  %``LanHEP: A package for automatic generation of Feynman rules in gauge  models,''

\bibitem{interface} A.~S.~Belyaev et al.,  To appear in the proceedings of the  Seventh International Workshop on Advanced Computing and Analysis Technics in Physics Research (ACAT2000, Fermilab, October 16-20, 2000).

\bibitem{Pythia} T. Sj\"ostrand, Comp. Phys. Comm. {\bf 82}, 74 (1994); 
  %``QCD radiation off heavy particles,''
  %%CITATION = ;%%  
  S.~Mrenna, Computer Phys. Commun. {\bf 101}, 232 (1997); 
  T.~Sj\"ostrand, P.~Eden, C.~Friberg, L.~Lonnblad, G.~Miu, S.~Mrenna and E.~Norrbin, Comp. Phys. Comm.  {\bf 135}, 238 (2001).
  %%CITATION = ;%
  %Title: QCD Radiation off Heavy Particles
  %``High-energy-physics event generation with PYTHIA 6.1,''

\bibitem{CTEQ4L} H.~L.~Lai \etal. \Journal{\PRD}{55}{1280}{1997}.

\bibitem{K-signal} H.~Baer \etal, \Journal{\PRD}{57}{5871}{1998}.

\bibitem{DZero Diboson Studies} D$\O$ Collaboration, F.~Abachi \etal, \Journal{\PRL}{75}{1023}{1995}.
  %``Search for W boson pair production in p anti-p collisions at s**(1/2) = 1.8-TeV,''
  %%CITATION = ;%%

\bibitem{K-background} J.~Smith \etal, \Journal{\ZPC}{44}{267}{1989}.

\bibitem{SHW} The Supersymmetry/Higgs Workshop, Version 2.3, Fermilab, Batavia IL, November 19, 1998. http://fnth37.fnal.gov/susy.html

\bibitem{CDF} CDF Collaboration, F.~Abe \etal, \Journal{\NIMA}{271}{387}{1988}.

\bibitem{DZero} D$\O$ Collaboration, F.~Abachi \etal, \Journal{\NIMA}{338}{185}{1994}.

\bibitem{Frequentest} G.~Zech, \Journal{\NIMA}{227}{608}{1989}; 
T. Huber \etal, \Journal{\PRD}{41}{2709}{1990}.

\end{thebibliography}

%***************************************************************************** 


\end{document}



