diff options
| author | Thomas White <taw@physics.org> | 2011-11-09 14:57:55 +0100 |
|---|---|---|
| committer | Thomas White <taw@physics.org> | 2012-02-22 15:27:39 +0100 |
| commit | e21d69a9141e3f64147aa2f68f9c3634d66a493f (patch) | |
| tree | 6f231e0aa0427c1a42c350e3908ae664fdb959e2 /src/partialator.c | |
| parent | 667c06f12ca19986acfc00b0bc3adc2ca7f58f58 (diff) | |
partialator: Use an I/sig(I) cutoff for determining scalability
Diffstat (limited to 'src/partialator.c')
| -rw-r--r-- | src/partialator.c | 5 |
1 files changed, 3 insertions, 2 deletions
diff --git a/src/partialator.c b/src/partialator.c index 4edf7bf4..ff03c039 100644 --- a/src/partialator.c +++ b/src/partialator.c @@ -155,14 +155,15 @@ static int select_scalable_reflections(RefList *list, RefList *reference) refl = next_refl(refl, iter) ) { int sc = 1; - double v; + double v, esd; /* This means the reflection was not found on the last check */ if ( get_redundancy(refl) == 0 ) sc = 0; if ( get_partiality(refl) < 0.1 ) sc = 0; v = fabs(get_intensity(refl)); - if ( v < 0.1 ) sc = 0; /* FIXME: Yuk */ + esd = get_esd_intensity(refl); + if ( v < 0.5*esd ) sc = 0; /* If we are scaling against a reference set, we additionally * require that this reflection is in the reference list. */ |