\documentstyle[fleqn,12pt]{article}
% put your own definitions here:
%   \newcommand{\cZ}{\cal{Z}}
%   \newtheorem{def}{Definition}[section] 
%   ...
\input epsf
% \documentstyle[12pt]{article}
 
% Redefinition of standard commands:
\let\un=\b
\let\ced=\c
\let\du=\d
\let\um=\H
\let\sll=\l
\let\Sll=\L
\let\slo=\o
\let\Slo=\O
\let\tie=\t
\let\br=\u
\let\SEC=\S 
\def\cA{{\cal A}}
\def\cB{{\cal B}}
\def\cC{{\cal C}}
\def\cD{{\cal D}}
\def\cE{{\cal E}}
\def\cF{{\cal F}}
\def\cG{{\cal G}}
\def\cH{{\cal H}}
\def\cI{{\cal I}}
\def\cJ{{\cal J}}
\def\cK{{\cal K}}
\def\cL{{\cal L}}
\def\cM{{\cal M}}
\def\cN{{\cal N}}
\def\cO{{\cal O}}
\def\cP{{\cal P}}
\def\cQ{{\cal Q}}
\def\cR{{\cal R}}
\def\cS{{\cal S}}
\def\cT{{\cal T}}
\def\cU{{\cal U}}
\def\cV{{\cal V}}
\def\cW{{\cal W}}
\def\cX{{\cal X}}
\def\cY{{\cal Y}}
\def\cZ{{\cal Z}}
% Greek alphabet
\def\a{\alpha}
\def\b{\beta}
% \def\c{\chi}
\def\d{\delta}      \def\D{\Delta}
\def\e{\epsilon}
\def\eps{\varepsilon}
\def\f{\phi}        \def\F{\Phi}
\def\vf{\varphi}
\def\g{\gamma}      \def\G{\Gamma}
\def\h{\eta}
\def\i{\iota}
\def\j{\psi}        \def\J{\Psi}
\def\k{\kappa}
\def\l{\lambda}   \def\L{\Lambda}
\def\m{\mu}
\def\n{\nu}
\def\o{\omega}   \def\O{\Omega}
\def\p{\pi}      \def\P{\Pi}
\def\q{\theta}   \def\Q{\Theta}
\def\r{\rho}
\def\s{\sigma}   \def\S{\Sigma}
\def\t{\tau}
\def\u{\upsilon}  \def\U{\Upsilon}
\def\x{\xi}      \def\X{\Xi}
\def\z{\zeta}
% Bars on quarks:
\def\qb{\bar{q}}
\def\ub{\bar{u}}
\def\db{\bar{d}}
\def\cb{\bar{c}}
\def\sb{\bar{s}}

% Bra-Kets
\def\bra{\langle}
\def\ket{\rangle}
% KM parameters
\def\Vcd{V_{cd}}
\def\Vus{V_{us}}
\def\Vcs{V_{cs}}
\def\Vud{V_{ud}}
\def\Vbc{V_{cb}}
\def\Vbu{V_{ub}}
\def\Vtd{V_{td}}
\def\Vts{V_{ts}}
\def\Vtb{V_{tb}}
\def\Im{\mbox{Im}}
\def\bg#1{\mbox{\boldmath $ #1 $}}     % FOR BOLD GREEK LETTERS
 
%Some groups
 
\def\gws{{SU(2)_L \times U(1)_Y}}
\def\glrs{{SU(2)_L \times SU(2)_R \times U(1)_{B-L}}}
\def\gps{{SU(2)_L \times SU(2)_R \times SU(4)_C}}
\def\gsm{{SU(2)_L \times U(1)_Y \times SU(3)_c}}
\def\gr{SU(2)_R \times U(1)_{B-L}}
 
\def\bo{{\raise.15ex\hbox{\large$\Box$}}}
\def\del{\partial}
\def\de{\nabla}
\def\dag{^{\dagger}{}}
\def\rarr{\rightarrow}
\def\larr{\leftarrow}
\def\iff{\leftrightarrow}
\def\conj{^{\raisebox{-.4cm}{\normalsize *}}}
 
\def\Bar{\overline}
\def\Tilde#1{\widetilde {#1}}
\def\Hat#1{\widehat {#1}}
\def\bra#1{\left\langle #1\right|}
\def\ket#1{\left| #1\right\rangle}
\def\vev#1{\left\langle #1\right\rangle}
\def\abs#1{\left| #1\right|}
\def\mag#1{\abs{\vec{#1}}}
\def\chain#1{\begin{array}{c} \\ \longrightarrow \\ {\scriptstyle #1}
    \end{array}}
\def\ordless{{\lower2mm\hbox{$\,\stackrel{\textstyle <}{\sim}\, $}}}
\def\ordmore{{\lower2mm\hbox{$\,\stackrel{\textstyle >}{\sim}\, $}}}
\def\matelem#1#2#3{\vev{#1 \left| \vphantom{#1} #2 \vphantom{#3} \right| #3}}
\def\momint#1{\int {d^4 #1 \over (2\p)^4}}
 
\newtoks\slashfraction
\slashfraction={.13}
\def\slash#1{\setbox0\hbox{$\, #1$}
\setbox0\hbox to \the\slashfraction\wd0{\hss \box0}/\box0}
\def\Slash#1{\rlap/#1}
 
\def\leftrightarrowfill{$\mathsurround=0pt \mathord\leftarrow \mkern-6mu
        \cleaders\hbox{$\mkern-2mu \mathord- \mkern-2mu$}\hfill
        \mkern-6mu \mathord\rightarrow$}
