aboutsummaryrefslogtreecommitdiff
path: root/scripts/peak-intensity
blob: 46f5be9edbfe84b70819e2ec9725c8312f8603a1 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
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);