system_functions_ALA4-21x21grids_msstate-cut-inc_2msstates.py

import unittest
import warnings
import numpy as np
from copy import copy, deepcopy
# from PyEmma
from pyemma.msm.analysis import pcca

def ALA4_msstates( T_AA ):

    # Define the metastable states from analysis of T_AA_mle
    membership=pcca(T_AA, 2)
    # the precise definitions were fine-tuned by hand, I should plot the PCCA states to double check!
    H = np.where(membership[:,0]>0.80)[0]
    E = np.where(membership[:,1]>0.80)[0]
    # make sure there is no overlap
    #print "checking the overlap of metastable states:"
    check = [val for val in H if val in E]
    if (len(check) != 0 ):  raise ValueError('Some metastable states are overlapping!')

    return H, E

# shortcut function for msstates, should always be linked to the particular calculation (above) for each system
get_msstates = lambda T, lcc, bins: ALA4_msstates( T )