TestEm3/geant3/src/guout.F

Go to the documentation of this file.

      SUBROUTINE guout 
*
*     User routine called at the end of each event
*     
#include "geant321/gcflag.inc"
#include "calor.inc"
#include "celoss.inc"
#include "histo.inc"
*
      CHARACTER*4 unit

*
*     histo: energy deposit per tracking medium
      do k=1,nbabsor
    if (histo(k)) call hfill (k,edepabs(k)/histunit(k),0.,1.)
      enddo
*
*     cumul energy deposited, charged track length and energy leakage
      do k=1,nbabsor
        sumedep(k) = sumedep(k) + edepabs(k)
        sumtrck(k) = sumtrck(k) + trckabs(k)
*
        su2edep(k) = su2edep(k) + edepabs(k)**2
        su2trck(k) = su2trck(k) + trckabs(k)**2
      enddo      
*
* *** print
      IF (iswit(1).NE.0) THEN
        do k=1,nbabsor
          CALL gevkev (edepabs(k),eunit,unit)
          print 751,k,eunit,unit,trckabs(k)
    enddo
    print 749
      ENDIF
*      
#ifndef batch      
*                                                                               
* *** drawing
*                                                          
      IF (iswit(1).NE.0) THEN
         CALL gdhead (110110,'TestEm3',0.)
         CALL gdshow (3)
         CALL gdxyz  (0)
      END IF
#endif                
  749 FORMAT(/) 
  751 FORMAT(1x,'Absorber',i2,' : Energy deposit: ',f8.3,a4,
     &                        '   Track length: '  ,f8.3,' cm')   
*
      END