fix for i4Ds#432

Author: nicHoch

stixcore/products/level0/quicklookL0.py

@@ -656,6 +656,8 @@ class EnergyCalibration(QLProduct):
     In level 0 format.
     """
 
+    PRODUCT_PROCESSING_VERSION = 3
+
     def __init__(
         self, *, service_type, service_subtype, ssid, control, data, idb_versions=defaultdict(SCETimeRange), **kwargs
     ):
@@ -680,6 +682,7 @@ def from_levelb(cls, levelb, parent=""):
         control.add_basic(name="integration_time", nix="NIX00122", packets=packets, dtype=np.uint32, attr="value")
         control.add_basic(name="quiet_time", nix="NIX00123", packets=packets, dtype=np.uint16, attr="value")
         control.add_basic(name="live_time", nix="NIX00124", packets=packets, dtype=np.uint32, attr="value")
+
         control.add_basic(name="average_temperature", nix="NIX00125", packets=packets, dtype=np.uint16, attr="value")
         control.add_data("detector_mask", _get_detector_mask(packets))
         control.add_data("pixel_mask", _get_pixel_mask(packets))
@@ -720,7 +723,8 @@ def from_levelb(cls, levelb, parent=""):
             channels.append(list(chain(*[ch.tolist() for ch in sub_channels])))
         control["num_channels"] = [len(c) for c in channels]
 
-        duration = SCETimeDelta(packets.get_value("NIX00122").astype(np.uint32))
+        # fix for wrong calibration in IDB https://github.com/i4Ds/STIXCore/issues/432
+        duration = SCETimeDelta(packets.get_value("NIX00122") * 10)
         time = SCETime(control["scet_coarse"], control["scet_fine"]) + duration / 2
 
         dids = packets.get_value("NIXD0155")

stixcore/products/level1/quicklookL1.py

@@ -6,6 +6,8 @@
 
 import numpy as np
 
+from stixcore.idb.manager import IDBManager
+from stixcore.processing import engineering
 from stixcore.products.level0.quicklookL0 import QLProduct
 from stixcore.products.product import L1Mixin
 from stixcore.time import SCETimeRange
@@ -191,6 +193,8 @@ class EnergyCalibration(QLProduct, L1Mixin):
     In level 1 format.
     """
 
+    PRODUCT_PROCESSING_VERSION = 3
+
     def __init__(
         self, *, service_type, service_subtype, ssid, control, data, idb_versions=defaultdict(SCETimeRange), **kwargs
     ):
@@ -212,6 +216,32 @@ def __init__(
     def is_datasource_for(cls, *, service_type, service_subtype, ssid, **kwargs):
         return kwargs["level"] == "L1" and service_type == 21 and service_subtype == 6 and ssid == 41
 
+    @classmethod
+    def from_level0(cls, l0product, parent=""):
+        l1 = cls(
+            service_type=l0product.service_type,
+            service_subtype=l0product.service_subtype,
+            ssid=l0product.ssid,
+            control=l0product.control,
+            data=l0product.data,
+            idb_versions=l0product.idb_versions,
+            comment=l0product.comment,
+            history=l0product.history,
+        )
+
+        l1.control.replace_column("parent", [parent] * len(l1.control))
+        l1.level = "L1"
+        engineering.raw_to_engineering_product(l1, IDBManager.instance)
+
+        # fix for wrong calibration in IDB https://github.com/i4Ds/STIXCore/issues/432
+        # nix00122 was wrong assumed to be in ds but it is plain s
+        l1.control["integration_time"] = l1.control["integration_time"] * 10
+        # nix00124 was wrong assumed to be in ds but it is unscaled ms
+        l1.control["live_time"] = (l1.control["live_time"] / 100.0).to("ms").astype(np.uint32)
+        # nix00124 was wrong assumed to be in s but it is us
+        l1.control["quiet_time"] = (l1.control["quiet_time"] / 100000.0).to("us")
+        return l1
+
 
 class TMStatusFlareList(QLProduct, L1Mixin):
     """Quick Look TM Management status and Flare list data product.

stixcore/products/tests/test_quicklook.py

@@ -48,7 +48,7 @@
         qll1.EnergyCalibration,
         "energy",
         "0659318520f00000",
-        "0659326919f00000",
+        "0659402519f00000",
         1,
     ),
 ]