function [freq,sp] = import_1124data(KOACH,ErrFB,varargin)
Freq = containers.Map;
Abs = containers.Map;
Phase = containers.Map;
if strcmp(KOACH,'on')
    disp('KOACH on')
    if strcmp(ErrFB,'ERR_IMC')
        for ii = 1:5
            mnimportdata(['E1124_' num2str(2*ii) '.txt'],Freq,Abs,Phase...
                ,['ERR_IMC' num2str(ii)],'sp');
        end
        mncombine(Freq,Abs,'ERR_IMC',5);
        freq = Freq('ERR_IMC');
        % divide by the SR560 gain
        SR560 = 100;
        sp = Abs('ERR_IMC')/SR560;
    elseif strcmp(ErrFB,'FB_AOM')
        for ii = 1:5
        mnimportdata(['E1124_' num2str(2*ii-1) '.txt'],Freq,Abs,Phase...
            ,['FB_IMC' num2str(ii)],'sp');
        end
        mncombine(Freq,Abs,'FB_IMC',5);
        freq = Freq('FB_IMC');
        % divide by the SR560 gain
        sp = Abs('FB_IMC');
    elseif strcmp(ErrFB,'TTFSS_OUT2')
        for ii = 1:5
        mnimportdata(['E1124_' num2str(ii+32) '.txt'],Freq,Abs,Phase...
            ,['ERR_FSS' num2str(ii)],'sp');
        end
        mncombine(Freq,Abs,'ERR_FSS',5);
        freq = Freq('ERR_FSS');
        % divide by the SR560 gain
        SR560 = 100;
%         if isempty(varargin)
%             error('Put TTFSS_Elec container as input argument')
%         end
%         TTFSS_elec = varargin{1};
        sp = Abs('ERR_FSS')/SR560;%/TTFSS_elec('common_gain_1');
    elseif strcmp(ErrFB,'FB_PZT')
        for ii = 1:5
            mnimportdata(['E1124_' num2str(2*ii+42) '.txt'],Freq,Abs,Phase...
                ,['FB_PZT' num2str(ii)],'sp');
        end
        mncombine(Freq,Abs,'FB_PZT',5);
        freq = Freq('FB_PZT');
        sp = Abs('FB_PZT');
    elseif strcmp(ErrFB,'RTM PZT MON')
        mnimportdata('F1124_58_PZTFB.txt',Freq,Abs,Phase,'RTM PZT MON','sp','diaggui')
        [temp_ff1, temp_abs1] = mntrim(0,2,Freq('RTM PZT MON'),Abs('RTM PZT MON'));
        [temp_ff2, temp_abs2] = mntrim(2.1,1e9,Freq('RTM PZT MON'),Abs('RTM PZT MON')); 
        temp_ff2 = downsample(temp_ff2,50);
        temp_abs2 = downsample(temp_abs2,50);
        freq = [temp_ff1; temp_ff2];
        sp = [temp_abs1; temp_abs2];
        
    elseif strcmp(ErrFB,'RTM MIXER MON')
        mnimportdata('E1124_58_FSSERR.txt',Freq,Abs,Phase,'RTM MIXER MON','sp','diaggui')
        [temp_ff1, temp_abs1] = mntrim(0,2,Freq('RTM MIXER MON'),Abs('RTM MIXER MON'));
        [temp_ff2, temp_abs2] = mntrim(2.1,1e9,Freq('RTM MIXER MON'),Abs('RTM MIXER MON')); 
        temp_ff2 = downsample(temp_ff2,50);
        temp_abs2 = downsample(temp_abs2,50);
        freq = [temp_ff1; temp_ff2];
        sp = [temp_abs1; temp_abs2];
        
    elseif strcmp(ErrFB,'ACC')
        mnimportdata('N1124_58_ACC.txt',Freq,Abs,Phase,'ACC','sp','diaggui')
        [temp_ff1, temp_abs1] = mntrim(0,2,Freq('ACC'),Abs('ACC'));
        [temp_ff2, temp_abs2] = mntrim(2.1,1e9,Freq('ACC'),Abs('ACC')); 
        temp_ff2 = downsample(temp_ff2,50);
        temp_abs2 = downsample(temp_abs2,50);
        freq = [temp_ff1; temp_ff2];
        sp = [temp_abs1; temp_abs2];
    else
        error('1124 IMC sus FB signal has unkown noies')
        mnimportdata('F1124_13.txt',Freq,Abs,Phase,'FB_IMC','sp','diaggui')
        [temp_ff1, temp_abs1] = mntrim(0,2,Freq('FB_IMC'),Abs('FB_IMC'));
        [temp_ff2, temp_abs2] = mntrim(2.1,1e9,Freq('FB_IMC'),Abs('FB_IMC')); 
        temp_ff2 = downsample(temp_ff2,50);
        temp_abs2 = downsample(temp_abs2,50);
        freq = [temp_ff1; temp_ff2];
        sp = [temp_abs1; temp_abs2];
    end
