#!/usr/bin/gnuplot
infile1 = 'tx_refl_varytemp.dat'
set terminal postscript eps color enhanced solid font "Helvetica,22" size 12.6cm,10cm
set output "../tx_refl_varytemp_vert.eps"
#set output "vert.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.02
set size 1,0.55
#set tmargin 2
set lmargin 10
#set bmargin 4
#set rmargin 0
set xtics nomirror
set ytics 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 xtics 25
set ytics 0.05
set xrange [25:132]
set yrange [0.42:0.63]

set xlabel 'temperature ({/Symbol m}K)'
set ylabel 'peak reflection (%)'
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

#################################
#######
set origin 0,0.5
set size 1,0.45
unset key
set border 3
set xtics nomirror

#set title 'transmission extinction'
set ylabel 'extinction (%)' offset -1.7
set key right bottom
set yrange [13.1:17.8]
unset key

set ytics 1
set xlabel ''
set xtics format ''

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

unset multiplot   
unset output
pause -1