Science/psdlag-agn
1
0
Fork 0
mirror of https://asciireactor.com/otho/psdlag-agn.git synced 2024-11-22 09:25:06 +00:00
Code Issues Packages Projects Releases Wiki Activity

fixed bug in averaging

Browse Source
This commit is contained in:
caes 2016-06-30 03:04:21 -04:00
parent ed7a5d8920
commit 2d9d4ebf44
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
scripts/heasoft_to_psdlag.pl
View File

@ -16,7 +16,7 @@ $keep_digits = -1 * log10($Δt) + 2;
$throw_digits = 5 - $keep_digits; $throw_digits = 5 - $keep_digits;
# Skip first line -- shouldn't need this for a single filter # Skip first line -- shouldn't need this for a single filter
# <>; <>;
#Start file with number of light curves (1) and Δt #Start file with number of light curves (1) and Δt
my $linetoprint="# 1 $Δt"; my $linetoprint="# 1 $Δt";
@ -47,11 +47,11 @@ while (<>) {
$this_t = "$t𝓃.$t𝜀"; $this_t = "$t𝓃.$t𝜀";
if ($this_t == $last_t) { if ($this_t == $last_t) {
if ( $linetoprint =~ /([0-9\.]+)\s+([0-9e\-\.]+)\s+([0-9e\-\.]+)/ ) { if ( $linetoprint =~ /([0-9\.]+)\s+([0-9e\-\.]+)\s+([0-9e\-\.]+)/ ) {
$new_flux_μ = (flux_μ + $2)/2; $new_flux_μ = ($flux_μ + $2)/2;
# Could just take the max error # Could just take the max error
# $new_flux_σ = max($flux_σ,$3); # $new_flux_σ = max($flux_σ,$3);
# Average error seems fair # Average error seems fair
$new_flux_σ = (flux_σ + $3)/2; $new_flux_σ = ($flux_σ + $3)/2;
$linetoprint="$this_t $new_flux_μ $new_flux_σ"; $linetoprint="$this_t $new_flux_μ $new_flux_σ";
$num_avg++; $num_avg++;
} }