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
|
#!/usr/bin/env python3
#
# Copyright (c) 2016 Aaron LI
# MIT license
#
# Created: 2016-07-11
#
"""
Make the configuration file for `fit_tprofile.py`.
"""
import argparse
import glob
from info import get_redshift
sample_config = """
# redshift of the object (for pixel to distance conversion)
redshift = %(redshift).6f
[model_params]
# name = initial, lower, upper, variable (FIXED/False to fix the parameter)
A = %(A_val).1f, %(A_min).1f, %(A_max).1f, %(A_vary)s
n = %(n_val).1f, %(n_min).1f, %(n_max).1f, %(n_vary)s
xi = %(xi_val).1f, %(xi_min).1f, %(xi_max).1f, %(xi_vary)s
a2 = %(a2_val).1f, %(a2_min).1f, %(a2_max).1g, %(a2_vary)s
a3 = %(a3_val).1f, %(a3_min).1f, %(a3_max).1f, %(a3_vary)s
beta = %(beta_val).1f, %(beta_min).1f, %(beta_max).1f, %(beta_vary)s
T0 = %(T0_val).1f, %(T0_min).1f, %(T0_max).1f, %(T0_vary)s
"""
def parse_wang2012_param(pfile):
"""
Parse the original configuration for `fit_wang2012_model`.
"""
params = {}
for name, value, minimum, maximum, vary in map(str.split,
open(pfile).readlines()):
par = [
name,
float(value),
float(minimum),
float(maximum),
{"T": True, "F": False}.get(vary),
]
params[name] = par
return params
def make_config(params, redshift):
cfg_data = {"redshift": redshift}
for name, value, minimum, maximum, vary in params.values():
cfg_data.update({
"%s_val" % name: value,
"%s_min" % name: minimum,
"%s_max" % name: maximum,
"%s_vary" % name: vary,
})
cfg = sample_config % cfg_data
return cfg
def main():
parser = argparse.ArgumentParser(
description="Make configuration for 'fit_tprofile.py'")
parser.add_argument("-j", "--json", dest="json", required=False,
help="the *_INFO.json file " +
"(default: find ../*_INFO.json)")
parser.add_argument("-i", "--infile", dest="infile",
default="wang2012_param.txt",
help="original wang2012 model parameter file " +
"(default: wang2012_param.txt")
parser.add_argument("-o", "--outfile", dest="outfile",
default="tprofile.conf",
help="output filename")
args = parser.parse_args()
# default "*_INFO.json"
info_json = glob.glob("../*_INFO.json")[0]
if args.json:
info_json = args.json
redshift = get_redshift(info_json)
params = parse_wang2012_param(args.infile)
config = make_config(params, redshift)
open(args.outfile, "w").write(config)
if __name__ == "__main__":
main()
|
