# use this to histogram a column - especially distribution of scattering angles
#filename = 'DownstreamInflectorScattering.dat'
max=10. #max value
min=-10. #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 "px[mrad]"
set ylabel "Frequency"
#set ytics 50
#count and plot
plot "fort.53" u (hist(($7)*1.e3,width)):(1.0) smooth freq w boxes lc rgb"green" t "px[mrad]",\
"" u (hist(($9)*1.e3,width)):(1.0) smooth freq w boxes lc rgb"blue" t "py[mrad]"