def alfredo(ra, dec):
"""Alfredo's eRosita survey extension"""
l, b = cel2gal(ra, dec)
sel = ((b > 20) & (b < 30))
sel &= ((l < 360) & (l > 240))
#sel &= ((ra > 135) & (ra < 245))
sel &= ((ra > 135) & (ra < 180))
sel &= ((dec > -30) & (dec < -10))
return sel
def ligo_mw(ra, dec):
"""The high-probability region for LIGO/MW dwarfs"""
l, b = cel2gal(ra, dec)
# Starting to the north of the Galactic plane
sel = (ra > 120) & (ra < 270)
# Upper limit in declination
sel &= (dec < -30)
# 10 deg above the Galactic plane
sel &= (b > 10.)
return sel
def footprint(ra, dec):
l, b = cel2gal(ra, dec)
angsep_lmc = angsep(constants.RA_LMC, constants.DEC_LMC, ra, dec)
angsep_smc = angsep(constants.RA_SMC, constants.DEC_SMC, ra, dec)
sel = (np.fabs(b) > 10.) \
& ((angsep_lmc < 30.) | (angsep_smc < 30.)) \
& (dec < -55.) & (ra > 100.) & (ra < 300.)
#sel = sel | ((dec < -65.) & (angsep_lmc > 5.) & (angsep_smc > 5.))
sel = sel | ((dec < -65.) & (ra > 300.) & (ra < 360.)) # SMC
sel = sel | (dec < -80.)
return sel
def footprint(ra, dec):
l, b = cel2gal(ra, dec)
# Starting to the north of the Galactic plane
sel = ((ra > 120) & (ra < 360)) | (ra < 10)
# Band in declination
sel &= (dec < -30) & (dec > -40)
# LIGO/DSPHS high probability region south of nominal stripe
sel |= ((ra > 140) & (ra < 185) & (dec < -30) & (dec > -60))
# Alfredo's eRosita survey
sel |= BlissSurvey.alfredo(ra, dec)
# 10 deg from the Galactic plane
sel &= (np.fabs(b) > 10.)
return sel
def footprintMilkyWay(ra,dec,angsep=10.):
""" Selecting fields close to the Milky Way plane.
Parameters:
-----------
ra : Right ascension (deg)
dec: Declination (deg)
angsep : Angular separation (deg)
Returns:
--------
sel : Selection of fields within the footprint
"""
glon,glat = cel2gal(ra,dec)
return (np.fabs(glat) < angsep)