Thursday, April 15, 2010

Reductions to QSO Catalog

At the suggestion of David and Joe I am reducing the number of objects in the QSO Catalog to speed up the likelihood calculation, and then multiplying the likelihoods by the amount we reduced them. Here is the code:
___________________
file2 = "/home/jessica/boss/R06New0.5-5.0.fits"
objs2 = mrdfits(file2, 1)

minz = 0.5
dz = 0.1
nzbin = 16
reduction = 10.

.com subsample.pro
reducedlow = subsample(minz, dz, nzbin, objs2, reduction)

minz = 2.1
dz = 0.1
nzbin = 29
reduction = 2.
reducedhigh = subsample(minz, dz, nzbin, objs2, reduction)

newQSOcatalog = [reducedlow,reducedhigh]

outfile = "./reducedQSOCatalog.fits"
mwrfits, newQSOcatalog, outfile, /create
___________________

Also in the logfile: ../logs/100415log.pro with this code. The above code uses the function subsample.pro which is in the likelihood directory.

No comments:

Post a Comment