diff options
author | Thomas White <taw@physics.org> | 2011-12-12 13:17:30 +0100 |
---|---|---|
committer | Thomas White <taw@physics.org> | 2012-02-22 15:27:42 +0100 |
commit | 8c40aed1b3ab66fbab2a0b457ccf8286bc1a3b05 (patch) | |
tree | 3a491e2cd4b63b9b80ead4fbe01e81d03ad05e18 /libcrystfel/src/peaks.c | |
parent | 2318ae07fa5baaead4c1bc55bc4e63694c1942dc (diff) |
indexamajig: Optional I/sigma(I) cutoff during integration
Diffstat (limited to 'libcrystfel/src/peaks.c')
-rw-r--r-- | libcrystfel/src/peaks.c | 5 |
1 files changed, 4 insertions, 1 deletions
diff --git a/libcrystfel/src/peaks.c b/libcrystfel/src/peaks.c index 0d72e514..34e0b8ca 100644 --- a/libcrystfel/src/peaks.c +++ b/libcrystfel/src/peaks.c @@ -564,7 +564,7 @@ int peak_sanity_check(struct image *image) /* Integrate the list of predicted reflections in "image" */ void integrate_reflections(struct image *image, int polar, int use_closer, - int bgsub) + int bgsub, double min_snr) { Reflection *refl; RefListIterator *iter; @@ -606,6 +606,9 @@ void integrate_reflections(struct image *image, int polar, int use_closer, &intensity, &bg, &max, &sigma, polar, 0, bgsub); + /* I/sigma(I) cutoff */ + if ( intensity/sigma < min_snr ) r = 0; + /* Record intensity and set redundancy to 1 on success */ if ( r == 0 ) { set_int(refl, intensity); |