'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)