[u/mrichter/AliRoot.git] / ISAJET / isasusy / sswzf5.F

#include "isajet/pilot.h"
        REAL FUNCTION SSWZF5(SS)
C-----------------------------------------------------------------------
C          SSWZBF: wiss -> zjss f fbar
C          Baer's XI2FUN
C-----------------------------------------------------------------------
#if defined(CERNLIB_IMPNONE)
      IMPLICIT NONE
#endif
#include "isajet/sssm.inc"
#include "isajet/sspar.inc"
#include "isajet/sstmp.inc"
C
      REAL MW,PI,SS
      DOUBLE PRECISION M1,M2,M3,EQ,Q,XMUS,D,XLOG,S
      DATA PI/3.14159265/,MW/80./
C
      S=SS
      M1=TMP(1)
      M2=TMP(2)
      M3=TMP(3)
      MW=AMW
C
      EQ=(S+M1**2-M3**2)/2./M1
      Q=DSQRT(MAX(0.D0,EQ**2-S))
      XMUS=M2**2+S-M3**2
      D=(M1*(EQ+Q)-XMUS)/(M1*(EQ-Q)-XMUS)
      XLOG=DLOG(D)
      SSWZF5=PI**2/2./M1*M3*S/4./(SS-MW**2)*XLOG
      RETURN
      END
Commit	Line	Data
0795afa3	1	#include "isajet/pilot.h"
	2	REAL FUNCTION SSWZF5(SS)
	3	C-----------------------------------------------------------------------
	4	C SSWZBF: wiss -> zjss f fbar
	5	C Baer's XI2FUN
	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 MW,PI,SS
	15	DOUBLE PRECISION M1,M2,M3,EQ,Q,XMUS,D,XLOG,S
	16	DATA PI/3.14159265/,MW/80./
	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+M12-M32)/2./M1
	25	Q=DSQRT(MAX(0.D0,EQ**2-S))
	26	XMUS=M22+S-M32
	27	D=(M1(EQ+Q)-XMUS)/(M1(EQ-Q)-XMUS)
	28	XLOG=DLOG(D)
	29	SSWZF5=PI*2/2./M1M3S/4./(SS-MW2)XLOG
	30	RETURN
	31	END