def coincidences(data, max_tau=(5*11.6e6 + 1e6)):
taus = T2_tools.get_coincidence_intervals(data, max_tau)
# make a histogram with one bin per 10us (interval between fake photons)
# such that expected coincidence intervals are centered in the bin
h,e = np.histogram(taus * 1e-3, bins=np.arange(6)*11.6e3 - 11.6e3/2)
print
print 'Coincidence intervals'
print 'bin edges (ns): '
print e
print 'events: '
print h
print
