-
Notifications
You must be signed in to change notification settings - Fork 0
/
pup_WASP18b_light_curve_compare.cfg
124 lines (102 loc) · 3.5 KB
/
pup_WASP18b_light_curve_compare.cfg
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
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
[MODEL]
# Output file to print results:
output = WASP18b_lightcurve_comparison
# Choose time units [phase jd mjd]:
timeunits = phase
# Apply an offset (in phase units) in case the values don't suit you
timeoffset = -1.0
# Least-squares fitting settings:
leastsq = True
# Select method: lm (Levenberg-Marquardt) or trf (this one repects boundries)
optimizer = trf
# Adjust data uncerts so reduced chisq = 1 (cause we don't trust Spitzer errors)
chiscale = True
# Posterior sampler setting:
walk = snooker
nsamples = 1e5
burnin = 1000
nchains = 16
ncpu = 16
thinning = 1
# Enfore constraining the transit/eclipse within first and last frames
boundecl = True
# Number of temperature calculations
numcalc = 50000
# Disply all (True) or only some (False) plots and how many bins in light curve
plots = True
nbins = 100
# Input dataset (end point of photometry step):
input = WASP18b_irac1/badpix/gauss/aper3500714/WASP18b_irac1.p
# Setup the astrophysical and systematics models:
# Concatenate as many models as you want to include in the fit
# Run multiple fits (eg. for BIC comparison) by adding lines with diff. models
model = mandelecl linramp bliss
mandelecl bliss
eclipse bliss
# Name of file(s) containing initial parameter values:
# This one sets the parameters' starting points, boundaries, and
# define free, fixed, and shared parameters.
modelfile = fit_values_WASP18b_irac1.txt
# Normalize flux at each position (set False for posflux)
normflux = True
# Iterative sigma-rejection cutoff:
sigrej = 5.0 5.0 5.0
# Clip out (remove) data points from modeling:
preclip = 0
postclip = 0
# Intermadiate ranges (one line per bracket to exclude):
#interclip = lower1 upper1
# lower2 upper2
# ...
# List of model-variables with priors (applies to each model):
#priorvars = cosi ars
# Tuple [prior, lower sigma, upper sigma] of (asymmetric) Gaussian prior:
#priorvals = 0.0849 0.002 0.002
# 9.1027 0.067 0.060
# BLISS mapping y,x axis bin sizes:
ystep = 0.015
xstep = 0.015
# Minimum number of points per BLISS bin:
minpt = 5
# Fix BLISS map to best-fit values after burn-in
isfixipmap = False
# To compare datasets, add new sections with sequential numbered names
# E.g., here we will apply the same fitting setting to all raw light curves
# Minimum variables to define are input, model, and modelfile
# But most of the previous setting can be customized for each entry
[MODEL2]
input = WASP18b_irac1/badpix/gauss/aper3750714/WASP18b_irac1.p
model = mandelecl linramp bliss
mandelecl bliss
eclipse bliss
modelfile = fit_values_WASP18b_irac1.txt
[MODEL3]
input = WASP18b_irac1/badpix/gauss/aper4000714/WASP18b_irac1.p
model = mandelecl linramp bliss
mandelecl bliss
eclipse bliss
modelfile = fit_values_WASP18b_irac1.txt
[MODEL4]
input = WASP18b_irac1/badpix/gauss/aper4250714/WASP18b_irac1.p
model = mandelecl linramp bliss
mandelecl bliss
eclipse bliss
modelfile = fit_values_WASP18b_irac1.txt
[MODEL5]
input = WASP18b_irac1/badpix/gauss/aper4500714/WASP18b_irac1.p
model = mandelecl linramp bliss
mandelecl bliss
eclipse bliss
modelfile = fit_values_WASP18b_irac1.txt
[MODEL6]
input = WASP18b_irac1/badpix/gauss/aper4750714/WASP18b_irac1.p
model = mandelecl linramp bliss
mandelecl bliss
eclipse bliss
modelfile = fit_values_WASP18b_irac1.txt
[MODEL7]
input = WASP18b_irac1/badpix/gauss/aper5000714/WASP18b_irac1.p
model = mandelecl linramp bliss
mandelecl bliss
eclipse bliss
modelfile = fit_values_WASP18b_irac1.txt