gukine.F

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file gukine.F

      SUBROUTINE gukine
*  
*    Generates Kinematics for primary track
*
*    Data card Kine : Itype  Ekine x0  y0  z0  ux  uy  uz rbeam
*
*    rbeam : size of the random beam extension in fraction of 0.5*calorYZ
*
#include "geant321/gcbank.inc"
#include "geant321/gcflag.inc"
#include "geant321/gckine.inc"
#include "calor.inc"
#include "histo.inc"                                                          
*                                                                               
      dimension vertex(3),plab(3)
      dimension rndm(2)
*
      data amev /1.e-3/      
*
* *** set/reset default kinematic
      IF ((ikine.LT.0).or.(abs(pkine(2)).GT.0.5*worldx)) THEN
          ikine = abs(ikine)
          pkine(2) = -0.5*worldx
          pkine(3) = 0.
          pkine(4) = 0.
          pkine(5) = 1.
          pkine(6) = 0.
          pkine(7) = 0.
      ENDIF
*
*     random in YZ ?
      if ((pkine(8).lt.0.).or.(pkine(8).gt.1.)) pkine(8) = 0.
      rbeam = pkine(8)*0.5*caloryz
      call grndm(rndm,2)
*        
      vertex(1) = pkine(2)
      vertex(2) = pkine(3) + (2*rndm(1)-1.)*rbeam
      vertex(3) = pkine(4) + (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))
      plab(1) = pmom*pkine(5)
      plab(2) = pmom*pkine(6)
      plab(3) = pmom*pkine(7)
*
      CALL gskine(plab,ikine,nvert,0,0,nt)
*
* *** Kinematics debug
      IF (ievent.EQ.1.OR.idebug.NE.0) CALL gprint('KINE',0)
*                                                                               
      END