diff options
| author | Thomas White <taw@physics.org> | 2013-11-29 16:06:32 +0100 |
|---|---|---|
| committer | Thomas White <taw@physics.org> | 2013-11-29 16:06:32 +0100 |
| commit | ec1fce9b2bfa07458f2c0d17bbbe86c449e7743d (patch) | |
| tree | 03f951a2a2c0a43678aa2a621400f2d7da5fbe4f | |
| parent | 15a538e5a0b66b590a8857b27e8e79d11d6174d9 (diff) | |
indexamajig: Add --int-diag=strong
| -rw-r--r-- | libcrystfel/src/integration.c | 7 | ||||
| -rw-r--r-- | libcrystfel/src/integration.h | 3 | ||||
| -rw-r--r-- | src/indexamajig.c | 4 |
3 files changed, 13 insertions, 1 deletions
diff --git a/libcrystfel/src/integration.c b/libcrystfel/src/integration.c index 042fcf7b..b51ef7f7 100644 --- a/libcrystfel/src/integration.c +++ b/libcrystfel/src/integration.c @@ -1310,6 +1310,13 @@ static int get_int_diag(struct intcontext *ic, Reflection *refl) return i < -5.0*sigi; } + if ( ic->int_diag == INTDIAG_STRONG ) { + double i, sigi; + i = get_intensity(refl); + sigi = get_esd_intensity(refl); + return i > 3.0*sigi; + } + if ( ic->int_diag == INTDIAG_INDICES ) { signed int h, k, l; get_indices(refl, &h, &k, &l); diff --git a/libcrystfel/src/integration.h b/libcrystfel/src/integration.h index ec674128..e023318a 100644 --- a/libcrystfel/src/integration.h +++ b/libcrystfel/src/integration.h @@ -41,7 +41,8 @@ typedef enum { INTDIAG_ALL, INTDIAG_INDICES, INTDIAG_NEGATIVE, - INTDIAG_IMPLAUSIBLE + INTDIAG_IMPLAUSIBLE, + INTDIAG_STRONG } IntDiag; diff --git a/src/indexamajig.c b/src/indexamajig.c index 8806d4c2..7615872a 100644 --- a/src/indexamajig.c +++ b/src/indexamajig.c @@ -556,6 +556,10 @@ int main(int argc, char *argv[]) iargs.int_diag = INTDIAG_IMPLAUSIBLE; } + if ( strcmp(int_diag, "strong") == 0 ) { + iargs.int_diag = INTDIAG_STRONG; + } + r = sscanf(int_diag, "%i,%i,%i", &h, &k, &l); if ( r == 3 ) { iargs.int_diag = INTDIAG_INDICES; |