GECAM_FLT and GECAM_GND conversions

gcn_classic_to_json/notices/GECAM_FLT/__init__.py

@@ -0,0 +1,22 @@
+import numpy as np
+
+from ..GECAM_GND import parse as parse_gecam_gnd
+
+
+def parse(bin):
+    bin[15]  # Unused. According to docs: '4 bytes for the future'
+    bin[27:39]  # Unused. According to docs: '32 bytes for the future'
+    bin[21]  # Temporarily Unused. This is trigger_duration; currently no keyword
+    bin[22:24]  # Unused. Equivalent to rate_energy_range
+
+    detector_options = ["on", "triggered"]
+    detector_bits = np.flip(np.unpackbits(bin[26:27].view(dtype="u1")))
+    detector_status = [detector_options[bit] for bit in detector_bits[:25]]
+    detectors = dict(zip(np.arange(25.0) + 1, detector_status))
+
+    return {
+        **parse_gecam_gnd(bin),
+        "rate_snr": bin[20] * 1e-2,
+        "rate_energy_range": [bin[24], bin[25]],
+        "detector_status": detectors,
+    }

gcn_classic_to_json/notices/GECAM_FLT/example.json

@@ -0,0 +1,58 @@
+{
+  "mission": "GECAM",
+  "id": [
+    368
+  ],
+  "messenger": "EM",
+  "mission_type": "B",
+  "trigger_time": "2024-05-17T13:01:22.550Z",
+  "trigger_type": "rate",
+  "net_count_rate": 550,
+  "rate_duration": 10.0,
+  "ra": 57.830000000000005,
+  "dec": 17.400000000000002,
+  "ra_dec_error": 1.4899,
+  "containment_probability": 0.6826894921370859,
+  "instrument_phi": 7.46,
+  "instrument_theta": 51.27,
+  "latitude": 11.09,
+  "longitude": -8.55,
+  "classification": [
+    {
+      "SFLARE": 1
+    }
+  ],
+  "additional_info": null,
+  "rate_snr": 5.09,
+  "rate_energy_range": [
+    30,
+    140
+  ],
+  "detector_status": {
+    "1.0": "triggered",
+    "2.0": "on",
+    "3.0": "on",
+    "4.0": "on",
+    "5.0": "on",
+    "6.0": "on",
+    "7.0": "on",
+    "8.0": "on",
+    "9.0": "on",
+    "10.0": "on",
+    "11.0": "triggered",
+    "12.0": "on",
+    "13.0": "on",
+    "14.0": "on",
+    "15.0": "on",
+    "16.0": "on",
+    "17.0": "triggered",
+    "18.0": "on",
+    "19.0": "on",
+    "20.0": "on",
+    "21.0": "on",
+    "22.0": "on",
+    "23.0": "on",
+    "24.0": "on",
+    "25.0": "on"
+  }
+}

gcn_classic_to_json/notices/GECAM_GND/__init__.py

@@ -0,0 +1,41 @@
+import numpy as np
+from scipy.stats import norm
+
+from ... import utils
+
+src_class_vals = ["GRB", "SFLARE", "KNOWN_SOURCE", "GENERIC_SOURCE"]
+
+
+def parse(bin):
+    bin[15]  # Unused. According to docs: '4 bytes for the future'
+    bin[20:39]  # Unused. According to docs: '76 bytes for the future'
+
+    soln_status_bits = np.flip(np.unpackbits(bin[18:19].view(dtype="u1")))
+
+    containment_prob = norm().cdf(1) - norm.cdf(-1)
+
+    return {
+        "mission": "GECAM",
+        "id": [bin[4]],
+        "messenger": "EM",
+        "mission_type": chr(bin[19] + 64),
+        # There is an error here; You are supposed to divide by 8640000 to get the correct value
+        # But it only seems to work if you divide by 864000000; I'm assuming it's some error in encoding the packets
+        "trigger_time": utils.datetime_to_iso8601(bin[5], bin[6] * 1e-2),
+        "trigger_type": "rate",
+        "net_count_rate": bin[9],
+        "rate_duration": bin[14] * 1e-4,
+        "ra": bin[7] * 1e-4,
+        "dec": bin[8] * 1e-4,
+        "ra_dec_error": bin[10] * 1e-4,
+        "containment_probability": containment_prob,
+        "instrument_phi": bin[12] * 1e-2,
+        "instrument_theta": bin[13] * 1e-2,
+        "latitude": bin[16] * 1e-2,
+        "longitude": bin[17] * 1e-2,
+        "classification": ({src_class_vals[bin[11] - 1]: 1},),
+        # "classification": ({str(bin[11]) : 1}, ),
+        "additional_info": (
+            "This is a test notice." if soln_status_bits[0] == 1 else None
+        ),
+    }

gcn_classic_to_json/notices/GECAM_GND/example.json

@@ -0,0 +1,26 @@
+{
+  "mission": "GECAM",
+  "id": [
+    368
+  ],
+  "messenger": "EM",
+  "mission_type": "B",
+  "trigger_time": "2024-05-17T13:01:22.550Z",
+  "trigger_type": "rate",
+  "net_count_rate": 1306,
+  "rate_duration": 27.0,
+  "ra": 56.25,
+  "dec": 16.21,
+  "ra_dec_error": 1.57,
+  "containment_probability": 0.6826894921370859,
+  "instrument_phi": 5.98,
+  "instrument_theta": 52.83,
+  "latitude": 11.09,
+  "longitude": -8.55,
+  "classification": [
+    {
+      "SFLARE": 1
+    }
+  ],
+  "additional_info": null
+}