def local_standard_deviations_target_per_site( unit_cell, density_map, weight_map, weight_map_scale_factor, sites_cart, site_radii): if (weight_map is None): return maptbx.standard_deviations_around_sites( unit_cell=unit_cell, density_map=density_map, sites_cart=sites_cart, site_radii=site_radii) d = maptbx.real_space_target_simple_per_site( unit_cell=unit_cell, density_map=density_map, sites_cart=sites_cart) w = flex.pow2(maptbx.standard_deviations_around_sites( unit_cell=unit_cell, density_map=weight_map, sites_cart=sites_cart, site_radii=site_radii)) w_min = 0.01 w.set_selected((w < w_min), w_min) w = 1. / w if (weight_map_scale_factor is not None): assert weight_map_scale_factor > 0 w *= weight_map_scale_factor return d / w
def get(): return maptbx.standard_deviations_around_sites( unit_cell=unit_cell, density_map=density_map, sites_cart=sites_cart, site_radii=site_radii)