import numpy as np

f=open('PopII_haloes_0038.txt','r')

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

fo = open('PopII_haloes.txt','w')
line = f.readline()
j=0
while line !='':
   mass = line.split()[1]
   line = f.readline()
   pn = line.split()[0] 
   m = []
   for i in xrange(int(pn)):
     line=f.readline()
     starmass = float(line.split()[1])
     m.append(starmass)
   if len(m) == 0: 
      pm = 0; pms =0
   else:
      pm=np.array(m).mean()
      pms=np.array(m).std()
   if j not in overlayed:
      fo.write(str(mass)+" "+str(len(m))+" "+str(pm)+" "+str(pms)+"\n")
   f.readline();f.readline()
   line = f.readline()
   j=j+1

f.close()
fo.close()