function runbKAGRAp1_LSC()
% Sweep plots and sensing matrix for bKAGRA Phase 1 LSC
% uses Optickle 2 and Tools directory in KAGRAOptickle SVN
% needs (1) parambKAGRAp1.m          for parameter setting
%       (2) optbKAGRAp1.m            for constructing Optickle model
%       (3) plotSweep.m              for sweep plots (lives in Tools directory)
%       (4) showSensingMatrix.m      for showing sensing matrix (lives in Tools directory)
% Initiated by: Yuta Michimura

%% Clear everything
close all
clear all

%% Set parameters and save directory
% set parameters
p = parambKAGRAp1();

% set results directory
saveDir=['./results/',p.modelName,'-',datestr(now,'yyyy-mm-dd'),'/'];
mkdir(saveDir);

%% Construct optickle model and set some parameters, lists
opt = optbKAGRAp1(p);
opt2dot(opt, [saveDir,p.modelName,'.pdf']); % GraphViz needs to be installed

% some parameters
lambda = opt.lambda(1);

% DOF definitions
probeNames={'REFL_I', 'REFL_Q','AS_I', 'AS_Q'};  % probes of interest
driveNames = {'ETMX','ETMY','BS','PRM'};  % drives of interest
dofNames = {'MICH','PRCL'};   % DOF of interest
mDrv = [1,-1,0,0; % MICH  Conversion matrix for optic motion to DOF motion
        0,0,0,1]; % PRCL
    
%% Check power at each field, each probe
[fDC, sigDC] = opt.tickle([], []);
fprintf('DC fields (fDC matrix):\n');
showfDC(opt, fDC);
fprintf('\nProbes (sigDC matrix):\n');
showsigDC(opt, sigDC);

%% Sweep in length
% PRCL sweep
sweepDOF = 'PRCL';
range = lambda;
plotSweep(opt,sweepDOF,probeNames,range,driveNames,dofNames,mDrv)
saveas(gcf,[saveDir,sweepDOF,'sweep.png'])

% MICH sweep
sweepDOF = 'MICH';
range = lambda;
plotSweep(opt,sweepDOF,probeNames,range,driveNames,dofNames,mDrv)
saveas(gcf,[saveDir,sweepDOF,'sweep.png'])

%% Sensing Matrix for LSC
[fDC, sigDC, sigAC] = opt.tickle([], p.ftickle);
fprintf('\nLSC Sensing Matrix [W/m] for %s:\n',p.modelName);
diary([saveDir,'/LSCSensingMatrix.txt'])
diary on
showSensingMatrix(opt,sigAC,probeNames,driveNames,dofNames,mDrv);
%showSensingMatrix(opt,sigAC,probeNames,driveNames)
%showSensingMatrix(opt,sigAC,{opt.probe.name},driveNames)
%showSensingMatrix(opt,sigAC)
diary off