Hooks for estimation of maximum resolution for each image

1 files changed, 28 insertions, 0 deletions

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)