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

#include "isajet/pilot.h"
      REAL FUNCTION SSXLAM(A,B,C)
C-----------------------------------------------------------------------
C          Kinematic function
C-----------------------------------------------------------------------
#if defined(CERNLIB_IMPNONE)
      IMPLICIT NONE
#endif
      REAL A,B,C
C          Rewrite SSXLAM=A**2+B**2+C**2-2*A*B-2*A*C-2*B*C
      IF(A.GE.B.AND.A.GE.C) THEN
        SSXLAM=(A-B-C)**2-4*B*C
      ELSEIF(B.GE.A.AND.B.GE.C) THEN
        SSXLAM=(B-A-C)**2-4*A*C
      ELSE
        SSXLAM=(C-A-B)**2-4*A*B
      ENDIF
      RETURN
      END
Commit	Line	Data
0795afa3	1	#include "isajet/pilot.h"
	2	REAL FUNCTION SSXLAM(A,B,C)
	3	C-----------------------------------------------------------------------
	4	C Kinematic function
	5	C-----------------------------------------------------------------------
	6	#if defined(CERNLIB_IMPNONE)
	7	IMPLICIT NONE
	8	#endif
	9	REAL A,B,C
	10	C Rewrite SSXLAM=A2+B2+C*2-2AB-2AC-2B*C
	11	IF(A.GE.B.AND.A.GE.C) THEN
	12	SSXLAM=(A-B-C)*2-4B*C
	13	ELSEIF(B.GE.A.AND.B.GE.C) THEN
	14	SSXLAM=(B-A-C)*2-4A*C
	15	ELSE
	16	SSXLAM=(C-A-B)*2-4A*B
	17	ENDIF
	18	RETURN
	19	END