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
