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()