Don't calculate pobs if reference reflection is weak

2 files changed, 6 insertions, 0 deletions

diff --git a/src/partialator.c b/src/partialator.c
index f826bf83..0142a628 100644
--- a/src/partialator.c
+++ b/src/partialator.c
@@ -648,6 +648,9 @@ static void write_to_pgraph(FILE *fh, RefList *list, RefList *full, Crystal *cr,
 		match = find_refl(full, h, k, l);
 		if ( match == NULL ) continue;
 
+		/* Don't calculate pobs if reference reflection is weak */
+		if ( fabs(get_intensity(match)) / get_esd_intensity(match) < 3.0 ) continue;
+
 		/* Calculated partiality */
 		pcalc = get_partiality(refl);
 
diff --git a/src/post-refinement.c b/src/post-refinement.c
index 007216b9..79ea187a 100644
--- a/src/post-refinement.c
+++ b/src/post-refinement.c
@@ -617,6 +617,9 @@ void write_specgraph(Crystal *crystal, const RefList *full,
 		match = find_refl(full, h, k, l);
 		if ( match == NULL ) continue;
 
+		/* Don't calculate pobs if reference reflection is weak */
+		if ( fabs(get_intensity(match)) / get_esd_intensity(match) < 3.0 ) continue;
+
 		corr = G * exp(B*res*res) * get_lorentz(refl);
 		Ipart = get_intensity(refl) * corr;
 		Ifull = get_intensity(match);