From 46a290bfc8de65bef253f8516e13b3ca09698caa Mon Sep 17 00:00:00 2001 From: Thomas White Date: Tue, 19 Nov 2013 17:44:48 +0100 Subject: Move debugging output --- libcrystfel/src/integration.c | 13 +++++++------ 1 file changed, 7 insertions(+), 6 deletions(-) (limited to 'libcrystfel/src/integration.c') diff --git a/libcrystfel/src/integration.c b/libcrystfel/src/integration.c index 7046e109..a118384a 100644 --- a/libcrystfel/src/integration.c +++ b/libcrystfel/src/integration.c @@ -1564,12 +1564,6 @@ static void integrate_rings_once(Reflection *refl, struct image *image, intensity = tentative_intensity(ic, bx); mean_var_area(ic, bx, BM_BG, &bgmean, &sig2_bg); - mean_var_area(ic, bx, BM_PK, &pkmean, &sig2_pk); - if ( bx->verbose ) { - STATUS("bg mean, var = %.2f, %.2f\n", bgmean, sig2_bg); - STATUS("pk mean, var = %.2f, %.2f\n", pkmean, sig2_pk); - } - aduph = bx->p->adu_per_eV * ph_lambda_to_eV(ic->image->lambda); sig2_poisson = aduph * intensity; @@ -1595,6 +1589,13 @@ static void integrate_rings_once(Reflection *refl, struct image *image, sigma = sqrt(sig2_poisson + bx->m*sig2_bg); + mean_var_area(ic, bx, BM_PK, &pkmean, &sig2_pk); + if ( bx->verbose ) { + STATUS("bg mean, var = %.2f, %.2f\n", bgmean, sig2_bg); + STATUS("pk mean, var = %.2f, %.2f\n", pkmean, sig2_pk); + STATUS("intensity = %.2f +/- %.2f\n", intensity, sigma); + } + if ( intensity < -5.0*sigma ) { delete_box(ic, bx); ic->n_implausible++; -- cgit v1.2.3