aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rwxr-xr-xscripts/peak-intensity35
1 files changed, 35 insertions, 0 deletions
diff --git a/scripts/peak-intensity b/scripts/peak-intensity
new file mode 100755
index 00000000..46f5be9e
--- /dev/null
+++ b/scripts/peak-intensity
@@ -0,0 +1,35 @@
+#!/usr/bin/perl -w
+
+use strict;
+
+open(FH, $ARGV[0]);
+
+my $line;
+my $total_i;
+my $n = 0;
+my $n_patt = 0;
+
+while ( $line = <FH> ) {
+
+ if ( $line =~ /^-----\ Begin chunk\ -----$/ ) {
+ $n_patt++;
+ }
+
+ if ( $line =~ /^\s*([\d\.]+)\s+([\d\.]+)\s+([\d\.]+)\s+([\d\.]+)$/ ) {
+
+ my $fs = $1;
+ my $ss = $2;
+ my $one_over_d = $3;
+ my $i = $4;
+
+ $total_i += $i;
+ $n++;
+
+ }
+
+}
+
+printf("%i patterns, %i peaks, %.2f total ADU\n", $n_patt, $n, $total_i);
+printf("Mean %i peaks per hit\n", $n / $n_patt);
+printf("Mean %.2f ADU per peak\n", $total_i / $n);
+printf("Mean %.2f ADU per hit\n", $total_i / $n_patt);