import pyhande
import matplotlib.pyplot as plt
(metadata, qmc_data) = pyhande.extract.extract_data('calcs/shoulder/h2o_plat.out')[0]
(metadata, qmc_data2) = pyhande.extract.extract_data('calcs/shoulder/h2o_plat_bigstartpop.out')[0]
plt.loglog(qmc_data['# H psips'], (qmc_data['# H psips'] / qmc_data['N_0']), label='initial population = 200')
plt.loglog(qmc_data2['# H psips'], (qmc_data2['# H psips'] / qmc_data2['N_0']), label='initial population = 800')
plt.legend(loc="best")
plt.xlabel('# particles')
plt.ylabel('# particles / # particles on reference')