ISAJET/isasusy/ssxlam.F

   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=A**2+B**2+C**2-2*A*B-2*A*C-2*B*C
  11       IF(A.GE.B.AND.A.GE.C) THEN
  12         SSXLAM=(A-B-C)**2-4*B*C
  13       ELSEIF(B.GE.A.AND.B.GE.C) THEN
  14         SSXLAM=(B-A-C)**2-4*A*C
  15       ELSE
  16         SSXLAM=(C-A-B)**2-4*A*B
  17       ENDIF
  18       RETURN
  19       END