ISAJET/isasusy/sszzf1.F

   1 #include "isajet/pilot.h"
   2       REAL FUNCTION SSZZF1(X)
   3 C-----------------------------------------------------------------------
   4 C          SSWZBF: ziss -> zjss f fbar
   5 C          Baer's TFUNC
   6 C-----------------------------------------------------------------------
   7 #if defined(CERNLIB_IMPNONE)
   8       IMPLICIT NONE
   9 #endif
  10 #include "isajet/sslun.inc"
  11 #include "isajet/sssm.inc"
  12 #include "isajet/sspar.inc"
  13 #include "isajet/sstmp.inc"
  14 C
  15       REAL X
  16       DOUBLE PRECISION N,S,Q,D,SQBRKT,TERM1,QS
  17       DOUBLE PRECISION MZ2,MZ1,M,DSN
  18 C
  19       MZ2=TMP(1)
  20       MZ1=TMP(2)
  21       M=TMP(3)
  22       DSN=TMP(4)
  23 C
  24       N=MZ1**2/MZ2**2
  25       S=M**2/MZ2**2
  26       Q=X*(1.D0-N)/2.
  27       QS=Q**2
  28       D=(S-2*S*Q-N)/((1.D0-2*Q)*(S-2*Q-N))
  29       IF(D.LE.0.) THEN
  30         WRITE(LOUT,*) 'ERROR IN SSZZF1: D,S,Q,N=',D,S,Q,N
  31         SSZZF1=0
  32         RETURN
  33       END IF
  34       SQBRKT=-Q*(1.D0-2*Q-N)/(1.D0-2*Q)-(2*Q-S+N)/2.D0*DLOG(D)
  35       TERM1=QS*(1.D0-2*Q-N)**2/(1.D0-2*Q-S)**2/(1.D0-2*Q)
  36      $+DSQRT(N)/2.D0/(1.D0-2*Q-S)*SQBRKT*DSN
  37       SSZZF1=(1.D0-N)/2.D0*TERM1
  38       RETURN
  39       END