# use this to histogram calo_hits.dat - number vs time
max=52.e6 #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 = 'fort.112'
set title filename
e=104000000./2
plot filename u (hist(($1),width)):(1.0) smooth freq w boxes lc rgb"green" t 'Energy
pause -1

plot filename u (hist(($6),width)):(1.0) smooth freq w boxes lc rgb"green" t 'pz'
pause -1
reset
width = 2./n
plot filename u (hist((cos($2)),width)):($1 > 4.e7 ?  1.0/1:0) smooth freq w boxes lc rgb"green" t 'cos(theta)'
pause -1
plot filename u (hist(($3),width)):(1.0) smooth freq w boxes lc rgb"green" t 'phi'