#filename = 'Fiber-015.00ns.dat' generated by analyze_harp from harp_plane_hit.dat and by analyze_fibers_harp from fiber_hits.dat
#time_bin = 10.e-9
min(a,b)=(a< b) ? a:b


time_scale = time_bin/1.e-6
do for [harp =1:4] {
set label 'plotting_scripts/harp_fiber_hits.gnu' at graph 0.02,-0.035 left font 'Verdana,6' noenhance
set terminal pdf enhanced fontscale 0.5 size 6.0in, 6.0in

set xlabel 'Time [micro-s]'

harp_num=sprintf("%i",harp)
tbin=sprintf("%i",time_bin*1.e9)

output_file = 'Fibers'.harp_num.'-'.tbin.'-hits.pdf'
set out output_file
print 'Harp '.harp_num
#set label 'Harp '.harp_num at graph 0.1, 1.04
title = 'harp '.harp_num
set multiplot title title  noenhance


unset ylabel
set xrange [0.1:100]
set xtics 0,20
unset ytics
#set timestamp font 'Verdana,6'
min=3
do for [fiber=-min:min] {
if(fiber == -min) {set size 1,0.18; set origin 0.,0.} else {set size 1,0.14; set origin 0.,((fiber+min-1)*0.14+0.18)}
fib = sprintf("%i",fiber)
print 'harp = ', harp
print 'fiber = ', fiber
plot filename u ($3*time_bin*1.e6):($1 == harp && $2 == fiber ? $5:1/0) w l t 'fiber '.fib
unset timestamp
unset xlabel
unset xtics
unset label
}
unset multiplot
unset terminal
print 'write '.output_file
}