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.
Thursday, April 15, 2010
Subscribe to:
Post Comments (Atom)
No comments:
Post a Comment