def explore_ob150211(): data = munge.getdata('ob150211', use_astrom_phot=True) mod_base = '/u/jlu/work/microlens/OB150211/a_2019_06_26/model_fits/' mod_base += '8_fit_multiphot_astrom_parallax2/tmp/aa_' mod_fit = model_fitter.PSPL_multiphot_astrom_parallax2_Solver( data, outputfiles_basename=mod_base) #mod_fit.summarize_results_modes() tab_all = mod_fit.load_mnest_modes() tab = tab_all[0] tab['mL'] = 10**tab['log_mL'] tab['piE'] = 10**tab['log_piE'] tab['thetaE'] = 10**tab['log_thetaE'] plt.figure(1) plt.clf() plt.hist(tab['log_thetaE'], weights=tab['weights'], bins=50) plt.xlabel('log(thetaE)') plt.figure(3) plt.clf() plt.hist(tab['log_mL'], weights=tab['weights'], bins=50) plt.xlabel('log(mL)') plt.figure(4) plt.clf() plt.hist(tab['log_piE'], weights=tab['weights'], bins=50) plt.xlabel('log(piE)') plt.figure(5) plt.clf() plt.hist(tab['thetaE'], weights=tab['weights'], bins=50) plt.xlabel('thetaE') plt.figure(6) plt.clf() plt.hist(tab['mL'], weights=tab['weights'], bins=50) plt.xlabel('mL') plt.figure(7) plt.clf() plt.hist(tab['piE'], weights=tab['weights'], bins=50) plt.xlabel('piE') plt.figure(8) plt.clf() plt.plot(tab['logLike'], tab['log_mL'], 'k.') summarize_results(tab) pdb.set_trace() return
def plot_ob140613_phot_ast(): data = munge.getdata('ob140613', use_astrom_phot=True) mod_base = '/u/jlu/work/microlens/OB140613/a_2019_06_26/model_fits/' mod_base += '8_fit_multiphot_astrom_parallax2/aa_' mod_fit = model_fitter.PSPL_multiphot_astrom_parallax2_Solver(data, outputfiles_basename=mod_base) #mod_fit.separate_modes() mod_fit.summarize_results_modes() mod_all = mod_fit.get_best_fit_modes_model(def_best='maxL') tab_all = mod_fit.load_mnest_modes() plot_4panel(data, mod_all[0], tab_all[0], 'ob140613_phot_astrom.png', r_min_k=4.0, mass_max_lim=2, log=False) return
def plot_ob150211_phot_ast(): data = munge.getdata('ob150211', use_astrom_phot=True) mod_base = '/u/jlu/work/microlens/OB150211/a_2019_06_26/model_fits/' mod_base += '8_fit_multiphot_astrom_parallax2/tmp/aa_' mod_fit = model_fitter.PSPL_multiphot_astrom_parallax2_Solver(data, outputfiles_basename=mod_base) #mod_fit.separate_modes() #mod_fit.summarize_results_modes() mod_all = mod_fit.get_best_fit_modes_model(def_best='maxL') tab_all = mod_fit.load_mnest_modes() tab_all[0]['weights'] = tab_all[0]['weights'] / tab_all[0]['weights'].sum() plot_4panel(data, mod_all[0], tab_all[0], 'ob150211_phot_astrom.png', mass_max_lim=10, log=True) return