# use this to histogram a column - especially distribution of scattering angles
width=44.686/24.
max=24
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)/6,max
set mxtics 4
set boxwidth 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
set title filename
plot filename u (hist(($4+1.533+.1)/width,1)):(1.0) smooth freq w boxes lc rgb"green" t "Calorimeter"