gukine.F

Go to the documentation of this file.

      SUBROUTINE gukine
*                                                                               
*    Generates Kinematics for primary track
*
*    Data card Kine : Itype  Ekine                                       
*                                                                               
#include "geant321/gconst.inc"
#include "geant321/gcflag.inc"
#include "geant321/gckine.inc"                                                                                                                           
*                                                                               
      dimension vertex(3),plab(3)
      dimension rndm(2)
*                                                     
      DATA vertex/3*0./                                                                                                                  
*
*                                                                  
      CALL gsvert(vertex,0,0,0,0,nvert)
*
      CALL grndm(rndm,2)
      costheta = 2*rndm(1) - 1.
      sintheta = sqrt(1. - costheta**2)
      phi   = twopi*rndm(2)                                
*                                                         
      plab(1)  = pkine(1)*sintheta*cos(phi)
      plab(2)  = pkine(1)*sintheta*sin(phi)
      plab(3)  = pkine(1)*costheta                                       
*                                                          
      CALL gskine(plab,ikine,nvert,0,0,nt)                           
*                                                                   
* *** Kinematics debug
      IF (ievent.EQ.1.OR.idebug.NE.0) CALL gprint('KINE',0)
*                                                                               
      END