TestEm1/geant3/src/uglast.F

Go to the documentation of this file.

      SUBROUTINE uglast                                                         
*                                                                               
*     Termination routine to print histograms and statistics                    
*                                                                                                                              
#include "geant321/gcflag.inc"
#include "geant321/gctrak.inc"
#include "celoss.inc"        
*
* *** nb of steps and tracks per event
      totevt = ievent
      fntrk0 = ntrak0/totevt
      fntrk1 = ntrak1/totevt
      fnstp0 = nstep0/totevt
      fnstp1 = nstep1/totevt
      print 750, fntrk0,fntrk1,fnstp0,fnstp1
*
* *** frequency of processes call
      CALL uctoh('MUNU',namec(12),4,4)
      print 760,(namec(i),i=1,12)
      print 761,(nbcall(i)/totevt,i=1,12)      
*                                                                                                                        
* *** geant termination                                                                               
      CALL glast
*
* *** print  selected histo                                                      

*
* *** close HIGZ file
      CALL hplend
*                                                                               
* *** Save selected histograms                                                  
      IF (iswit(2).EQ.1) THEN                                                   
         CALL hrput(0,'testem1.hbook','N')                                                                        
      ENDIF                                                                     
*
* *** formats            
  750 FORMAT(//,1x,'nb tracks/event neutral: ',f9.2,'  charged: ',f9.2,
     +        /,1x,'nb steps/event  neutral: ',f9.2,'  charged: ',f9.2)
     
  760 FORMAT(//,1x,'nb of process calls per event: ',
     +        /,3x,12a12)
  761 FORMAT(   3x,12f12.4,//)
*
      END