\documentclass[12pt]{article}
\usepackage{cite,epsfig}

\input paperdef

\hyphenation{ele-gant ligh-test re-nor-ma-li-za-ti-on pa-ra-me-ter
  pa-ra-me-ters ap-pro-xi-ma-tion} 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


\oddsidemargin -0.5cm
\evensidemargin -0.1cm
\marginparwidth 68pt
\marginparsep 10pt
\topmargin 0cm
\headheight 0pt
\headsep 0pt
\footskip 25pt
\textheight 22cm
\textwidth 16.5cm
\columnsep 10pt
\columnseprule 0pt



\begin{document}

\thispagestyle{empty}
\setcounter{page}{0}
\def\thefootnote{\fnsymbol{footnote}}

\begin{flushright}
CERN-TH/2000-55\\
DESY 00-020\\
KA-TP-3-2000\\
 \\
%\date{\today}
\end{flushright}

\vspace{1cm}

\begin{center}

{\large\sc {\bf \fhf: a program for a fast calculation}}

\vspace*{0.4cm} 

{\large\sc {\bf of masses and mixing angles}}

\vspace*{0.4cm}

{\large\sc {\bf in the Higgs Sector of the MSSM}}
\footnote{Extended version of the contribution to the 
          {\em Report of the HIGGS working group for the
                 Workshop ``Physics at TeV Colliders''}, 
                 Les Houches, France, June 1999.}
 
\vspace{1cm}

{\sc 
S.~Heinemeyer$^{1}$%
\footnote{email: Sven.Heinemeyer@desy.de}%
, W.~Hollik$^{2}$%
\footnote{email: Wolfgang.Hollik@physik.uni-karlsruhe.de}%
~and G.~Weiglein$^{3}$%
\footnote{email: Georg.Weiglein@cern.ch}
}

\vspace*{1cm}

{\sl
$^1$ 
DESY Theorie, Notkestr. 85, 22603 Hamburg, Germany

\vspace*{0.4cm}

$^2$ Institut f\"ur Theoretische Physik, Universit\"at Karlsruhe, \\
D--76128 Karlsruhe, Germany

\vspace*{0.4cm}

$^3$ CERN, TH Division, CH-1211 Geneva 23, Switzerland
}


\end{center}


\vspace*{1cm}

\begin{abstract}
\fhf\ is a Fortran code for the calculation of the masses and the
mixing angle of the
neutral $\cp$-even Higgs bosons in the MSSM up to \twol\ order. It is based
on a compact analytical approximation formula of the complete
diagrammatic \onel\ and the dominant \twol\ contributions.
At the \onel\ level a
leading logarithmic result is used, taking into account all sectors of
the MSSM. At the \twol\ level at $\oaas$ the leading logarithmic and
non-logarithmic contributions are taken into account. 
The approximation formula is valid for 
arbitrary choices of the parameters in the Higgs sector of the model.
Comparing its quality to the full diagrammatic
result, we find agreement better than $2 \gev$ for
most parts of the MSSM parameter space.
\end{abstract}



\def\thefootnote{\arabic{footnote}}
\setcounter{page}{0}
\setcounter{footnote}{0}

\newpage



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


\section{Introduction}

One of the most striking phenomenological implications of Supersymmetry
(SUSY) is the prediction of a relatively light Higgs boson, which is 
common to all supersymmetric models whose couplings remain in the
perturbative regime up to a very high energy
scale~\cite{susylighthiggs}. The search for the lightest Higgs boson thus 
allows a crucial test of SUSY and especially the Minimal
Supersymmetric Standard Model (MSSM)~\cite{mssm}. Therefore it is one
of the main 
goals at the present and the next generation of colliders. 
A precise knowledge of the dependence of the masses and mixing angles
of the Higgs sector of the MSSM
on the relevant SUSY parameters is necessary for a
detailed analysis of SUSY phenomenology at LEP2, the upgraded Tevatron,
and also for the LHC and a future linear $e^+e^-$ collider, where
high-precision measurements in the Higgs sector might become possible.

The mass of the lightest $\cp$-even MSSM Higgs boson, $\mh$, is
bounded from above by about $\mh \lsim 135 \gev$, including radiative
corrections up to the \twol\ level~%
\cite{mhiggsletter,mhiggslong,mhiggsRG1a,mhiggsRG1b,mhiggsRG2}.
The mixing angle $\aeff$ that diagonalizes the neutral $\cp$-even Higgs
sector receives the same kind of corrections. By incorporating
$\aeff$ into the Higgs decay widths, the leading electroweak corrections
to the decay of the neutral $\cp$-even Higgs bosons are taken into account.

In this note we present the Fortran code \fhf. Using low
energy MSSM parameters as input, it evaluates the masses of the neutral
$\cp$-even Higgs bosons, $\mh$ and $\mH$, as well as the corresponding
mixing angle, $\aeff$, at the \twol\ level. In addition the mass of
the charged Higgs boson, $\mhp$, is evaluated at the \onel\
level. The \rp, leading to constraints in the scalar fermion
sector of the MSSM, is evaluated up to $\oaas$, taking into account the
gluon exchange contribution at the \twol\ level~\cite{drhosuqcd}.
\fhf\ is based on a compact analytical approximation formula,
containing at the \twol\ level the dominant corrections in
$\oaas$ obtained in the Feynman-diagrammatic (FD)
approach~\cite{mhiggsletter,mhiggslong} and subdominant corrections of
$\ogmzmts$ obtained with renormalization group (RG)
methods~\cite{mhiggsRG1a,mhiggsRG1b,mhiggsRG2}.
In comparison with the FD result, consisting of the complete
\onel\ and the \twol\ contributions as given in \citere{mhiggslong}
(incorporated into the Fortran code \fh~\cite{feynhiggs}), the
approximation formula is much shorter. Thus \fhf\ is about $3 \times 10^4$
times faster than \fh. Agreement between these two codes of
better than $2 \gev$ is found for most parts of the MSSM parameter
space. 

The following sections are organized as follows: In Sect.~2 we 
shortly summarize the analytical approximation formula and perform a
comparison with the full FD result.
A description of how to use \fhf\ is given in
Sect.~3. The conclusions can be found in Sect.~4.


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\section{The analytical approximation}

In order to fix our notations, 
we first list the conventions for the
MSSM scalar top sector:
the mass matrix in the basis of the current eigenstates $\StopL$ and
$\StopR$ is given by (neglecting contributions $\sim \MZ^2$):
\BE
\label{stopmassmatrixwithdt}
{\cal M}^2_{\Stop} =
  \ML \MstL^2 + \mt^2  &
      \mt \Xt \\
      \mt \Xt &
      \MstR^2 + \mt^2  
  \MR,
\EE
where 
\BE
\mt \Xt = \mt (A_t - \mu \CTb)~.
\label{eq:mtlr}
\EE
Furtheron the following simplification is used:
\BE
\MstL = \MstR := \msq .
\EE
In this simplified case we define%
%
\footnote{
The more general case with $\MstL \neq \MstR$ can be found in
\refeq{msfkt}. 
}
%
\BE
\ms^2 := \msq^2 + \mt^2 .
\EE
For further details see \citere{mhiggslle}.


