#!/usr/bin/gnuplot
#Plots two time-resolved things: the transmission for three incoming photons as well as the ratio of the transmission.

#set terminal qt size 350,262 enhanced font 'Arial,15' persist
# Line width of the axes
set border linewidth 1.5

## input file ##
infile_long_dec='data_fig_3/non_reversed_with_alternating_09_26_onwards_combined_histo_tx_2ns'
infile_5ns_dec='data_fig_3/2015_11_03-11_06_non_reversed_combined_histo_tx_1ns'
infile_long_ris='data_fig_3/reversed_with_alternating_09_26_onwards_combined_histo_tx_2ns'

infile_long_dec_theory ='data_fig_3/trans_ratio_long_dec_w_bg'
infile_5ns_dec_theory ='data_fig_3/trans_ratio_5ns_dec_w_bg'
infile_long_ris_theory ='data_fig_3/trans_ratio_long_ris_w_bg'
infile_long_dec_theory_wo_bg ='data_fig_3/trans_ratio_long_dec_wo_bg'
infile_5ns_dec_theory_wo_bg ='data_fig_3/trans_ratio_5ns_dec_wo_bg'
infile_long_ris_theory_wo_bg ='data_fig_3/trans_ratio_long_ris_wo_bg'
#################


set terminal postscript eps color enhanced dashed "Helvetica,22" dl 2 size 25cm, 11cm
set output "fig3.eps"

#set border 3
set nobar
set multiplot 


###################
### TOP PLOTS  ####
###################

set origin 0,0.5
set size 0.37,0.45
set key at graph 1.01,0.95 right spacing 1.3 samplen 0.5 width 0.5 

#set tmargin 2
set lmargin 9
#set rmargin 0
#set bmargin 0.5
set mytics 10
set xtics nomirror
set ytics scale 2.5,1 nomirror
set label 1 "{/Symbol t}_p = 5.6ns" at graph 0.5,1.1 center
#set label 2 "1ns time-bin" at graph 0.47,0.8

#set xr [-1:5]
#set xtics (-1,0,1,2,3,4,5)
set xr [-1*5.61:5.05*5.61]

set xtics format " " 
set ytics format "%g" 
set yr [0.4:100]
#set ytics (1,10,100)
unset xlabel
set ylabel "Coincidence\n probability (x 10^{-5})"
set logscale y 10
unset grid
unset bars
plot infile_5ns_dec u (($2-909.5)):($9*1e5):($10*1e5)  with yerrorbars pt 7 ps 0.9 lc rgb  "black" notitle,\
'' u (($2-909.5)):($9*1e5) w p pt 7 ps 0.9 lc rgb "black" t "G_{f}, with atom",\
infile_5ns_dec u (($2-909.5)):($5*1e5):($6*1e5) with yerrorbars pt 13 ps 0.9 lc rgb "red" notitle,\
'' u (($2-909.5)):($5*1e5) w p pt 13 ps 0.9 lc rgb "red" t "G_{f,0}, w/o atom"

####################################################################
#set lmargin 3
#set rmargin 3

set origin 0.31,0.5
set size 0.37,0.45

set label 1 "{/Symbol t}_p = 13.2ns" at graph 0.5,1.1 center
#set label 2 "2ns time-bin" at graph 0.485,0.8
unset key

####
set xr [-1*13.19:5.05*13.19]
set xtics format " " 
set ytics format " " 
set yr [0.4:100]
unset xlabel
unset ylabel
set logscale y 10
unset grid
unset bars
plot infile_long_dec u (($2-881.5)):($9*1e5):($10*1e5)  with yerrorbars pt 7 ps 0.9 lc rgb "black", infile_long_dec u (($2-881.5)):($5*1e5):($6*1e5) with yerrorbars pt 13 ps 0.9 lc rgb "red"

####################################################################
# set lmargin 3
# set rmargin 3
set origin 0.62,0.5
set size 0.37,0.45

set label 1 "{/Symbol t}_p = 13.2ns" at graph 0.5,1.1 center
#set label 2 "2ns time-bin" at graph 0.255,0.8

