mod = StarModel(dar, obs=t)
mod.fit_multinest(n_live_points=1000,
basename='test{}_bound'.format(i))
#if rank == 0:
f1 = open('evidence_bound2.txt','a')
evi = mod.evidence
evi = str(evi)
f1.write('case{}: '.format(i) + evi + '\n')
#mod.corner(['mass_0_0','mass_0_1','distance_0','AV_0'])
fig = mod.corner_physical(props=['mass', 'distance', 'AV'])
fig.savefig('test{}_bound_corner_physical.png'.format(i))
plt.close(fig)
fig = mod.corner_observed()
fig.savefig('test{}_bound_corner_observed.png'.format(i))
plt.close(fig)
#unassociated case
dar = Dartmouth_Isochrone()
t = ObservationTree.from_df(df, name='test{}'.format(i))
t.define_models(dar, index = [0,1])
mod = StarModel(dar, obs=t)
mod.fit_multinest(n_live_points=1000,
basename='test{}_unassociated'.format(i))
#if rank == 0:
f2 = open('evidence_unassociated2.txt','a')
evi = mod.evidence
