addpath(genpath('../../Optickle'));
addpath(genpath('eLIGOcfg'));
addpath(genpath('LCGTcfg'));
addpath(genpath('lib_060504b'));
addpath(genpath('pickleCode'));
%Open used files

% edit DARMcouplBudget12DATA.m
% edit shotPropDATA.m
% edit pickleEligo.m
% edit paramEligo.m
% edit paramPower.m
% edit WFS1.m
% edit analGouyAllCOOL.m

% Configure LCGT

% (Po, Rprm, Rpr2, Rpr3, Rsrm, Rsr2, Rsr3)
par = paramPowerLCGT(82, 300.624, -3.251, 27.36, 300.624, -3.251, 27.36);
%par = paramPowerLCGT(82, 300.624, -3.251, 20, 300.624, -3.251, 20);
par = paramLCGT(par); %uses param_null for the estimate function
opt = optLCGT(par);
%opt = probesEligo_01(opt, par); A and B 90dg apart
opt = probeSensLCGT(opt, par);

setupLCGT;
pickle = pickleLCGT(opt);

% Frequency space

nP = 1000;
Ndof = pickle.param.Ndof;
Nmirr = pickle.param.Nmirr;

f = logspace(-1, 2.5, nP);

%TO BE DONE EVERY TIME THERE IS A CHANGE IN THE PLANT
[fDC, sigDC] = tickle(opt, [], []);

fprintf('REFL A and B power is %g and %g\n', sigDC(nREFLA_DC), sigDC(nREFLB_DC))
fprintf('POP A and B power is %g and %g\n', sigDC(nPOPA_DC), sigDC(nPOPB_DC))
fprintf('POX A and B power is %g and %g\n', sigDC(nPOXA_DC), sigDC(nPOXB_DC))
fprintf('AS A and B power is %g and %g\n', sigDC(nASA_DC), sigDC(nASB_DC))
fprintf('OMCtr A and B power is %g and %g\n', sigDC(nOMCtA_DC), sigDC(nOMCtB_DC))
fprintf('OMCref A and B power is %g and %g\n', sigDC(nOMCrA_DC), sigDC(nOMCrB_DC))
fprintf('TRX power is %g\n', sigDC(nTRX_DC))
fprintf('TRY power is %g\n', sigDC(nTRY_DC))
fprintf('PRC power is %g\n', sigDC(nPRM_DC))
fprintf('SRC power is %g\n', sigDC(nSRM_DC))
fprintf('Intracavity power of arms is %g and %g\n', sigDC(nIX_DC), sigDC(nIY_DC))

nl.REFL = getLinkNum(opt, 'IN2', 'AttREFL');
P_REFL = abs(fDC(nl.REFL, :)*conj(fDC(nl.REFL, :))');
fprintf('REFL power is %g\n', P_REFL);

nl.PRC = getLinkNum(opt, 'PRM', 'PR2');
P_PRC = abs(fDC(nl.PRC, par.Laser.vFrf==0))^2;
P_PRCs1 = abs(fDC(nl.PRC, par.Laser.vFrf == par.Mod.f1))^2;
P_PRCs2 = abs(fDC(nl.PRC, par.Laser.vFrf == par.Mod.f2))^2;
fprintf('Carrier power at PRC is %g\n', P_PRC);
fprintf('Sideband power f1 and f2 at PRC is %g an %g\n', P_PRCs1, P_PRCs2);

nl.armX = getLinkNum(opt, 'IX', 'EX');
nl.armY = getLinkNum(opt, 'IY', 'EY');
fprintf('Intracavity power of arms is %g and %g\n', abs(fDC(nl.armX, par.Laser.vFrf==0)).^2, abs(fDC(nl.armY, par.Laser.vFrf==0)).^2)

%save DATA/tickle_LASTnewConfig fDC sigDC
% Load old data instead of running tickle01

%[sigAC, mMech] = tickle01(opt, [], f);
%save DATA/tickle01_LASTnewConfig sigAC mMech
%save C:\Users\Agatsuma\Documents\MATLAB\Pickle_June2010\Pickle_Last\DATA/tickle01_LCGT_G90d0 sigAC mMech
%save C:\Users\Agatsuma\Documents\MATLAB\Pickle_June2010\Pickle_Last\DATA/tickle01_LCGT_110222 sigAC mMech


%load('tickle01_LCGT_110222');
%plotSensingMatrixASC_bLCGT(opt, sigAC);


% Useful constant values

hP = opt.h; % Planck's constant
e_charge = 1.60217646e-19;  % Coulombs
c_light = opt.c;
l_lambda = opt.lambda;
v_light = c_light / l_lambda;