def test_create_roi_from_source():
rm = ROIModel.create_from_source('3FGL J2021.0+4031e',
{'catalogs': ['3FGL'], 'src_radius': 2.0})
assert len(rm.sources) == 9
src = rm.sources[0]
assert src.name == '3FGL J2021.0+4031e'
assert src['SpatialType'] == 'SpatialMap'
assert src['SourceType'] == 'DiffuseSource'
assert src['SpectrumType'] == 'PowerLaw'
assert_allclose(src['ra'], 305.26999, rtol=1E-5)
assert_allclose(src['dec'], 40.52, rtol=1E-5)
assert_allclose(src.spectral_pars['Index']['value'], 1.53, rtol=1E-4)
assert_allclose(src.spectral_pars['Prefactor']['value'] *
src.spectral_pars['Prefactor']['scale'], 0.4003659112E-12, rtol=1E-4)
assert_allclose(src.spatial_pars['Prefactor']['value'], 1.0, rtol=1E-4)
assert_allclose(src.spatial_pars['Prefactor']['scale'], 1.0, rtol=1E-4)
def test_create_roi_from_source():
rm = ROIModel.create_from_source('3FGL J2021.0+4031e',
{'catalogs': ['3FGL'], 'src_radius': 2.0})
assert len(rm.sources) == 9
src = rm.sources[0]
assert src.name == '3FGL J2021.0+4031e'
assert src['SpatialType'] == 'SpatialMap'
assert src['SourceType'] == 'DiffuseSource'
assert src['SpectrumType'] == 'PowerLaw'
assert_allclose(src['ra'], 305.26999, rtol=1E-5)
assert_allclose(src['dec'], 40.52, rtol=1E-5)
assert_allclose(src.spectral_pars['Index']['value'], 1.53, rtol=1E-4)
assert_allclose(src.spectral_pars['Prefactor']['value'] *
src.spectral_pars['Prefactor']['scale'], 0.4003659112E-12, rtol=1E-4)
assert_allclose(src.spatial_pars['Prefactor']['value'], 1.0, rtol=1E-4)
assert_allclose(src.spatial_pars['Prefactor']['scale'], 1.0, rtol=1E-4)
def test_create_roi_from_source_4fgldr2():
rm = ROIModel.create_from_source('4FGL J2021.0+4031e',
{'catalogs': ['4FGL-DR2'], 'src_radius': 2.0})
assert len(rm.sources) == 6
src = rm.sources[0]
assert src.name == '4FGL J2021.0+4031e'
assert src['SpatialType'] == 'RadialDisk'
assert src['SourceType'] == 'DiffuseSource'
assert src['SpectrumType'] == 'LogParabola'
assert_allclose(src['ra'], 305.27, rtol=1E-5)
assert_allclose(src['dec'], 40.52, rtol=1E-5)
assert_allclose(src.spectral_pars['alpha']['value'], 1.8747, rtol=1E-4)
assert_allclose(src.spectral_pars['beta']['value'] *
src.spectral_pars['beta']['scale'], 0.0606, rtol=1E-4)
assert_allclose(src.spectral_pars['norm']['value'] *
src.spectral_pars['norm']['scale'], 2.0858e-13, rtol=1E-4)
assert_allclose(src.spatial_pars['Radius']['value'], 0.63, rtol=1E-4)
assert_allclose(src.spatial_pars['DEC']['scale'], 1.0, rtol=1E-4)