\def\overleftrightarrow#1{\vbox{\ialign{##\crcr
        \leftrightarrowfill\crcr\noalign{\kern-1pt\nointerlineskip}
        $\hfil\displaystyle{#1}\hfil$\crcr}}}
\def\startarray{\left( \begin{array}}
\def\finarray{\end{array} \right)}
\def\starteq{%\vspace{.1in}
\begin{eqnarray}}
\def\fineq{\end{eqnarray}
%\newline
}
 
\def\half{{1 \over 2}}
\def\partder#1#2{{\partial #1\over\partial #2}}
\def\secder#1#2#3{{\partial^2 #1\over\partial #2 \partial #3}}
 
\catcode`@=11
\def\underline#1{\relax\ifmmode\@@underline#1\else
$\@@underline{\hbox{#1}}$\relax\fi}
\catcode`@=12
 
% mboxes in math mode

\def\to{\mbox{-}} 
\def\tr{\mbox{tr}\,}
\def\Tr{\mbox{Tr}\,}
\def\csch{\mbox{csch}}
\def\sech{\mbox{sech}}
\def\hc{\mbox{h.c.}}
\def\cc{\mbox{c.c.}}
\def\re{\mbox{Re\,}}
\def\im{\mbox{Im\,}}
\def\sub#1{_{{\rm #1}}}
\def\sup#1{^{{\rm #1}}}
 
%       the stuff below defines \eqalign and \eqalignno in such a
%       way that they will run on Latex
\newskip\humongous \humongous=0pt plus 1000pt minus 1000pt
\def\caja{\mathsurround=0pt}
\def\eqalign#1{\,\vcenter{\openup2\jot \caja
        \ialign{\strut \hfil$\displaystyle{##}$&$
        \displaystyle{{}##}$\hfil\crcr#1\crcr}}\,}
\newif\ifdtup
 
% Referencing
 
\def\NP#1{{\em Nucl. Phys.} {\bf B{#1}}}
\def\PL#1{{\em Phys. Lett.\ }{\bf B{#1}}}
\def\PRL#1{{\em Phys. Rev. Lett.\ }{\bf {#1}}}
\def\PRD#1{{\em Phys. Rev.\ }{\bf D{#1}}}
\def\Astr#1{{\em Astron. and Astrophys.\ }{\bf{#1}}}
\def\ApJ#1{{\em Astrophys.~J.} {~\bf{#1}}}
\def\PRep#1{{\em Phys. Rep.\ }{\bf{#1}}}
\def\japan#1{{\em Prog. Theor. Phys.\ }{\bf{#1}}}
\def\italy#1{{\em Nuovo Cimento\ }{\bf {#1}C}}
\def\zph#1{{\em Z.~Phys.\ }{\bf C{#1}}}
\def\annphys#1{{\em Ann. Phys. (NY)\ }{\bf #1}}
\def\jetp#1{{\em Sov. Phys. JETP} {\bf #1}}
\def\IJMP#1{{\em Int. J.~Mod. Phys.\ }{\bf A{#1}}}
\def\amjphys#1{{\em Am.~J. Phys.\ }{\bf #1}}
 
\def\pref#1{{$^{\hbox{\scriptsize \cite{#1}}}$}}
\def\textcite#1{Ref.~{\cite{#1}}}
\def\Ref#1{(\ref{#1})}
\def\Eqn#1{Eq. \Ref{#1}}
\def\Sec#1{Sec.~\ref{#1}}
\def\next{\nonumber \\}
\def\thefootnote{\fnsymbol{footnote}}
 
% Authors, Frequent collaborators
 
\def\dept{Department of Physics}
\def\phast{Department of Physics and Astronomy}
\def\author#1#2{{\bf #1} \\ {\em #2}\vspace{5mm}}
\def\ecolepoly{Centre de Physique Th\'eorique, \\
Centre National de la Recherche Scientifique, UPR A0014, \\  
Ecole Polytechnique, 91128 Palaiseau Cedex, France}
\def\pham{\author {T. N. Pham} {\ecolepoly}}
%
 
\def\sect#1{\bigskip\medskip \goodbreak \noindent {\bf #1}
    \nobreak \medskip}
\def\Section#1{\begin{center}\section{\hspace{-.3cm}{: #1}}\end{center}}
\def\letsect#1{\vspace{.1in} \noindent\underline{\em #1} $\bullet$}
\def\Appendix#1{\addtocounter{section}{1} \setcounter{equation}{0}
\begin{center} {\Large \bf Appendix
    \thesection : {#1}}\end{center}\nopagebreak}
\def\singleappendix#1{\appendix \addtocounter{section}{1}
    \setcounter{equation}{0} \vspace{.6in}
    \begin{center} {\Large \bf Appendix : {#1}}\end{center}\nopagebreak}
\def\INSERT#1{{\samepage
    \begin{center} \rule{13cm}{.4mm} \\
    {\Large \bf $\Longrightarrow$ [\{(Insert here)\}]
    $\Longleftarrow$} \\ {\Large \em #1} \\ \rule{13cm}{.4mm}
    \end{center}}}
 
\def\bold#1{\setbox0=\hbox{$#1$}%
     \kern-.025em\copy0\kern-\wd0
     \kern.05em\copy0\kern-\wd0
     \kern-.025em\raise.0433em\box0 }

%\typein {      >>> input either PPT or JOUR styles  <<<}


\font\mysevenrm=cmr7
\font\mysevenmi=cmmi7
                   
%%%%%%%%%%%%%%%%%%%%%%%% letterhead %%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%% comments %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Style commands
\topmargin=0in
\headheight=0in
\headsep=0in
\textheight=9in \textwidth=6in
\footheight=2ex \footskip=5ex
\oddsidemargin=.25in \evensidemargin=.25in
\hsize=6in
\parskip=0pt
\lineskip=0pt
\abovedisplayskip=3mm plus.3em minus.5em
\belowdisplayskip=3mm plus.3em minus.5em
\abovedisplayshortskip=2mm plus.2em minus.4em
\belowdisplayshortskip=2mm plus.2em minus.4em
\baselineskip=14pt
\def\baselinestretch{1.2}
 
% TitlePage
\def\title#1#2#3#4#5{\thispagestyle{empty}
        \begin{center} \vspace*{1cm} { \bf #3} \\[.5in] {#4{}}
        \end{center} \vfill \centerline{ ABSTRACT}
   {\nopagebreak \noindent\begin{quotation}\noindent {\small #5}
   \end{quotation}} \vfill {#2} \hfill\begin{tabular}{r} {#1} 
        \end{tabular}  \newpage
        \def\thefootnote{\arabic{footnote}}}
%        \def\thefootnote{\sharp\arabic{footnote}}}
%
% Different styles for the bibliography: --------------------------------
%
\def\prefer{\section*{}
    \list{[\arabic{enumi}]}{\usecounter{enumi}\settowidth\labelwidth{[000]}
      \leftmargin\labelwidth\advance\leftmargin\labelsep \rightmargin=0pt}
        \small \sfcode`\.=1000\relax}
\let\endprefer=\endlist
  %
\def\refer#1{\section*{\large \sc {#1}}
    \list{\arabic{enumi}.}{\usecounter{enumi}\settowidth\labelwidth{[000]}
      \leftmargin\labelwidth\advance\leftmargin\labelsep \rightmargin=0pt}
        \raggedright \small \sfcode`\.=1000\relax}
\let\endrefer=\endlist
  %
\def\ReFer#1#2{\section*{\large\sc#1}
    \list{[\arabic{enumi}]}{\usecounter{enumi}\settowidth\labelwidth{#2}
      \leftmargin\labelwidth\advance\leftmargin\labelsep \rightmargin=0pt}
        \raggedright \small \sfcode`\.=1000\relax}
\let\endReFer=\endlist
  %
\def\REFER#1#2{\section*{\large\sc#1}
    \list{#2 {enumi}.}{\usecounter{enumi}\settowidth\labelwidth{[000]}
      \leftmargin\labelwidth\advance\leftmargin\labelsep \rightmargin=0pt}
        \raggedright \small \sfcode`\.=1000\relax}
\let\endREFER=\endlist
% ------------------------------------------------------------------------
%\def\startbib{\vspace{1in}\begin{thebibliography}{000}
\def\startbib{\vspace{1in}\begin{refer}{References}
\small\frenchspacing\nopagebreak}
\def\endbib{\end{refer} \normalsize \nonfrenchspacing}
\def\startfig{\newpage \centerline{{\sl Figure captions}} \begin{itemize}}
\def\figitem#1{\item[{\bf Fig.~{#1}:}\ ]}
\def\endfig{\end{itemize}}
 
%\begin{document}







\newcommand{\be}{\begin{equation}}
\newcommand{\ee}{\end{equation}}
\newcommand{\bea}{\begin{eqnarray}}
\newcommand{\eea}{\end{eqnarray}}
\newcommand{\nn}{\nonumber}
\newcommand{\dd}{\displaystyle}
%\newcommand{\bra}[1]{\left\langle #1 \right|}
%\newcommand{\ket}[1]{\left| #1 \right\rangle}
%\newcommand{\qq}{<0|{\bar q} q|0>}
%\newcommand{\spur}[1]{\not\! #1 \,}
\newcommand{\ttbs}{\char'134}
\newcommand{\AmS}{{\protect\the\textfont2
  A\kern-.1667em\lower.5ex\hbox{M}\kern-.125emS}} 
%\pagestyle{empty}
% add words to TeX's hyphenation exception list
\hyphenation{author another created financial paper re-commend-ed}
\begin{document}
% declarations for front matter
\title{September 2000} {PC084.0900}
{ $B \rightarrow K\pi$ DECAYS  AND THE  WEAK PHASE ANGLE $\gamma\dag$  
\footnotetext{${\dag}$  
Talk given at the QCD Euroconference 00, Montpellier 6-12 July 2000}}
{\author{T. N. Pham} {Centre de Physique Th\'eorique, \\
Centre National de la Recherche Scientifique, UMR 7644, \\  
Ecole Polytechnique, 91128 Palaiseau Cedex, France}}   
{The large 
branching ratios  for  $B \rightarrow K\pi$ decays as observed  by the CLEO
Collaboration  indicate that penguin interactions
contribute a major part to the decay rates and provide an interference
between the Cabibbo-suppressed tree and penguin contributions resulting in
a CP-asymmetry between the $B \rightarrow K\pi$ and its charge conjugate
mode. The CP-averaged decay rates  depend also on  the weak phase
$\gamma$ and give us a  determination of this phase. In this talk, I
would like to report on a recent analysis of $B \rightarrow K\pi$ decays
using factorisation model with final state interaction 
phase shift included. We find that factorisation
seems to describe qualitatively the latest CLEO data. We also obtain   
 a relation for the branching
ratios independent of the strength of the strong penguin interactions.
This relation gives a central value of $0.60 \times 10^{-5}$ for 
${\mathcal B}(\bar{B}^{0} \rightarrow \bar{K}^{0}\pi^{0})$, somewhat
smaller than the latest CLEO measurement. 
We also find that a ratio obtained from the 
CP-averaged $B \rightarrow K\pi$ decay rates could be used
to test the factorisation model and to determine the weak angle $\gamma$
with more precise data, though the latest CLEO data seem to favor
$\gamma$ in the range  $90^{\circ}-120^{\circ}$.}



With the measurement of all the four $B \rightarrow K\pi$ branching ratios, 
we seem to have a qualitative understanding of the $B \rightarrow K\pi$ decays. 
The  measured CP-averaged branching ratios(${\mathcal B}$)   by 
the CLEO Collaboration \cite{cleo2} show  that the penguin interactions
dominates the $B \rightarrow K\pi$ decays, as predicted by factorisation. 
The strong penguin
amplitude, because of the large CKM factors, becomes much larger than the
tree-level terms which are Cabibbo-suppressed and the non-leptonic
interaction for $B \rightarrow K\pi$ is dominated by an $ I = 1/2 $ amplitude.
This is borne out by the CLEO data which give
$ {\mathcal B}(B^{-} \rightarrow \bar{K}^{0}\pi^{-}\simeq 2{\mathcal B}(B^{-} \rightarrow
K^{-}\pi^{0}) $ and 
${\mathcal B}(B^{-} \rightarrow \bar{K}^{0}\pi^{-}) \simeq {\mathcal B}(\bar{B}^{0}
\rightarrow K^{-}\pi^{+})$~:
\begin{eqnarray}  
{\mathcal B} (B^{+} \rightarrow K^{+} \pi^0)
&=&(11.6^{+3.0+1.4}_{-2.7-1.3}) \times 10^{-6},  \nonumber  \\
{\mathcal B} (B^{+} \rightarrow K^0 \pi^{+})  &=& 
(18.2^{+4.6}_{-4.0}\pm 1.6) \times 10^{-6},   \nonumber  \\
{\mathcal B} (B^0\rightarrow K^{+} \pi^{-})  &=& 
(17.2^{+2.5}_{-2.4}\pm 1.2) \times 10^{-6}, \nonumber \\ 
{\mathcal B} (B^0  \rightarrow K^0\pi^0)\, \  &=& \
(14.6^{+5.9+2.4}_{-5.1-3.3}) \times 10^{-6}.
\end{eqnarray}
If the strength of the  interference between the tree-level and penguin
contributions is known, a determination of the weak phase $\gamma$ 
could be done in principle. Previous works \cite{Deshpande,Ali} shows that
factorisation model produces sufficient $B \rightarrow K\pi$ decay rates, in   
qualitative agreement with the CLEO measured values. Also, as argued  
in \cite{Neubert2}, for these very 
energetic decays,  because of color transparency, factorisation should be 
a good approximation for $B \rightarrow K\pi$ decays if the Wilson coefficients 
are evaluated at a scale $\mu =O(m_{b})$. Infact, recent hard scattering
calculations with perturbative QCD shows that factorisation is valid
up to corrections of the order $\Lambda_{\rm QCD}/m_{b}$ \cite{Beneke}. 
It is thus 
encouraging to use factorisation to analyse the $B \rightarrow K\pi$ decays,
bearing in mind that there are important 
theoretical uncertainties in the long-distance hadronic matrix elements, 
as the heavy to light form factors for the vector current and the value of
the current 
$s$ quark mass are currently not determined with good accuracy. 
In this talk, I
would like to report on a recent work \cite{Isola} on the $B \rightarrow K\pi$ 
decays as a possible way to measure the angle $\gamma$ and to see direct 
CP violation. 

In the standard model, the effective Hamiltonian for $B \rightarrow K\pi$ decays
are given by \cite{Buras,Ciuchini,Fleischer1,Kramer,Deshpande1},
\bea
\kern -0.7cm H_{\rm eff} = {G_F\over \sqrt{2}}[V_{ub}V^*_{us}(
c_1O_{1}^{u} + c_2O_{2}^{u})  + V_{cb}V^*_{cs}(c_1O_{1}^{c} 
 + \  c_2 O_{2}^{c}) -\sum_{i=3}^{10}V_{tb}V^*_{ts} c_iO_i]  
+ {\rm h.c.}\,.
\label{hw}
\eea
in standard notation.
At next-to-leading logarithms, $c_{i}$ take the form of an effective
Wilson coefficients $c_{i}^{\rm eff}$ which  contain also the penguin 
contribution from the $c$ quark loop and 
are given in \cite{Fleischer1,Deshpande1}.

The parameters  $V_{ub}$ etc. are the flavor-~changing charged current 
couplings of the weak gauge boson $W^{\pm}$ with the quarks as given by 
the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix 
$V$. $V$ is usually defined as the unitary transformation
relating the the weak interaction eigenstate of quarks to their mass 
eigenstate \cite{PDG}:
\be
\pmatrix{d' \cr s' \cr b'}= \pmatrix{V_{ud}& V_{us} & V_{ub} \cr
V_{cd} & V_{cs} & V_{cb} \cr 
V_{td} & V_{ts} & V_{tb}} \pmatrix{d \cr s \cr
b}
\label{Vckm} 
\ee
where $d,s,b$ and $d',s',b'$ are respectively the  mass eigenstates
and weak interaction eigenstates for the charge $Q=-1/3$ quarks. 
Since the neutral current is not affected by the unitary transformation
on the quark fields, 
flavor-changing neutral current is absent at the tree-level as implied
by the GIM mechanism. The unitarity condition $VV^{\dag}=1$ 
gives, for the (db) elements relevant to $B$ decays \cite{PDG}
\be
V_{ud}V_{ub}^{*} + V_{cd}V_{cb}^{*} + V_{td}V_{tb}^{*} =0 
\label{db}  
\ee
This can be represented by a triangle \cite{PDG} with the three 
angles $\alpha$, $\beta$
and $\gamma$ expressed in terms of the CKM matrix elements as \cite{Nir}:
\bea
&& \alpha = arg(-V_{td}V_{tb}^{*}/V_{ud}V_{ub}^{*})  \nonumber\\
&&  \beta = arg(-V_{cd}V_{cb}^{*}/V_{td}V_{tb}^{*})  \nonumber\\
&&  \gamma = arg(-V_{ud}V_{ub}^{*}/V_{cd}V_{cb}^{*})
\label{arg}
\eea
The angle $\gamma$  enters the $B \rightarrow K\pi$ decay amplitudes  through
the factor 
$V_{ub}V_{us}^*/V_{tb}V_{ts}^* $ which can be approximated by
$-(|V_{ub}|/|V_{cb}|)\times (|V_{cd}|/|V_{ud}|)\exp(-i\,\gamma)$
after neglecting terms of the order $O(\lambda^{5})$ in the (bs)
unitarity triangle, $\lambda $ being the Cabibbo angle in the 
Wolfenstein parametrisation of the CKM quark mixing matrix. The $B \rightarrow
K\pi$ decay amplitudes, expressed in terms of the $I=1/2$ and $I=3/2$
isospin amplitudes are given by \cite{Deshpande},

\bea
A_{K^-\pi^0} &=& 
 {2\over 3} B_3e^{i\delta_3} + \sqrt{{1\over 3}} 
(A_1+B_1)e^{i\delta_1}, \ \ \  \nonumber\\
A_{\bar K^0\pi^-} &=&
{\sqrt{2}\over 3} B_3e^{i\delta_3} - \sqrt{{2\over 3}} 
(A_1+B_1)e^{i\delta_1},\nonumber\\
A_{K^-\pi^+} &=&
{\sqrt{2}\over 3} B_3e^{i\delta_3} + \sqrt{{2\over 3}} 
(A_1-B_1)e^{i\delta_1},\ \ \ \nonumber\\
A_{\bar K^0\pi^0} &=& 
{2\over 3} B_3e^{i\delta_3} - \sqrt{{1\over 3}} 
(A_1-B_1)e^{i\delta_1},
\label{amplitude}
\eea
$A_{1}$ is the sum of  the strong penguin $A_{1}^{\rm S} $
and the $I=0$  tree-level $A_{1}^{\rm T}$ as well as the $I=0$ 
electroweak penguin $A_{1}^{\rm W}$ contributions 
to the $B \rightarrow K\pi$ $I=1/2$ amplitude;  similarly  
$ B_{1}$ is the sum of the
$I=1$ tree-level $B^{T}_{1}$ and electroweak penguin $B^{W}_{1}$
contribution to the $I=1/2$ amplitude, and $B_{3}$ is the sum of the
$I=1$  tree-level $B^{T}_{3}$ and electroweak penguin $B^{W}_{3}$
contribution to the $I=3/2$ amplitude. $\delta_1 $ and $ \delta_3$ are,
respectively, the  elastic  $\pi K \rightarrow \pi K$ $I=1/2$ and $I=3/2$ final 
state interaction (FSI) phase shift at the $B$ mass. 
The inelastic FSI contributions 
are also included through the internal quark loop contributions to the
penguin operators, for which the Wilson coefficients now have 
an absorptive  part and are given in \cite{Fleischer1,Deshpande1,Hou}.
The $B \rightarrow K\pi$ isospin amplitudes in the factorisation model 
are given by \cite{Deshpande},
\bea
&&  A_1^{\rm T} = i{\sqrt{3}\over 4}\, V_{ub}V_{us}^*\, r\, 
a_2,\nonumber\\
&& B_1^{\rm T} = i {1\over 2\sqrt{3}}\, V_{ub}V_{us}^*\, r 
\left[-{1\over 2}a_2 
+a_1 X\right],\nonumber\\
&& B_3^{\rm T} = i{1\over 2}\, V_{ub}V_{us}^*\, r 
\left[a_2 + a_1 X\right],
\nonumber\\
&& A_1^{\rm S} = -i{\sqrt{3}\over 2}\, V_{tb}V_{ts}^* r
\left[ a_4 + a_6 Y\right],
  \ B_1^{\rm S} = B_3^{\rm S} = 0,\nonumber\\
&& A_1^{\rm W} = -i{\sqrt{3}\over 8}\, V_{tb}V_{ts}^*r
\left[a_8Y +  a_{10}\right],
\nonumber\\
&& B_1^{\rm W} = i{\sqrt{3}\over 4}\, V_{tb}V_{ts}^*r
\left[{1\over 2}a_8Y 
+ {1\over 2}a_{10} 
+\left(a_7 - a_9\right)X\right],
\nonumber\\
&& B_3^{\rm W} = - i{3\over 4}\, V_{tb}V_{ts}^* r
\left[ a_8Y + a_{10} 
- (a_7 - a_9)X\right] 
\eea
where $r = G_F\, f_K F^{B\pi}_0(m^2_K) (m_B^2-m_\pi^2)$ ,\qquad
$X= (f_\pi/f_K)[F^{BK}_0(m^2_\pi)/F^{B\pi}_0(m^2_K)] 
[(m_B^2-m^2_K)/(m_B^2-m_\pi^2)]$,
$Y = 2m^2_K/[(m_s+ m_q)(m_b-m_q)]$ with $q=u,\ d$ for 
$\pi^{\pm,0}$ final states, respectively. In this analysis, 
$f_{\pi} =133 \,\rm MeV$, $f_{K} =158\,\rm MeV$, 
$ F_{0}^{B\pi}(0) = 0.33$, $ F_{0}^{BK}(0) = 0.38$ \cite{Ali,Bauer}~;
$|V_{cb}|=0.0395, |V_{cd}|=0.224 $ and $|V_{ub}|/|V_{cb}|=0.08$ \cite{PDG}.
The value of $m_{s}$ is not known to a good accuracy, but a value 
around $(100-120\rm )\, MeV$
inferred from $m_{K^{*}} - m_{\rho}$, $m_{D_{s}^{+}}- m_{D^{+}}$ 
and $m_{B_{s}^{0}}- m_{B^{0}}$ mass differences
\cite{Colangelo} seems not unreasonable and in this work, 
we use $m_{s} =120\,\rm MeV$.
$a_{j}$ are  the effective Wilson coefficients after Fierz reordering
in factorisation model and are given by \cite{Isola}
\bea
&& a_1 = \ 0.07    , \ \ \ a_2=  \ 1.05  , \nonumber\\
&& a_4 = -0.043 - 0.016i , \ \ a_6= -0.054 - 0.016i  ,
\label{coeff}
\eea
for the contributions from the
tree-level and  the strong penguin operators at $N_c =3$ and 
$m_{b} = 5.0\, \rm GeV$. 
The strong penguin contribution   $P = a_4 + a_6 Y$, 
as obtained from Eq.(\ref{coeff}) is enhanced by the
charm quark loop which increases the amplitude by $30\%$ as pointed out 
in \cite{Fleischer1}. This enhancement brings the predicted branching
ratios closer to the CLEO measured values, as shown in Fig.1. where
the CP-averaged $B \rightarrow K\pi$ branching ratios obtained for 
$\gamma = 70^\circ$ \cite{Ali}, are plotted against
the rescattering phase difference $\delta = \delta_{3} - \delta_{1}$. 
%\vspace{-2cm}
\begin{figure}[hbp]
\centering
\leavevmode
\epsfxsize=4.8cm
\epsffile{fig10.eps}
%\centerline{ \DESepsf(fig10.eps width 6 cm)}
%\vspace{-1cm}
\caption{${\mathcal B}(B\rightarrow K\pi)$ vs. $\delta$ for $\gamma = 70^\circ$.
The curves  
(a), (b), (c), (d)  are for the CP-averaged branching ratios
$B^-\rightarrow K^-\pi^0,\ \bar K^0 \pi^-$ and 
$\bar B^0 \rightarrow K^-\pi^+,\ \bar K^0\pi^0$, respectively.}
\label{fig:BR}
\end{figure}

For a determination of $\gamma$, two 
quantities obtained from the sum
of the two CP-averaged decay rates 
 $\Gamma_{B^-} = \Gamma(B^{-} \rightarrow K^- {\pi}^0) + 
\Gamma(B^{-} \rightarrow \bar{K}^{0} {\pi}^{-})$ and 
$\Gamma_{B^0} = \Gamma(\bar{B}^{0} \rightarrow K^- {\pi}^+) + 
\Gamma(\bar{B}^{0} \rightarrow \bar{K}^{0} {\pi}^{0})$ which are independent of
$\delta$ could be used \cite{Isola}.
As the CP-averaged $B \rightarrow K\pi$ decay rates depend  on $\gamma$,
the  computed partial rates $\Gamma_{B^{-}} $ 
and $\Gamma_{B^0} $ would now lie
between the upper and lower limit corresponding to $\cos(\gamma) =1 $ 
and $\cos(\gamma) =-1 $, respectively.
As shown in Fig.2, where the corresponding CP-averaged branching
ratios (${\mathcal B}_{B^0} $ and
${\mathcal B}_{B^-}$)
for $\Gamma_{B^{-}} $ and $\Gamma_{\bar{B}^{0}} $ are plotted against 
$\gamma$, the factorisation model values with the 
BWS form factors \cite{Bauer}
seem somewhat smaller
than the CLEO central values by about $10-20\% $. Also,  
${\mathcal B}_{B^{-}} > {\mathcal B}_{\bar{B}^{0}} $
while the data give
${\mathcal B}_{B^{-}} < {\mathcal B}_{\bar{B}^{0}} $
by a small amount which could be due to a large
measured $\bar{B}^{0} \rightarrow \bar{K}^{0}\pi^{0} $ branching ratio.    
%\vspace{-2cm}
\begin{figure}[hbp]
\centering
\leavevmode
\epsfxsize=4.8cm
\epsffile{fig130.eps}
%\vspace{-0.5cm}
\caption{${\mathcal B}_{B^-}$ (a), ${\mathcal B}_{\bar{B}^0}$ (b), $\Delta$ (c) vs. $\gamma$}
\label{fig:facttest}
\end{figure}
%\vspace{-0.5cm}

Note that smaller values for the form factors 
could easily accommodate the latest CLEO measured values, if a smaller
value for $m_{s}$, e.g, in the range $(80-100)\,\rm MeV$ is used. What
one learns from this analysis is that $B \rightarrow K\pi$ decays are 
penguin-dominated and the strength of the penguin interactions as obtained by
perturbative QCD, produce sufficient $B \rightarrow K\pi$ decay rates and that
factorisation seems to work with an accuracy better than a factor of 2,
considering large uncertainties from the form factors and possible
non-factorisation terms inherent in the factorisation model
and the uncertainties in the penguin amplitude which is sensitive to the
current $s$ quark mass. Since  the four $B \rightarrow K\pi$ decay rates 
depend on only three amplitudes $A_{1} $, $B_{1} $ and $B_{3} $, it is
possible to derive a relation
between the decay rates independent of $A_{1}$. Thus, the
quantity $\Delta $ obtained from the decay rates,
\bea
&& \kern -0.7cm \Delta =\left\{ \Gamma(B^{-} \rightarrow \bar{K}^{0} {\pi}^{-})+\Gamma(\bar{B}^{0}
\rightarrow K^- {\pi}^+) \right. \nonumber\\
&& \kern -0.7cm  \ \  \ \ \ \ -  \left. 2\left[\Gamma(B^{-} \rightarrow K^- {\pi}^0)
+\Gamma(\bar{B}^{0} \rightarrow
     \bar{K}^{0} {\pi}^{0})\right]\right\}\tau_{B^0}\nonumber\\
&&\kern -0.7cm \ \ \ =  \left[-{4\over 3}|B_3|^2-{8\over {\sqrt 3}}{\rm Re} (B_3^* B_1 {\rm  e}^{i\delta})\right](C\tau_{B^0}) 
\label{delta}
\eea
is independent of  the strong penguin term. It is given by the
tree-level and electroweak penguin contributions. As can be seen 
from Fig.2, where its values for  $\delta =0$ 
are plotted against $\gamma $. $\Delta$ is of the order $O(10^{-6})$
compared with  ${\mathcal B}_{B^{-}} $ and ${\mathcal
  B}_{\bar{B}^{0}} $ which 
are in the range  $(1.6-3.0)\times 10^{-5}$. Thus, to this
level of accuracy, we can put $\Delta \simeq 0$ and obtain the relation
($r_b=\tau_{\bar{B}^0}/ \tau_{B^-}$).
\be
{ r_b}{\mathcal B}_{\bar{K}^{0}
    {\pi}^{-}}+{\mathcal B}_{K^- {\pi}^+} = 
2\,\left[{\mathcal B}_{\bar{K}^{0} {\pi}^{0}}+
{ r_b}{\mathcal B}_{ K^- {\pi}^0}\right]\,.
\label{kpi0}
\ee
which can be used to test factorisation or to predict 
${\mathcal B}(\bar B^0\rightarrow \bar K^0 {\pi}^0)$ in terms of the 
other measured branching ratios. Eq.(\ref{kpi0}) then predicts a central value 
${\mathcal B}(\bar{B}^{0} \rightarrow\bar{K}^{0} 
{\pi}^{0})= 0.60\times 10^{-5}$.

Since  ${\mathcal B}_{\bar{K}^{0} {\pi}^{0}} $ is not known with good
accuracy at the moment, it is useful to use another quantity, defined as
\be
  { r_b}{\mathcal B}_{\bar{K}^{0}
    {\pi}^{-}}+{\mathcal B}_{K^- {\pi}^+} = (C_1\tau_{\bar{B}^{0}})
  \left[{1\over 3} |B_3|^2 +(|A_1|^2 + 
|B_1|^2)-{2\over \sqrt 3}{\rm Re}(B_3^* B_1 {e}^{i\delta})\right]
\label{B23}
\ee
which contains a negligible $\delta$-dependent term of the order 
$O(10^{-7})$. The quantity $R$ defined as
\be
R={{{\mathcal B}(B^-\rightarrow K^{-} {\pi}^{0})
+{\mathcal B}( B^-\rightarrow \bar{K}^{0} {\pi}^{-})}\over{
{\mathcal B}(B^-\rightarrow\bar{K}^{0}
    {\pi}^{-})+{{\mathcal B}(\bar B^0\rightarrow K^- {\pi}^+)}/{r_b}}}\,\, .
\label{R}
\ee
is thus essentially independent of $\delta$ and could also be used to
obtain $\gamma$, as it does not suffer from large uncertainties in the
form factors and in the CKM parameters. As can be seen in Fig.3, 
it is not possible to deduce a value for $\gamma$ with the present
data which give $R=(0.80\pm 0.25) $ 
as the  prediction for $R$ lies within
the experimental errors. If we could reduce the experimental
uncertainties to a level of less than $10\%$, we might be able to 
give a value for $\gamma$. Thus it is important to measure 
$B\rightarrow K\pi$ branching ratios to a high precision. Also shown in
Fig.3 are two other quantities more sensitive to $\gamma$, but involved 
${\mathcal B}(\bar B^0\rightarrow \bar K^0 {\pi}^0) $ and are given as \cite{Isola} 
\be
R_{1} = {\Gamma_{B^{-}} \over \Gamma_{\bar B^0}} \, \ ,\
R_{2} = {\Gamma_{B^{-}} \over (\Gamma_{B^{-}} +\Gamma_{\bar B^0})}
\label{ratios}
\ee
%\vspace{-2.0cm}
\begin{figure}[hbp]
\centering
\leavevmode
\epsfxsize=4.8cm
\epsffile{fig150.eps}
%\vspace{-1cm}
\caption{The curves (a), (b), (c) are for $R$, $R2$, $R1$ respectively.}
\label{fig:R}
\end{figure}

Thus a better way to obtain $\gamma$ would be to use $R_{1}$ when a 
precise value for ${\mathcal B}(\bar B^0\rightarrow \bar K^0 {\pi}^0) $ will be
available. The  central
value of $0.80$ for $R$ corresponds to $\gamma = 110^{\circ}$, close to 
the   value  of $(113^{+25}_{-23})^{\circ}$ found
by the CLEO Collaboration in an analysis of all known charmless two-body
$B$ decays with the factorisation model \cite{cleo1}. It seems that 
the CLEO data favor a large  
$\gamma$ in the range $90^{\circ}-120^{\circ}$. With a large $\gamma$, 
for example, with the central value of $110^{\circ} $, as shown in Fig.4,
the predicted $B \rightarrow K\pi$ branching ratios are larger than that for
$\gamma = 70^{\circ}$ and are closer to the data. The data also show that
 $B^-\rightarrow \bar K^0 \pi^-$ and $\bar B^0 \rightarrow K^-\pi^+$ 
are the two largest modes with near-equal
branching ratios in qualitative agreement with factorisation. However,
for $\gamma = 70^{\circ}$, 
Fig.1 shows that these two largest branching ratios are quite apart,
except for $\delta < 50^{\circ}$ while for $\gamma = 110^{\circ}$,
Fig.4 suggests that these two
branching ratios are closer to each other only for
$\delta $ in the range $40^{\circ}-70^{\circ}$. 
With  $\gamma < 110^{\circ}$ and some adjustment of form 
factors, the current $s$ quark mass and CKM parameters, 
it might be possible to accommodate these two largest
branching ratios with  $\delta < 50^{\circ}$.

The  CP-asymmetries, as shown in \cite{Isola}, for $\gamma = 110^\circ$,    
are in the range $\pm(0.04)$ to $\pm(0.3)$
for the preferred values of $\delta$ in the range
$(40-70)^{\circ}$, but could be smaller for $\delta < 50^{\circ}$. 
The CLEO measurements \cite{cleo3} however, do not show any 
large CP-asymmetry in $B \rightarrow K\pi$ decays, but the errors are still
too large to draw any conclusion at the moment. 

In conclusion, factorisation with enhancement of the strong penguin
contribution  seems to describe qualitatively
the $B \rightarrow K\pi$ decays, 
Further measurements will allow 
a more precise test of factorisation and a determination of
the weak  angle $\gamma $ from the FSI phase-independent
relations shown above.
\vspace{-1cm}
\begin{figure}[hbp]
\centering
\leavevmode
\epsfxsize=4.8cm
\epsffile{fig110.eps}
\vspace{-1cm}
\caption{${\mathcal B}(B\rightarrow K\pi)$ vs. $\delta$ for $\gamma = 110^\circ$.
The curves  
(a), (b), (c), (d)  are for the CP-averaged branching ratios
$B^-\rightarrow K^-\pi^0,\ \bar K^0 \pi^-$ and 
$\bar B^0 \rightarrow K^-\pi^+,\ \bar K^0\pi^0$, respectively.}
\label{fig:BRc}
\end{figure}



I would like to thank S. Narison and the organisers of QCD00 for the
warm hospitality extended to me at Montpellier.





 

 
 












\begin{thebibliography}{9}

\bibitem{cleo2} CLEO Collaboration, D. Cronin-Hennessy {\em et al.}, 
CLEO 99-18, .

\bibitem{Deshpande} N. G. Deshpande, X. G. He, W. S. Hou, and S. Pakvasa,
Phys. Rev. Lett. {\bf 82}  (1999) 2240. 

\bibitem{Ali} A.~Ali, G.~Kramer, and C.~L{\"u}, Phys. Rev. D {\bf 58}  
(1998) 094009.

\bibitem{Neubert2} M. Neubert and B. Stech,  to appear in
{\em Heavy Flavors II}, eds.  A. J. Buras and H. Lindner,  World
Scientific, Singapore.

\bibitem{Beneke} M. Beneke, G. Buchalla, M. Neubert, and 
C. T. Sachrajda,  
Phys. Rev. Lett. {\bf 83} (1999) 1914~; \\ 
M. Beneke, Talk given at the Heavy Flavor 8, 
Southampton, UK, 1999, .

\bibitem{Isola} Claudia Isola and  T. N. Pham,  
to appear in Phys. Rev. D.

\bibitem{Buras} A. Buras, M. Jamin. M. Lautenbacher, and P. Weisz, 
Nucl. Phys. B {\bf 400}, 37 (1993); A. Buras, M. Jamin, and 
M. Lautenbacher, ibid. B {\bf 400} (1993) 75.


\bibitem{Ciuchini} M. Ciuchini, E. Franco, G. Martinelli, L. Reina, and
L. Silvestrini, Phys. Lett. B {\bf 316}  (1993) 127; 
M. Ciuchini, E. Franco, G. Martinelli, and L. Reina,
Nucl. Phys. B {\bf 415}  (1994) 403.

\bibitem{Fleischer1} R. Fleischer, Z. Phys. C {\bf 58} (1993) 483~; 
R. Fleischer and T. Mannel, Phys. Rev. D {\bf 57}  (1998) 2752.

\bibitem{Kramer} G. Kramer, W. Palmer, and H. Simma, 
Nucl. Phys. B {\bf 428} (1994) 77.

\bibitem{Deshpande1} N. G. Deshpande and X. G. He, 
Phys. Lett. B {\bf 336} (1994) 471. 



\bibitem{PDG}  For a review of the CKM model, see, e.g,\\
Particle Data Group, C. Caso {\em et al.},
Eur. Phys. Jour. C{\bf 3} (1998) 103.

\bibitem{Nir}  Y. Nir and H. R. Quinn, Ann. Rev. Nucl. Sci. {\bf 42} 
(1992) 212.


\bibitem{Hou} J.-M. G\' erard and W.S. Hou, 
Phys. Rev. Lett. {\bf 62}  (1989) 855; Phys. Rev. D {\bf 43} (1991) 2909.

\bibitem{Bauer} M. Bauer, B. Stech and M. Wirbel, Z. Phys. C
{\bf 34}  (1987) 103.

\bibitem{Colangelo} P. Colangelo, M. Ladisa, G. Nardulli and  T. N. Pham,
Phys.Lett B {\bf 416} (1998) 208.

\bibitem{cleo1} CLEO Collaboration, Y. Kwon {\em et al.}, 
CLEO-CONF 99-14.

\bibitem{cleo3} CLEO Collaboration, S. Chen {\em et al.}, CLEO 99-17, 
.




\end{thebibliography}


\end{document}

%%% ====================================================================
%%%   This file is freely redistributable and placed into the
%%%   public domain by Tomas Rokicki.
%%%  @TeX-file{
%%%     author          = "Tom Rokicki",
%%%     version         = "2.7",
%%%     date            = "25 October 1996",
%%%     time            = "10:00:05 MDT",
%%%     filename        = "epsf.tex",
%%%     address         = "Tom Rokicki
%%%                        Box 2081
%%%                        Stanford, CA 94309
%%%                        USA",
%%%     telephone       = "+1 415 855 9989",
%%%     email           = "rokicki@cs.stanford.edu (Internet)",
%%%     codetable       = "ISO/ASCII",
%%%     keywords        = "PostScript, TeX",
%%%     supported       = "yes",
%%%     abstract        = "This file contains macros to support the inclusion
%%%                        of Encapsulated PostScript files in TeX documents.",
%%%     docstring       = "This file contains TeX macros to include an
%%%                        Encapsulated PostScript graphic.  It works
%%%                        by finding the bounding box comment,
%%%                        calculating the correct scale values, and
%%%                        inserting a vbox of the appropriate size at
%%%                        the current position in the TeX document.
%%%
%%%                        To use, simply say
%%%
%%%                        \input epsf % somewhere early on in your TeX file
%%%
%%%                        % then where you want to insert a vbox for a figure:
%%%                        \epsfbox{filename.ps}
%%%
%%%                        Alternatively, you can supply your own
%%%                        bounding box by
%%%
%%%                        \epsfbox[0 0 30 50]{filename.ps}
%%%
%%%                        This will not read in the file, and will
%%%                        instead use the bounding box you specify.
%%%
%%%                        The effect will be to typeset the figure as
%%%                        a TeX box, at the point of your \epsfbox
%%%                        command. By default, the graphic will have
%%%                        its `natural' width (namely the width of
%%%                        its bounding box, as described in
%%%                        filename.ps). The TeX box will have depth
%%%                        zero.
%%%
%%%                        You can enlarge or reduce the figure by
%%%                        saying
%%%
%%%                          \epsfxsize=<dimen> \epsfbox{filename.ps}
%%%                        or
%%%                          \epsfysize=<dimen> \epsfbox{filename.ps}
%%%
%%%                        instead. Then the width of the TeX box will
%%%                        be \epsfxsize and its height will be scaled
%%%                        proportionately (or the height will be
%%%                        \epsfysize and its width will be scaled
%%%                        proportionately).
%%%
%%%                        The width (and height) is restored to zero
%%%                        after each use, so \epsfxsize or \epsfysize
%%%                        must be specified before EACH use of
%%%                        \epsfbox.
%%%
%%%                        A more general facility for sizing is
%%%                        available by defining the \epsfsize macro.
%%%                        Normally you can redefine this macro to do
%%%                        almost anything.  The first parameter is
%%%                        the natural x size of the PostScript
%%%                        graphic, the second parameter is the
%%%                        natural y size of the PostScript graphic.
%%%                        It must return the xsize to use, or 0 if
%%%                        natural scaling is to be used.  Common uses
%%%                        include:
%%%
%%%                           \epsfxsize  % just leave the old value alone
%%%                           0pt         % use the natural sizes
%%%                           #1          % use the natural sizes
%%%                           \hsize      % scale to full width
%%%                           0.5#1       % scale to 50% of natural size
%%%                           \ifnum #1>\hsize\hsize\else#1\fi
%%%                                       % smaller of natural, hsize
%%%
%%%                        If you want TeX to report the size of the
%%%                        figure (as a message on your terminal when
%%%                        it processes each figure), say
%%%                        `\epsfverbosetrue'.
%%%
%%%                        If you only want to get the bounding box
%%%                        extents, without producing any output boxes
%%%                        or \special{}, then say
%%%                        \epsfgetbb{filename}.  The extents will be
%%%                        saved in the macros \epsfllx \epsflly
%%%                        \epsfurx \epsfury in PostScript units of
%%%                        big points.
%%%
%%%                        Revision history:
%%%
%%%                        ---------------------------------------------
%%%                        epsf.tex macro file:
%%%                        Originally written by Tomas Rokicki of
%%%                        Radical Eye Software, 29 Mar 1989.
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Don Knuth, 3 Jan 1990.
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Tomas Rokicki, 18 Jul 1990.
%%%                        Accept bounding boxes with no space after
%%%                        the colon.
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Nelson H. F. Beebe
%%%                        <beebe@math.utah.edu>, 03 Dec 1991 [2.0].
%%%                        Add version number and date typeout.
%%%
%%%                        Use \immediate\write16 instead of \message
%%%                        to ensure output on new line.
%%%
%%%                        Handle nested EPS files.
%%%
%%%                        Handle %%BoundingBox: (atend) lines.
%%%
%%%                        Do not quit when blank lines are found.
%%%
%%%                        Add a few percents to remove generation of
%%%                        spurious blank space.
%%%
%%%                        Move \special output to
%%%                        \epsfspecial{filename} so that other macro
%%%                        packages can input this one, then change
%%%                        the definition of \epsfspecial to match
%%%                        another DVI driver.
%%%
%%%                        Move size computation to \epsfsetsize which
%%%                        can be called by the user; the verbose
%%%                        output of the bounding box and scaled width
%%%                        and height happens here.
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Nelson H. F. Beebe
%%%                        <beebe@math.utah.edu>, 05 May 1992 [2.1].
%%%                        Wrap \leavevmode\hbox{} around \vbox{} with
%%%                        the \special so that \epsffile{} can be
%%%                        used inside \begin{center}...\end{center}
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Nelson H. F. Beebe
%%%                        <beebe@math.utah.edu>, 09 Dec 1992 [2.2].
%%%                        Introduce \epsfshow{true,false} and
%%%                        \epsfframe{true,false} macros; the latter
%%%                        suppresses the insertion of the PostScript,
%%%                        and instead just creates an empty box,
%%%                        which may be handy for rapid prototyping.
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Nelson H. F. Beebe
%%%                        <beebe@math.utah.edu>, 14 Dec 1992 [2.3].
%%%                        Add \epsfshowfilename{true,false}.  When
%%%                        true, and \epsfshowfalse is specified, the
%%%                        PostScript file name will be displayed
%%%                        centered in the figure box.
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Nelson H. F. Beebe
%%%                        <beebe@math.utah.edu>, 20 June 1993 [2.4].
%%%                        Remove non-zero debug setting of \epsfframemargin,
%%%                        and change margin handling to preserve EPS image
%%%                        size and aspect ratio, so that the actual
%%%                        box is \epsfxsize+\epsfframemargin wide by
%%%                        \epsfysize+\epsfframemargin high.
%%%                        Reduce output of \epsfshowfilenametrue to
%%%                        just the bare file name.
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Nelson H. F. Beebe
%%%                        <beebe@math.utah.edu>, 13 July 1993 [2.5].
%%%                        Add \epsfframethickness for control of
%%%                        \epsfframe frame lines.
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Nelson H. F. Beebe
%%%                        <beebe@math.utah.edu>, 02 July 1996 [2.6]
%%%                        Add missing initialization \epsfatendfalse;
%%%                        the lack of this resulted in the wrong
%%%                        BoundingBox being picked up, mea culpa, sigh...
%%%                        ---------------------------------------------
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Nelson H. F. Beebe
%%%                        <beebe@math.utah.edu>, 25 October 1996 [2.7]
%%%                        Update to match changes in from dvips 5-600
%%%                        distribution: new user-accessible macros:
%%%                        \epsfclipon, \epsfclipoff, \epsfdrafton,
%%%                        \epsfdraftoff, change \empty to \epsfempty.
%%%                        ---------------------------------------------
%%%  }
%%% ====================================================================

 \immediate\write16{This is `epsf.tex' v2.7 <25 October 1996>}%
%
\newread\epsffilein    % file to \read
\newif\ifepsfatend     % need to scan to LAST %%BoundingBox comment?
\newif\ifepsfbbfound   % success?
\newif\ifepsfdraft     % use draft mode?
\newif\ifepsffileok    % continue looking for the bounding box?
\newif\ifepsfframe     % frame the bounding box?
\newif\ifepsfshow      % show PostScript file, or just bounding box?
\epsfshowtrue          % default is to display PostScript file
\newif\ifepsfshowfilename % show the file name if \epsfshowfalse specified?
\newif\ifepsfverbose   % report what you're making?
\newdimen\epsfframemargin % margin between box and frame
\newdimen\epsfframethickness % thickness of frame rules
\newdimen\epsfrsize    % vertical size before scaling
\newdimen\epsftmp      % register for arithmetic manipulation
\newdimen\epsftsize    % horizontal size before scaling
\newdimen\epsfxsize    % horizontal size after scaling
\newdimen\epsfysize    % vertical size after scaling
\newdimen\pspoints     % conversion factor
%
\pspoints = 1bp        % Adobe points are `big'
\epsfxsize = 0pt       % default value, means `use natural size'
\epsfysize = 0pt       % ditto
\epsfframemargin = 0pt % default value: frame box flush around picture
\epsfframethickness = 0.4pt % TeX's default rule thickness
%
\def\epsfbox#1{\global\def\epsfllx{72}\global\def\epsflly{72}%
   \global\def\epsfurx{540}\global\def\epsfury{720}%
   \def\lbracket{[}\def\testit{#1}\ifx\testit\lbracket
   \let\next=\epsfgetlitbb\else\let\next=\epsfnormal\fi\next{#1}}%
%
% We use \epsfgetlitbb if the user specified an explicit bounding box,
% and \epsfnormal otherwise.  Because \epsfgetbb can be called
% separately to retrieve the bounding box, we move the verbose
% printing the bounding box extents and size on the terminal to
% \epsfstatus.  Therefore, when the user provided the bounding box,
% \epsfgetbb will not be called, so we must call \epsfsetsize and
% \epsfstatus ourselves.
%
\def\epsfgetlitbb#1#2 #3 #4 #5]#6{%
   \epsfgrab #2 #3 #4 #5 .\\%
   \epsfsetsize
   \epsfstatus{#6}%
   \epsfsetgraph{#6}%
}%
%
\def\epsfnormal#1{%
    \epsfgetbb{#1}%
    \epsfsetgraph{#1}%
}%
%
\def\epsfgetbb#1{%
%
%   The first thing we need to do is to open the
%   PostScript file, if possible.
%
    \openin\epsffilein=#1
    \ifeof\epsffilein
        \errmessage{Could not open file #1, ignoring it}%
    \else                       %process the file
        {%                      %start a group to contain catcode changes
            % Make all special characters, except space, to be of type
            % `other' so we process the file in almost verbatim mode
            % (TeXbook, p. 344).
            \chardef\other=12
            \def\do##1{\catcode`##1=\other}%
            \dospecials
            \catcode`\ =10
            \epsffileoktrue         %true while we are looping
            \epsfatendfalse     %[02-Jul-1996]: add forgotten initialization
            \loop               %reading lines from the EPS file
                \read\epsffilein to \epsffileline
                \ifeof\epsffilein %then no more input
                \epsffileokfalse %so set completion flag
            \else                %otherwise process one line
                \expandafter\epsfaux\epsffileline:. \\%
            \fi
            \ifepsffileok
            \repeat
            \ifepsfbbfound
            \else
                \ifepsfverbose
                    \immediate\write16{No BoundingBox comment found in %
                                    file #1; using defaults}%
                \fi
            \fi
        }%                      %end catcode changes
        \closein\epsffilein
    \fi                         %end of file processing
    \epsfsetsize                %compute size parameters
    \epsfstatus{#1}%
}%
%
% Clipping control:
\def\epsfclipon{\def\epsfclipstring{ clip}}%
\def\epsfclipoff{\def\epsfclipstring{\ifepsfdraft\space clip\fi}}%
\epsfclipoff % default for dvips is OFF
%
% The special that is emitted by \epsfsetgraph comes from this macro.
% It is defined separately to allow easy customization by other
% packages that first \input epsf.tex, then redefine \epsfspecial.
% This macro is invoked in the lower-left corner of a box of the
% width and height determined from the arguments to \epsffile, or
% from the %%BoundingBox in the EPS file itself.
%
% This version is for dvips:
\def\epsfspecial#1{%
     \epsftmp=10\epsfxsize
     \divide\epsftmp\pspoints
     \ifnum\epsfrsize=0\relax
       \special{PSfile=\ifepsfdraft psdraft.ps\else#1\fi\space
		llx=\epsfllx\space
		lly=\epsflly\space
		urx=\epsfurx\space
		ury=\epsfury\space
		rwi=\number\epsftmp
		\epsfclipstring
               }%
     \else
       \epsfrsize=10\epsfysize
       \divide\epsfrsize\pspoints
       \special{PSfile=\ifepsfdraft psdraft.ps\else#1\fi\space
		llx=\epsfllx\space
		lly=\epsflly\space
		urx=\epsfurx\space
		ury=\epsfury\space
		rwi=\number\epsftmp
		rhi=\number\epsfrsize
		\epsfclipstring
               }%
     \fi
}%
%
% \epsfframe macro adapted from the TeXbook, exercise 21.3, p. 223, 331.
% but modified to set the box width to the natural width, rather
% than the line width, and to include space for margins and rules
\def\epsfframe#1%
{%
  \leavevmode                   % so we can put this inside
                                % a centered environment
  \setbox0 = \hbox{#1}%
  \dimen0 = \wd0                                % natural width of argument
  \advance \dimen0 by 2\epsfframemargin         % plus width of 2 margins
  \advance \dimen0 by 2\epsfframethickness      % plus width of 2 rule lines
  \vbox
  {%
    \hrule height \epsfframethickness depth 0pt
    \hbox to \dimen0
    {%
      \hss
      \vrule width \epsfframethickness
      \kern \epsfframemargin
      \vbox {\kern \epsfframemargin \box0 \kern \epsfframemargin }%
      \kern \epsfframemargin
      \vrule width \epsfframethickness
      \hss
    }% end hbox
    \hrule height 0pt depth \epsfframethickness
  }% end vbox
}%
%
\def\epsfsetgraph#1%
{%
   %
   % Make the vbox and stick in a \special that the DVI driver can
   % parse.  \vfil and \hfil are used to place the \special origin at
   % the lower-left corner of the vbox.  \epsfspecial can be redefined
   % to produce alternate \special syntaxes.
   %
   \leavevmode
   \hbox{% so we can put this in \begin{center}...\end{center}
     \ifepsfframe\expandafter\epsfframe\fi
     {\vbox to\epsfysize
     {%
        \ifepsfshow
            % output \special{} at lower-left corner of figure box
            \vfil
            \hbox to \epsfxsize{\epsfspecial{#1}\hfil}%
        \else
            \vfil
            \hbox to\epsfxsize{%
               \hss
               \ifepsfshowfilename
               {%
                  \epsfframemargin=3pt % local change of margin
                  \epsfframe{{\tt #1}}%
               }%
               \fi
               \hss
            }%
            \vfil
        \fi
     }%
   }}%
   %
   % Reset \epsfxsize and \epsfysize, as documented above.
   %
   \global\epsfxsize=0pt
   \global\epsfysize=0pt
}%
%
%   Now we have to calculate the scale and offset values to use.
%   First we compute the natural sizes.
%
\def\epsfsetsize
{%
   \epsfrsize=\epsfury\pspoints
   \advance\epsfrsize by-\epsflly\pspoints
   \epsftsize=\epsfurx\pspoints
   \advance\epsftsize by-\epsfllx\pspoints
%
%   If `epsfxsize' is 0, we default to the natural size of the picture.
%   Otherwise we scale the graph to be \epsfxsize wide.
%
   \epsfxsize=\epsfsize{\epsftsize}{\epsfrsize}%
   \ifnum \epsfxsize=0
      \ifnum \epsfysize=0
	\epsfxsize=\epsftsize
        \epsfysize=\epsfrsize
	\epsfrsize=0pt
%
%   We have a sticky problem here:  TeX doesn't do floating point arithmetic!
%   Our goal is to compute y = rx/t. The following loop does this reasonably
%   fast, with an error of at most about 16 sp (about 1/4000 pt).
%
      \else
	\epsftmp=\epsftsize \divide\epsftmp\epsfrsize
	\epsfxsize=\epsfysize \multiply\epsfxsize\epsftmp
	\multiply\epsftmp\epsfrsize \advance\epsftsize-\epsftmp
	\epsftmp=\epsfysize
	\loop \advance\epsftsize\epsftsize \divide\epsftmp 2
	\ifnum \epsftmp>0
	   \ifnum \epsftsize<\epsfrsize
           \else
	      \advance\epsftsize-\epsfrsize \advance\epsfxsize\epsftmp
           \fi
	\repeat
	\epsfrsize=0pt
      \fi
   \else
     \ifnum \epsfysize=0
       \epsftmp=\epsfrsize \divide\epsftmp\epsftsize
       \epsfysize=\epsfxsize \multiply\epsfysize\epsftmp
       \multiply\epsftmp\epsftsize \advance\epsfrsize-\epsftmp
       \epsftmp=\epsfxsize
       \loop \advance\epsfrsize\epsfrsize \divide\epsftmp 2
       \ifnum \epsftmp>0
	  \ifnum \epsfrsize<\epsftsize
          \else
	     \advance\epsfrsize-\epsftsize \advance\epsfysize\epsftmp
          \fi
       \repeat
       \epsfrsize=0pt
     \else
       \epsfrsize=\epsfysize
     \fi
   \fi
}%
%
% Issue some status messages if the user requested them
%
\def\epsfstatus#1{% arg = filename
   \ifepsfverbose
     \immediate\write16{#1: BoundingBox:
                  llx = \epsfllx\space lly = \epsflly\space
                  urx = \epsfurx\space ury = \epsfury\space}%
     \immediate\write16{#1: scaled width = \the\epsfxsize\space
                  scaled height = \the\epsfysize}%
   \fi
}%
%
%   We still need to define the tricky \epsfaux macro. This requires
%   a couple of magic constants for comparison purposes.
%
{\catcode`\%=12 \global\let\epsfpercent=%\global\def\epsfbblit{%BoundingBox}}%
\global\def\epsfatend{(atend)}%
%
%   So we're ready to check for `%BoundingBox:' and to grab the
%   values if they are found.
%
%   If we find a line
%
%   %%BoundingBox: (atend)
%
%   then we ignore it, but set a flag to force parsing all of the
%   file, so the last %%BoundingBox parsed will be the one used.  This
%   is necessary, because EPS files can themselves contain other EPS
%   files with their own %%BoundingBox comments.
%
%   If we find a line
%
%   %%BoundingBox: llx lly urx ury
%
%   then we save the 4 values in \epsfllx, \epsflly, \epsfurx, \epsfury.
%   Then, if we have not previously parsed an (atend), we flag completion
%   and can stop reading the file.  Otherwise, we must keep on reading
%   to end of file so that we find the values on the LAST %%BoundingBox.
\long\def\epsfaux#1#2:#3\\%
{%
   \def\testit{#2}%             % save second character up to just before colon
   \ifx#1\epsfpercent           % then first char is percent (quick test)
       \ifx\testit\epsfbblit    % then (slow test) we have %%BoundingBox
            \epsfgrab #3 . . . \\%
            \ifx\epsfllx\epsfatend % then ignore %%BoundingBox: (atend)
                \global\epsfatendtrue
            \else               % else found %%BoundingBox: llx lly urx ury
                \ifepsfatend    % then keep parsing ALL %%BoundingBox lines
                \else           % else stop after first one parsed
                    \epsffileokfalse
                \fi
                \global\epsfbbfoundtrue
            \fi
       \fi
   \fi
}%
%
%   Here we grab the values and stuff them in the appropriate definitions.
%
\def\epsfempty{}%
\def\epsfgrab #1 #2 #3 #4 #5\\{%
   \global\def\epsfllx{#1}\ifx\epsfllx\epsfempty
      \epsfgrab #2 #3 #4 #5 .\\\else
   \global\def\epsflly{#2}%
   \global\def\epsfurx{#3}\global\def\epsfury{#4}\fi
}%
%
%   We default the epsfsize macro.
%
\def\epsfsize#1#2{\epsfxsize}%
%
%   Finally, another definition for compatibility with older macros.
%
\let\epsffile=\epsfbox
\endinput

