uhinit.F

Go to the documentation of this file.

      SUBROUTINE uhinit
*
*     To book the user's histograms
*
#include "histo.inc"
#include "calor.inc"
*
      character*50 title
*
      if (idhist.gt.maxhist) return
*      
      if (histo(idhist)) call hdelet(idhist)      
*
      if (idhist.lt.maxabs)
     +  title = 'Edep in absorber '//char(ichar('0')+(idhist))
      if ((idhist.gt.maxabs).and.(idhist.lt.2*maxabs))
     +  title = 'Edep long. profile in absorber '
     +          //char(ichar('0')+(idhist-maxabs))//' (MeV/event)'
      if (idhist.eq.(2*maxabs+1))
     +  title = 'energy flow (MeV/event)'
      if (idhist.eq.(2*maxabs+2))
     +  title = 'lateral energy leak (MeV/event)'               
*      
      call hbook1(idhist,title,nbbins,valmin,valmax,0.)
*      
      histo(idhist) = .true.
      binwidth(idhist) = (valmax-valmin)/nbbins
      if (valunit.le.0.) valunit = 1.                
      histunit(idhist) = valunit                   
*                                                                               
      END