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 "calor.inc"
#include "celoss.inc"
#include "histo.inc"
*      
      CHARACTER*20 parnam,matnam
      CHARACTER*4  unit1,unit2
*
* *** run summary
      print 749
      call gfpart(ikine,parnam,ityp,dum,duc,dut,udu,idu)      
      call gfmate(imat ,matnam,dua,duz,dud,dur,dui,udu,idu)
      CALL gevkev(pkine(1),ekine,unit1)
      print 750, ievent,parnam,ekine,unit1,sizex,matnam
      print 749
*
      fnorm    = 1./ievent      
*
* *** total energy deposit
      etot  = fnorm*etot
      etot2 = fnorm*etot2
      rms   = sqrt(abs(etot2-etot**2))
*
      CALL gevkev(etot,avee,unit1)
      CALL gevkev(rms ,rmse ,unit2)
      print 752, avee,unit1,rmse,unit2
*
* *** track length of primary track
      trlen  = fnorm*trlen
      trlen2 = fnorm*trlen2
      rms    = sqrt(abs(trlen2-trlen**2))            
      print 753, trlen,rms
*
* *** compare with csda range from table
      call gftmat(imat,ikine,'RANG',1,pkine(1),csdarang,dum,ixst)
      print 763, csdarang      
*     
* *** projected range of primary track
      xtrlen  = fnorm*xtrlen
      xtrlen2 = fnorm*xtrlen2
      rms     = sqrt(abs(xtrlen2-xtrlen**2))            
      print 754, xtrlen,rms
*     
* *** nb of steps and step size of primary track
      fnbstep  = fnorm*nbstep
      fnbstep2 = fnorm*nbstep2
      rms1     = sqrt(abs(fnbstep2-fnbstep**2))
*      
      steplen  = fnorm*steplen
      steplen2 = fnorm*steplen2
      rms2     = sqrt(abs(steplen2-steplen**2))                  
      print 755, fnbstep,rms1,steplen,rms2
*
* *** transmission coefficients
      absorbed  = fnorm*100*kcoef(1)
      transmit  = fnorm*100*kcoef(2)
      reflected = fnorm*100*kcoef(3)
      print 756, absorbed, transmit, reflected
*                  
      print 749
*                                             
* *** geant termination
      CALL glast
*
* *** close HIGZ
      CALL hplend
*                                                                               
* *** Save histo
      if (iswit(2).eq.1) then
        write(6,99) filename
99      format(/,15x,'histogram file --> Name: ',a25)                                                  
        CALL hrput(0,filename,'N')
      endif      
*
* *** formats
  749 FORMAT(/, 60(1h-),/)            
  750 FORMAT(  1x,'The run is',i8,1x,a10,' of',f8.2,a4,' through ',
     &             e12.4,' cm of ',a10)
  752 FORMAT(  1x,'Total Energy deposited        = ',
     &             f8.3, a4,' +- ', f8.3, a4)
  753 FORMAT(/,1x,'Track length of primary track = ',
     &             e12.5,' cm +- ', e12.5,' cm')
  763 FORMAT(  1x,'csda Range from table         = ',
     &             e12.5,' cm')     
  754 FORMAT(/,1x,'Projected range               = ',
     &             e12.5,' cm +- ', e12.5,' cm')      
  755 FORMAT(/,1x,'Nb of steps of primary track  = ',
     &             f6.2,' +- ', f6.2,
     &         3x,'Step size = ', e12.5,' cm +- ', e12.5,' cm')
  756 FORMAT(/,1x,'absorbed = ',f6.2,' %',3x,'transmit = ', f6.2,' %',
     &         3x,'reflected = ', f6.2,' %')              
*              
      END