aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorThomas White <taw@physics.org>2014-03-06 12:07:45 +0100
committerThomas White <taw@physics.org>2014-03-06 12:07:45 +0100
commite7ed8c50e67cbc35fe63df9003a83b331373f471 (patch)
tree27a9bbe098b75ca224dcd7f99083cb9c54319006
parent533c8aad7b256cceeb1ae9cae346ce00827cc3b1 (diff)
Add scripts/fg-graph
-rw-r--r--Makefile.am2
-rwxr-xr-xscripts/fg-graph61
2 files changed, 62 insertions, 1 deletions
diff --git a/Makefile.am b/Makefile.am
index 0a187d94..bdbbb829 100644
--- a/Makefile.am
+++ b/Makefile.am
@@ -176,7 +176,7 @@ script_DATA = scripts/alternate-stream scripts/cell-please \
scripts/random-image scripts/README scripts/sequence-image \
scripts/split-indexed scripts/stream_grep scripts/zone-axes \
scripts/Rsplit_surface scripts/Rsplit_surface.py \
- scripts/clean-stream.py
+ scripts/clean-stream.py scripts/fg-graph
EXTRA_DIST += $(script_DATA)
diff --git a/scripts/fg-graph b/scripts/fg-graph
new file mode 100755
index 00000000..701277c3
--- /dev/null
+++ b/scripts/fg-graph
@@ -0,0 +1,61 @@
+#!/bin/bash
+
+# Copyright © 2014 Wolfgang Brehm
+# Copyright © 2014 Deutsches Elektronen-Synchrotron DESY,
+# a research centre of the Helmholtz Association.
+#
+# Authors:
+# 2013-2014 Wolfgang Brehm <wolfgang.brehm@gmail.com>
+# 2014 Thomas White <taw@physics.org>
+#
+# This file is part of CrystFEL.
+#
+# CrystFEL is free software: you can redistribute it and/or modify
+# it under the terms of the GNU General Public License as published by
+# the Free Software Foundation, either version 3 of the License, or
+# (at your option) any later version.
+#
+# CrystFEL is distributed in the hope that it will be useful,
+# but WITHOUT ANY WARRANTY; without even the implied warranty of
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+# GNU General Public License for more details.
+#
+# You should have received a copy of the GNU General Public License
+# along with CrystFEL. If not, see <http://www.gnu.org/licenses/>.
+#
+
+FILE=fg.dat
+NPATT=4000
+NITER=4
+CMIN=-0.1
+CMAX=0.3
+
+gnuplot -p << EOF
+set terminal pngcairo size 1600,1000 enhanced monochrome font 'Helvetica,20' linewidth 2
+set output "correlation.png"
+set xlabel "Iteration i"
+set ylabel "Correlation"
+rnd(x) = x - floor(x) < 0.5 ? floor(x) : ceil(x)
+round(x1,x2)=rnd(10**x2*x1)/10.0**x2 # funktion mit zwei Argumenten, x1 wird gerundet x2 anzahl der Stellen nach dem Komma
+set xzeroaxis lc rgb "black" lt 1
+set key top left
+set xrange [0:${NITER}*${NPATT}]
+set yrange [${CMIN}:${CMAX}]
+set ytics nomirror
+set xtics nomirror
+set x2tics nomirror
+
+set x2range [0:${NITER}]
+set x2label "Iteration i gemessen in n"
+set arrow from ${NPATT},${CMIN} to ${NPATT},${CMAX} nohead lc rgb "black"
+
+plot\
+ 1 lc rgb "black" lt 1 notitle,\
+ "${FILE}" u 0:1 ps 0.2 pt 7 lc rgb "orange" notitle,\
+ "${FILE}" u 0:2 ps 0.2 pt 7 lc rgb "#0088FF" notitle,\
+ "${FILE}" u 0:(rand(0)>0.5?\$1:1/0) ps 0.2 pt 7 lc rgb "orange" notitle ,\
+ "${FILE}" u (round(\$0, -3)):(\$1>\$2?\$1:\$2) lw 3 lc rgb "black" smooth unique notitle,\
+ "${FILE}" u (round(\$0, -3)):(\$1<\$2?\$1:\$2) lw 3 lc rgb "black" smooth unique notitle,\
+ "${FILE}" u (round(\$0, -3)):1 lw 3 lc rgb "red" smooth unique notitle,\
+ "${FILE}" u (round(\$0, -3)):2 lw 3 lc rgb "blue" smooth unique notitle
+EOF