import numpy as np
import glob
import re

filelist1 = glob.glob("XrayLuminosity_512_??.??.dat")
filelist2 = glob.glob("XrayLuminosity_WV_??.??.dat")
filelist5 = glob.glob("XrayLuminosity_??.??.dat")
filelist7 = glob.glob("XrayLuminosity_??.??_density.dat")
filelist1.sort(reverse=True)
filelist2.sort(reverse=True)
filelist5.sort(reverse=True)
filelist7.sort(reverse=True)
filelist3 = glob.glob("XrayLuminosity_512_?.??.dat")
filelist4 = glob.glob("XrayLuminosity_WV_?.??.dat")
filelist6 = glob.glob("XrayLuminosity_?.??.dat")
filelist8 = glob.glob("XrayLuminosity_?.??_density.dat")
filelist3.sort(reverse=True)
filelist4.sort(reverse=True)
filelist6.sort(reverse=True)
filelist8.sort(reverse=True)

filelist1+=filelist3
filelist2+=filelist4
filelist5+=filelist6
filelist7+=filelist8

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,filelist5,filelist7):

  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()