From e21d69a9141e3f64147aa2f68f9c3634d66a493f Mon Sep 17 00:00:00 2001 From: Thomas White Date: Wed, 9 Nov 2011 14:57:55 +0100 Subject: partialator: Use an I/sig(I) cutoff for determining scalability --- src/partialator.c | 5 +++-- 1 file 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. */ -- cgit v1.2.3