function [noise, f] = noise1_sd_hp()

% load data
clear
load advirgo1_sd_hp.mat
opt = adVirgo1_sd_hp;

% get photo-diodes numbers
nXP_SB1_P   = getProbeNum(opt, 'XP_SB1_P');
nXP_SB1_Q   = getProbeNum(opt, 'XP_SB1_Q');
nXP_SB2_P   = getProbeNum(opt, 'XP_SB2_P');
nXP_SB2_Q   = getProbeNum(opt, 'XP_SB2_Q');
nXP_SB3_P   = getProbeNum(opt, 'XP_SB3_P');
nXP_SB3_Q   = getProbeNum(opt, 'XP_SB3_Q');

nYP_SB1_P   = getProbeNum(opt, 'YP_SB1_P');
nYP_SB1_Q   = getProbeNum(opt, 'YP_SB1_Q');
nYP_SB2_P   = getProbeNum(opt, 'YP_SB2_P');
nYP_SB2_Q   = getProbeNum(opt, 'YP_SB2_Q');
nYP_SB3_P   = getProbeNum(opt, 'YP_SB3_P');
nYP_SB3_Q   = getProbeNum(opt, 'YP_SB3_Q');

nSYM_SB1_P   = getProbeNum(opt, 'SYM_SB1_P');
nSYM_SB1_Q   = getProbeNum(opt, 'SYM_SB1_Q');
nSYM_SB2_P   = getProbeNum(opt, 'SYM_SB2_P');
nSYM_SB2_Q   = getProbeNum(opt, 'SYM_SB2_Q');
nSYM_SB3_P   = getProbeNum(opt, 'SYM_SB3_P');
nSYM_SB3_Q   = getProbeNum(opt, 'SYM_SB3_Q');

nASYp_SB1_P   = getProbeNum(opt, 'ASYp_SB1_P');
nASYp_SB1_Q   = getProbeNum(opt, 'ASYp_SB1_Q');
nASYp_SB2_P   = getProbeNum(opt, 'ASYp_SB2_P');
nASYp_SB2_Q   = getProbeNum(opt, 'ASYp_SB2_Q');
nASYp_SB3_P   = getProbeNum(opt, 'ASYp_SB3_P');
nASYp_SB3_Q   = getProbeNum(opt, 'ASYp_SB3_Q');

nPOBS_SB1_P   = getProbeNum(opt, 'POBS_SB1_P');
nPOBS_SB1_Q   = getProbeNum(opt, 'POBS_SB1_Q');
nPOBS_SB2_P   = getProbeNum(opt, 'POBS_SB2_P');
nPOBS_SB2_Q   = getProbeNum(opt, 'POBS_SB2_Q');
nPOBS_SB3_P   = getProbeNum(opt, 'POBS_SB3_P');
nPOBS_SB3_Q   = getProbeNum(opt, 'POBS_SB3_Q');

nPOX_SB1_P   = getProbeNum(opt, 'POX_SB1_P');
nPOX_SB1_Q   = getProbeNum(opt, 'POX_SB1_Q');
nPOX_SB2_P   = getProbeNum(opt, 'POX_SB2_P');
nPOX_SB2_Q   = getProbeNum(opt, 'POX_SB2_Q');
nPOX_SB3_P   = getProbeNum(opt, 'POX_SB3_P');
nPOX_SB3_Q   = getProbeNum(opt, 'POX_SB3_Q');

nPOY_SB1_P   = getProbeNum(opt, 'POY_SB1_P');
nPOY_SB1_Q   = getProbeNum(opt, 'POY_SB1_Q');
nPOY_SB2_P   = getProbeNum(opt, 'POY_SB2_P');
nPOY_SB2_Q   = getProbeNum(opt, 'POY_SB2_Q');
nPOY_SB3_P   = getProbeNum(opt, 'POY_SB3_P');
nPOY_SB3_Q   = getProbeNum(opt, 'POY_SB3_Q');

nASY_DC      = getProbeNum(opt, 'ASY_DC');

% extract quantum noise (W)

NOISE_ASYp_SB1_P = noiseAC1hp(nASYp_SB1_P, :);
NOISE_ASYp_SB1_Q = noiseAC1hp(nASYp_SB1_Q, :);
NOISE_ASYp_SB2_P = noiseAC1hp(nASYp_SB2_P, :);
NOISE_ASYp_SB2_Q = noiseAC1hp(nASYp_SB2_Q, :);
NOISE_ASYp_SB3_P = noiseAC1hp(nASYp_SB3_P, :);
NOISE_ASYp_SB3_Q = noiseAC1hp(nASYp_SB3_Q, :);

NOISE_SYM_SB1_P = noiseAC1hp(nSYM_SB1_P, :);
NOISE_SYM_SB1_Q = noiseAC1hp(nSYM_SB1_Q, :);
NOISE_SYM_SB2_P = noiseAC1hp(nSYM_SB2_P, :);
NOISE_SYM_SB2_Q = noiseAC1hp(nSYM_SB2_Q, :);
NOISE_SYM_SB3_P = noiseAC1hp(nSYM_SB3_P, :);
NOISE_SYM_SB3_Q = noiseAC1hp(nSYM_SB3_Q, :);

