#!/usr/bin/gnuplot

infile1 = 'data/rabi_2019-12-16.dat'
set terminal postscript eps color enhanced solid font "Helvetica,30" size 12cm,10cm
set output "rabi.eps"
set samples 1000
set colorsequence classic

a0=-0.405
#t0 = 12
y0 = 0.926
tau = 700
f0 = 70
rabi(x) = a0 - exp(-x/tau)*a0*cos(2*3.14*x*f0/1000) + y0

fit rabi(x) infile1 u 1:5:6 via a0,f0,y0,tau

#set label 1 at 5, 5
#set label 1 "Rabi oscillation" tc lt 1

#set title "MW-Rabi, lin pol trap"
set xlabel "Pulse length ( {/Symbol m}s)"
set ylabel 'Probability in F=2 (%)'
set key right bottom
set pointsize 1.0
set xrange [0:65]
set yrange [0:100]
set ytics 25
set xtics 15
#set label 1 sprintf("Lorentz+Dispersion fit:") at graph 0.02,0.45

plot 100*rabi(x) lc 1 lw 1.5 not, infile1 u 1:(100*$5):(100*$6) w yerr ps 1.2 pt 7 lc 1 lw 2 not

unset output
#pause -1