power_per_band = [np.sum(psd[np.bitwise_and(freqs >= low, freqs<up)])