ISAJET/isasusy/sswzf3.F

   1 #include "isajet/pilot.h"
   2         REAL FUNCTION SSWZF3(QQ)
   3 C-----------------------------------------------------------------------
   4 C          SSWZBF: wiss -> zjss f fbar
   5 C          Baer's PHIFUN
   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 QQ
  15       DOUBLE PRECISION M1,M2,M3,T,B,XLOG,SQBKT,Q,PI
  16       DATA PI/3.14159265D0/
  17 C
  18       Q=QQ
  19       M1=TMP(1)
  20       M2=TMP(2)
  21       M3=TMP(3)
  22 C
  23       T=M2**2*(M1-2*Q)-M1*M3**2
  24       B=(M1-2*Q)*(M2**2-2*M1*Q-M3**2)
  25       XLOG=DLOG(T/B)
  26       SQBKT=-Q*(M1**2-M3**2-2*M1*Q)/M1/(M1-2*Q)-
  27      $(2*M1*Q-M2**2+M3**2)*XLOG/2./M1
  28       SSWZF3=.5*PI**2*M1*M3*SQBKT/(M1**2-M2**2-2*M1*Q)
  29       RETURN
  30       END