gustep.F

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file gustep.F

      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/gconst.inc"
#include "geant321/gctrak.inc"
#include "geant321/gctmed.inc"
#include "geant321/gcvolu.inc"
#include "calor.inc"
#include "celoss.inc"
#include "histo.inc"
*
      data amev /1.e-3/
*
* *** 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
*
* *** Something generated ?
      IF(ngkine.GT.0) CALL gsking(0)
*
* *** Energy deposited and charged track length
      IF (numed.LE.nbabsor) THEN
        if (destep.NE.0.)      edepabs(numed) = edepabs(numed) + destep
        if (charge*step.NE.0.) trckabs(numed) = trckabs(numed) + step
*       longitudinal profile of Edep
        layer = number(nlevel-1)
  ih = maxabs + numed
  if (histo(ih)) call hfill(ih,0.+layer,0.,destep/(amev*nevent))  
      ENDIF
*
* *** Energy flow : leaving an absorber ?
      IF ((numed.LE.nbabsor).and.(inwvol.eq.2)) THEN
        idnow = nbabsor*(layer-1) + numed
        eflow = gekin
        if (ipart.eq.2) eflow = gekin + 2*emass
*         
        if (max(abs(vect(2)),abs(vect(3))).ge.(0.5*caloryz)) then
    eleaklat(idnow)   = eleaklat(idnow)   + eflow
        elseif (vect(4).ge.0.) then
    enerflow(idnow+1) = enerflow(idnow+1) + eflow
        else
    enerflow(idnow)   = enerflow(idnow)   - eflow   
        endif
      ENDIF      
*
      END