paramagpy.fit.gridsearch_fit_atom_restrain_distance_cutoff

paramagpy.fit.gridsearch_fit_atom_restrain_distance_cutoff(densityMapA, densityMapB, distUpper, distLower, cutoffValue)[source]

Given two RMSD density maps, this function will compare all points pairwise and return only those within the bounds of a given distance cutoff and within a given RMSD cutoff. This might be useful if two density maps for separate atoms in a molecule are known to be constrained w.r.t. one another and you would like to use that restraint to further restrict the space of PCS RMSD points. The calculation first selects points with an RMSD less than the given cutoff value and then compares each point pariwise to fulfill the distance condition. It then returns those points from both maps. Unfortunately there is no correlation data available between these two maps.

Parameters

  • densityMapA (paramagpy.fit.DensityMap) – a density map of PCS RMSD values.

  • densityMapB (paramagpy.fit.DensityMap) – a second density map of PCS RMSD values.

  • distUpper (float) – The upper limit for distance. Any pairwise distances larger than this value will be rejected from the final space of points

  • distLower (float) – The lower distance limit for distance Any pairwise distance smaller than this value will be rejectet from the final space of points

  • cutoffValue (float) – The RMSD threshold for which points are taken for the pairwise distance comparison. Note that the total number of points returned is significnalty influenced by this parameter and should be chosen carefully

Returns

tuple – two lists of [x,y,z] coordinates are returned associated with the inputs <densityMapA> and <densityMapB>. The returned coordinates are taken from the original grids and represent points that have another associated point in the other grid which is within the distance bounds and have an RMSD below the cutoff threshold.

Return type

np.ndarray of position coordinates