Check B1 law
Contents
Check B1 law#
With this notebook, you can check whether a source is in B1 law or not. The code uses the source name to search its coordinates in the SIMBAD catalog. If the source is not found, the user can enter the coordinates manually.
import os
if not os.path.isfile("B1_law.txt"):
os.system("wget https://raw.githubusercontent.com/fcangemi/gp-tools-svom/main/B1_law.txt")
!pip install astroquery
Requirement already satisfied: astroquery in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (0.4.6)
Requirement already satisfied: numpy>=1.16 in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from astroquery) (1.21.5)
Requirement already satisfied: astropy>=4.0 in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from astroquery) (5.1)
Requirement already satisfied: requests>=2.4.3 in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from astroquery) (2.28.1)
Requirement already satisfied: beautifulsoup4>=4.3.2 in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from astroquery) (4.11.1)
Requirement already satisfied: html5lib>=0.999 in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from astroquery) (1.1)
Requirement already satisfied: keyring>=4.0 in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from astroquery) (23.4.0)
Requirement already satisfied: pyvo>=1.1 in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from astroquery) (1.4)
Requirement already satisfied: PyYAML>=3.13 in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from astropy>=4.0->astroquery) (6.0)
Requirement already satisfied: pyerfa>=2.0 in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from astropy>=4.0->astroquery) (2.0.0)
Requirement already satisfied: packaging>=19.0 in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from astropy>=4.0->astroquery) (21.3)
Requirement already satisfied: soupsieve>1.2 in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from beautifulsoup4>=4.3.2->astroquery) (2.3.1)
Requirement already satisfied: six>=1.9 in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from html5lib>=0.999->astroquery) (1.16.0)
Requirement already satisfied: webencodings in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from html5lib>=0.999->astroquery) (0.5.1)
Requirement already satisfied: importlib-metadata>=3.6 in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from keyring>=4.0->astroquery) (4.11.3)
Requirement already satisfied: certifi>=2017.4.17 in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from requests>=2.4.3->astroquery) (2022.9.24)
Requirement already satisfied: urllib3<1.27,>=1.21.1 in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from requests>=2.4.3->astroquery) (1.26.11)
Requirement already satisfied: charset-normalizer<3,>=2 in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from requests>=2.4.3->astroquery) (2.0.4)
Requirement already satisfied: idna<4,>=2.5 in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from requests>=2.4.3->astroquery) (3.3)
Requirement already satisfied: zipp>=0.5 in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from importlib-metadata>=3.6->keyring>=4.0->astroquery) (3.8.0)
Requirement already satisfied: pyparsing!=3.0.5,>=2.0.2 in /Users/fcangemi/opt/anaconda3/lib/python3.9/site-packages (from packaging>=19.0->astropy>=4.0->astroquery) (3.0.9)
from astroquery.simbad import Simbad
from astroquery.vizier import Vizier
from astroquery.ipac.ned import Ned
from astropy.coordinates import SkyCoord, Galactic
import astropy.units as u
import matplotlib.pyplot as plt
import numpy as np
import warnings
warnings.filterwarnings('ignore')
def which_source(source):
result_table = Simbad.query_object(source)
if(result_table == None):
print("Unable to find", source.rstrip("\n"), ", please enter the source coordinates:")
ra = float(input("ra (degrees):"))
dec = float(input("dec (degrees):"))
else:
ra = result_table["RA"]
dec = result_table["DEC"]
return source, ra, dec
def read_B1():
b1_law = np.genfromtxt("B1_law.txt", unpack = True)
b1_ras = b1_law[1]
b1_decs = b1_law[2]
return b1_ras, b1_decs
def check_inB1_name(source):
source, ra_source, dec_source = which_source(source)
b1_ras, b1_decs = read_B1()
c_source = SkyCoord(ra_source, dec_source, frame = "icrs", unit = (u.hourangle, u.deg))
count = 0
all_sep = []
for b1_ra, b1_dec in zip(b1_ras, b1_decs):
c_b1 = SkyCoord(b1_ra, b1_dec, frame = "icrs", unit = "deg")
sep = c_source.separation(c_b1)
all_sep.append(sep.value)
if sep.value <= 10:
count += 1
else:
continue
if count == 0:
print(source.rstrip('\n'), "is not in B1 law, minimal separation =", min(all_sep)*u.degree)
else:
print(source.rstrip('\n'), "is in B1 law")
def check_inB1_list(list_of_sources):
#list_of_sources = open(list_of_sources_file, "r")
for source in list_of_sources:
check_inB1_name(source)
Using the source name#
You can directly write the source name. Click on the rocket at the top of this page, and then click on the "Live Code" button to edit the cell below.
Here an example for Cygnus X-1; write your source name and then click on "run" (the first time you are running, you may need to click on "restart & run all").
check_inB1_name("Cygnus X-1")
Cygnus X-1 is not in B1 law, minimal separation = [62.11449208] deg
Using a list of sources#
Alternatively, you can provide a list of sources:
list_of_sources = ["Mrk 501",
"Mrk 421",
"1ES 1959+650",
"1ES 2344+514",
"M87",
"PG 1553+113",
"NGC 1246",
"IC 310",
"1ES 1011+496",
"1ES 1215+303"
]
check_inB1_list(list_of_sources)
Mrk 501 is not in B1 law, minimal separation = [27.62596468] deg
Mrk 421 is in B1 law
1ES 1959+650 is not in B1 law, minimal separation = [58.38624248] deg
1ES 2344+514 is not in B1 law, minimal separation = [47.73709681] deg
M87 is in B1 law
PG 1553+113 is not in B1 law, minimal separation = [21.05383455] deg
NGC 1246 is not in B1 law, minimal separation = [37.30023443] deg
IC 310 is not in B1 law, minimal separation = [48.45962333] deg
1ES 1011+496 is not in B1 law, minimal separation = [17.64474591] deg
1ES 1215+303 is not in B1 law, minimal separation = [14.47156307] deg