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
|
#!/usr/bin/env python3
#
# Copyright (c) 2016 Aaron LI
# MIT license
#
# Make the configuration file for the `sbp_fit.py`.
# The name, obsid, r500 information is extracted from the '*_INFO.json' file.
#
# Created: 2016-04-21
#
# Change logs:
# 2016-07-12:
# * Rewrite according to `make_tprofile_config.py`,
# and make the sample config built-in
# 2016-04-26:
# * Minor update to output file write
#
import glob
import argparse
from info import get_name, get_obsid, get_r500
sample_config = """
## Configuration for `fit_sbp.py`
name = %(name)s
obsid = %(obsid)d
r500_pix = %(r500_pix)f
r500_kpc = %(r500_kpc)f
# sbp model: "sbeta" or "dbeta"
model = sbeta
[sbeta]
outfile = sbpfit_sbeta.json
imgfile = sbpfit_sbeta.png
#ignore = 0.0-20.0,
ignore_r500 = 0.0-0.1,
[[params]]
# name = initial, lower, upper, variable (FIXED/False to fix the parameter)
s0 = 1.0e-8, 0.0, 1.0e-6
rc = 30.0, 5.0, 1.0e4
beta = 0.7, 0.3, 1.1
bkg = 1.0e-10, 0.0, 1.0e-8
"""
def make_config(info):
cfg_data = {
"name": get_name(info)["name"],
"obsid": get_obsid(info),
"r500_pix": get_r500(info)["r500_pix"],
"r500_kpc": get_r500(info)["r500_kpc"],
}
cfg = sample_config % cfg_data
return cfg
def main():
parser = argparse.ArgumentParser(
description="Make the config for 'fit_sbp.py'")
parser.add_argument("-j", "--json", dest="json", required=False,
help="the *_INFO.json file " +
"(default: find ../*_INFO.json)")
parser.add_argument("outfile", nargs="?", default="sbpfit.conf",
help="filename of the output sbpfit config " +
"(default: sbpfit.conf)")
args = parser.parse_args()
# default "*_INFO.json"
info_json = glob.glob("../*_INFO.json")[0]
if args.json:
info_json = args.json
config = make_config(info_json)
open(args.outfile, "w").write(config)
if __name__ == "__main__":
main()
|
