names=('ObsId', 'RevolutionId', 'ExposureID', 'Start', 'End',
'Duration_Obs', 'EPIC_rate_soft', 'EPIC_erate_soft',
'EPIC_rate_hard', 'EPIC_erate_hard', 'hr', 'fvar_soft',
'efvar_soft', 'xs_soft', 'exs_soft', 'nxs_soft', 'enxs_soft',
'VA_soft', 'eVA_soft', 'fvar_hard', 'efvar_hard', 'xs_hard',
'exs_hard', 'nxs_hard', 'enxs_hard', 'VA_hard', 'eVA_hard'),
dtype=('i', 'i', 'U9', 'U30', 'U30', 'd', 'd', 'd', 'd', 'd', 'd', 'd',
'd', 'd', 'd', 'd', 'd', 'd', 'd', 'd', 'd', 'd', 'd', 'd', 'd',
'd', 'd'))
i = 0
for obsid in epic_observations:
obs = Observation(obsid=obsid, target_dir=target_dir)
print('-----------------------------------------------------------')
#Process each observation
obs.cifbuild()
obs.odfingest()
obs.epproc()
obs.filter_epic(pileup=True)
obs.epiclccorr(pileup=True) #always correct for pile-up
obs.epic_lightcurve(mjdref=mjdref)
obs.fracvartest(instrument='epic')
obs.pn_spectrum(pileup=True)
obs.pn_xspec(target_REDSHIFT)
#Save attributes of observation into the EPIC_table (to edit)
EPIC_obs_table.add_row(
(str(obs.obsid), str(obs.revolution), f"{obs.epic_expid}",
str(obs.starttime), str(obs.endtime), str(int(
obs.duration)), obs.epicrate[0], obs.epic_erate[0],
obs.epicrate[1], obs.epic_erate[1], obs.mean_hr,
"/Users/kym/PycharmProjects/astrostat/data/A1589/0149900301/ODF",
"ccf_path": "/Users/kym/Desktop/software/ccf/",
"mos1_prefix": "1S001",
"mos2_prefix": "2S002",
"pn_prefix": "S003",
"mos1_ccds": (1, 1, 1, 0, 1, 1, 1),
"mos2_ccds": (1, 1, 1, 1, 1, 1, 1),
"pn_ccds": (1, 1, 1, 1),
"ra": "190.32455",
"dec": "18.5725719"
}
A1589 = Observation(A1589_info_dict)
A1589.set_folders()
A1589.info()
A1589.cifbuild()
A1589.odfingest()
A1589.chain()
#A1589.copy_events()
#A1589.epic_filter()
A1589.convert_coords()
#A1589.cheese()
A1589.create_region(mos1=True,
mos2=True,
pn=True,
ra=None,
dec=None,
inner_radius=0,
outer_radius=300)
A1589.list_regions()