from deproject import Deproject
import logging
logger = logging.getLogger("sherpa")
logging.basicConfig(level=logging.INFO, filename="deproject.log", filemode="w")
import numpy
redshift = 0.004233 # M87 redshift
arcsec_per_pixel = 0.492 # ACIS plate scale
angdist = 4.9e25 # M87 distance (cm) (16 Mpc)
radii = numpy.arange(30., 640., 30) * arcsec_per_pixel
dep = Deproject(radii, theta=75, angdist=angdist)
# Load datasets for each annulus
for annulus in range(len(radii)-1):
dep.load_pha('m87/r%dgrspec.pha' % (annulus+1), annulus)
# Subtract background
dep.subtract()
# Set source model and ignore specified energy ranges
dep.set_source('xswabs*xsmekal')
dep.ignore(None, 0.5)
dep.ignore(1.8, 2.2)
dep.ignore(7, None)
# Specify Galactic absorption
dep.set_par('xswabs.nh', 0.0255)
dep.freeze("xswabs.nh")
from deproject import Deproject
from cosmocalc import cosmocalc
redshift = 0.004233
from math import pi
dep = Deproject(numpy.arange(30., 640., 30)*0.492)
dep.theta = 75
dep.angdist = cosmocalc(redshift)['DA_cm'] * 0.892
set_method("levmar")
set_stat("chi2gehrels")
for ann in range(dep.nshell):
dep.load_pha('m87/r%dgrspec.pha' % (ann+1), annulus=ann)
dep.set_source('xswabs*xsmekal')
dep.ignore(None, 0.5)
dep.ignore(1.8, 2.2)
dep.ignore(7, None)
dep.set_par('xswabs.nh', 0.0255)
dep.freeze("xswabs.nh")
dep.set_par('xsmekal.abundanc', 0.5)
dep.thaw('xsmekal.abundanc')
dep.set_par('xsmekal.redshift', redshift)
dep.subtract()
dep.fit()
