uglast.F

Go to the documentation of this file.

      SUBROUTINE uglast
*
#include "geant321/gcflag.inc"
*
#include "runinfo.inc"
#include "histo.inc"
*
      CHARACTER*4 unitm,unitr
      dimension transmit(2),reflect(2)
*
      CALL glast
*
      IF(ievent.le.0) return
*      
      print 749
      print 750,ievent
*             
      totevent = float(ievent)
*
*     energy deposit      
      fmean = sedep/totevent
      rms2 = sedep2/totevent-fmean**2
      rms = 0.
      if (rms2.gt.0.) rms = sqrt(rms2/totevent)
*
      call gevkev(fmean,emean,unitm)
      call gevkev(rms  ,erms ,unitr)      
      print 752, emean,unitm,erms,unitr           
*
*     track length (charged)      
      fmean = strkch/totevent
      rms2 = strkch2/totevent-fmean**2
      rms = 0.
      if (rms2.gt.0.) rms = sqrt(rms2/totevent)
      print 753, fmean,rms
*
*     track length (neutral)      
      fmean = strkne/totevent
      rms2 = strkne2/totevent-fmean**2
      rms = 0.
      if (rms2.gt.0.) rms = sqrt(rms2/totevent)
      print 754, fmean,rms
*
*     number of steps (charged)      
      fmean = sstpch/totevent
      rms2 = sstpch2/totevent-fmean**2
      rms = 0.
      if (rms2.gt.0.) rms = sqrt(rms2/totevent)
      print 755, fmean,rms
*
*     number of steps (neutral)      
      fmean = sstpne/totevent
      rms2 = sstpne2/totevent-fmean**2
      rms = 0.
      if (rms2.gt.0.) rms = sqrt(rms2/totevent)
      print 756, fmean,rms            
*
*     number of secondaries
      gamma = ngamma/totevent
      elect = nelect/totevent
      posit = nposit/totevent
      print 757, gamma,elect,posit
*
*     transmission & reflection coefficients
      transmit(1) = 100.*nbtransm(1)/totevent
      transmit(2) = 100.*nbtransm(2)/totevent
      reflect(1) = 100.*nbreflec(1)/totevent
      reflect(1) = 100.*nbreflec(2)/totevent                  
      print 761, transmit(2)
      print 762, transmit(1)
      print 763, reflect(2)
      print 764, reflect(1)
*
      print 749      
*
*
* *** histograms
*
      call hplend      
      call hrput(0,filename,'N')
*      
* *** formats
*
  749 FORMAT(/ ,80(1h=),/)            
  750 FORMAT(' Number of events = ', i7)
*  
  752 FORMAT(/,' Total energy deposit in absorber per event = ',
     +       f8.3,a4,' +- ',f6.2,a4)
*     
  753 FORMAT(/,' Total track length (charged) in absorber per event = ',
     +       f8.4,' cm +- ',f8.4, ' cm')
  754 FORMAT(  ' Total track length (neutral) in absorber per event = ',
     +       f8.4,' cm +- ',f8.4, ' cm')
*     
  755 FORMAT(/,' Number of steps (charged) in absorber per event = ',
     +       f7.3,' +- ',f7.3)
  756 FORMAT(  ' Number of steps (neutral) in absorber per event = ',
     +       f7.3,' +- ',f7.3)
*     
  757 FORMAT(/,' Number of secondaries per event : ',
     +  'gammas = ',f8.3, '; electrons = ',f8.3,'; positrons =',f8.3)
*     
  761 FORMAT(/,' Nb of events with the primary particle transmitted = ',
     +       f6.2, ' %')
  762 FORMAT(  ' Nb of events with at least 1 particle transmitted ',
     +       '(same charge as primary) = ',f6.2, ' %')
  763 FORMAT(/,' Nb of events with the primary particle reflected = ',
     +       f6.2, ' %')
  764 FORMAT(  ' Nb of events with at least 1 particle reflected ',
     +       '(same charge as primary) = ',f6.2, ' %')                                                           
*
      END