aboutsummaryrefslogtreecommitdiff
path: root/src/scaling-report.c
diff options
context:
space:
mode:
authorThomas White <taw@physics.org>2011-10-06 10:51:10 +0200
committerThomas White <taw@physics.org>2012-02-22 15:27:38 +0100
commit6a39404f3500f5c0475030219d5b1d6fc93731ca (patch)
treead1d76da3c4f0377cea40784bedc1c093f4105b0 /src/scaling-report.c
parent08fe40822656e153e300f4cd54747bde3888c6ea (diff)
Show sample standard deviation on intensity histograms
Not the standard deviation of the mean, i.e. multiply by sqrt(N) for display
Diffstat (limited to 'src/scaling-report.c')
-rw-r--r--src/scaling-report.c8
1 files changed, 4 insertions, 4 deletions
diff --git a/src/scaling-report.c b/src/scaling-report.c
index 1aadad2c..805a131d 100644
--- a/src/scaling-report.c
+++ b/src/scaling-report.c
@@ -526,17 +526,17 @@ static void intensity_histogram(cairo_t *cr, const struct image *images,
const double g_height = 55.0;
char tmp[64];
Reflection *f;
- double Ifull, esd_Ifull, pos, mI, bit;
+ double Ifull, dsd_Ifull, pos, mI, bit;
int have_full;
f = find_refl(full, h, k, l);
if ( f != NULL ) {
Ifull = get_intensity(f);
- esd_Ifull = get_esd_intensity(f);
+ dsd_Ifull = get_esd_intensity(f) * sqrt(get_redundancy(f));
have_full = 1;
} else {
Ifull = 0.0;
- esd_Ifull = 0.0;
+ dsd_Ifull = 0.0;
have_full = 0;
}
@@ -657,7 +657,7 @@ static void intensity_histogram(cairo_t *cr, const struct image *images,
if ( have_full ) {
- double eW = g_width*esd_Ifull/mI;
+ double eW = g_width*dsd_Ifull/mI;
cairo_new_path(cr);
cairo_rectangle(cr, 0.0, g_height+bit*2.0,