notebook.community

Edit and run

Reading FRB Cat [v1.3]

v1.2 -- Update in 2018
v1.3 -- Update on 07 Dec 2019

http://frbcat.org/




In [1]:



    


# imports
from imp import reload
import numpy as np
from matplotlib import pyplot as plt

from astropy.table import Table
from astropy.time import Time

from frb import frbcat as ffcat
from frb import experiment

Load with 2018 format





In [2]:



    


reload(ffcat)
obs = ffcat.FRBCat()



    


















    






Using /home/xavier/Projects/FRB/frb/data/FRBs/frbcat_2019-12-07.csv for the FRB catalog

















In [3]:



    


obs.frbcat.sort('frb_name')



    











In [4]:



    


obs.frbcat[-5:]



    


















    
Out[4]:






Table length=5



frb_nameutctelescoperop_rajrop_decjrop_glrop_gbrmp_dmrmp_widthrmp_snrRADEC


str18str23str9str12str12float64float64str19float64float64float64float64


FRB1812282018-12-28 13:48:50.100UTMOST06:09:23.64-45:58:02.4253.3915-26.0633354.2&plusmn0.91.2412.092.47723726394653-45.98330613177865


FRB190116.J1249+272019-01-16 13:07:33.833CHIME/FRB12:4927:09210.589.5444&plusmn0.64.00.0192.2980758081284627.14835236872834


FRB190209.J0937+772019-02-09 08:20:20.977CHIME/FRB09:3777:40134.234.8424.6&plusmn0.63.70.0144.1099325933291577.67748003429321


FRB190222.J2052+692019-02-22 18:46:01.367CHIME/FRB20:5269:50104.915.9460.6&plusmn0.12.970.0313.046679135227269.85418471743013


FRB1905232019-05-23 06:05:55.815DSA-1013:48:15.6+72:28:11117.0344.0760.8&plusmn0.60.4211.5207.07557574047872.47075196238896



















In [5]:



    


obs.frbcat[0:20]



    


















    
Out[5]:






Table length=20



frb_nameutctelescoperop_rajrop_decjrop_glrop_gbrmp_dmrmp_widthrmp_snrRADEC


str18str23str9str12str12float64float64str19float64float64float64float64


FRB0101252001-01-25 00:29:15.790parkes19:06:53-40:37:14356.641-20.0206790&plusmn39.417.0286.72072305094264-40.62064517309942


FRB0103122001-03-12 11:06:47.980parkes05:26:54.9-64:56:19.2274.72-33.31187&plusmn1424.311.081.71890591818621-64.93631423273322


FRB0106212001-06-21 13:02:11.299parkes18:52:05-08:29:3525.434-4.00381745&plusmn107.016.3283.0208291465951-8.493049586958135


FRB0107242001-07-24 19:50:01.690parkes01:18:06-75:12:18300.653-41.80513755.023.019.525065904947255-75.20496347123685


FRB0906252009-06-25 21:53:51.379parkes03:07:47-29:55:36226.444-60.0303899.55&plusmn0.011.9230.046.94584178165078-29.926709866601552


FRB1102142011-02-14 07:14:10.353parkes01:21:17-49:47:11290.7-66.6168.9&plusmn0.51.913.020.351647596049716-49.81046938710202


FRB1102202011-02-20 01:55:48.096parkes22:34:38-12:23:4550.829-54.7663944.38&plusmn0.055.649.0338.65833967919855-12.395829957098739


FRB1105232011-05-23 15:06:19.700GBT21:45:12-00:09:3756.12-37.82623.3&plusmn0.061.7342.0326.2965238534899-0.15926998703343367


FRB1106262011-06-26 21:33:17.477parkes21:03:43-44:44:19355.862-41.7522723&plusmn0.31.411.0315.9291924469609-44.73847978976432


FRB1107032011-07-03 18:59:40.607parkes23:30:51-02:52:2480.9978-59.01911103.6&plusmn0.74.316.0352.71250345649395-2.8733513425079926


