In [1]:
%matplotlib inline
import numpy as np
import dmsky.roster
import dmsky.skymap
/Users/alex/fermi/software/dmsky/dmsky/data/targets
In [2]:
# Build the library of pre-defined rosters
library = dmsky.roster.RosterLibrary()
/Users/alex/fermi/software/dmsky/dmsky/data/rosters
In [5]:
# The roster name often follows from a literature reference
name = 'martinez2015_nfw'
# To get a list of the available rosters
print library.library.keys()
['test2015', 'dsphs_gc', 'ackermann2015_dsphs', 'martinez2015_nfw', 'ackermann2013_gc']
In [6]:
# Build the roster of targets
dsphs = library.create_roster(name)
In [7]:
# Fill the skymap with the jfactors from the target list
skymap = dmsky.skymap.Skymap(dsphs, nside=128)
skymap.fill()
Dwarf
name : bootes_I
ra : 210.025
dec : 14.5
distance : 66.0
density : <class 'dmsky.density.NFWProfile'>
Parameters:
rs : Parameter(8.331336e+20, [None, None], [None, None], False)
rhos : Parameter(1.78718119032e-23, [None, None], [None, None], False)
rmin : Parameter(0.0, [None, None], [None, None], False)
rmax : Parameter(3.08568e+21, [None, None], [None, None], False)
rhomax : Parameter(inf, [None, None], [None, None], False)
covar : None
Dwarf
name : carina
ra : 100.4029
dec : -50.966
distance : 105.0
density : <class 'dmsky.density.NFWProfile'>
Parameters:
rs : Parameter(6.479928e+20, [None, None], [None, None], False)
rhos : Parameter(2.05796621916e-23, [None, None], [None, None], False)
rmin : Parameter(0.0, [None, None], [None, None], False)
rmax : Parameter(3.08568e+21, [None, None], [None, None], False)
rhomax : Parameter(inf, [None, None], [None, None], False)
covar : None
Dwarf
name : canes_venatici_II
ra : 194.29167
dec : 34.32083
distance : 160.0
density : <class 'dmsky.density.NFWProfile'>
Parameters:
rs : Parameter(5.862792e+20, [None, None], [None, None], False)
rhos : Parameter(2.47091338813e-23, [None, None], [None, None], False)
rmin : Parameter(0.0, [None, None], [None, None], False)
rmax : Parameter(3.08568e+21, [None, None], [None, None], False)
rhomax : Parameter(inf, [None, None], [None, None], False)
covar : None
Dwarf
name : coma_berenices
ra : 186.74583
dec : 23.90417
distance : 44.0
density : <class 'dmsky.density.NFWProfile'>
Parameters:
rs : Parameter(5.554224e+20, [None, None], [None, None], False)
rhos : Parameter(2.58599702539e-23, [None, None], [None, None], False)
rmin : Parameter(0.0, [None, None], [None, None], False)
rmax : Parameter(3.08568e+21, [None, None], [None, None], False)
rhomax : Parameter(inf, [None, None], [None, None], False)
covar : None
Dwarf
name : draco
ra : 260.05167
dec : 57.91528
distance : 76.0
density : <class 'dmsky.density.NFWProfile'>
Parameters:
rs : Parameter(1.079988e+21, [None, None], [None, None], False)
rhos : Parameter(1.55701391581e-23, [None, None], [None, None], False)
rmin : Parameter(0.0, [None, None], [None, None], False)
rmax : Parameter(3.08568e+21, [None, None], [None, None], False)
rhomax : Parameter(inf, [None, None], [None, None], False)
covar : None
Dwarf
name : fornax
ra : 39.997
dec : -34.449
distance : 147.0
density : <class 'dmsky.density.NFWProfile'>
Parameters:
rs : Parameter(1.4502696e+21, [None, None], [None, None], False)
rhos : Parameter(9.0036022088e-24, [None, None], [None, None], False)
rmin : Parameter(0.0, [None, None], [None, None], False)
rmax : Parameter(3.08568e+21, [None, None], [None, None], False)
rhomax : Parameter(inf, [None, None], [None, None], False)
covar : None
Dwarf
name : hercules
ra : 247.75833
dec : 12.79167
distance : 132.0
density : <class 'dmsky.density.NFWProfile'>
Parameters:
rs : Parameter(6.17136e+20, [None, None], [None, None], False)
rhos : Parameter(2.47091338813e-23, [None, None], [None, None], False)
rmin : Parameter(0.0, [None, None], [None, None], False)
rmax : Parameter(3.08568e+21, [None, None], [None, None], False)
rhomax : Parameter(inf, [None, None], [None, None], False)
covar : None
Dwarf
name : leo_II
ra : 168.37
dec : 22.15167
distance : 233.0
density : <class 'dmsky.density.NFWProfile'>
Parameters:
rs : Parameter(9.565608e+20, [None, None], [None, None], False)
rhos : Parameter(1.23884150693e-23, [None, None], [None, None], False)
rmin : Parameter(0.0, [None, None], [None, None], False)
rmax : Parameter(3.08568e+21, [None, None], [None, None], False)
rhomax : Parameter(inf, [None, None], [None, None], False)
covar : None
Dwarf
name : leo_IV
ra : 173.2375
dec : -0.53333
distance : 154.0
density : <class 'dmsky.density.NFWProfile'>
Parameters:
rs : Parameter(6.17136e+20, [None, None], [None, None], False)
rhos : Parameter(2.47091338813e-23, [None, None], [None, None], False)
rmin : Parameter(0.0, [None, None], [None, None], False)
rmax : Parameter(3.08568e+21, [None, None], [None, None], False)
rhomax : Parameter(inf, [None, None], [None, None], False)
covar : None
Dwarf
name : sculptor
ra : 15.039
dec : -33.709
distance : 86.0
density : <class 'dmsky.density.NFWProfile'>
Parameters:
rs : Parameter(1.2034152e+21, [None, None], [None, None], False)
rhos : Parameter(1.13052749539e-23, [None, None], [None, None], False)
rmin : Parameter(0.0, [None, None], [None, None], False)
rmax : Parameter(3.08568e+21, [None, None], [None, None], False)
rhomax : Parameter(inf, [None, None], [None, None], False)
covar : None
Dwarf
name : segue_1
ra : 151.76667
dec : 16.08194
distance : 23.0
density : <class 'dmsky.density.NFWProfile'>
Parameters:
rs : Parameter(4.62852e+20, [None, None], [None, None], False)
rhos : Parameter(3.74360302366e-23, [None, None], [None, None], False)
rmin : Parameter(0.0, [None, None], [None, None], False)
rmax : Parameter(3.08568e+21, [None, None], [None, None], False)
rhomax : Parameter(inf, [None, None], [None, None], False)
covar : None
Dwarf
name : segue_2
ra : 34.81667
dec : 20.17528
distance : 35.0
density : <class 'dmsky.density.NFWProfile'>
Parameters:
rs : Parameter(4.937088e+20, [None, None], [None, None], False)
rhos : Parameter(2.97186569148e-23, [None, None], [None, None], False)
rmin : Parameter(0.0, [None, None], [None, None], False)
rmax : Parameter(3.08568e+21, [None, None], [None, None], False)
rhomax : Parameter(inf, [None, None], [None, None], False)
covar : None
Dwarf
name : sextans
ra : 153.2625
dec : -1.61472
distance : 86.0
density : <class 'dmsky.density.NFWProfile'>
Parameters:
rs : Parameter(6.17136e+20, [None, None], [None, None], False)
rhos : Parameter(2.58599702539e-23, [None, None], [None, None], False)
rmin : Parameter(0.0, [None, None], [None, None], False)
rmax : Parameter(3.08568e+21, [None, None], [None, None], False)
rhomax : Parameter(inf, [None, None], [None, None], False)
covar : None
Dwarf
name : ursa_major_II
ra : 132.875
dec : 63.13
distance : 32.0
density : <class 'dmsky.density.NFWProfile'>
Parameters:
rs : Parameter(5.862792e+20, [None, None], [None, None], False)
rhos : Parameter(2.83647317706e-23, [None, None], [None, None], False)
rmin : Parameter(0.0, [None, None], [None, None], False)
rmax : Parameter(3.08568e+21, [None, None], [None, None], False)
rhomax : Parameter(inf, [None, None], [None, None], False)
covar : None
Dwarf
name : ursa_minor
ra : 227.2854
dec : 67.2225
distance : 76.0
density : <class 'dmsky.density.NFWProfile'>
Parameters:
rs : Parameter(1.4502696e+21, [None, None], [None, None], False)
rhos : Parameter(9.81595729531e-24, [None, None], [None, None], False)
rmin : Parameter(0.0, [None, None], [None, None], False)
rmax : Parameter(3.08568e+21, [None, None], [None, None], False)
rhomax : Parameter(inf, [None, None], [None, None], False)
covar : None
Dwarf
name : willman_1
ra : 162.3375
dec : 51.05
distance : 38.0
density : <class 'dmsky.density.NFWProfile'>
Parameters:
rs : Parameter(4.937088e+20, [None, None], [None, None], False)
rhos : Parameter(3.11402783162e-23, [None, None], [None, None], False)
rmin : Parameter(0.0, [None, None], [None, None], False)
rmax : Parameter(3.08568e+21, [None, None], [None, None], False)
rhomax : Parameter(inf, [None, None], [None, None], False)
covar : None
/Users/alex/fermi/software/dmsky/dmsky/jcalc.py:380: RuntimeWarning: divide by zero encountered in log10
return np.log10([_dhalo, _psi, _jval])
In [20]:
# Some plotting (not the best projection)
import healpy
healpy.mollview(np.log10(skymap.values))
/usr/local/lib/python2.7/site-packages/IPython/kernel/__main__.py:3: RuntimeWarning: divide by zero encountered in log10
app.launch_new_instance()
Content source: kadrlica/dmsky
Similar notebooks: