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,
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'),
#Process each observation
obs.cifbuild()
obs.odfingest()
obs.rgsproc()
obs.create_pairs_exposures()
if obs.skipobs == True:
print(
f'The observation {obsid} is tricky. Please analyse individually. Moving forward to next observation.'
)
print('----------------')
continue
obs.bkg_lightcurve()
obs.check_flaring_particle_bkgr()
obs.rgslccorr()
obs.lightcurve(mjdref=mjdref)
obs.fracvartest(screen=True, instrument='rgs')
obs.vaughan_panel(N=30,
M=20,
timescale=timescale_fvar,
timebinsize=25)
obs.divide_spectrum()
obs.xspec_divided_spectra_average(target_REDSHIFT)
obs.xspec_divided_spectra(target_REDSHIFT)
#Save attributes of observable into a table
if len(obs.rgsrate) == 0:
obs_table.add_row(
(str(obs.obsid), str(obs.revolution), None,
str(obs.starttime), str(obs.endtime),
str(int(obs.duration)), None, None, None, None, None,
None, None, None, None, None, None))
