
% frequency bins
f = logspace(-1, 3.7, 1000);

disp('AdVirgo1 - Low power - Single demodulation - define interferometer');
opt1 = adVirgo1_sd_lp;
disp('AdVirgo1 - Low power - Single demodulation - run simulation');
[fDC1lp, sigDC1lp, sigAC1lp, mMech1lp, noiseAC1lp] = tickle(opt1, [], f);
disp('AdVirgo1 - Low power - Single demodulation - saving');
save advirgo1_sd_lp.mat f fDC1lp sigDC1lp sigAC1lp mMech1lp noiseAC1lp

disp('AdVirgo1 - High power - Single demodulation - define interferometer');
opt2 = adVirgo1_sd_hp;
disp('AdVirgo1 - High power - Single demodulation - run simulation');
[fDC1hp, sigDC1hp, sigAC1hp, mMech1hp, noiseAC1hp] = tickle(opt2, [], f);
disp('AdVirgo1 - High power - Single demodulation - saving');
save advirgo1_sd_hp.mat f fDC1hp sigDC1hp sigAC1hp mMech1hp noiseAC1hp

disp('AdVirgo1 - Low power - Double demodulation - define interferometer');
opt1 = adVirgo1_dd_lp;
disp('AdVirgo1 - Low power - Double demodulation - run simulation');
[fDC1lp, sigDC1lp, sigAC1lp, mMech1lp, noiseAC1lp] = tickle(opt1, [], f);
disp('AdVirgo1 - Low power - Double demodulation - saving');
save advirgo1_dd_lp.mat f fDC1lp sigDC1lp sigAC1lp mMech1lp noiseAC1lp

disp('AdVirgo1 - High power - Double demodulation - define interferometer');
opt2 = adVirgo1_dd_hp;
disp('AdVirgo1 - High power - Double demodulation - run simulation');
[fDC1hp, sigDC1hp, sigAC1hp, mMech1hp, noiseAC1hp] = tickle(opt2, [], f);
disp('AdVirgo1 - High power - Double demodulation - saving');
save advirgo1_dd_hp.mat f fDC1hp sigDC1hp sigAC1hp mMech1hp noiseAC1hp

extractTF1_sd_lp
extractTF1_sd_hp
extractTF1_dd_lp
extractTF1_dd_hp

tune_phases1_sd_lp
tune_phases1_sd_hp
tune_phases1_dd_lp
tune_phases1_dd_hp

save_text1_sd_lp
save_text1_sd_hp
save_text1_dd_lp
save_text1_dd_hp

noise1_sd_lp;
noise1_sd_hp;
noise1_dd_lp;
noise1_dd_hp;
