# use this to histogram a column - especially distribution of scattering angles
#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 xtics min,(max-min)/10,max
set boxwidth width*0.9
set style fill solid 0.5 #fillstyle
set tics out nomirror
set xlabel xlab
set ylabel "Frequency"
plot filename u (hist(((column(index_col))),width)):(1.0) smooth freq w boxes lc rgb"green" 
#pause -1
#reset
#plot filename u (hist(((column(index_col))),width)):(1.0) smooth freq w boxes lc rgb"green"