PyTables example

import tables
from tables import *
import tables
import numpy
class Particle(IsDescription):
     name      = StringCol(16)  
     idnumber  = Int64Col()     
     ADCcount  = UInt16Col()    
     TDCcount  = UInt8Col()     
     grid_i    = Int32Col()     
     grid_j    = Int32Col()     
     pressure  = Float32Col()   
     energy    = Float64Col()

h5file = open_file("tutorial1.h5", mode = "w", title = "Test file")
group = h5file.create_group("/", 'detector', 'Detector information')
table = h5file.create_table(group, 'readout', Particle, "Readout example")
print(h5file)
h5file
particle = table.row

for i in xrange(10):
     particle['name']  = 'Particle: %6d' % (i)
     particle['TDCcount'] = i % 256
     particle['ADCcount'] = (i * 256) % (1 << 16)
     particle['grid_i'] = i
     particle['grid_j'] = 10 - i
     particle['pressure'] = float(i*i)
     particle['energy'] = float(particle['pressure'] ** 4)
     particle['idnumber'] = i * (2 ** 34)
     # Insert a new particle record
     particle.append()

table.flush()

table = h5file.root.detector.readout
pressure = [x['pressure'] for x in table.iterrows() if x['TDCcount'] > 3 and 20 <= x['pressure'] < 50]
pressure