# 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