diff options
| author | Thomas White <taw@physics.org> | 2021-01-07 11:29:23 +0100 |
|---|---|---|
| committer | Thomas White <taw@physics.org> | 2021-01-08 10:52:46 +0100 |
| commit | d6bfd45c3cd48023532fcc40021a8c98d18e05bf (patch) | |
| tree | 7b1656bf63606a572e002c3dff515c522be82a85 /libcrystfel/src/integration.c | |
| parent | 845c338b3032837ce37e02fbed0a82879030537d (diff) | |
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.
Diffstat (limited to 'libcrystfel/src/integration.c')
| -rw-r--r-- | libcrystfel/src/integration.c | 17 |
1 files changed, 0 insertions, 17 deletions
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; q<ic->w; 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) |