/
summarize-timestamps
32 lines (28 loc) · 1.05 KB
/
summarize-timestamps
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
#!/usr/bin/env python3
# -*- coding: utf-8 -*-
import argparse
import numpy as np
from photon_tools.io import read_photons
import matplotlib.pyplot as pl
from photon_tools.bin_photons import bin_photons
from scipy.stats import linregress
parser = argparse.ArgumentParser()
parser.add_argument('file', nargs='+', type=argparse.FileType('r'), help='Timestamp file')
args = parser.parse_args()
for f in args.file:
ts = read_photons.open(f.name)
print('%s:' % (f.name))
print(' jiffy = %g' % ts.jiffy)
for c in ts.valid_channels:
d = ts.channel(c)
N = len(d)
if N == 0:
continue
dur = d[-1] - d[0]
dur *= ts.jiffy
bins = bin_photons(d, 1. / ts.jiffy)
m, b, r, p, stderr = linregress(bins['start_t'], bins['count'])
print(' channel %s: T=%1.1f seconds, N=%1.2e' % (c, dur, N))
print(' rate = %1.1f Hz, variance = %1.1f Hz²' % (N / dur, np.var(bins['count'])))
print(' dI/dt = %g ± %1.1g Hz / second' % (m / ts.jiffy, stderr / ts.jiffy))
print()