industry-website/htdocs/examples/run_analysis.sh

47 lines
1.7 KiB
Bash
Raw Normal View History

2022-07-11 04:07:03 +00:00
#!/usr/bin/env bash
# This script is meant to encapsulate all functions in proper order to produce
# the analyses. Parameters that this script should control: Δt, σ type, input
# dataset(lightcurve directory), more?
ref_band=$(cat ref_band)
ref_curve="lightcurves/${ref_band}.lc"
err_str=$(cat err_type)
mkdir -p analyses/tables
for echo_curve in lightcurves/*
do
# Grab and determine labels of analyses, skip if over the same band.
echo_band=$(basename $echo_curve|sed 's|\(.*\)\.lc|\1|')
if [[ $ref_band == $echo_band ]]; then continue; fi
echo "Analysing $echo_band$ref_band."
if [[ -e "analyses/${ref_band}_≺_${echo_band}/" ]]; then
echo "Results already exists. Create tables from stored results."
cp analyses/${ref_band}_≺_${echo_band}/*.out .
else
time python scripts/analyze_lightcurve.py $ref_curve $echo_curve >> ${echo_curve}.log
mkdir -p "analyses/${ref_band}_≺_${echo_band}"
cp *.out analyses/${ref_band}_≺_${echo_band}/
fi
echo "Tabling PSD and time lags referred to ${ref_band} for $echo_band{${err_str}}."
# Propagate tables into analyses/tables
echoPSD_tabfile=analyses/tables/PSD_${echo_band}_\{${err_str}\}.tab
refPSD_tabfile=analyses/tables/PSD_${ref_band}_\{${err_str}\}.tab
timelag_tabfile=analyses/tables/timelag_${ref_band}_≺_${echo_band}_\{${err_str}\}.tab
# Output curves to temporary files using perl script, move tables to
# permanent location. This just assumes there are no conflicts.
scripts/extract_tables.pl
mv tmp.echoPSD $echoPSD_tabfile
mv tmp.refPSD $refPSD_tabfile
mv tmp.timelag $timelag_tabfile
rm *.out
done
rm tmp.*