From 269badfbb6be0a144b9a1662daca3db9f379502e Mon Sep 17 00:00:00 2001 From: Thomas White Date: Mon, 29 Oct 2012 07:06:41 -0700 Subject: Add new script: scripts/peak-intensity --- scripts/peak-intensity | 35 +++++++++++++++++++++++++++++++++++ 1 file changed, 35 insertions(+) create mode 100755 scripts/peak-intensity 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 = ) { + + 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); -- cgit v1.2.3