From 27a61c97bb72b67e44c3d27d27a25114a9c5c437 Mon Sep 17 00:00:00 2001
From: caes <xaw5719@wmich.edu>
Date: Thu, 30 Jun 2016 02:51:19 -0400
Subject: [PATCH] created script for converting heasoft lightcurves to psdlag
 lightcurves for a single filter

---
 scripts/heasoft_to_psdlag.pl | 74 ++++++++++++++++++++++++++++++++++++
 1 file changed, 74 insertions(+)
 create mode 100755 scripts/heasoft_to_psdlag.pl

diff --git a/scripts/heasoft_to_psdlag.pl b/scripts/heasoft_to_psdlag.pl
new file mode 100755
index 0000000..7e00167
--- /dev/null
+++ b/scripts/heasoft_to_psdlag.pl
@@ -0,0 +1,74 @@
+#!/usr/local/bin/perl
+use feature say;
+use utf8;
+
+# reads (stdin) a heasoft data set (single filter) as compiled circa 2016 by Cackett @ WSU
+# outputs (stdout) a lightcurve in psdlag circa 2016 Zoghbi @ UMY input format
+
+# smallest time interval.
+my $Δt = .01;
+# normalization factor based on inspection of data
+my $flux_norm = 1e14;
+
+# This sets up the capture widths for rounding
+# Add two because need to read more digits for rounding
+$keep_digits = -1 * log10($Δt) + 2;
+$throw_digits = 5 - $keep_digits;
+
+# Skip first line -- shouldn't need this for a single filter
+# <>;
+
+#Start file with number of light curves (1) and Δt
+my $linetoprint="# 1 $Δt";
+
+my $num_avg = 0;
+while (<>) {
+    $_ =~ /^.{8}\s+([0-9]{4})\s+([0-9]{5})\.([0-9]{$keep_digits})[0-9]{$throw_digits}(.*)/;
+    # wavelength/label, characteristic, mantissa, everything else
+    ($λ, $t𝓃, $t𝜀, $vals) = ($1,$2,$3,$4);
+
+    # Round t according to Δt
+    if (!(chop($t𝜀) == 0) + chop($t𝜀) > 5) {$t𝜀++;}
+    if ($t𝜀 >= (1/$Δt)) {
+        $t𝓃 += 1;
+        $t𝜀 = 0;
+    }
+
+    # Print all values for this line
+    # say "$λ\t$t𝓃.$t𝜀$vals";
+
+    # Compile variables for this line
+    ($flux_μ,$flux_σ,$mag_μ,$mag_σ,$telescope_id) = split ' ',$vals;
+    $flux_μ *= $flux_norm;
+    $flux_σ *= $flux_norm;
+
+    # Print light curve for psdlag
+    # This mess to average measurements occuring at the same time coordinate
+    $this_t = "$t𝓃.$t𝜀";
+    if ($this_t == $last_t) {
+        if ( $linetoprint =~ /([0-9\.]+)\s+([0-9e\-\.]+)\s+([0-9e\-\.]+)/ ) {
+            $new_flux_μ = (flux_μ + $2)/2;
+            # Could just take the max error
+            # $new_flux_σ = max($flux_σ,$3);
+            # Average error seems fair
+            $new_flux_σ = (flux_σ + $3)/2;
+            $linetoprint="$this_t  $new_flux_μ  $new_flux_σ";
+            $num_avg++;
+        }
+        else { die "Malformed data"; }
+    }
+    else {
+        say $linetoprint;
+        $linetoprint = "$t𝓃.$t𝜀  $flux_μ  $flux_σ";
+    }
+    $last_t = $this_t;
+}
+say $linetoprint;
+
+
+sub max ($$) { $_[$_[0] < $_[1]] }
+sub min ($$) { $_[$_[0] > $_[1]] }
+sub log10 {
+  my $n = shift;
+  return log($n)/log(10);
+}
\ No newline at end of file