TEvtGen/PHOTOS/phoan1.F

   1       FUNCTION PHOAN1(X,Y)
   2 C.----------------------------------------------------------------------
   3 C.
   4 C.    PHOTOS:   PHOton radiation in decays calculation of ANgle '1'
   5 C.
   6 C.    Purpose:  Calculate angle from X and Y
   7 C.
   8 C.    Input Parameters:  X, Y
   9 C.
  10 C.    Output Parameter:  Function value
  11 C.
  12 C.    Author(s):  S. Jadach                       Created at:  01/01/89
  13 C.                B. van Eijk                     Last Update: 02/01/90
  14 C.
  15 C.----------------------------------------------------------------------
  16       IMPLICIT NONE
  17       DOUBLE PRECISION PHOAN1
  18       REAL*8 X,Y
  19       REAL*8 PI,TWOPI
  20       COMMON/PHPICO/PI,TWOPI
  21       IF (ABS(Y).LT.ABS(X)) THEN
  22         PHOAN1=ATAN(ABS(Y/X))
  23         IF (X.LE.0.D0) PHOAN1=PI-PHOAN1
  24       ELSE
  25         PHOAN1=ACOS(X/SQRT(X**2+Y**2))
  26       ENDIF
  27       IF (Y.LT.0.D0) PHOAN1=TWOPI-PHOAN1
  28       RETURN
  29       END