set term epslatex color 12 size 14.9cm,8cm 
set encoding utf8
set output "g2X_laser.tex"

# raw data
raw_data_1="g2X_laser.dat"
fitted_data_1="g2X_laser_fitted.dat"


symbols_color_1 = "#0000FF"
# symbols_color_2 = "#cc79a7"
# symbols_color_3 = "#009e73"


time_offset = 157.27 #bin number
# time_offset_2 = 242 #bin number
bin_size = 2 #ns


lw_number = 3

laser_lower_range = 0.45
laser_upper_range = 1.05
laser_normalisation_factor = 587155661792/(22544776*21350362*16)

laser_lower_y2_range = laser_lower_range/laser_normalisation_factor
laser_upper_y2_range = laser_upper_range/laser_normalisation_factor

# subthreshold
# set size 1,0.33
# set origin 0,0.66
#set title "30.0 mA"
set xrange [-500:500]
set yrange [laser_lower_range:laser_upper_range]
set xlabel "time difference between photoevents $\\tau$ (ns)"
#set style histogram clustered
set bars large
set ylabel "$g^{(2X)}$"
set ytics nomirror 0.5,0.1,1.0 format "%.2f"
set y2tics nomirror 7000,2000,13000 format "%.0f"
set y2label "coincidences" 
set y2range [laser_lower_y2_range:laser_upper_y2_range]
# set label 1 "30.0mA" at -285,0.55
set key right bottom 

plot raw_data_1 using ($1-time_offset)*bin_size:($3):(sqrt($2)/$2*$3) with yerrorbars linestyle 1 lw lw_number ps 0 linetype rgb symbols_color_1 notitle,fitted_data_1 using ($1):($3) with lines lw lw_number linetype rgb "#000000" notitle
# # # near threshold
# set size 1,0.33
# set origin 0,0.365
# #set title "38.4 mA"
# set ylabel "g^{(2X)}"
# set label 1 "38.4mA" at -285,0.55

# plot fitted_data_2 using ($1-(time_offset))*bin_size:3 with lines linetype rgb "#000000" notitle, raw_data_2 using ($1-time_offset)*bin_size:3:(sqrt($2)/$2*$3) with yerrorbars linestyle 1 pointsize 0 linetype rgb symbols_color notitle


# # above threshold
# set size 1,0.40
# set origin 0,0
# #set title "60.0 mA"
# set xlabel "t_{2} - t_{1} (ns)"
# set ylabel "g^{(2X)}"
# set label 1 "60.0mA" at -285,0.55
# unset format x

# plot fitted_data_3 using ($1-(time_offset))*bin_size:3 with lines linetype rgb "#000000" notitle, raw_data_3 using ($1-time_offset)*bin_size:3:(sqrt($2)/$2*$3) with yerrorbars linestyle 1 ps 0 linetype rgb symbols_color notitle





# unset multiplot