# use this to histogram calo_hits.dat - number vs time
#max=700. #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 "time [us]"
#set ylabel "Frequency"
#set ytics 50
#count and plot
#filename = 'muon_decay.dat'
#set title filename
#Emin = (1.8-3.1)/3.1
#tmin=35.e-6
#tmin=0.
#pmin=1.8/3.1
#plot filename u (hist(($3)*1.e6,width)):(($10+1)**2> pmin**2 && $3>tmin ? 1.0/1:0) smooth freq w boxes lc rgb"green" not
#stats filename u (hist(($3)*1.e6,width)):($10> Emin && $3>tmin ? 1.0/1:0)