1. First I make a catalog of all the BOSS and otherwise known QSOs which aren't already in the inputs for which we also have coadded fluxes by spherematching them with the varcat files. The code to do this is in the following log file ../logs/100602_2log.pro. This makes a file: ~/boss/coaddedALLQSOs.fits which has the coadded fluxes (from varcat files) and redshifts of these QSOs. These fluxes are not de-reddened. I've also added the ra/dec to this structure because Schlegel just reminded me of a function of the likelihood_compute that excludes any object in the QSO Catalog when doing the calculation of the likelihoods that has the same ra/dec as the test object. There are 6,063 QSOs in this file.
Here is a plot of these quasars
2. Combine this BOSS catalog with SDSS DR5 quasars that are included in the old version of the Monte Carlo. This is done by running Joe Hennawi's qso_photosamp.pro program (to get the SDSS DR5 QSOs), reading in coaddedBOSSQSOs file, de-reddening them both, and then merging them into one file, which is then saved as ~/boss/allbrightQSOMCInput.fits. The code to do this is in the following log file ../logs/100602_3log.pro. (I tried to do this with the SDSS DR7 and was having some problems converting from magnitudes to psffluxes)
3. Change qso_fakephoto.pro to read in the ~/boss/allbrightQSOMCInput.fits file instead of just running qso_photosamp.pro. Save this modified program as ../likelihood/qsocatalog/qso_fakephoto_jess.pro
4. Change hiz_kde_numerator.pro to call qso_fakephoto_jess.pro, instead of qso_fakephoto.pro and save as ../likelihood/qsocatalog/hiz_kde_numerator_jess.pro
5. Follow instructions from this post for how to change/create the luminosity function.
6. Enter the command :
.run hiz_kde_numerator_jess.pro into IDL to generate QSO Catalog.
7. Change likelihood_compute.pro so that it points to this QSO Catalog that you just generated as the input QSO Catalog:
; Read in the QSO
filefile2 = '/home/jessica/repository/ccpzcalib/Jessica/likelihood/qsocatalog/QSOCatalog-Mon-Jun-7-13:56:32-2010.fits'
Code to generate plots is in the log file: ../logs/100602_4log.pro
2. Combine this BOSS catalog with SDSS DR5 quasars that are included in the old version of the Monte Carlo. This is done by running Joe Hennawi's qso_photosamp.pro program (to get the SDSS DR5 QSOs), reading in coaddedBOSSQSOs file, de-reddening them both, and then merging them into one file, which is then saved as ~/boss/allbrightQSOMCInput.fits. The code to do this is in the following log file ../logs/100602_3log.pro. (I tried to do this with the SDSS DR7 and was having some problems converting from magnitudes to psffluxes)
allbrightQSOMCInput.fits QSOs
3. Change qso_fakephoto.pro to read in the ~/boss/allbrightQSOMCInput.fits file instead of just running qso_photosamp.pro. Save this modified program as ../likelihood/qsocatalog/qso_fakephoto_jess.pro
4. Change hiz_kde_numerator.pro to call qso_fakephoto_jess.pro, instead of qso_fakephoto.pro and save as ../likelihood/qsocatalog/hiz_kde_numerator_jess.pro
5. Follow instructions from this post for how to change/create the luminosity function.
6. Enter the command :
.run hiz_kde_numerator_jess.pro into IDL to generate QSO Catalog.
7. Change likelihood_compute.pro so that it points to this QSO Catalog that you just generated as the input QSO Catalog:
; Read in the QSO
filefile2 = '/home/jessica/repository/ccpzcalib/Jessica/likelihood/qsocatalog/QSOCatalog-Mon-Jun-7-13:56:32-2010.fits'
QSO Catalog (including the BOSS quasars)
Code to generate plots is in the log file: ../logs/100602_4log.pro
No comments:
Post a Comment