efb929c891fb73a5aaedfc3816c6df919ec4bf60.svn-base

function [dac,total,stageNames] = plothxdsactuatornoise(freq,susModel,driverConfigs,printFigs)

%%

susType   = 'HXDS';
svnDir.sus = '/ligo/svncommon/SusSVN/sus/trunk/';
svnDir.sei = '/ligo/svncommon/SeiSVN/seismic/';
svnDir.anb = '/ligo/svncommon/NbSVN/aligonoisebudget/trunk/';
toolsDir.sus = [svnDir.sus 'Common/MatlabTools/'];
toolsDir.anb = [svnDir.anb 'Dev/Utils/NoiseModel/'];
modelDir.sus = [svnDir.sus 'Common/MatlabTools/TripleModel_Production/'];
modelDir.isi = [svnDir.sei 'Common/MatlabTools/'];
actNoiseDir.sus = [svnDir.anb 'Dev/SusElectronics/matlab/NoiseData/'];

resultsDir = [svnDir.sus 'HXDS/Common/Data/'];

figFileTag = [datestr(now,29) '_' susType '_ActNoiseComp'];

% Modeled noises are
% (1) Top Driver (as originally designed)

stageNames.model   = {'M1'};
driverNames  = {'HAM-A Driver -v4 (D1100117)'};
legendNames.(stageNames.model{1}) = {'UL','LL','UR','LR'};

nOSEMs.M1 = 4;

nStates.M1 = 2;

stateNames.M1 = {'acq';'lp'};

             %L T V R P Y         
nActs.htts = [4 0 0 0 4 4]; % M1
              
                  
                %%L T V R    P     Y          % Distances in T2100275-v3
leverArms.htts = [1 0 0 0 40.64e-3 48.26e-3]; % [m] M1        
          
dofs = {'L';'T';'V';'R';'P';'Y'};

ctrlDOFs = [1 5 6];

%%

% Use when using as a function, comment out when running as a script
buildType = susModel.paramsName;
xLimits = [freq(1) freq(end)];
%%
switch driverConfigs
    case 'Acquire'
        cdState.M1 = 1;
    case 'LowNoise'
        cdState.M1 = 2;
        isoStage = 'M1';
end

stageNames.cdStateNum = cdState;
for iStage = 1
    stageNames.cdStateName.(stageNames.model{iStage}) = stateNames.(stageNames.model{iStage}){cdState.(stageNames.model{iStage})};
    
    addOn = '';
    
    stageNames.driverName.(stageNames.model{iStage}) = [addOn driverNames{iStage}(1:(find(driverNames{iStage} == ' ',1,'first')-1))];
end

figFileTag = [figFileTag '_' driverConfigs];

%%
addpath(toolsDir.sus);
addpath(modelDir.isi);
addpath(modelDir.sus);
addpath(genpath(toolsDir.anb));

%% ADC / DAC  
% DAC
dac.dcGain_Vpct = 2^16 / 20; % [ct / V]
dac.noise_VprtHz = SUS_actuator_noise('GS_16bit_DAC',freq);

% ADC 
adc.dcGain_ctpV = 2^16 / 40; % [ct / V]
adc.noise_VprtHz = SEI_sensor_noise('GS_16bit_ADC',freq); % [V/rtHz]


%% Get BSC and SUS Model
hamReqs = HAM_req(freq);

