TestEm14/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/gckine.inc"
#include "geant321/gctrak.inc"
*
#include "detector.inc"  
#include "process.inc"
#include "histo.inc"        
*
      character*20 material, particle
      character*4  unit
      real massMFP, massAC
*
* *** run conditions
      call gfmate(imat,material,dum,dum,density,dum,dum,dum,ndum)
      call gfpart(ikine,particle,ndum,dum,dum,dum,dum,ndum)
      call gevkev(pkine(1),energy,unit)
      print 750, ievent,particle,energy,unit,boxsize,material,density
*      
* *** frequency of processes call
      CALL uctoh('MUNU',namec(12),4,4)
      print 760,(namec(i),i=1,12)
      print 761,(nbcall(i),i=1,12)
*
* *** compute mean free path and related quantities
      avtrak  = sumtrak /nbtot
      avtrak2 = sumtrak2/nbtot
      rms = sqrt(abs(avtrak2 - avtrak*avtrak))
      attecoef = 1./avtrak
*
      massmfp = avtrak*density
      massac  = 1./massmfp
      print 770, avtrak,rms,massmfp,attecoef,massac                              
*                                                                                                                        
* *** geant termination                                                                               
      CALL glast
*
* *** close HIGZ file
      CALL hplend
*                                                                               
* *** histograms                                                                                                  
      CALL hrput(0,filename,'N')
*
* *** formats            
*
  750 FORMAT(/,1x,'The run consists of ',i7,1x,a8,' of ',f7.2,a4,
     +      ' throught',f10.2,' cm of ',a8,' (density: ',f7.3,' g/cm3)')
          
  760 FORMAT(/,1x,'Frequency of process calls: ',
     +       /,1x,12a8)
  761 FORMAT(  1x,12i8,/)
  
  770 FORMAT(/,1x,'MeanFreePath:',f12.5,' cm +- ',f12.5,' cm',
     +          5x,'massic:',f12.5,' g/cm2', 
     +        /,1x,'CrossSection:',f12.5,' cm^-1 ',15x,
     +          5x,'massic:',f12.5,' cm2/g',/)  
*
      END