import pyhande
import matplotlib.pyplot as plt
ifciqmc_data = pyhande.extract.extract_data('calcs/ifciqmc/hubbard_ifciqmc.out')
for (md, data) in ifciqmc_data[:3]:
    plt.plot(data['iterations'], data['\sum H_0j N_j']/data['N_0'], label=r'$N_{\mathrm{target}} = %s$' % int(md['qmc']['target_particles']))
plt.ylim(-0.40, -0.25)
plt.xlabel('iteration')
plt.ylabel(r'$E(\tau) = \sum_j H_{0j} N_j(\tau)/N_0(\tau)$ / $t$')
plt.legend()
plt.tight_layout()