Python script to separate the channel data from the auxiliary data in the HDF5 files.

separate_data_aux.py

ns = 30000    # number of samples

f = h5py.File('/home/chrono/sng/data/take5.h5')
dset = f['wired-dataset'] # this has length 30000, where each element is a non-rectangular tuple

data = np.zeros((1024,ns), dtype='int16')      # 1024 channels at 30 kHz
aux = np.zeros((32,45,ns/15), dtype='int16')   # 32 chips, 45 aux values, samples at 2 kHz

for i in range(ns):
    data[:,i] = dset[i][3][:1024]
    for j in range(96):
        aux[j//3, (i%15)*3+j%3, i//15] = dset[i][3][1024+j]