In [1]:

    

%matplotlib inline
from __future__ import division
import om10,os
import numpy as np
import matplotlib.pyplot as plt
import triangle


    

In [5]:

    

db = om10.DB(catalog='$OM10_DIR/data/qso_mock.fits')
db.lenses


    

    
Out[5]:




LENSID
FLAGTYPE
NIMG
ZLENS
VELDISP
...
SIGCRIT
DSLUM
L_I
REFF
REFF_T
14428
0
2
0.062
264.2409
...
0.0
35956.6
0.0
0.0
0.0
21703
0
2
0.072
165.1962
...
0.0
8193.9
0.0
0.0
0.0
31711
0
2
0.082
155.9553
...
0.0
35956.6
0.0
0.0
0.0
34331
0
2
0.084
172.5838
...
0.0
14498.4
0.0
0.0
0.0
51810
0
2
0.096
244.9063
...
0.0
26239.2
0.0
0.0
0.0
83414
0
2
0.114
154.17
...
0.0
11904.1
0.0
0.0
0.0
115074
0
2
0.128
113.5011
...
0.0
7101.16
0.0
0.0
0.0
125027
0
2
0.13
247.1724
...
0.0
16302.5
0.0
0.0
0.0
147709
0
2
0.138
200.9093
...
0.0
17748.4
0.0
0.0
0.0
180815
0
2
0.148
185.3532
...
0.0
14029.4
0.0
0.0
0.0
202550
0
4
0.154
181.134
...
0.0
8794.67
0.0
0.0
0.0
...
...
...
...
...
...
...
...
...
...
...
113167635
0
2
1.936
250.6109
...
0.0
47637.1
0.0
0.0
0.0
115614728
0
4
1.96
242.1029
...
0.0
38304.5
0.0
0.0
0.0
118666840
0
2
1.99
207.9697
...
0.0
27265.9
0.0
0.0
0.0
120546345
0
2
2.008
281.8383
...
0.0
28090.1
0.0
0.0
0.0
123598414
0
2
2.038
199.5262
...
0.0
31000.1
0.0
0.0
0.0
129002676
0
2
2.09
236.0478
...
0.0
39270.7
0.0
0.0
0.0
138366386
0
2
2.18
227.5097
...
0.0
36169.7
0.0
0.0
0.0
157075130
0
2
2.358
258.8213
...
0.0
29955.8
0.0
0.0
0.0
168229053
0
2
2.464
200.9093
...
0.0
48954.0
0.0
0.0
0.0
185790887
0
2
2.63
215.7744
...
0.0
42184.9
0.0
0.0
0.0
193208710
0
2
2.7
259.4179
...
0.0
37235.8
0.0
0.0
0.0

In [6]:

    

test_lens = db.get_lens(5924799)
test_lens


    

    
Out[6]:




LENSID
FLAGTYPE
NIMG
ZLENS
VELDISP
...
SIGCRIT
DSLUM
L_I
REFF
REFF_T
5924799
0
4
0.516
235.5049
...
0.0
12636.8
0.0
0.0
0.0

In [7]:

    

om10.plot_lens(test_lens)


    

    




om10.plot_lens: plotting image configuration of lens ID  5924799
om10.plot_lens: lens, image magnitudes: 19.9553566 [ 19.55916786  20.40661523  21.98854659  19.53562427]






    

In [10]:

    

test_lens['XIMG'], test_lens['YIMG']


    

    
Out[10]:





(<Column name='XIMG' unit=None format=None description=None>
 array([[-0.2291,  0.5897, -1.1676,  0.3106]]),
 <Column name='YIMG' unit=None format=None description=None>
 array([[-0.7459, -0.6133,  0.3543,  0.4498]]))

In [8]:

    

input_cat = om10.DB(catalog='./input_cat_for_eric.fits')


    

In [10]:

    

test_sim = input_cat.get_lens(5924799)
test_sim


    

    
Out[10]:




LENSID
RA
DEC
...
I
Z
5924799.0
34.64855907
-7.843399404
...
18.92477
18.41967
5924799.0
34.6485409841
-7.84370237622
...
17.2676578609
17.2377178609
5924799.0
34.6482323468
-7.84338776511
...
18.1151052333
18.0851652333
5924799.0
34.6486761094
-7.84328101511
...
19.6970365894
19.6670965894
5924799.0
34.6484111865
-7.84366326511
...
17.2441142721
17.2141742721

In [11]:

    

plt.plot(3600*test_sim['RA'][0],3600*test_sim['DEC'][0],color='Orange',marker='+')
plt.plot(3600*test_sim['RA'][1:],3600*test_sim['DEC'][1:],color='Blue',marker='+',linestyle='None')


    

    
Out[11]:





[<matplotlib.lines.Line2D at 0x1092b6090>]






    

In [ ]: