From a9d6c673d16ea9275becaa0aaa42a1d1b347a03a Mon Sep 17 00:00:00 2001 From: othocaeS Date: Mon, 20 Jun 2016 14:26:27 -0400 Subject: [PATCH] lc_plot should have really been power_lag_plot, lc_plot will be light curve plot --- scripts/power_lag_plot.pl | 205 ++++++++++++++++++++++++++++++++++++++ 1 file changed, 205 insertions(+) create mode 100644 scripts/power_lag_plot.pl diff --git a/scripts/power_lag_plot.pl b/scripts/power_lag_plot.pl new file mode 100644 index 0000000..e165487 --- /dev/null +++ b/scripts/power_lag_plot.pl @@ -0,0 +1,205 @@ +#!/usr/local/bin/perl + +use utf8; +use Encode qw(encode decode); +use feature 'say'; +use locale; +use Switch; + + +# Enables debug output. +our $debug=0; + + +# This section locate the output data of interest in a +# psdlab output file. +if (!${^UTF8LOCALE}) { + say encode($charset,"You are not using UTF-8 encoding. :("); +} +my $charset=$ENV{LANG}; +our $outputfilename=decode($charset,$ARGV[0]); +open $outputfile,'<',$outputfilename or die $!; + +my $star_linenum_1=0; +my $star_linenum_2=0; +my $star_bytenum_1=0; +my $star_bytenum_2=0; + +my $linenum=0; +while(<$outputfile>) { + $linenum++; + if ($_ =~ /^\*+$/) { + $star_linenum_1=$star_linenum_2; + $star_linenum_2 = $linenum; + $star_bytenum_1=$star_bytenum_2; + #$star_bytenum_2 = $outputfile.tell(); + $star_bytenum_2 = tell($outputfile); + } +} +if ($debug) { + say encode($charset,"Final set found between lines "), + $star_linenum_1, + " and ", + $star_linenum_2; + + say encode($charset,"Final set found between bytes "), + $star_bytenum_1, + " and ", + $star_bytenum_2; +} +seek($outputfile,$star_bytenum_1,0); +<$outputfile>; +$bin_bounds_line = <$outputfile>; +#say $bin_bounds_line; +$bin_bounds_line =~ s/^#\s*(.*)$/$1/; +#$line =~ / +@bin_bounds = split /\s/,$bin_bounds_line; +if ($debug) { + print encode($charset,"Found bin boundaries: "); + foreach (@bin_bounds) {print encode($charset,"$_ ");} + say encode($charset," "); +} + + +# These arrays store the various quantities of interest. +@freq_coords_avg = []; +@freq_coords_err = []; +@power_curve_source_avg = []; +@power_curve_source_err = []; +@power_curve_reprocessed_avg = []; +@power_curve_reprocessed_err = []; +@cross_correlation_power_curve_avg = []; +@cross_correlation_power_curve_err = []; +@phase_difference_avg = []; +@phase_difference_err = []; +@time_lag_avg = []; +@time_lag_err = []; + +# This captures the frequency bins. +my $upper_bound = 0; +my $lower_bound = 0; +foreach (@bin_bounds) { + $lower_bound = $upper_bound; + $upper_bound = $_; + if ($lower_bound == 0) {next} + my $avg = ($upper_bound + $lower_bound)/2; + my $err = ($upper_bound - $lower_bound)/2; + #say ($avg,":",$err); + push(@freq_coords_avg,$avg); + push(@freq_coords_err,$err); +} + +$numbins = scalar @freq_coords_avg; +say encode($charset,"$numbins frequency boundaries captured."); + + +# This section collects the various quantities. The mode counter +# increments to designate the data being captured. +my $mode=0; +while (<$outputfile>) { + next if $_ =~ /\*+/; + switch ($mode) { + case 0 { + my $avg = my $err = $_; + $avg =~ s/^([\-\+e0-9\.]+)+\s+[\-\+e0-9\.]+\s*$/$1/; + $err =~ s/^[\-\+e0-9\.]+\s+([\-\+e0-9\.]+)\s*$/$1/; + if ($debug) { + say encode($charset,""); + say encode($charset," New Bin"); + say encode($charset,"─────────────────────────────────────────────────"); + print encode($charset,"Driving light curve power from $_"); + say encode($charset,"Average: $avg, Err: $err"); + } + push(@power_curve_source_avg,$avg); + push(@power_curve_source_err,$err); + $mode++; + } + case 1 { + my $avg = my $err = $_; + $avg =~ s/^([\-\+e0-9\.]+)+\s+[\-\+e0-9\.]+\s*$/$1/; + $err =~ s/^[\-\+e0-9\.]+\s+([\-\+e0-9\.]+)\s*$/$1/; + if ($debug) { + print encode($charset,"Reprocessed light curve power from $_"); + say encode($charset,"Average: $avg, Err: $err"); + } + push(@power_curve_reprocessed_avg,$avg); + push(@power_curve_reprocessed_err,$err); + $mode++; + } + case 2 { + my $avg = my $err = $_; + $avg =~ s/^([\-\+e0-9\.]+)+\s+[\-\+e0-9\.]+\s*$/$1/; + $err =~ s/^[\-\+e0-9\.]+\s+([\-\+e0-9\.]+)\s*$/$1/; + if ($debug) { + print encode($charset,"Cross-Correlation from $_"); + say encode($charset,"Average: $avg, Err: $err"); + } + push(@cross_correlation_power_curve_avg,$avg); + push(@cross_correlation_power_curve_err,$err); + $mode++; + } + case 3 { + my $avg = my $err = $_; + $avg =~ s/^([\-\+e0-9\.]+)\s+[\-\+e0-9\.]+\s*$/$1/; + $err =~ s/^[\-\+e0-9\.]+\s+([\-\+e0-9\.]+)\s*$/$1/; + if ($debug) { + print encode($charset,"Phase different from $_"); + say encode($charset,"Average: $avg, Err: $err"); + } + push(@phase_difference_avg,$avg); + push(@phase_difference_err,$err); + $mode = 0; + } + } # End switch +} + +$numvals = scalar @power_curve_source_avg; +say encode($charset,"Saved $numvals sets of quantities."); + + + + + + + + + +=pod + Tried to set up the program to find the boundaries on its own. + Should return to complete the program this way; just a syntax + problem. + +#while (<$outputfile>) {say $_;} +#while (<$outputfile> =~ /^[^#]+(.*)$/) {} +#if (eof($outputfile)) { +# say encode($charset,"Could not scrape values of energy bin boundaries."); +# exit 9; +#} +=cut + + + + + + +#$thing = <$outputfile>; +#say encode($charset,"$thing"); + + + + + + + +# switch ($val) { +# case 1 { print encode($charset,"number 1" )} +# case "a" { print encode($charset,"string a" )} +# case [1..10,42] { print encode($charset,"number in list" )} +# case (@array) { print encode($charset,"number in list" )} +# case /\w+/ { print encode($charset,"pattern" )} +# case qr/\w+/ { print encode($charset,"pattern" )} +# case (%hash) { print encode($charset,"entry in hash" )} +# case (\%hash) { print encode($charset,"entry in hash" )} +# case (\&sub) { print encode($charset,"arg to subroutine" )} +# else { print encode($charset,"previous case not true" )} +# } \ No newline at end of file