diff options
author | Thomas White <taw@physics.org> | 2011-10-06 10:51:10 +0200 |
---|---|---|
committer | Thomas White <taw@physics.org> | 2012-02-22 15:27:38 +0100 |
commit | 6a39404f3500f5c0475030219d5b1d6fc93731ca (patch) | |
tree | ad1d76da3c4f0377cea40784bedc1c093f4105b0 /src/scaling-report.c | |
parent | 08fe40822656e153e300f4cd54747bde3888c6ea (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.c | 8 |
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, |