SWIFT_BAT_GRB_POS_ACK conversion

gcn_classic_to_json/notices/SWIFT_BAT_GRB_POS_ACK/__init__.py

@@ -1,10 +1,79 @@
+import numpy as np
+
 from ... import utils
 
 
 def parse(bin):
+    bin[15]  # Unused. According to docs: '4 bytes for the future'
+    bin[19]  # Unused. Flags are either internal or equivalent to bin[18]
+    bin[26:36]  # Unused. According to docs: '40 bytes for the future'
+    bin[36]  # Unused. Flags Equivalent to bin[18]
+    bin[38]  # Unused. Sun/Moon parameters
+
+    integ_time = bin[14] * 16 / 1000
+
+    lat, lon = bin[16:17].view(">i2")
+
+    soln_status_bits = np.flip(np.unpackbits(bin[18:19].view(dtype="u1")))
+
+    flag_descriptions = {
+        0: "A point source was found.\n",
+        1: "It is a GRB.\n",
+        2: "It is an interesting src.\n",
+        3: "It is in the flight catalog.\n",
+        5: "It is definitely not a GRB.\n",
+        6: "It is probably not a GRB or Transient(hi bkg level).\n",
+        7: "It is probably not a GRB or Transient(low image significance; < 7).\n",
+        8: "It is in the ground catalog.\n",
+        9: "It is probably not a GRB or Transient(negative bkg slop).\n",
+        10: "StraTracker not locked so trigger porbably bogus.\n",
+        11: "It is probably not a GRB or Transient(very low image significance; < 6.5).",
+        12: "It is the catalog of sources to be blocked.\n",
+        13: "There is a bright star nearby.\n",
+        14: "This was orginally a SubTresh, but it is now converted to a real BAT_POS.\n",
+        15: "This is a source that has purposefully been removed from on-board catalog.\n",
+        16: "This matched a Nearby_Galaxy in the on-board catalog.\n",
+        28: "There was a temporal coincidence with another event.\n",
+        29: "There was a spatial coincidence with another event.\n",
+        30: "This is a test submission",
+    }
+    comments = "".join(
+        [
+            val
+            for (key, val) in flag_descriptions.items()
+            if (soln_status_bits[key] == 1)
+        ]
+    )
+
+    energy_ranges = [[15, 25], [15, 50], [25, 100], [50, 350]]
+    energy_range_idx = np.flip(bin[37:38].view(dtype="i1"))[0]
+    energy_range = energy_ranges[energy_range_idx]
+
     return {
+        "mission": "SWIFT",
+        "instrument": "BAT",
+        "id": [bin[4]],
         "trigger_time": utils.datetime_to_iso8601(bin[5], bin[6]),
+        "trigger_type": "image" if soln_status_bits[4] == 1 else "rate",
+        "image_duration": integ_time if soln_status_bits[4] == 1 else None,
+        "image_energy_range": energy_range if soln_status_bits[4] == 1 else None,
+        "rate_duration": integ_time if soln_status_bits[4] == 0 else None,
+        "rate_energy_range": integ_time if soln_status_bits[4] == 0 else None,
         "ra": 1e-4 * bin[7],
         "dec": 1e-4 * bin[8],
         "ra_dec_error": 1e-4 * bin[11],
+        "instrument_phi": 1e-2 * bin[12],
+        "instrument_theta": 1e-2 * bin[13],
+        "latitude": lat * 1e-2,
+        "longitude": lon * 1e-2,
+        "rate_snr": bin[21] * 1e-2,
+        "image_snr": bin[20] * 1e-2,
+        "n_events": bin[9],
+        "image_peak": bin[10],
+        "background_events": bin[22],
+        "background_start_time": utils.datetime_to_iso8601(bin[5], bin[23]),
+        "backgroun_duration": bin[24] * 1e-2,
+        "trigger_index": bin[17],
+        "catalog_number": bin[25] if soln_status_bits[3] == 1 else None,
+        "additional_info": comments if comments else None,
     }

gcn_classic_to_json/notices/SWIFT_BAT_GRB_POS_ACK/example.json

@@ -1,6 +1,30 @@
 {
+  "mission": "SWIFT",
+  "instrument": "BAT",
+  "id": [
+    1227767
+  ],
   "trigger_time": "2024-05-11T18:06:53.220Z",
+  "trigger_type": "rate",
+  "image_duration": null,
+  "image_energy_range": null,
+  "rate_duration": 8.192,
+  "rate_energy_range": 8.192,
   "ra": 336.66450000000003,
   "dec": 8.5135,
-  "ra_dec_error": 0.05
+  "ra_dec_error": 0.05,
+  "instrument_phi": -173.54,
+  "instrument_theta": 28.38,
+  "latitude": -11.14,
+  "longitude": 151.31,
+  "rate_snr": 14.89,
+  "image_snr": 8.52,
+  "n_events": 7538,
+  "image_peak": 281,
+  "background_events": 74677,
+  "background_start_time": "2024-05-11T18:06:22.570Z",
+  "backgroun_duration": 24.0,
+  "trigger_index": 262,
+  "catalog_number": null,
+  "additional_info": "A point source was found.\nIt is a GRB.\nThere was a spatial coincidence with another event.\n"
 }