Friday, August 6, 2010

HRH07 Likelihood Results

logfile is here:
../log/100806log.pro and here ../log/100809log.pro

The likelihood results for the Hopkins, Richards, Hernquist 2007 (HRH07) luminosity function:

The way that lratio is calculated (optimized with redshift range for best results):


;Read in computed likelihoods from threshold file
eps = 1e-30
bossqsolike = total(likelihood.L_QSO_Z[16:35],1) ;quasars between redshift 2.1 and 3.5
qsolcut = where(alog10(bossqsolike) 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[qsolcut] = 0 ; eliminate objects with low L_QSO value
likelihood.l_ratio = newratio

Targets/dec^2------v1 threshold----new threshold
-----0------------1.00000--------1.00000
-----10-----------0.749085-------0.641645
-----20-----------0.490932-------0.390440
-----30-----------0.347334-------0.273209
-----40-----------0.253020-------0.195684
-----50-----------0.198630-------0.156788
-----60-----------0.160074-------0.128583
-----70-----------0.128912-------0.105899
-----80-----------0.107913-------0.0899322
-----90-----------0.0917588------0.0786345

Targets/dec^2--# Targeted QSOs:----v1----new version
-----0------------------------1---------1
-----10----------------------443-------426.000
-----20----------------------656-------664.000
-----30----------------------781-------797.000
-----40----------------------897-------902.000
-----50----------------------965-------975.000
-----60---------------------1035-------1030.00
-----70---------------------1089-------1084.00
-----80---------------------1123-------1115.00
-----90---------------------1157-------1146.00

If you look at color-color plots of what was targeted by both version (white), with version 1 (cyan) versus the new version (magenta):


We do seem to be better at targeting QSOs at redshift ~ 3.0 with this version of the method.


No comments:

Post a Comment