import pandas as pd
import matplotlib.pyplot as plt
data = pd.read_csv('calcs/dmqmc/heisenberg_reweighted_block.out', sep=r'\s+')
data = data[::20]
plt.errorbar(data['Beta'], data['Tr[Hp]/Tr[p]']/36., yerr=data['Tr[Hp]/Tr[p]_error']/36., fmt='s')
plt.xlabel(r'$\beta J$')
plt.ylabel(r'$U/(JN)$')