gukine.F

Go to the documentation of this file.

      SUBROUTINE gukine                                                         
*                                                                               
*    Generates Kinematics for primary track
*
*    Data card Kine : Itype  Ekine
*                     (pkine(3) is used internaly to store Etot)                                       
*                                                                               
#include "geant321/gcbank.inc"
#include "geant321/gcflag.inc"
#include "geant321/gckine.inc"
#include "pvolum.inc"                                                                     
*                                                                               
      dimension vertex(3),plab(3)                                               
      DATA vertex/3*0./                                                         
      DATA plab  /3*0./                                                         
*                                                            
      vertex(3) = - z1     + 0.0001                         
      CALL gsvert(vertex,0,0,0,0,nvert)                          
*                                                         
      jpa = lq(jpart-ikine)                              
      xmass = q(jpa+7)                                     
      pkine(3) = xmass + pkine(1)                          
      plab(3)  = sqrt(pkine(1)*(pkine(3)+xmass))                                 
*                                                          
      CALL gskine(plab,ikine,nvert,0,0,nt)                           
*                                                                   
* *** Kinematics debug                      
      IF(ievent.EQ.1.OR.idebug.NE.0) CALL gprint('KINE',0)    
*                                                                               
      END