GEANT321/gxuser/gukine.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:21:52  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.35  by  S.Giani
  11 *-- Author :
  12       SUBROUTINE GUKINE
  13 *
  14 ************************************************************************
  15 *                                                                      *
  16 *             Generates Kinematics for primary tracks                  *
  17 *                                                                      *
  18 ************************************************************************
  19 *
  20 #include "geant321/gcflag.inc"
  21 #include "geant321/gckine.inc"
  22 #include "geant321/gconsp.inc"
  23 #include "geant321/gcscan.inc"
  24 #include "geant321/gcomis.inc"
  25 *
  26       DIMENSION VERTEX(3),PLAB(3)
  27       DIMENSION RNDM(2)
  28 *
  29 *     -----------------------------------------------------------------
  30 *
  31 #if defined(CERNLIB_COMIS)
  32       IF(JUKINE.NE.0)THEN
  33          CALL CSJCAL(JUKINE,0,X,X,X,X,X,X,X,X,X,X)
  34          GO TO 99
  35       ENDIF
  36 #endif
  37
  38       VERTEX(1)=VSCAN(1)
  39       VERTEX(2)=VSCAN(2)
  40       VERTEX(3)=VSCAN(3)
  41       IF(IKINE.GT.100)THEN
  42          IK=IKINE-100
  43          THETA=PKINE(2)*DEGRAD
  44          PHI=PKINE(3)*DEGRAD
  45       ELSE
  46          IK=IKINE
  47          CALL GRNDM(RNDM,2)
  48          THETA=PI*RNDM(1)
  49          PHI=TWOPI*RNDM(2)
  50       ENDIF
  51 C
  52       PLAB(1) = PKINE(1)*SIN(THETA)*COS(PHI)
  53       PLAB(2) = PKINE(1)*SIN(THETA)*SIN(PHI)
  54       PLAB(3) = PKINE(1)*COS(THETA)
  55 C
  56       CALL GSVERT(VERTEX,0,0,0,0,NVERT)
  57       CALL GSKINE(PLAB,IK,NVERT,0,0,NT)
  58 *
  59 *              Kinematic debug (controled by ISWIT(1))
  60 *
  61       IF(IDEBUG.EQ.1.AND.ISWIT(1).EQ.1) THEN
  62         CALL GPRINT('VERT',0)
  63         CALL GPRINT('KINE',0)
  64       ENDIF
  65 *
  66   99  END