TestEm5/geant3/src/gukine.F

Go to the documentation of this file.

      SUBROUTINE gukine
*  
*    Generates Kinematics for primary track
*
*    Data card Kine : Itype  Ekine theta phi rbeam
*
*    rbeam : size of the random beam extension in fraction of 0.5*calorYZ      
*
#include "geant321/gcbank.inc"
#include "geant321/gcflag.inc"
#include "geant321/gconst.inc"
#include "geant321/gckine.inc"
*
#include "geomate.inc"
#include "runinfo.inc"
*
      dimension vertex(3),plab(3)
      dimension rndm(2)
*
*     random in YZ ?
      if ((pkine(4).lt.0.).or.(pkine(4).gt.1.)) pkine(4) = 0.
      rbeam = pkine(4)*0.5*yzabso
      call grndm (rndm,2)            
*
      vertex(1)=xvert     
      vertex(2)=yvert + (2*rndm(1)-1.)*rbeam     
      vertex(3)=zvert + (2*rndm(2)-1.)*rbeam
*      
      CALL gsvert(vertex,0,0,0,0,nvert)           
*
      jpa     = lq(jpart-ikine)
      xmass   =  q(jpa+7)
      pmom    = sqrt(pkine(1)*(pkine(1)+2*xmass))
      theta   = pkine(2)*degrad
      phi     = pkine(3)*degrad   
C
      plab(1)=pmom*sin(theta)*cos(phi)
      plab(2)=pmom*sin(theta)*sin(phi)
      plab(3)=pmom*cos(theta)
*      
*     store the charge of primary particle in a free location
      pkine(10) = q(jpa+8)         
C
      CALL gskine(plab,ikine,nvert,0,0,nt)
*                                                                   
* *** Kinematics debug                     
      IF (ievent.EQ.1.OR.idebug.NE.0) CALL gprint('KINE',0)            
*
      END