Now I am going to see how the likelihood compares if we use the Richards luminosity function and we also include the (0.5 < z < 5.0). The code to run this is in the following log file: ../log/100427_3log.pro. This works even better than the last test:

eps = 1e-30

bossqsolike = total(likelihood.L_QSO_Z[16:32],1) ;quasars between redshift 2.1 and 3.7

qsolcut1 = where( alog10(total(likelihood.L_QSO_Z[16:32],1)) LT -9.0)

den = total(targets.L_EVERYTHING_ARRAY[0:4],1) + total(likelihood.L_QSO_Z[0:44],1) + eps

num = bossqsolike + eps

NEWRATIO = num/den

NEWRATIO[qsolcut1] = 0 ; eliminate objects with low L_QSO value

IDL> print, n_elements(ql) ; number quasars

461

IDL> print, n_elements(ql)*1.0/(n_elements(sl)+n_elements(ql)) ;percent accuracy

0.256111

IDL> print, n_elements(ql) + n_elements(sl) ; total targeted

1800

IDL>

IDL> print, n_elements(nql) ; number quasars

575

IDL> print, n_elements(nql)*1.0/(n_elements(nsl)+n_elements(nql)) ;percent accuracy

0.319444

IDL> print, n_elements(nql) + n_elements(nsl) ; total targeted

We now go from 25.6% accuracy to 31.9% accuracy! And we find 114 more quasars!

Congratulations on the improvements!

ReplyDelete(And I love the blog, of course.)