antares-devkit 0.6.0__tar.gz → 0.6.2__tar.gz

{antares_devkit-0.6.0 → antares_devkit-0.6.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: antares-devkit
-Version: 0.6.0
+Version: 0.6.2
 Summary: Provides tools and utility methods to write and debug ANTARES filters with Python.
 Keywords: antares,devkit,filter development
 Author-Email: NSF NOIRLab ANTARES Team <antares@noirlab.edu>

{antares_devkit-0.6.0 → antares_devkit-0.6.2}/antares_devkit/filters/rubin_NovTest_soraisam/__init__.py

@@ -70,8 +70,17 @@ class rubin_NovTest_soraisam(BaseFilter):
         import os
 
         antares_environment = os.getenv("ANTARES_ENVIRONMENT_URL", "")
+
+        ## check locus aggregation
+        C0 = False
+        if 'lsst' in locus.properties['survey'].keys():
+            if len(locus.properties['survey']['lsst']['dia_object_id']) > 1:
+                C0 = True 
+        if C0:
+            return        
+        
         
-        C1 = (locus.alert.properties['lsst_diaSource_reliability'] < 0.9) or (locus.alert.properties['lsst_diaSource_snr'] < 10.0)
+        C1 = (locus.alert.properties['lsst_diaSource_reliability'] < 0.9)
         if C1:
             return
 
@@ -95,8 +104,8 @@ class rubin_NovTest_soraisam(BaseFilter):
         TS_all = locus.timeseries('all')
         TS_lsst = locus.timeseries('lsst')
 
-        ## want more than two measurements across all filters and time span less than 15 days (young)
-        C4 = (len(TS_all) > 2) and (max(TS_lsst["ant_mjd"].data) - min(TS_lsst["ant_mjd"].data) < 15.0)
+        ## want more than two measurements across all filters 
+        C4 = (len(TS_lsst) > 2) and (max(TS_lsst["ant_mjd"].data) - min(TS_lsst["ant_mjd"].data) < 45.0)
         if not C4:
             return
 
@@ -184,6 +193,6 @@ class rubin_NovTest_soraisam(BaseFilter):
 
         locus.tag("SN_candies")
         
-        msg = f"Young SN candidate: triggering alert={alert_id} for locus <{antares_environment}/loci/{ant_id}|{ant_id}> {known_to_tns},  RA={ra_str} DEC={dec_str},  {current_pb}_mag={current_mag:.2f},  #dets in LSST={len(TS_lsst)},  UTC={current.iso}"
-        if current_mag < 22.5:
-            self.send_message_to_slack(msg)
+        msg = f"SN candidate: triggering alert={alert_id} for locus <{antares_environment}/loci/{ant_id}|{ant_id}> {known_to_tns},  RA={ra_str} DEC={dec_str},  {current_pb}_mag={current_mag:.2f},  #dets in LSST={len(TS_lsst)},  UTC={current.iso}"
+        
+        self.send_message_to_slack(msg)

{antares_devkit-0.6.0 → antares_devkit-0.6.2}/antares_devkit/filters/young_rubin_transients_soraisam/__init__.py

@@ -1,4 +1,4 @@
-## v1
+## v2
 from antares_devkit.models import BaseFilter
 
 
@@ -70,7 +70,17 @@ class young_rubin_transients_soraisam(BaseFilter):
         import os
 
         antares_environment = os.getenv("ANTARES_ENVIRONMENT_URL", "")
-        C1 = (locus.alert.properties['lsst_diaSource_reliability'] < 0.9) or (locus.alert.properties['lsst_diaSource_snr'] < 10.0)
+
+        ## check locus aggregation
+        C0 = False
+        if 'lsst' in locus.properties['survey'].keys():
+            if len(locus.properties['survey']['lsst']['dia_object_id']) > 1:
+                C0 = True 
+        if C0:
+            return
+        
+        
+        C1 = (locus.alert.properties['lsst_diaSource_reliability'] < 0.9)
         if C1:
             return
 
@@ -106,7 +116,7 @@ class young_rubin_transients_soraisam(BaseFilter):
         TS_all = locus.timeseries('all')
         TS_lsst = locus.timeseries('lsst')
 
-        C4 = (len(TS_all) > 2) and (max(TS_lsst["ant_mjd"].data) - min(TS_lsst["ant_mjd"].data) < 15.0)
+        C4 = (len(TS_lsst) > 2) and (max(TS_lsst["ant_mjd"].data) - min(TS_lsst["ant_mjd"].data) < 45.0)
         if not C4:
             return
 
@@ -126,6 +136,20 @@ class young_rubin_transients_soraisam(BaseFilter):
         C7 = np.all(TS_lsst['lsst_diaSource_isNegative'].data==False)
         if not C7:
             return
+
+        ## check for MO, if all detections within 2 hr            
+        if max(TS_lsst["ant_mjd"].data) - min(TS_lsst["ant_mjd"].data) < 2.0/24.0:
+            C8 = False
+            for PB in np.unique(TS_lsst["ant_passband"].data):
+                y = TS_lsst["ant_mag"].data[TS_lsst["ant_passband"].data==PB]
+                yerr = TS_lsst["ant_magerr"].data[TS_lsst["ant_passband"].data==PB]
+                if np.std(y) > 0.5:
+                    C8 = True
+                    break
+        else:
+            C8 = True
+        if not C8:
+            return
         
         
         ant_id = locus.locus_id
@@ -143,5 +167,4 @@ class young_rubin_transients_soraisam(BaseFilter):
         locus.tag("young_rubin_transients")
         
         msg = f"Young extragalactic transient: triggering alert={alert_id} for locus <{antares_environment}/loci/{ant_id}|{ant_id}> {known_to_tns},  RA={ra_str} DEC={dec_str},  {current_pb}_mag={current_mag:.2f},  #dets in LSST={len(TS_lsst)},  UTC={current.iso}"
-        if current_mag<22.5:
-            self.send_message_to_slack(msg)
+        self.send_message_to_slack(msg)

{antares_devkit-0.6.0 → antares_devkit-0.6.2}/pyproject.toml

@@ -6,7 +6,7 @@ build-backend = "pdm.backend"
 
 [project]
 name = "antares-devkit"
-version = "0.6.0"
+version = "0.6.2"
 description = "Provides tools and utility methods to write and debug ANTARES filters with Python."
 readme = "README.md"
 requires-python = ">=3.9"