At the tree level the mass matrix of the neutral $\cp$-even Higgs bosons
in the $\phi_1,\phi_2$ basis can be expressed as follows:
\BEA
M_{\rm Higgs}^{2, {\rm tree}} &=& \ML \mpe^2 & \mpez^2 \\ 
                           \mpez^2 & \mpz^2 \MR \non\\
&=& \ML \MA^2 \SQb + \MZ^2 \CQb & -(\MA^2 + \MZ^2) \Sb \Cb \\
    -(\MA^2 + \MZ^2) \Sb \Cb & \MA^2 \CQb + \MZ^2 \SQb \MR.
\label{higgsmassmatrixtree}
\EEA
$\MA$ denotes the mass of the $\cp$-odd Higgs boson, $\MZ$ is the mass
of the $Z$ boson, and $\tb = v_2/v_1$ is the ratio of the two vacuum
expectation values of the two Higgs doublets in the MSSM, see \citere{hhg}.

At higher orders, the Higgs mass matrix (\ref{higgsmassmatrixtree}) is
supplemented by the renormalized self-energies 
$\hSi_s(q^2)\,,\;s = \Pe,\Pz,\PePz$. Here we use the approximation of
neglected external 
momentum: $\hSi_s \equiv \hSi_s(0)$. The masses $\mh$ and
$\mH$ are then 
obtained by diagonalizing the higher order corrected mass matrix with
the mixing angle $\aeff$:
\BEA
&&
\label{higgsmassmatrix}
\ML \mpe^2 - \hSi_{\Pe} & \mpez^2 - \hSi_{\PePz} \\ 
    \mpez^2 - \hSi_{\PePz} & \mpz^2 - \hSi_{\Pz} \MR 
   \stackrel{\aeff}{\longrightarrow}
   \ML \mH^2 & 0 \\ 0 &  \mh^2 \MR~, \\
\label{alphaeff}
&& \aeff = {\rm arctan}\KKL 
  \frac{-(\MA^2 + \MZ^2) \Sb \Cb - \hSi_{\PePz}}
       {\MZ^2 \CQb + \MA^2 \SQb - \hSi_{\Pe} - \mh^2} \KKR~,~~
  -\frac{\pi}{2} < \aeff < \frac{\pi}{2}~.
\EEA


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%\subsection{Calculation of the renormalized self-energies}

Here we only give a very brief outline of the calculation of the
renormalized Higgs boson self-energies. A detailed description can be
found in \citere{mhiggslle}.

The mass of the lightest Higgs boson receives contributions from all
sectors of the MSSM, but not all are numerically of equal relevance.
The dominant corrections arise from the $t-\Stop$ sector of the MSSM: 
The results for the renormalized self-energies $\hSi_s$ have been
derived analytically in the FD approach in
\citeres{mhiggsletter,mhiggslong}. These results, however, are rather
lengthy. In order to derive a compact analytical expression, %the following 
several approximations have been made as described in \citere{mhiggslle}.
The main step of the approximation consists of a Taylor
expansion in 
\BE
\dst = \frac{|\mt\,\Xt|}{\ms^2}
\label{deltastop}
\EE
of the $\hSi_s(0)$.
For the \onel\ correction the expansion has been performed up to 
${\cal O}(\dst^8)$; 
all three renormalized \onel\ Higgs-boson self-energies give a contribution.
%In the \onel\ self-energies we have kept terms up to
%${\cal O}(\MZ^4/\mt^4)$ while terms of ${\cal O}(\MZ^2/\ms^2)$ have
%been neglected. We have checked that the numerical effect of the
%latter terms is insignificant.
Concerning the \twol\ self-energies,
the expansion has been carried out
up to ${\cal O}(\dst^4)$; 
in the approximation considered here only $\hSiz_{\Pz}$ gives a
non-zero contribution. 
%Further details can be found in \citere{mhiggslle}.

Leading contributions beyond $\oaas$ have been taken into account by
incorporating the leading \twol\ Yukawa 
correction of ${\cal O}(\gf^2\mt^6)$
\cite{mhiggsRG1a,mhiggsRG1b,mhiggsRG2}. Furthermore the result has
been expressed in terms of the \msbar\ top-quark mass
\BE
\mtms = \mtms(\mt) \approx \frac{\mt}{1 + \frac{4}{3\,\pi} \als(\mt)}~
\label{mtrun}
\EE
instead of the pole mass $\mt$.
This leads to an additional contribution in 
${\cal O}(\al\als^2)$.

The analytical expressions arising from the $t-\Stop$ sector are
given at the \onel\ level as follows:
%
\BEA
\label{p1seol}
\hSie_{\Pe}(0) &=&
    \frac{\gf\wz}{\pi^2} \MZ^4 \Lambda \CQb
    \log\lmtmsms~, \\
\label{p1p2seol}
\hSie_{\PePz}(0) &=&
    -\frac{\gf\wz}{\pi^2} \MZ^2 \CTb 
    \KKL -\frac{3}{8} \mtms^2 + \MZ^2 \Lambda \SQb \KKR \log\lmtmsms~, \\
\label{p2seol}
\hSie_{\Pz}(0) &=&
    \frac{\gf\wz}{\pi^2} \frac{\mtms^4}{8 \SQb} 
      \Bigg\{ -2 \frac{\MZ^2}{\mtms^2} 
              + \frac{11}{10}\frac{\MZ^4}{\mtms^4} \non \\
 && {} + \KKL 12 - 6 \frac{\MZ^2}{\mtms^2} \SQb +
                 8 \frac{\MZ^4}{\mtms^4} \Lambda \SVb \KKR \log\lmtmsms \non \\
 && {} + \frac{\Xt^2}{\ms^2}\, 
      \KKL -12 + 4 \frac{\MZ^2}{\mtms^2}
          + 6 \frac{\mtms^2}{\ms^2} \KKR 
    + \frac{\Xt^4}{\ms^4}\,
      \KKL 1 - 4\,\frac{\mtms^2}{\ms^2} 
           +3 \frac{\mtms^4}{\ms^4} \KKR \non \\
 && {} + \frac{\Xt^6}{\ms^6}\,
      \KKL  \frac{3}{5}\frac{\mtms^2}{\ms^2}
           -\frac{12}{5}\frac{\mtms^4}{\ms^4}
           +2\frac{\mtms^6}{\ms^6} \KKR \non \\
 && {} + \frac{\Xt^8}{\ms^8}\,
      \Bigg[  \frac{3}{7}\frac{\mtms^4}{\ms^4}
             -\frac{12}{7}\frac{\mtms^6}{\ms^6}
             +\frac{3}{2}\frac{\mtms^8}{\ms^8} \Bigg]
    \Bigg\}~, %\\[1em]
\EEA
with
\BE
\Lambda = \left(\frac{1}{8} - \frac{1}{3} \sw^2 + \frac{4}{9} 
                 \sw^4 \right), \qquad \sw^2 = 1 - \frac{\MW^2}{\MZ^2} . 
\label{eq:lambda}
\EE
The \twol\ contributions in $\oaas$ read:
\BEA
\hSiz_{\Pe}(0) &=& 0~, \non \\
\hSiz_{\PePz}(0) &=& 0~, \non \\
\label{p2setl}
\hSiz_{\Pz}(0) &=&
    \frac{\gf\wz}{\pi^2} \frac{\als}{\pi} \frac{\mtms^4}{\SQb}
      \Biggl[ 4 + 3 \log^2\lmtmsms + 2 \log\lmtmsms 
             -6 \frac{\Xt}{\ms} \non \\
 && {}  %\hspace{7em}
     - \frac{\Xt^2}{\ms^2} \KKKL 3 \log\lmtmsms +8 \KKKR
     +\frac{17}{12} \frac{\Xt^4}{\ms^4} \Biggr]~.
