sn_lookup now using tns data.

CheerfulUser · Apr 23, 2024 · d36b14b · d36b14b
1 parent baefe2f
commit d36b14b
Show file tree

Hide file tree

Showing 3 changed files with 663 additions and 18 deletions.
diff --git a/tessreduce/catalog_tools.py b/tessreduce/catalog_tools.py
@@ -31,10 +31,8 @@ def Get_Catalogue(tpf, Catalog = 'gaia'):
 		Gmag 	array 	Gmags of sources
 	"""
 	c1 = SkyCoord(tpf.ra, tpf.dec, frame='icrs', unit='deg')
-	# Use pixel scale for query size
-	pix_scale = 4.0  # arcseconds / pixel for Kepler, default
-	if tpf.mission == 'TESS':
-		pix_scale = 21.0
+	# Use pixel scale for query size	
+	pix_scale = 21.0
 	# We are querying with a diameter as the radius, overfilling by 2x.
 	from astroquery.vizier import Vizier
 	Vizier.ROW_LIMIT = -1
@@ -50,10 +48,10 @@ def Get_Catalogue(tpf, Catalog = 'gaia'):
 		raise ValueError("{} not recognised as a catalog. Available options: 'gaia', 'dist','ps1'")
 	if Catalog == 'gaia':
 		result = Vizier.query_region(c1, catalog=[catalog],
-                             		 radius=Angle(1890, "arcsec"),column_filters={'Gmag':'<19'})
+                             		 radius=Angle(np.max(tpf.shape[1:]) * pix_scale + 60, "arcsec"),column_filters={'Gmag':'<19'})
 	else:
 		result = Vizier.query_region(c1, catalog=[catalog],
-									 radius=Angle(np.max(tpf.shape[1:]) * pix_scale, "arcsec"))
+									 radius=Angle(np.max(tpf.shape[1:]) * pix_scale + 60, "arcsec"))
 
 	no_targets_found_message = ValueError('Either no sources were found in the query region '
 										  'or Vizier is unavailable')