Exemplo n.º 1
0
print('Assigned science: {}'.format(len(_assigned)))

##  Load one tile, currently. 
Tile       = 74057
sv_mtl     = Table(fits.open(scratch + '/BGS/SV-ASSIGN/mtls/svmtl_{:06d}.fits'.format(Tile))[1].data)

##  Limit to only BGS.
sv_mtl     = sv_mtl[sv_mtl['SV1_BGS_TARGET'].quantity.value > 0]
print('Total targets in Tile {}: {}.'.format(Tile, len(sv_mtl)))

##  
_bitdefs   = load_mask_bits("sv1")
bgs_mask   = BitMask('sv1_bgs_mask', _bitdefs)

types      =  bgs_mask.names()
bits       = [bgs_mask.bitnum(x) for x in types]

##  Those that were assigned in this tile. 
assigned   = join(_assigned, sv_mtl, keys=['BRICKID', 'BRICK_OBJID'], join_type='inner')

##  SV Target class.
splits     = {}
title      = ''

##
for _type in ['BGS_FAINT', 'BGS_BRIGHT', 'BGS_FAINT_EXT', 'BGS_LOWQ', 'BGS_FIBMAG']:
  sample        = _type.split('_')
  sample.pop(0)
  sample        = ''.join(x for x in sample)

  splits[_type] = assigned[(assigned['SV1_BGS_TARGET'] & 2 ** bgs_mask.bitnum(_type)) != 0]
Exemplo n.º 2
0
##  Retain only the science targets; Dropping standards etc.
_assigned = _assigned[_assigned['FA_TYPE'].quantity == 1]

##  Load one tile, currently.
sv_mtl = Table(
    fits.open(scratch + '/BGS/SV-ASSIGN/mtls/svmtl_074057.fits')[1].data)

##  Limit to only BGS.
sv_mtl = sv_mtl[sv_mtl['SV1_BGS_TARGET'].quantity.value > 0]

##
_bitdefs = load_mask_bits("sv1")
bgs_mask = BitMask('sv1_bgs_mask', _bitdefs)

types = bgs_mask.names()
bits = [bgs_mask.bitnum(x) for x in types]

##
assigned = join(_assigned,
                sv_mtl,
                keys=['BRICKID', 'BRICK_OBJID'],
                join_type='inner')
assert len(_assigned) == len(assigned)

splits = {}
title = ''

for _type in [
        'BGS_FAINT', 'BGS_BRIGHT', 'BGS_FAINT_EXT', 'BGS_LOWQ', 'BGS_FIBMAG'
]:
    sample = _type.split('_')
Exemplo n.º 3
0
scratch = os.environ['CSCRATCH']

tiles = Table(
    fits.open(
        '/global/cscratch1/sd/mjwilson/BGS/SV-ASSIGN/tiles/BGS_SV_30_3x_superset60_Sep2019.fits'
    )[1].data)
utiles = np.unique(tiles['TILEID'].quantity)

##
_bitdefs = load_mask_bits("sv1")

desi_mask = BitMask('sv1_desi_mask', _bitdefs)
bgs_mask = BitMask('sv1_bgs_mask', _bitdefs)

types = bgs_mask.names()
bits = [bgs_mask.bitnum(x) for x in types]

##  L428 of https://github.com/desihub/desimodel/blob/master/py/desimodel/focalplane/geometry.py
params = load_desiparams()
fiber_dia = params['fibers']['diameter_um']

#- Platescales in um/arcsec
ps = load_platescale()

##  Add in GAMA labels.
_fits = fits.open(
    '/global/cscratch1/sd/mjwilson/BGS/SV-ASSIGN/truth/legacy/ls-GAMA-south.fits'
)
gama = Table(_fits[1].data)

for tile in utiles:
Exemplo n.º 4
0
    )[1].data)
utiles = np.unique(tiles['TILEID'].quantity)

result = Table(names=[
    'TARGETID', 'FA_TYPE', 'BRICK_OBJID', 'TARGET_RA', 'TARGET_DEC', 'BRICKID'
],
               dtype=('i8', 'i1', 'i4', 'f8', 'f8', 'i4'))

##
_bitdefs = load_mask_bits("sv1")

desi_mask = BitMask('sv1_desi_mask', _bitdefs)
bgs_mask = BitMask('sv1_bgs_mask', _bitdefs)

types = bgs_mask.names()
bits = [bgs_mask.bitnum(x) for x in types]

##  L428 of https://github.com/desihub/desimodel/blob/master/py/desimodel/focalplane/geometry.py
params = load_desiparams()
fiber_dia = params['fibers']['diameter_um']

#- Platescales in um/arcsec
ps = load_platescale()

for tile in utiles:
    print('Solving for Tile {}.'.format(tile))

    ##  Scrape from the tile picker site.
    ##  cmd = 'wget http://www.astro.utah.edu/~u6022465/SV/tiles/SV_BGS/fits_files/tile-{:06}.fits -O /global/cscratch1/sd/mjwilson/BGS/SV-ASSIGN/fiberassign/tile-{:06}.fits'.format(tile, tile)
    ##  os.system(cmd)