ISAJET/isasusy/sswzf4.F

   1 #include "isajet/pilot.h"
   2         REAL FUNCTION SSWZF4(SS)
   3 C-----------------------------------------------------------------------
   4 C          SSWZBF: wiss -> zjss f fbar
   5 C          Baer's XI1FUN
   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 SS,PI
  15       DOUBLE PRECISION M1,M2,M3,EQ,Q,XMUS,XLOG,TERM,S,MW
  16       DATA PI/3.14159265/
  17 C
  18       S=SS
  19       M1=TMP(1)
  20       M2=TMP(2)
  21       M3=TMP(3)
  22       MW=AMW
  23 C
  24       EQ=(S+M1**2-M3**2)/2./M1
  25       IF (EQ**2.GE.S) THEN
  26         Q=DSQRT(EQ**2-S)
  27       ELSE
  28         Q=0.D0
  29       END IF
  30       XMUS=M2**2+S-M3**2
  31       XLOG=DLOG((M1*(EQ+Q)-XMUS)/(M1*(EQ-Q)-XMUS))
  32       TERM=-.5*M1*EQ*Q-.5*(M2**2-M1**2-S)*Q-
  33      $.25/M1*(M2**2-M3**2)*(M2**2-M1**2)*XLOG
  34 C          SS can stay single precision below
  35       SSWZF4=PI**2/2./M1/(SS-MW**2)*TERM
  36       RETURN
  37       END