'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,
obs.fracvardict[0].get('Fractional Variability'),
obs.fracvardict[0].get('Fractional Variability Error'),
obs.fracvardict[0].get('Excess variance'),
obs.fracvardict[0].get('Excess variance error'),
obs.fracvardict[0].get('Normalized excess variance'),
obs.fracvardict[0].get('Normalized excess variance error'),
obs.fracvardict[0].get('Variability Amplitude'),
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()
mos1_spectrum_region = A1589.regions["mos1"][0]
mos2_spectrum_region = A1589.regions["mos2"][0]
pn_spectrum_region = A1589.regions["pn"][0]
A1589.mos1_spectrum(region_object=mos1_spectrum_region, mask=1)
A1589.mos2_spectrum(region_object=mos2_spectrum_region, mask=1)
A1589.pn_spectrum(region_object=pn_spectrum_region, mask=1)
A1589.rename_mos1_spectra_output(region_object=mos1_spectrum_region)
A1589.rename_mos2_spectra_output(region_object=mos2_spectrum_region)
A1589.rename_pn_spectra_output(region_object=pn_spectrum_region)
A1589.grppha_mos1(region_object=mos1_spectrum_region, bin_value=30)
A1589.grppha_mos2(region_object=mos2_spectrum_region, bin_value=30)
A1589.grppha_pn(region_object=pn_spectrum_region, bin_value=30)
A1589.proton_scale_mos1(region_object=mos1_spectrum_region)
A1589.proton_scale_mos2(region_object=mos2_spectrum_region)
A1589.proton_scale_pn(region_object=pn_spectrum_region)
