Add scripts/frequency and scripts/histogram

2 files changed, 66 insertions, 0 deletions

diff --git a/scripts/frequency b/scripts/frequency
new file mode 100755
index 00000000..1f9c66e6
--- /dev/null
+++ b/scripts/frequency
@@ -0,0 +1,51 @@
+#!/usr/bin/perl -w
+
+use strict;
+
+my $bins = 200;
+open(FH, $ARGV[0]);
+
+my $min = 99999999999999999999.0;
+my $max = 0.0;
+my $line;
+
+while ( $line = <FH> ) {
+
+	chomp($line);
+	if ( $line < $min ) {
+		$min = $line;
+	}
+	if ( $line > $max ) {
+		$max = $line;
+	}
+
+}
+
+printf("%f -> %f\n", $min, $max);
+
+seek(FH, 0, 0);
+
+my $binsize = ($max+1-$min)/$bins;
+my @hist;
+my $i;
+
+for ( $i=0; $i<$bins; $i++ ) {
+	$hist[$i] = 0;
+}
+
+while ( $line = <FH> ) {
+
+	my $bin;
+
+	chomp($line);
+
+	$bin = int(($line-$min)/$binsize);
+	$hist[$bin]++;
+
+}
+
+for ( $i=0; $i<$bins; $i++ ) {
+	printf("%f\t%i\n", $min+(($i+0.5)*$binsize), $hist[$i]);
+}
+
+close(FH);
diff --git a/scripts/histogram b/scripts/histogram
new file mode 100755
index 00000000..329b775e
--- /dev/null
+++ b/scripts/histogram
@@ -0,0 +1,15 @@
+#!/bin/sh
+
+./frequency.pl 103.hist > hist.dat
+
+gnuplot << EOF
+set term postscript enhanced font "Helvetica,20"
+set output "histo.ps"
+
+plot [-2000:60000] "hist.dat" u 1:2 w l t "103 reflection"
+EOF
+
+ps2pdf histo.ps
+rm -f histo.ps hist.dat
+
+evince histo.pdf