antares-devkit 0.5.1__tar.gz → 0.6.1__tar.gz

{antares_devkit-0.5.1 → antares_devkit-0.6.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: antares-devkit
-Version: 0.5.1
+Version: 0.6.1
 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.1/antares_devkit/filters/lsst_ddf/__init__.py

@@ -0,0 +1,43 @@
+"""Description
+This filter finds locus located within a radius of 1.75 deg from the center of the LSST Deep Drilling fields.
+"""
+from antares_devkit.models import BaseFilter
+
+class In_LSSTDDF(BaseFilter):
+    NAME = 'in_LSSTDDF' #for this default filter, name is same as its tag
+
+    SLACK_CHANNEL = "#filter-in-lsst-ddf"
+
+    ## bare absolute minimum alert properties proving it's good to go
+    REQUIRED_ALERT_PROPERTIES = [
+        'ant_ra',
+        'ant_dec',
+        'ant_survey'
+    ]
+
+    OUTPUT_TAGS = [
+        {
+            'name': 'in_LSSTDDF',
+            'description': 'Locus is along the line of sight to the LSST Deep Drilling Fields',
+        }
+    ]
+
+    def _run(self, locus):
+        import numpy as np
+        from astropy import units as u
+        from astropy.coordinates import SkyCoord
+
+        if locus.alert.properties['ant_survey'] == 2: #be sure we don't work with upper limits
+            return
+
+        ra = locus.alert.properties['ant_ra']
+        dec = locus.alert.properties['ant_dec']
+        #ELAISS1, XMM_LSS, ECDFS, COSMOS, EDFS_a, EDFS_b
+        #From https://survey-strategy.lsst.io/baseline/ddf.html
+        LSST_DDF_ra =   np.array([9.45 , 35.57, 52.98, 150.11, 58.9, 63.6])
+        LSST_DDF_dec = np.array([-44.02,-4.82, -28.12, 2.23, -49.32,-47.6])
+        field_centers = SkyCoord(ra=LSST_DDF_ra*u.deg, dec=LSST_DDF_dec*u.deg)
+        c = SkyCoord(ra*u.deg, dec*u.deg, frame='icrs')
+        sep = c.separation(field_centers)
+        if (sep.deg < 1.75).any(): #LSST FOV = 3.5 deg diameter
+            locus.tag('in_LSSTDDF')

{antares_devkit-0.5.1 → antares_devkit-0.6.1}/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)
         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.5.1 → antares_devkit-0.6.1}/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.5.1 → antares_devkit-0.6.1}/antares_devkit/utils.py

@@ -23,6 +23,7 @@ def get_filters_path():
         "HighSNR": "antares_devkit.filters.highsnr",
         "IF_anomaly_detection": "antares_devkit.filters.if_anomaly_detection",  # ["numpy", "pandas", "sklearn", "astropy"]
         "IMBH": "antares_devkit.filters.imbh",
+        "In_LSSTDDF": "antares_devkit.filters.lsst_ddf", #["astropy", "numpy"]
         "InM31": "antares_devkit.filters.inm31",
         "LAISS_RFC_AD_filter": "antares_devkit.filters.laiss_rfc_ad_filter",  # ["astropy", "timeout_decorator", "alerce", "numpy", "pandas", "matplotlib", "astro_ghost", "sklearn", "scipy"]
         "AnomalyTransientDmDt": "antares_devkit.filters.malanchev_anomaly_transient_dmdt",  # ["numpy", "astropy", "light_curve"]

{antares_devkit-0.5.1 → antares_devkit-0.6.1}/pyproject.toml

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