import numpy as np

f=open('PopIII_list_0171_type1.txt','r')
f1=open('PopIII_list_0171_type5.txt','r')

fo = open('PopIII_histogram.txt','w')
starlist = []
line = f.readline()
line = f.readline()
while line !='':
  starmass = float(line.split()[1])
  #if starmass > 40:
  starlist.append(starmass)   
  line = f.readline()

line = f1.readline()
line = f1.readline()
while line !='':
  starmass = float(line.split()[1])
  if starmass > 1e20: starmass= starmass/1e20
  starlist.append(starmass)
  line = f1.readline()

starcount = np.zeros((300-5)/5)
for starmass in starlist:
    starcount[int((starmass-5)/5)] = starcount[int((starmass-5)/5)]+1

bin = 7.5
for starnumber in starcount:
  fo.write(str(bin)+" "+str(starnumber)+"\n")
  bin = bin + 5   


f.close()
f1.close()
fo.close()