amps.py

na = np.loadtxt('../nanoamps')
_,bins,_ = plt.hist(na, bins=40, color='gray', alpha=0.2, density=False);
plt.hist(na.reshape(4,-1).T, bins=bins, histtype='step', density=False);
na_ = na.reshape(4,-1)
plt.hist(np.mean(na_, axis=0), histtype='step', bins=30);
_,bins,_ = plt.hist(na, bins=40, color='gray', alpha=0.2, density=False);

plt.hist(np.std(na_, axis=0, ddof=1), bins=bins);

np.argmax(np.std(na_,axis=0,ddof=1))
np.argwhere(na_ == np.max(na_))