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"     
                                                                    
*                                                                                                                                         
*
* *** initialisation track per track
      nstrk0 = 0
      nstrk1 = 0
*                                                                                     
      CALL gtrack
*
* *** count nb of tracks
      IF (charge.NE.0.) THEN
         ntrak1 = ntrak1 + 1
      ELSE 
         ntrak0 = ntrak0 + 1
      ENDIF
*
* *** track length of primary particle
      IF ((itra.EQ.1).AND.(istak.EQ.0)) THEN
     CALL hfill(1,sleng*10,0.,1.)
     nstrk = nstrk0
     if (charge.NE.0.) nstrk = nstrk1
         CALL hfill(2,REAL(nstrk),0.,1.)
      ENDIF                                                                                     
*                                                                               
      END