# 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)