FRB1201272012-01-27 08:11:21.725parkes23:15:06-18:25:3849.2871-66.2037553.3&plusmn0.31.111.0348.7750050379304-18.42723215173697


FRB1210022012-10-02 13:09:18.436parkes18:14:47-85:11:53308.22-26.26471629.18&plusmn0.025.4416.0273.6967687919654-85.19777502573858


FRB1211022012-11-02 06:35:53.244arecibo05:32:0933:05:13174.95-0.225138557&plusmn23.014.083.0373314739749933.08715440393463


FRB1306262013-06-26 14:55:59.771parkes16:27:06-07:27:487.4500327.4203952.4&plusmn0.11.9821.0246.77495842935215-7.4633132463924605


FRB1306282013-06-28 03:58:00.178parkes09:03:0203:26:16225.95530.6556469.88&plusmn0.010.6429.0135.75811865081833.438101042303216


FRB1307292013-07-29 09:01:51.190parkes13:41:21-05:59:43324.78854.7446861&plusmn215.6114.0205.33766462069678-5.995193975709084


FRB1311042013-11-04 18:04:11.200parkes06:44:10-51:16:40260.55-21.9253779&plusmn12.0830.0101.0403428601284-51.271538619174756


FRB1405142014-05-14 17:14:11.060parkes22:34:06-12:18:4650.8413-54.612562.7&plusmn0.62.816.0338.5250223506901-12.312768823204856


FRB1411132014-11-13 07:42:55.220arecibo06:13:00.118:47:11.2191.90.36400.32.08.493.2484012452370618.799947101323646


FRB1502152015-02-15 20:41:41.714parkes18:17:27-04:54:1524.66285.280921105.6&plusmn0.82.8819.0274.36250785471003-4.9041513275113635

Simple histogram





In [6]:



    


times = Time(obs.frbcat['utc'])



    











In [7]:



    


times.jyear



    


















    
Out[7]:








array([2001.06644852, 2001.19360179, 2001.47034411, 2001.56146861,
       2009.48230637, 2011.12060646, 2011.13642825, 2011.38981354,
       2011.48363619, 2011.50250908, 2012.07074941, 2012.75304074,
       2012.83716611, 2013.48424975, 2013.48847441, 2013.57392549,
       2013.843266  , 2014.36610043, 2014.86604099, 2015.12487964,
       2015.29277723, 2015.43799336, 2015.59820886, 2015.79341737,
       2015.92816541, 2015.9950106 , 2016.00233603, 2016.20773559,
       2016.27339277, 2016.43439238, 2016.71903893, 2017.01940287,
       2017.29080389, 2017.3230713 , 2017.42893652, 2017.51338043,
       2017.5278772 , 2017.65415678, 2017.68116766, 2017.72477671,
       2017.7540638 , 2017.75671848, 2017.79893275, 2017.80133086,
       2017.87576917, 2017.93937636, 2017.94962103, 2017.95824622,
       2018.02550495, 2018.05069716, 2018.07403537, 2018.0744799 ,
       2018.0799595 , 2018.0827916 , 2018.1176992 , 2018.18375837,
       2018.18939066, 2018.20044398, 2018.22559088, 2018.29173039,
       2018.32694615, 2018.36937684, 2018.395998  , 2018.40296603,
       2018.5322841 , 2018.56331194, 2018.56683412, 2018.57420394,
       2018.57230269, 2018.57536207, 2018.58142681, 2018.59573552,
       2018.6070593 , 2018.60765212, 2018.61188217, 2018.61770819,
       2018.61765199, 2018.62443786, 2018.70753542, 2018.73013767,
       2018.78898954, 2018.79242646, 2018.79391246, 2018.82731238,
       2018.86442174, 2018.88350645, 2018.90719388, 2018.98993999,
       2019.04188068, 2019.10704303, 2019.14382467, 2019.38878609])

















In [8]:



    


bins = np.arange(2001, 2020,1)
bins



    


















    
Out[8]:








