\documentclass[12pt,twoside]{article}

\usepackage{epsf,cite} 
\usepackage{colordvi}

\textwidth       16.0cm
\textheight      25.0cm
\evensidemargin   0.0cm 
\oddsidemargin    0.0cm 
\topmargin       -1.5cm 
\headsep          1.0cm 
\hoffset          0.0cm 
\voffset          0.0cm 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%  Input macros for SUSY and such 
\input{SH_macros.tex}
\input{SA_macros.tex}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%% Beginning of text %%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{document}
\pagestyle{myheadings} 
\markboth
{\it Aspects of GMSB Phenomenology at TeV Colliders}
{\it Aspects of GMSB Phenomenology at TeV Colliders}

\begin{titlepage}
\thispagestyle{empty} 
\null~\vspace{-3.0cm} 
\begin{flushright} 
                           CERN-TH/2000-054 \\
     \hfill  DESY 00/021 \\ 
\end{flushright} 
\hrule\hfill

\vspace{-0.1cm}

\hrule\hfill

\vspace{0.2cm} 
                                        
\begin{center}
{\BrickRed{\Large \bf 
Aspects of GMSB Phenomenology at TeV Colliders}} \\

\vspace{0.2cm} 

\parbox{0.75\textwidth}{\hrule\hfill}

\vspace{0.2cm}

{\Large \sl
Report of the GMSB SUSY Working Group} \\
\vspace{0.2cm} 
{\Large \sl 
Workshop ``Physics at TeV Colliders''} \\
\vspace{0.2cm}
{\Large \sl
Les Houches, 7-18 June 1999} \\ 

\vspace{0.2cm} 

\parbox{0.75\textwidth}{\hrule\hfill}

\vspace{0.5cm}

{\large {\bf S.~Ambrosanio}~$^{a,\star}$ (convener),  
        {\bf S.~Heinemeyer}~$^{b}$,
        {\bf B.~Mele}~$^{c}$,
        {\bf S.~Petrarca}~$^{c}$,
        {\bf G.~Polesello}~$^{d}$,
        {\bf A.~Rimoldi}~$^{d,e}$,
        {\bf G.~Weiglein}~$^{a}$} \\
~\\
$^a$ CERN -- {\it Theory Division,
     CH-1211 Geneva 23, Switzerland}\\
~\\
$^b$ DESY -- {\it Theory Group, 
     Notkestr. 85, D-22603 Hamburg, Germany}\\
