def setUp(self):
mock.set_loglevel(3)
mock.cosmology_set(0.31)
self.hod = mock.Hod()
self.n = len(self.hod)
self.c = [12.0, 0.5, 0.1, 0.05, 0.95, 0.12, 13.0, 1.0, 1.5, 0.15]
for i, c in enumerate(self.c):
self.hod[i] = c
def setUp(self):
mock.set_loglevel(3)
mock.cosmology_set(0.31)
u = [2, 2, 1, 1, 0, 0, 0, 1, 0]
boxsize = 600.0
self.remap = mock.Remap(u, boxsize)
# sky
ra_range = [30.1, 38.9]
dec_range = [-6.0, -4.1]
z_range = [0.38889, 1.21239]
self.sky = mock.Sky(ra_range, dec_range, z_range)
# slice
self.slice = mock.Slice(self.remap, self.sky)
def setUp(self):
mock.set_loglevel(3)
mock.cosmology_set(0.31)
u = [2, 2, 1, 1, 0, 0, 0, 1, 0]
boxsize = 600.0
self.remap = mock.Remap(u, boxsize)
# sky
ra_range = [30.1, 38.9]
dec_range = [-6.0, -4.1]
z_range = [0.38889, 1.21239]
self.sky = mock.Sky(ra_range, dec_range, z_range)
# slice
self.slice = mock.Slice(self.remap, self.sky)
def setUp(self):
mock.set_loglevel(3)
self.ps = mock.PowerSpectrum('../data/planck1_matterpower.dat')
def setUp(self):
self.z_max = 1.5
mock.set_loglevel(3)
mock.cosmology_set(0.31)
self.distance = mock.Distance(self.z_max)
#parser.add_argument('--ssr-quad', default='vipers_quad_ssr.txt',
# help='ssr quad file')
arg = parser.parse_args()
#
# Read parameter file
#
# print('Parameter file: %s' % arg.param)
with open(arg.param, 'r') as f:
param = json.load(f)
omega_m = param['omega_m']
omega_l = 1.0 - omega_m
mock.set_loglevel(7)
mock.cosmology.set(omega_m)
mock.distance.init(1.21)
sys.stderr.write('Setting cosmology: omega_m= %.4f\n' % omega_m)
zmin = float(arg.zmin)
zmax = float(arg.zmax)
#rmin = mock.cosmology.compute_comoving_distance(1.0/(1.0 + arg.zmin))
#rmax = mock.cosmology.compute_comoving_distance(1.0/(1.0 + arg.zmax)))
#sys.stderr.write('rmin, max= %.4f %.4f\n' % (rmin, rmax))
def read_ssr_quad(filename):
"""read ssr file vipers_quad_ssr.txt and set to ssr_quad dictionary"""
ssr_quad = {}
#
# Read parameter file
#
print('# Parameter file: %s' % arg.param)
with open(arg.param, 'r') as f:
param = json.load(f)
omega_m = param['omega_m']
print('# Setting cosmology: omega_m= %.4f' % omega_m)
print('# redshift-range %f %f' % (arg.zmin, arg.zmax))
#
# Initilise
#
mock.set_loglevel(0)
mock.cosmology.set(omega_m)
mock.distance.init(1.2)
def read_catalogues(filebase, irange):
cats = mock.Catalogues()
for i in range(int(irange[0]), int(irange[1]) + 1):
filename = '%s%05d.txt' % (filebase, i)
a = np.loadtxt(filename, delimiter=' ', usecols=[1,2,3,6,7])
cats.append(a, z_min=arg.zmin, z_max= arg.zmax)
return cats
galaxies = read_catalogues('../mocks/%s/mock_%s_' % (arg.reg, arg.reg), igalaxies)
randoms = read_catalogues('../rands/%s/rand_%s_' % (arg.reg, arg.reg), irandoms)
with open(arg.param, 'r') as f:
param = json.load(f)
# find region name w1 or w4
for reg in param['reg']:
if reg['name'] == arg.reg:
break
if reg['name'] != arg.reg:
raise Exception('region %s not found in parameter file %s' %
(arg.reg, arg.param))
print('Region: %s' % reg['name'])
mock.set_loglevel(0)
# omega_m
omega_m = param['omega_m']
print('Setting cosmology: omega_m= %.4f' % omega_m)
mock.cosmology.set(omega_m)
# power spectrum
mock.power.init(param['power_spectrum'])
mock.sigma.init(M_min=9.0e9)
# Random
if arg.random:
random = True
out_dir = 'rand_lightcone'
print('random')
#parser.add_argument('--ssr-quad', default='vipers_quad_ssr.txt',
# help='ssr quad file')
arg = parser.parse_args()
#
# Read parameter file
#
# print('Parameter file: %s' % arg.param)
with open(arg.param, 'r') as f:
param = json.load(f)
omega_m = param['omega_m']
omega_l = 1.0 - omega_m
mock.set_loglevel(7)
mock.cosmology.set(omega_m)
mock.distance.init(1.21)
sys.stderr.write('Setting cosmology: omega_m= %.4f\n' % omega_m)
zmin = float(arg.zmin)
zmax = float(arg.zmax)
#rmin = mock.cosmology.compute_comoving_distance(1.0/(1.0 + arg.zmin))
#rmax = mock.cosmology.compute_comoving_distance(1.0/(1.0 + arg.zmax)))
#sys.stderr.write('rmin, max= %.4f %.4f\n' % (rmin, rmax))
def read_ssr_quad(filename):
"""read ssr file vipers_quad_ssr.txt and set to ssr_quad dictionary"""
ssr_quad = {}
for line in open(filename, 'r'):
def setUp(self):
mock.set_loglevel(3)
mock.cosmology_set(0.31)
def setUp(self):
mock.set_loglevel(3)
mock.cosmology_set(0.31)
self.ps = mock.PowerSpectrum('../data/planck1_matterpower.dat')
self.s = mock.Sigma(self.ps)
self.mf = mock.MassFunction(self.s)
def setUp(self):
mock.set_loglevel(3)
mock.cosmology_set(0.31)
self.ps = mock.PowerSpectrum('../data/planck1_matterpower.dat')
self.s = mock.Sigma(self.ps)
self.mf = mock.MassFunction(self.s)
def setUp(self):
mock.set_loglevel(3)
self.ps= mock.PowerSpectrum('../data/planck1_matterpower.dat')
def setUp(self):
mock.set_loglevel(3)
self.lightcones = mock.LightCones()
def setUp(self):
mock.set_loglevel(3)
self.lightcones = mock.LightCones()