array([2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008, 2009, 2010, 2011,
       2012, 2013, 2014, 2015, 2016, 2017, 2018, 2019])

















In [9]:



    


plt.clf()
ax = plt.gca()
ax.hist(times.jyear, bins=bins)
#
plt.show()

Load an old one





In [10]:



    


reload(ffcat)
obs = ffcat.FRBCat(frbcat_file='frbcat_2017-04-06.csv')
obs



    


















    






Using /home/xavier/Projects/FRB/frb/data/FRBs/frbcat_2017-04-06.csv for the FRB catalog










    
Out[10]:








<FRBCat:  FRB Catalog with 21 sources
>

















In [11]:



    


obs.frbcat[0:5]



    


















    
Out[11]:






Table masked=True length=5



NameTelescopeTypeUTCBeamReceiverBackendRAJDECJPointing ErrorFWHMSampling TimeBandwidthCentre FrequnecyBits per SampleGainSystem TemperatureNE2001 DM LimitDMDM errorSNRWidthWidth error lowerWidth error upperFluxFlux error lowerFlux error upperDM IndexDM Index ErrorScattering IndexScattering Index ErrorScattering TimeScattering Time ErrorLinear Polarization FractionLinear Polarization Fraction ErrorCircular Polarization FractionCircular Polarization Fraction ErrorPhotometric RedshiftPhotometric Redshift ErrorSpectroscopic RedshiftSpectroscopic Redshift ErrorReferenceRADEC


str9str7str5str19int64str6str18str8str9int64float64float64float64int64int64float64float64float64float64float64float64float64float64float64float64float64float64float64float64float64float64float64float64float64float64int64float64float64float64float64float64str97float64float64


FRB010125parkesradio2001-01-25 00:29:145MB20AFB19:06:53-40:37:141115.00.125288.0--10.6928.0110.0790.03.017.09.40.20.20.3----2.00.01-4.21.2--------------------Burke-Spolaor S., Bannister K. W., 2014, ApJ, 792, 19286.7208333333333-40.620555555555555


FRB010125parkesradio2001-01-25 00:29:145MB20AFB19:06:53-40:37:141115.00.125288.0--10.6928.0110.0790.33.025.010.62.52.80.540.070.11----------------------------Burke-Spolaor S., Bannister K. W., 2014, ApJ, 792, 19286.7208333333333-40.620555555555555


FRB010621parkesradio2001-06-21 13:02:0910MB20AFB18:52:05-08:29:351115.00.25288.0--10.58128.0523.0748.03.018.08.02.254.00.530.090.26----------------------------Keane E. F., Kramer M., Lyne A. G., Stappers B. W., McLaughlin M. A., 2011, MNRAS, 415, 3065283.0208333333333-8.493055555555555


FRB010621parkesradio2001-06-21 13:02:0910MB20AFB18:52:05-08:29:351115.00.25288.0--10.58128.0523.0745.010.0--7.0----0.41----2.020.01------------------------Keane E. F., Kramer M., Lyne A. G., Stappers B. W., McLaughlin M. A., 2011, MNRAS, 415, 3065283.0208333333333-8.493055555555555


FRB010724parkesradio2001-07-24 19:50:006MB20AFB01:18:06-75:12:181115.01.0288.0--10.6928.044.58375.03.0100.020.00.00.01.5740.00.0----------------------------Lorimer D. R., Bailes M., McLaughlin M. A., Narkevic, D. J., Crawford F., 2007, Science, 318, 77719.525-75.205

For Sky





In [15]:



    


obs = ffcat.FRBCat()



    


















    






Using /home/xavier/Projects/FRB/frb/data/FRBs/frbcat_2019-12-07.csv for the FRB catalog

















In [19]:



    


obs.frbcat.sort('utc')



    











In [21]:



    


obs.frbcat.write('sky_tbl.fits')



    











In [ ]:

Content source: FRBs/FRB

Similar notebooks:

notebook.community | gallery | about