In [1]:
%matplotlib inline
import matplotlib.pyplot as plt
import numpy as np
from vespa import transit_basic as tr
/u/tdm/anaconda/lib/python2.7/site-packages/setuptools-7.0-py2.7.egg/pkg_resources.py:1045: UserWarning: /u/tdm/.python-eggs is writable by group/others and vulnerable to attack when used with get_resource_filename. Consider a more secure location (set with .set_extraction_path or the PYTHON_EGG_CACHE environment variable).
In [2]:
mafn = tr.MAInterpolationFunction()
In [3]:
zs = np.linspace(0,1.2,500)
plt.plot(zs, mafn(0.1,zs))
Out[3]:
[<matplotlib.lines.Line2D at 0x7f69ea91a650>]
In [6]:
import logging
rootLogger = logging.getLogger()
rootLogger.setLevel(logging.WARNING)
In [9]:
p0, b, aR = tr.eclipse_pars(10, 1.0, 0.5, 1.0, 0.5)
ts, fs = tr.eclipse(p0,b,aR, P=10, MAfn=mafn)
plt.plot(ts,fs)
dur, depth, slope = tr.eclipse_tt(p0, b, aR, P=10, MAfn=mafn)
plt.plot(ts,tr.traptransit(ts,np.array([dur,depth,slope,0])))
Out[9]:
[<matplotlib.lines.Line2D at 0x7f69ea8158d0>]
In [10]:
print dur, depth, slope
0.191708658273 0.289214147953 2.93380430212
In [ ]:
In [1]:
from vespa.kepler import KeplerTransitSignal
In [2]:
s = KeplerTransitSignal(87.01)
In [3]:
s.df
Out[3]:
TIME
TIMECORR
CADENCENO
SAP_FLUX
SAP_FLUX_ERR
SAP_BKG
SAP_BKG_ERR
PDCSAP_FLUX
PDCSAP_FLUX_ERR
SAP_QUALITY
PSF_CENTR1
PSF_CENTR1_ERR
PSF_CENTR2
PSF_CENTR2_ERR
MOM_CENTR1
MOM_CENTR1_ERR
MOM_CENTR2
MOM_CENTR2_ERR
POS_CORR1
POS_CORR2
0
120.539007
0.000855
568
301338.75000
14.351400
1732.762695
0.795363
314778.31250
14.928895
0
NaN
NaN
NaN
NaN
179.907726
0.000051
615.855901
0.000047
-0.001460
0.003533
1
120.559442
0.000856
569
301317.34375
14.350985
1738.260620
0.796874
314749.90625
14.942488
0
NaN
NaN
NaN
NaN
179.907888
0.000051
615.855985
0.000047
-0.001193
0.003698
2
120.579876
0.000856
570
301328.18750
14.352235
1738.363281
0.795001
314776.43750
14.936597
0
NaN
NaN
NaN
NaN
179.908033
0.000051
615.855737
0.000047
-0.001002
0.003532
3
120.600310
0.000857
571
301302.81250
14.351120
1737.150513
0.798276
314716.71875
14.892340
128
NaN
NaN
NaN
NaN
179.908521
0.000051
615.855697
0.000047
-0.000561
0.003343
4
120.620745
0.000858
572
301303.90625
14.351894
1737.046997
0.796383
314750.46875
14.961715
0
NaN
NaN
NaN
NaN
179.908486
0.000051
615.855425
0.000047
-0.000296
0.003350
5
120.641179
0.000859
573
301323.93750
14.351185
1735.657349
0.798350
314744.43750
14.912598
0
NaN
NaN
NaN
NaN
179.908748
0.000051
615.855266
0.000047
-0.000082
0.002904
6
120.661613
0.000859
574
301302.12500
14.350945
1737.417969
0.796073
314732.96875
14.912435
128
NaN
NaN
NaN
NaN
179.908803
0.000051
615.855461
0.000047
0.000100
0.003117
7
120.682048
0.000860
575
301318.18750
14.351001
1734.431763
0.795743
314747.81250
14.932611
0
NaN
NaN
NaN
NaN
179.908870
0.000051
615.855271
0.000047
0.000091
0.002942
8
120.702482
0.000861
576
301329.71875
14.351317
1734.204834
0.798131
314755.12500
14.897735
0
NaN
NaN
NaN
NaN
179.908820
0.000051
615.855325
0.000047
0.000077
0.003002
9
120.722916
0.000861
577
301289.18750
14.350432
1734.150879
0.796400
314713.28125
14.890948
0
NaN
NaN
NaN
NaN
179.909054
0.000051
615.855066
0.000047
0.000244
0.002669
10
120.743351
0.000862
578
301315.56250
14.351839
1735.950806
0.794790
314756.65625
14.935699
0
NaN
NaN
NaN
NaN
179.909052
0.000051
615.855378
0.000047
0.000297
0.003181
11
120.763785
0.000863
579
301296.53125
14.350765
1735.640503
0.794837
314716.06250
14.937016
128
NaN
NaN
NaN
NaN
179.908892
0.000051
615.855278
0.000047
0.000180
0.002853
12
120.784219
0.000864
580
301293.56250
14.351567
1733.614136
0.797433
314735.75000
14.948170
0
NaN
NaN
NaN
NaN
179.908848
0.000051
615.855224
0.000047
0.000151
0.003064
13
120.804653
0.000864
581
301296.00000
14.351581
1733.131348
0.796645
314730.90625
14.936604
0
NaN
NaN
NaN
NaN
179.908948
0.000051
615.855453
0.000047
0.000232
0.003219
14
120.825088
0.000865
582
301308.59375
14.351114
1735.676636
0.795741
314739.46875
14.912021
128
NaN
NaN
NaN
NaN
179.908948
0.000051
615.855088
0.000047
0.000174
0.002847
15
120.845522
0.000866
583
301301.43750
14.351683
1736.945068
0.794803
314752.00000
14.989100
0
NaN
NaN
NaN
NaN
179.908783
0.000051
615.855314
0.000047
0.000067
0.003089
16
120.865956
0.000867
584
301321.65625
14.352068
1736.778931
0.796134
314753.59375
14.916309
0
NaN
NaN
NaN
NaN
179.908937
0.000051
615.855483
0.000047
0.000224
0.003166
17
120.886391
0.000867
585
301282.59375
14.351349
1736.359009
0.795323
314729.00000
14.962685
8192
NaN
NaN
NaN
NaN
179.909156
0.000051
615.855102
0.000047
0.000356
0.003081
18
120.906825
0.000868
586
301293.37500
14.351609
1742.116577
0.793769
314739.21875
14.944017
0
NaN
NaN
NaN
NaN
179.909132
0.000051
615.855239
0.000047
0.000448
0.003255
19
120.927259
0.000869
587
301288.43750
14.351261
1735.490356
0.795046
314725.40625
14.906545
0
NaN
NaN
NaN
NaN
179.909403
0.000051
615.855282
0.000047
0.000711
0.003160
20
120.947694
0.000869
588
301274.31250
14.351209
1733.981323
0.796705
314727.75000
14.936511
0
NaN
NaN
NaN
NaN
179.909226
0.000051
615.855444
0.000047
0.000610
0.003302
21
120.968128
0.000870
589
301282.34375
14.349618
1735.573364
0.796179
314734.37500
14.917085
0
NaN
NaN
NaN
NaN
179.909351
0.000051
615.855838
0.000047
0.000776
0.003796
22
120.988562
0.000871
590
301310.34375
14.350084
1733.979858
0.794329
314758.28125
14.939182
0
NaN
NaN
NaN
NaN
179.909502
0.000051
615.855969
0.000047
0.000859
0.003950
23
121.008997
0.000872
591
301263.75000
14.349236
1737.160400
0.794501
314730.46875
14.935127
0
NaN
NaN
NaN
NaN
179.909668
0.000051
615.856142
0.000047
0.001077
0.004249
24
121.029431
0.000872
592
301308.81250
14.350004
1733.037842
0.795205
314754.75000
14.880700
0
NaN
NaN
NaN
NaN
179.909813
0.000051
615.856045
0.000047
0.001219
0.003925
25
121.049865
0.000873
593
301283.71875
14.349568
1732.783325
0.796115
314756.90625
14.916881
0
NaN
NaN
NaN
NaN
179.909947
0.000051
615.856305
0.000047
0.001448
0.004523
26
121.070300
0.000874
594
301294.56250
14.349793
1736.393433
0.794646
314756.15625
14.942303
0
NaN
NaN
NaN
NaN
179.909832
0.000051
615.855991
0.000047
0.001258
0.004013
27
121.090734
0.000874
595
301268.18750
14.349350
1732.714111
0.796738
314732.31250
14.918449
0
NaN
NaN
NaN
NaN
179.909851
0.000051
615.856039
0.000047
0.001296
0.004131
28
121.111168
0.000875
596
301264.43750
14.349226
1737.358154
0.794555
314730.21875
14.936155
0
NaN
NaN
NaN
NaN
179.909994
0.000051
615.856135
0.000047
0.001436
0.004276
29
121.131603
0.000876
597
301239.53125
14.348791
1737.773438
0.795076
314715.68750
14.943715
0
NaN
NaN
NaN
NaN
179.909974
0.000051
615.856267
0.000047
0.001484
0.004458
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
1526
1590.408872
0.000294
72502
290612.00000
14.117769
1891.195435
0.796728
303680.00000
14.736322
0
179.913901
0.000038
615.869673
0.000032
179.945870
0.000053
615.879837
0.000049
0.000214
0.015142
1527
1590.429306
0.000295
72503
290614.62500
14.117785
1890.302490
0.795970
303702.12500
14.735714
0
179.913978
0.000038
615.869866
0.000032
179.946029
0.000053
615.879960
0.000049
0.000369
0.015356
1528
1590.449740
0.000295
72504
290618.56250
14.117774
1888.268677
0.795751
303701.15625
14.738446
0
179.914003
0.000038
615.869712
0.000032
179.945964
0.000053
615.879746
0.000049
0.000205
0.015235
1529
1590.470174
0.000296
72505
290632.90625
14.118023
1888.238770
0.796032
303697.90625
14.726659
0
179.913895
0.000038
615.869574
0.000032
179.945930
0.000053
615.879740
0.000049
0.000221
0.015094
1530
1590.490609
0.000297
72506
290647.03125
14.118378
1890.244507
0.795827
303734.25000
14.759251
0
179.913789
0.000038
615.869641
0.000032
179.945869
0.000053
615.879758
0.000049
0.000174
0.015215
1531
1590.511043
0.000297
72507
290641.96875
14.118310
1887.171997
0.796601
303713.34375
14.752502
0
179.914010
0.000038
615.869623
0.000032
179.945984
0.000053
615.879792
0.000049
0.000262
0.015143
1532
1590.531477
0.000298
72508
290650.09375
14.118516
1888.597168
0.796697
303716.93750
14.749522
0
179.913835
0.000038
615.869619
0.000032
179.945875
0.000053
615.879789
0.000049
0.000156
0.015114
1533
1590.551912
0.000299
72509
290620.84375
14.117930
1888.846436
0.795699
303706.00000
14.729953
0
179.913806
0.000038
615.869693
0.000032
179.945819
0.000053
615.879758
0.000049
0.000174
0.015228
1534
1590.572346
0.000300
72510
290639.28125
14.118202
1886.013672
0.796123
303702.46875
14.721685
0
179.913716
0.000038
615.869410
0.000032
179.945799
0.000053
615.879703
0.000049
-0.000030
0.014934
1535
1590.592780
0.000300
72511
290607.90625
14.117559
1888.696533
0.794004
303681.25000
14.734468
8192
179.913789
0.000038
615.869731
0.000032
179.945820
0.000053
615.879884
0.000049
0.000194
0.015219
1536
1590.613215
0.000301
72512
290627.00000
14.117975
1887.402954
0.795935
303699.78125
14.738303
0
179.913704
0.000038
615.869576
0.000032
179.945724
0.000053
615.879785
0.000049
0.000125
0.015078
1537
1590.633649
0.000302
72513
290658.81250
14.118552
1889.432373
0.795500
303714.40625
14.725591
0
179.913690
0.000038
615.869381
0.000032
179.945687
0.000053
615.879701
0.000049
-0.000009
0.014920
1538
1590.654083
0.000302
72514
290647.25000
14.118498
1887.824463
0.798430
303724.78125
14.739362
0
179.913877
0.000038
615.869596
0.000032
179.945899
0.000053
615.879714
0.000049
0.000307
0.015165
1539
1590.674518
0.000303
72515
290625.62500
14.117902
1890.046143
0.796249
303695.28125
14.755366
0
179.913924
0.000038
615.869645
0.000032
179.945962
0.000053
615.879792
0.000049
0.000291
0.015145
1540
1590.694952
0.000304
72516
290632.03125
14.117958
1887.885742
0.794197
303685.46875
14.735544
0
179.913559
0.000038
615.869349
0.000032
179.945611
0.000053
615.879682
0.000049
-0.000125
0.014868
1541
1590.715386
0.000304
72517
290636.03125
14.118173
1888.309082
0.795350
303709.59375
14.749816
0
179.913740
0.000038
615.869599
0.000032
179.945782
0.000053
615.879758
0.000049
0.000127
0.015136
1542
1590.735820
0.000305
72518
290663.62500
14.118673
1886.197876
0.796545
303716.75000
14.718997
8192
179.913895
0.000038
615.869569
0.000032
179.945783
0.000053
615.879802
0.000049
0.000229
0.015058
1543
1590.756255
0.000306
72519
290655.06250
14.118562
1886.248901
0.797058
303718.37500
14.725267
0
179.913785
0.000038
615.869460
0.000032
179.945782
0.000053
615.879659
0.000049
0.000120
0.014998
1544
1590.776689
0.000307
72520
290631.06250
14.118109
1886.391235
0.797105
303694.65625
14.732637
0
179.913626
0.000038
615.869552
0.000032
179.945695
0.000053
615.879711
0.000049
0.000007
0.015049
1545
1590.797123
0.000307
72521
290620.03125
14.117822
1886.728882
0.798193
303666.62500
14.723521
0
179.913920
0.000038
615.869453
0.000032
179.945834
0.000053
615.879729
0.000049
0.000245
0.014956
1546
1590.817558
0.000308
72522
290622.09375
14.117957
1883.789062
0.800048
303688.65625
14.747580
0
179.913761
0.000038
615.869530
0.000032
179.945788
0.000053
615.879750
0.000049
0.000144
0.015071
1547
1590.837992
0.000309
72523
290634.40625
14.118080
1889.025513
0.796464
303694.78125
14.711802
0
179.913811
0.000038
615.869508
0.000032
179.945806
0.000053
615.879732
0.000049
0.000200
0.015026
1548
1590.858426
0.000309
72524
290630.59375
14.117905
1884.987427
0.795485
303673.84375
14.710948
0
179.913931
0.000038
615.869368
0.000032
179.945896
0.000053
615.879668
0.000049
0.000224
0.014923
1549
1590.878861
0.000310
72525
290614.15625
14.117755
1885.464233
0.798310
303673.18750
14.733409
0
179.913742
0.000038
615.869471
0.000032
179.945780
0.000053
615.879585
0.000049
0.000047
0.015008
1550
1590.899295
0.000311
72526
290635.03125
14.118193
1883.331665
0.798422
303683.06250
14.732606
0
179.913709
0.000038
615.869394
0.000032
179.945721
0.000053
615.879640
0.000049
0.000129
0.014936
1551
1590.919729
0.000311
72527
290616.59375
14.117698
1890.242310
0.796420
303661.34375
14.725244
8192
179.913913
0.000038
615.869520
0.000032
179.945836
0.000053
615.879783
0.000049
0.000265
0.015031
1552
1590.940164
0.000312
72528
290627.50000
14.118028
1886.663574
0.797689
303686.62500
14.730895
0
179.913748
0.000038
615.869509
0.000032
179.945834
0.000053
615.879687
0.000049
0.000108
0.015009
1553
1590.960598
0.000313
72529
290641.71875
14.118085
1885.665771
0.795836
303683.12500
14.711903
0
179.913966
0.000038
615.869451
0.000032
179.945843
0.000053
615.879769
0.000049
0.000189
0.014963
1554
1590.981032
0.000314
72530
290601.40625
14.118630
1886.195190
0.796545
303685.03125
14.722792
128
179.911801
0.000038
615.868524
0.000032
179.944361
0.000053
615.879106
0.000049
-0.000820
0.014272
1555
1591.001467
0.000314
72531
290626.43750
14.118395
1886.976440
0.798337
303706.06250
14.736184
0
179.913053
0.000038
615.869359
0.000032
179.945443
0.000053
615.879514
0.000049
0.000105
0.014914
69778 rows × 20 columns
In [ ]:
Content source: timothydmorton/VESPA
Similar notebooks: