From d6bfd45c3cd48023532fcc40021a8c98d18e05bf Mon Sep 17 00:00:00 2001 From: Thomas White Date: Thu, 7 Jan 2021 11:29:23 +0100 Subject: Remove "Bragg proximity check" during integration This test is meant to avoid integrating overlapping spots. It does a very bad job of this, because it doesn't take into account whether the neighbouring reflection is actually excited, or even allowed by the lattice. It's even worse when we think about wide bandwidth. --- libcrystfel/src/integration.c | 17 ----------------- 1 file changed, 17 deletions(-) (limited to 'libcrystfel/src/integration.c') diff --git a/libcrystfel/src/integration.c b/libcrystfel/src/integration.c index d859b45d..0ae19cf0 100644 --- a/libcrystfel/src/integration.c +++ b/libcrystfel/src/integration.c @@ -836,9 +836,6 @@ static int check_box(struct intcontext *ic, struct peak_box *bx, int *sat) for ( q=0; qw; q++ ) { int fs, ss; - double hd, kd, ld; - signed int h, k, l; - double dv[3]; float lsat; fs = bx->cfs + p; @@ -894,20 +891,6 @@ static int check_box(struct intcontext *ic, struct peak_box *bx, int *sat) if ( sat != NULL ) *sat = 1; } - /* Ignore if this pixel is closer to the next reciprocal lattice - * point */ - detgeom_transform_coords(bx->p, fs, ss, - 1.0/ic->k, dv); - hd = dv[0] * adx + dv[1] * ady + dv[2] * adz; - kd = dv[0] * bdx + dv[1] * bdy + dv[2] * bdz; - ld = dv[0] * cdx + dv[1] * cdy + dv[2] * cdz; - h = lrint(hd); - k = lrint(kd); - l = lrint(ld); - if ( (h != hr) || (k != kr) || (l != lr) ) { - bx->bm[p+ic->w*q] = BM_BH; - } - /* Find brightest pixel */ if ( (bx->bm[p+ic->w*q] != BM_IG) && (bx->bm[p+ic->w*q] != BM_BH) -- cgit v1.2.3