From 940c8f05cde9b2df0ccbcb970e6959c7d3c043fb Mon Sep 17 00:00:00 2001
From: Anna Bergamaschi <anna.bergamaschi@psi.ch>
Date: Fri, 24 Jan 2020 15:32:48 +0100
Subject: [PATCH] interpolation without quads

---
 .../moenchExecutables/moench03Interpolation.cpp           | 8 +++++++-
 1 file changed, 7 insertions(+), 1 deletion(-)

diff --git a/slsDetectorCalibration/moenchExecutables/moench03Interpolation.cpp b/slsDetectorCalibration/moenchExecutables/moench03Interpolation.cpp
index 76ce75739..1032256a8 100644
--- a/slsDetectorCalibration/moenchExecutables/moench03Interpolation.cpp
+++ b/slsDetectorCalibration/moenchExecutables/moench03Interpolation.cpp
@@ -5,7 +5,7 @@
 //#include "moench03T1ZmqData.h"
 //#define DOUBLE_SPH
 //#define MANYFILES
-#define WRITE_QUAD
+//#define WRITE_QUAD
 #ifdef DOUBLE_SPH
 #include "single_photon_hit_double.h"
 #endif
@@ -177,6 +177,12 @@ int main(int argc, char *argv[]) {
 // 	    if (cl.x>210 && cl.x<240 && cl.y>210 && cl.y<240) {
 // #endif
 	    if (cl.x<0 || cl.y<0 || cl.x>NC || cl.y>NR) {
+	      if (cl.x<-1 || cl.y<-1 || cl.x>NC+1 || cl.y>NR+1) {
+		cout <<"**************"<< endl;
+		cout << cl.x << " " << cl.y << " " << sum << endl;
+		cl.print();
+		cout <<"**************"<< endl;
+	      }
 	    } else  {
 #ifndef FF
 	      // interp->getInterpolatedPosition(cl.x,cl.y, cl.get_cluster(),int_x, int_y);