1 #include "isajet/pilot.h"
2 REAL FUNCTION SSWZF4(SS)
3 C-----------------------------------------------------------------------
4 C SSWZBF: wiss -> 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"
13 C-----------------------------------------------------------------------
15 DOUBLE PRECISION M1,M2,M3,EQ,Q,XMUS,XLOG,TERM,S,MW
24 EQ=(S+M1**2-M3**2)/2./M1
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