import os
from dessn.models.d_simple_stan.load import plot_quick
if __name__ == "__main__":
dir_name = os.path.dirname(__file__)
std = dir_name + "/stan_output"
plot_quick(std, "stan_mc", include_sep=False)
c = ChainConsumer()
if do_weight:
c.add_chain(chain, weights=w, name="calib")
else:
c.add_chain(chain, name="calib")
if load_second:
res2 = load_stan_from_folder(std + "_calib_data_no_calib_model_and_bias", num=0, merge=False, cut=False)
chain, posterior, _, p, f, l, w, ow = res2[0]
if do_weight:
c.add_chain(chain, weights=w, name="nocalib")
else:
c.add_chain(chain, name="nocalib")
c.plot(filename="output.png", truth=t, figsize=0.75)
# c = ChainConsumer()
# c.add_chain(chain, name="approx")
# c.add_chain(chain, weights=w, name="full")
# c.plot(filename="output.png", truth=t)
if __name__ == "__main__":
dir_name = os.path.dirname(__file__)
std = dir_name + "/stan_output"
plot_quick(std, "approx", include_sep=False)
# plot_all_weight(std, dir_name + "/plot_approx_weight.png")
#for i in range(20):
# plot_single_cosmology_weight(std, dir_name + "/plot_approx_single_weight_%d.png" % i, i=i)
# debug_plots(std)