GEANT321/geocad/gckpar.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1999/05/18 15:55:17  fca
   6 * AliRoot sources
   7 *
   8 * Revision 1.1.1.1  1995/10/24 10:20:47  cernlib
   9 * Geant
  10 *
  11 *
  12 #include "geant321/pilot.h"
  13 *CMZ :  3.21/02 29/03/94  15.41.32  by  S.Giani
  14 *-- Author :
  15       SUBROUTINE GCKPAR (ISH,NPAR,PAR)
  16 *
  17 #include "geant321/gconsp.inc"
  18 *
  19       DIMENSION PAR(100)
  20 *
  21       IF (ISH.EQ.4) THEN
  22 * ***    trap
  23          NPAR = 11
  24          PH = 0.
  25          IF (PAR(2).NE.0.) PH = ATAN2(PAR(3),PAR(2))*RADDEG
  26          TT = SQRT(PAR(2)**2+PAR(3)**2)
  27          PAR(2) = ATAN(TT)*RADDEG
  28          IF (PH.LT.0.0) PH = PH + 360.0
  29          PAR(3) = PH
  30          PAR(7) = ATAN(PAR(7))*RADDEG
  31          IF (PAR(7).GT.90.0) PAR(7) = PAR(7)-180.0
  32          PAR(11)= ATAN(PAR(11))*RADDEG
  33          IF (PAR(11).GT.90.0) PAR(11) = PAR(11)-180.0
  34 *
  35       ELSE IF (ISH.EQ.10) THEN
  36 * ***    para
  37          PH = 0.
  38          IF (PAR(5).NE.0.) PH = ATAN2(PAR(6),PAR(5))*RADDEG
  39          TT = SQRT(PAR(5)**2+PAR(6)**2)
  40          PAR(4) = ATAN(PAR(4))*RADDEG
  41          IF (PAR(4).GT.90.0) PAR(4) = PAR(4)-180.0
  42          PAR(5) = ATAN(TT)*RADDEG
  43          IF (PH.LT.0.0) PH = PH + 360.0
  44          PAR(6) = PH
  45 *
  46       ELSE IF (ISH.EQ.28) THEN
  47 * ***    gtra
  48          NPAR = 12
  49       ENDIF
  50 *
  51       END