-
Notifications
You must be signed in to change notification settings - Fork 1
/
bestfit.py
129 lines (127 loc) · 4.85 KB
/
bestfit.py
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
124
125
126
127
128
# Best-fit params from optical(g,z) + radio, no SMF, 2016-10-19
params_bf_1 = {'sfr_sfms_sigma': 0.37189501386500001,
'mass_mstar_min': 6.0,
'opt_bands': ['u', 'g', 'r', 'i', 'z'],
'fpass_alpha0': 10.7816242649,
'fpass_alpha1': 0.5,
'sfr_pass_pthres': 0.05,
'fpass_zeta': -1.5124955364999999,
'ms_cen_nu': -0.72,
'mass_mstar_max': 14.,
'opt_pdf_sigma': [0.28122634653850337, 0.25232918833346546, 0.2468073941298409, 0.25273681573440887, 0.2724513351999828],
'extinction_alpha': -7.57886,
'extinction_beta': 5.4443420979599999,
'sfr_pass_beta': 0.37,
'extinction_diskfac': 4.9942409330900004,
'mass_mhalo_max': 16.0,
'ms_cen_sigma1': 0.57533203575699998,
'opt_offset': [-27.3176081, -25.9445857, -25.2836569, -24.981922, -24.8096689],
'opt_cross_amp': [-3.61194266, -2.70027825, -2.01635745, -1.72575487, -1.56393268],
'opt_mstar_beta': -0.000242854097,
'mhi_omegaHI': 0.0001,
'ms_cen_mu': 0.019,
'nsamp_mstar': 300,
'ms_cen_norm': 0.017996514026499998,
'opt_cross_beta': -0.262636097,
'ms_cen_gamma0': 0.558604689313,
'ms_cen_gamma1': -0.26,
'extinction_amp': 0.615135404291,
'mass_mhalo_min': 7.0,
'ms_cen_xi': 4.25,
'opt_cross_gamma': 0.290177366,
'opt_mstar_amp': 3876.55879,
'ms_cen_logM2': 11.8,
'ms_cen_logM1': 12.4992783145,
'sfr_pass_gamma': 1.9,
'sfr_min': -9.0,
'sfr_pass_mshift': 0.0024012941371300001,
'ms_cen_sigmainf': 0.022, #-0.031528396942999998,
'sfr_pass_sfrmin': 1e-07,
'sfr_sfms_gamma': 1.9,
'sfr_pass_mscale': 1.0,
'sfr_sfms_mscale': 1.0,
'sfr_sfms_beta': 1.09014789144,
'nsamp_mhalo': 250,
'sfr_sfms_alpha1': 1.07,
'sfr_sfms_alpha0': -0.060352680060800003,
'opt_pdf_mean': [-0.06587481168919591, -0.053777765775930214, -0.01854712885192855, -0.008538656095465969, -0.008732300503716532],
'nsamp_sfr': 300,
'sfr_pass_sigma': 0.080482132361400005,
'ms_cen_beta0': 1.2865496082700001,
'ms_cen_beta1': 0.17,
'opt_mstar_c': -1.0108709,
'sfr_max': 4.0,
'fpass_beta': -2.5515276576699999,
'sfr_pass_type': 'shifted sfms',
'sfr_pass_alpha1': 0.0,
'sfr_pass_alpha0': -5.0,
'mhi_vc1': 200.0,
'mhi_vc0': 50.0,
'extinction_tau0': 4.52232780e-01,
'extinction_kappa': 3.29006914e-03,
'extinction_lambda0': 4.96662013e+03,
}
# Best-fit params from optical(g,z) + radio, no SMF, better priors, 2016-10-26
params_bf = {
'sfr_sfms_sigma': 0.39059251672099998,
'mass_mstar_min': 6.0,
'opt_bands': ['u', 'g', 'r', 'i', 'z'],
'fpass_alpha0': 10.804369278199999,
'fpass_alpha1': 0.5,
'sfr_pass_pthres': 0.05,
'fpass_zeta': -1.62086498436,
'ms_cen_nu': -0.72,
'mass_mstar_max': 14.0,
'opt_pdf_sigma': [0.28122634653850337, 0.25232918833346546, 0.2468073941298409, 0.25273681573440887, 0.2724513351999828],
'extinction_kappa': 0.003534, #0.00329006914,
'extinction_alpha': -7.57886,
'extinction_beta': 5.4783395970699997,
'extinction_lambda0': 4977.189036, #4966.62013,
'sfr_pass_beta': 0.37,
'extinction_diskfac': 1.4294834244600001,
'mass_mhalo_max': 16.0,
'ms_cen_sigma1': 0.55697362786500004,
'opt_offset': [-27.3176081, -25.9445857, -25.2836569, -24.981922, -24.8096689],
'ms_cen_logM1': 12.498161809000001,
'mhi_omegaHI': 0.0001,
'ms_cen_mu': 0.019,
'nsamp_mstar': 300,
'ms_cen_norm': 0.017739115493300001,
'nsamp_mhalo': 250,
'opt_cross_beta': -0.262636097,
'ms_cen_gamma0': 0.59559404090000001,
'ms_cen_gamma1': -0.26,
'extinction_amp': 0.615135404291,
'mass_mhalo_min': 7.0,
'ms_cen_xi': 4.25,
'opt_cross_gamma': 0.290177366,
'opt_mstar_amp': 3876.55879,
'ms_cen_logM2': 11.8,
'opt_mstar_beta': -0.000242854097,
'sfr_pass_type': 'shifted sfms',
'sfr_min': -9.0,
'sfr_pass_mshift': 0.00105527163519,
'ms_cen_sigmainf': 0.030873590974500001,
'sfr_pass_alpha1': 0.0,
'sfr_sfms_gamma': 1.9,
'nsamp_sfr': 300,
'sfr_sfms_mscale': 1.0,
'sfr_sfms_beta': 1.0373848138599999,
'opt_cross_amp': [-3.61194266, -2.70027825, -2.01635745, -1.72575487, -1.56393268],
'sfr_sfms_alpha1': 1.07,
'sfr_sfms_alpha0': -0.077393675287999994,
'opt_pdf_mean': [-0.06587481168919591, -0.053777765775930214, -0.01854712885192855, -0.008538656095465969, -0.008732300503716532],
'sfr_pass_mscale': 1.0,
'extinction_tau0': 1.040598, #0.74590281490999999,
'sfr_pass_sigma': 0.028559680849299999,
'ms_cen_beta0': 1.3211623293100001,
'ms_cen_beta1': 0.17,
'opt_mstar_c': -1.0108709,
'sfr_max': 4.0,
'fpass_beta': -2.4362458372,
'sfr_pass_gamma': 1.9,
'sfr_pass_sfrmin': 1e-07,
'sfr_pass_alpha0': -5.0,
'mhi_vc1': 200.0,
'mhi_vc0': 50.0
}