\EEA
The \twol\ Yukawa correction in this approximation reads:
\BE
\label{mh2yuk}
\hSiz_{\Pz}(0) = -\frac{9}{16\pi^4} \frac{\gf^2 \mtms^6}{\SQb}
               \KKL \log^2\lmtmsms 
                    - 2 \frac{\Xt^2}{\ms^2} \log\lmtmsms
                    + \ed{6} \frac{\Xt^4}{\ms^2} \log\lmtmsms \KKR~.
\EE
$\ms$ has to be chosen according to
\BE
\label{msfkt}
\ms = \KKKL \begin{array}{l@{\quad:\quad}l}
            \sqrt{\msq^2 + \mtms^2} & \MstL = \MstR = \msq \\
            \KKL \MstL^2 \MstR^2 + \mtms^2 (\MstL^2 + \MstR^2) +
                     \mtms^4 \KKR^\frac{1}{4} &
                                    \MstL \neq \MstR
            \end{array} \right. 
\EE

For the \onel\ corrections from the other sectors of the MSSM 
the logarithmic approximation given in \citere{mhiggs1lrest} has been
used for \fhf, see also \citere{mhiggslle}. 
In order to obtain the radiatively corrected Higgs boson masses and
the mixing angle, the renormalized Higgs boson self-energies have to be
inserted into \refeq{higgsmassmatrix}, 
and the corresponding diagonalization has to be performed.


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%\subsection{Calculation of the \rp}
\bigskip

We have also implemented the result of the MSSM contributions to
$\De\rho$~\cite{drhosuqcd,precobssusyqcd}. Here the corrections 
arising from $\Stop/\Sbot$-loops up to $\oaas$
have been taken into account. The result is valid for arbitrary
parameters in  
the $\Stop$- and $\Sbot$-sector, also taking into account the mixing
in the $\Sbot$-sector which can have a non-negligible
effect in the large $\tb$ scenario~\cite{precobssusyqcd}.  

The \twol\ result can be separated into the pure gluon-exchange contribution,
which can be expressed by a very compact formula, allowing a very
fast evaluation, and the pure
gluino-exchange contribution, which is given by a rather lengthy 
expression. The latter correction goes to zero with increasing gluino mass
and can thus be discarded for a heavy gluino.
Concerning the implementation of the \rp\ into \fhf, we have neglected the
gluino exchange contribution.
The \rp\ can be used as an additional constraint (besides
the experimental bounds) on the squark masses. 
A value of $\De\rho$ outside the experimentally preferred region of 
$\dr^{\SU} \approx 10^{-3}$~\cite{delrhoexp} indicates experimentally
disfavored $\Stop$- and $\Sbot$-masses.



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%\section{Discussion of the compact expression for $\mh$ and $\aeff$}

\bigskip
For illustration of the quality of the compact approximation,
we compare $\mh$ and $\aeff$ with the results from the complete
%In this section we analyze the quality of our
%compact approximation formula for $\mh$ and the mixing angle $\aeff$
%with respect to the full 
Feynman-diagrammatic calculation.
The FD calculation contains the full diagrammatic \onel\
contribution~\cite{mhiggs1lfull}, the complete leading \twol\
corrections in $\oaas$~\cite{mhiggsletter,mhiggslong},
and the two contributions beyond $\oaas$ (see \refeq{mtrun} 
and \refeq{mh2yuk}) without approximation. 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[ht!]
\begin{center}
 \mbox{
\psfig{figure=mh_mq_MA500_mgl500_tb1p6.eps,width=7cm,height=5.4cm}}
\hspace{1em}
\mbox{
\psfig{figure=mh_mq_MA500_mgl500_tb40.eps,width=7cm,height=5.4cm}}
\end{center}
\vspace{-1.5em}
\caption[]{
$\mh$ as a function of $\msq$, calculated from the full formula and
from the approximation formula for $\MA = 500 \gev, 
\mgl = 500 \gev$ and $\tb = 1.6$ or $40$.
}
\label{fig:mh_mq}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[ht!]
%\vspace{2cm}
\begin{center}
\mbox{
\psfig{figure=mh_MA_mq1000_mgl500_tb1p6.eps,width=7cm,height=5.4cm}}
\hspace{1em}
\mbox{
\psfig{figure=mh_MA_mq1000_mgl500_tb40.eps,width=7cm,height=5.4cm}}
\end{center}
\vspace{-1.5em}
\caption[]{
$\mh$ as a function of $\MA$, calculated from the full formula and
from the approximation formula for $\msq = 1000 \gev, 
\mgl = 500 \gev$ and $\tb = 1.6$ or $40$.
}
\label{fig:mh_MA}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[ht!]
%\vspace{2cm}
\begin{center}
\mbox{
\psfig{figure=ghRC2l07b.bw.eps,width=7cm,height=5.4cm}}
\hspace{1em}
\mbox{
\psfig{figure=ghRC2l08b.bw.eps,width=7cm,height=5.4cm}}
\end{center}
\vspace{-1.5em}
\caption[]{
The relative difference $(\sin\aeff - \sin\aeffapprox)/\sin\aeff$  
is shown as a function
of $\MA$ for three values of $\tb$ in the no mixing and the maximal
mixing scenario. The other parameters are $\mu = -100 \gev$, 
$M = \msq$, $\mgl = 500 \gev$ and $A_b = A_t$.
}
\label{fig:saeff_MA}
\vspace{-2em}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


For the Standard Model parameters we use 
$\MZ = 91.187 \gev,~\MW = 80.39 \gev,~
G_F = 1.16639\times 10^{-5} \gev^{-2},~\als(\mt) = 0.1095$, and
$\mt = 175 \gev$. In the numerical evaluation we
have furthermore chosen the trilinear couplings in the scalar top and
bottom sector to be $A_b = A_t$. This fixes together with the
choice of $\mu$ (the Higgs mixing parameter) the mixing in the scalar
bottom sector. 
The parameter $M$ appearing in the plots is the $SU(2)$ gaugino mass
parameter, it enters as an independent parameter in the full result
only, see \citere{mhiggslong}. 

In \citeres{mhiggslong,mhiggslle} it has been shown that the lightest
Higgs boson mass $\mh$ as a function of $\Xt$ reaches a maximum at
around $|\Xt/\msq| \approx 2$. This case we refer to as 'maximal
mixing'. A minimum of $\mh$ is reached for $\Xt \approx 0$, which we
refer to as 'no mixing'.

\smallskip
\reffi{fig:mh_mq} displays the dependence of $\mh$ on $\msq$ for the cases
of no mixing and maximal mixing, and we have set $\MA = 500 \gev$.
For $\tb$ we have chosen two typical values: $\tb = 1.6$ as a low and
$\tb = 40$ as a typical high value.
Very good agreement is found  in the
no-mixing scenario as well as in the maximal-mixing scenario,
the deviations lie below $2 \gev$.

