#!/usr/bin/gnuplot
infile1 = 'tx_refl_varytemp.dat'
set terminal postscript eps color enhanced solid font "Helvetica,22" size 12.6cm,8cm
set output "../tx_refl_varytemp.eps"

set border linewidth 1.5
set border 3
set samples 5000
set multiplot 
set nobar
unset key

##### tx diff plots ####

set origin 0,0
set size 0.47,1
set tmargin 2
set lmargin 8
set bmargin 4
set rmargin 0
set xtics nomirror
set ytics nomirror

#set title 'transmission extinction'
set xlabel 'temperature ({/Symbol m}K)'
set ylabel 'extinction (%)'
set key right bottom
#set yrange [82:101]
set xrange [15:111]
set xtics 25
set ytics 1
unset key

plot infile1 u 11:((1-$8)*100):12:($9*100) w xye lt 7 lc 1 lw 2 not

#################################
#######
set origin 0.50,0
set size 0.47,1
unset key
set border 3
set xtics nomirror

f1='tx_refl_varytemp.dat'

attenuation = 0.0073 #0.012
fibercoupling = 0.62
apdeff = 0.52
opticalpath = 0.9
col_attn = 1.0 #0.726
bg_red = 180*20e-6 # bg per atom bin
bg_blue = 180*20e-6 # bg per atom bin

gamma = 1/26.24e-9
h = 6.626e-34
c = 3e8
lambda = 780.24e-9

eff_correction = 1/(fibercoupling*attenuation*opticalpath) 

set ytics 0.05
#set yrange [0:0.62]

set xlabel 'temperature ({/Symbol m}K)'
set ylabel 'peak reflection (%)' offset 1
unset grid

plot f1 u 11:(($4-bg_blue*$6)/($5*eff_correction)*100):12:(sqrt ( (sqrt($4)/($5*eff_correction) )**2 + (sqrt($5)*$4/($5**2*eff_correction))**2 )*100) w xyerr lc 1 pt 7 lw 2 not

unset multiplot   
unset output
pause -1