# Check the distribution j = np.zeros(6) jj =[] for jjj inrange(len(j)): jj.append(np.round(pow(0.1, jjj+1),8)) for i inrange(len(f_nan_0)): for ii inrange(len(j)): if f_nan_0[i]< jj[ii]: j[ii]+=1 print('Energy cut is:') print(' '.join(map(str, jj))) print('Labels below energy cut:') print(' '.join(map(str, j))) print('Percentage below energy cut:') print(' '.join(map(str, np.round(j/len(f_nan_0),4))))