I need to calculate the free energy from my MD simulations which respect contact frequency and RMSD as my reaction coordinates. I have used the pyWham library, and I want to know more flexible python based code for free energy calculations which I can use for my analysis code. Thank you!!!