def frb_190523(): """Ravi+19, Nature -- https://ui.adsabs.harvard.edu/abs/2019Natur.572..352R/abstract """ fname = 'FRB190523' frb190523 = frb.FRB(fname, 'J134815.6+722811', 760.8 * units.pc / units.cm**3, z_frb=0.660, repeater=False) # Error ellipse [REQUIRED] frb190523.set_ee( 4, 1.5, cl=68., theta=340.) # Halfing the 3x8'' at 95% reported by Ravi et al. # Error in DM frb190523.DM_err = 0.6 * units.pc / units.cm**3 # Fluence frb190523.fluence = 280 * units.Jy * units.ms #frb180924.fluence_err = 1 * units.Jy * units.ms # NE2001 frb190523.set_DMISM() # RM #frb190102.RM = 10 * units.rad / units.m**2 #frb190102.RM_err = 1 * units.rad / units.m**2 # References frb190523.refs = ['Ravi2019'] # Write path = resource_filename('frb', 'data/FRBs') frb190523.write_to_json(path=path)
def frb_190102(): """Bhandari+20, ApJL -- """ fname = 'FRB190102' wv_oiii = 6466.48 z_OIII = wv_oiii / 5008.239 - 1 frb190102 = frb.FRB(fname, 'J212939.7-792832.3', # Pulled from Slack 26-03-2019 362 * units.pc / units.cm**3, z_frb=z_OIII) # Error ellipse [REQUIRED] frb190102.set_ee(0.1, 0.1, 0., 95.) # Error in DM frb190102.DM_err = 1 * units.pc / units.cm**3 # NE2001 frb190102.set_DMISM() # RM #frb190102.RM = 10 * units.rad / units.m**2 #frb190102.RM_err = 1 * units.rad / units.m**2 # References frb190102.refs = ['Bhandari2019'] # Write path = resource_filename('frb', 'data/FRBs') frb190102.write_to_json(path=path)
def frb_190711(): """MacQuart+20, Day+2020 https://ui.adsabs.harvard.edu/abs/2020MNRAS.497.3335D/abstract """ fname = 'FRB190711' frb190711 = frb.FRB( fname, 'J215740.68-802128.8', # MacQuarter+2020, Day+2020 587.9 * units.pc / units.cm**3, # Day+2020 z_frb=0.52172, repeater=True) # Error ellipse frb190711.set_ee(0.12, 0.07, theta=90., cl=68.) # Statistical (Day+2020) frb190711.set_ee(0.38, 0.3, theta=90., cl=68., stat=False) # Systematic # Error in DM frb190711.DM_err = 1 * units.pc / units.cm**3 # NE2001 frb190711.set_DMISM() # RM -- Day+2020 frb190711.RM = 9 * units.rad / units.m**2 frb190711.RM_err = 2 * units.rad / units.m**2 frb190711.fluence = 34 * units.Jy * units.ms frb190711.fluence_err = 3 * units.Jy * units.ms # References frb190711.refs = ['MacQuart2020', 'Day2020'] # Write path = resource_filename('frb', 'data/FRBs') frb190711.write_to_json(path=path)
def frb_121102(): """ FRB 121102 All of the data currently comes from Tendulkar et al. 2017 https://ui.adsabs.harvard.edu/abs/2017ApJ...834L...7T/abstract or Chatterjee et al. 2017 Update for Marcote et al. 2017 05h31m58.7013s +/- 4 mas +33d08m52.5536s +/- 2.3 mas """ frb121102 = frb.FRB('FRB20121102', 'J053158.7013+330852.5536', 558.1*units.pc/units.cm**3, z_frb=0.19273, repeater=True) # NE2001 frb121102.set_DMISM() # RM frb121102.RM = 1.e5 * units.rad / units.m**2 # Pulse properties frb121102.set_pulse(1*units.GHz, Wi=3.0*units.ms, Wi_err=0.5*units.ms, tscatt=0.024*units.ms) # no errors given # Error ellipse frb121102.set_ee(0.004, 0.002, theta=90., cl=68.) frb121102.set_ee(a=0.0015, b=0.0015, theta=0., cl=68.,stat=False) # Marcote priv. corr # References frb121102.refs = ['Tendulkar2017', 'Marcote2017'] # Write path = resource_filename('frb', 'data/FRBs') frb121102.write_to_json(path=path)
def frb_190611(): """ Macquart al. Nature Day et al. 2020 """ FRB_190611_coord = SkyCoord('J212258.91-792351.3', # Day+2020 unit=(units.hourangle, units.deg)) frb190611 = frb.FRB('FRB190611', FRB_190611_coord, 332.6 * units.pc / units.cm**3) # Error ellipse [REQUIRED] frb190611.set_ee(0.7, 0.7, 0., 68.) # Error in DM frb190611.DM_err = 1 * units.pc / units.cm**3 # NE2001 frb190611.set_DMISM() # RM #frb190611.RM = 20 * units.rad / units.m**2 #frb190611.RM_err = 4 * units.rad / units.m**2 # References frb190611.refs = ['MacQuart2019', 'Day2020'] # Write path = resource_filename('frb', 'data/FRBs') frb190611.write_to_json(path=path)
def frb_190608(): """Bhandari+20, ApJL, Day+20 https://ui.adsabs.harvard.edu/abs/2020MNRAS.497.3335D/abstract Macquart al. Nature """ fname = 'FRB190608' frb190608 = frb.FRB( fname, "J221604.77-075353.7", # Pulled from Slack on 2020 Mar 18 340.05 * units.pc / units.cm**3, z_frb=0.1177805) # Taken from the SDSS table # Error ellipse [REQUIRED] frb190608.set_ee(0.19, 0.18, theta=90., cl=68.) # Statistsical (Day+2020) frb190608.set_ee(0.18, 0.18, theta=90., cl=68., stat=False) # Systematic # Error in DM frb190608.DM_err = 0.6 * units.pc / units.cm**3 # NE2001 frb190608.set_DMISM() # FRB properties frb190608.fluence = 26 * units.Jy * units.ms # From Bhandari+20 frb190608.fluence_err = 4 * units.Jy * units.ms # -- //-- frb190608.RM = 353 * units.rad / units.m**2 frb190608.RM_err = 2 * units.rad / units.m**2 frb190608.tau = 3.3 * units.ms frb190608.tau_err = 0.2 * units.ms # References frb190608.refs = ['Bhandari2020', 'Day2020'] # Write path = resource_filename('frb', 'data/FRBs') frb190608.write_to_json(path=path)
def frb_190711(): """MacQuart+20, Day+2020 """ fname = 'FRB190711' frb190711 = frb.FRB(fname, 'J215740.68-802128.8', # MacQuarter+2020, Day+2020 587.9 * units.pc / units.cm ** 3, # Day+2020 z_frb=0.52172) # Error ellipse frb190711.set_ee(0.3, 0.3, 0., 68.) # (Statistical) frb190711.set_ee(0.38, 0.3, 0., 68., stat=False) # Systematic # Error in DM frb190711.DM_err = 1 * units.pc / units.cm ** 3 # NE2001 frb190711.set_DMISM() # RM -- Day+2020 frb190711.RM = 9 * units.rad / units.m**2 frb190711.RM_err = 2 * units.rad / units.m**2 # References frb190711.refs = ['MacQuart2019', 'Day+2020'] # Write path = resource_filename('frb', 'data/FRBs') frb190711.write_to_json(path=path)
def frb_180916(): """ FRB 180916.J0158+65 All of the data currently comes from Marcote et al. 2020 https://ui.adsabs.harvard.edu/abs/2020Natur.577..190M/abstract """ coord = SkyCoord("01h58m00.75017s 65d43m00.3152s", frame='icrs') name = 'FRB180916' frb180916 = frb.FRB(name, coord, 348.76 * units.pc / units.cm**3, z_frb=0.0337) # Error ellipse frb180916.set_ee(0.0023, 0.0023, 0., 68.) # Error in DM frb180916.DM_err = 0.10 * units.pc / units.cm**3 # NE2001 frb180916.set_DMISM() # #frb180924.fluence = 16 * units.Jy * units.ms #frb180924.fluence_err = 1 * units.Jy * units.ms # # References frb180916.refs = ['Marcote2020'] # Write path = resource_filename('frb', 'data/FRBs') frb180916.write_to_json(path=path)
def frb_190523(): """Ravi+19, Nature -- https://ui.adsabs.harvard.edu/abs/2019Natur.572..352R/abstract """ fname = 'FRB190523' frb190523 = frb.FRB(fname, 'J134815.6+722811', 760.8 * units.pc / units.cm**3, z_frb=0.660) # Error ellipse [REQUIRED] frb190523.set_ee(5, 2, 0., 340) # JXP eyeball # Error in DM frb190523.DM_err = 0.6 * units.pc / units.cm**3 # Fluence frb190523.fluence = 280 * units.Jy * units.ms #frb180924.fluence_err = 1 * units.Jy * units.ms # NE2001 frb190523.set_DMISM() # RM #frb190102.RM = 10 * units.rad / units.m**2 #frb190102.RM_err = 1 * units.rad / units.m**2 # References frb190523.refs = ['Ravi2019'] # Write path = resource_filename('frb', 'data/FRBs') frb190523.write_to_json(path=path)
def frb_171020(): frbname = 'FRB20171020' FRB_171020_coord = SkyCoord("22h15m18.55s -19d40m11.23s", frame='icrs') # From Shannon+18 frb171020 = frb.FRB(frbname, FRB_171020_coord, 114.1 * units.pc / units.cm ** 3,z_frb=0.00867,repeater=False) # Error ellipse (Statistical) frb171020.set_ee(600,600, 0., 68.) # Error ellipse (Systematic) frb171020.set_ee(0., 0., 0., 68., stat=False) # Error in DM frb171020.DM_err = 0.2 * units.pc / units.cm ** 3 # NE2001 frb171020.set_DMISM() # RM # frb190102.RM = 10 * units.rad / units.m**2 # frb190102.RM_err = 1 * units.rad / units.m**2 # References frb171020.refs = ['Shannon2018','Mahony2018'] # Write path = resource_filename('frb', 'data/FRBs') frb171020.write_to_json(path=path)
def frb_200430(): """ Heintz+2020 Bhandari+2021 Updated FRB localization + error from Day+2021 """ frbname = 'FRB20200430' FRB_200430_coord = SkyCoord('15h18m49.54s +12d22m36.3s', frame='icrs') frb200430 = frb.FRB(frbname, FRB_200430_coord, 380 * units.pc / units.cm ** 3, # First Slack posting z_frb = 0.161, repeater=False) # Error ellipse (Statistical) frb200430.set_ee(0.24, 0.17, 0., 68.) # Error ellipse (Systematic) frb200430.set_ee(0.98, 0.251, 0., 68., stat=False) # Error in DM #frb191001.DM_err = 1 * units.pc / units.cm ** 3 # NE2001 frb200430.set_DMISM() # RM # frb190102.RM = 10 * units.rad / units.m**2 # frb190102.RM_err = 1 * units.rad / units.m**2 # References frb200430.refs = ['Bhandari2021', 'Day2021'] # Write path = resource_filename('frb', 'data/FRBs') frb200430.write_to_json(path=path)
def frb_191228(): """ Bhandari+2021 Updated FRB localization + error from Day+2021 """ frbname = 'FRB20191228' FRB_191228_coord = SkyCoord('22h57m43.30s -29d35m38.7s', frame='icrs') # Taken from Slack on 11 Oct 2020 frb191228 = frb.FRB(frbname, FRB_191228_coord, 298 * units.pc / units.cm ** 3, repeater=False) # First Slack posting # Error ellipse : Day+2021 frb191228.set_ee(0.34, 0.34, 0., 68.) frb191228.set_ee(0.830, 0.823, 0., 68., stat=False) # Error in DM frb191228.DM_err = 0.05 * units.pc / units.cm ** 3 # NE2001 frb191228.set_DMISM() # RM # frb190102.RM = 10 * units.rad / units.m**2 # frb190102.RM_err = 1 * units.rad / units.m**2 # References frb191228.refs = ['Bhandari2021', 'Day2021'] # Write path = resource_filename('frb', 'data/FRBs') frb191228.write_to_json(path=path)
def frb_20200120E(): """M81 + globular cluster """ frbname = 'FRB20200120E' FRB_20200120E_coord = SkyCoord("09h57m54.68s +68d49m08.0s", frame='icrs') # From EVN localisation frb20200120E = frb.FRB(frbname, FRB_20200120E_coord, 87.818 * units.pc / units.cm ** 3,z_frb=0.0008,repeater=True) # Error ellipse (Statistical) frb20200120E.set_ee(0.4,0.4, 0., 68.) # Error ellipse (Systematic) frb20200120E.set_ee(0., 0., 0., 68., stat=False) # Error in DM frb20200120E.DM_err = 0.007 * units.pc / units.cm ** 3 # RM frb20200120E.RM = -29.8 * units.rad / units.m**2 frb20200120E.RM_err = 0.5 * units.rad / units.m**2 # NE2001 frb20200120E.set_DMISM() # RM # frb190102.RM = 10 * units.rad / units.m**2 # frb190102.RM_err = 1 * units.rad / units.m**2 # References frb20200120E.refs = ['Bhardwaj2021','Kirsten2021'] # Write path = resource_filename('frb', 'data/FRBs') frb20200120E.write_to_json(path=path)
def frb_180301(): """ Bhandari+2021 """ frbname = 'FRB20180301' FRB_180301_coord = SkyCoord("06h12m54.44s +04d40m15.8s", frame='icrs') # Bhandari+2021 frb180301 = frb.FRB(frbname, FRB_180301_coord, 536 * units.pc / units.cm ** 3, z_frb=0.33044, repeater=True) # Slack posting # Error ellipse (Statistical) # frb180301.set_ee(0.01142,0.00825, 0., 68.) # Error ellipse (Systematic) frb180301.set_ee(0.619, 0.603, 0., 68., stat=False) # Error in DM frb180301.DM_err = 13 * units.pc / units.cm ** 3 # NE2001 frb180301.set_DMISM() # RM # frb190102.RM = 10 * units.rad / units.m**2 # References frb180301.refs = ['Bhandari2021'] # Write path = resource_filename('frb', 'data/FRBs') frb180301.write_to_json(path=path)
def frb_190714(): """ Day+2020 https://ui.adsabs.harvard.edu/abs/2020MNRAS.497.3335D/abstract Heintz+2020 Updated FRB localization + error from Day+2021 """ # Taken from Slack; Cherie posted on 19 June 2020 # Updated as per Day 2021 FRB_190714_coord = SkyCoord('J121555.13-130115.6', unit=(units.hourangle, units.deg)) frb190714 = frb.FRB('FRB20190714', FRB_190714_coord, 504.13 * units.pc / units.cm ** 3, z_frb=0.2365, repeater=False) # Error ellipse (Day+2021) frb190714.set_ee(0.17, 0.10, theta=90., cl=68.) # Statistical frb190714.set_ee(0.5191, 0.376, theta=90., cl=68., stat=False) # Systematic # Error in DM frb190714.DM_err = 0.1 * units.pc / units.cm ** 3 # NE2001 frb190714.set_DMISM() # Fluence frb190714.fluence = 12 * units.Jy * units.ms # From Cherie on Slack (2/10 - 2020) frb190714.fluence_err = 2 * units.Jy * units.ms # References frb190714.refs = ['Heintz2020', 'Day2021'] # Write path = resource_filename('frb', 'data/FRBs') frb190714.write_to_json(path=path)
def frb_181112(): """ Generate the JSON file for FRB 181112 All of the data comes from Prochaska+2019, Science, in press Returns: """ frb181112 = frb.FRB('FRB181112', 'J214923.63-525815.33', 589.27 * units.pc / units.cm**3, z_frb=0.4755) # Error in DM frb181112.DM_err = 0.03 * units.pc / units.cm**3 # Error ellipse frb181112.set_ee(a=555.30/1e3, b=152.93/1e3, theta=120.15, cl=68.) # RM frb181112.RM = 10.9 * units.rad / units.m**2 frb181112.RM_err = 0.9 * units.rad / units.m**2 # NE2001 frb181112.set_DMISM() # References frb181112.refs = ['Prochaska2019'] # Write path = resource_filename('frb', 'data/FRBs') frb181112.write_to_json(path=path)
def frb_190714(): """ Day+2020 https://ui.adsabs.harvard.edu/abs/2020MNRAS.497.3335D/abstract Heintz+2020 """ # Taken from Slack; Cherie posted on 19 June 2020 FRB_190714_coord = SkyCoord('J121555.12-130115.7', unit=(units.hourangle, units.deg)) frb190714 = frb.FRB('FRB190714', FRB_190714_coord, 504.13 * units.pc / units.cm**3, z_frb=0.2365, repeater=False) # Error ellipse frb190714.set_ee(0.16, 0.1, theta=90., cl=68.) # Statistical frb190714.set_ee(0.32, 0.2, theta=90., cl=68., stat=False) # Systematic # Error in DM frb190714.DM_err = 0.1 * units.pc / units.cm**3 # NE2001 frb190714.set_DMISM() # RM # References frb190714.refs = ['Heintz2020'] # Write path = resource_filename('frb', 'data/FRBs') frb190714.write_to_json(path=path)
def frb_200906(): """ Bhandari+2021 Updated FRB localization + error from Day+2021 """ frbname = 'FRB20200906' FRB_200906_coord = SkyCoord("03h33m59.08s -14d04m59.46s", frame='icrs') # Day+2021 frb200906 = frb.FRB(frbname, FRB_200906_coord, 577.84 * units.pc / units.cm**3, z_frb=0.36879, repeater=False) # Slack posting # Error ellipse (Statistical) frb200906.set_ee(0.11, 0.102, 0., 68.) # Error ellipse (Systematic) frb200906.set_ee(0.55, 0.340, 0., 68., stat=False) # Error in DM frb200906.DM_err = 0.02 * units.pc / units.cm ** 3 # NE2001 frb200906.set_DMISM() # RM # frb190102.RM = 10 * units.rad / units.m**2 # frb190102.RM_err = 1 * units.rad / units.m**2 # References frb200906.refs = ['Bhandari2021', 'Day2021'] # Write path = resource_filename('frb', 'data/FRBs') frb200906.write_to_json(path=path)
def frb_190608(): """Bhandari+20, ApJL, Day+20, in prep. -- Macquart al. Nature """ fname = 'FRB190608' frb190608 = frb.FRB( fname, "J221604.77-075353.7", # Pulled from Slack on 2020 Mar 18 339.8 * units.pc / units.cm**3, z_frb=0.1177805) # Taken from the SDSS table # Error ellipse [REQUIRED] frb190608.set_ee(0.19315, 0.18, 0., 68.) # Statistsical frb190608.set_ee(0.178292, 0.18, 0., 68., stat=False) # Systematic # Error in DM frb190608.DM_err = 1 * units.pc / units.cm**3 # NE2001 frb190608.set_DMISM() # RM #frb190102.RM = 10 * units.rad / units.m**2 #frb190102.RM_err = 1 * units.rad / units.m**2 # References frb190608.refs = ['Bhandari2020', 'Day2020'] # Write path = resource_filename('frb', 'data/FRBs') frb190608.write_to_json(path=path)
def frb_180924(): """ FRB 180924 All of the data currently comes from Bannister et al. 2019 https://ui.adsabs.harvard.edu/abs/2019Sci...365..565B/abstract """ frb180924 = frb.FRB('FRB180924', 'J214425.26-405400.1', 361.4*units.pc / units.cm**3, z_frb=0.3212) # Error in DM (Bannister 2019) frb180924.DM_err = 0.06 * units.pc / units.cm**3 # NE2001 frb180924.set_DMISM() # Bannister 2019 frb180924.fluence = 16 * units.Jy * units.ms frb180924.fluence_err = 1 * units.Jy * units.ms frb180924.RM = 14 * units.rad / units.m**2 frb180924.RM_err = 1 * units.rad / units.m**2 frb180924.lpol = 80. # % frb180924.lpol_err = 10. # Error ellipse frb180924.set_ee(a=100./1e3, b=100./1e3, theta=0., cl=68.) # References frb180924.refs = ['Bannister2019'] # Write path = resource_filename('frb', 'data/FRBs') frb180924.write_to_json(path=path)
def frb_190608(): """Bhandari+20, ApJL -- Macquart al. Nature """ fname = 'FRB190608' frb190608 = frb.FRB(fname, "J221604.74-075353.6", # Pulled from Slack on 2019 Sep 20 339.8 * units.pc / units.cm**3, z_frb=0.1177805) # Taken from the SDSS table # Error ellipse [REQUIRED] frb190608.set_ee(0.1, 0.1, 0., 95.) # Error in DM frb190608.DM_err = 1 * units.pc / units.cm**3 # NE2001 frb190608.set_DMISM() # RM #frb190102.RM = 10 * units.rad / units.m**2 #frb190102.RM_err = 1 * units.rad / units.m**2 # References frb190608.refs = ['Bhandari2019'] # Write path = resource_filename('frb', 'data/FRBs') frb190608.write_to_json(path=path)
def frb_190614(): """ Law+2020 https://ui.adsabs.harvard.edu/abs/2020ApJ...899..161L/abstract """ fname = 'FRB190614' FRB_190614_coord = SkyCoord(ra=65.07552, dec=73.70674, unit='deg') frb190614 = frb.FRB(fname, FRB_190614_coord, 959 * units.pc / units.cm**3, repeater=False) # Error ellipse [REQUIRED] frb190614.set_ee(a=0.8, b=0.4, theta=67., cl=68.) # Error in DM frb190614.DM_err = 1 * units.pc / units.cm**3 # NE2001 frb190614.set_DMISM() # RM # frb190102.RM = 10 * units.rad / units.m**2 # frb190102.RM_err = 1 * units.rad / units.m**2 # References frb190614.refs = ['Law2020'] # Write path = resource_filename('frb', 'data/FRBs') frb190614.write_to_json(path=path)
def main(pargs): """ Run """ import warnings import numpy as np from matplotlib import pyplot as plt from astropy.io import fits from astropy import units from frb import frb from frb.figures import galaxies as ffgal from frb.figures import utils as ffutils from linetools.scripts.utils import coord_arg_to_coord # Load up hdu = fits.open(pargs.fits_file) icoord = coord_arg_to_coord(pargs.frb_coord) # Parse if pargs.vmnx is not None: tstr = pargs.vmnx.split(',') vmnx = (float(tstr[0]), float(tstr[1])) else: vmnx = (None, None) # Dummy FRB object FRB = frb.FRB('TMP', icoord, 0.) FRB.set_ee(1.0, 1.0, 0., 95.) fig = plt.figure(figsize=(7, 7)) ffutils.set_mplrc() ffgal.sub_image(fig, hdu, FRB, vmnx=vmnx, cmap='gist_heat', frb_clr='white', imsize=pargs.imsize) #img_center=HG190608.coord, # Layout and save plt.tight_layout(pad=0.2, h_pad=0.1, w_pad=0.1) plt.savefig(pargs.outfile, dpi=300) plt.close() print('Wrote {:s}'.format(pargs.outfile))
def frb_180924(): """ FRB 180924 All of the data currently comes from Bannister et al. 2019 https://ui.adsabs.harvard.edu/abs/2019Sci...365..565B/abstract Now including updated values on FRB position, DM, pulse width and scattering time from Day et al. 2020 Uncertinties in FRB localization updated by Day et al. 2021 """ frb180924 = frb.FRB('FRB20180924', 'J214425.255-405400.10', 362.16*units.pc / units.cm**3, z_frb=0.3212, repeater=False) # Error ellipse frb180924.set_ee(a=0.07, b=0.07, theta=0., cl=68.) # Statistical (Day+2020) frb180924.set_ee(a=0.1611, b=0.1611, theta=0., cl=68.,stat=False) # Systematic (Day+2021) # Error in DM (Day 2019) frb180924.DM_err = 0.01 * units.pc / units.cm**3 # NE2001 frb180924.set_DMISM() # RM, fluence and polarization frb180924.fluence = 16 * units.Jy * units.ms # From Bhandari+20 frb180924.fluence_err = 1 * units.Jy * units.ms # -- //-- frb180924.RM = 22 * units.rad / units.m**2 frb180924.RM_err = 2 * units.rad / units.m**2 frb180924.lpol = 80. # % frb180924.lpol_err = 10. # Pulse properties frb180924.set_pulse(1.2725*units.GHz, Wi=0.09*units.ms, Wi_err=0.04*units.ms, tscatt=0.68*units.ms, tscatt_err=0.03*units.ms) # References frb180924.refs = ['Bannister2019', 'Day2020', 'Day2021'] # Write path = resource_filename('frb', 'data/FRBs') frb180924.write_to_json(path=path)
def frb_201124(): """ ATELs only so far """ # ATEL 14603 (VLBI) FRB_201124_coord = SkyCoord("05h08m03.5077s 26d03m38.504s", frame='icrs') frb201124 = frb.FRB('FRB20201124', FRB_201124_coord, 411. * units.pc / units.cm ** 3, z_frb=0.0982, repeater=True) # Error ellipse [REQUIRED] frb201124.set_ee(0.004, 0.004, theta=0., cl=68.) # ATEL # Error in DM #frb191001.DM_err = 0.07 * units.pc / units.cm ** 3 # NE2001 frb201124.set_DMISM() # RM (Kumar+21) frb201124.RM = -613 * units.rad / units.m**2 frb201124.RM_err = 2 * units.rad / units.m**2 # RM and fluence (Bhandari+20b) #frb201124.RM = 55.5 * units.rad / units.m**2 #frb201124.RM_err = 0.9 * units.rad / units.m**2 #frb201124.fluence = 143 * units.Jy * units.ms #frb201124.fluence_err = 15 * units.Jy * units.ms # Pulse properties #frb191001.set_pulse(0.920*units.GHz, # Wi=0.22*units.ms, # Wi_err=0.03*units.ms, # tscatt=3.3*units.ms, # tscatt_err=0.2*units.ms) # References #frb191001.refs = ['Bhandari2020b'] # Write path = resource_filename('frb', 'data/FRBs') frb201124.write_to_json(path=path)
def frb_121102(): """ FRB 121102 All of the data currently comes from Tendulkar et al. 2017 https://ui.adsabs.harvard.edu/abs/2017ApJ...834L...7T/abstract """ frb121102 = frb.FRB('FRB121102', 'J053158.7+330852.5', 558.1*units.pc/units.cm**3, z_frb=0.19273) # NE2001 frb121102.set_DMISM() # Error ellipse frb121102.set_ee(0.1, 0.1, 0., 95.) # References frb121102.refs = ['Tendulkar2017'] # Write frb121102.write_to_json() # Test frb121102.from_json('FRB121102.json')
def frb_121102(): """ FRB 121102 All of the data currently comes from Tendulkar et al. 2017 https://ui.adsabs.harvard.edu/abs/2017ApJ...834L...7T/abstract """ frb121102 = frb.FRB('FRB121102', 'J053158.7+330852.5', 558.1 * units.pc / units.cm**3, z_frb=0.19273, repeater=True) # NE2001 frb121102.set_DMISM() # Error ellipse frb121102.set_ee(0.1, 0.1, theta=0., cl=95.) # References frb121102.refs = ['Tendulkar2017'] # Write path = resource_filename('frb', 'data/FRBs') frb121102.write_to_json(path=path)
def frb_181112(): """ Generate the JSON file for FRB 181112 All of the data comes from Prochaska+2019, Science, in press Returns: FRB position and localization updated by Day+2021 """ frb181112 = frb.FRB('FRB20181112', 'J214923.63-525815.4', 589.27 * units.pc / units.cm**3, z_frb=0.4755, repeater=False) # Error in DM frb181112.DM_err = 0.03 * units.pc / units.cm**3 # Error ellipse -- Updated by Day+2021 frb181112.set_ee(a=555.30/1e3, b=152.93/1e3, theta=120.15, cl=68.) # Statistical (Prochaska+2019) frb181112.set_ee(a=3.2*1.79, b=0.8*1.79, theta=120.15, cl=68.,stat=False) # Systematic # RM and fluence frb181112.RM = 10.5 * units.rad / units.m**2 frb181112.RM_err = 0.4 * units.rad / units.m**2 frb181112.fluence = 20.2 * units.Jy * units.ms frb181112.fluence_err = 0.1 * units.Jy * units.ms # Pulse properties frb181112.set_pulse(1.2725*units.GHz, Wi=0.016*units.ms, Wi_err=0.001*units.ms, tscatt=0.021*units.ms, tscatt_err=0.001*units.ms) # NE2001 frb181112.set_DMISM() # References frb181112.refs = ['Prochaska2019', 'Day2021'] # Write path = resource_filename('frb', 'data/FRBs') frb181112.write_to_json(path=path)
def frb_190611(): """ Macquart al. Nature Day et al. 2020 https://ui.adsabs.harvard.edu/abs/2020MNRAS.497.3335D/abstract Update systematic uncertainty in FRB localization from Day+2021 """ FRB_190611_coord = SkyCoord('J212258.94-792351.3', # Day+2021 unit=(units.hourangle, units.deg)) frb190611 = frb.FRB('FRB20190611', FRB_190611_coord, 332.63 * units.pc / units.cm**3, z_frb=0.3778, repeater=False) # Bright # Error ellipse frb190611.set_ee(0.34, 0.32, theta=0., cl=68.) # Statistical (Day+2020) frb190611.set_ee(1.12, 1.07, theta=0., cl=68., stat=False) # Systematic (Day+2021) # Error in DM frb190611.DM_err = 0.04 * units.pc / units.cm**3 # NE2001 frb190611.set_DMISM() # RM and fluence frb190611.RM = 20 * units.rad / units.m**2 # From Day+20 frb190611.RM_err = 4 * units.rad / units.m**2 frb190611.fluence = 10 * units.Jy * units.ms # From Macquart+20 frb190611.fluence_err = 2 * units.Jy * units.ms # Pulse properties frb190611.set_pulse(1.2725*units.GHz, Wi=0.09*units.ms, Wi_err=0.02*units.ms, tscatt=0.18*units.ms, tscatt_err=0.02*units.ms) # References frb190611.refs = ['MacQuart2020', 'Day2020', 'Day2021'] # Write path = resource_filename('frb', 'data/FRBs') frb190611.write_to_json(path=path)
def frb_190102(): """Bhandari+20, ApJL -- Now including updated values on FRB position, DM, pulse width and scattering time from Day et al. 2020 Update systematic uncertainty in FRB localization from Day+2021 """ fname = 'FRB20190102' #wv_oiii = 6466.48 #z_OIII = wv_oiii / 5008.239 - 1 frb190102 = frb.FRB(fname, 'J212939.76-792832.5', # Day+20 364.545 * units.pc / units.cm**3, z_frb=0.2912, repeater=False) # Updated redshift # Error ellipse [REQUIRED] frb190102.set_ee(0.21, 0.17, theta=0., cl=68.) # Statistical (Day+2020) frb190102.set_ee(0.936, 0.7876, theta=0., cl=68., stat=False) # Systematic (Day+2021) # Error in DM frb190102.DM_err = 0.004 * units.pc / units.cm**3 # NE2001 frb190102.set_DMISM() # RM and fluence frb190102.fluence = 14 * units.Jy * units.ms # From Bhandari+20 frb190102.fluence_err = 1 * units.Jy * units.ms # -- //-- frb190102.RM = -105 * units.rad / units.m**2 frb190102.RM_err = 1 * units.rad / units.m**2 # Pulse properties frb190102.set_pulse(1.2725*units.GHz, Wi=0.053*units.ms, Wi_err=0.002*units.ms, tscatt=0.041*units.ms, tscatt_err=0.003*units.ms) # References frb190102.refs = ['Bhandari2020', 'Day2020', 'Day2021'] # Write path = resource_filename('frb', 'data/FRBs') frb190102.write_to_json(path=path)