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()