#set multiplot layout 2,1

filename="processed_histogram_with_yerror_bars/results.sorted"
rootdir="/mnt/buffalo/"
square="sq_pulse/"
exp="exp_pulse/vary_photon_number_perkin_elmer_"

#### Columns for square pulse ####
data_filename=1
p_tran=2
p_tran_err=3
od_used=4
photon_number=5
photon_err=6
p_refl_flou=7
p_refl_flou_err=8
p_refl_flou_bg=9
p_refl_flou_bg_err=10
p_refl_scat=11
p_refl_scat_err=12
flou_lifetime=13
flou_lifetime_err=14
max_p_refl=15
max_p_refl_err=16

#### Columns for exponential pulse ####
exp_risetime=17
exp_risetime_err=18
fit_lambda=19
fit_lambda_err=20
fit_photon_number=21
fit_photon_err=22

#### End Columns ####

currentdir="`pwd`"
cd rootdir
unset logscale xy
set logscale x
set key left
# plot square.'5ns/'.filename u photon_number:p_refl_flou:photon_err:p_refl_flou_err w xyerrorl,\
#      square.'5ns/'.filename u photon_number:max_p_refl:photon_err:max_p_refl_err w xyerror axes x1y2,\
#      square.'15ns/'.filename u photon_number:p_refl_flou:photon_err:p_refl_flou_err w xyerrorl,\
#      square.'15ns/'.filename u photon_number:max_p_refl:photon_err:max_p_refl_err w xyerror axes x1y2,\
#      square.'20ns/'.filename u photon_number:p_refl_flou:photon_err:p_refl_flou_err w xyerrorl,\
#      square.'20ns/'.filename u photon_number:max_p_refl:photon_err:max_p_refl_err w xyerror axes x1y2,\
#      square.'60ns/'.filename u photon_number:p_refl_flou:photon_err:p_refl_flou_err w xyerrorl,\
#      square.'60ns/'.filename u photon_number:max_p_refl:photon_err:max_p_refl_err w xyerror axes x1y2,\
#      square.'250ns/'.filename u photon_number:p_refl_flou:photon_err:p_refl_flou_err w xyerrorl,\
#      square.'250ns/'.filename u photon_number:max_p_refl:photon_err:max_p_refl_err w xyerror axes x1y2

#set xrange restore ; #[1:1000]

# plot exp.'5ns/'.filename u 7:9:8:10 w xyerrorl,\
#      exp.'5ns/'.filename u 7:17:8:18 axes x1y2 w xyerror,\
#      exp.'15ns/'.filename u 7:9:8:10 w xyerrorl,\
#      exp.'15ns/'.filename u 7:17:8:18 axes x1y2 w xyerror,\
#      exp.'25ns/'.filename u 7:9:8:10 w xyerrorl,\
#      exp.'25ns/'.filename u 7:17:8:18 axes x1y2 w xyerror,\
#      exp.'60ns/'.filename u 7:9:8:10 w xyerrorl,\
#      exp.'60ns/'.filename u 7:17:8:18 axes x1y2 w xyerror,\
#      exp.'250ns/'.filename u 7:9:8:10 w xyerrorl,\
#      exp.'250ns/'.filename u 7:17:8:18 axes x1y2 w xyerror


set title '5ns'
plot square.'5ns/'.filename u photon_number:max_p_refl:photon_err:max_p_refl_err w xyerror axes x1y1 t 'Square pulse',\
     exp.'5ns/'.filename u photon_number:max_p_refl:photon_err:max_p_refl_err axes x1y1 w xyerror t 'Exponential pulse' ls 3

pause -1

set title '15ns'
plot square.'15ns/'.filename u photon_number:max_p_refl:photon_err:max_p_refl_err w xyerror axes x1y1 t 'Square pulse',\
     exp.'15ns/'.filename u photon_number:max_p_refl:photon_err:max_p_refl_err axes x1y1 w xyerror t 'Exponential pulse' ls 3

pause -1

set title '60ns'
plot square.'60ns/'.filename u photon_number:max_p_refl:photon_err:max_p_refl_err w xyerror axes x1y1 t 'Square pulse',\
     exp.'60ns/'.filename u photon_number:max_p_refl:photon_err:max_p_refl_err axes x1y1 w xyerror t 'Exponential pulse' ls 3

pause -1

set title '250ns'
plot square.'250ns/'.filename u photon_number:max_p_refl:photon_err:max_p_refl_err w xyerror axes x1y1 t 'Square pulse',\
     exp.'250ns/'.filename u photon_number:max_p_refl:photon_err:max_p_refl_err axes x1y1 w xyerror t 'Exponential pulse' ls 3

pause -1

unset multiplot

cd currentdir