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/gconst.inc"
#include "geant321/gckine.inc"
#include "geant321/gcking.inc"
#include "geant321/gctmed.inc"
#include "geant321/gctrak.inc"
#include "geant321/gcvolu.inc"
#include "pvolum.inc"
#include "celoss.inc"                                                             
*
*                                                                               
* *** Debug event and strore track for drawing
      IF (idebug.NE.0) CALL gpcxyz
      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                                                          
      IF (destep.GT.0.)THEN                                                     
         radius = sqrt(vect(1)**2 + vect(2)**2)
     nr = 1 + radius/(drx0*x0)
     nl = 1 + (z1 + vect(3))/(dlx0*x0)
         dedr(nr) = dedr(nr) + destep                                                
         dedl(nl) = dedl(nl) + destep                                                                                    
      ENDIF
*
* *** track length and total energy deposit                                                                         
      IF (charge.NE.0.) THEN
         strch  = strch  + step
     edepch = edepch + destep 
      ELSE
         strne  = strne  + step
     edepne = edepne + destep
      ENDIF                                                                                                         
*
* *** energy of particles contributing to edep
      IF ((destep).GT.0.)THEN
         elog = log10((gekin+destep)/pkine(3))
     call hfill(21,elog,0.,destep)
      ENDIF                                                                                      
      END