From 2318ae07fa5baaead4c1bc55bc4e63694c1942dc Mon Sep 17 00:00:00 2001 From: Thomas White Date: Sun, 28 Aug 2011 04:29:43 -0700 Subject: Hooks for estimation of maximum resolution for each image --- libcrystfel/src/peaks.c | 28 ++++++++++++++++++++++++++++ 1 file changed, 28 insertions(+) (limited to 'libcrystfel/src/peaks.c') diff --git a/libcrystfel/src/peaks.c b/libcrystfel/src/peaks.c index 82ce1cf3..0d72e514 100644 --- a/libcrystfel/src/peaks.c +++ b/libcrystfel/src/peaks.c @@ -290,6 +290,34 @@ int integrate_peak(struct image *image, int cfs, int css, } +void estimate_resolution(RefList *list, UnitCell *cell, + double *min, double *max) +{ + Reflection *refl; + RefListIterator *iter; + + for ( refl = first_refl(list, &iter); + refl != NULL; + refl = next_refl(refl, iter) ) + { + double one_over_d; + signed int h, k, l; + + get_indices(refl, &h, &k, &l); + one_over_d = 2.0 * resolution(cell, h, k, l); + + if ( one_over_d > *max ) *max = one_over_d; + if ( one_over_d < *min ) *min = one_over_d; + + /* FIXME: Implement this */ + + } + + *min = 0.0; + *max = 0.0; +} + + static void search_peaks_in_panel(struct image *image, float threshold, float min_gradient, float min_snr, struct panel *p) -- cgit v1.2.3