options.multiple_input(parser)
options.single_output(parser)
options.merge(parser)
(opts, args) = parser.parse_args()
if not opts.input_files and not opts.input_file_list:
parser.error('Input filename(s) not provided.')
if not opts.output_file:
parser.error('Output filename not provided.')
if not opts.bg_expect:
parser.error('Expected background density not provided.')
# Read List of Files:
if opts.input_file_list:
errors.file_name_error(opts.input_file_list)
file_list = np.genfromtxt(opts.input_file_list, dtype="S", unpack=True)
elif opts.input_files:
file_list = opts.input_files
# Read Files and Store Elements in Clusters:
clusters = []
cluster_count = 0
for file in file_list:
errors.file_name_error(file)
print 'Reading: ', file
c_id, g_num, g_id, g_ra, g_dec, g_z = read_file(file)
cluster_list = np.unique(c_id)
for clt in cluster_list:
default=range(0, 6),
help=c_col_help)
(opts, args) = parser.parse_args()
if not opts.cluster_file:
parser.error('Cluster file name not provided.')
if not opts.mock_file:
parser.error('Mock file name not provided.')
#############
# READ DATA #
#############
#make sure files exit
errors.file_name_error(opts.mock_file)
errors.file_name_error(opts.cluster_file)
print '================================================================='
#read mock catalogue
print 'Reading file:', opts.mock_file,
mock = np.genfromtxt(opts.mock_file, dtype="S", unpack=True)[opts.m_cols]
dtypes = [('id', 'S22'), ('cen', 'i4'), ('ra', 'f8'), ('dec', 'f8'),
('z', 'f8'), ('rich', 'f8'), ('mass', 'f8'), ('r200', 'f8'),
('minra', 'f8'), ('maxra', 'f8'), ('mindec', 'f8'), ('maxdec', 'f8')]
mock = np.core.records.fromarrays(mock, dtype=dtypes)
mock.rich = np.log10(mock.rich)
print '\tComplete:', len(mock)
#read cluster catalogue
options.multiple_input(parser)
options.single_output(parser)
options.merge(parser)
(opts, args) = parser.parse_args()
if not opts.input_files and not opts.input_file_list:
parser.error('Input filename(s) not provided.')
if not opts.output_file:
parser.error('Output filename not provided.')
if not opts.bg_expect:
parser.error('Expected background density not provided.')
# Read List of Files:
if opts.input_file_list:
errors.file_name_error(opts.input_file_list)
file_list = np.genfromtxt(opts.input_file_list, dtype="S", unpack = True)
elif opts.input_files:
file_list = opts.input_files
# Read Files and Store Elements in Clusters:
clusters = []
cluster_count = 0
for file in file_list:
errors.file_name_error(file)
print 'Reading: ', file
c_id, g_num, g_id, g_ra, g_dec, g_z = read_file(file)
cluster_list = np.unique(c_id)
for clt in cluster_list:
mass_bin_size = float(arg)
elif opt in ("-i", "--min_z_bin"): #set minimum for binning in redshift
min_z_bin = float(arg)
elif opt in ("-j", "--max_z_bin"): #set maximum for binning in redshift
max_z_bin = float(arg)
elif opt in ("-k", "--z_bin_size"): #set size of bins in redshift
z_bin_size = float(arg)
else:
print 'Something is wrong!'
sys.exit(2)
#############
# READ DATA #
#############
errors.file_name_error(mock_file) #make sure files exit
errors.file_name_error(cluster_file)
print '================================================================='
print 'Reading file:', mock_file,
mock_data = np.genfromtxt(mock_file, dtype="S", unpack = True) #mock haloes
print '\tComplete'
print 'Reading file:', cluster_file,
cluster_data = np.genfromtxt(cluster_file, dtype="S", unpack = True) #clusters
print '\tComplete'
####################
# DEFINE VARIABLES #
####################
parser = optparse.OptionParser()
parser.add_option("-c", "--cluster_file", dest = "cluster_file", help = "Input cluster file name.")
parser.add_option("-m", "--mock_file", dest = "mock_file", help = "Input mock file name.")
parser.add_option("-z", "--delta_z", dest = "delta_z", type = "float", help = "Delta_z value. [Default 0.03]")
(opts, args) = parser.parse_args()
if not opts.cluster_file:
parser.error('Cluster file name not provided.')
if not opts.mock_file:
parser.error('Mock file name not provided.')
#############
# READ DATA #
#############
errors.file_name_error(opts.mock_file) #make sure files exit
errors.file_name_error(opts.cluster_file)
print '================================================================='
print 'Reading file:', opts.mock_file,
mock_data = np.genfromtxt(opts.mock_file, dtype="S", unpack = True) #mock haloes
print '\tComplete'
print 'Reading file:', opts.cluster_file,
cluster_data = np.genfromtxt(opts.cluster_file, dtype="S", unpack = True) #clusters
print '\tComplete'
####################
# DEFINE VARIABLES #
####################
