from sedbot.probf import LnUniform, LnNormal
from sedbot.ensemble_multipix.gibbsbg import MultiPixelGibbsBgModeller
from sedbot.ensemble_multipix.threeparam_lnprob import ThreeParamLnProb, \
GlobalThreeParamLnProb
N_PIXELS = 5
BANDS = ['sdss_u', 'sdss_g', 'sdss_r', 'sdss_i', '2mass_J', '2mass_Ks']
D0 = 785. * 1000. # distance in parsecs (McConnachie 2005)
D0_SIGMA = 25. * 1000.
PIX_AREA = 1. # arcsec-sq
MSUN_I = fsps.get_filter('sdss_i').msun_ab
BKG_AMPLITUDE = 1e-10 # units of janksy per sq arcsec
# Hard limits on parameters
LIMITS = {'logmass': (sb_to_mass(30., MSUN_I, PIX_AREA, 0.4, D0),
sb_to_mass(12., MSUN_I, PIX_AREA, 0.4, D0)),
'logtau': (-1., 2.),
'const': (0., 1.0),
'sf_start': (0.5, 10.),
'tburst': (10.0, 13.8),
'fburst': (0., 1.0),
'logZZsol': (-1.98, 0.2),
'd': (D0 - 3. * D0_SIGMA, D0 + 3. * D0_SIGMA),
'dust1': (0., 5.),
'dust2': (0., 3.),
'ml': (-0.5, 1.5)}
# Prior probability functions
PRIOR_FUNCS = {"logtau": LnUniform(*LIMITS['logtau']),
def __init__(self, logml_min, logml_max, sb, D_pc, area, msun):
self._low_mass = sb_to_mass(sb, msun, logml_min, area, D_pc)
self._high_mass = sb_to_mass(sb, msun, logml_max, area, D_pc)
super(LnUniformMass, self).__init__(self._low_mass, self._high_mass)