gutrak.F

Go to the documentation of this file.

      SUBROUTINE gutrak
* 
*     User routine to control tracking of one track
*     Called by GTREVE
*
#include "geant321/gckine.inc"
#include "geant321/gconst.inc"
#include "geant321/gctrak.inc" 
#include "geant321/gctmed.inc"
#include "geant321/gcvolu.inc"
#include "calor.inc" 
#include "celoss.inc" 
*
* *** initialisation of Eflow for primary track
*
      if ((itra.eq.1).and.(istak.eq.0)) then
*       unique Id for absorber+layer
        idnow = 1
        if (numed.le.nbabsor) then
          layer = number(nlevel-1)
      idnow = nbabsor*(layer-1) + numed 
        endif 
*
        eflow = gekin
        if (ipart.eq.2) eflow = gekin + 2*emass
*
        do k=1,idnow
          enerflow(idnow) = enerflow(idnow) + eflow
        enddo
      endif       
                                                                                         
      CALL gtrack
*
      END