ISAJET/isasusy/ssdlam.F

   1 #include "isajet/pilot.h"
   2         DOUBLE PRECISION FUNCTION SSDLAM(A,B,C)
   3 C-----------------------------------------------------------------------
   4 C          Kinematic function
   5 C-----------------------------------------------------------------------
   6 #if defined(CERNLIB_IMPNONE)
   7       IMPLICIT NONE
   8 #endif
   9         DOUBLE PRECISION A,B,C
  10 C            Rewrite SSDLAM=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           SSDLAM=(A-B-C)**2-4*B*C
  13         ELSEIF(B.GE.A.AND.B.GE.C) THEN
  14           SSDLAM=(B-A-C)**2-4*A*C
  15         ELSE
  16           SSDLAM=(C-A-B)**2-4*A*B
  17         ENDIF
  18         RETURN
  19         END