~\\
$^c$ INFN -- {\it Sezione di Roma I e Dipartimento di Fisica, 
     Universit\`a ``La Sapienza'',}\\
     {\it p.~le Aldo Moro 2, I-00185 Roma, Italy}\\
~\\
$^d$ INFN -- {\it Sezione di Pavia e Dip.~Fis.~Nucl.~Teor.,
     Universit\`a di Pavia,} \\
     {\it via Bassi 6, I-27100 Pavia, Italy}\\
~\\
$^e$ CERN -- {\it EP Division,
     CH-1211 Geneva 23, Switzerland} \\
\end{center}


\vspace*{\fill}  

\begin{center} 
{\bf Abstract} \\
\end{center}
{\small 
The status of two on-going studies concerning important aspects of the 
phenomenology of gauge-mediated supersymmetry breaking (GMSB) 
models at TeV colliders is reported. \\
The first study deals with the characteristics of the light Higgs 
boson spectrum allowed by the (minimal and non-minimal) GMSB framework.
Today's most accurate GMSB model generation and two-loop 
Feynman-diagrammatic calculation of $\mh$ have been combined. 
The Higgs masses are shown in dependence of various model parameters
at the messenger and electroweak scales. In the minimal model, an upper
limit on $\mh$ of about 124 GeV is found for $\mt = 175$~GeV. \\
The second study is focused on the measurement of the fundamental
SUSY breaking scale $\sqrt{F}$ at the LHC in the GMSB scenario where
a stau is the next-to-lightest SUSY particle (NLSP) and decays into a 
gravitino with $c\tau_{\rm NLSP}$ in the range 0.5~m to 1~km.
This implies the measurement of mass and lifetime of long lived sleptons.
The identification is performed by determining the time of flight in the 
ATLAS muon chambers. Accessible range and precision on $\sqrt{F}$ 
achievable with a counting method are assessed.
}

\vspace*{\fill}  

\noindent 
\parbox{0.4\textwidth}{\hrule\hfill} \\ 
{\small 
$^\star$\ E-mail: {\tt ambros@mail.cern.ch} 
}  

\end{titlepage}

\setcounter{page}{0} 

\thispagestyle{empty} 
~\\ 
\newpage 
\thispagestyle{plain} 

\section{Introduction to Gauge-Mediated SUSY Breaking}
\label{sec:intro}
\noindent 
Since no superpartners have been detected at collider experiments so far,
supersymmetry (SUSY) cannot be an exact symmetry of Nature.
The requirement of ``soft'' supersymmetry breaking~\cite{StevePrimer} 
alone is not sufficient to reduce the free parameters
to a number suitable for predictive phenomenological studies. 
Hence, motivated theoretical hypotheses on the nature of SUSY breaking 
and the mechanism through which it is transmitted to the visible sector 
of the theory [here assumed to be the one predicted by the minimal SUSY 
extension of the standard model (MSSM)] are highly desirable. 
If SUSY is broken at energies of the order of the Planck mass and the SUSY 
breaking sector communicates with the MSSM sector through gravitational 
interactions only, one falls in the supergravity-inspired (SUGRA) scheme.
The most recognised alternative to SUGRA is based instead on the hypothesis
that SUSY breaking occurs at relatively low energy scales and it is mediated 
mainly by gauge interactions (GMSB)~\cite{oldGMSB,newGMSB,GR-GMSB}. 
A good theoretical reason to consider such a possibility is that it 
provides a natural, automatic suppression of the SUSY contributions 
to flavour-changing neutral current and $\cp$-violating processes. 
A pleasant consequence is that, at least in the simplest versions of GMSB, 
the MSSM spectrum and other observables depend on just a handful of 
parameters, typically

\be
M_{\rm mess}, \; N_{\rm mess}, \; \Lambda, \; \tb, \; {\rm sign}(\mu),
\label{eq:pars}
\ee

\noindent
where $M_{\rm mess}$ is the overall messenger scale; $N_{\rm mess}$ is the 
so-called messenger index, parameterising the structure of the messenger
sector; $\Lambda$ is the universal soft SUSY breaking scale felt by the
low-energy sector; $\tb$ is the ratio of the vacuum expectation 
values of the two Higgs doublets; sign($\mu$) is the ambiguity
left for the SUSY higgsino mass after conditions for correct 
electroweak symmetry breaking (EWSB) are imposed (see e.g.  
Refs.~\cite{GMSBmodels1,GMSBmodels2,AKM-LEP2,AB-LC}).

The phenomenology of GMSB (and more in general of any theory with low-energy
SUSY breaking) is characterised by the presence of a very light gravitino 
$\G$~\cite{Fayet},

\be 
m_{3/2} \equiv m_{\G} = \frac{F}{\sqrt{3}M'_P} \simeq 
\left(\frac{\sqrt{F}}{100 \; {\rm TeV}}\right)^2 2.37 \; {\rm eV},  
\label{eq:Gmass}
\ee

\noindent
where $\sqrt{F}$ is the fundamental scale of SUSY breaking, 100 TeV is a 
typical value for it, and $M'_P = 2.44 \times 10^{18}$ GeV is the reduced 
Planck mass.
Hence, the $\G$ is always the lightest SUSY particle (LSP) in these theories. 
If $R$-parity is assumed to be conserved, any produced MSSM particle will 
finally decay into the gravitino. Depending on $\sqrt{F}$, the interactions 
of the gravitino, although much weaker than gauge and Yukawa interactions, 
can still be strong enough to be of relevance for collider physics. 
As a result, in most cases the last step of any SUSY decay chain is 
the decay of the next-to-lightest SUSY particle (NLSP), which can  
occur outside or inside a typical detector or even close to the interaction 
point. The pattern of the resulting spectacular signatures is determined 
by the identity of the NLSP and its lifetime before decaying into the $\G$,

\be 
c \tau_{\rm NLSP} \simeq \frac{1}{100 {\cal B}}
\left(\frac{\sqrt{F}}{100 \; {\rm TeV}}\right)^4 
\left(\frac{m_{\rm NLSP}} {100 \; {\rm GeV}}\right)^{-5},
\label{eq:NLSPtau}
\ee 

\noindent
where ${\cal B}$ is a number of order unity depending on the nature
of the NLSP.

The identity of the NLSP [or, to be more precise, the identity of the 
sparticle(s) having a large branching ratio (BR) for decaying into the 
gravitino and the relevant SM partner] determines four main scenarios 
giving rise to qualitatively different phenomenology:

\begin{description} 

\item[Neutralino NLSP scenario:] Occurs whenever 
$m_{\NI} < (m_{\tauu} - m_{\tau})$. Here typically a decay of the $\NI$  
to $\G\gamma$ is the final step of decay chains following
any SUSY production process. As a consequence, the main inclusive signature 
at colliders is prompt or displaced photon pairs + X + missing energy. 
$\NI$ decays to $\G \Z$ and other minor channels may also be relevant 
at TeV colliders.    

\item[Stau NLSP scenario:] Defined by 
$m_{\tauu} < {\rm Min}[m_{\NI}, m_{\lR}] - m_{\tau}$, 
features $\tauu \to \G \tau$ decays, producing $\tau$ pairs or 
charged semi-stable $\tauu$ tracks or decay kinks + X + missing energy.
Here and in the following, $\ell$ stands for $e$ or $\mu$. 

\item[Slepton co-NLSP scenario:] When 
$m_{\lR} < {\rm Min}[m_{\NI}, m_{\tauu} + m_{\tau}]$, 
$\lR \to \G \ell$ decays are also open with large BR. In addition to the 
signatures of the stau NLSP scenario, one also gets $\ell^+\ell^-$
pairs or $\lR$ tracks or decay kinks. 

\item[Neutralino-stau co-NLSP scenario:] If 
$| m_{\tauu} - m_{\NI} | < m_{\tau}$ and $m_{\NI} < m_{\lR}$, 
both signatures of the neutralino NLSP and stau NLSP scenario are present
at the same time, since $\NI \leftrightarrow \tauu$ decays are not allowed
by phase space. 

\end{description}

Note that in the GMSB parameter space the relation $m_{\lR} > m_{\tauu}$ 
always holds. 
Also, one should keep in mind that the classification above is only 
valid as an indicative scheme in the limit $m_e$, $m_\mu \to 0$, neglecting 
also those cases where a fine-tuned choice of $\sqrt{F}$ and the sparticle 
masses may give rise to competition between phase-space suppressed decay 
channels from one ordinary sparticle to another and sparticle decays to the 
gravitino~\cite{AKM2}. 

In this report, we treat two important aspects of the GMSB phenomenology 
at TeV colliders:
\begin{description}
\item[(A)] The consequences of the GMSB hypotheses on the light Higgs spectrum
      using the most accurate tools available today for model generation and 
      $\mh$ calculation;
\item[(B)] Studies and possible measurements at the LHC with the ATLAS detector
      in the stau NLSP or slepton co-NLSP scenarios, with focus 
      on determining the fundamental SUSY breaking scale $\sqrt{F}$.
\end{description}

For this purpose, we generated about 30000 GMSB models under
well defined hypotheses, using the program {\tt SUSYFIRE}~\cite{SUSYFIRE}, 
as described in the following section. 

\section{GMSB Models}
\label{sec:models}

\noindent 
In the GMSB framework, the pattern of the MSSM spectrum is simple,
as all sparticle masses are generated in the same way and scale 
approximately with a single parameter $\Lambda$, which sets
the amount of soft SUSY breaking felt by the visible sector. 
As a consequence, scalar and gaugino masses are related to each
other at a high energy scale, which is not the case in other SUSY 
frameworks, e.g.\ SUGRA. Also, it is possible to impose other 
conditions at a lower scale to achieve EWSB and  
further reduce the dimension of the parameter space. 

To build our GMSB models, we adopt the usual phenomenological approach,
in particular following Ref.~\cite{AB-LC}, where problems relevant 
for GMSB physics at TeV colliders were also approached.
We do not specify the origin of the SUSY higgsino mass $\mu$, nor do we 
assume that the analogous soft SUSY breaking parameter $B\mu$ vanishes at the 
messenger scale. Instead, we impose correct EWSB to trade $\mu$
and $B\mu$ for $M_Z$ and $\tb$, leaving the sign of $\mu$ undetermined. 
However, we are aware that to build a satisfactory GMSB model one should also 
solve the latter problem in a more fundamental way, perhaps by providing a 
dynamical mechanism to generate $\mu$ and $B\mu$, possibly with values
of the same order of magnitude. This might be accomplished radiatively 
through some new interactions. However, in this case the other
soft terms in the Higgs potential, namely $m^2_{H_{1,2}}$, will be also  
affected and this will in turn change the values of $|\mu|$ and $B\mu$ 
coming from EWSB conditions~\cite{GR-GMSB,GMSBmodels1,GMSBmodels2}. 
Within the study {\bf (A)}, we are currently considering some ``non-minimal'' 
possibilities for GMSB models that to some extent take this problem into 
account, and we are trying to assess the impact on the light Higgs mass. 
We do not treat this topic here, but refer to~\cite{AHW} for further details.

To determine the MSSM spectrum and low-energy parameters, we solve 
the renormalisation group equation (RGE) evolution with the following 
boundary conditions at the $M_{\rm mess}$ scale,

\bea
M_a & = & N_{\rm mess} \Lambda 
g\left(\frac{\Lambda}{M_{\rm mess}}\right) \frac{\alpha_a}{4\pi}, 
\; \; \; (a=1, 2, 3) \nonumber \\
\tilde{m}^2 & =  & 2 N_{\rm mess} \Lambda^2 
f\left(\frac{\Lambda}{M_{\rm mess}}\right) 
\sum_a \left(\frac{\alpha_a}{4\pi}\right)^2 C_a,
\label{eq:bound}
\eea 

\noindent
respectively for the gaugino and the scalar masses. 
In Eq.~(\ref{eq:bound}), $g$ and $f$ are the one-loop and two-loop functions 
whose exact expressions can be found e.g. in Ref.~\cite{AKM-LEP2}, 
and $C_a$ are the quadratic Casimir invariants for the scalar fields.
As usual, the scalar trilinear couplings $A_f$ are assumed to vanish
at the messenger scale, as suggested by the fact that they (and not
their squares) are generated via gauge interactions with the messenger 
fields at the two loop-level only. 

 To single out the interesting region of the GMSB parameter space, 
we proceed as follows.
Barring the case where a neutralino is the NLSP and decays outside
the detector (large $\sqrt{F}$), the GMSB signatures are very spectacular
and are generally free from SM background. Keeping this in mind and being 
interested in GMSB phenomenology at future TeV colliders, we consider 
only models where the NLSP mass is larger than 100 GeV, assuming that 
searches at LEP and Tevatron, if unsuccessful, will in the end exclude 
a softer spectrum in most cases. 
We require that $M_{\rm mess} > 1.01 \Lambda$, to prevent an excess
of fine-tuning of the messenger masses, and that the mass of the 
lightest messenger scalar be at least 10 TeV. We also impose 
$M_{\rm mess} > M_{\rm GUT} \; {\rm exp}(-125/N_{\rm mess})$,
to ensure perturbativity of gauge interactions up to the 
GUT scale. Further, we do not consider models with $M_{\rm mess} \gtap 10^{5}
\Lambda$. As a result of this and other constraints, the messenger index 
$N_{\rm mess}$, which we assume to be an integer independent of the gauge 
group, cannot be larger than~8. To prevent the top Yukawa coupling from 
blowing up below the GUT scale, we require $\tb > 1.2$ (and in some 
cases $> 1.5$). This is also motivated by the current bounds from SUSY Higgs 
searches at LEP~II~\cite{mhiggsmSUGRA}.
Models with $\tb \gtap 55$ (with a mild dependence on $\Lambda$) are 
forbidden by the EWSB requirement and typically fail in giving $m_A^2 > 0$.

To calculate the NLSP lifetime relevant to our study {\bf (B)}, one needs 
to specify the value of the fundamental SUSY breaking scale $\sqrt{F}$ on a
model-by-model basis. Using perturbativity arguments, for each given 
set of GMSB parameters, it is possible to determine a lower bound
according to Ref.~\cite{AKM-LEP2}, 

\be 
\sqrt{F} \ge \sqrt{F_{\rm mess}} \equiv \sqrt{\Lambda M_{\rm mess}} > \Lambda.
\label{eq:sqrtFmin}
\ee

On the contrary, no solid arguments can be used to set an upper limit 
on $\sqrt{F}$ of relevance for collider physics, although some 
semi-qualitative cosmological arguments are sometimes evoked.

To generate our model samples using {\tt SUSYFIRE}, we used logarithmic
steps for $\Lambda$ (between about 45 TeV/$N_{\rm mess}$ and about 220 
TeV/$\sqrt{N_{\rm mess}}$, which corresponds to excluding models with
sparticle masses above $\sim 4$ TeV), $M_{\rm mess}/\Lambda$ (between
1.01 and $10^5$) and $\tb$ (between 1.2 and about 60), subject
to the constraints described above. 
{\tt SUSYFIRE} starts from the values of SM particle masses and gauge
couplings at the weak scale and then evolves up to the messenger
scale through RGE's. At the messenger scale, it imposes the boundary 
conditions (\ref{eq:bound}) for the soft sparticle masses and then 
evolves the full set of RGE's back to the weak scale. The decoupling 
of each sparticle at the proper threshold is taken into account. 
Two-loop RGE's are used for gauge couplings, third generation Yukawa
couplings and gaugino soft masses. The other RGE's are taken at the
one-loop level. At the scale $\sqrt{m_{\stopu}m_{\stopd}}$, EWSB conditions 
are imposed by means of the one-loop effective potential approach, 
including corrections from stops, sbottoms and staus. 
The program then evolves up again to $M_{\rm mess}$ and so on.
Three or four iterations are usually enough to get a good approximation
for the MSSM spectrum. 

\vspace{1.0cm} 

\hrule\hfill

\vspace{0.5cm} 

\noindent 
{\Large {\bf (A) The Light Higgs Boson Spectrum}} \\
{\Large {\bf \null~~~~~~in GMSB Models}} \\

\vspace{0.2cm} 

\hrule\hfill 

\vspace{0.2cm} 

\noindent
{\large Contribution by:} \\
{\large S.~Ambrosanio, S.~Heinemeyer, G.~Weiglein}

\vspace{0.5cm}

\section*{A.1~~Introduction}
\label{sec:A1}

\noindent
Within the MSSM, the masses of the $\cp$-even neutral Higgs bosons are 
calculable in terms of the other low-energy parameters. The mass of the
lightest Higgs boson, $\mh$, has been of particular interest, as
it is bounded to be smaller than the $\Z$~boson mass at the tree level. 
The \onel\ results~\cite{mhiggs1l,mhiggsf1l,mhiggsf1ldab,pierce} 
for $\mh$ have been supplemented in the
last years with the leading \twol\ corrections, performed in the
renormalisation group (RG)
approach~\cite{mhiggsRG1,mhiggsRG2}, in the effective
potential approach~\cite{mhiggsEP} and most recently in
the Feynman-diagrammatic (FD)
approach~\cite{mhiggsletter,mhiggslong,mhiggslle}. 
The \twol\ corrections have turned out to be sizeable. They can
lower the \onel\ results by up to 20\%.
These calculations predict an upper bound on $\mh$ of about 
$\mh \le 130$ GeV for an unconstrained MSSM with $\mt = 175$ GeV 
and a common SUSY mass scale $\msusy \le 1$ TeV.

As discussed in Sec.~\ref{sec:intro}, the GMSB scenario provides a 
relatively simple set of constraints and thus constitutes a very 
predictive and readily testable realization of the MSSM.
The main goal of the present analysis is to study the spectrum of the
lightest neutral $\cp$-even Higgs boson, $\mh$, within the GMSB framework. 
Particular emphasis is given to the maximal value of $\mh$ achievable in 
GMSB after an exhaustive scanning of the parameter space. 
Our results are discussed in terms of the GMSB constraints on the 
low-energy parameters and compared to the cases of a SUGRA-inspired or an 
unconstrained MSSM. 

\section*{A.2~~Calculation of $\mh$}
\label{sec:A2}

\noindent
To evaluate $\mh$, we employ the currently most accurate calculation 
based on the FD approach~\cite{mhiggsletter,mhiggslong,mhiggslle}. 
The most important radiative corrections to $\mh$ arise from the top and
scalar top sector of the MSSM, with the input parameters $\mt$, 
the masses of the scalar top quarks, $\mste$, $\mstz$, and the
$\Stop$-mixing angle, $\tst$. Here we adopt the conventions
introduced in Ref.~\cite{mhiggslong}.
The complete diagrammatic \onel\ result~\cite{mhiggsf1ldab} has been 
combined with the dominant \twol\ corrections of 
$\oaas$~\cite{mhiggsletter,mhiggslong} and with the subdominant 
corrections of $\ogmzmts$~\cite{mhiggsRG1,mhiggsRG2}. 
%and the leading
%contributions of \order{\al\als^2}~\cite{mhiggslong}. 
GMSB models are generated with the program {\tt SUSYFIRE}, according 
to the discussion of Sec.~\ref{sec:models}. For this study, we consider 
only models with $\tb > 1.5$~\cite{mhiggsmSUGRA} 
and $m_A > 80$ GeV~\cite{lepc}. 
In addition, we always use $\mt = 175$ GeV. A change of 1 GeV in $\mt$ 
translates roughly into a shift of 1 GeV (with the same sign) in $\mh$ 
as well. Thus, changing $\mt$ affects our results on $\mh$ in an 
easily predictable way. 

The results of the $\mh$ calculation have been implemented in the 
program \fh\ \cite{feynhiggs}. This {\tt Fortran} code has been combined 
with {\tt SUSYFIRE}, which has been used to calculate the low 
energy parameters $\mste$, $\mstz$, $\tst$, $\mu$, $M_1$, $M_2$, $\mgl$, 
$\dots$ for each of the $\sim$30000 GMSB models generated. 
These have then been passed to \fh\ for the $\mh$ evaluation in a
coherent way. Indeed, we transform the \msbar\ parameters in the 
{\tt SUSYFIRE} output into on-shell parameters before feeding
them into \fh. Compact expressions for the relevant transition 
formulas can be found in Refs.~\cite{m20,bse}.

Compared to an existing analysis in the GMSB framework~\cite{kaeding},
we use a more complete evaluation of $\mh$. This leads in particular
to smaller values of $\mh$ for a given set of input parameters in our 
analysis. Also, in Ref.~\cite{kaeding} although some GMSB scenarios with 
generalised messenger sectors were considered, the parameter space for the 
``minimal'' case with a unique, integer messenger index 
$N_{\rm mess} = N_1 = N_2 = N_3$ was not fully explored. 
Indeed, $\Lambda$ was in most cases limited to values smaller 
than 100 TeV and $M_{\rm mess}$ was fixed to $10^5$ TeV. Furthermore, 
partly as a consequence of the above assumptions, the authors did not 
consider models with $N_{\rm mess} > 4$, i.e. their requirements  
for perturbativity of the MSSM gauge couplings up to the GUT scale were 
stronger than ours. We will see in the following section that maximal 
$\mh$ values in our analysis are instead obtained for larger values of 
the messenger scale and the messenger index. 

\section*{A.3~~The Light Higgs Spectrum in GMSB}
\label{sec:A3}

\noindent 
In the following, we give some results in the form of scatter plots showing 
the pattern in GMSB for $\mh$, $m_A$ as well as other low-energy parameters
of relevance for the light Higgs spectrum. 

\begin{figure} 
\hspace{-0.5cm}
\epsfxsize=0.55\textwidth 
\epsffile{mhvstb.eps}
\epsfxsize=0.55\textwidth 
\epsffile{mavstb.eps}
\caption{\sl Scatter plots for the light scalar (a) and pseudoscalar (b)
Higgs masses as functions of $\tb$. Only GMSB models with 
$\tb > 1.5$, $m_A > 80$ GeV and $m_{\rm NLSP} > 100$ GeV are considered.
}
\label{fig:mh-mavstb}
\end{figure}

In Fig.~\ref{fig:mh-mavstb}(a), we show the dependence of $\mh$ on $\tb$,
where only models with $\tb > 1.5$, $m_A > 80$ GeV and
$m_{\rm NLSP} > 100$ GeV are considered, while $m_t$ is fixed to 175 GeV.
The dependence is strong for small $\tb \ltap 10$, while for larger 
$\tb$ the increase of the lightest Higgs mass is rather mild.
The maximum values for $\mh \simeq 124$ GeV are achieved for $\tb > 50$.  
It should be noted that for very large $\tb \gtap 52$, we also find a few
models with relatively small $\mh \ltap 100$ GeV. This is due
to the fact that in this case EWSB conditions tend to drive $m_A$ toward
very small values (cfr. Sec.~\ref{sec:models}). 
This is made visible by the scatter plot in Fig.~\ref{fig:mh-mavstb}(b), 
where the pseudoscalar Higgs mass is shown as a function of $\tb$. 
For such small values of $m_A$ and for large $\tb$, the relation 
$\mh \approx m_A$ holds. Thus small $\mh$ values are quite natural in 
this region of the parameter space. On the other hand, one can see that 
extremely large values of $m_A \gtap 2$ TeV can only be obtained for small 
or moderate $\tb \ltap 10$ GeV.
A comparison between Fig~\ref{fig:mh-mavstb}(a) and (b) reveals that the 
largest $\mh$ values $\gtap 123$ GeV correspond in GMSB to $m_A$ values in 
the 300--800 GeV range. 
Indeed, it has been checked that such large $\mh$ values are in general 
obtained in the FD calculation for 
$300 \ltap m_A \ltap 1000$ GeV, see Ref.~\cite{mhiggslong}.

\begin{figure} 
\hspace{3.0cm}
\epsfxsize=0.6\textwidth 
\epsffile{mhvsxt.eps}
\caption{\sl The light $\cp$-even Higgs boson mass is given as a function of 
$x_{\rm top}$ for $\tb = 3, 50$, $m_A = 1000$ GeV, a common soft SUSY 
breaking scale for the squarks, $\msq = 1000$ GeV, and a gluino 
mass $\mgl = 800$ GeV.
}
\label{fig:mhvsxt}
\end{figure}


\begin{figure} 
\hspace{-0.5cm}
\epsfxsize=0.55\textwidth 
\epsffile{atvsms.eps}
\epsfxsize=0.55\textwidth 
\epsffile{tbvsxt.eps}
\caption{\sl Scatter plots of $A_{\rm top}$ vs. $m_S$, the
mass scale appearing in the denominator of Eq.~(\ref{eq:stmix}) (a)
and $\tb$ vs. $x_{\rm top}$ (b).
}
\label{fig:atvsmsxtvstb}
\end{figure}

In Fig.~\ref{fig:mhvsxt}, we show the dependence of the lightest Higgs 
boson mass on the stop mixing parameter $x_{\rm top}$ defined by 

\be 
x_{\rm top} \equiv \frac{A_{\rm top} - \mu/\tb}{m_S}, \; \; \; \;
{\rm where} \; \; \; \; m_S = {\sqrt{(m_{\stopu}^2 + m_{\stopd}^2)/2}}.
\label{eq:stmix}
\ee

For equal soft SUSY breaking parameters in the stop sector with
the $D$-terms neglected, $x_{\rm top}$ corresponds to the ratio 
$X_t/M_S$ of the off-diagonal and diagonal entries in the stop mixing
matrix, see e.g. Ref.~\cite{bse}.

Maximal $\mh$ values are obtained for $x_{\rm top} \approx \pm 2$, a
minimum is reached around $x_{\rm top} \approx 0$. Thus, for large
$\mh$ values a large numerator in Eq.~(\ref{eq:stmix}) is required.
From Fig.~\ref{fig:atvsmsxtvstb}(a), one can see that in GMSB only negative 
values of $A_{\rm top}$ are allowed at the electroweak scale, as a 
consequence of the fact that the trilinear couplings are negligible at 
the messenger scale. Due to the logarithmic dependence of $\mh$ on the stop 
masses, relatively large values of $|A_{\rm top}|$ are needed for large
$\mh$. In addition, large $\tb$ is also required. 
From Fig.~\ref{fig:atvsmsxtvstb}(b) 
one can check that this leads to values of $x_{\rm top} \approx -0.95$, 
which can only be achieved for positive $\mu$. 
Fig.~\ref{fig:at-mhvsmu}(a) shows the dependence of $A_{\rm top}$
on $\mu$. Large values of $|A_{\rm top}|$ are only reached for large 
$|\mu|$ values. Therefore maximal $h$ masses are obtained for relatively 
large and positive $\mu$, as can be seen in 
Fig.~\ref{fig:at-mhvsmu}(b).\footnote{In general, for large values of 
$|\mu|$ and $\tb$ the effects of the corrections from the $b$--$\tilde{b}$ 
sector can become important, leading to a decrease in $\mh$. For the 
GMSB models under consideration, however, this is not the case as a 
consequence of the relatively large $\tilde{b}$ masses.}

\newpage

All these arguments about the combination of low energy parameters 
needed for large $\mh$ in GMSB are summarised in Tab.~\ref{tab:one},   
where we report the 10 models in our sample that give 
rise to the highest $\mh$ values. Together with $\mh$, Tab.~\ref{tab:one} 
shows the corresponding input GMSB parameters [cfr. Eq.~(\ref{eq:pars})] 
as well as the values of the low energy parameters mentioned above.

It is interesting to note that all the models shown in Tab.~1 feature
a large messenger index and values of the messenger scale not far from 
the maximum we allowed while generating GMSB models. We could not construct
a single model with $\mh \gtap 122.5$ GeV having $N_{\rm mess} < 6$ or
$M_{\rm mess} < 10^5$ TeV, for $m_t = 175$ GeV. It is hence worth mentioning
here that our choice of imposing $M_{\rm mess}/\Lambda < 10^5 \Rightarrow 
M_{\rm mess} \ltap 2 \times 10^{10}$ GeV does not correspond to any solid 
theoretical prejudice. On the other hand it is true that 
$M_{\rm mess} \gtap 3 \times 10^8$ GeV always corresponds to gravitino masses 
larger than $\sim 1$ keV, due to Eqs.~(\ref{eq:Gmass}) and (\ref{eq:sqrtFmin}).
The latter circumstance might be disfavoured by cosmological 
arguments~\cite{cosmo}. A curious consequence is that the GMSB models with
the highest $\mh$ belong always to the stau NLSP or slepton co-NLSP scenarios.

Note also that restricting ourselves to GMSB models with 
$\Lambda < 100$ TeV, $M_{\rm mess} < 10^5$ TeV and $N_{\rm mess} \le 4$,
we find a maximal $\mh$ value of 122.2 GeV, for $\mt = 175$ GeV and 
$\tb \sim 52$. This is to be compared with the one-loop result of 
Ref.~\cite{kaeding}, $\mh$(max) = 131.7, for $\tb$ around 30 (the
assumed value of $\mt$ is not quoted).

Values for $\mh$ slightly larger than those we found here may also arise 
from non-minimal contributions to the Higgs potential, 
in connection with a dynamical generation of $\mu$ and $B\mu$ 
(cfr. Sec.~\ref{sec:models}). 
A treatment of this problem can be found in Ref.~\cite{AHW}.

\newpage

One should also keep in mind that our analysis still suffers from 
uncertainties due to unknown higher order corrections both in the RGE's
for GMSB model generation and in the evaluation of $\mh$ from low energy 
parameters. A rough estimate of these effects leads to shifts in $\mh$
not larger than 3 to 5 GeV.

\section*{A.4~~Conclusions}
\label{sec:A4}

\noindent
We conclude that in the minimal GMSB framework described above, values 
of $\mh \gtap 124.2$ GeV are not allowed for $\mt = 175$ GeV. This is 
almost 6 GeV smaller than the maximum value for $\mh$ one can achieve 
in the MSSM without any constraints or assumptions about the structure 
of the theory at high energy scales~\cite{mhiggslong,tbexcl,bench}.
On the other hand, the alternative mSUGRA framework allows values of 
$\mh$ that are $\sim 3$ GeV larger than in GMSB~\cite{mhiggsmSUGRA}. 
This makes the GMSB scenario slightly easier to explore via Higgs boson 
search. This result was expected in the light of the rather strong GMSB 
requirements, such as the presence of a unique soft SUSY breaking scale, 
the relative heaviness of the squarks and the gluino compared to 
non-strongly interacting sparticles, and the fact that the soft SUSY
breaking trilinear couplings $A_f$ get nonzero values at the electroweak 
scale only by RGE evolution.
Nevertheless, once the whole parameter space is explored, it is not true 
that mGMSB gives rise to $\mh$ values that are considerably smaller than in 
mSUGRA. Even smaller differences in the maximal $\mh$ might be present 
when considering non-minimal, complex messenger sectors~\cite{kaeding}
or additional contributions to the Higgs potential~\cite{GMSBmodels1,AHW}.
In any case, as for mSUGRA, current LEP~II or Tevatron data on Higgs boson 
searches are far from excluding mGMSB, and the upgraded Tevatron and 
the LHC will certainly be needed to deeply test any realistic SUSY model. 

\newpage 

\begin{figure}[ht]
\hspace{-0.5cm}
\epsfxsize=0.55\textwidth 
\epsffile{atvsmu.eps}
\epsfxsize=0.55\textwidth 
\epsffile{mhvsmu.eps}
\caption{\sl Scatter plots for $A_{\rm top}$ (a) and 
the light scalar Higgs mass (b) vs. the SUSY higgsino 
mass $\mu$ evaluated at the electroweak scale.
}
\label{fig:at-mhvsmu}
\end{figure}

\begin{table}[ht]
\begin{tabular}{|r||c|c|c|c|c|c|c|c|c|c|c|}  \hline
ID & $\mh$ & $N_{\rm mess}$ & $M_{\rm mess}$    & $\Lambda$ &
$\mu$ & $\tb$ & $m_A$ & $m_{\stopu}$ 
& $m_{\stopd}$ & $A_{\rm top}$ & $x_{\rm top}$ \\
& GeV &        & $10^6$ TeV    & TeV  & GeV &             
& GeV & GeV  & GeV    & GeV    &  \\ \hline \hline
A1    & 124.2 & 7 & 1.00 & 72.7 & 1470 & 53.4 & 367 & 2320
      & 2510  & -2150 & -0.90 \\ \hline 

A2    & 124.2 & 8 & 4.48 & 66.7 & 1570 & 53.1 & 436 & 2400
      & 2600  & -2310 & -0.94 \\ \hline

A3    & 123.7 & 6 & 2.07 & 87.0 & 1580 & 52.8 & 485 & 2420
      & 2630  & -2240 & -0.90 \\ \hline 

A4    & 123.7 & 8 & 4.67 & 52.4 & 1270 & 52.9 & 373 & 1930 
      & 2100  & -1850 & -0.93 \\ \hline

A5    & 123.5 & 8 & 4.89 & 51.1 & 1250 & 52.7 & 388 & 1880
      & 2050  & -1810 & -0.93 \\ \hline

A6    & 123.5 & 6 & 2.54 & 67.1 & 1260 & 53.0 & 349 & 1910 
      & 2080  & -1760 & -0.89 \\ \hline

A7    & 123.4 & 8 & 4.62 & 61.6 & 1470 & 51.9 & 549 & 2230
      & 2430  & -2160 & -0.94 \\ \hline

A8    & 123.4 & 6 & 4.15 & 88.1 & 1630 & 51.9 & 609 & 2450  
      & 2670  & -2300 & -0.91 \\ \hline
 
A9    & 123.3 & 7 & 3.77 & 70.3 & 1490 & 51.8 & 567 & 2260
      & 2460  & -2150 & -0.92 \\ \hline

A10   & 123.3 & 8 & 3.74 & 72.1 & 1677 & 50.7 & 756 & 2580 
      & 2800  & -2500 & -0.94 \\ \hline
\end{tabular}
\caption{\sl The 10 GMSB models giving rise to the highest $\mh$ values
in our sample. For each model, together with the light Higgs mass, we show 
the values of the GMSB input parameters and other low energy parameters of
interest for calculating $\mh$.}
\label{tab:one}
\end{table}

\bigskip
\subsection*{Acknowledgements}
S.~A. and S.~H. thank the organisers of the Workshop ``Physics at TeV
Colliders'', for the hospitality and the pleasant and productive atmosphere
in Les Houches. 

\newpage

\hrule\hfill

\vspace{0.5cm} 

\noindent 
{\Large {\bf (B) Measuring the SUSY Breaking Scale at the LHC}} \\
{\Large {\bf \null~~~~~~in the Slepton NLSP Scenario of GMSB Models}} \\

\vspace{0.2cm} 

\hrule\hfill

\vspace{0.2cm}

\noindent
{\large Contribution by:} \\
{\large S.~Ambrosanio, B.~Mele, S.~Petrarca, G.~Polesello, A.~Rimoldi}

\vspace{0.5cm}

\section*{B.1~~Introduction}
\label{sec:B1}

\noindent
The fundamental scale of SUSY breaking $\sqrt{F}$ is perhaps the most 
important quantity to determine from phenomenology in a SUSY theory.  
In the mSUGRA framework, the gravitino mass sets the scale of the 
soft SUSY breaking masses in the MSSM ($\sim 0.1-1$ TeV), so that $\sqrt{F}$ 
is typically large $\sim 10^{10-11}$ GeV [cfr. Eq.~(\ref{eq:Gmass})]. 
As a consequence, the interactions of the $\G$ with the other MSSM 
particles $\sim F^{-1}$ are too weak for the gravitino to be of relevance 
in collider physics and there is no direct way to access $\sqrt{F}$ 
experimentally. In GMSB theories, the situation is completely different. 
The soft SUSY breaking scale of the MSSM and the sparticle masses are set by
gauge interactions between the messenger and low energy sectors  
to be $\sim \alpha_{\rm SM}\Lambda$ [cfr. Eq.~(\ref{eq:bound})], so that 
typical $\Lambda$ values are $\sim 10-100$ TeV. On the other hand, 
$\sqrt{F}$ is subject to the lower bound (\ref{eq:sqrtFmin}) only, 
which tells us that values well below $10^{10}$ GeV and even as low as 
several tens of TeV are perfectly reasonable.
The $\G$ is in this case the LSP and its interactions are strong enough to
allow NLSP decays to the $\G$ inside a typical detector size.
The latter circumstance gives us a chance for extracting $\sqrt{F}$ 
experimentally through a measurement of the NLSP mass and lifetime
[cfr. Eq.~(\ref{eq:NLSPtau})]. 

Furthermore, the possibility of determining $\sqrt{F}$ with good 
precision opens a window on the physics of the SUSY breaking (the so-called 
``secluded'') sector and the way this SUSY breaking is transmitted to the 
messenger sector. Indeed, the characteristic scale of SUSY breaking felt 
by the messengers (and hence the MSSM sector) given by 
$\sqrt{F_{\rm mess}}$ in Eq.~(\ref{eq:sqrtFmin}) can be also determined
once the MSSM spectrum is known. By comparing the measured values 
of $\sqrt{F}$ and $\sqrt{F_{\rm mess}}$ it might well be possible to
get information on the way the secluded and messenger sector communicate
to each other. For instance, if it turns out that $\sqrt{F_{\rm mess}}
\ll \sqrt{F}$, then it is very likely that the communication occurs
radiatively and the ratio $\sqrt{F_{\rm mess}/F}$ is given by some loop 
factor. On the contrary, if the communication occurs via a direct
interaction, this ratio is just given by a Yukawa-type coupling 
constant, with values $\ltap 1$, see Refs.~\cite{GR-GMSB,AKM-LEP2}. 

An experimental method to determine $\sqrt{F}$ at a TeV scale $\epem$
collider through the measurement of the NLSP mass and lifetime was 
presented in Ref.~\cite{AB-LC}, in the neutralino NLSP scenario. 
Here, we are concerned about the same problem, but at a hadron collider,
the LHC, and in the stau NLSP or slepton co-NLSP scenarios. 
These scenarios provide a great opportunity at the LHC, since 
the characteristic signatures with semi-stable charged tracks 
are muon-like, but come from massive sleptons with a $\beta$ 
significantly smaller than 1. In particular, we perform
our simulations in the ATLAS muon detector, whose large size and 
excellent time resolution~\cite{TDR} allow a precision measurement of 
the slepton time of flight from the production vertex out to the muon 
chambers and hence of the slepton velocity. 
Moreover, in the stau NLSP or slepton co-NLSP scenarios, the 
knowledge of the NLSP mass and lifetime is sufficient to 
determine $\sqrt{F}$, since the factor ${\cal B}$ in 
Eq.~(\ref{eq:NLSPtau}) is exactly equal to 1. This is not the
case in the neutralino NLSP scenario, where ${\cal B}$ depends 
at least on the neutralino physical composition, and more information 
and measurements are needed for extracting a precise value of $\sqrt{F}$.

\section*{B.2~~Choice of the Sample Models \\ \null~~~~~~~and Event Simulation}
\label{sec:B2}

\noindent 
The two main parameters affecting the experimental measurement at the
LHC of the slepton NLSP properties are the slepton mass and momentum 
distribution. Indeed, at a hadron collider most of the NLSP's come
from squark and gluino production, followed by cascade decays. 
Thus, the momentum distribution is in general a function of the whole 
MSSM spectrum. However, one can approximately assume that most of the 
information on the NLSP momentum distribution is provided by the squark 
mass scale $m_{\tilde q}$ only (in the stau NLSP scenario or slepton 
co-NLSP scenarios of GMSB, one generally finds 
$m_{\tilde{g}} \gtap m_{\tilde q}$). 
To perform detailed simulations, we select a representative set of 
GMSB models generated by {\tt SUSYFIRE}. We limit ourselves to models 
with $m_{\rm NLSP} > 100$ GeV, motivated by the discussion in 
Sec.~\ref{sec:models}, and $m_{\tilde q} < 2$ TeV, in order to yield 
an adequate event statistics after a three-year low-luminosity run 
(corresponding to 30 fb$^{-1}$) at the LHC.
Within these ranges, we choose eight extreme points (four in the stau 
NLSP scenario and four in the slepton co-NLSP scenario) allowed by GMSB
in the ($m_{\rm NLSP}$, $m_{\tilde q}$) plane, in order to cover the 
various possibilities. 

\begin{table}
\begin{center}
\begin{tabular}{|r||r|r|r|r|c|} \hline
ID & $M_{\rm mess}$ (TeV) & $N_{\rm mess}$ 
   & $\Lambda$ (TeV)  & $\tb$ & sign($\mu$)  \\ \hline\hline
B1 & 1.79$\times 10^4$ & 3 &  26.6  &  7.22  & --   \\ \hline
B2 & 5.28$\times 10^4$ & 3 &  26.0  &  2.28  & --   \\ \hline
B3 & 4.36$\times 10^2$ & 5 &  41.9  & 53.7~  & +    \\ \hline
B4 & 1.51$\times 10^2$ & 4 &  28.3  &  1.27  & --   \\ \hline
B5 & 3.88$\times 10^4$ & 6 &  58.6  & 41.9~  & +    \\ \hline
B6 & 2.31$\times 10^5$ & 3 &  65.2  &  1.83  & --   \\ \hline
B7 & 7.57$\times 10^5$ & 3 & 104~~  &  8.54  & --   \\ \hline
B8 & 4.79$\times 10^2$ & 5 &  71.9  &  3.27  & --   \\ \hline
\end{tabular}
\caption{\sl Input parameters of the sample GMSB models chosen for our
study.}
\label{tab:tbg}
\end{center}
\end{table}

In Tab.~\ref{tab:tbg}, we list the input GMSB parameters we used to generate
these eight points, while in Tab.~\ref{tab:tbg1} we report the corresponding 
values of the stau mass, the squark mass scale and the gluino mass. 
The ``NLSP'' column indicates whether the model belongs to the stau NLSP 
or slepton co-NLSP scenario. The last column gives the total cross section 
in pb for producing any pairs of SUSY particles at the LHC. 

\begin{table}
\begin{center}
\begin{tabular}{|r||r|c|r|r|c|} \hline
ID &  $m_{\tilde\tau_1}$ (GeV) 
   & ``NLSP'' & $m_{\tilde q}$ (GeV) & $m_{\tilde g}$ (GeV) 
   & $\sigma$ (pb) \\ \hline \hline
B1 &  100.1 & $\tilde\tau$  &  577 &  631 & 42~~~~~ \\ \hline
B2 &  100.4 & $\tilde\ell$  &  563 &  617 & 50~~~~~ \\ \hline
B3 &  101.0 & $\tilde\tau$  & 1190 & 1480 & ~0.59~  \\ \hline
B4 &  103.4 & $\tilde\ell$  &  721 &  859 & 10~~~~~ \\ \hline
B5 &  251.2 & $\tilde\tau$  & 1910 & 2370 & ~0.023  \\ \hline
B6 &  245.3 & $\tilde\ell$  & 1290 & 1410 & ~0.36~  \\ \hline
B7 &  399.2 & $\tilde\tau$  & 2000 & 2170 & ~0.017  \\ \hline
B8 &  302.9 & $\tilde\ell$  & 1960 & 2430 & ~0.022  \\ \hline 
\end{tabular}
\caption{\sl Features of the sample GMSB model points studied.}
\label{tab:tbg1}
\end{center}
\end{table}

For each model, the events were generated with the {\tt ISAJET} 
Monte Carlo \cite{isajet} that incorporates the calculation of 
the SUSY mass spectrum and branching fraction 
using the GMSB parameters as input. We have checked that for 
the eight model points considered the sparticle masses 
calculated with {\tt ISAJET} are in good agreement with the output  
of {\tt SUSYFIRE}.

The generated events were then passed through {\tt ATLFAST} \cite{ATLFAST}, 
a fast particle-level simulation of the ATLAS detector. The {\tt ATLFAST}
package, however, was only used to evaluate the efficiency of the calorimetric
trigger that selects the GMSB events. The detailed response of the detector 
to the slepton NLSP has been parametrised for this work using the results  
of a full simulation study, as described in the next section.

\section*{B.3~~Slepton detection}
\label{sec:B3} 

The experimental signatures of heavy long-lived charged 
particles at a hadron collider have already been studied both 
in the framework of GMSB and in more general scenarios 
\cite{leandro,drtata, femoroi,marthom}.
The two main observables one can use to separate
these particles from muons are the high specific ionisation
and the time of flight in the detector.

We concentrate here on the measurement of the time of flight, 
made possible by the timing precision ($\ltap 1$~ns) 
and the size of the ATLAS muon spectrometer.

It was demonstrated with a full simulation of the ATLAS muon 
detector~\cite{gpar} that the $\beta$ of a particle can be measured 
with a resolution that can be approximately parameterised as 
$\sigma(\beta)/\beta^2=0.028$.
The resolution on the transverse momentum measurement for heavy
particles is found to be comparable to the one expected for muons.
We have therefore simulated the detector response to NLSP sleptons 
by smearing the slepton momentum and $\beta$ according 
to the parameterisations in Ref.~\cite{gpar}. 

An important issue is the online selection of the SUSY events.
We have not made any attempt to evaluate whether the heavy sleptons 
can be selected using the muon trigger. For the event selection, we 
rely on the calorimetric $E_{T}^{\rm miss}$ trigger, 
consisting in the requirement of at least a hadronic jet with $p_T>50$~GeV,
and a transverse momentum imbalance calculated only from the energy 
deposition in the calorimeter larger than 50~GeV. We checked that 
this trigger has an efficiency in excess of 80\% for all the considered
models. 

A detailed discussion of the experimental assumptions underlying
the results presented here is given in Ref.~\cite{AMPPR}.

\section*{B.4~~Event Selection and Slepton Mass Measurement} 
\label{sec:B4} 

\noindent
In order to select a clean sample of sleptons, we apply the
following requirements: 

\begin{itemize}
\item
at least a hadronic jet with $P_T>50$~GeV and a calorimetric \\
\mbox{$E_T^{\rm miss}>50$~GeV} (trigger requirement);
\item
at least one candidate slepton satisfying the following cuts: 
\begin{itemize}
\item
$|\eta|<$2.4 to ensure that the particle is in the
acceptance of the muon trigger chamber, and therefore 
both coordinates can be measured;
\item
$\beta_{\rm meas}<0.91$, where $\beta_{\rm meas}$ is the $\beta$ of the
particle measured with the time of flight in the precision chambers;
\item
The $P_T$ of the slepton candidate, after the energy loss in the
calorimeters has been taken into account, must be larger than
10~GeV, to ensure that the particle traverse all of the muon stations.
\end{itemize}
\end{itemize} 

Considering an integrated luminosity of 30~fb$^{-1}$, 
a number of events ranging from a few hundred for the models with
2~TeV squark mass scale to a few hundred thousand for a 500~GeV
mass scale survive these cuts and can be used for measuring the 
NLSP properties.

From the measurements of the slepton momentum and of particle $\beta$,
the mass can be determined using the standard relation 
$m=p\frac{\sqrt{1-\beta^2}}{\beta}$. 
For each value of $\beta$ and momentum, the measurement error is
known and it is given by the parametrisations in Ref.~\cite{gpar}.
Therefore, the most straightforward way of measuring the mass 
is just to use the weighted average of all the masses 
calculated with the above formula. 

In order to perform this calculation, the particle momentum is needed,
which implies measuring the $\eta$ coordinate. In fact, with the precision 
chambers only one can only measure the momentum components  
transverse to the beam axis. 

The measurement of the second coordinate must be provided by the 
trigger chambers, for which only a limited time window around the beam 
crossing is read out, therefore restricting the $\beta$ range where 
this measurement is available. Hence, we have evaluated the achieved 
measurement precision for two different $\beta$ intervals: $0.6<\beta<0.91$ 
and $0.8<\beta<0.91$ for the eight sample points.  
We found a statistical error well below the 0.1\% level for those model
points having $m_{\tilde q} < 1300$ GeV. Even for the three models 
(B5, B7, B8) with lower statistics ($m_{\tilde q} \simeq 2$ TeV), the 
error stays below the 0.4\% level. 

Many more details, tables and figures about this part of our study 
can be found in Ref.~\cite{AMPPR}. 

\section*{B.5~~Slepton Lifetime Measurement}
\label{sec:B5}

\noindent 
The measurement of the NLSP lifetime at a high energy $e^+e^-$ collider
was studied in detail in Ref.~\cite{AB-LC} for the neutralino NLSP case.
Similar to that study, the measurement of the slepton NLSP lifetime we are 
interested in here can be performed by exploiting the fact that two NLSP's 
are produced in each event. 
One can therefore select $N_1$ events where a slepton is detected through 
the time-of-flight measurement described above, count the number of 
times $N_2$ when a second slepton is observed and use this information  
to measure the lifetime. 
Although in principle very simple, in practice this method requires 
an excellent control on all possible sources of inefficiency for 
detecting the second slepton.

\begin{figure} 
\hspace{3.0cm}
\epsfxsize=0.7\textwidth
\epsffile{ratiovsctau.eps}
\caption{\sl The ratio $R=N_2/N_1$ defined in the text as a function 
of the slepton lifetime $c\tau$. Only the curves corresponding to 
the model points B1, B5, B6, B8 are shown.} 
\label{fig:ratio}
\end{figure} 

We give here the basis of the method, without mentioning the experimental 
details. We provide an estimate of the achievable statistical error for 
the models considered and a parametrisation of the effect on the lifetime
measurement of a generic systematic uncertainty on the slepton efficiency.
In case the sparticle spectrum and BR's can be measured from the SUSY events, 
as e.g. shown in Ref.~\cite{ihfp}, an accurate simulation of all the SUSY 
production processes can be performed, and the results from this section are 
representative of the measurement precision achievable in a real experiment.

Another method based on the same principles, but assuming minimal knowledge 
of the SUSY spectrum, is described in Ref.~\cite{AMPPR}, where 
a detailed estimate of the achievable systematic precision is given.

We define $N_1$ starting from the event sample 
defined by the cuts discussed in Sec.~\ref{sec:B4}, 
with the additional requirement that, 
for a given value of the slepton lifetime, 
at least one of the produced sleptons decays at 
at a distance from the interaction vertex $>10$~m,
and is therefore reconstructed in the muon system.
For the events thus selected, we define $N_2$ as the 
subsample where a second particle with a transverse momentum 
$>10$~GeV is identified in the muon system. The search for the
second particle should be as inclusive as possible, in order to minimise
the corrections to the ratio. In particular, the cut $\beta_{\rm meas}<0.91$ 
is not applied, but  particles with a mass measured from $\beta$ and momentum 
incompatible with the measured slepton mass are rejected. This leaves
a background of high momentum muons in the sample that can be statistically
subtracted using the momentum distribution of electrons. 
The ratio

\be
R=\frac{N_2}{N_1}
\label{eq:ratio}
\ee

\noindent
is a function of the slepton lifetime. Its dependence on the
NLSP lifetime $c\tau$ in metres in shown in Fig.~\ref{fig:ratio}
for four among our eight sample models. The curves for the model
points not shown are either very similar to one of the curves we 
show or are mostly included between the external curves corresponding to 
points B1 and B8, thus providing no essential additional information.
Note that the curve for model 6 starts from $c\tau = 2.5$~m and not
from $c\tau = 50$~cm, as for the other models. This is due to the 
large value of $M_{\rm mess}$ (cfr. Tab.~\ref{tab:tbg}), determining
a minimum NLSP lifetime allowed by theory which is macroscopic in 
this case [cfr. Eqs.~(\ref{eq:NLSPtau}) and (\ref{eq:sqrtFmin})].

\begin{figure}[ht] 
\epsfxsize = \textwidth 
\epsffile{sigctau1.eps}
\caption{\sl Fractional error on the measurement of the 
slepton lifetime $c\tau$, for model sample points B1 to B4.
We assume an integrated luminosity of 30~fb$^{-1}$. 
The curves are shown for three different assumptions on the 
fractional systematic error on the $R$ measurement: 
statistical error only (full line), 1\% systematic error (dashed line),
5\% systematic error (dotted line).
}
\label{fig:sigctau1}
\end{figure}

The probability for a particle of mass $m$, momentum $p$ and proper 
lifetime $\tau$ to travel a distance $L$ before decaying is
given by the expression

\be
P(L)=e^{-mL/pc\tau}.
\label{eq:prob} 
\ee

$N_2$ is therefore a function of the momentum distribution
of the slepton, which is determined by the details of the SUSY spectrum.
One needs therefore to be able to simulate the full SUSY cascade decays
in order to construct the $c\tau$--$R$ relationship. 

\begin{figure} 
\epsfxsize = \textwidth
\epsffile{sigctau2.eps}
\caption{\sl The same as in Fig.~\ref{fig:sigctau2}, but for 
the model sample points B5 to B8.
}
\label{fig:sigctau2}
\end{figure}

The statistical error on $R$ can be evaluated as

\be
\sigma(R)=\sqrt{\frac{R(1-R)}{N_1}}.
\label{eq:errat} 
\ee

Relevant for the precision with which the SUSY breaking scale can be 
measured is instead the error on the measured $c\tau$. This can be 
extracted from the curves shown in Fig.~\ref{fig:ratio} and can be 
evaluated as

\be
\sigma({\mathrm c}\tau)=\sigma(R)/
\left[\frac{\partial R(c\tau)}{\partial c\tau}\right].
\label{eq:erctau}
\ee

The measurement precision calculated according to this formula is shown 
in Figs.~\ref{fig:sigctau1} and \ref{fig:sigctau2} for the eight sample 
points, for an integrated luminosity of 30~fb$^{-1}$. 
The full line in the plots is the error on $c\tau$ considering the statistical
error on $R$ only. The available statistics is a function 
of the strongly interacting sparticles' mass scale. 
Even if a precise $R$--$c\tau$ relation can be built from the 
knowledge of the model details, there will be a systematic uncertainty 
in the evaluation of the losses in $N_2$, because of sleptons produced outside
the $\eta$ acceptance, or absorbed in the calorimeters, or escaping the 
calorimeter with a transverse momentum below the cuts. 
The full study of these uncertainties is in progress. 

\begin{figure}[ht]
\epsfxsize = \textwidth
\epsffile{sigsqrtf1.eps}
\caption{\sl Fractional error on the measurement of the 
SUSY breaking scale $\sqrt{F}$ for model sample points B1 to B4.
We assume an integrated luminosity of 30~fb$^{-1}$. 
The curves are shown for the three different assumptions on the fractional 
systematic error used in Figs.~\ref{fig:sigctau1} and \ref{fig:sigctau2}.
}
\label{fig:sigsqrtf1}
\end{figure}

At this level, we just parameterise the systematic error as a term 
proportional to $R$, added in quadrature to the statistical error. 
We choose two values, $1\% R$ and $5\% R$, and propagate the error to 
the $c\tau$ measurement. 
The results are represented by the dashed and dotted lines in
Figs.~\ref{fig:sigctau1} and \ref{fig:sigctau2}. 

For the models with squark mass scales up to 1200~GeV,
assuming a 1\% systematic error on the measured ratio, 
a precision better than 10\% on the $c\tau$ measurement 
can be obtained for lifetimes between 0.5--1~m and 50--80~m. 
If the systematic uncertainty grows to 5\%, the 10\% precision
can only be achieved in the range 1--10~m. If the mass scale goes up to 2~TeV,
even considering a pure statistical error only, a 10\% precision
is not achievable. However a 20\% precision is possible  
over $c\tau$ ranges between 5 and 100~m, assuming a 1\% systematic error.

\begin{figure}
\epsfxsize = \textwidth
\epsffile{sigsqrtf2.eps}
\caption{\sl The same as in Fig.~\ref{fig:sigsqrtf1}, but for model
sample points B5 to B8.
}
\label{fig:sigsqrtf2}
\end{figure}

Note that the curves corresponding to the model points B2, B6 and B7 
do not start from $c\tau = 50$ cm, but from the theoretical lower 
limit on $c\tau$ of 1.8, 2.5 and 6.1 metres, respectively.

\section*{B.6~~Determining the SUSY Breaking Scale $\sqrt{F}$}
\label{sec:B6}

\noindent 
Using the measured values of $c\tau$ and the NLSP mass,
the SUSY breaking scale $\sqrt{F}$ can be calculated from 
Eq.~(\ref{eq:NLSPtau}), where ${\cal B} = 1$ for the case where 
the NLSP is a slepton.  
From simple error propagation, the fractional uncertainty on the 
$\sqrt{F}$ measurement can be obtained adding in quadrature one 
fourth of the fractional error in $c\tau$ and five fourths of the 
fractional error on the slepton mass. 

In Figs.~\ref{fig:sigsqrtf1} and \ref{fig:sigsqrtf2}, we show the fractional 
error on the $\sqrt{F}$ measurement as a function of $\sqrt{F}$ for 
our three different assumptions on the $c\tau$ error. 
The uncertainty is dominated by $c\tau$ for the higher part of the 
$\sqrt{F}$ range and grows quickly when approaching
the lower limit on $\sqrt{F}$. This is because very few sleptons 
survive and the statistical error on both $m_{\tilde \ell}$ and 
$c\tau$ gets very large. 
If we assume a 1\% systematic error on the ratio $R$ from which $c\tau$
is measured (dashed lines in Figs.~\ref{fig:sigsqrtf1} and 
\ref{fig:sigsqrtf2}), the error on $\sqrt{F}$ is better than 10\% 
for $1000 \ltap \sqrt{F} \ltap 4000$~TeV for model points B1--B4 with
higher statistics. For points B5--B8, in general one can explore 
a range of higher $\sqrt{F}$ values with a small relative error, 
essentially due to the heaviness of the decaying NLSP in these models.
Note also that the theoretical lower limit (\ref{eq:sqrtFmin}) on $\sqrt{F}$ 
is equal to about 1200, 1500, 3900, 8900 TeV respectively in model 
points B2, B5, B6, B7, while it stays well below 1000 TeV for the other models.

\section*{B.7~~Conclusions}
\label{sec:B7}

\noindent 
We have discussed a simple method to measure at the LHC with the ATLAS 
detector the fundamental SUSY breaking scale $\sqrt{F}$ in the GMSB scenarios 
where a slepton is the NLSP and decays to the gravitino with a lifetime
in the range 0.5~m $\ltap c\tau_{\rm NLSP} \ltap 1$~km. This method requires
the measurement of the time of flight of long lived sleptons 
and is based on counting events with one or two identified NLSP's.  
It relies on the assumptions that a good knowledge of the MSSM sparticle 
spectrum and BR's can be extracted from the observation of the SUSY events 
and that the systematic error in evaluating the slepton losses can be kept
below the few percent level.
We performed detailed, particle level simulations for eight representative 
GMSB models, some of them being particularly hard due to low statistics. 
We found that a level of precision of a few 10's \% on the SUSY breaking 
scale measurement can be achieved in significant parts of the 
$1000 \ltap \sqrt{F} \ltap 30000$~TeV range, for all models considered. 
More details as well as a full study of the systematics associated with this 
procedure and another less ``model-dependent'' method to measure $\sqrt{F}$ 
is presented in detail in Ref.~\cite{AMPPR}. 

\bigskip
\subsection*{Acknowledgements}
S.~A. and G.~P. thank the organisers of the Workshop ``Physics at TeV
Colliders'', for the hospitality and the pleasant and productive atmosphere
in Les Houches.  

\vspace{0.7cm}

\begin{thebibliography}{10}

\bibitem{StevePrimer}
  For a recent pedagogical review of supersymmetry and supersymmetry 
  breaking, see S.~P.~Martin, ``A Supersymmetry Primer'', in ``Perspectives
  on Supersymmetry'', G.~L.~Kane ed., World Scientific 1998,  
  and references therein. 

\bibitem{oldGMSB}
    M.~Dine, W.~Fischler, M.~Srednicki, \NPB{189}{1981}{575};
    S.~Dimopoulos, S.~Raby, \NPB{192}{1981}{353};
    M.~Dine, W.~Fischler, \PLBold{110}{1982}{227};
    M.~Dine,  M.~Srednicki, \NPB{202}{1982}{238};
    M.~Dine, W.~Fischler, \NPB{204}{1982}{346};
    L.~Alvarez-Gaum\'e, M.~Claudson, M.~B.~Wise, \NPB{207}{1982}{96};
    C.~R.~Nappi, B.~A.~Ovrut, \PLBold{113}{1982}{175};
    S.~Dimopoulos, S.~Raby, \NPB{219}{1983}{479}.

\bibitem{newGMSB}
    M.~Dine, A.~E.~Nelson, \PRD{48}{1993}{1277};
    M.~Dine, A.~E.~Nelson, Y.~Shirman, \PRD{51}{1995}{1362};
    M.~Dine, A.~E.~Nelson, Y.~Nir, Y.~Shirman, \PRD{53}{1996}{2658}.

\bibitem{GR-GMSB}
  For a recent review, see G.~F.~Giudice, R.~Rattazzi, ``Theories with 
  Gauge-Mediated Supersymmetry Breaking'',  submitted to 
  Phys. Rep.

\bibitem{GMSBmodels1}
    S.~Dimopoulos, S.~Thomas, J.~D.~Wells, \PRD{54}{1996}{3283};
    \NPB{488}{1997}{39}.

\bibitem{GMSBmodels2}
    J.~A.~Bagger, K.~Matchev, D.~M.~Pierce, R.~Zhang, 
    \PRD{55}{1997}{3188}. 

\bibitem{AKM-LEP2} 
    S.~Ambrosanio, G.~D.~Kribs, S.~P.~Martin, \PRD{56}{1997}{1761}. 

\bibitem{AB-LC}
    S.~Ambrosanio, G.~A.~Blair, \EPJC{12}{2000}{287--321}.

\bibitem{Fayet}
    P.~Fayet, \PLBold{70}{1977}{461}; \PLBold{86}{1979}{272};
    \PLB{175}{1986}{471} and in ``Unification of the fundamental 
    particle interactions", eds.~S.~Ferrara, J.~Ellis,   
    P.~van Nieuwenhuizen (Plenum, New York, 1980) p.~587.

\bibitem{AKM2} 
  A relevant example is discussed in:
  S.~Ambrosanio, G.~D.~Kribs, S.~P.~Martin, \NPB{516}{1998}{55}.

\bibitem{SUSYFIRE}
  An updated, generalised and {\tt Fortran}-linked version of the 
  program used in Ref.~\cite{AKM-LEP2}. It generates minimal and 
  non-minimal GMSB and SUGRA models.
  For inquiries about this software package, please send e-mail to  
 {\tt ambros@mail.cern.ch}. 

\bibitem{AHW}
  S.~Ambrosanio, S.~Heinemeyer, G.~Weiglein, in preparation.

\bibitem{mhiggsmSUGRA} 
  A.~Dedes, S.~Heinemeyer, P.~Teixeira-Dias, G.~Weiglein,
   to appear in Jour. Phys. G.

\bibitem{mhiggs1l} 
  H.~Haber, R.~Hempfling, \PRL{66}{1991}{1815};
  J.~Ellis, G.~Ridolfi, F.~Zwirner, \PLB{257}{1991}{83};
                                       \PLB{262}{1991}{477}.

\bibitem{mhiggsf1l} 
  P.~Chankowski, S.~Pokorski, J.~Rosiek, \NPB{423}{1994}{437}.

\bibitem{mhiggsf1ldab} 
  A.~Dabelstein, \NPB{456}{1995}{25}; \ZPC{67}{1995}{495}.

\bibitem{pierce} 
  J.~Bagger, K.~Matchev, D.~Pierce, R.~Zhang, \NPB{491}{1997}{3}.

\bibitem{mhiggsRG1}  
  J.~Casas, J.~Espinosa, M.~Quir\'os, A.~Riotto, \NPB{436}{1995}{3}, 
  E: ibid {\bf B439} (1995) 466;
  M.~Carena, J.~Espinosa, M.~Quir\'os, C.~Wagner, \PLB{355}{1995}{209};
  M.~Carena, M.~Quir\'os, C.~Wagner, \NPB{461}{1996}{407}.

\bibitem{mhiggsRG2} 
  H.~Haber, R.~Hempfling, A.~Hoang, \ZPC{75}{1997}{539}.

\bibitem{mhiggsEP} 
  R.~Hempfling, A.~Hoang, \PLB{331}{1994}{99};
  R.-J.~Zhang, \PLB{447}{1999}{89}.
%  J.~Espinosa, R.-J.~Zhang, .

\bibitem{mhiggsletter} 
  S.~Heinemeyer, W.~Hollik, G.~Weiglein, \PRD{58}{1998}{091701};
  \PLB{440}{1998}{296}.

\bibitem{mhiggslong} 
  S.~Heinemeyer, W.~Hollik, G.~Weiglein, \EPJC{9}{1999}{343}.

\bibitem{mhiggslle} 
  S.~Heinemeyer, W.~Hollik, G.~Weiglein, \PLB{455}{1999}{179}.

\bibitem{lepc} 
  P.~McNamara, ``Combined LEP Higgs search results up to 
  $\sqrt{s} = 196$ GeV'', talk given at the LEPC meeting, September 7, 1999.

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

\bibitem{m20} 
  S.~Heinemeyer, W.~Hollik, G.~Weiglein, .

\bibitem{bse} 
  M.~Carena, H.~Haber, S.~Heinemeyer, W.~Hollik, C.~Wagner, 
  G.~Weiglein, DESY 197-99, .

\bibitem{kaeding} 
  T.~A.~Kaeding, S.~Nandi, .

\bibitem{cosmo}
  See e.g.: H.~Pagels, J.~R.~Primack, \PRL{48}{1982}{223}. 

\bibitem{tbexcl} 
  S.~Heinemeyer, W.~Hollik, G.~Weiglein, .

\bibitem{bench} 
  M.~Carena, S.~Heinemeyer, C.~Wagner, G.~Weiglein, .

\bibitem{TDR}
  The ATLAS Collaboration, ``ATLAS Detector and Physics Performance Technical
  Design Report'', ATLAS TDR 15, CERN/LHCC/99-15 (1999).

\bibitem{isajet}
  H.~Baer, F.~E.~Paige, S.~D.~Protopopescu, X. Tata, ; 
  .

\bibitem{ATLFAST}
  E.~Richter-Was, D.~Froidevaux, L.~Poggioli, ``{\tt ATLFAST 2.0}: A Fast
  Simulation Package for ATLAS'', ATLAS Internal Note ATL-PHYS-98-131 (1998).

\bibitem{leandro}
  A.~Nisati, S.~Petrarca, G.~Salvini, Mod. Phys. Lett. {\bf A12} (1997) 2213.

\bibitem{drtata}
  M.~Drees, X.~Tata, \PLB{252}{1990}{695}.

\bibitem{femoroi}
  J.~L.~Feng, T.~Moroi, \PRD{58}{1998}{035001}.

\bibitem{marthom}
  S.~P.~Martin, J.~D.~Wells, \PRD{59}{1999}{035008}.

\bibitem{gpar}
  G.~Polesello, A.~Rimoldi, ``Reconstruction of Quasi-stable Charged
  Sleptons in the ATLAS Muon Spectrometer'', 
  ATLAS Internal Note ATL-MUON-99-06.

\bibitem{AMPPR}
  S.~Ambrosanio, B.~Mele, S.~Petrarca, G.~Polesello, A.~Rimodi, in preparation.

%\bibitem{leandro1}
%  A.~Nisati, ``Preliminary Timing Studies of the Barrel Muon Trigger 
%  System'', ATLAS Internal Note ATL-DAQ-98-083 (1998).

\bibitem{ihfp}
  I.~Hinchliffe, F.~E.~Paige, \PRD{60}{1999}{095002}.

\end{thebibliography}

\end{document}

%
%  SA_macros.tex:  A useful set of macros to be \input{} into .tex files.
%
%  25 Jan 98  -SA

%  Equation lines
\newcommand{ \be  }{\begin{equation}}
\newcommand{ \ee  }{\end{equation}} 
\newcommand{ \bea }{\begin{eqnarray}}
\newcommand{ \eea }{\end{eqnarray}}
\newcommand{ \beas }{\begin{eqnarray*}}
\newcommand{ \eeas }{\end{eqnarray*}}

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

%  References to main physics journals
\newcommand{ \rNPB    }[3]{Nucl. Phys. \bfseries{B#1}, #3 (#2)}
\newcommand{ \rPLB    }[3]{Phys. Lett. B \bfseries{#1}, #3 (#2)}
\newcommand{ \rPLBold }[3]{Phys. Lett. \bfseries{#1B}, #3 (#2)}
\newcommand{ \rPRD    }[3]{Phys. Rev. \bfseries{D#1}, #3 (#2)}
\newcommand{ \rPRold  }[3]{Phys. Rev. \bfseries{#1}, #3 (#2)}
\newcommand{ \rPRL    }[3]{Phys. Rev. Lett. \bfseries{#1}, #3 (#2)}
\newcommand{ \rPREP   }[3]{Phys. Rep. \bfseries{#1}, #3 (#2)}
\newcommand{ \rZPC    }[3]{Z. Phys. C \bfseries{#1}, #3 (#2)}

%\newcommand{ \NPB    }[3]{Nucl. Phys.      {\bf B#1}, #2 (#3)}
%\newcommand{ \PLB    }[3]{Phys. Lett. B    {\bf #1},  #2 (#3)}
%\newcommand{ \PLBold }[3]{Phys. Lett.      {\bf #1B}, #2 (#3)}
%\newcommand{ \PRD    }[3]{Phys. Rev. D     {\bf #1},  #2 (#3)}
%\newcommand{ \PRold  }[3]{Phys. Rev.       {\bf #1},  #2 (#3)}
%\newcommand{ \PRL    }[3]{Phys. Rev. Lett. {\bf #1},  #2 (#3)}
%\newcommand{ \PREP   }[3]{Phys. Rep.       {\bf #1},  #2 (#3)}
%\newcommand{ \ZPC    }[3]{Z. Phys. C       {\bf #1},  #2 (#3)}
%\newcommand{ \EPJC   }[3]{Eur. Phys. J. C  {\bf #1},  #2 (#3)}

\newcommand{ \NPB    }[3]{Nucl. Phys.      {\bf B#1} (#2) #3}
\newcommand{ \PLB    }[3]{Phys. Lett. B    {\bf #1}  (#2) #3}
\newcommand{ \PLBold }[3]{Phys. Lett.      {\bf #1B} (#2) #3}
\newcommand{ \PRD    }[3]{Phys. Rev. D     {\bf #1}  (#2) #3}
\newcommand{ \PRold  }[3]{Phys. Rev.       {\bf #1}  (#2) #3}
\newcommand{ \PRL    }[3]{Phys. Rev. Lett. {\bf #1}  (#2) #3}
\newcommand{ \PREP   }[3]{Phys. Rep.       {\bf #1}  (#2) #3}
\newcommand{ \ZPC    }[3]{Z. Phys. C       {\bf #1}  (#2) #3}
\newcommand{ \EPJC   }[3]{Eur. Phys. J. C  {\bf #1}  (#2) #3}


\newcommand{ \nNPB    }[3]{Nucl. Phys. B #1, #2 (#3)}
\newcommand{ \nPLB    }[3]{Phys. Lett. B #1, #2 (#3)}
\newcommand{ \nPLBold }[3]{Phys. Lett. #1B, #2 (#3)}
\newcommand{ \nPRD    }[3]{Phys. Rev. D#1, #2 (#3)}
\newcommand{ \nPRold  }[3]{Phys. Rev. #1, #2 (#3)}
\newcommand{ \nPRL    }[3]{Phys. Rev. Lett. #1, #2 (#3)}
\newcommand{ \nPREP   }[3]{Phys. Rep. #1, #2 (#3)}
\newcommand{ \nZPC    }[3]{Z. Phys. C #1, #2 (#3)}

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

%  \gsim and \lsim provide >= and <= signs.
\newcommand{ \centeron }[2]{{\setbox0=\hbox{#1}\setbox1=\hbox{#2}\ifdim
                             \wd1>\wd0\kern.5\wd1\kern-.5\wd0\fi \copy0
                             \kern-.5\wd0\kern-.5\wd1\copy1\ifdim\wd0>\wd1
                             \kern.5\wd0\kern-.5\wd1\fi}}
\newcommand{ \ltap }{\;\centeron{\raise.35ex\hbox{$<$}}
                     {\lower.65ex\hbox{$\sim$}}\;}
\newcommand{ \gtap }{\;\centeron{\raise.35ex\hbox{$>$}}
                     {\lower.65ex\hbox{$\sim$}}\;}
\newcommand{ \gsim }{\mathrel{\gtap}}
\newcommand{ \lsim }{\mathrel{\ltap}}

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

%  Slash character...
\newcommand{ \slashchar }[1]{\setbox0=\hbox{$#1$}   % set a box for #1
   \dimen0=\wd0                                     % and get its size 
   \setbox1=\hbox{/} \dimen1=\wd1                   % get size of /
   \ifdim\dimen0>\dimen1                            % #1 is bigger 
      \rlap{\hbox to \dimen0{\hfil/\hfil}}          % so center / in box 
      #1                                            % and print #1
   \else                                            % / is bigger 
      \rlap{\hbox to \dimen1{\hfil$#1$\hfil}}       % so center #1
      /                                             % and print /
   \fi}                                             %

%%EXAMPLE:  $\slashchar{E}$ or $\slashchar{E}_{t}$

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

\def\doublespaced{\baselineskip=\normalbaselineskip\multiply
                             \baselineskip by 200\divide\baselineskip by 100}
\def\singleandhalfspaced{\baselineskip=\normalbaselineskip\multiply
                             \baselineskip by 150\divide\baselineskip by 100}
\def\singlespaced{\baselineskip=\normalbaselineskip}

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

%  Various particle physics / supersymmetry commands.
\newcommand{ \ph         }{ \gamma }
\newcommand{ \textfrac   }[2]{ {\textstyle\frac{#1}{#2}} }

\newcommand{ \Ni         }{ {\tilde \chi}^{0}_i }
\newcommand{ \Nj         }{ {\tilde \chi}^{0}_j }
\newcommand{ \NI         }{ \tilde{N}_1 }
\newcommand{ \NII        }{ \tilde{N}_2 }
\newcommand{ \NIII       }{ \tilde{N}_3 }
\newcommand{ \NIIII      }{ \tilde{N}_4 }
\newcommand{ \Ci         }{ \tilde{C}_i }
\newcommand{ \Cj         }{ \tilde{C}_j }
\newcommand{ \CI         }{ \tilde{C}_1 }
\newcommand{ \CII        }{ \tilde{C}_2 }
\newcommand{ \CIp        }{ \tilde{C}^{+}_1 }
\newcommand{ \CIm        }{ \tilde{C}^{-}_1 }
\newcommand{ \CIpm       }{ \tilde{C}^{\pm}_1 }
\newcommand{ \CIIp       }{ \tilde{C}^{+}_2 }
\newcommand{ \CIIm       }{ \tilde{C}^{-}_2 }
\newcommand{ \CIIpm      }{ \tilde{C}^{\pm}_2} 
\newcommand{ \CIImp      }{ \tilde{C}^{\mp}_2} 
\def\G{ \tilde{G} }

\newcommand{ \lRp        }{ \tilde{\ell}_R^+ }
\newcommand{ \lRm        }{ \tilde{\ell}_R^- }
\newcommand{ \lRpm       }{ \tilde{\ell}_R^{\pm} }

\newcommand{ \XI         }{ {\tilde X}_1 }
\newcommand{ \XII        }{ {\tilde X}_2 }
\newcommand{ \eL         }{ {\tilde e}_{\scriptscriptstyle L} }
\newcommand{ \eR         }{ {\tilde e}_{\scriptscriptstyle R} }
\newcommand{ \eLR        }{ {\tilde e}_{\scriptscriptstyle L/R} }
\newcommand{ \veL        }{ {\tilde \nu} }
\newcommand{ \SG         }{ {\tilde \gamma} }
\newcommand{ \SZ         }{ {\tilde Z} }
%
\newcommand{ \gmu        }{ \gamma^{\mu} }
\newcommand{ \gnu        }{ \gamma^{\nu} }
\newcommand{ \gfive      }{ \gamma_{5} }
\newcommand{ \PL         }{ P_{L} }
\newcommand{ \PR         }{ P_{R} }
\newcommand{ \epsiloni   }{ \epsilon_{i} }
\newcommand{ \epsilonj   }{ \epsilon_{j} }
\newcommand{ \OLij       }{ O^{L}_{ij} }
\newcommand{ \ORij       }{ O^{R}_{ij} }
\newcommand{ \OLji       }{ O^{L}_{ji} }
\newcommand{ \ORji       }{ O^{R}_{ji} }
\newcommand{ \OLijc      }{ O^{L*}_{ij} }
\newcommand{ \ORijc      }{ O^{R*}_{ij} }
\newcommand{ \OLjic      }{ O^{L*}_{ji} }
\newcommand{ \ORjic      }{ O^{R*}_{ji} }
\newcommand{ \dL         }{ \tilde d_L }
\newcommand{ \dR         }{ \tilde d_R }
\newcommand{ \uL         }{ \tilde u_L }
\newcommand{ \uR         }{ \tilde u_R }
\newcommand{ \slepton    }{ \widetilde l }
\newcommand{ \M          }{ {\cal M} }
\newcommand{ \cO         }{ {\cal O} }
%
\newcommand{ \ra         }{ \rightarrow }
\newcommand{ \lra        }{ \longleftrightarrow }

%
\newcommand{ \sle        }{ {\tilde \ell} }
\newcommand{ \snu        }{ {\tilde \nu} }
\newcommand{ \lL         }{ { \tilde \ell}_L }
\newcommand{ \lLstar     }{ { \tilde \ell}_L^* }
\newcommand{ \lR         }{ { \tilde \ell}_R }
\newcommand{ \lRstar     }{ { \tilde \ell}_R^* }
\newcommand{ \snuL       }{ { \tilde \nu}_L }
\newcommand{ \snuLst     }{ { \tilde \nu}_L^* }
\newcommand{ \nubar      }{ \overline{ \nu } }
\newcommand{ \muL        }{ { \tilde \mu}_L }
\newcommand{ \muR        }{ { \tilde \mu}_R }
\newcommand{ \tauL       }{ { \tilde \tau}_L }
\newcommand{ \tauR       }{ { \tilde \tau}_R }
\newcommand{ \tauu       }{ { \tilde \tau}_1 }
\newcommand{ \taud       }{ { \tilde \tau}_2 }
\newcommand{ \h          }{ { h^0 } }
%
\newcommand{ \E          }{ {\slashchar{E}} }
\newcommand{ \Et         }{ { \slashchar{E}_{\scriptscriptstyle T} } }
\newcommand{ \Etcut      }{ { \Et^{\rm cut} } }
\newcommand{ \sigbreff   }{ \sigma \times {\rm BR} \times {\rm EFF} }
%
\newcommand{ \eegg       }{ {ee\gamma\gamma + \Et} }
\newcommand{ \GeV        }{ {\rm GeV} }
%
%  Obsolete
\newcommand{ \neuti      }{ {\tilde \chi}^{0}_i }
\newcommand{ \neutj      }{ {\tilde \chi}^{0}_j }
\newcommand{ \neutI      }{ {\tilde \chi}^{0}_1 }
\newcommand{ \neutII     }{ {\tilde \chi}^{0}_2 }
\newcommand{ \neutIII    }{ {\tilde \chi}^{0}_3 }
\newcommand{ \neutIIII   }{ {\tilde \chi}^{0}_4 }
\newcommand{ \chari      }{ {\tilde \chi}^{\pm}_i }
\newcommand{ \charj      }{ {\tilde \chi}^{\pm}_j }
\newcommand{ \charI      }{ {\tilde \chi}^{\pm}_1 }
\newcommand{ \charII     }{ {\tilde \chi}^{\pm}_2 }
\newcommand{ \charIplus  }{ {\tilde \chi}^{+}_1 }
\newcommand{ \charIminus }{ {\tilde \chi}^{-}_1 }

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

%  Sandro's def's

\newcommand{\susy}{{\sc SuSy}}
\newcommand{\sugra}{{\sc SuGra}} 
\newcommand{\sss}{\scriptscriptstyle}
\newcommand{\rar}{\rightarrow} 
\newcommand{\Rar}{\Rightarrow} 

\newcommand{\Z}{Z^{\sss 0}}  
\newcommand{\ep}{e^{\sss +}} 
\newcommand{\eme}{e^{\sss -}}  
\newcommand{\epm}{e^{\sss \pm}}  
\newcommand{\emp}{e^{\sss \mp}}  
\newcommand{\epem}{e^{\sss +}e^{\sss -}}   
\newcommand{\ve}{\nu_{\sss e}} 
\newcommand{\veb}{\bar{\nu}_{\sss e}}  

\newcommand{\sel}{\tilde{e}_{\sss L}}
\newcommand{\ser}{\tilde{e}_{\sss R}}
\newcommand{\selr}{\tilde{e}_{\sss L,R}}
\newcommand{\sepl}{\tilde{e}^{\sss +}_{\sss L}}
\newcommand{\sepr}{\tilde{e}^{\sss +}_{\sss R}}
\newcommand{\seplr}{\tilde{e}^{\sss +}_{\sss L,R}}
\newcommand{\seml}{\tilde{e}^{\sss -}_{\sss L}}
\newcommand{\semr}{\tilde{e}^{\sss -}_{\sss R}}
\newcommand{\semlr}{\tilde{e}^{\sss -}_{\sss L,R}}
\newcommand{\sepml}{\tilde{e}^{\sss \pm}_{\sss L}}
\newcommand{\sempl}{\tilde{e}^{\sss \mp}_{\sss L}}
\newcommand{\sepmr}{\tilde{e}^{\sss \pm}_{\sss R}}
\newcommand{\sempr}{\tilde{e}^{\sss \mp}_{\sss R}}
\newcommand{\sepmlr}{\tilde{e}^{\sss \pm}_{\sss L,R}}
\newcommand{\semplr}{\tilde{e}^{\sss \mp}_{\sss L,R}}
\newcommand{\semprl}{\tilde{e}^{\sss \mp}_{\sss R,L}}
\newcommand{\smul}{\tilde{\mu}_{\sss L}}
\newcommand{\smur}{\tilde{\mu}_{\sss R}}
\newcommand{\smulr}{\tilde{\mu}_{\sss L,R}}
\newcommand{\smupl}{\tilde{\mu}^{\sss +}_{\sss L}}
\newcommand{\smupr}{\tilde{\mu}^{\sss +}_{\sss R}}
\newcommand{\smuplr}{\tilde{\mu}^{\sss +}_{\sss L,R}}
\newcommand{\smuml}{\tilde{\mu}^{\sss -}_{\sss L}}
\newcommand{\smumr}{\tilde{\mu}^{\sss -}_{\sss R}}
\newcommand{\smumlr}{\tilde{\mu}^{\sss -}_{\sss L,R}}
\newcommand{\smupml}{\tilde{\mu}^{\sss \pm}_{\sss L}}
\newcommand{\smumpl}{\tilde{\mu}^{\sss \mp}_{\sss L}}
\newcommand{\smupmr}{\tilde{\mu}^{\sss \pm}_{\sss R}}
\newcommand{\smumpr}{\tilde{\mu}^{\sss \mp}_{\sss R}}
\newcommand{\smupmlr}{\tilde{\mu}^{\sss \pm}_{\sss L,R}}
\newcommand{\smumplr}{\tilde{\mu}^{\sss \mp}_{\sss L,R}}
\newcommand{\smumprl}{\tilde{\mu}^{\sss \mp}_{\sss R,L}}
\newcommand{\staul}{\tilde{\tau}_{\sss L}}
\newcommand{\staur}{\tilde{\tau}_{\sss R}}
\newcommand{\staulr}{\tilde{\tau}_{\sss L,R}}
\newcommand{\staupl}{\tilde{\tau}^{\sss +}_{\sss L}}
\newcommand{\staupr}{\tilde{\tau}^{\sss +}_{\sss R}}
\newcommand{\staupu}{\tilde{\tau}^{\sss +}_{\sss 1}}
\newcommand{\staupd}{\tilde{\tau}^{\sss +}_{\sss 2}}
\newcommand{\stauplr}{\tilde{\tau}^{\sss +}_{\sss L,R}}
\newcommand{\stauml}{\tilde{\tau}^{\sss -}_{\sss L}}
\newcommand{\staumr}{\tilde{\tau}^{\sss -}_{\sss R}}
\newcommand{\staumu}{\tilde{\tau}^{\sss -}_{\sss 1}}
\newcommand{\staumd}{\tilde{\tau}^{\sss -}_{\sss 2}}
\newcommand{\staumlr}{\tilde{\tau}^{\sss -}_{\sss L,R}}
\newcommand{\staupml}{\tilde{\tau}^{\sss \pm}_{\sss L}}
\newcommand{\staumpl}{\tilde{\tau}^{\sss \mp}_{\sss L}}
\newcommand{\staupmr}{\tilde{\tau}^{\sss \pm}_{\sss R}}
\newcommand{\staumpr}{\tilde{\tau}^{\sss \mp}_{\sss R}}
\newcommand{\staupmu}{\tilde{\tau}^{\sss \pm}_{\sss 1}}
\newcommand{\staupmd}{\tilde{\tau}^{\sss \pm}_{\sss 2}}
\newcommand{\staumpu}{\tilde{\tau}^{\sss \mp}_{\sss 1}}
\newcommand{\staumpd}{\tilde{\tau}^{\sss \mp}_{\sss 2}}
\newcommand{\staupmlr}{\tilde{\tau}^{\sss \pm}_{\sss L,R}}
\newcommand{\staumplr}{\tilde{\tau}^{\sss \mp}_{\sss L,R}}
\newcommand{\staumprl}{\tilde{\tau}^{\sss \mp}_{\sss R,L}}
\newcommand{\svl}{\tilde{\nu}_{\sss \ell,L}}
\newcommand{\svlb}{\bar{\tilde{\nu}}_{\sss \ell,L}} 
\newcommand{\sve}{\tilde{\nu}_{\sss e,L}}
\newcommand{\sveb}{\bar{\tilde{\nu}}_{\sss e,L}} 
\newcommand{\phino}{\tilde{\gamma}}
\newcommand{\Zino}{\tilde{Z}}
\newcommand{\Winop}{\tilde{W}^{\sss +}} 
\newcommand{\Winom}{\tilde{W}^{\sss -}}
\newcommand{\Winopm}{\tilde{W}^{\sss \pm}} 
\newcommand{\gino}{\tilde{g}} 
\newcommand{\hinop}{\tilde{H}^{\sss +}}
\newcommand{\hinom}{\tilde{H}^{\sss -}}
\newcommand{\hinopm}{\tilde{H}^{\sss \pm}} 
\newcommand{\lsp}{\tilde{\chi}}  

\newcommand{\mz}{M_{\sss Z}}
\newcommand{\mw}{M_{\sss W}}
\newcommand{\msel}{m_{\tilde{e}_L}}
\newcommand{\mser}{m_{\tilde{e}_R}}
\newcommand{\mselr}{m_{\tilde{e}_{L,R}}}
\newcommand{\msve}{m_{\tilde{\nu}_{e,L}}} 

\newcommand{\stopu}{\tilde{t}_1}
\newcommand{\stopd}{\tilde{t}_2}

\newcommand{\sthw}{\sin\theta_{\sss W}}
\newcommand{\cthw}{\cos\theta_{\sss W}}
\newcommand{\tgw}{\tan\theta_{\sss W}} 
\newcommand{\tgwq}{\tan^2\!\theta_{\sss W}} 
\newcommand{\sthwq}{\sin^2\!\theta_{\sss W}}
\newcommand{\cthwq}{\cos^2\!\theta_{\sss W}}
\newcommand{\sthwqq}{\sin^4\!\theta_{\sss W}}
\newcommand{\cthwqq}{\cos^4\!\theta_{\sss W}}
\newcommand{\sina}{\sin \alpha} 
\newcommand{\sinaq}{\sin^2\!\alpha} 
\newcommand{\sinda}{\sin 2\alpha}
\newcommand{\sinqa}{\sin^2\!\alpha} 
\newcommand{\cosa}{\cos \alpha} 
\newcommand{\cosaq}{\cos^2\!\alpha}
\newcommand{\cosda}{\cos 2\alpha} 
\newcommand{\cosqa}{\cos^2\!\alpha}
\newcommand{\tgb}{\tan \beta}
\newcommand{\tgbq}{\tan^2\!\beta}
\newcommand{\ctgb}{\cot \beta}
\newcommand{\ctgbq}{\cot^2\!\beta}
\newcommand{\sinb}{\sin \beta} 
\newcommand{\sinbq}{\sin^2\!\beta} 
\newcommand{\sindb}{\sin 2\beta}
\newcommand{\sinqb}{\sin^2\!\beta} 
\newcommand{\cosb}{\cos \beta} 
\newcommand{\cosbq}{\cos^2\!\beta}
\newcommand{\cosdb}{\cos 2\beta} 
\newcommand{\cosqb}{\cos^2\!\beta}
\newcommand{\tgthv}{\tan\theta_v}
\newcommand{\sthv}{\sin\theta_v}
\newcommand{\cthv}{\cos\theta_v}
\newcommand{\sdthv}{\sin 2\theta_v} 
\newcommand{\cdthv}{\cos 2\theta_v} 

% Macro-instructions
%
\newcommand{ \neut   }[1]{\tilde{\chi}^{\sss 0}_{#1}}  
\newcommand{ \hino   }[1]{\tilde{H}^{\sss 0}_{#1}}  
\newcommand{ \charp  }[1]{\tilde{\chi}^{\sss +}_{#1}}
\newcommand{ \charm  }[1]{\tilde{\chi}^{\sss -}_{#1}} 
%\newcommand{ \char   }[1]{\tilde{\chi}^{\sss \pm}_{#1}} 
\newcommand{ \charmp }[1]{\tilde{\chi}^{\sss \mp}_{#1}} 

\newcommand{ \mn     }[1]{m_{\tilde{\chi}^0_{#1}}} 
\newcommand{ \mc     }[1]{m_{\tilde{\chi}^{\pm}_{#1}}} 
\newcommand{ \mhino  }[1]{m_{\tilde{H}_{#1}}} 

\newcommand{ \N }[1]{ {\tilde N}_{#1} }
\newcommand{ \C }[1]{ {\tilde C}_{#1} }
\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{\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_{\rm 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{\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}{{\tt FeynHiggs}}
\newcommand{\fhf}{{\tt FeynHiggsFast}}
\newcommand{\suf}{{\tt SusyFIRE}}
\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{\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}}


