diff options
| author | Thomas White <taw@physics.org> | 2023-08-30 17:14:21 +0200 |
|---|---|---|
| committer | Thomas White <taw@physics.org> | 2023-08-30 17:14:21 +0200 |
| commit | 63111bd3c29a20657f301601330acf70b3de42a1 (patch) | |
| tree | 2193e20f133349fe9643f677407b88dd25a36d40 /libcrystfel/src/crystfel-mille.c | |
| parent | 4b1ad5c62bfd203503a68b136453c810ba1bb11f (diff) | |
Add excitation error to Millepede data (local-only measurement)
Diffstat (limited to 'libcrystfel/src/crystfel-mille.c')
| -rw-r--r-- | libcrystfel/src/crystfel-mille.c | 7 |
1 files changed, 7 insertions, 0 deletions
diff --git a/libcrystfel/src/crystfel-mille.c b/libcrystfel/src/crystfel-mille.c index fd957693..6a80e323 100644 --- a/libcrystfel/src/crystfel-mille.c +++ b/libcrystfel/src/crystfel-mille.c @@ -186,6 +186,7 @@ void write_mille(Mille *mille, int n, UnitCell *cell, float local_gradients_fs[nl]; float local_gradients_ss[nl]; + float local_gradients_r[nl]; float global_gradients_fs[ng*max_hierarchy_levels]; float global_gradients_ss[ng*max_hierarchy_levels]; int labels[ng*max_hierarchy_levels]; @@ -199,6 +200,8 @@ void write_mille(Mille *mille, int n, UnitCell *cell, Minvs[rps[i].peak->pn], 0, 0, 0, &local_gradients_fs[j], &local_gradients_ss[j]); + local_gradients_r[j] = EXC_WEIGHT * r_gradient(rvl[j], rps[i].refl, + cell, image->lambda); } /* Global gradients for each hierarchy level, starting at the @@ -243,6 +246,10 @@ void write_mille(Mille *mille, int n, UnitCell *cell, nl, local_gradients_ss, j, global_gradients_ss, labels, ss_dev(&rps[i], image->detgeom), 0.22); + + /* Add excitation error "measurement" (local-only) */ + mille_add_measurement(mille, nl, local_gradients_r, + 0, NULL, NULL, r_dev(&rps[i])*EXC_WEIGHT, 1.0); } } |