function [noise, f] = noise1_dd_lp()

% load data
clear
load advirgo1_dd_lp.mat
opt = adVirgo1_dd_lp;

% get photo-diodes numbers
nXP_SB1_P   = getProbeNum(opt, 'XP_SB1_P');
nXP_SB1_Q   = getProbeNum(opt, 'XP_SB1_Q');
nXP_2M1_P   = getProbeNum(opt, 'XP_2M1_P');
nXP_2M1_Q   = getProbeNum(opt, 'XP_2M1_Q');
nXP_3P1_P   = getProbeNum(opt, 'XP_3P1_P');
nXP_3P1_Q   = getProbeNum(opt, 'XP_3P1_Q');

nYP_SB1_P   = getProbeNum(opt, 'YP_SB1_P');
nYP_SB1_Q   = getProbeNum(opt, 'YP_SB1_Q');
nYP_2M1_P   = getProbeNum(opt, 'YP_2M1_P');
nYP_2M1_Q   = getProbeNum(opt, 'YP_2M1_Q');
nYP_3P1_P   = getProbeNum(opt, 'YP_3P1_P');
nYP_3P1_Q   = getProbeNum(opt, 'YP_3P1_Q');

nSYM_SB1_P   = getProbeNum(opt, 'SYM_SB1_P');
nSYM_SB1_Q   = getProbeNum(opt, 'SYM_SB1_Q');
nSYM_2M1_P   = getProbeNum(opt, 'SYM_2M1_P');
nSYM_2M1_Q   = getProbeNum(opt, 'SYM_2M1_Q');
nSYM_3P1_P   = getProbeNum(opt, 'SYM_3P1_P');
nSYM_3P1_Q   = getProbeNum(opt, 'SYM_3P1_Q');

nASYp_SB1_P   = getProbeNum(opt, 'ASYp_SB1_P');
nASYp_SB1_Q   = getProbeNum(opt, 'ASYp_SB1_Q');
nASYp_2M1_P   = getProbeNum(opt, 'ASYp_2M1_P');
nASYp_2M1_Q   = getProbeNum(opt, 'ASYp_2M1_Q');
nASYp_3P1_P   = getProbeNum(opt, 'ASYp_3P1_P');
nASYp_3P1_Q   = getProbeNum(opt, 'ASYp_3P1_Q');

nPOBS_SB1_P   = getProbeNum(opt, 'POBS_SB1_P');
nPOBS_SB1_Q   = getProbeNum(opt, 'POBS_SB1_Q');
nPOBS_2M1_P   = getProbeNum(opt, 'POBS_2M1_P');
nPOBS_2M1_Q   = getProbeNum(opt, 'POBS_2M1_Q');
nPOBS_3P1_P   = getProbeNum(opt, 'POBS_3P1_P');
nPOBS_3P1_Q   = getProbeNum(opt, 'POBS_3P1_Q');

nPOX_SB1_P   = getProbeNum(opt, 'POX_SB1_P');
nPOX_SB1_Q   = getProbeNum(opt, 'POX_SB1_Q');
nPOX_2M1_P   = getProbeNum(opt, 'POX_2M1_P');
nPOX_2M1_Q   = getProbeNum(opt, 'POX_2M1_Q');
nPOX_3P1_P   = getProbeNum(opt, 'POX_3P1_P');
nPOX_3P1_Q   = getProbeNum(opt, 'POX_3P1_Q');

nPOY_SB1_P   = getProbeNum(opt, 'POY_SB1_P');
nPOY_SB1_Q   = getProbeNum(opt, 'POY_SB1_Q');
nPOY_2M1_P   = getProbeNum(opt, 'POY_2M1_P');
nPOY_2M1_Q   = getProbeNum(opt, 'POY_2M1_Q');
nPOY_3P1_P   = getProbeNum(opt, 'POY_3P1_P');
nPOY_3P1_Q   = getProbeNum(opt, 'POY_3P1_Q');

nASY_DC      = getProbeNum(opt, 'ASY_DC');

% extract quantum noise (W)

NOISE_ASYp_SB1_P = noiseAC1lp(nASYp_SB1_P, :);
NOISE_ASYp_SB1_Q = noiseAC1lp(nASYp_SB1_Q, :);
NOISE_ASYp_2M1_P = noiseAC1lp(nASYp_2M1_P, :);
NOISE_ASYp_2M1_Q = noiseAC1lp(nASYp_2M1_Q, :);
NOISE_ASYp_3P1_P = noiseAC1lp(nASYp_3P1_P, :);
NOISE_ASYp_3P1_Q = noiseAC1lp(nASYp_3P1_Q, :);

NOISE_SYM_SB1_P = noiseAC1lp(nSYM_SB1_P, :);
NOISE_SYM_SB1_Q = noiseAC1lp(nSYM_SB1_Q, :);
NOISE_SYM_2M1_P = noiseAC1lp(nSYM_2M1_P, :);
NOISE_SYM_2M1_Q = noiseAC1lp(nSYM_2M1_Q, :);
NOISE_SYM_3P1_P = noiseAC1lp(nSYM_3P1_P, :);
NOISE_SYM_3P1_Q = noiseAC1lp(nSYM_3P1_Q, :);

