gustep.F

Go to the documentation of this file.

      SUBROUTINE gustep                                                         
*                                                                               
*     User routine called at the end of each tracking step           
*                                                               
#include "geant321/gcflag.inc"
#include "geant321/gckine.inc"
#include "geant321/gcking.inc"
#include "geant321/gctrak.inc"
#include "celoss.inc"                                                             
*                                                                               
*                                                                               
* *** Debug event and store tracks for drawing
      IF (idebug.NE.0)   CALL gpcxyz
      IF (idebug.NE.0)   CALL gpgkin            
      IF ((iswit(1).EQ.1).AND.(charge.NE.0.)) CALL gsxyz
      IF  (iswit(1).EQ.2)                     CALL gsxyz                        
*                                                                               
* *** Something generated ?                                                     
      IF(ngkine.GT.0) CALL gsking(0)                                            
*                                                                               
* *** Energy deposited                                                          
      etot = etot + destep
*
* *** count number of steps
      IF (charge.NE.0.) THEN
           nstep1 = nstep1 + 1
       nstrk1 = nstrk1 + 1
      ELSE 
           nstep0 = nstep0 + 1
       nstrk0 = nstrk0 + 1
      ENDIF                                                          
*
* *** count nb of invoked processes
      IF (nmec.GT.0) THEN
        DO im = 1,nmec
          iproc = lmec(im)
          IF (iproc.EQ.21) iproc = 12
          IF (iproc.LE.12) nbcall(iproc) = nbcall(iproc)+1
        ENDDO
      ENDIF
*
* *** step size of primary particle
      IF ((itra.EQ.1).AND.(istak.EQ.0)) CALL hfill(3,step*10,0.,1.)
*                
      END