Reading FRB Cat [v1.2]

v1.2 -- Update in 2018

In [53]:
# 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 [45]:
reload(ffcat)
obs = ffcat.FRBCat()


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

In [46]:
obs.frbcat[-5:]


Out[46]:
Table length=5
frb_nameutctelescoperop_rajrop_decjrop_glrop_gbrmp_dmrmp_widthrmp_snrRADEC
str11str23str7str11str11float64float64str18float64float64float64float64
FRB1102142011-02-14 07:14:10.353parkes01:21:17-49:47:11290.7-66.6168.9&plusmn0.51.913.020.351647596049716-49.81046938710202
FRB0906252009-06-25 21:53:51.379parkes03:07:47-29:55:36226.444-60.0303899.55&plusmn0.011.9230.046.94584178165078-29.92670986660156
FRB0107242001-07-24 19:50:01.690parkes01:18:06-75:12:18300.653-41.80513755.023.019.525065904947255-75.20496347123685
FRB0106212001-06-21 13:02:11.299parkes18:52:05-08:29:3525.434-4.00381745&plusmn107.016.3283.0208291465951-8.493049586958135
FRB0101252001-01-25 00:29:15.790parkes19:06:53-40:37:14356.641-20.0206790&plusmn39.417.0286.72072305094264-40.62064517309942

Simple histogram


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

In [50]:
times.jyear


Out[50]:
array([2018.395998  , 2018.22559088, 2018.20044398, 2018.1176992 ,
       2018.0827916 , 2018.0799595 , 2018.0744799 , 2018.07403537,
       2018.05069716, 2018.02550495, 2017.95824622, 2017.94962103,
       2017.87576917, 2017.80133086, 2017.79893275, 2017.75671848,
       2017.7540638 , 2017.68116766, 2017.65415678, 2017.5278772 ,
       2017.51338043, 2017.3230713 , 2017.29080389, 2017.01940287,
       2016.43439238, 2016.27339277, 2016.20773559, 2016.00233603,
       2015.9950106 , 2015.92816541, 2015.59820886, 2015.43799336,
       2015.29277723, 2015.12487964, 2014.86604099, 2014.36610043,
       2013.843266  , 2013.57392549, 2013.48847441, 2013.48424975,
       2012.83716611, 2012.75304074, 2012.07074941, 2011.50250908,
       2011.48363619, 2011.38981354, 2011.13642825, 2011.12060646,
       2009.48230637, 2001.56146861, 2001.47034411, 2001.06644852])

In [55]:
bins = np.arange(2001, 2020,1)
bins


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

In [60]:
plt.clf()
ax = plt.gca()
ax.hist(times.jyear, bins=bins)
#
plt.show()


Load an old one


In [3]:
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[3]:
<FRBCat:  FRB Catalog with 21 sources
>

In [4]:
obs.frbcat[0:5]


Out[4]:
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


In [7]:
from astropy.coordinates import SkyCoord

In [9]:
c = SkyCoord(b=23.2, l=232.2, frame='galactic', unit='deg')

In [10]:
c


Out[10]:
<SkyCoord (Galactic): (l, b) in deg
    (232.2, 23.2)>

In [ ]: