-
Notifications
You must be signed in to change notification settings - Fork 2
/
LYR_input_parameters.py
122 lines (88 loc) · 3.08 KB
/
LYR_input_parameters.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
############################
### LYR input parameters ###
############################
import numpy as np
print('\n')
print('Reading catalogues and parameters..')
# f(r) searching radius:
r_fr = 5 # in arcsec
# total(m) searching radius:
r_in_tm = 2 # in arcsec
# LR searching radius for FR
r_lr = 5 #in arcsec
# n(m) annulus searching annulus:
r_min_nm = 5 # in arcsec
r_max_nm = 30
# distributions bins:
distrib_bins = 21
# results:
path_LR = '/path/to/the/code/'
path_output = path_LR + 'output/'
path_images = path_LR + 'images/'
# output catalogue errors moltiplicative factor
plus_erf = 1 # wavdetect errors factor, put 1 if there isn't
# sigma finder ('all' == every single source; '1sigma' == 1 sigma mean errror):
sigma_out_finder = 'all'
# input catalog type(0 = static pos. error (put 0 below is there are no errors); 1 = usual errors; 2 = two catalogue with no errors):
input_cat_type = 2
# catalogues parameters:
noxy = -99. # coords of undetected sources
nomag = -99. # mag/fluxes of undetected sources
# use of fluxes of IO catalogue:
delta_f2 = False
delta_f2_fluxerr = 0 # 0: static error; 1: every single error for the input cat
# python interactions
iwp = True
cursors = False
save_images = True
save_output = True
# sdding string to all output images
add_str = 'Final_XJ'
# adding string to images titles
add_title = add_str
# out catalog and info
filename_LR = 'LR_' + add_str + '.txt'
filename_outinfo = 'LR_outinfo_' + add_str + '.txt'
############################################################### import data ########
# coordinates in degrees
######################################################################################################## ONIR - X:
# Input catalogue:
file_input = "/path/to/the/catalog/catalog.txt"
data_input = np.genfromtxt(file_input)
ID_input = np.array(data_input[:,0])
ra_input = np.array(data_input[:,1])
dec_input = np.array(data_input[:,2])
mag_input = np.array(data_input[:,7])
# Output catalogue:
file_output = "/path/to/the/catalog/catalog.txt"
data_output = np.genfromtxt(file_output)
ID_output = np.array(data_output[:,0])
ra_output = np.array(data_output[:,1])
dec_output = np.array(data_output[:,2])
ra_err_output = np.array(data_output[:,4])
dec_err_output = np.array(data_output[:,5])
#---------------------------------------------------------------------- for sigma selection ------
# Sigma input (ONIR):
file1 = "/path/to/the/catalog/catalog.txt"
data1 = np.genfromtxt(file1)
ID1_input_s = np.array(data1[:,0])
ra1_input_s = np.array(data1[:,1])
dec1_input_s = np.array(data1[:,2])
mag1_input_s = np.array(data1[:,3])
if input_cat_type == 0:
sigma_input = 0 # 0.3 per cosmos
elif input_cat_type == 1:
ra1_input_s_err = np.array(data1[:,30])
dec1_input_s_err = np.array(data1[:,30])
elif input_cat_type == 2:
mag2_input_s = np.array(data1[:,30])
ID2_input_s = np.array(data1[:,34])
ra2_input_s = np.array(data1[:,35])
dec2_input_s = np.array(data1[:,36])
# Sigma output (X-ray):
ID_output_s = ID_output
ra_output_s = ra_output
dec_output_s = dec_output
ra_err_output_s = ra_err_output
dec_err_output_s = dec_err_output
print('... done.\n')