else
    disp('KOACH off')
    if strcmp(ErrFB,'ERR_IMC')
        for ii = 1:5
            mnimportdata(['E1124_' num2str(2*ii+13) '.txt'],Freq,Abs,Phase...
                ,['ERR_IMC' num2str(ii)],'sp');
        end
        mncombine(Freq,Abs,'ERR_IMC',5);
        freq = Freq('ERR_IMC');
        % divide by the SR560 gain
        SR560 = 100;
        sp = Abs('ERR_IMC')/SR560;
    elseif strcmp(ErrFB,'FB_AOM')
        for ii = 1:5
        mnimportdata(['E1124_' num2str(2*ii+12) '.txt'],Freq,Abs,Phase...
            ,['FB_IMC' num2str(ii)],'sp');
        end
        mncombine(Freq,Abs,'FB_IMC',5);
        freq = Freq('FB_IMC');
        % divide by the SR560 gain
        sp = Abs('FB_IMC');
    elseif strcmp(ErrFB,'TTFSS_OUT2')
        SR560 = 1000;
        for ii = 1:5
            mnimportdata(['E1124_' num2str(ii+57) '.txt'],Freq,Abs,Phase...
            ,['ERR_FSS' num2str(ii)],'sp');
            if ii == 5
                Abs('ERR_FSS5') = Abs('ERR_FSS5')*10;
            end
        end
        mncombine(Freq,Abs,'ERR_FSS',5);
        freq = Freq('ERR_FSS');
        sp = Abs('ERR_FSS')/SR560;
    elseif strcmp(ErrFB,'FB_PZT')
        for ii = 1:5
            mnimportdata(['E1124_' num2str(2*ii+67) '.txt'],Freq,Abs,Phase...
                ,['FB_PZT' num2str(ii)],'sp');
        end
        mncombine(Freq,Abs,'FB_PZT',5);
        freq = Freq('FB_PZT');
        sp = Abs('FB_PZT');
        
    elseif strcmp(ErrFB,'RTM PZT MON')
        mnimportdata('F1124_59_PZTFB.txt',Freq,Abs,Phase,'RTM PZT MON','sp','diaggui')
        [temp_ff1, temp_abs1] = mntrim(0,2,Freq('RTM PZT MON'),Abs('RTM PZT MON'));
        [temp_ff2, temp_abs2] = mntrim(2.1,1e9,Freq('RTM PZT MON'),Abs('RTM PZT MON')); 
        temp_ff2 = downsample(temp_ff2,50);
        temp_abs2 = downsample(temp_abs2,50);
        freq = [temp_ff1; temp_ff2];
        sp = [temp_abs1; temp_abs2];
        
    elseif strcmp(ErrFB,'RTM MIXER MON')
        mnimportdata('E1124_59_FSSERR.txt',Freq,Abs,Phase,'RTM MIXER MON','sp','diaggui')
        [temp_ff1, temp_abs1] = mntrim(0,2,Freq('RTM MIXER MON'),Abs('RTM MIXER MON'));
        [temp_ff2, temp_abs2] = mntrim(2.1,1e9,Freq('RTM MIXER MON'),Abs('RTM MIXER MON')); 
        temp_ff2 = downsample(temp_ff2,50);
        temp_abs2 = downsample(temp_abs2,50);
        freq = [temp_ff1; temp_ff2];
        sp = [temp_abs1; temp_abs2];
        
    
    elseif strcmp(ErrFB,'ACC')
        mnimportdata('N1124_59_ACC.txt',Freq,Abs,Phase,'ACC','sp','diaggui')
        [temp_ff1, temp_abs1] = mntrim(0,2,Freq('ACC'),Abs('ACC'));
        [temp_ff2, temp_abs2] = mntrim(2.1,1e9,Freq('ACC'),Abs('ACC')); 
        temp_ff2 = downsample(temp_ff2,50);
        temp_abs2 = downsample(temp_abs2,50);
        freq = [temp_ff1; temp_ff2];
        sp = [temp_abs1; temp_abs2];
    else
        mnimportdata('F1124_25.txt',Freq,Abs,Phase,'FB_IMC','sp','diaggui')
        [temp_ff1, temp_abs1] = mntrim(0,2,Freq('FB_IMC'),Abs('FB_IMC'));
        [temp_ff2, temp_abs2] = mntrim(2.1,1e9,Freq('FB_IMC'),Abs('FB_IMC')); 
        temp_ff2 = downsample(temp_ff2,50);
        temp_abs2 = downsample(temp_abs2,50);
        freq = [temp_ff1; temp_ff2];
        sp = [temp_abs1; temp_abs2];
    end
end
end