set terminal postscript  eps color enhanced solid font "Helvetica,12" size 9cm, 4cm
#set terminal postscript eps color solid enhanced font "Helvetica,10" size 8cm,10cm
set border 3
set bars 0
set isosamples 40,40
set sample 1000,1000
set output 'plot_phase.eps'
set multiplot
#set grid
#unset xlabel
#unset label
#set format x ""
f(x)=(x>=0?bg+((a1**2)*exp(-2*(x-x0)/t1))+((a2**2)*exp(-2*(x-x0)/t2))+(2*a1*a2*exp(-(x-x0)*((1/t1)+(1/t2)))*cos((2*3.14*0.266*(x-x0))+p)):1)
set origin 0.01,0.035
set size 1,0.51
set ytics 200 nomirror
set xtics 5 nomirror
x0=0;a1=30;a2=60;t1=10.3;t2=26;p=-1;bg=19
fit [1:] f(x) 'pump_HV_out_LL_processed.dat' u ($1-459):2 via a1,a2,p,bg
set label font "Helvetica,11" "Coincidences" at first -14.5, first 510 rotate by 90
set label font "Helvetica,11" "{/Symbol D}t_{is} (ns)" at first 12, first -180
set label font "Helvetica,12" "signal: {/Symbol \367} L {/Symbol \361}" at first 15, first 450
set label font "Helvetica,12" "idler: {/Symbol \367} L {/Symbol \361}" at first 15, first 330
plot [-10:40] 'pump_HV_out_LL_processed.dat' u ($1-459):2:(sqrt($2))w yerrorbars notitle pt 7 ps 0.5 lc rgb "red", f(x) lc rgb "blue" notitle

unset xlabel
unset label
set format x ""
set size 1,0.5
set origin 0.01,0.48
x0=0;a1=300;a2=300;t1=10.3;t2=26;p=0.1;bg=13
fit [2:] f(x) 'pump_HV_out_VH_processed.dat' u ($1-459.2):2 via a1,a2,p,bg
set label font "Helvetica,12" "signal: {/Symbol \367} V {/Symbol \361}" at first 15, first 450
set label font "Helvetica,12" "idler:{/Symbol \367} H {/Symbol \361}" at first 15, first 330
plot [-10:40]'pump_HV_out_VH_processed.dat' u ($1-459.2):2:(sqrt($2))w yerrorbars notitle pt 7 ps 0.5 lc rgb "red", f(x) lc rgb "blue"  notitle
#set origin 0.01,0.5
#set size 1,0.5


#set format y "%4.0f"
#set format x "%2.0f"
#set label font "Helvetica,12" "Coincidence /s" at first -14.5, first 400 rotate by 90
#set xlabel "Time delay (ns)" font "Helvetica,18”
#f(x)=(x>=1?bg+((a1**2)*exp(-2*(x-x0)/t1))+((a2**2)*exp(-2*(x-x0)/t2))+(2*a1*a2*exp(-(x-#x0)*((1/t1)+(1/t2)))*cos((2*3.14*0.266*(x-x0))+p)):1)
#x0=1;a1=1000;a2=0.001;t1=15;t2=30;p=-10;bg=10
#fit [2:] f(x) 'pump_HV_out_HV.dat' u 1:2 via a1,a2,p
#plot [-10:40][0:550] 'pump_HV_out_HV.dat' u 1:2:(sqrt($2))w yerrorbars notitle pt 7 ps #0.5,f(x) w l lc rgb "black" lw 2 lt 1 notitle
unset multiplot