uglast.F

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file uglast.F

      SUBROUTINE uglast                                                         
*                                                                               
*     Termination routine to print histograms and statistics                    
*
#include "geant321/gcflag.inc" 
#include "geant321/gctrak.inc"
#include "calor.inc"
#include "celoss.inc"
#include "histo.inc"
*
      dimension edeptot(maxabs)
*      
      CHARACTER*20 matnam
      CHARACTER*4 unit1, unit2
*      
      data amev /1.e-3/      
*
* *** mean energy deposit and track length
      print 749
      print 750,ievent
      print 751 
      fnorm    = 1./ievent
      do k=1,nbabsor
        call gfmate(materabs(k),matnam,dua,duz,dud,dur,dui,udu,idu)      
        aveedep  = fnorm*sumedep(k)
  aveedep2 = fnorm*su2edep(k)
        rmsedep  = sqrt(abs(aveedep2-aveedep**2))
        avetrck  = fnorm*sumtrck(k)
  avetrck2 = fnorm*su2trck(k)
        rmstrck  = sqrt(abs(avetrck2-avetrck**2))
*     
        CALL gevkev(aveedep,edmean,unit1)
        CALL gevkev(rmsedep,edrms ,unit2)
        print 752,k,matnam,edmean,unit1,edrms,unit2,avetrck,rmstrck
      enddo
      print 749
*
* *** Energy flow
      iplmax = nbabsor*nblayer + 1
      do k=1,iplmax
        fk = k
  enerflow(k) = enerflow(k)*fnorm
  eleaklat(k) = eleaklat(k)*fnorm
        ih = 2*maxabs + 1
        if (histo(ih)) call hfill(ih, fk, 0., enerflow(k)/amev)
        ih = 2*maxabs + 2
        if (histo(ih)) call hfill(ih, fk, 0., eleaklat(k)/amev) 
      enddo 
*
* *** Energy deposit from energy flow balance
      do k=1,nbabsor
        edeptot(k)  = 0.
      enddo
*
      idmax = nbabsor*nblayer
      do k=1,idmax
        iabsor = mod(k,nbabsor)
        if (iabsor.eq.0) iabsor = nbabsor
  ebalance = enerflow(k) - enerflow(k+1) - eleaklat(k)
  edeptot(iabsor) = edeptot(iabsor) + ebalance
      enddo
      
      print 761
      do k=1,nbabsor
        call gfmate(materabs(k),matnam,dua,duz,dud,dur,dui,udu,idu)
  CALL gevkev(abs(edeptot(k)),edtot,unit1)        
  if (edeptot(k).lt.0.) edtot = -edtot
        print 762,k,matnam,edtot,unit1
      enddo
      print 749             
*                       
*                                             
* *** geant termination
      CALL glast
*
* *** close HIGZ
      CALL hplend
*                                                                               
* *** Save histo                                                  
      CALL hrput(0,filename,'N')

*      
*
* *** formats
  749 FORMAT(/, 60(1h-),/)            
  750 FORMAT(1x,'Nb of events:',i6)
  751 FORMAT(1x,'AbsNo  Material',5x,'Energy deposit',14x,'Trck length')
  752 FORMAT(i5,3x,a10,f8.3,a4,' +- ',f6.2,a4,f9.3,' cm +- ',f6.2,' cm')
  761 FORMAT('  Energy deposition from energy flow balance:',/,/,
     &       1x,'AbsNo  Material',5x,'Total Edep')
  762 FORMAT(i5,3x,a10,f10.3,a4)     
*              
      END