Scripts: fitsearch.py

# Get fit search results in nogui mode.

def print_fit_search_results(map1_id = 1, map2_id = 0, search = 10):

    from chimera import openModels as om, selection
    m1 = om.list(id = map1_id)[0]
    m2 = om.list(id = map2_id)[0]
    s1 = selection.ItemizedSelection([m1])
    
    from FitMap.fitcmd import fitmap
    fit_list = fitmap(s1, m2, search = search, listFits = False)

    print
    print '%d fits' % len(fit_list)
    import Matrix
    for fit in fit_list:
        print 'correlation =', fit.correlation()
        print 'ave map value = ', fit.average_map_value()
        print 'points inside = ', fit.points_inside_contour()
        print 'times found', fit.hits()
        print Matrix.transformation_description(fit.transforms[0])


print_fit_search_results(2,0,20)