#!/bin/csh
#              H B O O K   D E M O N S T R A T I O N
#
echo Copy source code to demo_hbook.f
#
cat >demo_hbook.f <<EOF
      PROGRAM DEMO_HBOOK
      INTEGER PAW_MEMORY_SIZE, NBIN, NBINX, NBINY
      REAL LOW, UPP, LOWX, UPPX, LOWY, UPPY
      PARAMETER (PAW_MEMORY_SIZE=100000)
      COMMON /PAWC/ PAW_SPACE(PAW_MEMORY_SIZE)
*
      CALL HLIMIT(PAW_MEMORY_SIZE)      ! Tell HBOOK the size of work space 
*
      NBIN = 100
      LOW = 0.
      UPP = 1.
*
* Initialize 1D histogram with ID 10 and name 'One dim...'
*
      CALL HBOOK1(10,'One dimensional histo',NBIN,LOW,UPP,0.)
*
* Initialize 2D histogram with ID 100 and name 'Two D ...'
*
      NBINX = 60 
      NBINY = 40
      LOWX  = -1. 
      LOWY= 0.  
      UPPX = 1. 
      UPPY = 2.
      CALL HBOOK2(100,'Two D histo',NBINX,LOWX,UPPX,NBINY,LOWY,UPPY,0.)
*
* Loop to fill histograms entries
*
      DO I=1,10000
*
* Generate 3 random numbers (0,1) using RNDM, a CERN library function
*
        R1 = RNDM(I)
        R2 = RNDM(I+1)
        R3 = RNDM(I+2)
*
* Put one entry to 1D histogram
*
        WEIGHT = 1.
        CALL HFILL(10,R1,0,WEIGHT)
*
* Put one entry to 2D histogram
*
        X = -1. + 2.0*R2
        Y =       2.0*R3
        CALL HFILL (100,X,Y,WEIGHT)
      ENDDO
*
* Now the histogram outputs
*
      CALL HISTDO                         ! On program listing
      CALL HRPUT(0,'histos_for_paw','N')   ! On a file for PAW
*
      STOP
      END
EOF
#
echo Compile and link
#
f77 demo_hbook.f -o demo_hbook.exe \
`cernlib`
#
echo Execute
#
demo_hbook.exe > Program.Listing 
#