from scipy import *
from pylab import *

# load
darm = load('opt1_dd_hp_darm.txt', comments='%')
carm = load('opt1_dd_hp_carm.txt', comments='%')
prcl = load('opt1_dd_hp_prcl.txt', comments='%')
mich = load('opt1_dd_hp_mich.txt', comments='%')
srcl = load('opt1_dd_hp_srcl.txt', comments='%')

f = darm[:,0]


# PLOTS BY SIGNALS

# ASYp SB1
figure()
loglog(f, darm[:,13]*1e6, 'r-')
loglog(f, carm[:,13]*1e6, 'g-')
loglog(f, prcl[:,13]*1e6, 'b-')
loglog(f, mich[:,13]*1e6, 'c-')
loglog(f, srcl[:,13]*1e6, 'k-')
loglog(f, darm[:,15]*1e6, 'r--')
loglog(f, carm[:,15]*1e6, 'g--')
loglog(f, prcl[:,15]*1e6, 'b--')
loglog(f, mich[:,15]*1e6, 'c--')
loglog(f, srcl[:,15]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('ASYp SB1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_asyp_sb1.eps')
savefig('plots/advirgo1_dd_hp_asyp_sb1.png')

# ASYp 2M1
figure()
loglog(f, darm[:,17]*1e6, 'r-')
loglog(f, carm[:,17]*1e6, 'g-')
loglog(f, prcl[:,17]*1e6, 'b-')
loglog(f, mich[:,17]*1e6, 'c-')
loglog(f, srcl[:,17]*1e6, 'k-')
loglog(f, darm[:,19]*1e6, 'r--')
loglog(f, carm[:,19]*1e6, 'g--')
loglog(f, prcl[:,19]*1e6, 'b--')
loglog(f, mich[:,19]*1e6, 'c--')
loglog(f, srcl[:,19]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('ASYp 2M1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_asyp_2m1.eps')
savefig('plots/advirgo1_dd_hp_asyp_2m1.png')

# ASYp 3P1
figure()
loglog(f, darm[:,21]*1e6, 'r-')
loglog(f, carm[:,21]*1e6, 'g-')
loglog(f, prcl[:,21]*1e6, 'b-')
loglog(f, mich[:,21]*1e6, 'c-')
loglog(f, srcl[:,21]*1e6, 'k-')
loglog(f, darm[:,23]*1e6, 'r--')
loglog(f, carm[:,23]*1e6, 'g--')
loglog(f, prcl[:,23]*1e6, 'b--')
loglog(f, mich[:,23]*1e6, 'c--')
loglog(f, srcl[:,23]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('ASYp 3P1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_asyp_3p1.eps')
savefig('plots/advirgo1_dd_hp_asyp_3p1.png')

# SYM SB1
figure()
loglog(f, darm[:,1]*1e6, 'r-')
loglog(f, carm[:,1]*1e6, 'g-')
loglog(f, prcl[:,1]*1e6, 'b-')
loglog(f, mich[:,1]*1e6, 'c-')
loglog(f, srcl[:,1]*1e6, 'k-')
loglog(f, darm[:,3]*1e6, 'r--')
loglog(f, carm[:,3]*1e6, 'g--')
loglog(f, prcl[:,3]*1e6, 'b--')
loglog(f, mich[:,3]*1e6, 'c--')
loglog(f, srcl[:,3]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('SYM SB1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_sym_sb1.eps')
savefig('plots/advirgo1_dd_hp_sym_sb1.png')

# SYM 2M1
figure()
loglog(f, darm[:,5]*1e6, 'r-')
loglog(f, carm[:,5]*1e6, 'g-')
loglog(f, prcl[:,5]*1e6, 'b-')
loglog(f, mich[:,5]*1e6, 'c-')
loglog(f, srcl[:,5]*1e6, 'k-')
loglog(f, darm[:,7]*1e6, 'r--')
loglog(f, carm[:,7]*1e6, 'g--')
loglog(f, prcl[:,7]*1e6, 'b--')
loglog(f, mich[:,7]*1e6, 'c--')
loglog(f, srcl[:,7]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('SYM 2M1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_sym_2m1.eps')
savefig('plots/advirgo1_dd_hp_sym_2m1.png')

# SYM 3P1
figure()
loglog(f,  darm[:,9]*1e6, 'r-')
loglog(f,  carm[:,9]*1e6, 'g-')
loglog(f,  prcl[:,9]*1e6, 'b-')
loglog(f,  mich[:,9]*1e6, 'c-')
loglog(f,  srcl[:,9]*1e6, 'k-')
loglog(f, darm[:,11]*1e6, 'r--')
loglog(f, carm[:,11]*1e6, 'g--')
loglog(f, prcl[:,11]*1e6, 'b--')
loglog(f, mich[:,11]*1e6, 'c--')
loglog(f, srcl[:,11]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('SYM 3P1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_sym_3p1.eps')
savefig('plots/advirgo1_dd_hp_sym_3p1.png')

# POBS SB1
figure()
loglog(f, darm[:,25]*1e6, 'r-')
loglog(f, carm[:,25]*1e6, 'g-')
loglog(f, prcl[:,25]*1e6, 'b-')
loglog(f, mich[:,25]*1e6, 'c-')
loglog(f, srcl[:,25]*1e6, 'k-')
loglog(f, darm[:,27]*1e6, 'r--')
loglog(f, carm[:,27]*1e6, 'g--')
loglog(f, prcl[:,27]*1e6, 'b--')
loglog(f, mich[:,27]*1e6, 'c--')
loglog(f, srcl[:,27]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('POBS SB1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_pobs_sb1.eps')
savefig('plots/advirgo1_dd_hp_pobs_sb1.png')

# POBS 2M1
figure()
loglog(f, darm[:,29]*1e6, 'r-')
loglog(f, carm[:,29]*1e6, 'g-')
loglog(f, prcl[:,29]*1e6, 'b-')
loglog(f, mich[:,29]*1e6, 'c-')
loglog(f, srcl[:,29]*1e6, 'k-')
loglog(f, darm[:,31]*1e6, 'r--')
loglog(f, carm[:,31]*1e6, 'g--')
loglog(f, prcl[:,31]*1e6, 'b--')
loglog(f, mich[:,31]*1e6, 'c--')
loglog(f, srcl[:,31]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('POBS 2M1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_pobs_2m1.eps')
savefig('plots/advirgo1_dd_hp_pobs_2m1.png')

# POS 3P1
figure()
loglog(f, darm[:,33]*1e6, 'r-')
loglog(f, carm[:,33]*1e6, 'g-')
loglog(f, prcl[:,33]*1e6, 'b-')
loglog(f, mich[:,33]*1e6, 'c-')
loglog(f, srcl[:,33]*1e6, 'k-')
loglog(f, darm[:,35]*1e6, 'r--')
loglog(f, carm[:,35]*1e6, 'g--')
loglog(f, prcl[:,35]*1e6, 'b--')
loglog(f, mich[:,35]*1e6, 'c--')
loglog(f, srcl[:,35]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('POBS 3P1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_pobs_3p1.eps')
savefig('plots/advirgo1_dd_hp_pobs_3p1.png')

# XP SB1
figure()
loglog(f, darm[:,37]*1e6, 'r-')
loglog(f, carm[:,37]*1e6, 'g-')
loglog(f, prcl[:,37]*1e6, 'b-')
loglog(f, mich[:,37]*1e6, 'c-')
loglog(f, srcl[:,37]*1e6, 'k-')
loglog(f, darm[:,39]*1e6, 'r--')
loglog(f, carm[:,39]*1e6, 'g--')
loglog(f, prcl[:,39]*1e6, 'b--')
loglog(f, mich[:,39]*1e6, 'c--')
loglog(f, srcl[:,39]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('XP SB1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_xp_sb1.eps')
savefig('plots/advirgo1_dd_hp_xp_sb1.png')

# XP 2M1
figure()
loglog(f, darm[:,41]*1e6, 'r-')
loglog(f, carm[:,41]*1e6, 'g-')
loglog(f, prcl[:,41]*1e6, 'b-')
loglog(f, mich[:,41]*1e6, 'c-')
loglog(f, srcl[:,41]*1e6, 'k-')
loglog(f, darm[:,43]*1e6, 'r--')
loglog(f, carm[:,43]*1e6, 'g--')
loglog(f, prcl[:,43]*1e6, 'b--')
loglog(f, mich[:,43]*1e6, 'c--')
loglog(f, srcl[:,43]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('XP 2M1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_xp_2m1.eps')
savefig('plots/advirgo1_dd_hp_xp_2m1.png')

# XP 3P1
figure()
loglog(f, darm[:,45]*1e6, 'r-')
loglog(f, carm[:,45]*1e6, 'g-')
loglog(f, prcl[:,45]*1e6, 'b-')
loglog(f, mich[:,45]*1e6, 'c-')
loglog(f, srcl[:,45]*1e6, 'k-')
loglog(f, darm[:,47]*1e6, 'r--')
loglog(f, carm[:,47]*1e6, 'g--')
loglog(f, prcl[:,47]*1e6, 'b--')
loglog(f, mich[:,47]*1e6, 'c--')
loglog(f, srcl[:,47]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('XP 3P1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_xp_3p1.eps')
savefig('plots/advirgo1_dd_hp_xp_3p1.png')

# YP SB1
figure()
loglog(f, darm[:,49]*1e6, 'r-')
loglog(f, carm[:,49]*1e6, 'g-')
loglog(f, prcl[:,49]*1e6, 'b-')
loglog(f, mich[:,49]*1e6, 'c-')
loglog(f, srcl[:,49]*1e6, 'k-')
loglog(f, darm[:,51]*1e6, 'r--')
loglog(f, carm[:,51]*1e6, 'g--')
loglog(f, prcl[:,51]*1e6, 'b--')
loglog(f, mich[:,51]*1e6, 'c--')
loglog(f, srcl[:,51]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('YP SB1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_yp_sb1.eps')
savefig('plots/advirgo1_dd_hp_yp_sb1.png')

# YP 2M1
figure()
loglog(f, darm[:,53]*1e6, 'r-')
loglog(f, carm[:,53]*1e6, 'g-')
loglog(f, prcl[:,53]*1e6, 'b-')
loglog(f, mich[:,53]*1e6, 'c-')
loglog(f, srcl[:,53]*1e6, 'k-')
loglog(f, darm[:,55]*1e6, 'r--')
loglog(f, carm[:,55]*1e6, 'g--')
loglog(f, prcl[:,55]*1e6, 'b--')
loglog(f, mich[:,55]*1e6, 'c--')
loglog(f, srcl[:,55]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('YP 2M1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_yp_2m1.eps')
savefig('plots/advirgo1_dd_hp_yp_2m1.png')

# YP 3P1
figure()
loglog(f, darm[:,57]*1e6, 'r-')
loglog(f, carm[:,57]*1e6, 'g-')
loglog(f, prcl[:,57]*1e6, 'b-')
loglog(f, mich[:,57]*1e6, 'c-')
loglog(f, srcl[:,57]*1e6, 'k-')
loglog(f, darm[:,59]*1e6, 'r--')
loglog(f, carm[:,59]*1e6, 'g--')
loglog(f, prcl[:,59]*1e6, 'b--')
loglog(f, mich[:,59]*1e6, 'c--')
loglog(f, srcl[:,59]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('YP 3P1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_yp_3p1.eps')
savefig('plots/advirgo1_dd_hp_yp_3p1.png')


# POX SB1
figure()
loglog(f, darm[:,61]*1e6, 'r-')
loglog(f, carm[:,61]*1e6, 'g-')
loglog(f, prcl[:,61]*1e6, 'b-')
loglog(f, mich[:,61]*1e6, 'c-')
loglog(f, srcl[:,61]*1e6, 'k-')
loglog(f, darm[:,63]*1e6, 'r--')
loglog(f, carm[:,63]*1e6, 'g--')
loglog(f, prcl[:,63]*1e6, 'b--')
loglog(f, mich[:,63]*1e6, 'c--')
loglog(f, srcl[:,63]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('POX SB1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_pox_sb1.eps')
savefig('plots/advirgo1_dd_hp_pox_sb1.png')

# POX 2M1
figure()
loglog(f, darm[:,65]*1e6, 'r-')
loglog(f, carm[:,65]*1e6, 'g-')
loglog(f, prcl[:,65]*1e6, 'b-')
loglog(f, mich[:,65]*1e6, 'c-')
loglog(f, srcl[:,65]*1e6, 'k-')
loglog(f, darm[:,67]*1e6, 'r--')
loglog(f, carm[:,67]*1e6, 'g--')
loglog(f, prcl[:,67]*1e6, 'b--')
loglog(f, mich[:,67]*1e6, 'c--')
loglog(f, srcl[:,67]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('POX 2M1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_pox_2m1.eps')
savefig('plots/advirgo1_dd_hp_pox_2m1.png')

# POX 3P1
figure()
loglog(f, darm[:,69]*1e6, 'r-')
loglog(f, carm[:,69]*1e6, 'g-')
loglog(f, prcl[:,69]*1e6, 'b-')
loglog(f, mich[:,69]*1e6, 'c-')
loglog(f, srcl[:,69]*1e6, 'k-')
loglog(f, darm[:,71]*1e6, 'r--')
loglog(f, carm[:,71]*1e6, 'g--')
loglog(f, prcl[:,71]*1e6, 'b--')
loglog(f, mich[:,71]*1e6, 'c--')
loglog(f, srcl[:,71]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('POX 3P1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_pox_3p1.eps')
savefig('plots/advirgo1_dd_hp_pox_3p1.png')

# POY SB1
figure()
loglog(f, darm[:,73]*1e6, 'r-')
loglog(f, carm[:,73]*1e6, 'g-')
loglog(f, prcl[:,73]*1e6, 'b-')
loglog(f, mich[:,73]*1e6, 'c-')
loglog(f, srcl[:,73]*1e6, 'k-')
loglog(f, darm[:,75]*1e6, 'r--')
loglog(f, carm[:,75]*1e6, 'g--')
loglog(f, prcl[:,75]*1e6, 'b--')
loglog(f, mich[:,75]*1e6, 'c--')
loglog(f, srcl[:,75]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('POY SB1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_poy_sb1.eps')
savefig('plots/advirgo1_dd_hp_poy_sb1.png')

# POY 2M1
figure()
loglog(f, darm[:,77]*1e6, 'r-')
loglog(f, carm[:,77]*1e6, 'g-')
loglog(f, prcl[:,77]*1e6, 'b-')
loglog(f, mich[:,77]*1e6, 'c-')
loglog(f, srcl[:,77]*1e6, 'k-')
loglog(f, darm[:,79]*1e6, 'r--')
loglog(f, carm[:,79]*1e6, 'g--')
loglog(f, prcl[:,79]*1e6, 'b--')
loglog(f, mich[:,79]*1e6, 'c--')
loglog(f, srcl[:,79]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('POY 2M1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_poy_2m1.eps')
savefig('plots/advirgo1_dd_hp_poy_2m1.png')

# POY 3P1
figure()
loglog(f, darm[:,81]*1e6, 'r-')
loglog(f, carm[:,81]*1e6, 'g-')
loglog(f, prcl[:,81]*1e6, 'b-')
loglog(f, mich[:,81]*1e6, 'c-')
loglog(f, srcl[:,81]*1e6, 'k-')
loglog(f, darm[:,83]*1e6, 'r--')
loglog(f, carm[:,83]*1e6, 'g--')
loglog(f, prcl[:,83]*1e6, 'b--')
loglog(f, mich[:,83]*1e6, 'c--')
loglog(f, srcl[:,83]*1e6, 'k--')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('POY 3P1 phase and quadrature')
savefig('plots/advirgo1_dd_hp_poy_3p1.eps')
savefig('plots/advirgo1_dd_hp_poy_3p1.png')




# ASY DC signal
figure()
loglog(f, darm[:,85]*1e6, 'r-')
loglog(f, carm[:,85]*1e6, 'g-')
loglog(f, prcl[:,85]*1e6, 'b-')
loglog(f, mich[:,85]*1e6, 'c-')
loglog(f, srcl[:,85]*1e6, 'k-')
legend(('DARM','CARM','PRCL','MICH','SRCL'), 'lower left')
grid(True)
xlabel('Frequency [Hz]')
ylabel('Transfer function [W/um]')
title('ASY_DC signal')
savefig('plots/advirgo1_dd_hp_asy_dc.eps')
savefig('plots/advirgo1_dd_hp_asy_dc.png')