def __init__(self, **kwargs):
""" C'tor
Keyword arguments
-----------------
basedir : str
Top level directory for finding files
"""
self._name_factory = NameFactory(**kwargs)
self._catalog_comp_info_dicts = {}
self._split_comp_info_dicts = {}
def __init__(self, **kwargs):
""" C'tor
Keyword arguments
-----------------
name_policy : str
Name of yaml file contain file naming policy definitions
basedir : str
Top level directory for finding files
"""
self._name_factory = NameFactory(basedir=kwargs.get('basedir'))
self._diffuse_comp_info_dict = {}
def __init__(self, **kwargs):
""" C'tor
Keyword arguments
-----------------
projtype : str [healpix]
Used to define path to gasmap files
basedir : str
Top level directory for finding files
"""
self.projtype = kwargs.get('projtype', 'healpix')
self._name_factory = NameFactory(basedir=kwargs.get('basedir'))
self._ring_dicts = {}
self._diffuse_comp_info_dicts = {}
def __init__(self, **kwargs):
"""
Keyword arguments
-------------------
basedir : str
Top level directory for finding files
"""
self._name_factory = NameFactory(**kwargs)
self._dmm = kwargs.get('DiffuseModelManager',
DiffuseModelManager(**kwargs))
self._gmm = kwargs.get('GalpropMapManager',
GalpropMapManager(**kwargs))
self._csm = kwargs.get('CatalogSourceManager',
CatalogSourceManager(**kwargs))
self._library = {}
self._models = {}
self._spec_lib = SpectralLibrary({})
def test_name_factory():
names = dict(data_pass='******',
data_ver='P302',
evclass='source',
data_time='8years',
zcut='zmax105',
psftype='PSF3',
mktime='none',
ebin='E3',
coordsys='GAL',
irf_ver='V6',
sourcekey='ptsrc')
name_fact = NameFactory(**names)
test_dict = name_fact.make_filenames(**names)
test_irfs = name_fact.irfs(**names)
assert (test_dict['ltcube'] == 'lt_cubes/ltcube_8years_none_zmax105.fits')
assert (test_dict['ft1file'] == 'P8_P302_8years_source_zmax105.lst')
assert (
test_dict['bexpcube'] ==
'bexp_cubes/bexcube_P8_P302_8years_source_none_zmax105_E3_PSF3_GAL_V6.fits'
)
assert (
test_dict['srcmaps'] ==
'srcmaps/srcmaps_ptsrc_P8_P302_8years_source_none_zmax105_E3_PSF3_GAL_V6.fits'
)
assert (
test_dict['ccube'] ==
'counts_cubes/ccube_P8_P302_8years_source_none_zmax105_E3_PSF3_GAL.fits'
)
assert (
test_dict['mcube'] ==
'model_cubes/mcube_ptsrc_P8_P302_8years_source_none_zmax105_E3_PSF3_GAL_V6.fits'
)
assert (test_irfs == 'P8R2_SOURCE_V6')
import math
from fermipy.jobs.file_archive import FileFlags
from fermipy.jobs.chain import Link
from fermipy.jobs.gtlink import Gtlink
from fermipy.jobs.scatter_gather import ConfigMaker, build_sg_from_link
from fermipy.jobs.lsf_impl import make_nfs_path, get_lsf_default_args, LSF_Interface
from fermipy.diffuse.utils import create_inputlist
from fermipy.diffuse.name_policy import NameFactory
from fermipy.diffuse.binning import Component
from fermipy.diffuse.diffuse_src_manager import make_ring_dicts,\
make_diffuse_comp_info_dict
from fermipy.diffuse.catalog_src_manager import make_catalog_comp_dict
from fermipy.diffuse import defaults as diffuse_defaults
NAME_FACTORY = NameFactory()
def make_input_file_list(binnedfile, num_files):
"""Make the list of input files for a particular energy bin X psf type """
outdir_base = os.path.dirname(binnedfile)
outbasename = os.path.basename(binnedfile)
filelist = ""
for i in range(num_files):
split_key = "%06i" % i
output_dir = os.path.join(outdir_base, split_key)
filepath = os.path.join(
output_dir, outbasename.replace('.fits',
'_%s.fits.gz' % split_key))
filelist += ' %s' % filepath
return filelist
from fermipy import fits_utils
from fermipy.utils import load_yaml
from fermipy.jobs.file_archive import FileFlags
from fermipy.jobs.scatter_gather import ScatterGather
from fermipy.jobs.slac_impl import make_nfs_path
from fermipy.jobs.link import Link
from fermipy.jobs.chain import Chain
from fermipy.diffuse.binning import Component
from fermipy.diffuse.name_policy import NameFactory
from fermipy.diffuse import defaults as diffuse_defaults
from fermipy.diffuse.gt_split_and_mktime import SplitAndMktimeChain
NAME_FACTORY = NameFactory()
NAME_FACTORY_CLEAN = NameFactory()
NAME_FACTORY_DIRTY = NameFactory()
class ResidualCR(Link):
"""Small class to analyze the residual cosmic-ray contaimination.
"""
appname = 'fermipy-residual-cr'
linkname_default = 'residual-cr'
usage = '%s [options]' % (appname)
description = "Compute the residual cosmic-ray contamination"
default_options = dict(ccube_dirty=diffuse_defaults.residual_cr['ccube_dirty'],
ccube_clean=diffuse_defaults.residual_cr['ccube_clean'],
bexpcube_dirty=diffuse_defaults.residual_cr['bexpcube_dirty'],