diff options
| author | Thomas White <taw@physics.org> | 2015-03-18 18:01:56 +0100 |
|---|---|---|
| committer | Thomas White <taw@physics.org> | 2015-04-20 15:50:39 +0200 |
| commit | 4745fc71282e7383063898ddad8248dd98dd6c17 (patch) | |
| tree | 36bb0c9e426934838646132ab211f8409409e4a2 | |
| parent | 6c677e891d1da131c6f985e84e6cc8dd7a9ebde5 (diff) | |
Add plot-predict-refine and plot-radius-resolution
| -rw-r--r-- | Makefile.am | 3 | ||||
| -rwxr-xr-x | scripts/plot-predict-refine | 11 | ||||
| -rwxr-xr-x | scripts/plot-radius-resolution | 12 |
3 files changed, 25 insertions, 1 deletions
diff --git a/Makefile.am b/Makefile.am index b340dc86..aa479254 100644 --- a/Makefile.am +++ b/Makefile.am @@ -186,7 +186,8 @@ script_DATA = scripts/alternate-stream scripts/cell-please \ scripts/clean-stream.py scripts/fg-graph scripts/truncate-stream \ scripts/gen-sfs-expand scripts/add-beam-params \ scripts/find-pairs scripts/plot-cc-and-scale.R \ - scripts/ave-resolution scripts/crystal-frame-number + scripts/ave-resolution scripts/crystal-frame-number \ + scripts/plot-radius-resolution scripts/plot-predict-refine EXTRA_DIST += $(script_DATA) diff --git a/scripts/plot-predict-refine b/scripts/plot-predict-refine new file mode 100755 index 00000000..68baffdc --- /dev/null +++ b/scripts/plot-predict-refine @@ -0,0 +1,11 @@ +#!/bin/sh + +INFILE=$1 + +grep "predict_refine/R" $INFILE > plotme.dat +gnuplot -persist << EOF +set xlabel "Profile radius before refinement / nm^-1" +set ylabel "Profile radius after refinement / nm^-1" +plot "plotme.dat" using 4:7 +replot x +EOF diff --git a/scripts/plot-radius-resolution b/scripts/plot-radius-resolution new file mode 100755 index 00000000..da7a047d --- /dev/null +++ b/scripts/plot-radius-resolution @@ -0,0 +1,12 @@ +#!/bin/sh + +INFILE=$1 + +grep "profile_radius" $INFILE > plotme1.dat +grep "diffraction_resolution_limit" $INFILE > plotme2.dat +paste plotme1.dat plotme2.dat > plotme.dat +gnuplot -persist << EOF +set xlabel "Estimated crystal diameter (= 2/profile radius) / nm" +set ylabel "Resolution limit / nm^-1" +plot "plotme.dat" using (2/\$3):7 +EOF |