From b2c5b05c7e8da5a20610c7a7b098a76c9ad9c4e0 Mon Sep 17 00:00:00 2001 From: Thomas White Date: Mon, 30 Mar 2015 16:51:13 +0200 Subject: Add scripts/ave-resolution --- scripts/ave-resolution | 35 +++++++++++++++++++++++++++++++++++ 1 file changed, 35 insertions(+) create mode 100755 scripts/ave-resolution (limited to 'scripts') diff --git a/scripts/ave-resolution b/scripts/ave-resolution new file mode 100755 index 00000000..ff93ed68 --- /dev/null +++ b/scripts/ave-resolution @@ -0,0 +1,35 @@ +#!/usr/bin/env python +# -*- coding: utf-8 -*- + +# +# Find mean diffracting resolution +# +# Copyright © 2014-2015 Deutsches Elektronen-Synchrotron DESY, +# a research centre of the Helmholtz Association. +# +# Author: +# 2014-2015 Thomas White +# + +import sys +import numpy + +f = open(sys.argv[1]) +a = [] + +while True: + fline = f.readline() + if not fline: + break + if fline.find("diffraction_resolution_limit") != -1: + res = float(fline.split('= ')[1].split(' ')[0].rstrip("\r\n")) + a.append(res) + continue + +f.close() + +b = numpy.array(a) +print " Mean: %.2f nm^-1 = %.2f A" % (numpy.mean(b),10.0/numpy.mean(b)) +print " Best: %.2f nm^-1 = %.2f A" % (numpy.max(b),10.0/numpy.max(b)) +print "Worst: %.2f nm^-1 = %.2f A" % (numpy.min(b),10.0/numpy.min(b)) +print "Std deviation: %.2f nm^-1" % (numpy.std(b)) -- cgit v1.2.3