diff options
| author | Thomas White <taw@physics.org> | 2019-10-24 15:29:27 +0200 |
|---|---|---|
| committer | Thomas White <taw@physics.org> | 2019-10-24 15:29:27 +0200 |
| commit | 43147db0b3c5232068f0cec70008fa217fa97476 (patch) | |
| tree | 52b7d1bcde402b0b32ea39d79ff695abe9ae83e6 | |
| parent | 049c3eb48c31dfa006dd78e2612babf6d6b3b6db (diff) | |
polarisation_correction: Scale the sigma as well
| -rw-r--r-- | libcrystfel/src/geometry.c | 4 |
1 files changed, 3 insertions, 1 deletions
diff --git a/libcrystfel/src/geometry.c b/libcrystfel/src/geometry.c index fc42a8de..3d6cd699 100644 --- a/libcrystfel/src/geometry.c +++ b/libcrystfel/src/geometry.c @@ -1042,7 +1042,7 @@ void polarisation_correction(RefList *list, UnitCell *cell, refl = next_refl(refl, iter) ) { double pol; - double intensity; + double intensity, sigma; double xl, yl, zl; double phi, tt, kpred; signed int h, k, l; @@ -1062,7 +1062,9 @@ void polarisation_correction(RefList *list, UnitCell *cell, + (1.0-P)*(1.0 - sin(phi)*sin(phi)*sin(tt)*sin(tt)); intensity = get_intensity(refl); + sigma = get_esd_intensity(refl); set_intensity(refl, intensity / pol); + set_esd_intensity(refl, sigma / pol); } } |