def registerOptions(parser):
parser.add_option("--wf_type",dest="wftype",type="string",
help="window func type [none,radecbcut,healpix, depthvar]",
default="depthvar")
WindowDecBcut.registerOptions(parser)
WindowHealpix.registerOptions(parser)
WindowDepthVariations.registerOptions(parser)
def registerOptions(parser):
parser.add_option(
"--wf_type",
dest="wftype",
type="string",
help="window func type [none,radecbcut,healpix, depthvar]",
default="depthvar")
WindowDecBcut.registerOptions(parser)
WindowHealpix.registerOptions(parser)
WindowDepthVariations.registerOptions(parser)
def getWindowFunc(o):
if o.wftype == "none":
return WindowBase()
elif o.wftype == "radecbcut":
return WindowDecBcut(options=o)
elif o.wftype == "healpix":
return WindowHealpix(options=o)
elif o.wftype == "depthvar":
return WindowDepthVariations(options=o)
else:
print "Bad WF type:", o.wftype
stop()
def readWindowH5(dataset):
## note that humna map is just a wrapper and gets saved/realoaded from healpix map
name=dataset.attrs['type']
## loop over possible types
## (why can't i do 'for t in [WindowBase, WindowDecBcut]:'?)
toret=WindowBase.readH5(dataset)
if toret is not None: return toret
toret=WindowDecBcut.readH5(dataset)
if toret is not None: return toret
toret=WindowHealpix.readH5(dataset)
if toret is not None: return toret
print "Unknown window type!"
stop()
def readWindowH5(dataset):
## note that humna map is just a wrapper and gets saved/realoaded from healpix map
name = dataset.attrs['type']
## loop over possible types
## (why can't i do 'for t in [WindowBase, WindowDecBcut]:'?)
toret = WindowBase.readH5(dataset)
if toret is not None: return toret
toret = WindowDecBcut.readH5(dataset)
if toret is not None: return toret
toret = WindowHealpix.readH5(dataset)
if toret is not None: return toret
print "Unknown window type!"
stop()
