with open(HUGSRoot, 'rb') as f:
    HUGSdata = pkl.load(f)