function [act_aom, sigma_aom] = Fit_Aaom( varargin)
%UNTITLED9 Summary of this function goes here
%   Detailed explanation goes here
% plot or not
if mnismember(varargin,'Plot')
    doplot = 'Plot';
else
    doplot = 'noPlot';
end

% set containers
temp_Sigma = containers.Map;
temp_Act = containers.Map;

% measurement: see klog3172
data = importdata('P0727.txt');
vv = data.data(:,1);
freq = data.data(:,2);

mnlinearFit(vv,freq,temp_Sigma,temp_Act,'AOM',doplot,'title','AOM driver input voltage .vs output RF frequency',...
    'xlabel','input voltage [V]','ylabel','output frequency [MHz]');
act_aom = 2*temp_Act('AOM_p1')*1e6;
sigma_aom = 2*temp_Sigma('AOM_p1')*1e6;

end