NOISE_POBS_SB1_P = noiseAC1hp(nPOBS_SB1_P, :);
NOISE_POBS_SB1_Q = noiseAC1hp(nPOBS_SB1_Q, :);
NOISE_POBS_SB2_P = noiseAC1hp(nPOBS_SB2_P, :);
NOISE_POBS_SB2_Q = noiseAC1hp(nPOBS_SB2_Q, :);
NOISE_POBS_SB3_P = noiseAC1hp(nPOBS_SB3_P, :);
NOISE_POBS_SB3_Q = noiseAC1hp(nPOBS_SB3_Q, :);

NOISE_POX_SB1_P = noiseAC1hp(nPOX_SB1_P, :);
NOISE_POX_SB1_Q = noiseAC1hp(nPOX_SB1_Q, :);
NOISE_POX_SB2_P = noiseAC1hp(nPOX_SB2_P, :);
NOISE_POX_SB2_Q = noiseAC1hp(nPOX_SB2_Q, :);
NOISE_POX_SB3_P = noiseAC1hp(nPOX_SB3_P, :);
NOISE_POX_SB3_Q = noiseAC1hp(nPOX_SB3_Q, :);

NOISE_POY_SB1_P = noiseAC1hp(nPOY_SB1_P, :);
NOISE_POY_SB1_Q = noiseAC1hp(nPOY_SB1_Q, :);
NOISE_POY_SB2_P = noiseAC1hp(nPOY_SB2_P, :);
NOISE_POY_SB2_Q = noiseAC1hp(nPOY_SB2_Q, :);
NOISE_POY_SB3_P = noiseAC1hp(nPOY_SB3_P, :);
NOISE_POY_SB3_Q = noiseAC1hp(nPOY_SB3_Q, :);

NOISE_XP_SB1_P = noiseAC1hp(nXP_SB1_P, :);
NOISE_XP_SB1_Q = noiseAC1hp(nXP_SB1_Q, :);
NOISE_XP_SB2_P = noiseAC1hp(nXP_SB2_P, :);
NOISE_XP_SB2_Q = noiseAC1hp(nXP_SB2_Q, :);
NOISE_XP_SB3_P = noiseAC1hp(nXP_SB3_P, :);
NOISE_XP_SB3_Q = noiseAC1hp(nXP_SB3_Q, :);

NOISE_YP_SB1_P = noiseAC1hp(nYP_SB1_P, :);
NOISE_YP_SB1_Q = noiseAC1hp(nYP_SB1_Q, :);
NOISE_YP_SB2_P = noiseAC1hp(nYP_SB2_P, :);
NOISE_YP_SB2_Q = noiseAC1hp(nYP_SB2_Q, :);
NOISE_YP_SB3_P = noiseAC1hp(nYP_SB3_P, :);
NOISE_YP_SB3_Q = noiseAC1hp(nYP_SB3_Q, :);

NOISE_ASY_DC = noiseAC1hp(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_SB2_P';
noise(:, 5) = NOISE_ASYp_SB2_Q';
noise(:, 6) = NOISE_ASYp_SB3_P';
noise(:, 7) = NOISE_ASYp_SB3_Q';
	 
noise(:, 8) = NOISE_SYM_SB1_P';
noise(:, 9) = NOISE_SYM_SB1_Q';
noise(:,10) = NOISE_SYM_SB2_P';
noise(:,11) = NOISE_SYM_SB2_Q';
noise(:,12) = NOISE_SYM_SB3_P';
noise(:,13) = NOISE_SYM_SB3_Q';

noise(:,14) = NOISE_POBS_SB1_P';
noise(:,15) = NOISE_POBS_SB1_Q';
noise(:,16) = NOISE_POBS_SB2_P';
noise(:,17) = NOISE_POBS_SB2_Q';
noise(:,18) = NOISE_POBS_SB3_P';
noise(:,19) = NOISE_POBS_SB3_Q';

noise(:,20) = NOISE_POX_SB1_P';
noise(:,21) = NOISE_POX_SB1_Q';
noise(:,22) = NOISE_POX_SB2_P';
noise(:,23) = NOISE_POX_SB2_Q';
noise(:,24) = NOISE_POX_SB3_P';
noise(:,25) = NOISE_POX_SB3_Q';

noise(:,26) = NOISE_POY_SB1_P';
noise(:,27) = NOISE_POY_SB1_Q';
noise(:,28) = NOISE_POY_SB2_P';
noise(:,29) = NOISE_POY_SB2_Q';
noise(:,30) = NOISE_POY_SB3_P';
noise(:,31) = NOISE_POY_SB3_Q';

noise(:,32) = NOISE_XP_SB1_P';
noise(:,33) = NOISE_XP_SB1_Q';
noise(:,34) = NOISE_XP_SB2_P';
noise(:,35) = NOISE_XP_SB2_Q';
noise(:,36) = NOISE_XP_SB3_P';
noise(:,37) = NOISE_XP_SB3_Q';

noise(:,38) = NOISE_YP_SB1_P';
noise(:,39) = NOISE_YP_SB1_Q';
noise(:,40) = NOISE_YP_SB2_P';
noise(:,41) = NOISE_YP_SB2_Q';
noise(:,42) = NOISE_YP_SB3_P';
noise(:,43) = NOISE_YP_SB3_Q';

save advirgo1_noise_sd_hp.txt noise -ASCII
