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/gctrak.inc"
#include "celoss.inc" 
#include "histo.inc"
*
* *** per track initialisations    
*
* *** process a track                                                                                     
      CALL gtrack
*
* *** count track length of primary particle
      if (istak.eq.0) then 
        trlen  = trlen + sleng
    trlen2 = trlen2 + sleng*sleng
*   
    xlen    = vect(1) - vert(1)
    xtrlen  = xtrlen  + xlen
    xtrlen2 = xtrlen2 + xlen*xlen
*
        nbstep  = nbstep  + nstep
    nbstep2 = nbstep2 + nstep*nstep
*
        stepsize = sleng/nstep
    steplen  = steplen  + stepsize
    steplen2 = steplen2 + stepsize*stepsize         
      endif
*
* *** absorbed, transmited, reflected ?
      if (istak.eq.0) then 
        i = 1
    if (inwvol.eq.3) then
      i = 2
      if (vect(4).le.0.) i = 3
    endif
        kcoef(i) = kcoef(i) + 1
      endif
*
* *** histo track length
      ih = 0
      if (istak.eq.0) ih = 3
      if (istak*charge.ne.0.) ih = 6
      if ((ih.gt.0).and.(histo(ih)))
     &   call hfill(ih,sleng/histunit(ih),0.,1.)
*
* *** histo projected range of primary
      if (istak.eq.0) then 
        ih = 5
        if (histo(ih)) call hfill(ih,xlen/histunit(ih),0.,1.)
      endif                                                                                             
*                                                                               
      END