In [13]:
from astropy.io.fits import getdata
from astropy.time import Time
from nustar_lunar_pointing.tracking import get_epoch_tle, get_moon_j2000
from astropy.coordinates import SkyCoord
import astropy.units as u

In [3]:
infile ='../data/numaster.fits'
tlefile = '../data/NusTAR.tle'

numaster = getdata(infile)

In [30]:
ctr = 0
for ind, t in enumerate(numaster['Time']):
    obs_time = Time(t, format='mjd')
    obs_ra = numaster[ind]['RA']
    obs_dec = numaster[ind]['DEC']
    if (obs_dec == -99):
        continue

    mindt, line1, line2 = get_epoch_tle(obs_time, tlefile)
    ra_moon, dec_moon = get_moon_j2000(obs_time, line1, line2)
    moon = SkyCoord(ra_moon, dec_moon)
        
    obs = SkyCoord(obs_ra*u.deg, obs_dec*u.deg)
    doff = obs.separation(moon).deg
    if (doff < 30):
        print(ind, numaster[ind]['NAME'], doff)
        ctr += 1


10 1E1740d7m2942 17.486514286115753
62 G21d5m0d9 28.266906976014354
63 G21d5m0d9 15.675594109083082
64 G21d5m0d9 15.127668288270966
238 Crab_cal12 19.984006178866675
319 COSMOS_MOS007 29.043495988766757
320 COSMOS_MOS008 22.71217503930299
322 B2_1023p25 22.611100057931587
328 COSMOS_MOS009 14.911181792336775
329 COSMOS_MOS010 22.619724166841834
392 MCGm6m30m15 21.580873093277535
393 MCGm6m30m15 21.49685300309463
431 3C105 26.85700265312411
433 NGC0513 24.00647386619671
441 Ark120 24.152417460753355
583 Crab_cal18 24.489628462694554
593 MCG_01_24_012 29.25485627113909
615 gcmagnetar 29.90101008773521
670 COSMOS_MOS061 16.722398907424804
671 COSMOS_MOS062 23.90525730246475
690 gcmagnetar 15.086872263158325
728 Norma_MOS014 27.751291529035118
729 Norma_MOS015 26.761766029108017
730 Norma_MOS016 28.48654543864051
837 GX339_4 29.78788325096706
841 GalSurvey_B_MOS001 23.93618339336372
893 1E1841m045 13.838149749790201
935 GalSurvey_B_MOS013 26.939467041899857
942 GalSurvey_B_MOS014 17.86200133584996
972 XTEJ1908p094 27.114607450923646
1054 NGC4388 29.188368899429847
1094 COSMOS_MOS111 26.11296838388768
1104 COSMOS_MOS112 29.555674064747805
1132 Mrk728 19.148666384648962
1134 PG_1211p143 18.516879927235117
1358 HE1136m2304 29.580367783222957
1405 SDSSJ1243m0232 29.068762104109112
1440 NGC1068 29.740873949178262
1531 F2M0830p3759 29.986409545093256
1553 GalSurvey_B_MOS019 22.90997406115764
1612 SS433 23.04580932272511
1641 2MASXJ04532576p0403416 19.87462841486597
1642 2MASXJ04532576p0403416 21.39566520279529
1663 SS433 23.439005079943886
1784 NGC5135 22.020957517460154
1785 NGC5135 23.01296792364984
1886 HESSJ1640m465 29.864304537798564
1940 IGR_J17511m3057 22.99167635605468
1946 HESSJ1834m087_TeV 14.757737845201088
2032 CXOUJ171405d7m381031 25.79074325852613
2119 SGR_1806m20 29.217787383274693
2121 HESSJ1640m465 29.283247895890252
2125 4U_1909p07 28.96966804877288
2177 H1743m322 22.923980543195082
2246 IGR_J18027m2016 24.210784382731912
2248 PG2112p059 21.12092025469217
2394 Arches_cluster 14.366616287255273
2455 2MASXJ18305065p0928414 27.624207536576073
2473 2MASXJ01392400p2924067 28.81950811843348
2475 1ES_0414p009 15.286413828421383
2493 NGC2992 24.130415656597986
2543 1RXSJ073308d7p455511 28.74898640802971
2561 J133331d15m012653d3 27.356523595688998
2641 UDS_MOS011 26.849432274428853
2642 UDS_MOS012 16.300914251274577
2682 4U_1700m37 24.590318987765322
2719 IGRJ17091m3624 20.855744569643647
2721 RXJ1713d7m3946_NW_P2 22.678403442183505
2783 IGRJ17402m3656 19.758364681482217
2906 G11d2m0d3 28.447978507121732
2918 2MASXJ01485967p2145343 14.448246133276347
2920 J0113p0106 25.95859934780647
2956 GRS_1915p105 28.894174237426697
3060 1RXSJ010713d0m113554 26.563897537030208
3142 TCP_J18102829m2729590 17.837761355714335
3160 UDS_MOS042 25.19589369792038
3305 TON_396 15.006708399952515
3313 M87 20.36258715389459
3317 NGC4939 24.599524217250845
3331 SwiftJ1753d5m0127 16.96682961368693
3337 Swift_J1816d7m1613 20.501491681338734
3525 M87 26.632391684963245
3541 PSR_B1821m24 14.580686622421064
WARNING: ErfaWarning: ERFA function "d2dtf" yielded 1 of "dubious year (Note 5)" [astropy._erfa.core]
WARNING: ErfaWarning: ERFA function "dtf2d" yielded 1 of "dubious year (Note 6)" [astropy._erfa.core]
WARNING: ErfaWarning: ERFA function "utctai" yielded 1 of "dubious year (Note 3)" [astropy._erfa.core]
WARNING: ErfaWarning: ERFA function "taiutc" yielded 1 of "dubious year (Note 4)" [astropy._erfa.core]
WARNING: ErfaWarning: ERFA function "epv00" yielded 1 of "warning: date outsidethe range 1900-2100 AD" [astropy._erfa.core]
WARNING: Tried to get polar motions for times before IERS data is valid. Defaulting to polar motion from the 50-yr mean for those. This may affect precision at the 10s of arcsec level [astropy.coordinates.builtin_frames.utils]
WARNING: ErfaWarning: ERFA function "utcut1" yielded 1 of "dubious year (Note 3)" [astropy._erfa.core]
3581 2MASXJ08555426p0051110 22.351948374201736
3583 IC2515 25.293372384777392
3588 SDSSJ104208d36p004206d1 13.792001911132301
3593 2MASXJ11324928p1017473 19.677711192466123
3595 CGCG041_020 27.185436982941507
3597 NGC4074 27.82939783779156
3787 2MASXJ08420557p0759253 22.79506654244454
3788 2MASXJ08451850p1420345 21.40117987186948
3792 1RXSJ090915d6p035453 17.919115694823958
3798 Mrk705 11.455815171342158
3799 NGC2885 15.152773193356612
3801 2MASXJ09352707p2617093 16.410383144372688
3803 CGCG122m055 13.415899332281828
3804 NGC3035 19.972211177551003
3806 NGC3080 3.2408294904606714
3808 3C234d0 16.541450967878653
3811 2MASXJ10084862m0954510 22.38843388784685
3812 ARK241 16.04927083858927
3813 NGC3227 7.832126008064174
3816 2MASSiJ1031543m141651 27.13811605592872
3819 FBQSJ110340d2p372925 27.482915242587914
3826 Mrk423 28.253963033735406
3827 2MASXJ11271632p1909198 19.066226692450336
3831 HB89_1130p106 19.704869691843008
3833 2MASXJ11383367p2523532 23.963112589283153
3834 SDSSJ113915d13p253557d9 24.196614777621157
3836 PG1138p222 23.07997319943807
3840 2MASXJ11475508p0902284 23.55241350736741
3841 MCGp05m28m032 28.0791585393877
4057 ABELL754 27.235306941927647
4059 2MASXJ10171680m0404558 16.566493132276822
4138 Mrk1239 14.969794277153493
4151 CGCG_004m040 29.770416473092155

In [32]:
print(ctr, ctr / ind)


116 0.027586206896551724

In [ ]: