SUBROUTINE GUKINE                                                         
*                                                                               
*    Generates Kinematics for primary track
*
*    Data card Kine : Itype  Ekine  xbeam  rbeam                                      
*                                                                               
#include "geant321/gcbank.inc"
#include "geant321/gcflag.inc"
#include "geant321/gckine.inc"
#include "calor.inc"                                                                                                                          
*                                                                               
      DIMENSION VERTEX(3),PLAB(3)
      dimension rndm(3)
*                                                     
      DATA VERTEX/3*0./
      DATA PLAB  /3*0./                                                                                                                  
*
      IF (PKINE(2).NE.0.) vertex(1) = -0.5*BoxSize
*      
*     random in YZ ?
      if ((pkine(3).lt.0.).or.(pkine(3).gt.1.)) pkine(3) = 0.
      rbeam = pkine(3)*0.5*BoxSize
      call GRNDM (rndm,2)            
*
      VERTEX(2) = (2*rndm(1)-1.)*rbeam     
      VERTEX(3) = (2*rndm(2)-1.)*rbeam      
*                                                                  
      CALL GSVERT(VERTEX,0,0,0,0,NVERT)
*
      JPA     = LQ(JPART-IKINE)
      XMASS   =  Q(JPA+7)
      PLAB(1) = SQRT(PKINE(1)*(PKINE(1)+2*XMASS))                                       
*                                                          
      CALL GSKINE(PLAB,IKINE,NVERT,0,0,NT)                           
*                                                                   
* *** Kinematics debug                     
      IF (IEVENT.EQ.1.OR.IDEBUG.NE.0) CALL GPRINT('KINE',0)                                                                     
*                                                                               
      END