! questo programma illustra l'utilizzo delle
! routines precedenti.
program histogram
use HistoRoutines
implicit none
real :: rnd
integer :: i
integer :: nBin
real :: xmin, xmax
character(len=128):: fileName, loadfileName, xlabel, ylabel, title

! dichiarazione dell'istogramma
type(Histo) :: h1

xmin = 0.0
xmax = 1.99
nBin = 40
! qui si costruisce l'istogramma, cioè si dichiarano
! il numero di bin, e l'intervallo [xmin, xmax]
call ConstructHisto(h1,nBin,xmin,xmax)


! riempiamo l'istogramma
do i=1, 10000
    call random_number(rnd)
    call FillHisto(h1,-log(rnd))
enddo

fileName = 'exp'
call printHisto(h1,fileName,100)
print*,"underflow: ", h1%underflow
print*,"overflow: ", h1%overflow
loadfileName = 'exp'
xlabel = "x"
ylabel = "counts"
title = "exp. distribution"
call drawHisto( fileName, loadfileName, 100 , xlabel, xmin, xmax, ylabel, 0.0, 1500.0, title)


! per distruggere l'istogramma quando si è finito di
! utilizzarlo
call DestructHisto(h1)





end program