#!/usr/bin/env python

import os
#import numpy as np
from pylab import *

#name = sys.argv[1]

#atten = 10 # 10dB

#for meas in ['9_24', '9_25', '9_26', '9_27', '9_28']:
#for meas in os.listdir(sys.argv[1]):
for meas in sorted(sys.argv[1:]):
    if not meas.find('.dat') > 0:
        continue

    #path = '%s/%s' % (sys.argv[1], meas)
    path = meas
    label = '/'.join(meas.split('/')[1].split('.dat')[0].split('_')[0:2])
    print label

    data = loadtxt(path, delimiter=',')

    loglog(data[:,0], data[:,1], label=label)
    #semilogy(data[:,0], data[:,1], label=label)
    hold('on')

axis('tight')
grid('on')
legend()
title('MC error point spectrum')
xlabel('Frequency (Hz)')
ylabel('Volts')

show()
# show(False)
# input("Enter to exit:")
# sys.exit()

#savefig('MCspectrum.pdf')