import numpy as np
import glob
import re

filelist1 = glob.glob("XrayLuminosity_512*")
filelist2 = glob.glob("XrayLuminosity_WV*")
filelist3 = glob.glob("XrayLuminosity_??.??.dat")
filelist4 = glob.glob("XrayLuminosity_??.??_density*")
filelist1.sort(reverse=True)
filelist2.sort(reverse=True)
filelist3.sort(reverse=True)
filelist4.sort(reverse=True)

fo = open("X-ray_luminosity.txt",'w')
fo.write("%-16s%-16s%-16s%-16s%-16s\n"%("redshift","refined volume","whole volume",'fitting','global density'))

for name1,name2,name3,name4 in zip(filelist1,filelist2,filelist3,filelist4):

  redshift=float(re.findall(r"[-+]?\d*\.\d+|\d+",name1)[1])
  Xray1=np.fromfile(name1)
  Xray2=np.fromfile(name2)
  Xray3=np.fromfile(name3)
  Xray4=np.fromfile(name4)
  fo.write("%-16.2f%-16e%-16e%-16e%-16e\n"%(redshift,Xray1.sum(),Xray2.sum(),Xray3.sum(),Xray4.sum()))

fo.close()