# sky sampler: cos(dec) uniform in [-1, 1), adaptive sampling
#
dec_sampler = mcsampler.dec_samp_vector
sampler.add_parameter("declination", \
pdf = dec_sampler, \
cdf_inv = None, \
left_limit = param_limits["declination"][0], \
right_limit = param_limits["declination"][1], \
prior_pdf = mcsampler.uniform_samp_dec, \
adaptive_sampling = not opts.no_adapt)
#
# Determine pinned and non-pinned parameters
#
pinned_params = common_cl.get_pinned_params(opts)
unpinned_params = common_cl.get_unpinned_params(opts, sampler.params)
print "{0:<25s} {1:>5s} {2:>5s} {3:>20s} {4:<10s}".format("parameter", "lower limit", "upper limit", "pinned?", "pin value")
plen = len(sorted(sampler.params, key=lambda p: len(p))[-1])
for p in sampler.params:
if pinned_params.has_key(p):
pinned, value = True, "%1.3g" % pinned_params[p]
else:
pinned, value = False, ""
if isinstance(p, tuple):
for subp, subl, subr in zip(p, sampler.llim[p], sampler.rlim[p]):
subp = subp + " "*min(0, plen-len(subp))
print "|{0:<25s} {1:>1.3g} {2:>1.3g} {3:>20s} {4:<10s}".format(subp, subl, subr, str(False), "")
else:
p = p + " "*min(0, plen-len(p))
# sky sampler: cos(dec) uniform in [-1, 1), adaptive sampling
#
dec_sampler = mcsampler.dec_samp_vector
sampler.add_parameter("declination", \
pdf = dec_sampler, \
cdf_inv = None, \
left_limit = param_limits["declination"][0], \
right_limit = param_limits["declination"][1], \
prior_pdf = mcsampler.uniform_samp_dec, \
adaptive_sampling = not opts.no_adapt)
#
# Determine pinned and non-pinned parameters
#
pinned_params = common_cl.get_pinned_params(opts)
unpinned_params = common_cl.get_unpinned_params(opts, sampler.params)
print "{0:<25s} {1:>5s} {2:>5s} {3:>20s} {4:<10s}".format(
"parameter", "lower limit", "upper limit", "pinned?", "pin value")
plen = len(sorted(sampler.params, key=lambda p: len(p))[-1])
for p in sampler.params:
if pinned_params.has_key(p):
pinned, value = True, "%1.3g" % pinned_params[p]
else:
pinned, value = False, ""
if isinstance(p, tuple):
for subp, subl, subr in zip(p, sampler.llim[p], sampler.rlim[p]):
subp = subp + " " * min(0, plen - len(subp))
print "|{0:<25s} {1:>1.3g} {2:>1.3g} {3:>20s} {4:<10s}".format(
subp, subl, subr, str(False), "")