# use this to histogram calo_hits.dat - number vs time
max=1 #max value microseconds
min=0. #min value
width=(max-min)/n #interval width
#function used to map a value to the intervals
hist(x,width)=width*floor(x/width)+width/2.0
set xrange [min:max]
#set yrange [0:]
#set noytics 
#to put an empty boundary around the
#data inside an autoscaled graph.
#set offset graph 0.05,0.05,0.05,0.0
#set xtics min,(max-min)/10,max
set boxwidth width*0.9
set style fill solid 0.5 #fillstyle
set tics out nomirror
set xlabel "energy"
set ylabel "Frequency"
#set ytics 50
#count and plot
filename = 'muon_decay.dat'
set title filename
e=104000000./2
plot filename u (hist(($10+1),width)):(1.0) smooth freq w boxes lc rgb"green" not