#
p.xlabel(qty) #'Observed fraction per '+pixel_area_str+r' deg$^2$ pixels')
p.ylabel('fraction')
p.grid()
p.ylim((0, 1))
#p.yscale('log')
p.legend(frameon=False, loc=4)
p.savefig(fig_out)
p.clf()
# NSIDE=2^5=32
#
nside_int = int(sys.argv[1])
nside = 2**nside_int
nside_str = 'hp' + str(nside_int)
pixel_area = healpy.nside2pixarea(nside, degrees=True)
pixel_area_str = str(np.round(pixel_area, 2))
uniq_all = np.unique(data[nside_str], return_counts=True)
uniq_observed = np.unique(data[nside_str][data['in_BOSS_v5_13_0']],
return_counts=True)
uniq_goodZ = np.unique(data[nside_str][(data['in_BOSS_v5_13_0'])
& (data['ZWARNING'] == 0)],
return_counts=True)
N_pixels_with_targets = len(uniq_all[0])
N_pixels_with_observed_targets = len(uniq_observed[0])
print('N_pixels_with_targets', N_pixels_with_targets)
print('N_pixels_with_observed_targets', N_pixels_with_observed_targets)
N_pixels_with_goodZ = len(uniq_goodZ[0])
print('N_pixels_with_goodZ', N_pixels_with_goodZ)
from astropy.table import Table, Column
import astropy.units as u
import numpy as n
import extinction
print('CREATES 4FS FITS FILES')
print('------------------------------------------------')
print('------------------------------------------------')
# astropy
#from lib_model_agn import *
#option = sys.argv[1]
#option = 'SNR3'
#option = 'eRASS8'
#option = 'eRASS3'
area_per_cat = healpy.nside2pixarea(8, degrees=True)
env = sys.argv[1]
# simulation setup
if env == "MD10" or env == "MD04":
from astropy.cosmology import FlatLambdaCDM
import astropy.units as u
cosmoMD = FlatLambdaCDM(H0=67.77 * u.km / u.s / u.Mpc,
Om0=0.307115) # , Ob0=0.048206)
h = 0.6777
L_box = 1000.0 / h
cosmo = cosmoMD
if env == "UNIT_fA1_DIR" or env == "UNIT_fA1i_DIR" or env == "UNIT_fA2_DIR":
from astropy.cosmology import FlatLambdaCDM
import astropy.units as u
from scipy.interpolate import interp1d
import matplotlib
matplotlib.use('Agg')
matplotlib.rcParams.update({'font.size': 14})
import matplotlib.pyplot as p
import astropy.io.fits as fits
#import h5py
import numpy as n
print('PLOT XLF and logNlogS')
print('------------------------------------------------')
print('------------------------------------------------')
t0 = time.time()
#import astropy.io.fits as fits
area = healpy.nside2pixarea(8, degrees=True)
matplotlib.use('Agg')
matplotlib.rcParams.update({'font.size': 14})
# import all pathes
env = sys.argv[1] # 'MD04'
ftyp = sys.argv[2] # "sat"
print(env, ftyp)
if env == "MD10" or env == "UNIT_fA1_DIR" or env == "UNIT_fA1i_DIR" or env == "UNIT_fA2_DIR":
z_maximum = 6.0
if env == "MD04":
z_maximum = 0.45
DZ = 0.2