rlwastro / mastcasjobs Goto Github PK

import mastcasjobs
import getpass
import os

if not os.environ.get('CASJOBS_USERID'):
    os.environ['CASJOBS_USERID'] = input('Enter Casjobs UserID:')
if not os.environ.get('CASJOBS_PW'):
    os.environ['CASJOBS_PW'] = getpass.getpass('Enter Casjobs password:')

query = """select o.objID, o.raMean, o.decMean,
o.nDetections, o.ng, o.nr, o.ni, o.nz, o.ny,
m.gMeanPSFMag, m.rMeanPSFMag, m.iMeanPSFMag, m.zMeanPSFMag, m.yMeanPSFMag
from fGetNearbyObjEq(187.706,12.391,50.0/60.0) nb
inner join ObjectThin o on o.objid=nb.objid and o.nDetections>1
inner join MeanObject m on o.objid=m.objid and o.uniquePspsOBid=m.uniquePspsOBid
"""

jobs = mastcasjobs.MastCasJobs(context="PanSTARRS_DR2")
results = jobs.quick(query, task_name="python cone search")
print(results)

---------------------------------------------------------------------------
ExpatError                                Traceback (most recent call last)
Input In [15], in <cell line: 11>()
      1 query = """select o.objID, o.raMean, o.decMean,
      2 o.nDetections, o.ng, o.nr, o.ni, o.nz, o.ny,
      3 m.gMeanPSFMag, m.rMeanPSFMag, m.iMeanPSFMag, m.zMeanPSFMag, m.yMeanPSFMag
   (...)
      6 inner join MeanObject m on o.objid=m.objid and o.uniquePspsOBid=m.uniquePspsOBid
      7 """
      9 jobs = mastcasjobs.MastCasJobs(context="PanSTARRS_DR2")
---> 11 results = jobs.quick(query, task_name="python cone search")
     12 print(results)

File ~\anaconda3\envs\stsci\lib\site-packages\mastcasjobs\__init__.py:152, in MastCasJobs.quick(self, q, context, task_name, system, astropy)
    131 def quick(self, q, context=None, task_name="quickie", system=False, astropy=True):
    132     """
    133     Run a quick job. Like CasJobs method but adds astropy option.
    134 
   (...)
    150 
    151     """
--> 152     results = super(MastCasJobs,self).quick(q, context=context, task_name=task_name, system=system)
    153     if astropy:
    154         return MastCasJobs.convert_quick_table(results)

File ~\anaconda3\envs\stsci\lib\site-packages\casjobs.py:176, in CasJobs.quick(self, q, context, task_name, system)
    173 params = {"qry": q, "context": context, "taskname": task_name,
    174         "isSystem": system}
    175 r = self._send_request("ExecuteQuickJob", params=params)
--> 176 return self._parse_single(r.text, "string")

File ~\anaconda3\envs\stsci\lib\site-packages\casjobs.py:148, in CasJobs._parse_single(self, text, tagname)
    134 def _parse_single(self, text, tagname):
    135     """
    136     A hack to get the content of the XML responses from the CAS server.
    137 
   (...)
    146 
    147     """
--> 148     return minidom.parseString(text)\
    149             .getElementsByTagName(tagname)[0].firstChild.data

File ~\anaconda3\envs\stsci\lib\xml\dom\minidom.py:1998, in parseString(string, parser)
   1996 if parser is None:
   1997     from xml.dom import expatbuilder
-> 1998     return expatbuilder.parseString(string)
   1999 else:
   2000     from xml.dom import pulldom

File ~\anaconda3\envs\stsci\lib\xml\dom\expatbuilder.py:925, in parseString(string, namespaces)
    923 else:
    924     builder = ExpatBuilder()
--> 925 return builder.parseString(string)

File ~\anaconda3\envs\stsci\lib\xml\dom\expatbuilder.py:223, in ExpatBuilder.parseString(self, string)
    221 parser = self.getParser()
    222 try:
--> 223     parser.Parse(string, True)
    224     self._setup_subset(string)
    225 except ParseEscape:

ExpatError: not well-formed (invalid token): line 14, column 132