%% Create Models of DAC noise and Coil Driver Self Noise
for iStage = 1
    iState = cdState.(stageNames.model{iStage});
    % Get Coil Driver Response for ith Stage in ith State
    [~,~,osemModel,coilDriverModel,~,~] = make_OSEM_filter_model(freq,susType,isoStage,stateNames.(stageNames.model{iStage}){iState});
    stage(iStage).forceCoeff = osemModel.NperA;
    driverDCTrans_ApV = coilDriverModel.dcTransConductance;
    
    for iDOF = ctrlDOFs
        % Convert everthing to Force noise in Euler Basis
        NperA = stage(iStage).forceCoeff .* leverArms.htts(iStage,iDOF) * sqrt(nActs.htts(iStage,iDOF)) ;
        
        dac.stage(iStage).dof(iDOF).force   = dac.noise_VprtHz(:) .* driverDCTrans_ApV .* abs(coilDriverModel.f(:)) .* NperA;
       
        % FIXME
        % JUST ASSUME THAT DAC NOISE DOMINATES WAY ABOVE SELF NOISE OF THE
        % DRIVER FOR NOW.
        total.stage(iStage).dof(iDOF).force = dac.stage(iStage).dof(iDOF).force;
        
        % Convert everything to displacement of optic
        mperN = abs(squeeze(susModel.dampedf(susModel.dampedout.m1.disp.(dofs{iDOF}),susModel.dampedin.(lower(stageNames.model{iStage})).drive.(dofs{iDOF}),:)));
        
        dac.stage(iStage).dof(iDOF).disp   =   dac.stage(iStage).dof(iDOF).force .* mperN;
        
        total.stage(iStage).dof(iDOF).disp = dac.stage(iStage).dof(iDOF).disp;
        
        sei.dof(iDOF).disp = hamReqs(:) .* abs(squeeze(susModel.dampedf(susModel.dampedout.m1.disp.(dofs{iDOF}),susModel.dampedin.gnd.disp.(dofs{iDOF}),:)));
        dispReqs = SUS_reqs(freq,susType);
    end                            
end

%% Plot it all!
figNum = 1;

driveUnits = {'N/rtHz';...
              'N/rtHz';...
              'N/rtHz';...
              'N.m/rtHz';...
              'N.m/rtHz';...
              'N.m/rtHz'}; 
          
dispUnits = {'m/rtHz';...
             'm/rtHz';...
             'm/rtHz';...
             'rad/rtHz';...
             'rad/rtHz';...
             'rad/rtHz'};           
         
legLocation = {'NorthEast';...
               'SouthEast';...
               'NorthEast';...
               'SouthEast'}; 
           
theLegend = {[driverNames{1} ', Model CD Self Noise'],...
             [driverNames{1} ', Model DAC Noise'],...
             [susType ' ' stageNames.model{1} ' Stage LowNoise State Requirements']};
    
