gukine.F

Go to the documentation of this file.

      SUBROUTINE gukine
*
*    Generates Kinematics for primary track
*
*    Data card Kine : Itype  Ekine
*                                                                               
#include "geant321/gcbank.inc"
#include "geant321/gcflag.inc"
#include "geant321/gckine.inc"
*
#include "detector.inc"
*
      dimension vertex(3),plab(3)
      dimension rndm(2)
* 
      DATA vertex/3*0./
      DATA plab  /3*0./                                                                                                                  
*
      vertex(1) = -0.5*boxsize
*      
*     random in YZ 
      beam = 0.4*boxsize
      call grndm(rndm,2)
*
      vertex(2) = (2*rndm(1)-1.)*beam     
      vertex(3) = (2*rndm(2)-1.)*beam      
*
      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