This is the standard primary on stripe 82:
http://data.sdss3.org/sas/
These are alternative runs on stripe 82:
"good seeing"
http://data.sdss3.org/sas/
"median seeing"
http://data.sdss3.org/sas/
"bad seeing"
http://data.sdss3.org/sas/
And it matched them to the variable catalogs (co-adds) with the following code:
erinfile2 = '/home/jessica/boss/boss-qso-stripe82median.fits'
confTT = mrdfits(erinfile2, 2)
varcat0 = '/home/schlegel/varcat4/varcat-ra0-30.fits.gz'
varcat = mrdfits(varcat0, 1)
spherematch, confTT.ra, confTT.dec, varcat.ra, varcat.dec, 2./3600, i1, i2, d12
truthmatch0 = confTT[i1]
varmatch0 = varcat[i2]
varcat1 = '/home/schlegel/varcat4/varcat-ra30-60.fits.gz'
varcat = mrdfits(varcat1, 1)
spherematch, confTT.ra, confTT.dec, varcat.ra, varcat.dec, 2./3600, i1, i2, d12
truthmatch1 = confTT[i1]
varmatch1 = varcat[i2]
varcat2 = '/home/schlegel/varcat4/varcat-ra60-90.fits.gz'
varcat = mrdfits(varcat2, 1)
spherematch, confTT.ra, confTT.dec, varcat.ra, varcat.dec, 2./3600, i1, i2, d12
truthmatch2 = confTT[i1]
varmatch2 = varcat[i2]
varcat3 = '/home/schlegel/varcat4/varcat-ra90-120.fits.gz'
varcat = mrdfits(varcat3, 1)
spherematch, confTT.ra, confTT.dec, varcat.ra, varcat.dec, 2./3600, i1, i2, d12
truthmatch3 = confTT[i1]
varmatch3 = varcat[i2]
varcat4 = '/home/schlegel/varcat4/varcat-ra300-330.fits.gz'
varcat = mrdfits(varcat4, 1)
spherematch, confTT.ra, confTT.dec, varcat.ra, varcat.dec, 2./3600, i1, i2, d12
truthmatch4 = confTT[i1]
varmatch4 = varcat[i2]
varcat5 = '/home/schlegel/varcat4/varcat-ra330-360.fits.gz'
varcat = mrdfits(varcat5, 1)
spherematch, confTT.ra, confTT.dec, varcat.ra, varcat.dec, 2./3600, i1, i2, d12
truthmatch5 = confTT[i1]
varmatch5 = varcat[i2]
truthmatch = [truthmatch0, truthmatch1, truthmatch4, truthmatch5]
varmatch = [varmatch0, varmatch1, varmatch4, varmatch5]
outfile = "./singleEpochMatch.fits"
splog, 'Writing file ', outfile
mwrfits, truthmatch, outfile, /create
outfile = "./coaddedMatch.fits"
splog, 'Writing file ', outfile
mwrfits, varmatch, outfile, /create
This creates two matched catalogs, each with 735,498 objects in Stripe-82. I'm sending these files to Jo and Joe to help with the "extreme deconvolution." (Man I love that name)
No comments:
Post a Comment