def stack_it(specList ): outfile = join(spec_dir, os.path.basename(specList)[:-4]+".stack") print(outfile) test_D = n.loadtxt(specList, unpack=True) print(len(test_D[0])) if len(test_D[0])>10: stack=sse.SpectraStackingEBOSS(specList, outfile ) stack.createStackMatrix()
def stack_it(specList):
outfile = join(spec_dir, os.path.basename(specList)[:-4] + ".stack")
print(outfile)
#test_D = n.loadtxt(specList, unpack=True)
#print(len(test_D[0]))
#if len(test_D[0])>10:
stack = sse.SpectraStackingEBOSS(specList,
outfile,
l_start=3.,
l_end=4.1,
csv_input=True)
stack.createStackMatrix()
stack.stackSpectra()
def stack_it(specList = join(spec_dir, "catalog_AGN.dat") ): outfile = join(spec_dir, os.path.basename(specList)[:-4]+".stack") print(specList, outfile) stack=sse.SpectraStackingEBOSS(specList, outfile, KZ_input=True, l_start=3.3, l_end=3.9 ) stack.createStackMatrix() stack.stackSpectra()
stack_it(join(spec_dir, 'catalog_lowz_liner.dat')) stack_it(join(spec_dir, 'catalog_lowz_sf.dat')) sys.exit() stack_it(join(spec_dir, 'catalog_AGN.ascii')) stack_it(join(spec_dir, 'catalog_liner.ascii')) stack_it(join(spec_dir, 'catalog_comp.ascii')) stack_it(join(spec_dir, 'catalog_sf.ascii')) spec_dir = join(os.environ['HOME'], "afs_comparat/sdss/elg/stacks/KZ/high_z/") specList = join(spec_dir, 'catalog_AGN.ascii') outfile = join(spec_dir, os.path.basename(specList)[:-6] + ".stack") print(specList, outfile) stack = sse.SpectraStackingEBOSS(specList, outfile, l_start=3.3, l_end=3.9) stack.stackSpectra() specList = join(spec_dir, 'catalog_liner.ascii') outfile = join(spec_dir, os.path.basename(specList)[:-6] + ".stack") print(specList, outfile) stack2 = sse.SpectraStackingEBOSS(specList, outfile, l_start=3.3, l_end=3.9) stack2.stackSpectra() specList = join(spec_dir, 'catalog_comp.ascii') outfile = join(spec_dir, os.path.basename(specList)[:-6] + ".stack") print(specList, outfile) stack = sse.SpectraStackingEBOSS(specList, outfile, l_start=3.3, l_end=3.9) stack.stackSpectra() specList = join(spec_dir, 'catalog_sf.ascii')
""" This script produces the stacks for emission line luminosity limited samples. """ import sys import os from os.path import join import glob import numpy as n import SpectraStackingEBOSS as sse import astropy.io.fits as fits spec_dir = join(os.environ['HOME'],"SDSS/stacks") specList = join(spec_dir, "eboss-elg_0.2_z_1.5.asc") outfile = join(spec_dir, os.path.basename(specList)[:-4]+".specMatrix") stack=sse.SpectraStackingEBOSS(specList, outfile) def getSpectra(path_to_spectrum): hdulist = fits.open(path_to_spectrum) wave = 10**hdulist[1].data['loglam'] ok=(wave>3740)&(wave<9604) flux = hdulist[1].data['flux'] hdulist.close() return wave[ok], flux[ok] for IDX_j in n.arange(0, len(stack.plates), 4096): IDX_min=IDX_j IDX_max=IDX_j+4096 IDX_str=str(IDX_min).zfill(6)+'-'+str(IDX_max).zfill(6) samp_plates, samp_mjds, samp_fiberids, samp_redshifts = stack.plates[IDX_min:IDX_max], stack.mjds[IDX_min:IDX_max], stack.fiberids[IDX_min:IDX_max], stack.redshifts[IDX_min:IDX_max]
import SpectraStackingEBOSS as sse
import time
t0 = time.time()
# create all stacks
stack_dir = join(os.environ['HOME'], "SDSS/stacks/v5")
dataList = n.array(glob.glob(join(stack_dir, "eboss-elg_*.asc")))
dataList.sort()
for specList in dataList:
print('considers', specList, time.time() - t0)
out_file = join(stack_dir, os.path.basename(specList)[:-4] + ".stack")
out_file_UV = join(stack_dir, os.path.basename(specList)[:-4] + ".UVstack")
print('starts working on', out_file, time.time() - t0)
stack = sse.SpectraStackingEBOSS(specList,
out_file,
dLambda=0.0001,
dV=-9999.99,
l_start=3.35,
l_end=3.579)
if os.path.isfile(out_file + '.specMatrix.dat') == False:
print('creates matrix', time.time() - t0)
stack.createStackMatrix()
if os.path.isfile(out_file) == False:
print('stacks', time.time() - t0)
stack.stackSpectra()
hdu = fits.open(out_file)
dd = hdu[1].data
wl = dd['wavelength']
s1 = (dd['NspectraPerPixel'] > 0.5 * n.max(dd['NspectraPerPixel']))
x, y = dd['wavelength'][s1], dd['medianStack'][s1]
yerr = dd['jackknifStackErrors'][s1]
pfit = stack.fit_UV_continuum(x, y, yerr, degree=5)
/home/comparat/software/linux/pySU/galaxy/python/SpectraStackingEBOSS.py
"""
import os
import sys
import SpectraStackingEBOSS as sse
gal_type = sys.argv[1]
zmin = float(sys.argv[2])
zmax = float(sys.argv[3])
#gal_type='qso_BL'
#zmin = 0.
#zmax = 0.5
#Nspec_max = 10000
Nspec_max = 200000
version = 'v1'
name = "sdss_" + gal_type + "_zmin_" + str(int(
10 * zmin)).zfill(2) + "_zmax_" + str(int(
10 * zmax)).zfill(2) + "_Nlt_" + str(int(Nspec_max))
out_dir = os.path.join(os.environ['OBS_REPO'], 'spectrastacks', version,
gal_type)
list_2_stack = os.path.join(out_dir, "stacklist_" + name + ".ascii")
output_file_name = os.path.join(out_dir, "spec_" + name + ".fits")
st = sse.SpectraStackingEBOSS(list_2_stack, output_file_name)
st.stackSpectra()
def stack_it(specList):
outfile = join(spec_dir, os.path.basename(specList)[:-6] + ".stack")
stack = sse.SpectraStackingEBOSS(specList, outfile, PBKT_input=True)
stack.createStackMatrix()
print(outfile)
stack.stackSpectra()
def stack_it(specList=join(spec_dir, "catalog.ascii")):
outfile = join(spec_dir, os.path.basename(specList)[:-4] + "red.stack")
stack = sse.SpectraStackingEBOSS(specList, outfile, l_end=3.9)
stack.createStackMatrix()
print(outfile)
stack.stackSpectra()