NOISE_POBS_SB1_P = noiseAC1lp(nPOBS_SB1_P, :);
NOISE_POBS_SB1_Q = noiseAC1lp(nPOBS_SB1_Q, :);
NOISE_POBS_2M1_P = noiseAC1lp(nPOBS_2M1_P, :);
NOISE_POBS_2M1_Q = noiseAC1lp(nPOBS_2M1_Q, :);
NOISE_POBS_3P1_P = noiseAC1lp(nPOBS_3P1_P, :);
NOISE_POBS_3P1_Q = noiseAC1lp(nPOBS_3P1_Q, :);

NOISE_POX_SB1_P = noiseAC1lp(nPOX_SB1_P, :);
NOISE_POX_SB1_Q = noiseAC1lp(nPOX_SB1_Q, :);
NOISE_POX_2M1_P = noiseAC1lp(nPOX_2M1_P, :);
NOISE_POX_2M1_Q = noiseAC1lp(nPOX_2M1_Q, :);
NOISE_POX_3P1_P = noiseAC1lp(nPOX_3P1_P, :);
NOISE_POX_3P1_Q = noiseAC1lp(nPOX_3P1_Q, :);

NOISE_POY_SB1_P = noiseAC1lp(nPOY_SB1_P, :);
NOISE_POY_SB1_Q = noiseAC1lp(nPOY_SB1_Q, :);
NOISE_POY_2M1_P = noiseAC1lp(nPOY_2M1_P, :);
NOISE_POY_2M1_Q = noiseAC1lp(nPOY_2M1_Q, :);
NOISE_POY_3P1_P = noiseAC1lp(nPOY_3P1_P, :);
NOISE_POY_3P1_Q = noiseAC1lp(nPOY_3P1_Q, :);

NOISE_XP_SB1_P = noiseAC1lp(nXP_SB1_P, :);
NOISE_XP_SB1_Q = noiseAC1lp(nXP_SB1_Q, :);
NOISE_XP_2M1_P = noiseAC1lp(nXP_2M1_P, :);
NOISE_XP_2M1_Q = noiseAC1lp(nXP_2M1_Q, :);
NOISE_XP_3P1_P = noiseAC1lp(nXP_3P1_P, :);
NOISE_XP_3P1_Q = noiseAC1lp(nXP_3P1_Q, :);

NOISE_YP_SB1_P = noiseAC1lp(nYP_SB1_P, :);
NOISE_YP_SB1_Q = noiseAC1lp(nYP_SB1_Q, :);
NOISE_YP_2M1_P = noiseAC1lp(nYP_2M1_P, :);
NOISE_YP_2M1_Q = noiseAC1lp(nYP_2M1_Q, :);
NOISE_YP_3P1_P = noiseAC1lp(nYP_3P1_P, :);
NOISE_YP_3P1_Q = noiseAC1lp(nYP_3P1_Q, :);

NOISE_ASY_DC = noiseAC1lp(nASY_DC, :);

% save all in a single table

noise(:, 1) = NOISE_ASY_DC';
	 
noise(:, 2) = NOISE_ASYp_SB1_P';
noise(:, 3) = NOISE_ASYp_SB1_Q';
noise(:, 4) = NOISE_ASYp_2M1_P';
noise(:, 5) = NOISE_ASYp_2M1_Q';
noise(:, 6) = NOISE_ASYp_3P1_P';
noise(:, 7) = NOISE_ASYp_3P1_Q';
	 
noise(:, 8) = NOISE_SYM_SB1_P';
noise(:, 9) = NOISE_SYM_SB1_Q';
noise(:,10) = NOISE_SYM_2M1_P';
noise(:,11) = NOISE_SYM_2M1_Q';
noise(:,12) = NOISE_SYM_3P1_P';
noise(:,13) = NOISE_SYM_3P1_Q';

noise(:,14) = NOISE_POBS_SB1_P';
noise(:,15) = NOISE_POBS_SB1_Q';
noise(:,16) = NOISE_POBS_2M1_P';
noise(:,17) = NOISE_POBS_2M1_Q';
noise(:,18) = NOISE_POBS_3P1_P';
noise(:,19) = NOISE_POBS_3P1_Q';

noise(:,20) = NOISE_POX_SB1_P';
noise(:,21) = NOISE_POX_SB1_Q';
noise(:,22) = NOISE_POX_2M1_P';
noise(:,23) = NOISE_POX_2M1_Q';
noise(:,24) = NOISE_POX_3P1_P';
noise(:,25) = NOISE_POX_3P1_Q';

noise(:,26) = NOISE_POY_SB1_P';
noise(:,27) = NOISE_POY_SB1_Q';
noise(:,28) = NOISE_POY_2M1_P';
noise(:,29) = NOISE_POY_2M1_Q';
noise(:,30) = NOISE_POY_3P1_P';
noise(:,31) = NOISE_POY_3P1_Q';

noise(:,32) = NOISE_XP_SB1_P';
noise(:,33) = NOISE_XP_SB1_Q';
noise(:,34) = NOISE_XP_2M1_P';
noise(:,35) = NOISE_XP_2M1_Q';
noise(:,36) = NOISE_XP_3P1_P';
noise(:,37) = NOISE_XP_3P1_Q';

noise(:,38) = NOISE_YP_SB1_P';
noise(:,39) = NOISE_YP_SB1_Q';
noise(:,40) = NOISE_YP_2M1_P';
noise(:,41) = NOISE_YP_2M1_Q';
noise(:,42) = NOISE_YP_3P1_P';
noise(:,43) = NOISE_YP_3P1_Q';

save advirgo1_noise_dd_lp.txt noise -ASCII
