#!/usr/bin/gnuplot

f1 = 'deep_trap/tmp'
f2 = 'shallow_trap/tmp'

set term pngcairo #size 800, 600
set output 'tmp.png'
#set term pdf
#set output 'tmp.pdf'

#fitting to gaussian
g(x)=a*exp(-x/b)+c
a=100
b= 10
c=60
#fit g(x) file u ($1/1000):($5*100):($6*100) via a,b,c


#set label 1 sprintf("repump tau = %g ms\n", b) at graph 0.6,0.35

#set yrange [0:100]
set title 'depolarisation'
set xlabel 'wait time (ms)'
set ylabel 'bright state population (%)'
set grid
set key

plot f2 u ($1):($5*100):($6*100) w yerr lc 1 pt 7 ps 0.5 t 'shallowtrap',\
f1 u ($1):($5*100):($6*100) w yerr pt 7 ps 0.5 t 'deeptrap'

unset output