reset
# set terminal postscript  eps color enhanced solid font "Helvetica,12" size 8.7cm, 6cm
# set output 'plot_phase.eps'
set sample 1000,1000
set border 3
set nobar
set multiplot
f(x)=(x>=0?bg+((g**2)*(r**2)*exp(-(x-x0)/t1))+((g**2)*exp(-(x-x0)/t2))+(2*(g**2)*r*exp(-(x-x0)*((1/(2*t1))+(1/(2*t2))))*cos((2*pi*0.266*(x-x0))+p)):1)
#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.0,0.63
set size 1,0.29
set ytics 200
set xtics 5
set format x ""
unset xlabel
unset label
set xtics nomirror
set ytics nomirror

#x0=0.1;a1=3.6;a2=0.05;t1=12.3;t2=24.8;p=0.1;bg=10 
x0=0;g=0.1;r=35;t1=6;t2=13.1;p=180;bg=15
#x0=0;g=0.1;r=35;t1=5.8;t2=13;p=180;bg=10

#x0=0;a1=30;a2=1;t1=8;t2=;p=-1;bg=0.4
#458.8,6
set label 1 font "Helvetica,12" "signal: {/Symbol \367} H {/Symbol \361}" at first 25, first 400
set label 2 font "Helvetica,12" "idler: {/Symbol \367} V {/Symbol \361}" at first 25, first 250
fit [0:] f(x) 'out_HV_new_4ns.dat' u ($1-459.2):2 via g,r,p
plot [-10:40][0:600] 'out_HV_new_4ns.dat' u ($1-459.2):2:(sqrt($2)) w yerrorbars notitle pt 7 ps 0.5,f(x) lc rgb "#007700" lw 2 lt 2 notitle



set size 1,0.29

set ytics 200
set xtics 5
set xtics nomirror
set ytics nomirror
set origin 0.0,0.36
set format x ""
x0=0;g=10;r=0.5;t1=5.6;t2=13.2;p=0.01;bg=0.01
set label 1 font "Helvetica,12" "signal: {/Symbol \367} V {/Symbol \361}" at first 25, first 400
set label 2 font "Helvetica,12" "idler: {/Symbol \367} H {/Symbol \361}" at first 25, first 250
fit [-0:] f(x) 'pump_HV_out_VH_processed.dat' u ($1-459.2):2 via g,p,r
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 "#007700" lw 2 lt 2 notitle

set size 1,0.37
set origin 0.0,0.01
unset xlabel
unset label
set ytics 200
set xtics 5
set xtics nomirror
set ytics nomirror
#set format x ""
set format y "%4.0f"
set format x "%2.0f"
#x0=0;g=10;r=1.5;t1=5.6;t2=13.1;p=180;bg=15
x0=0;g=10;r=1.5;t1=5.6;t2=13.1;p=180;bg=10
set xlabel "{/Symbol D}t (ns)" font "Helvetica,12"
#set ylabel "Coincidence rate" font "Helvetica,12"
set label 3 font "Helvetica,12" "Coincidence events" at first -15, first 600 rotate by 90
#set label font "Helvetica,18" "Count rate" at first -15, first 6 rotate by 90
set label 1 font "Helvetica,12" "signal: {/Symbol \367} L {/Symbol \361}" at first 25, first 400
set label 2 font "Helvetica,12" "idler: {/Symbol \367} L {/Symbol \361}" at first 25, first 250
fit [0:] f(x) 'pump_HV_out_LL_processed.dat' u ($1-459.2):2 via g,p,r
plot [-10:40] 'pump_HV_out_LL_processed.dat' u ($1-459.2):2:(sqrt($2)) w yerrorbars notitle pt 7 ps 0.5 lc rgb "red", f(x) lc rgb "#007700" lw 2 lt 2 notitle


unset multiplot