Vetting the DR5

v6.1



In [3]:

    
def lya(z,zem=None):
    wvobs = (1+z)*1215.67
    print('wvobs = {:g}'.format(wvobs))
    if zem is not None:
        wvrest = wvobs/(1+zem)
        print('wvrest = {:g}'.format(wvrest))
    return wvobs

Examples

493, 507

Double DLA
Missed the low z case



In [4]:

    
z=2.22131
zem = 2.569999
lya(z,zem)









    



wvobs = 3916.05
wvrest = 1096.93

604, 535 -- High NHI

Some residual flux
And a skyline spike
But this one *should* be found

1415,423

Double DLA
Both are weak and could be bogus..

1165, 582

Large NHI
Skyline spike
Should be found!

1226, 638

Large NHI
Unclear what happened here!
Too large?

1165, 426

Skyline

1408, 460

z=5
dodgy..

1936, 618

Large NHI
Seems pretty obvious..

787, 48

Non-zero flux
good miss

778, 38

Large NHI
high z
Should be found!

2010, 536

Lower NHI but still obvious
bad sky nearby



In [5]:

    
z=2.58369
zem = 3.092
lya(z,zem)









    



wvobs = 4356.58
wvrest = 1064.66






    Out[5]:





4356.5844223

1920, 273

Nearby skyline

2136, 438

Really strange..
Drawn away by Lya emission?

2099, 635

Near zem, but obvious

876, 346

Bad z??
Or have I screwed up??



In [4]:

    
z=3.14621
zem=3.1881999
lya(z,zem)









    



wvobs = 5040.42
wvrest = 1203.48






    Out[4]:





5040.4231107000005

v2_gensample

Same olds

787, 48

Ones for JXP to debug

836, 603
1433, 54

1431, 358

Double very close
Missed both, oddly

425, 278

Low NHI, but high S/N

397,333

Just missed (funny core)

517, 259

Double

770, 291

Performed but missed

1370, 335

Low S/N
Almost got it

2136, 438

Almost but no luck



In [ ]: