# Path: dart_board/constants.py
# G = 6.674e-8 # Gravitational constant in cgs
# GGG = 1.909e5 # Gravitational constant in Rsun * (km/s)^2 / Msun
# R_NS = 10.0  # NS radius in km


import sys
import numpy as np
import matplotlib.pyplot as plt
import scipy.optimize as so
    import mpl_toolkits.axisartist.angle_helper as angle_helper
    import matplotlib.gridspec as gridspec
    import mpl_toolkits.axisartist.angle_helper as angle_helper
    import matplotlib.gridspec as gridspec
from dart_board import constants as c
from matplotlib import colors
    from matplotlib.projections import PolarAxes
    from matplotlib.transforms import Affine2D
    from mpl_toolkits.axisartist import SubplotHost
    from mpl_toolkits.axisartist import GridHelperCurveLinear
    from matplotlib.projections import PolarAxes
    from matplotlib.transforms import Affine2D
    from mpl_toolkits.axisartist import SubplotHost
    from mpl_toolkits.axisartist import GridHelperCurveLinear





def get_theta_proj_degree(ra, dec, ra_b, dec_b):
    """ Return angular distance between two points

    Parameters
    ----------
    ra : float64
        Right ascension of first coordinate (degrees)
    dec : float64
        Declination of first coordinate (degrees)
    ra_b : float64
        Right ascension of second coordinate (degrees)
    dec_b : float64
        Declination of second coordinate (degrees)

    Returns
    -------
    theta : float64
        Angular distance (radians)
    """


# Path: dart_board/constants.py
# G = 6.674e-8 # Gravitational constant in cgs
# GGG = 1.909e5 # Gravitational constant in Rsun * (km/s)^2 / Msun
# R_NS = 10.0  # NS radius in km


# Path: dart_board/sf_history/sf_plotting.py
import sys
import numpy as np
import matplotlib.pyplot as plt
import scipy.optimize as so
    import mpl_toolkits.axisartist.angle_helper as angle_helper
    import matplotlib.gridspec as gridspec
    import mpl_toolkits.axisartist.angle_helper as angle_helper
    import matplotlib.gridspec as gridspec
from dart_board import constants as c
from matplotlib import colors
    from matplotlib.projections import PolarAxes
    from matplotlib.transforms import Affine2D
    from mpl_toolkits.axisartist import SubplotHost
    from mpl_toolkits.axisartist import GridHelperCurveLinear
    from matplotlib.projections import PolarAxes
    from matplotlib.transforms import Affine2D
    from mpl_toolkits.axisartist import SubplotHost
    from mpl_toolkits.axisartist import GridHelperCurveLinear






def get_theta_proj_degree(ra, dec, ra_b, dec_b):
    """ Return angular distance between two points

    Parameters
    ----------
    ra : float64
        Right ascension of first coordinate (degrees)
    dec : float64
        Declination of first coordinate (degrees)
    ra_b : float64
        Right ascension of second coordinate (degrees)
    dec_b : float64
        Declination of second coordinate (degrees)

    Returns
    -------
    theta : float64
        Angular distance (radians)
    """


    ra1 = c.deg_to_rad * ra