uhinit.F

Go to the documentation of this file.

      SUBROUTINE uhinit                                                         
*                                                                               
*     To book the user's histograms                                             
*                                                                               
#include "geant321/gckine.inc"
#include "pvolum.inc"                                                                     
*                                                                               
* *** Histograms for showers development                                                                                                                                                                                                                  
*                                                                                                                                         
      edmin=0.                                                              
      edmax=110.                                                             
      trkmx=110.*pkine(1)                                                                                                                         
*                                                                               
      CALL hbook1(1,'total energy deposition (in percent of E inc)'             
     *,110,edmin,edmax,0.0)                                                     
      CALL hbook1(2,'total charged tracklengh (in radl)'                       
     *,110, 0. , trkmx, 0.0)
      CALL hbook1(3,'total neutral tracklengh (in radl)'                       
     *,110, 0. , 10*trkmx, 0.0)                                                                                        
*                                                                               
      CALL hidopt(0,'STAT')                                                     
*
* *** Longitudinal profile
      zmax=nltot*dlx0
                                                                                           
      CALL hbprof(4,'longit energy profile    (in percent of E inc)'
     *, nltot, 0.,zmax, 0., 1000.,' ')
      CALL hbprof(5,'rms longit energy profile (in percent of E inc)'
     *, nltot, 0.,zmax, 0., 1000.,' ')
*      
      zmin = 0.5*dlx0
      zmax = zmin + nltot*dlx0
                                                                                                                          
      CALL hbprof(6,'cumul longit energy dep. (in percent of E inc)'          
     *, nltot, zmin,zmax, 0., 100.,' ')
      CALL hbook1(7,'resolution: cumul L energy dep. (% of E inc)'          
     *, nltot, zmin,zmax, 0.0)                                  
*
* *** Radial profile
      rmax=nrtot*drx0
      
      CALL hbprof(8,'radial energy profile    (in percent of E inc)'          
     *, nrtot, 0.,rmax, 0., 1000.,' ')
      CALL hbprof(9,'rms radial energy profile (in percent of E inc)'          
     *, nrtot, 0.,rmax, 0., 1000.,' ')     
*
      rmin = 0.5*drx0
      rmax = rmin + nrtot*drx0
      
      CALL hbprof(10,'cumul radial energy dep. (in percent of E inc)'          
     *, nrtot, rmin,rmax, 0., 100.,' ')
      CALL hbook1(11,'resolution: cumul R energy dep. (% of E inc)'          
     *, nrtot, rmin,rmax, 0.0)
     
*
* *** Origin of the deposited energy
      tminlog = -6.
      tmaxlog =  0.
      CALL hbook1(21,'edep (normalized 100 percent) versus log(ekin/e0)'          
     *, 100, tminlog, tmaxlog, 0.0)
      CALL hbook1(22,'cumul edep (100 percent) versus log(ekin/e0)'          
     *, 100, tminlog, tmaxlog, 0.0)
                                                                                                                                                     
*                                                                                                                                     
      END