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 "calor.inc"
#include "celoss.inc"
#include "histo.inc"
*
      common/local1/ fnorm
*      
      data amm, amev /1.e-1, 1.e-3/
*
      if (ievent.eq.1) fnorm = amm/(nevent*amev*binwidth(1)*histunit(1))                                                                         
*
* *** 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)
*
      if (destep.le.0.) return                                                   
*                                                                               
* *** Energy deposited per event                                                          
      eevent = eevent + destep
*
* *** histo edep longitudinal profile
      xlen = 0.5*sizex + vect(1)
      ih = 1
      if (histo(ih))
     &   call hfill(ih,xlen/histunit(ih),0.,fnorm*destep)
*
* *** histo step size
      ih = 0
      if (istak.eq.0) ih = 4
      if (istak*charge.ne.0.) ih = 7
      if ((ih.gt.0).and.(histo(ih)))
     &   call hfill(ih,step/histunit(ih),0.,1.)
*               
      END