[Chimera-users] representative fits in chimera fitmap search
Jan Kosinski
kosinski at embl.de
Fri Mar 18 10:19:59 PDT 2016
Hi,
In the fitmap command global search, how the unique fit is selected from
the cluster of similar fits?
Ideally, from each cluster, I would like to obtain the fit the gave the
best cross-correlation. But, I think, based on the fit_search code in
FitMap/search.py, what is happening is that the first unique fit ever
encountered is added to the fit list:
close = b.close_transforms(ptf)
if len(close) == 0:
transforms = [M.multiply_matrices(ptf, mtv) for mtv in
mtv_list]
stats['hits'] = 1
f = Fit(models, transforms, volume, stats)
f.ptf = ptf
flist.append(f)
b.add_transform(ptf)
fo[id(ptf)] = f
else:
s = fo[id(close[0])].stats
s['hits'] += 1
and any subsequent fit that is close to the first, would be discarded,
even if it gives better cross-correlation.
Do I understand correctly that the unique fit is not necessarily the
best scoring fit from the cluster?
Thanks in advance for clarification,
Jan
More information about the Chimera-users
mailing list