The dependence on $\MA$ is shown in \reffi{fig:mh_MA}. 
The quality of the approximation is typically better than $1 \gev$ for
the no-mixing case and better than $2 \gev$ for the maximal-mixing
case. Only for very small (and experimentally already excluded) values
of $\MA$ a deviation of $5 \gev$ occurs. 
The peaks in the plot for $\tb = 1.6$ in the full result are due to
the threshold  $\MA = 2\,\mt$ in the
\onel\ contribution, originating from the top-loop diagram in
the $A$~self-energy. This peak does not occur in the approximation
formula (where the momentum dependence of the $A$ self-energy
has been neglected) and can thus lead to a larger deviation around the
threshold. 

In \reffi{fig:saeff_MA} we display the quality of the short analytical
approximation formula for the effective mixing angle. The relative
difference 
between $\aeff$ obtained from the full calculation (i.e.\ calculating
the renormalized Higgs self-energies $\hSi_s(0)$ without any
approximation) and the angle obtained with the help of the
approximated Higgs self-energies, denoted as $\aeffapprox$ is shown as
a function of $\MA$. 
We use three values of $\tb$, $\tb = 3, 20, 40$.
Sizable deviations
occur only in the region $100 \gev \le \MA \le 150 \gev$, especially
for large $\tb$. In this region of parameter space the values of $\mh$
and $\mH$ are very close to each other. This results in a high
sensitivity to small deviations in the Higgs boson self-energies
entering the Higgs-boson mass
matrix~(\ref{higgsmassmatrix}),~(\ref{alphaeff}). Otherwise the
relative difference stays below~3\%.



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\section{The Fortran program \fhf}

%\subsection{The main structure}

