geant3/src/gutrak.F

0001       SUBROUTINE GUTRAK
0002 *
0003 *     User routine to control tracking of one track
0004 *     Called by GTREVE
0005 *
0006 #include "geant321/gckine.inc"
0007 #include "geant321/gconst.inc"
0008 #include "geant321/gctrak.inc"
0009 #include "geant321/gctmed.inc"
0010 #include "geant321/gcvolu.inc"
0011 #include "calor.inc"
0012 #include "celoss.inc"
0013 *
0014 * *** initialisation of Eflow for primary track
0015 *
0016       if ((itra.eq.1).and.(istak.eq.0)) then
0017 *       unique Id for absorber+layer
0018         Idnow = 1
0019         if (numed.le.NbAbsor) then
0020           layer = number(nlevel-1)
0021           Idnow = NbAbsor*(layer-1) + numed
0022         endif
0023 *
0024         Eflow = gekin
0025         if (ipart.eq.2) Eflow = gekin + 2*emass
0026 *
0027         do k=1,Idnow
0028           EnerFlow(Idnow) = EnerFlow(Idnow) + Eflow
0029         enddo
0030       endif
0031
0032       CALL GTRACK
0033 *
0034       END