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)