Source code for hendrics.timelags
# -*- coding: utf-8 -*-
from astropy import log
from .io import load_pds
from .io import save_as_qdp
from .base import hen_root
[docs]
def main(args=None):
import argparse
from .base import _add_default_args
description = (
"Read timelags from cross spectrum results and save them" " to a qdp file"
)
parser = argparse.ArgumentParser(description=description)
parser.add_argument("files", help="List of files", nargs="+")
_add_default_args(parser, ["loglevel", "debug"])
args = parser.parse_args(args)
if args.debug:
args.loglevel = "DEBUG"
log.setLevel(args.loglevel)
with log.log_to_file("HENlags.log"):
filelist = []
for fname in args.files:
cross = load_pds(fname)
lag, lag_err = cross.time_lag()
out = hen_root(fname) + "_lags.qdp"
save_as_qdp([cross.freq, lag], [None, lag_err], filename=out)
filelist.append(out)
return filelist