From ec97a9b53bc56ec7789a5449959a9a7722abf681 Mon Sep 17 00:00:00 2001
From: Valerio Mariani <valerio.mariani@desy.de>
Date: Mon, 26 Jan 2015 16:34:46 +0100
Subject: Add image_reflection_closest()

---
 libcrystfel/src/image.c | 49 +++++++++++++++++++++++++++++++++++++++++++++++++
 libcrystfel/src/image.h |  4 ++++
 2 files changed, 53 insertions(+)

diff --git a/libcrystfel/src/image.c b/libcrystfel/src/image.c
index 6caeeb3c..3cabef96 100644
--- a/libcrystfel/src/image.c
+++ b/libcrystfel/src/image.c
@@ -148,6 +148,55 @@ struct imagefeature *image_feature_closest(ImageFeatureList *flist,
 }
 
 
+Reflection *image_reflection_closest(RefList *rlist,
+                                     double fs, double ss,
+                                     struct detector *det,
+                                     double *d)
+{
+
+	double dmin = HUGE_VAL;
+	Reflection *closest = NULL;
+	struct panel *p1;
+	Reflection *refl;
+	RefListIterator *iter;
+
+	p1 = find_panel(det, fs, ss);
+
+	for ( refl = first_refl(rlist, &iter);
+	      refl != NULL;
+	      refl = next_refl(refl, iter) )
+	{
+		double ds;
+		struct panel *p2;
+		double rfs, rss;
+
+		get_detector_pos(refl, &rfs, &rss);
+
+		p2 = find_panel(det, rfs, rss);
+
+		if ( p1 != p2 ) {
+			continue;
+		}
+
+		ds = distance(rfs, rss, fs, ss);
+
+		if ( ds < dmin ) {
+			dmin = ds;
+			closest = refl;
+		}
+
+	}
+
+	if ( dmin < +HUGE_VAL ) {
+		*d = dmin;
+		return closest;
+	}
+
+	*d = +INFINITY;
+	return NULL;
+}
+
+
 int image_feature_count(ImageFeatureList *flist)
 {
 	if ( flist == NULL ) return 0;
diff --git a/libcrystfel/src/image.h b/libcrystfel/src/image.h
index 2498b5c2..e196d052 100644
--- a/libcrystfel/src/image.h
+++ b/libcrystfel/src/image.h
@@ -224,6 +224,10 @@ extern struct imagefeature *image_feature_closest(ImageFeatureList *flist,
                                                   double fs, double ss,
                                                   double *d, int *idx,
                                                   struct detector *det);
+extern Reflection *image_reflection_closest(RefList *rlist,
+                                            double fs, double ss,
+                                            struct detector *det,
+                                            double *d);
 
 extern int image_feature_count(ImageFeatureList *flist);
 extern struct imagefeature *image_get_feature(ImageFeatureList *flist, int idx);
-- 
cgit v1.2.3