import numpy as np

f=open('RedshiftOutput0028-halos.dat','r')
lines=f.readlines()
lines.pop(0)

sm = []
for line in lines:
    sm.append(float(line.split()[6]))

sm=np.array(sm)
sm=sm[sm>0]
a,b=np.histogram(np.log10(sm),bins=14,range=[2,9])

f=open("stellar_mass_function.txt",'w')
for i,c in enumerate(a):
    f.write(str((b[i]+b[i+1])/2.0)+" "+str(c/0.5/(6*6*6.125))+"\n")
f.close()