hist.py 284 B

123456789101112131415
  1. import sys
  2. import numpy as np
  3. data = np.loadtxt(sys.argv[1])
  4. bias = 7.8
  5. num_bins = 75
  6. data = data[data < bias]
  7. hist, bins = np.histogram(data, bins=num_bins)
  8. hist = hist.reshape((-1, 1))
  9. bins = bins.reshape((-1, 1))
  10. out = np.hstack((bins[:-1], hist))
  11. np.savetxt(sys.argv[2], out)