Skip to content
Snippets Groups Projects
Select Git revision
  • 9e55b84fc18dd26cf79c4e8bdef7f28f930a6f0b
  • develop default
  • Blockage
  • GuidoBasten
  • Anticipation_Model
  • tgf19ts
  • 313-waiting-behaviour
  • tgf19mc
  • wa_testing
  • trajectories
  • 287-waitingarea
  • 320-one-sided-closed-doors
  • kapakrit-anna
  • 311-improve-ff
  • 306-schedule-for-was
  • 294-events
  • 307-temp-close
  • 302-external-files
  • split_files
  • 298-sources
  • 293-statistics-crossings
  • v0.8.4
  • v0.8.3
  • v0.8.2
  • v0.8.1
  • v0.8
  • v0.7
  • v0.6
  • v0.5-alpha
  • v0.5-alpha1
  • v0.4
31 results

ffRouter.cpp

Blame
  • example_contact_mapping.py 779 B
    from src.BioHelpers_FABER.cmap import Cmap
    from src.BioHelpers_FABER.visualisation import show_contact_map
    import src.BioHelpers_FABER.bio_mod as bm
    import Bio.PDB as PDB
    import matplotlib.pyplot as plt
    from matplotlib.patches import Rectangle
    
    
    def main():
        test_pdb = Cmap()
        test_pdb.load_native_pdb("1l9a.pdb", start_from_one=False, added_tail=0)
    
        print(test_pdb.positioning, test_pdb.l, test_pdb.contacts["native"].shape)
        print(test_pdb.get_unmodelled_regions())
        test_pdb.load_contacts_from_restraints("simrna_1l9a_L0.5.res", "DCA")
        print(test_pdb.ppvs["DCA"])
        # print(test_pdb.get_unmodelled_regions(), test_pdb.ppvs["DCA"])
        fig, ax = show_contact_map(test_pdb, ["DCA"], mark_match=True)
        plt.show()
    
    
    if __name__ == "__main__":
        main()