TestEm5/geant3/src/gustep.F File Reference

#include "geant321/gcflag.inc"
#include "geant321/gckine.inc"
#include "geant321/gcking.inc"
#include "geant321/gconst.inc"
#include "geant321/gctmed.inc"
#include "geant321/gctrak.inc"
#include "runinfo.inc"
#include "histo.inc"

Include dependency graph for TestEm5/geant3/src/gustep.F:

Go to the source code of this file.

Functions/Subroutines
subroutine	gustep

Function/Subroutine Documentation

gustep()

subroutine gustep

(

)

Definition at line 3 of file TestEm5/geant3/src/gustep.F.

*
#include "geant321/gcflag.inc"
#include "geant321/gckine.inc"
#include "geant321/gcking.inc"
#include "geant321/gconst.inc"
#include "geant321/gctmed.inc"
#include "geant321/gctrak.inc"
*
#include "runinfo.inc"
#include "histo.inc"
*                                                                               
*                                                                               
* *** Debug event and store tracks 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                        
*
* *** collect informations in the absorber only
*
      if (numed.ne.1) return
*
*     energy deposit
      edep = edep + destep
*
*     number of steps and total track length
      if (charge.ne.0.) then
        stpch = stpch + 1.
    trkch = trkch + step
      else
        stpne = stpne + 1.
    trkne = trkne + step
      endif                     
*
*     manage secondaries
      if (ngkine.gt.0) then
        do k = 1,ngkine
      energy = gkin(4,k)
          id     = gkin(5,k)
          if (id.eq.1) then
        ngamma = ngamma + 1
        if (histo(3)) call hfill(3,log10(1000*energy),0.,1.)
      else if (id.eq.2) then
        nposit = nposit + 1
        energy = energy - emass
        if (histo(2)) call hfill(2,energy/histunit(2),0.,1.)
      else if (id.eq.3) then
        nelect = nelect + 1
        energy = energy - emass
        if (histo(2)) call hfill(2,energy/histunit(2),0.,1.)
      endif         
    enddo
*       keep or kill all secondaries    
    if (.not.kill) call gsking(0)
      endif 
*