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