def test_scwid_list():
import ddosa
import dataanalysis.core as da
da.debug_output()
fa = ddosa.IDScWList(use_scwid_list=["066500230010.001"])
print(fa.get())
print(fa.get_version())
assert len(fa.scwlistdata) == 1
fa = ddosa.IDScWList(
use_scwid_list=["066500230010.001", "066500250010.001"])
print(fa.get())
print(fa.get_version())
assert len(fa.scwlistdata) == 2
fa = ddosa.IDScWList(
use_scwid_list=["066500230010.001", "066600250010.001"])
print(fa.get())
print(fa.get_version())
assert len(fa.scwlistdata) == 2
def test_mosaic_ii_skyimage():
import ddosa
import gencat
cat = gencat.SourceCatalog()
cat.catalog = [dict(
NAME="TEST_SOURCE",
RA=83,
DEC=22,
)]
cat.promote()
assert cat.virtual
assert hasattr(cat, 'catalog')
grcat = gencat.GRcat()
grcat.serialize()
sc = ddosa.da.AnalysisFactory.byname("SourceCatalog")
assert cat == sc
print sc.get_version()
#ips=ddosa.ImageProcessingSummary().get()
#print "HASH:",ips.hash
##
#da.debug_output()
mf = gencat.GRcat()
c = mf.process(output_required=False, run_if_haveto=False)[0]
print("mf:", c)
mf = gencat.GRcat(assume=ddosa.ScWData(input_scwid="any"))
c = mf.process(output_required=False, run_if_haveto=False)[0]
print("mf:", c)
assert c[1][2].startswith("SourceCatalog.v0.1.")
##
mosaic = ddosa.mosaic_ii_skyimage(assume=[
ddosa.ii_skyimage(input_cat=gencat.CatForImage),
ddosa.ImageGroups(input_scwlist=ddosa.IDScWList(
use_scwid_list=["066500330010.001", "066500340010.001"])),
ddosa.ImageBins(use_ebins=[(25, 40)], use_version="onebin_25_40"),
ddosa.ImagingConfig(use_SouFit=0, use_version="soufit0")
])
mosaic.read_caches = []
mosaic.ii_NegModels = 1
mosaic.get()
d = pyfits.open(mosaic.srclres.get_path())[1].data
print d['NAME']
assert len(d) == 2
assert d[0]['NAME'] == cat.catalog[0]['NAME']
assert d[0]['RA_FIN'] == cat.catalog[0]['RA']
assert d[0]['DEC_FIN'] == cat.catalog[0]['DEC']
def test_mosaic_ii_skyimage():
#da.debug_output()
ig = ddosa.ImageGroups(input_scwlist=ddosa.IDScWList(
use_scwid_list=["066500230010.001", "066500240010.001"]))
fa = ddosa.mosaic_ii_skyimage(input_imagegroups=ig)
fa.read_caches = []
fa.get()
def test_mosaic_ii_skyimage():
import ddosa
mosaic=ddosa.mosaic_ii_skyimage(
assume=[
ddosa.ImageGroups(input_scwlist=ddosa.IDScWList(use_scwid_list=["066500330010.001","066500340010.001"])),
ddosa.ImageBins(use_ebins=[(25, 40)], use_version="onebin_25_40"),
ddosa.ImagingConfig(use_SouFit=0, use_version="soufit0")
]
)
mosaic.get()
def test_construct_many_images():
ig = ddosa.ImageGroups(input_scwlist=ddosa.IDScWList(
use_scwid_list=["066500230010.001", "066500240010.001"]))
ig = ig.get()
print ig
assert len(ig.members) == 2
assert len(ig.members[0]) == 5
assert hasattr(ig.members[0][1], 'skyima')
assert os.path.exists(ig.members[0][1].skyima.get_path())
ig.construct_og("og.fits")
def test_lc_pick():
#da.debug_output()
from dataanalysis import analysisfactory
analysisfactory.AnalysisFactory.WhatIfCopy(
"timebin", ddosa.LCTimeBin(use_time_bin_seconds=200))
ig = ddosa.LCGroups(input_scwlist=ddosa.IDScWList(
use_scwid_list=["066500230010.001", "066500240010.001"]))
fa = ddosa.lc_pick(input_lcgroups=ig, use_source_names=["Crab"])
fa.read_caches = []
fa.get()
# assert os.path.exi
def test_spectra_sum():
da.reset()
import process_isgri_spectra
import ddosa
import gencat
cat = gencat.SourceCatalog()
cat.catalog = [dict(
NAME="TEST_SOURCE",
RA=83,
DEC=22,
)]
cat.promote()
process_isgri_spectra.ScWSpectraList(input_scwlist=ddosa.IDScWList(
use_scwid_list=["066500330010.001", "066500340010.001"])).promote()
iss = process_isgri_spectra.ISGRISpectraSum(use_extract_all=True)
iss.get()