diff options
Diffstat (limited to 'deproject_sbp.py')
-rwxr-xr-x | deproject_sbp.py | 37 |
1 files changed, 2 insertions, 35 deletions
diff --git a/deproject_sbp.py b/deproject_sbp.py index 5dfbc2d..ce6ad8d 100755 --- a/deproject_sbp.py +++ b/deproject_sbp.py @@ -6,6 +6,7 @@ # # Change logs: # 2016-06-24: +# * Move class 'ChandraPixel' to module 'astro_params.py' # * Split class 'Projection' to a separate module 'projection.py' # * Move class 'DensityProfile' to tool 'calc_mass_potential.py' # * Split class 'AstroParams' to separate module 'astro_params.py' @@ -148,7 +149,6 @@ import json from collections import OrderedDict import astropy.units as au -from astropy.cosmology import FlatLambdaCDM import numpy as np import pandas as pd import scipy.optimize @@ -159,7 +159,7 @@ from matplotlib.backends.backend_agg import FigureCanvasAgg as FigureCanvas from matplotlib.figure import Figure from configobj import ConfigObj -from astro_params import AstroParams +from astro_params import AstroParams, ChandraPixel from projection import Projection plt.style.use("ggplot") @@ -809,39 +809,6 @@ class SBP: ###################################################################### -class ChandraPixel: - """ - Chandra pixel unit conversions. - """ - angle = 0.492 * au.arcsec - z = None - # cosmology calculator - cosmo = None - # angular diameter distance - D_A = None - # length of one pixel at the given redshift - length = None - - def __init__(self, z=None): - self.z = z - self.cosmo = FlatLambdaCDM(H0=AstroParams.H0, - Om0=AstroParams.OmegaM0) - if z is not None: - self.D_A = self.cosmo.angular_diameter_distance(z) - self.length = self.D_A * self.angle.to(au.radian).value - - def get_angle(self): - return self.angle - - def get_length(self, z=None): - if z is None: - length = self.length - else: - D_A = self.cosmo.angular_diameter_distance(z) - length = D_A * self.angle.to(au.radian).value - return length - - class BrightnessProfile: """ Calculate the electron number density and/or gas mass density profile |