---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
Input In [11], in <cell line: 4>()
      5     print("Retrieving table MyDB.{} (if it exists)".format(DBtable))
----> 6     tab = jobs.fast_table(DBtable, verbose=True)
      7 except ValueError:

File ~\anaconda3\envs\stsci\lib\site-packages\mastcasjobs\__init__.py:208, in MastCasJobs.fast_table(self, table, verbose)
    206 except Exception as e:
    207     # raise ValueError("table MyDB.{} not found".format(table)) from None
--> 208     raise_from(ValueError("table MyDB.{} not found".format(table)), None)
    209 # get table from the quick_casjobs.cgi service

File <string>:3, in raise_from(value, from_value)

ValueError: table MyDB.SWEEPS not found

During handling of the above exception, another exception occurred:

ExpatError                                Traceback (most recent call last)
Input In [11], in <cell line: 4>()
      8 print("Table MyDB.{} not found, running query to create it".format(DBtable))
     10 # drop table if it already exists
---> 11 jobs.drop_table_if_exists(DBtable)
     13 #get main information
     14 query = """
     15     select a.ObjID,  RA=a.raMean, Dec=a.decMean, RAerr=a.raMeanErr, Decerr=a.decMeanErr,
     16         c.NumFilters, c.NumVisits,
   (...)
     30      join AstromSumPropMagAper2Cat c on a.ObjID=c.ObjID
     31 """.format(DBtable)

File ~\anaconda3\envs\stsci\lib\site-packages\mastcasjobs\__init__.py:129, in MastCasJobs.drop_table_if_exists(self, table)
    120 def drop_table_if_exists(self, table):
    121     """
    122     Drop table from MyDB without an error if it does not exist
    123 
   (...)
    127 
    128     """   
--> 129     results = self.quick("DROP TABLE IF EXISTS {}".format(table),context="MYDB")

File ~\anaconda3\envs\stsci\lib\site-packages\mastcasjobs\__init__.py:152, in MastCasJobs.quick(self, q, context, task_name, system, astropy)
    131 def quick(self, q, context=None, task_name="quickie", system=False, astropy=True):
    132     """
    133     Run a quick job. Like CasJobs method but adds astropy option.
    134 
   (...)
    150 
    151     """
--> 152     results = super(MastCasJobs,self).quick(q, context=context, task_name=task_name, system=system)
    153     if astropy:
    154         return MastCasJobs.convert_quick_table(results)

File ~\anaconda3\envs\stsci\lib\site-packages\casjobs.py:176, in CasJobs.quick(self, q, context, task_name, system)
    173 params = {"qry": q, "context": context, "taskname": task_name,
    174         "isSystem": system}
    175 r = self._send_request("ExecuteQuickJob", params=params)
--> 176 return self._parse_single(r.text, "string")

File ~\anaconda3\envs\stsci\lib\site-packages\casjobs.py:148, in CasJobs._parse_single(self, text, tagname)
    134 def _parse_single(self, text, tagname):
    135     """
    136     A hack to get the content of the XML responses from the CAS server.
    137 
   (...)
    146 
    147     """
--> 148     return minidom.parseString(text)\
    149             .getElementsByTagName(tagname)[0].firstChild.data

File ~\anaconda3\envs\stsci\lib\xml\dom\minidom.py:1998, in parseString(string, parser)
   1996 if parser is None:
   1997     from xml.dom import expatbuilder
-> 1998     return expatbuilder.parseString(string)
   1999 else:
   2000     from xml.dom import pulldom

File ~\anaconda3\envs\stsci\lib\xml\dom\expatbuilder.py:925, in parseString(string, namespaces)
    923 else:
    924     builder = ExpatBuilder()
--> 925 return builder.parseString(string)

File ~\anaconda3\envs\stsci\lib\xml\dom\expatbuilder.py:223, in ExpatBuilder.parseString(self, string)
    221 parser = self.getParser()
    222 try:
--> 223     parser.Parse(string, True)
    224     self._setup_subset(string)
    225 except ParseEscape:

ExpatError: not well-formed (invalid token): line 14, column 132