1 #include "isajet/pilot.h"
2 REAL FUNCTION SSZZF3(SP)
3 C-----------------------------------------------------------------------
4 C SSWZBF: ziss -> zjss f fbar
6 C-----------------------------------------------------------------------
7 #if defined(CERNLIB_IMPNONE)
10 #include "isajet/sssm.inc"
11 #include "isajet/sspar.inc"
12 #include "isajet/sstmp.inc"
15 DOUBLE PRECISION S,PI,M2,M1,MS,E,QS,MUS,BKT
16 DOUBLE PRECISION THZ,TERM,MZ,DFI,Q,MF,QP
27 E=(S+M2**2-M1**2)/2.D0/M2
30 QP=Q*DSQRT(MAX(0.D0,1.D0-4*MF**2/S))
31 MUS=MS**2+S-M1**2-MF**2
32 BKT=(MS**2-M1**2-MF**2)*(MS**2-M2**2-MF**2)+THZ*M1*M2*(S-2*MF**2)
33 TERM=DLOG((M2*(E+QP)-MUS)/(M2*(E-QP)-MUS))
34 DFI=-.5D0*M2*E*QP-.5D0*(MS**2-M2**2-S-MF**2)*QP-BKT*TERM/4.D0/M2
35 DFI=DFI*PI**2/2.D0/M2/(S-MZ**2)