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)