ve/del *
 opt ndate
 opt fit
 opt nstat
 
ve/read dnu,gain,c0,mmcor cavity_data.dat ! ! -/*/

sigma c0=c0/mmcor
 
sigma signu=1000.0*dnu/2./sqrt(2.0*log(2)) 
 
 sigma siglas = sqrt(signu**2/c0**2-signu**2)
 sigma dlas = siglas*2.*sqrt(2.0*log(2))
 sigma eslas = 1/siglas*signu**2/c0**3*0.05
 sigma eslas =eslas*2.*sqrt(2.0*log(2))
 
np=$vlen(dnu)
ve/cr ex([np])
ve/cr ey([np]) r [np]*0.05


 
 zone 1 2
 
 null 0 600 0.0 1.1 
set hcol 2
set pmci 2
 hplot/err dnu c0 ex  ey [np] 20 ! s
 
 set hcol 1
 fun/plo convolution.f 1 1E3 c,s
 
 atit 'Cavity linewidth (kHz)' 'c0'
 
 sigma dlas=dlas/1000
 sigma eslas=eslas/1000
 
 null 0 600 0 400
 set hcol 4
 set pmci 4
 hplot/err dnu dlas ex  eslas [np] 20 ! s
 
 ve/fit dnu dlas eslas p0 s,e
 
 atit 'Cavity linewidth (kHz)' 'Verdi linewidth (kHz)'