Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Study efficiency - purity tradeoffs #2

Open
dkirkby opened this issue Mar 13, 2016 · 1 comment
Open

Study efficiency - purity tradeoffs #2

dkirkby opened this issue Mar 13, 2016 · 1 comment

Comments

@dkirkby
Copy link
Owner

dkirkby commented Mar 13, 2016

Define a boolean truth selector:

T(v1) =  { c | zfit - ztrue | / (1 + ztrue) < v1 }

and a boolean data selector, for example:

D(v2) = { c (z68hi - z68lo) / 2 / (1 + zbest) < v2 }

Then we can count the numbers of fits in four exclusive categories, which sum to N(total):

  • Positively selected in data: N(D)
    • True positives: N(D & T)
    • False positives: N(D & ~T)
  • Negatively selected in data: N(~D)
    • True negatives: N(~D & ~T)
    • False negatives: N(~D & T)

Focus on two metrics:

  • efficiency = N(D & T) / N(total)
  • purity = N(D & T) / N(D)

This issue is to:

  • Write code to calculate curves of efficiency and purity versus v2 given a definition of D(v2) and a fixed v1.
  • Find the optimum definition of D(v2) and choice of v2 given a minimum allowed efficiency and purity for some v1.

The optimum might be different for different target categories and can be used to set ZWARN = ~D(v2).

@fjaviersanchez
Copy link
Collaborator

The results for the different object types show a small improvement but we lose some efficiency, and at this point it might not be worth to make the cuts until we have more realistic simulations. The biggest improvement has been to include z_avg (average of the posterior) instead of z_best (maximum of posterior) as our redshift estimate.

ELGs

elg_nfail
elg_nwarn
elg_s68cut_dv
elg_s95cut_dv

LRG

lrg_nfail
lrg_nwarn
lrg_s68cut_dv
lrg_s95cut_dv

QSO

qso_nfail
qso_nwarn
qso_s68cut_dv
qso_s95cut_dv

STAR

star_nfail
star_nwarn
star_s68cut_dv
star_s95cut_dv

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants