% frequency bins
clear all
%f = logspace(log10(5), log10(2000), 100);
f = logspace(log10(5), log10(5000), 100);

lcgt
clc

disp('LCGT - Single demodulation - define interferometer');
opt2 = lcgt;
disp('LCGT - Single demodulation - run simulation');
[fDC1hp, sigDC1hp, sigAC1hp, mMech1hp, noiseAC1hp] = tickle(opt2, [], f);
%save f fDC1hp sigDC1hp sigAC1hp mMech1hp noiseAC1hp

extractTF
tune_phases
Sensitivity
IR