if printFigs         
    figNames(figNum) = {[figFileTag '_ICoil_' ...
                         stageNames.channel{1} 'State' num2str(cdState.(stageNames.channel{1})) '.pdf']};
    figure(figNum);figNum = figNum+1;
    ll=loglog(freq(:),[cd.stage(1).icoil(:),...
                       dac.stage(1).icoil(:)],...
              req.freq(:),[req.stage(1).icoil(:)]);
    title({[susType ' Coil Driver Noise Comparison: Current Noise across Coil'];...
           [stageNames.model{1} ': ' stateNames.(stageNames.channel{1}){cdState.(stageNames.channel{1})}]});
    set(ll([1]),'Color',[0.0 0.0 1.0])
    set(ll([2]),'Color',[0.0 0.7 0.0])
    set(ll([3]),'Color',[1.0 0.0 0.0])
    set(ll(1),'LineStyle','--')
    set(ll(3),'LineStyle','-.','LineWidth',5)
    leg = legend(theLegend);
    ylabel('Noise Amplitude [A/rtHz]')
    xlabel('Frequency [Hz]')
    set(leg,'FontSize',10)
    ylim([1e-12 1e-8])
    set(gca,'YTick',10.^(-12:-8))
    xlim(xLimits)
    
    for iDOF = 1:length(dofs)
        figNames(figNum) = {[figFileTag '_OpticForce_' stageNames.channel{1} 'State' num2str(cdState.(stageNames.channel{1})) '.pdf']};
        figure(figNum);figNum = figNum+1;
        ll=loglog(freq(:),[cd.stage(1).dof(iDOF).force(:),...
                           dac.stage(1).dof(iDOF).force(:)],...
                  req.freq(:),[req.stage(1).dof(iDOF).force(:)]);
        title({[susType ' Coil Driver Noise Comparison: Force Noise at Breadboard, ' dofs{iDOF} ];...
               [stageNames.model{1} ': ' stateNames.(stageNames.channel{1}){cdState.(stageNames.channel{1})}]});
        set(ll([1]),'Color',[0.0 0.0 1.0])
        set(ll([2]),'Color',[0.0 0.7 0.0])
        set(ll([3]),'Color',[1.0 0.0 0.0])
        set(ll(1),'LineStyle','--')
        set(ll(3),'LineStyle','-.','LineWidth',5)
        leg = legend(theLegend);
        ylabel(['Noise Amplitude [' driveUnits{iDOF} ']'])
        xlabel('Frequency [Hz]')
        set(leg,'FontSize',10)
        ylim([1e-14 1e-7])
        set(gca,'YTick',10.^(-14:-7))
        xlim(xLimits)
        
        figNames(figNum) = {[figFileTag '_OpticDispBudget_' stageNames.channel{1} 'State' num2str(cdState.(stageNames.channel{1})) '.pdf']};
        figure(figNum);figNum = figNum+1;
        ll=loglog(freq(:),[cd.stage(1).dof(iDOF).disp(:),...
                           dac.stage(1).dof(iDOF).disp(:)],...
                  req.freq(:),[req.stage(1).dof(iDOF).disp(:)],...
                  freq(:),sei.dof(iDOF).disp(:));
        title({[susType ' Coil Driver Noise Comparison: Displacement Noise at Breadboard, ' dofs{iDOF} ];...
               [stageNames.model{1} ': ' stateNames.(stageNames.channel{1}){cdState.(stageNames.channel{1})}]});
        set(ll([1]),'Color',[0.0 0.0 1.0])
        set(ll([2]),'Color',[0.0 0.7 0.0])
        set(ll([3]),'Color',[1.0 0.0 0.0])
        set(ll(1),'LineStyle','--')
        set(ll(3),'LineStyle','-.','LineWidth',5)
        set(ll(4),'LineWidth',4,'Color',[0.5 0.0 0.8])
        leg = legend(theLegend,...
                     'Location','SouthWest');
        ylabel(['Noise Amplitude [' dispUnits{iDOF} ']'])
        xlabel('Frequency [Hz]')
        set(leg,'FontSize',14)
        ylim([1e-20 1e-7])
        set(gca,'YTick',10.^(-20:-7))
        xlim(xLimits)
        
        figNames(figNum) = {[figFileTag '_OpticDispTotal_' stageNames.channel{1} 'State' num2str(cdState.(stageNames.channel{1})) '.pdf']};
        figure(figNum);figNum = figNum+1;
        ll=loglog(freq(:),[total.stage(1).dof(iDOF).disp(:)],...
                  req.freq(:),[req.stage(1).dof(iDOF).disp(:)],...
                  freq(:),sei.dof(iDOF).disp(:));
        title({[susType ' Coil Driver Noise Comparison: Displacement Noise at Breadboard, ' dofs{iDOF} ];...
               [stageNames.model{1} ': ' stateNames.(stageNames.channel{1}){cdState.(stageNames.channel{1})}]});
        set(ll([1]),'Color',[0.0 0.0 1.0])
        set(ll([2]),'Color',[1.0 0.0 0.0])
        set(ll(3),'LineWidth',4,'Color',[0.5 0.0 0.8])
        leg = legend([driverNames{1} ', Total Actuator Noise (DAC+CD)'],...
                     ['HAM-ISI Reqs -> ' susType ' Model'],'Location',legLocation{cdState.(stageNames.channel{3})},...
                     'Location','SouthWest');
        ylabel(['Noise Amplitude [' dispUnits{iDOF} ']'])
        xlabel('Frequency [Hz]')
        set(leg,'FontSize',14)
        ylim([1e-20 1e-7])
        set(gca,'YTick',10.^(-20:-6))
        xlim(xLimits)
    end
end
figNum = figNum-1;

%%
if printFigs
    mergeCommand = ['pdfmerge ' resultsDir figFileTag '.pdf'];

    
    for iFig = [1 2 5 8 11 14 17 3 6 9 12 15 18 4 7 10 13 16 19];
        figure(iFig)
        FillPage('w')
        IDfig(', J. Kissel')
        saveas(gcf,[resultsDir figNames{iFig}])
        
        mergeCommand = [mergeCommand ' ' resultsDir figNames{iFig}];
    end
    
    system(mergeCommand);
end