1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
|
#!/usr/bin/env python3
#
# Copyright (c) 2017 Weitian LI <weitian@aaronly.me>
# MIT license
#
"""
Reprocess Chandra ACIS raw "secondary" (level=1) data using the
CIAO contrib tool ``chandra_repro``, and build the ``manifest.yaml``.
"""
import os
import argparse
import subprocess
import logging
from glob import glob
from _context import acispy
from acispy.manifest import get_manifest
from acispy.pfiles import setup_pfiles
import update_manifest
logging.basicConfig(level=logging.INFO)
logger = logging.getLogger(__name__)
def repro_acis(indir=".", outdir="repro", clobber=False):
"""
Reprocess Chandra ACIS raw level=1 data using ``chandra_repro``.
"""
logger.info("Reprocess Chandra ACIS raw level=1 data ...")
clobber = "yes" if clobber else "no"
subprocess.check_call(["punlearn", "ardlib"])
subprocess.check_call(["punlearn", "chandra_repro"])
subprocess.check_call([
"chandra_repro", "indir=%s" % indir, "outdir=%s" % outdir,
"verbose=2", "clobber=%s" % clobber
])
logger.info("Fix asol.lis by striping the absolute directory ...")
asol_lis = glob("%s/acisf*_asol1.lis" % outdir)[0]
asol = [os.path.basename(fp) for fp in open(asol_lis).readlines()]
os.rename(asol_lis, asol_lis+".orig")
open(asol_lis, "w").writelines(asol)
def build_manifest(reprodir="repro", manifestfile="manifest.yaml"):
"""
Build the ``manifest.yaml`` with reprocessed products.
"""
logger.info("Build manifest with reprocessed products ...")
if os.path.exists(manifestfile):
logger.info("Use existing manifest file: %s" % manifestfile)
else:
logger.info("Create a new manifest file: %s" % manifestfile)
open(manifestfile, "w").close()
manifest = get_manifest(manifestfile)
update_manifest.add_directory(manifest)
update_manifest.add_repro(reprodir=reprodir, manifest=manifest)
def main():
parser = argparse.ArgumentParser(
description="Reprocess ACIS raw data & create manifest.yaml")
parser.add_argument("-C", "--clobber",
dest="clobber", action="store_true",
help="overwrite existing file")
args = parser.parse_args()
setup_pfiles(["chandra_repro", "ardlib"])
repro_acis(clobber=args.clobber)
build_manifest()
if __name__ == "__main__":
main()
|
