aboutsummaryrefslogtreecommitdiffstats
path: root/scripts/acis.py
blob: 91d9fcf0a50e53f99e1480f46dfed1501ca5778a (plain)
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
80
81
82
83
84
85
86
87
88
89
90
91
# Copyright (c) 2017 Weitian LI <liweitianux@live.com>
# MIT license

"""
Chandra ACIS utilities
"""

import math
import subprocess
import re


class ACIS:
    """
    Chandra ACIS detector properties and utilities.

    References
    ----------
    [1] CIAO - dictionary - ACIS (advanced camera for imaging and spectroscopy)
        http://cxc.harvard.edu/ciao/dictionary/acis.html
    [2] CIAO - Dictionary - PI (pulse invariant)
        http://cxc.harvard.edu/ciao/dictionary/pi.html
    """
    # Pixel size
    pixel2arcsec = 0.492  # [arcsec]
    # Number of channels
    nchannel = 1024
    # Channel energy width
    echannel = 14.6  # [eV]

    @classmethod
    def energy2channel(self, energy):
        """
        Convert energy [eV] to channel number.
        """
        return math.floor(energy/self.echannel + 1)

    @classmethod
    def get_type(self, filepath):
        """
        Determine the Chandra ACIS type (``I`` or ``S``) according the
        active ACIS chips.

        Parameters
        ----------
        filepath : str
            Path to the input FITS file

        Returns
        -------
        acis_type : str
            ``I`` if ACIS-I, ``S`` if ACIS-S;
            otherwise, ``ValueError`` raised.
        """
        subprocess.check_call(["punlearn", "dmkeypar"])
        detnam = subprocess.check_output([
            "dmkeypar", "infile=%s" % filepath, "keyword=DETNAM", "echo=yes"
        ]).decode("utf-8").strip()
        if re.match(r"^ACIS-0123", detnam):
            return "I"
        elif re.match(r"^ACIS-[0-6]*7", detnam):
            return "S"
        else:
            raise ValueError("unknown chip combination: %s" % detnam)

    @classmethod
    def get_chips_str(self, filepath, sep=":"):
        """
        Return a string of the chips of interest according to the
        active ACIS type.

        Parameters
        ----------
        filepath : str
            Path to the input FITS file
        sep : str, optional
            Separator to join the chip ranges, e.g., 0:3, 0-3

        Returns
        -------
        chips : str
            ``0:3`` if ACIS-I, ``7`` if ACIS-S;
            otherwise, ``ValueError`` raised.
        """
        acis_type = self.get_type(filepath)
        if acis_type == "I":
            return sep.join(["0", "3"])
        elif acis_type == "S":
            return "7"
        else:
            raise ValueError("unknown ACIS type")