#!/usr/local/bin/gnuplot -persist
set term postscript enhanced eps color font "Helvetica,16" size 8.6cm,5.5cm 
set encoding utf8
set output "g2-loops.eps"

point_size = 0.7
set pointinterval 0.0001
set multiplot

f(x) = 1+2*x*(1-x)
f_corr(x) = 1+ 2*x*(1-x)*0.886
# fabove(x) = x <   0.5 ? f(x) : f(0.5)
# fbelow(x) = x >=  0.5 ? f(x) : 0 

# gabove(x) = x < 0.5 ? 2*x :  f (0.5)
# gbelow(x) = x >= 1 ? 2*x :  fbelow(x) 
# g(x) = x <= 0.5 ? 2*x : 0
# ghigh(x) = x >= 1 ? 2*x : 0

# Avoid the cross-system mess
# raw_data_0=".\\g2_merged_data\\normal_alpha-0p000.dat"
raw_data_0="normal_alpha-0p000.dat"
#raw_data_1=".\\g2_merged_data\\normal_EYP-DFB-0780-100mA_25degC_2_emitter_400m_delay_ts_130ns_20240829_2234.g2_merged.dat"
raw_data_1="normal_EYP-DFB-0780-100mA_25degC_2_emitter_400m_delay_ts_130ns_20240829_2234.g2_merged.dat"
#raw_data_2=".\\g2_merged_data\\normal_EYP-DFB-0780-100mA_25degC_4_emitter_400m_186m_delay_ts_130ns_2out2in_20240830_0025.g2_merged.dat"
raw_data_2="normal_EYP-DFB-0780-100mA_25degC_4_emitter_400m_186m_delay_ts_130ns_2out2in_20240830_0025.g2_merged.dat"
#raw_data_3=".\\g2_merged_data\\normal_EYP-DFB-0780-100mA_25degC_8_emitter_400m_186m_100m_delay_ts_130ns_2out2in_20240830_0052.g2_merged.dat"
raw_data_3="normal_EYP-DFB-0780-100mA_25degC_8_emitter_400m_186m_100m_delay_ts_130ns_2out2in_20240830_0052.g2_merged.dat"

#fitted_data_1=".\\normal_g2_fit\\normal_EYP-DFB-0780-100mA_25degC_2_emitter_400m_delay_ts_130ns_20240829_2234.g2_merged_fitted.dat"
#fitted_data_2=".\\normal_g2_fit\\normal_EYP-DFB-0780-100mA_25degC_4_emitter_400m_186m_delay_ts_130ns_2out2in_20240830_0025.g2_merged_fitted.dat"
#fitted_data_3=".\\normal_g2_fit\\normal_EYP-DFB-0780-100mA_25degC_8_emitter_400m_186m_100m_delay_ts_130ns_2out2in_20240830_0052.g2_merged_fitted.dat"
fitted_data_1="normal_EYP-DFB-0780-100mA_25degC_2_emitter_400m_delay_ts_130ns_20240829_2234.g2_merged_fitted.dat"
fitted_data_2="normal_EYP-DFB-0780-100mA_25degC_4_emitter_400m_186m_delay_ts_130ns_2out2in_20240830_0025.g2_merged_fitted.dat"
fitted_data_3="normal_EYP-DFB-0780-100mA_25degC_8_emitter_400m_186m_100m_delay_ts_130ns_2out2in_20240830_0052.g2_merged_fitted.dat"


symbols_color_1 = "#000000"
symbols_color_2 = "#0072b2"
symbols_color_3 = "#D55e00"
symbols_color_4 = "#CC79A7"

lw_number  = 2
bin_size = 2
time_offset  = 67.5


# set lmargin screen 0.1
# set rmargin screen 0.9
# set size 1,0.5
# set origin 0,0.5
set ytics format "%.1f"
set xrange [-300:300]
set xtics 100
# set mxtics 4
set yrange [0.9:1.9]
set ytics 0.2
# set mytics 4
set errorbars small
set ylabel "g^{(2)}" 
set xlabel "two-photoevent separation {/Symbol t} (ns)"
# set style fill pattern 2
# set xtics (0.5, 0.75 ,1)
# set xlabel "Two-Photoevent Separation {/Symbol t} (ns)"
# #set ylabel "Two-Photon Coincidences"
# set ylabel "g^{(2)}({/Symbol t})" #rotate by -90

# set arrow from 0.5, graph 0 to 0.5, graph 1 nohead lc rgb 'black' dt 2 front
# set arrow from 2/3, graph 0 to 2/3, graph 1 nohead lc rgb 'black' dt 2 front

unset key
# set key top left samplen 0.5
set label "n=0" at 0, 1.08 center
set label "n=1" at 0, 1.27 center
set label "n=2" at -60,1.70 center
set arrow from -35,1.70 to 0, 1.66 nohead lc rgb 'black' dt 1 front
set label "n=3" at 0, 1.835 center



plot	raw_data_0 using ($1-time_offset)*bin_size:($3):(sqrt($2)/$2*$3) with yerrorbars linestyle 1 pt 12 lw lw_number pointsize point_size linetype rgb symbols_color_1 title "n=0",\
		raw_data_1 using ($1-time_offset)*bin_size:($3):(sqrt($2)/$2*$3) with yerrorbars linestyle 2 pt 8 lw lw_number pointsize point_size linetype rgb symbols_color_2 title "n=1",\
 		raw_data_2 using ($1-time_offset)*bin_size:($3):(sqrt($2)/$2*$3) with yerrorbars linestyle 3 pt 6 lw lw_number pointsize point_size linetype rgb symbols_color_3 title "n=2", \
		raw_data_3 using ($1-time_offset)*bin_size:($3):(sqrt($2)/$2*$3) with yerrorbars linestyle 4 pt 4 lw lw_number pointsize point_size linetype rgb symbols_color_4 title "n=3", \
		fitted_data_1 using ($1)*bin_size/2:($3) with lines lw lw_number linetype rgb "#000000" notitle,\
		fitted_data_2 using ($1)*bin_size/2:($3) with lines lw lw_number linetype rgb "#000000" notitle,\
	 	fitted_data_3 using ($1)*bin_size/2:($3) with lines lw lw_number linetype rgb "#000000" notitle, \


# # unset object 1
# # set size 0.55,0.55
# # set origin 0.4,0.1
# unset ytics 

# unset label 1
# set size 1,0.50
# set origin 0,0
# # unset arrow 5
# # unset arrow 6

# set ytics format "%.2f"
# set ytics 0, 0.1
# set xrange[0:0.51]
# set yrange[0.0:0.55]
# set label 1 "(b)" at graph 0.02, graph 0.90
# # set xzeroaxis ls -1

# set ylabel "$g^{(2)}(0)-1$"
# set xlabel "power ratio $\\alpha$"
# set key top left

# plot 	"g2_vs_alpha.txt" u 1:($2):3 with yerrorbars linestyle 1 lw 3 ps 0 linetype rgb "#FF0000" notitle,\
# 		f(x)-1 w l lw 2 lc rgb "black" notitle,\
# 		f_corr(x)-1 w l lw 2 dashtype 3 lc rgb "blue" notitle

# # plot 	g(x) w l lw 2 lc rgb "black" notitle,\
# #      	gbelow(x) w l lw 2 lc rgb "black" notitle
# #      	# f(x) w l lw 2 lc rgb "black" notitle,\
# #      	# (18*x -12)  w l lw 2 lc rgb "black" notitle

# set ytics border
# set xtics border