mirror of
https://asciireactor.com/otho/psdlag-agn.git
synced 2024-11-22 14:15:08 +00:00
73 lines
2.7 KiB
Bash
Executable File
73 lines
2.7 KiB
Bash
Executable File
#!/usr/bin/env bash
|
|
|
|
ywindow_of () {
|
|
max=0
|
|
min=-5
|
|
while read line
|
|
# for line in $(cat $1)
|
|
do
|
|
read xval yval xerr yerr <<< $line
|
|
echo $xval $yval
|
|
if (bc <<< "$yval > $max"); then max=$yval; fi
|
|
if (bc <<< "$yval < $min"); then min=$yval; fi
|
|
# done
|
|
done < $1
|
|
echo $(bc <<< "$min - 0.5") $(bc <<< "$max + 0.5")
|
|
}
|
|
|
|
mkdir -p analyses/tables
|
|
mkdir -p analyses/plots
|
|
|
|
|
|
for analysis in analyses/*
|
|
do
|
|
# Grab and determine labels of analyses, skip if over the same band.
|
|
ref_band=$(basename $analysis|sed 's@\([^≺]*\)_≺_[^≺_]*_[^_]*_[^_]*@\1@')
|
|
echo_band=$(basename $analysis|sed 's@[^≺]*_≺_\([^≺_]*\)_[^_]*_[^_]*@\1@')
|
|
if [[ $ref_band == $echo_band ]]; then continue; fi
|
|
err_type=$(basename $analysis|sed 's@[^≺]*_≺_[^≺_]*_[^_]*_\(σ∊[CLM][MFC]\)@\1@')
|
|
|
|
# Prepare files
|
|
echo "Plotting PSD and time lags for $echo_band, referred to ${ref_band}."
|
|
echoPSD_tabfile=analyses/tables/${echo_band}_${err_type}_PSD.tab
|
|
refPSD_tabfile=analyses/tables/${ref_band}_${err_type}_PSD.tab
|
|
timelag_tabfile=analyses/tables/${ref_band}_≺_${echo_band}_${err_type}_timelag.tab
|
|
PSD_plotfile=analyses/plots/${ref_band}_≺_${echo_band}_${err_type}_PSD.png
|
|
timelag_plotfile=analyses/plots/${ref_band}_≺_${echo_band}_${err_type}_timelag.png
|
|
|
|
# Output curves to temporary files using perl script, move tables to
|
|
# permanent location. This just assumes there are no conflicts.
|
|
scripts/create_tables.pl $analysis > /dev/null
|
|
mv tmp.echoPSD $echoPSD_tabfile
|
|
mv tmp.refPSD $refPSD_tabfile
|
|
mv tmp.timelag $timelag_tabfile
|
|
|
|
|
|
# Plot PSD and save using gnuplot
|
|
# read ymin ymax <<< $(ywindow_of $timelag_tabfile)
|
|
# echo $ymin $ymax
|
|
cat scripts/templates/psd_freq.gp|
|
|
sed "s|%TITLE|Power Spectrum for Lightcurves $echo_band and $ref_band|"|
|
|
sed "s|%SUBTITLE|as reported by Fausnaugh et. al, STORM III, 2016|"|
|
|
sed "s|%FILE1|$refPSD_tabfile|"|
|
|
sed "s|%LABEL1|${ref_band} PSD|"|
|
|
sed "s|%FILE2|$echoPSD_tabfile|"|
|
|
sed "s|%LABEL2|${echo_band} PSD|"|
|
|
sed "s|%YMIN|$ymin|"|sed "s|%YMAX|$ymax|"|
|
|
sed "s|%OUTPUTFILE|$PSD_plotfile|" > tmp.gp
|
|
gnuplot tmp.gp
|
|
|
|
# Plot time lags and save using gnuplot
|
|
# read ymin ymax <<< $(ywindow_of $timelag_tabfile)
|
|
cat scripts/templates/timelag_freq.gp|
|
|
sed "s|%TITLE|Time Lag for Lightcurve $echo_band relative to $ref_band|"|
|
|
sed "s|%SUBTITLE|as reported by Fausnaugh et. al, STORM III, 2016|"|
|
|
sed "s|%FILE1|$timelag_tabfile|"|
|
|
sed "s|%LABEL1|${echo_band} Lag|"|
|
|
sed "s|%YMIN|$ymin|"|sed "s|%YMAX|$ymax|"|
|
|
sed "s|%OUTPUTFILE|$timelag_plotfile|" > tmp.gp
|
|
gnuplot tmp.gp
|
|
done
|
|
|
|
rm tmp.*
|