ISAJET/isasusy/sszzf3.F

   1 #include "isajet/pilot.h"
   2         REAL FUNCTION SSZZF3(SP)
   3 C-----------------------------------------------------------------------
   4 C          SSWZBF: ziss -> zjss f fbar
   5 C          Baer's FI
   6 C-----------------------------------------------------------------------
   7 #if defined(CERNLIB_IMPNONE)
   8       IMPLICIT NONE
   9 #endif
  10 #include "isajet/sssm.inc"
  11 #include "isajet/sspar.inc"
  12 #include "isajet/sstmp.inc"
  13 C
  14       REAL SP
  15       DOUBLE PRECISION S,PI,M2,M1,MS,E,QS,MUS,BKT
  16       DOUBLE PRECISION THZ,TERM,MZ,DFI,Q,MF,QP
  17       DATA PI/3.14159265D0/
  18 C
  19       MZ=AMZ
  20       M2=TMP(1)
  21       M1=TMP(2)
  22       MS=TMP(3)
  23       THZ=TMP(4)
  24       MF=TMP(6)
  25 C
  26       S=SP
  27       E=(S+M2**2-M1**2)/2.D0/M2
  28       QS=E**2-S
  29       Q=DSQRT(MAX(0.D0,QS))
  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)
  36       SSZZF3=DFI
  37       RETURN
  38       END