#!/usr/bin/gnuplot

set terminal postscript eps color enhanced solid font "Helvetica,22" size 12.6cm,10cm
set termoption dashed
set output "../scat_tx_res.eps"
#set output 'test.eps'

f1 = '2016-08-08/scat_tx'
f2 = '2016-08-08/scat_resfreq'
f3 = 'model.dat'

set border linewidth 1.5
set border 3
set samples 5000
set multiplot 
set bar 2
unset key

#################################
#######
set origin 0,0
set size 1,0.55
unset key
set border 3
set lmargin 10
set xtics nomirror
set ytics nomirror

set label 1 "b)" at graph -0.2,1 center
scale = 0.458

set ylabel 'extinction (%)' #offset -1.7
set key right bottom
#set yrange [11.5:19.5]
set xrange [0:550]
unset key

set ytics 10
set xlabel 'number of scattered photons'

plot f1 u 1:2:3 w ye not lw 3 lt 7 lc rgb 'red', f3 u 1:4 w l lw 3 lc rgb 'blue' dt 2, f3 u 1:3 w l lw 3 lc rgb 'red'

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

set origin 0,0.55
set size 1,0.45
set border 3
#set tmargin 2
#set bmargin 4
#set rmargin 0

unset label 1
set label 2 "a)" at graph -0.2,1 center

set ylabel 'resonance shift (MHz)'
#set xtics 25
set ytics 1
set xrange [0:550]
set yrange [45:47]
set xlabel ''
set xtics format ''

unset grid

plot f3 u 1:2 w l lw 3 lc rgb 'red', f2 u 1:2:3 w ye not lw 3 lt 7 lc rgb 'red'
#################################
#######


unset multiplot   
unset output
pause -1