gukine.F

Go to the documentation of this file.

      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