The complete program \fhf\ consists of about 1300 lines Fortran code.
The executable file fills about 65 KB disk space.
The calculation for one set of parameters, including the $\De\rho$
constraint, takes about $2 \times 10^{-5}$ seconds on 
a Sigma station (Alpha processor, 600 MHz processing speed, 512 MB RAM).
The program can be obtained from the \fh\ home page:\\
{\tt http://www-itp.physik.uni-karlsruhe.de/feynhiggs}~.\\
Here the code is available,
together with a short instruction, information about bug fixes etc.

\fhf\ consists of a front-end ({\tt program FeynHiggsFast}) and the
main part where the evaluation is performed (starting with 
{\tt subroutine feynhiggsfastsub}). The front-end can be manipulated
by the user at will, whereas the main part should not be changed. In
this way \fhf\ can be accommodated as a subroutine to existing
programs%
\footnote{
This has been carried out, for example, for the program
HDECAY~\cite{hdecay,spirix}, into which \fhf\ has been incorporated as
a subroutine.
}
%
, thus providing an extremely fast and for many purposes sufficiently
accurate evaluation for the masses and
mixing angles in the MSSM Higgs sector. 

The input of \fhf\ are the low energy SUSY parameters, listed in
\refta{tab:mssmparameters}. Concerning the $\Stop$-sector, the user
has the option to enter either the physical parameters, i.e.\ the
masses and the mixing angle ($\mste$, $\mstz$ and $\sintt$) or the
unphysical parameters $\MstL$, $\MstR$ and $\Xt$. From these input
parameters \fhf\ calculates the masses and the mixing angle of the
neutral $\cp$-even Higgs sector, as well as the mass of the charged
Higgs boson and the \rp.

\begin{table}[ht!]
\renewcommand{\arraystretch}{1.5}
\begin{center}
\begin{tabular}{|c||c||c|} \hline
{\bf input for \fhf} & {\bf MSSM expr.} & 
{\bf internal expr. in \fhf}
\\ \hline \hline
{\tt tan(beta)}          & $\Tb$                     & {\tt ttb} \\
{\tt Msusy\_top\_L  }    & $\MstL$                   & {\tt msusytl} \\
{\tt Msusy\_top\_R}      & $\MstR$                   & {\tt msusytr} \\
{\tt MtLR }              & $\Xt \equiv \Mtlr$        & {\tt mtlr} \\
{\tt MSt2}               & $\mstz$                   & {\tt mst2} \\
{\tt delmst}             & $\delmst = \mstz - \mste$ & {\tt delmst} \\
{\tt sin(theta\_stop)}   & $\sintt$                  & {\tt stt} \\
{\tt MT}                 & $\mt$                     & {\tt mmt} \\
{\tt Mue}                & $\mu$                     & {\tt mmue} \\
{\tt MA}                 & $\MA$                     & {\tt mma} \\
\hline
\end{tabular}
\renewcommand{\arraystretch}{1}
\caption[]{The meaning of the different MSSM variables which have to
be entered into \fhf.}
\label{tab:mssmparameters}
\end{center}
\end{table}


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\section{Conclusions}

\fhf\ is a Fortran code for the calculation of masses and mixing
angles in the Higgs sector of the MSSM. In addition it evaluates the
SUSY corrections to the \rp, arising from the scalar top and bottom sector.
Concerning the evaluation in the neutral $\cp$-even Higgs sector, 
\fhf\ is based on a simple analytic formula, derived from a more
complete result obtained in 
the Feynman-diagrammatic approach. 
Beyond the \onel\ level, \fhf\ contains the dominant corrections in
$\oaas$ and further subdominant contributions.
The accuracy of the approximation compared to the full result is better
than $2 \gev$ for most parts of the MSSM parameter space.

The program is available via the WWW page\\
{\tt http://www-itp.physik.uni-karlsruhe.de/feynhiggs}.

The code consists of a front-end and a subroutine. The front-end can
be manipulated at the user's will. By accommodating the front-end,
\fhf\ can serve as a subroutine to existing programs, thus providing
an extremely fast and reasonably accurate evaluation of the masses and 
mixing angles in the MSSM Higgs sector. These 
can then be used as inputs for further computations.



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%



\bigskip
\subsection*{Acknowledgements}
S.H. thanks the organizers of the Les Houches workshop for the
inspiring and relaxed atmosphere.
%In addition we would like to express our deep gratitude to
%W.~Chopin, W.~Wyborowa and W.~Belv\'ed\`ere for constant encouragement.


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%\newpage

\begin{thebibliography}{00}  

\bibitem{susylighthiggs} G.~Kane, C.~Kolda and J.~Wells,
                         {\em Phys. Rev. Lett.} {\bf 70} (1993) 2686,
                         ;\\
                         J.~Espinosa and M.~Quir\'os,
                         {\em Phys. Lett.} {\bf B 302} (1993) 51,
                         ;
                         {\em Phys. Rev. Lett.} {\bf 81} (1998) 516,
                         .

\bibitem{mssm} H.~Haber and G.~Kane, 
               {\em Phys. Rep.} {\bf 117} (1985) 75. \\
               H.P.~Nilles, 
               {\em Phys. Rep.} {\bf 110} (1984) 1.


\bibitem{mhiggsletter} S.~Heinemeyer, W.~Hollik and G.~Weiglein,
                       {\em Phys. Rev.} {\bf D 58} (1998) 091701,
                       ;\\
                       {\em Phys. Lett.} {\bf B 440} (1998) 296,
                       .

\bibitem{mhiggslong} S.~Heinemeyer, W.~Hollik and G.~Weiglein,
                     {Eur. Phys. Jour.} {\bf C 9} (1999) 343,
                     .

\bibitem{mhiggsRG1a}   J.~Casas, J.~Espinosa, M.~Quir\'os and A.~Riotto,
                      {\em Nucl. Phys.} {\bf B 436} (1995) 3,
                      E: {\em ibid.} {\bf B 439} (1995) 466,
                      .

\bibitem{mhiggsRG1b}  M.~Carena, J.~Espinosa, M.~Quir\'os and C.~Wagner,
                      {\em Phys. Lett.} {\bf B 355} (1995) 209,
                      ;\\
                      M.~Carena, M.~Quir\'os and C.~Wagner,
                      {\em Nucl. Phys.} {\bf B 461} (1996) 407,
                      .

\bibitem{mhiggsRG2} H.~Haber, R.~Hempfling and A.~Hoang,
                    {\em Z. Phys.} {\bf C 75} (1997) 539,
                    .

%\bibitem{mhiggsEffPot} R.~Hempfling and A.~Hoang,
%                       {\em Phys. Lett.} {\bf B 331} (1994) 99,
%                       .

\bibitem{drhosuqcd} A.~Djouadi, P.~Gambino, S.~Heinemeyer, W.~Hollik,
                    C.~J\"unger and G.~Weiglein,
                    {\em Phys. Rev. Lett.} {\bf 78} (1997) 3626,
                    ;
                    {\em Phys. Rev.} {\bf D 57} (1998) 4179,
                    .

\bibitem{feynhiggs} S.~Heinemeyer, W.~Hollik and G.~Weiglein,
                    {\em Comp. Phys. Comm.} {\bf 124} (2000) 76,
                    .

\bibitem{mhiggslle}  S.~Heinemeyer, W.~Hollik and G.~Weiglein,
                     {\em Phys. Lett.} {\bf B 455} (1999) 179,
                     .

\bibitem{hhg} J.~Gunion, H.~Haber, G.~Kane and S.~Dawson, 
              {\em The Higgs Hunter's Guide}, 
              Addison-Wesley, 1990.

\bibitem{mhiggs1lrest} H.~Haber and R.~Hempfling,
                       {\em Phys. Rev.} {\bf D 48} (1993) 4280,
                       .

\bibitem{precobssusyqcd} S.~Heinemeyer, W.~Hollik and G.~Weiglein,
                         {in preparation.}                   

\bibitem{delrhoexp} G.~Altarelli, .

\bibitem{mhiggs1lfull}  A.~Dabelstein, 
                        {\em Z. Phys.} {\bf C 67} (1995) 495,
                        ;
                        {\em Nucl. Phys.} {\bf B 456} (1995) 25,
                        .

\bibitem{hdecay} A.~Djouadi, J.~Kalinowski and M.~Spira,
                 {\em Comp. Phys. Comm.} {\bf 108} (1998) 56,
                 .

\bibitem{spirix} A.~Djouadi et al. (conv.),
                 {\em Report of the HIGGS working group for the
                 Workshop ``Physics at TeV Colliders''}, 
                 Les Houches, France, June 1999.\\
                 B.~Allanach et al.,
                 Proceedings of {\em UK Phenomenology Workshop on
                 Collider Physics}, Durham, UK, September 1999,
                 .


\end{thebibliography}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%



\end{document}


\catcode`@=11
\def\citer{\@ifnextchar [{\@tempswatrue\@citexr}{\@tempswafalse\@citexr[]}}
 
\def\@citexr[#1]#2{\if@filesw\immediate\write\@auxout{\string\citation{#2}}\fi
  \def\@citea{}\@cite{\@for\@citeb:=#2\do
    {\@citea\def\@citea{--\penalty\@m}\@ifundefined
       {b@\@citeb}{{\bf ?}\@warning
       {Citation `\@citeb' on page \thepage \space undefined}}%
\hbox{\csname b@\@citeb\endcsname}}}{#1}}
\catcode`@=12

%\renewcommand{\topfraction}{1.0}
%\renewcommand{\bottomfraction}{0.3}
\renewcommand{\textfraction}{0}
\renewcommand{\floatpagefraction}{1.0}

% new commands for cross referencing
%\def\refeq#1{\mbox{(\ref{#1})}}
\def\refeq#1{\mbox{eq.~(\ref{#1})}}
\def\refeqs#1{\mbox{eqs.~(\ref{#1})}}
\def\reffi#1{\mbox{Fig.~\ref{#1}}}
\def\reffis#1{\mbox{Figs.~\ref{#1}}}
\def\refta#1{\mbox{Tab.~\ref{#1}}}
%\def\refta#1{\mbox{Table~\ref{#1}}}
\def\reftas#1{\mbox{Tabs.~\ref{#1}}}
%\def\reftas#1{\mbox{Tables~\ref{#1}}}
\def\refse#1{\mbox{Sect.~\ref{#1}}}
\def\refses#1{\mbox{Sects.~\ref{#1}}}
%\def\refse#1{\mbox{Section~\ref{#1}}}
\def\refapp#1{\mbox{App.~\ref{#1}}}
%\def\refapp#1{\mbox{Appendix~\ref{#1}}}
\def\citere#1{\mbox{Ref.~\cite{#1}}}
\def\citeres#1{\mbox{Refs.~\cite{#1}}}

\newcommand{\mst}{m_{\tilde{t}}}
\newcommand{\delmst}{\Delta\mst}
\newcommand{\mstr}{m_{\tilde{t}_R}}
\newcommand{\mstl}{m_{\tilde{t}_L}}
\newcommand{\mste}{m_{\tilde{t}_1}}
\newcommand{\mstz}{m_{\tilde{t}_2}}
\newcommand{\msti}{m_{\tilde{t}_i}}
\newcommand{\mstip}{m_{\tilde{t}_{i'}}}
\newcommand{\mstj}{m_{\tilde{t}_j}}
\newcommand{\msb}{m_{\tilde{b}}}
\newcommand{\msbr}{m_{\tilde{b}_R}}
\newcommand{\msbl}{m_{\tilde{b}_L}}
\newcommand{\msbe}{m_{\tilde{b}_1}}
\newcommand{\msbz}{m_{\tilde{b}_2}}
\newcommand{\msbi}{m_{\tilde{b}_i}}
\newcommand{\msbj}{m_{\tilde{b}_j}}
\newcommand{\MstL}{M_{\tilde{t}_L}}
\newcommand{\MstR}{M_{\tilde{t}_R}}
\newcommand{\MsbL}{M_{\tilde{b}_L}}
\newcommand{\MsbR}{M_{\tilde{b}_R}}
\newcommand{\MsqL}{M_{\tilde{q}_L}}
\newcommand{\MsqR}{M_{\tilde{q}_R}}
\newcommand{\Mtlr}{M_{t}^{LR}}
\newcommand{\Mtlrz}{\KL M_{t}^{LR}\KR^2}
\newcommand{\Mtlrd}{\KL M_{t}^{LR}\KR^3}
\newcommand{\Mtlrv}{\KL M_{t}^{LR}\KR^4}
\newcommand{\Mtlrf}{\KL M_{t}^{LR}\KR^5}
\newcommand{\Mtlrse}{\KL M_{t}^{LR}\KR^6}
\newcommand{\Mtlrsi}{\KL M_{t}^{LR}\KR^7}
\newcommand{\Mtlra}{\KL M_{t}^{LR}\KR^8}
\newcommand{\Mblr}{M_{b}^{LR}}
\newcommand{\Mblrz}{\KL M_{b}^{LR}\KR^2}
\newcommand{\Mqlr}{M_{q}^{LR}}
\newcommand{\Xt}{X_t}
\newcommand{\Xb}{X_b}
\newcommand{\At}{A_t}
\newcommand{\Ab}{A_b}

\newcommand{\mf}{m_f}
\newcommand{\mfp}{m_{f'}}
\newcommand{\mfi}{m_{f_i}}
\newcommand{\msfr}{m_{\tilde{f}_R}}
\newcommand{\msfl}{m_{\tilde{f}_L}}
\newcommand{\msf}{m_{\tilde{f}}}
\newcommand{\msfi}{m_{\tilde{f}_i}}
\newcommand{\msfe}{m_{\tilde{f}_1}}
\newcommand{\msfz}{m_{\tilde{f}_2}}

\newcommand{\ms}{M_S}
\newcommand{\msusy}{M_{SUSY}}
\newcommand{\msq}{m_{\tilde{q}}}
\newcommand{\msqi}{m_{\tilde{q}_i}}
\newcommand{\msqip}{m_{\tilde{q}_{i'}}}
\newcommand{\msqj}{m_{\tilde{q}_j}}
\newcommand{\msqe}{m_{\tilde{q}_1}}
\newcommand{\msqz}{m_{\tilde{q}_2}}

\newcommand{\sfn}{\tilde{f}}
\newcommand{\sfl}{\tilde{f}_L}
\newcommand{\sfr}{\tilde{f}_R}
\newcommand{\sfe}{\tilde{f}_1}
\newcommand{\sfz}{\tilde{f}_2}
\newcommand{\sfi}{\tilde{f}_i}
\newcommand{\sfj}{\tilde{f}_j}
\newcommand{\sfez}{\tilde{f}_{12}}
\newcommand{\sfze}{\tilde{f}_{21}}

\newcommand{\Pe}{\phi_1}
\newcommand{\Pz}{\phi_2}
\newcommand{\Pez}{\phi_{1,2}}
\newcommand{\PePz}{\phi_1\phi_2}
\newcommand{\mpe}{m_{\Pe}}
\newcommand{\mpz}{m_{\Pz}}
\newcommand{\mpez}{m_{\PePz}}

\newcommand{\smallz}{{\scriptscriptstyle Z}} %  a smaller Z
\newcommand{\smallD}{{\scriptscriptstyle D}} %
\newcommand{\smallQ}{{\scriptscriptstyle Q}} %
\newcommand{\smallU}{{\scriptscriptstyle U}} %

\newcommand{\SU}{\mathrm {SUSY}}
\newcommand{\SM}{\mathrm {SM}}
\newcommand{\msbar}{$\overline{\rm{MS}}$}
\newcommand{\oa}{{\cal O}(\alpha)}
\newcommand{\oas}{{\cal O}(\alpha_s)}
\newcommand{\oaas}{{\cal O}(\alpha\alpha_s)}
\newcommand{\oaasz}{{\cal O}(\alpha\alpha_s^2)}
\newcommand{\ogmzmts}{{\cal O}(\gf^2\mt^6)}
\def\order#1{${\cal O}(#1)$}
\newcommand{\cp}{{\cal CP}}
\newcommand{\wz}{\sqrt{2}}
\newcommand{\edz}{\frac{1}{2}}
\def\ed#1{\frac{1}{#1}}
\newcommand{\twol}{two-loop}
\newcommand{\onel}{one-loop}
\newcommand{\mma}{{\em Mathematica}}
\newcommand{\tc}{{\em TwoCalc}}
\newcommand{\fa}{{\em FeynArts}}
\newcommand{\fh}{{\em FeynHiggs}}
\newcommand{\fhf}{{\em FeynHiggsFast}}
\newcommand{\rp}{$\rho\,$-parameter}

\newcommand{\MW}{M_W}
\newcommand{\MZ}{M_Z}
\newcommand{\MA}{M_A}
\newcommand{\mh}{m_h}
\newcommand{\mH}{m_H}
\newcommand{\Mh}{M_h}
\newcommand{\MH}{M_H}
\newcommand{\mhp}{m_{H^\pm}}

\newcommand{\sudxsuzxue}{SU(3)_C \otimes SU(2)_L \otimes U(1)_Y}
\newcommand{\suzxue}{SU(2)_L \otimes U(1)_Y}

\newcommand{\dr}{\De\rho}
\newcommand{\mt}{m_{t}}
\newcommand{\mtms}{\overline{m}_t}
\newcommand{\mq}{m_{q}}
\newcommand{\mb}{m_{b}}
\newcommand{\gl}{\tilde{g}}
\newcommand{\Mgl}{m_{\tilde{g}}}
\newcommand{\mgl}{m_{\tilde{g}}}
\newcommand{\sq}{\tilde{q}}
\newcommand{\sqi}{\tilde{q}_i}
\newcommand{\sqj}{\tilde{q}_j}
\newcommand{\sql}{\tilde{q}_L}
\newcommand{\sqr}{\tilde{q}_R}
\newcommand{\sqe}{\tilde{q}_1}
\newcommand{\sqz}{\tilde{q}_2}
\newcommand{\sqez}{\tilde{q}_{12}}
\newcommand{\Stop}{\tilde{t}}
\newcommand{\StopL}{\tilde{t}_L}
\newcommand{\StopR}{\tilde{t}_R}
\newcommand{\Stope}{\tilde{t}_1}
\newcommand{\Stopz}{\tilde{t}_2}
\newcommand{\Stopi}{\tilde{t}_i}
\newcommand{\Stopj}{\tilde{t}_j}
\newcommand{\Sbot}{\tilde{b}}
\newcommand{\SbotL}{\tilde{b}_L}
\newcommand{\SbotR}{\tilde{b}_R}
\newcommand{\Sbote}{\tilde{b}_1}
\newcommand{\Sbotz}{\tilde{b}_2}
\newcommand{\Sboti}{\tilde{b}_i}
\newcommand{\Sbotj}{\tilde{b}_j}
\newcommand{\dst}{\Delta_{\tilde{t}}}
\newcommand{\tst}{\theta_{\tilde{t}}}
\newcommand{\tsb}{\theta_{\tilde{b}}}
\newcommand{\tsf}{\theta\kern-.20em_{\tilde{f}}}
\newcommand{\tsfp}{\theta\kern-.20em_{\tilde{f}\prime}}
\newcommand{\tsq}{\theta\kern-.15em_{\tilde{q}}}
\newcommand{\sw}{s_W}
\newcommand{\cw}{c_W}
\newcommand{\sweff}{\sin^2\theta_{\mathrm{eff}}}

\newcommand{\sintt}{\sin\tst}
\newcommand{\sinQtt}{\sin^2\tst}
\newcommand{\sinZtt}{\sin 2\tst}
\newcommand{\sinQZtt}{\sin^2 2\tst}
\newcommand{\sintb}{\sin\tsb}
\newcommand{\sinQtb}{\sin^2\tsb}
\newcommand{\sinZtb}{\sin 2\tsb}
\newcommand{\sintf}{\sin\tsf}
\newcommand{\sintfp}{\sin\tsfp}
\newcommand{\sinQtf}{\sin^2\tsf}
\newcommand{\sinZtf}{\sin 2\tsf}
\newcommand{\sintq}{\sin\tsq}
\newcommand{\sinQtq}{\sin^2\tsq}
\newcommand{\sinZtq}{\sin 2\tsq}

\newcommand{\costt}{\cos\tst}
\newcommand{\cosQtt}{\cos^2\tst}
\newcommand{\cosZtt}{\cos 2\tst}
\newcommand{\costb}{\cos\tsb}
\newcommand{\cosQtb}{\cos^2\tsb}
\newcommand{\cosZtb}{\cos 2\tsb}
\newcommand{\costf}{\cos\tsf}
\newcommand{\costfp}{\cos\tsfp}
\newcommand{\cosQtf}{\cos^2\tsf}
\newcommand{\cosZtf}{\cos 2\tsf}
\newcommand{\costq}{\cos\tsq}
\newcommand{\cosQtq}{\cos^2\tsq}
\newcommand{\cosZtq}{\cos 2\tsq}

\newcommand{\stt}{s_{\tilde{t}}}
\newcommand{\stb}{s_{\tilde{b}}}
\newcommand{\stf}{s_{\tilde{f}}}
\newcommand{\stq}{s_{\tilde{q}}}
\newcommand{\ctt}{c_{\tilde{t}}}
\newcommand{\ctb}{c_{\tilde{b}}}
\newcommand{\ctf}{c_{\tilde{f}}}
\newcommand{\ctq}{c_{\tilde{q}}}

\newcommand{\KL}{\left(}
\newcommand{\KR}{\right)}
\newcommand{\KKL}{\left[}
\newcommand{\KKR}{\right]}
\newcommand{\KKKL}{\left\{}
\newcommand{\KKKR}{\right\}}
\newcommand{\BL}{\lbrack}
\newcommand{\BR}{\rbrack}

\newcommand{\VL}{\left( \begin{array}{c}}
\newcommand{\VR}{\end{array} \right)}
\newcommand{\ML}{\left( \begin{array}{cc}}
\newcommand{\MLd}{\left( \begin{array}{ccc}}
\newcommand{\MLv}{\left( \begin{array}{cccc}}
\newcommand{\MR}{\end{array} \right)}
\newcommand{\dd}{\partial}
\newcommand{\dmu}{\partial_{\mu}}
\newcommand{\dmo}{\partial^{\mu}}
\newcommand{\mn}{{\mu\nu}}
\newcommand{\hc}{\mbox {h.c.}}
\newcommand{\re}{\mbox {Re}\,}
\newcommand{\OP}{\omega_+}
\newcommand{\OM}{\omega_-}
\newcommand{\Tb}{\tan \beta\hspace{1mm}}
\newcommand{\tb}{\tan \beta}
\newcommand{\TQb}{\tan^2 \beta\hspace{1mm}}
\newcommand{\CTb}{\cot \beta\hspace{1mm}}
\newcommand{\CTQb}{\cot^2 \beta\hspace{1mm}}
\newcommand{\Sb}{\sin \beta\hspace{1mm}}
\newcommand{\sbe}{\sin \beta}
\newcommand{\SQb}{\sin^2\beta\hspace{1mm}}
\newcommand{\SDb}{\sin^3\beta\hspace{1mm}}
\newcommand{\SVb}{\sin^4\beta\hspace{1mm}}
\newcommand{\Cb}{\cos \beta\hspace{1mm}}
\newcommand{\CQb}{\cos^2\beta\hspace{1mm}}
\newcommand{\CDb}{\cos^3\beta\hspace{1mm}}
\newcommand{\CVb}{\cos^4\beta\hspace{1mm}}
\newcommand{\Sa}{\sin \alpha\hspace{1mm}}
\newcommand{\SQa}{\sin^2\alpha\hspace{1mm}}
\newcommand{\Ca}{\cos \alpha\hspace{1mm}}
\newcommand{\CQa}{\cos^2\alpha\hspace{1mm}}
\newcommand{\Ta}{\tan\al}
\newcommand{\CTa}{\cot\al}
\newcommand{\Sab}{\sin (\alpha + \beta)\hspace{1mm}}
\newcommand{\Cab}{\cos (\alpha + \beta)\hspace{1mm}}
\newcommand{\Sba}{\sin (\beta - \alpha)\hspace{1mm}}
\newcommand{\Cba}{\cos (\beta - \alpha)\hspace{1mm}}
\newcommand{\SZa}{\sin 2\alpha\hspace{1mm}}
\newcommand{\SQZa}{\sin^2 2\alpha\hspace{1mm}}
\newcommand{\CZa}{\cos 2\alpha\hspace{1mm}}
\newcommand{\CQZa}{\cos^2 2\alpha\hspace{1mm}}
\newcommand{\SZb}{\sin 2\beta\hspace{1mm}}
\newcommand{\SQZb}{\sin^2 2\beta\hspace{1mm}}
\newcommand{\CZb}{\cos 2\beta\hspace{1mm}}
\newcommand{\CQZb}{\cos^2 2\beta\hspace{1mm}}
\newcommand{\CDZb}{\cos^3 2\beta\hspace{1mm}}
\newcommand{\CVZb}{\cos^4 2\beta\hspace{1mm}}
\newcommand{\CQba}{\cos^2 (\beta - \alpha)\hspace{1mm}}
\newcommand{\SQba}{\sin^2 (\beta - \alpha)\hspace{1mm}}
\newcommand{\SDba}{\sin^3 (\beta - \alpha)\hspace{1mm}}
\newcommand{\CQab}{\cos^2 (\alpha + \beta)\hspace{1mm}}
\newcommand{\SQab}{\sin^2 (\alpha + \beta)\hspace{1mm}}
\newcommand{\SDab}{\sin^3 (\alpha + \beta)\hspace{1mm}}
\newcommand{\SDea}{\sin \De\al}
\newcommand{\CDea}{\cos \De\al}
\newcommand{\CQDea}{\cos^2\De\al}
\newcommand{\TDea}{\tan \De\al}
\newcommand{\TQDea}{\tan^2\De\al}

\newcommand{\tev}{\,\, \mathrm{TeV}}
\newcommand{\gev}{\,\, \mathrm{GeV}}
\newcommand{\mev}{\,\, \mathrm{MeV}}

\newcommand{\BC}{\begin{center}}
\newcommand{\EC}{\end{center}}
\newcommand{\BE}{\begin{equation}}
\newcommand{\EE}{\end{equation}}
\newcommand{\BEA}{\begin{eqnarray}}
\newcommand{\BEAnn}{\begin{eqnarray*}}
\newcommand{\EEA}{\end{eqnarray}}
\newcommand{\EEAnn}{\end{eqnarray*}}
\newcommand{\non}{\nonumber}
\newcommand{\id}{{\rm 1\kern-.12em
\rule{0.3pt}{1.5ex}\raisebox{0.0ex}{\rule{0.1em}{0.3pt}}}}
\newcommand{\lsim}
{\;\raisebox{-.3em}{$\stackrel{\displaystyle <}{\sim}$}\;}
\newcommand{\gsim}
{\;\raisebox{-.3em}{$\stackrel{\displaystyle >}{\sim}$}\;}


\newcommand{\cf}{C_F}
\newcommand{\gf}{G_F}
\newcommand{\nf}{N_f}
% shorthands for greek letters
\def\al{\alpha}
\def\aeff{\al_{\rm eff}}
\def\aeffapprox{\al_{\rm eff}{\rm (approx)}}
\newcommand{\Saeff}{\sin\aeff}
\newcommand{\Caeff}{\cos\aeff}
\def\als{\alpha_s}
\def\altl{\alpha_{\rm 2-loop}}
\def\gs{g_s}
\def\be{\beta}
\def\ben{\be_0}
\def\bee{\be_1}
\def\bez{\be_2}
\def\ga{\gamma}
\def\gan{\ga_0}
\def\gae{\ga_1}
\def\gaz{\ga_2}
\def\de{\delta}
\def\deqed{\de\Ga_\ga}
\def\eps{\varepsilon}
\def\la{\lambda}
\def\si{\sigma}
\def\ie{i\epsilon}
\def\om{\omega}
\def\Ga{\Gamma}
\def\Gh{\Ga_h}
\def\GH{\Ga_H}
\def\Gz{\Ga_0}
\def\Ge{\Ga_1}
\def\Geeff{\Ga_{1,\rm eff}}
\def\Geg{\Ga_{1,g}}
\def\Gga{\Ga_\ga}
\def\Gg{\Ga_g}
\def\Ggl{\Ga_{\gl}}
\def\Gtot{\Ga_{\rm tot}}
\def\De{\Delta}
\def\La{\Lambda}
\def\Laqcd{\La_{\rm QCD}}
\def\Si{\Sigma}
\def\Sip{\Sigma'}
\def\Sie{\Sigma^{(1)}}
\def\Siz{\Sigma^{(2)}}
\def\Sii{\Sigma^{(i)}}
\def\Siep{\Sigma'^{(1)}}
\def\Sizp{\Sigma'^{(2)}}
\def\Siip{\Sigma'^{(i)}}

\def\hSi{\hat{\Sigma}}
\def\hSip{\hat{\Sigma}'}
\def\hSie{\hat{\Sigma}^{(1)}}
\def\hSiz{\hat{\Sigma}^{(2)}}
\def\hSiep{\hat{\Sigma}^{(1)}'}
\def\hSizp{\hat{\Sigma}^{(2)}'}
\def\hSiH{\hSi_{H}}
\def\hSih{\hSi_{h}}
\def\hSihH{\hSi_{hH}}
\def\hSipH{\hSip_{H}}
\def\hSiph{\hSip_{h}}
\def\hSiphH{\hSip_{hH}}
\def\hSieh{\hSie_{h}}
\def\hSieH{\hSie_{H}}
\def\hSiehH{\hSie_{hH}}
\def\hSizh{\hSiz_{h}}
\def\hSizH{\hSiz_{H}}
\def\hSizhH{\hSiz_{hH}}
\newcommand{\htad}{\hat{t}}
\newcommand{\tade}{t^{(1)}}
\newcommand{\tadz}{t^{(2)}}

\newcommand{\Ph}{\mathswitchr h}
\newcommand{\PH}{\mathswitchr H}

\newcommand{\deMZe}{\de \MZ^{2\,(1)}}
\newcommand{\deMZz}{\de \MZ^{2\,(2)}}
\newcommand{\deMZi}{\de \MZ^{2\,(i)}}
\newcommand{\deMAe}{\de \MA^{2\,(1)}}
\newcommand{\deMAz}{\de \MA^{2\,(2)}}
\newcommand{\deMAi}{\de \MA^{2\,(i)}}
\newcommand{\deZHie}{\de Z_{H_i}^{(1)}}
\newcommand{\deZHiz}{\de Z_{H_i}^{(2)}}
\newcommand{\deZHee}{\de Z_{H_1}^{(1)}}
\newcommand{\deZHez}{\de Z_{H_1}^{(2)}}
\newcommand{\deZHze}{\de Z_{H_2}^{(1)}}
\newcommand{\deZHzz}{\de Z_{H_2}^{(2)}}
\newcommand{\deTbe}{\de\Tb^{(1)}}
\newcommand{\deTbz}{\de\Tb^{(2)}}
\newcommand{\detie}{\de t_i^{(1)}}
\newcommand{\detiz}{\de t_i^{(2)}}
\newcommand{\detee}{\de t_1^{(1)}}
\newcommand{\detez}{\de t_1^{(2)}}
\newcommand{\detei}{\de t_1^{(i)}}
\newcommand{\detze}{\de t_2^{(1)}}
\newcommand{\detzz}{\de t_2^{(2)}}
\newcommand{\detzi}{\de t_2^{(i)}}
\newcommand{\deVie}{\de V_i^{(1)}}
\newcommand{\deViz}{\de V_i^{(2)}}
\newcommand{\deVee}{\de V_{\Pe}^{(1)}}
\newcommand{\deVez}{\de V_{\Pe}^{(2)}}
\newcommand{\deVei}{\de V_{\Pe}^{(i)}}
\newcommand{\deVze}{\de V_{\Pz}^{(1)}}
\newcommand{\deVzz}{\de V_{\Pz}^{(2)}}
\newcommand{\deVzi}{\de V_{\Pz}^{(i)}}
\newcommand{\deVeze}{\de V_{\PePz}^{(1)}}
\newcommand{\deVezz}{\de V_{\PePz}^{(2)}}
\newcommand{\deVezi}{\de V_{\PePz}^{(i)}}


\newcommand{\SLASH}[2]{\makebox[#2ex][l]{$#1$}/}
\newcommand{\Dslash}{\SLASH{D}{.5}\,}
\newcommand{\dslash}{\SLASH{\dd}{.15}}
\newcommand{\kslash}{\SLASH{k}{.15}}
\newcommand{\pslash}{\SLASH{p}{.2}}
\newcommand{\qslash}{\SLASH{q}{.08}}

\newcommand{\sqmsmt}{\sqrt{\ms^2-\mt^2}}
\newcommand{\sqmsmtbar}{\sqrt{\ms^2-\mtms^2}}
\newcommand{\lmtms}{\KL\frac{\mt^2}{\ms^2}\KR}
\newcommand{\lmtmsms}{\KL\frac{\mtms^2}{\ms^2}\KR}
\newcommand{\lmsmt}{\KL\frac{\ms^2 - \mt^2}{\mt^2}\KR}

\newcommand{\lc}{linear collider}
\newcommand{\epem}{$e^+e^-$}
\newcommand{\ipb}{$\mbox{pb}^{-1}$}
\newcommand{\ifb}{$\mbox{fb}^{-1}$}
\newcommand{\iab}{$\mbox{ab}^{-1}$}
\newcommand{\gvf}{g_V^f}
\newcommand{\gaf}{g_A^f}

\newcommand{\hff}{h \to f\bar{f}}
\newcommand{\Hff}{H \to f\bar{f}}
\newcommand{\hbb}{h \to b\bar{b}}
\newcommand{\Hbb}{H \to b\bar{b}}
\newcommand{\htautau}{h \to \tau^+\tau^-}
\newcommand{\hcc}{h \to c\bar{c}}
\newcommand{\rf}{R_f}
%\newcommand{\cff}{C_{ff}}
\newcommand{\cff}{C_f^{(u)}}
%\newcommand{\cfpfp}{C_{f'f'}}
\newcommand{\cfpfp}{C_f^{(d)}}
\newcommand{\Zh}{Z_h}
\newcommand{\wZh}{\sqrt{\Zh}}
\newcommand{\ZH}{Z_H}
\newcommand{\wZH}{\sqrt{\ZH}}
\newcommand{\ZhH}{Z_{hH}}
\newcommand{\ZHh}{Z_{Hh}}