######
set xr [-5.05*13.19:1*13.19]
set yr [0.4:100]
#set xtics (-5,-4,-3,-2,-1,0,1)
set ytics format " " 
set xtics format " " 
unset xlabel
unset ylabel
set logscale y 10
unset grid
unset bars
plot infile_long_ris u (($2-878.5)):($9*1e5):($10*1e5)  with yerrorbars pt 7 ps 0.9 lc rgb "black", infile_long_ris u (($2-878.5)):($5*1e5):($6*1e5) with yerrorbars pt 13 ps 0.9 lc rgb "red"

####################################################################

######################
### BOTTOM PLOTS  ####
######################

# set lmargin 8
# set bmargin 4
set tmargin 0
#set rmargin 0
set origin 0,0
set size 0.37,0.54

unset label 1
unset label 2
###########
unset log
set xr [-1*5.61:5.05*5.61]
set xtics format "%g" 
set ytics format "%g" 
#set xtics (-1,0,1,2,3,4,5)
set yr [0.8:1.1]
set mytics 2
set ytics 0.1 scale 2,1
#set xlabel "Detection time delay {/Symbol D}t_{hf} ({/Symbol t}_{p})"
set xlabel "Time (ns)"
#set ylabel "ratio: atom/no atom"
set ylabel "G_{f} / G_{f,0}"
unset grid
unset bars
plot infile_5ns_dec_theory_wo_bg u ($1):2 w lines lt 1 lc rgb "web-blue" lw 2, infile_5ns_dec u (($2-909.5)):($5/$9):(sqrt( ($6/$9)**2 + ($5**2*$10**2/$9**4))) with yerrorbars pt 6 ps 0.8 lt 1 lc rgb "red", infile_5ns_dec_theory u ($1):2 w lines lt 4 lc rgb "dark-green" lw 3
####################################################################
#set lmargin 5
#set rmargin 3
set origin 0.31,0
set size 0.37,0.54
unset log
set xr [-1*13.19:5.05*13.19]
set xtics format "%g" 
set ytics format " " 
#set xtics (-1,0,1,2,3,4,5)
set yr [0.8:1.1]
set mytics 2
set ytics 0.1 scale 2,1
#set xlabel "Detection time delay {/Symbol D}t_{hf} ({/Symbol t}_{p})"
set xlabel "Time (ns)"
unset ylabel
unset grid
unset bars
plot infile_long_dec_theory_wo_bg u ($1):2 w lines lt 1 lc rgb "web-blue" lw 2, infile_long_dec u (($2-881.5)):($5/$9):(sqrt( ($6/$9)**2 + ($5**2*$10**2/$9**4))) with yerrorbars pt 6 ps 0.8 lt 1  lc rgb "red", infile_long_dec_theory u ($1):2 w lines lt 4 lc rgb "dark-green" lw 3
####################################################################

####################################################################
#set lmargin 5
#set rmargin 3
set origin 0.62,0
set size 0.37,0.54
unset log
set xr [-5.05*13.19:1*13.19]
set yr [0.8:1.1]
set xtics format "%g" 
set ytics format " " 
#set xtics (-5,-4,-3,-2,-1,0,1)
set mytics 2
set ytics 0.1 scale 2,1
set ytics format " " 
#set xlabel "Detection time delay {/Symbol D}t_{hf} ({/Symbol t}_{p})"
set xlabel "Time (ns)"
unset ylabel
unset grid
unset bars
plot infile_long_ris_theory_wo_bg u ($1):2 w lines lt 1 lc rgb "web-blue" lw 2, infile_long_ris u (($2-878.5)):($5/$9):(sqrt( ($6/$9)**2 + ($5**2*$10**2/$9**4))) with yerrorbars pt 6 ps 0.8 lt 1 lc rgb "red", infile_long_ris_theory u ($1):2 w lines lt 4 lc rgb "dark-green" lw 3

unset multiplot                     # exit multiplot mode (prompt changes back to 'gnuplot')

unset output

pause -1