1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
|
import tables as tb
class Particle(tb.IsDescription):
name = tb.StringCol(16) # 16-character String
idnumber = tb.Int64Col() # Signed 64-bit integer
ADCcount = tb.UInt16Col() # Unsigned short integer
TDCcount = tb.UInt8Col() # Unsigned byte
grid_i = tb.Int32Col() # Integer
grid_j = tb.IntCol() # Integer (equivalent to Int32Col)
pressure = tb.Float32Col() # Float (single-precision)
energy = tb.FloatCol() # Double (double-precision)
with tb.open_file("tutorial.h5", mode="w", title="Test file") as h5file:
group = h5file.create_group("/", "detector", "Detector information")
table = h5file.create_table(group, "readout", Particle, "Readout example")
print(h5file)
particle = table.row
for i in range(10):
particle["name"] = f"Particle: {i:6d}"
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)
particle.append()
table.flush()
with tb.open_file("tutorial.h5", mode="r", title="Test file") as h5file:
table = h5file.root.detector.readout
pressure = [
x["pressure"]
for x in table.iterrows()
if x["TDCcount"] > 3 and 20 <= x["pressure"] < 50
]
print(pressure)
|
