def calc_wilson(observations_full, n_residues): """ Caculate isotropic Wilson G and B-factors """ if n_residues == 0: return 0, 0 from prime.postrefine.mod_util import mx_handler mxh = mx_handler() asu_contents = mxh.get_asu_contents(n_residues) try: observations_as_f = observations_full.as_amplitude_array() binner_template_asu = observations_as_f.setup_binner(auto_binning=True) wp = statistics.wilson_plot(observations_as_f, asu_contents, e_statistics=True) G = wp.wilson_intensity_scale_factor B = wp.wilson_b except Exception: G,B = (0,0) return G, B
def calc_wilson(observations_full, n_residues): """ Caculate isotropic Wilson G and B-factors """ if n_residues == 0: return 0, 0 from prime.postrefine.mod_util import mx_handler mxh = mx_handler() asu_contents = mxh.get_asu_contents(n_residues) try: observations_as_f = observations_full.as_amplitude_array() binner_template_asu = observations_as_f.setup_binner(auto_binning=True) wp = statistics.wilson_plot(observations_as_f, asu_contents, e_statistics=True) G = wp.wilson_intensity_scale_factor B = wp.wilson_b except Exception: G, B = (0, 0) return G, B