function [Noise_TTFSS,Noise_CMS] = import_noise(TTFSS_Elec,CMS_Elec,Fitting,freq,varargin)

data = importdata('N0728_7.txt');
ff_TTFSS_n1 = data(:,1);
abs_TTFSS_n1 = sqrt(data(:,2));
data = importdata('N0728_8.txt');
ff_TTFSS_n2 = data(:,1);
abs_TTFSS_n2 = sqrt(data(:,2));
data = importdata('N0728_9.txt');
ff_TTFSS_n3 = data(:,1);
abs_TTFSS_n3 = sqrt(data(:,2));
data = importdata('N0728_10.txt');
ff_TTFSS_n4 = data(:,1);
abs_TTFSS_n4 = sqrt(data(:,2));
data = importdata('N0728_11.txt');
ff_TTFSS_n5 = data(:,1);
abs_TTFSS_n5 = sqrt(data(:,2));

[ff_TTFSS_n2,abs_TTFSS_n2] = mntrim(max(ff_TTFSS_n1),max(ff_TTFSS_n2),ff_TTFSS_n2,abs_TTFSS_n2);
[ff_TTFSS_n3,abs_TTFSS_n3] = mntrim(max(ff_TTFSS_n2),max(ff_TTFSS_n3),ff_TTFSS_n3,abs_TTFSS_n3);
[ff_TTFSS_n4,abs_TTFSS_n4] = mntrim(max(ff_TTFSS_n3),max(ff_TTFSS_n4),ff_TTFSS_n4,abs_TTFSS_n4);
[ff_TTFSS_n5,~] = mntrim(max(ff_TTFSS_n4),max(ff_TTFSS_n5),ff_TTFSS_n5,abs_TTFSS_n5);
Freq_TTFSS = [ff_TTFSS_n1; ff_TTFSS_n2; ff_TTFSS_n3; ff_TTFSS_n4];
Output_Noise_TTFSS = [abs_TTFSS_n1; abs_TTFSS_n2; abs_TTFSS_n3; abs_TTFSS_n4];
% loglog(Freq_TTFSS,Output_Noise_TTFSS)

SR560 = 100;
TTFSS_PZT = -(TTFSS_Elec('common_gain1') * TTFSS_Elec('common_gain2') * SR560 ...
    * TTFSS_Elec('fast_gain1') * TTFSS_Elec('fast_LPF_10k') * TTFSS_Elec('fast_LPF_34k')...
    * TTFSS_Elec('fast_notch') * TTFSS_Elec('fast_LPF_out'))*Fitting('TTFSS_PZT')*10^(30/20)*10^(30/20);
servo_c = mnbode_list(TTFSS_PZT,Freq_TTFSS);
Input_Noise_TTFSS = Output_Noise_TTFSS ./ servo_c;
Noise_TTFSS = interp1(Freq_TTFSS,Input_Noise_TTFSS,freq);
for ii = 1:length(Noise_TTFSS)
    if isnan(Noise_TTFSS(ii))
        Noise_TTFSS(ii) = last;
    else
        last = Noise_TTFSS(ii);
    end
end
% loglog(freq,Noise_TTFSS)
%%
data = importdata('N0728_1.txt');
ff_CMS_n1 = data(:,1);
abs_CMS_n1 = sqrt(data(:,2));
data = importdata('N0728_6.txt');
ff_CMS_n2 = data(:,1);
abs_CMS_n2 = sqrt(data(:,2));
data = importdata('N0728_3.txt');
ff_CMS_n3 = data(:,1);
abs_CMS_n3 = sqrt(data(:,2));
data = importdata('N0728_4.txt');
ff_CMS_n4 = data(:,1);
abs_CMS_n4 = sqrt(data(:,2));
data = importdata('N0728_5.txt');
ff_CMS_n5 = data(:,1);
abs_CMS_n5 = sqrt(data(:,2));
[ff_CMS_n2,abs_CMS_n2] = mntrim(max(ff_CMS_n1),max(ff_CMS_n2),ff_CMS_n2,abs_CMS_n2);
[ff_CMS_n3,abs_CMS_n3] = mntrim(max(ff_CMS_n2),max(ff_CMS_n3),ff_CMS_n3,abs_CMS_n3);
[ff_CMS_n4,abs_CMS_n4] = mntrim(max(ff_CMS_n3),max(ff_CMS_n4),ff_CMS_n4,abs_CMS_n4);
[ff_CMS_n5,abs_CMS_n5] = mntrim(max(ff_CMS_n4),max(ff_CMS_n5),ff_CMS_n5,abs_CMS_n5);
Freq_CMS = [ff_CMS_n1; ff_CMS_n2; ff_CMS_n3; ff_CMS_n4;ff_CMS_n5];
Output_Noise_CMS = [abs_CMS_n1; abs_CMS_n2; abs_CMS_n3; abs_CMS_n4; abs_CMS_n5];



SR560 = 10;
CMS = CMS_Elec('common_comp') *  CMS_Elec('slow_LPF') * 10^(31/20) * Fitting('CMS_AOM') * SR560;
servo_c = mnbode_list(CMS,Freq_CMS);
Input_Noise_CMS = Output_Noise_CMS ./ servo_c;
Noise_CMS = interp1(Freq_CMS,Input_Noise_CMS,freq);
for ii = 1:length(Noise_CMS)
    if isnan(Noise_CMS(ii))
        Noise_CMS(ii) = last;
    else
        last = Noise_CMS(ii);
    end
end
% loglog(freq,Noise_CMS)
%%

% close all
% mnloglog(Freq_CMS,Output_Noise_CMS,Freq_CMS2,Output_Noise_CMS2,Freq_CMS,Output_Noise_CMS,Freq_CMS3,Output_Noise_CMS3,'title','Output Noise',...
%     'legend',{'Measured at IOO rack','Measured in PSL room'},'xlabel','frequency[Hz]','ylabel','noise [V/rtHz]')
end