#!/usr/bin/gnuplot

set terminal postscript eps color enhanced solid font 'Helvetica, 16pt' size 8.6cm, 6cm  dl .5

mean_offset = 410035.24218755

range_Dt = 600
unset bars
set output 'g2.eps'
set xrange [- range_Dt/2: range_Dt/2]
set ytics 1000
set yrange [0: 5700]

# set palette
load 'paired.pal'

set xlabel "{/Symbol t} - ~{/Symbol d&{.}}{.6/Symbol \55 } (ns)" #offset 0, -0.3
set ylabel'Normalized correlation g^{(2)}' #offset 1
# set format y '%.0s%c'

# zero line
set arrow from 0,0 to 0,4800 nohead lw 1 dt '.'

# set key title "Fiber length" samplen 0
unset key

fit_lw = 2
point_sz = .5
arr_d = '-'

fib1_a=-8.416
fib1_b=8.416
fib1_y=5300
set arrow from fib1_a,fib1_y to fib1_b,fib1_y ls 1 lw fit_lw heads
set arrow from fib1_a,0 to fib1_a,fib1_y  ls 1 lw fit_lw dt arr_d nohead
set arrow from fib1_b,0 to fib1_b,fib1_y  ls 1 lw fit_lw dt arr_d nohead
set label ' 1.7m' at fib1_b-10, fib1_y + 150

fib2_a=-33.615
fib2_b=33.615
fib2_y=4100
set arrow from fib2_a,fib2_y to fib2_b,fib2_y ls 3 lw fit_lw heads
set arrow from fib2_a,0 to fib2_a,fib2_y ls 3 lw fit_lw dt arr_d nohead
set arrow from fib2_b,0 to fib2_b,fib2_y ls 3 lw fit_lw dt arr_d nohead
set label ' 6.7m' at fib2_b, fib2_y + 150

fib3_a=-157.682
fib3_b=157.682
fib3_y=3300
set arrow from fib3_a,fib3_y to fib3_b,fib3_y ls 5 lw fit_lw heads
set arrow from fib3_a,0 to fib3_a,fib3_y  ls 5 lw fit_lw dt arr_d nohead
set arrow from fib3_b,0 to fib3_b,fib3_y  ls 5 lw fit_lw dt arr_d nohead
set label '31.7m' at fib3_b - 80, fib3_y + 200


fib4_a=-255.546
fib4_b=255.546
fib4_y=2500
set arrow from fib4_a,fib4_y to fib4_b,fib4_y ls 7 lw fit_lw heads
set arrow from fib4_a,0 to fib4_a,fib4_y  ls 7 lw fit_lw dt arr_d nohead
set arrow from fib4_b,0 to fib4_b,fib4_y  ls 7 lw fit_lw dt arr_d nohead
set label '51.7m' at fib4_b - 80, fib4_y + 200


plot [- range_Dt/2: range_Dt/2]\
	'fibre_00m_data.dat' u (($1-mean_offset)):2 every 10 ls 2 lw 1 pt 6 ps point_sz title " 1.7m", \
	'fibre_05m_data.dat' u (($1-mean_offset)):2 every 10 ls 4 lw 1 pt 6 ps point_sz title " 6.7m", \
	'fibre_30m_data.dat' u (($1-mean_offset)):2 every 10 ls 6 lw 1 pt 6 ps point_sz title "31.7m", \
	'fibre_50m_data.dat' u (($1-mean_offset)):2 every 10 ls 8 lw 1 pt 6 ps point_sz title "51.7m", \