import numpy as np

f=open('PopIII_haloes_0014_type5.txt','r')
f1=open('PopIII_haloes_0014_type1.txt','r')
f3=open('PopII_haloes_0014.txt','r')

#dir2='../../Rarepeak_LW1/Redshift18/'
#f4=open(dir2+'PopIII_haloes_0014_type5.txt','r')
#f5=open(dir2+'PopIII_haloes_0014_type1.txt','r')
#f6=open(dir2+'PopII_haloes_0014.txt','r')

f2=open('Overlaped_haloes_list.txt','r')
overlayed=[]
line=f2.readline()
while line !='':
   overlayed.append(int(line))
   line=f2.readline()

fo = open('star_haloes.txt','w')
line = f.readline()
line1 = f1.readline()
line3 = f3.readline()
j=0
while line !='':
   mass = line.split()[1]
   line = f.readline()
   line1 = f1.readline()
   line3 = f3.readline()
   pn = line.split()[0] 
   pn1 = line1.split()[0]
   pn3 = line3.split()[0]
   m0 = [];m1 = [];m3 = []
   for i in xrange(int(pn)):
     line=f.readline()
     starmass = float(line.split()[1])
     if starmass > 5:
        m0.append(starmass)
   for i in xrange(int(pn1)):
     line1=f1.readline()
     starmass = float(line1.split()[1])
     if starmass > 5:
       m1.append(starmass)
   for i in xrange(int(pn3)):
     line3=f3.readline()
     starmass = float(line3.split()[1])
     m3.append(starmass)
   if len(m0) == 0: 
      pm = 0; pms =0
   else:
      pm=np.array(m0).mean()
      pms=np.array(m0).std()
   if len(m1) == 0:
      pm1 = 0; pms1 =0
   else:
      pm1=np.array(m1).mean()
      pms1=np.array(m1).std()
   if len(m3) == 0:
      pm3 = 0; pms3 =0
   else:
      pm3=np.array(m3).mean()
      pms3=np.array(m3).std()
   if j not in overlayed:
      fo.write(str(mass)+" "+str(len(m0))+" "+str(pm)+" "+str(len(m1))+" "+str(pm1)+" "+str(len(m3))+" "+str(pm3)+"\n")
   f.readline();f.readline();f1.readline();f1.readline();f3.readline();f3.readline()
   line = f.readline()
   line1 = f1.readline()
   line3 = f3.readline()
   j=j+1

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