Wednesday, March 3, 2010

Single Epoch vs Coadded

Joe Hennawi asked me to give him targeting files that were single epoch versus coadded for Stripe 82. Erin Sheldon had created some files for me which had different seeings:

This is the standard primary on stripe 82:
http://data.sdss3.org/sas/bosswork/groups/boss/target/2009-12-05-chunk1-retrain/bosstarget-qso-2009-12-05-chunk1-retrain-collate-fullcache.fits

These are alternative runs on stripe 82:
"good seeing"
http://data.sdss3.org/sas/bosswork/groups/boss/target/2009-11-30-stripe82pairs/bosstarget-qso-2009-11-30-stripe82pairs-collate-good-fullcache.fits

"median seeing"
http://data.sdss3.org/sas/bosswork/groups/boss/target/2009-11-30-stripe82pairs/bosstarget-qso-2009-11-30-stripe82pairs-collate-median-fullcache.fits

"bad seeing"
http://data.sdss3.org/sas/bosswork/groups/boss/target/2009-11-30-stripe82pairs/bosstarget-qso-2009-11-30-stripe82pairs-collate-bad-fullcache.fits

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