sir - SIR spectrum manipulations

Single-dither and combined spectra from EUC_SIR_W-COMBSPEC*.fits.

Data repository:

  • repository: ousir:S1RD@ta@metis.lambrate.inaf.it
  • path: home/ousir/ftp/data/euclid/[prodname]/data/
  • FTP command: sftp://[repository]/[path]/EUC_SIR_W-COMBSPEC*.fits

Internal structure:

  • ext[0].header: TILE_ID and N_OBJ
  • for each object NNN from 0 to N_OBJ-1,
    • ext['N_META'].header: N_DITH, OBJECT_ID, RA_OBJ, DEC_OBJ
    • ext['N_COMBINED1D_SIGNAL'].data: columns WAVELENGTH, SIGNAL, QUALITY and VAR. QUALITY=0 for good data.
    • for each dither X from 1 to N_DITH:
      • ext['N_DITH1D_X_SIGNAL'], with X`=`DITH_ID (1-4)
      • data: same columns as combined extension
      • header: GWA_POS (str) and GWA_TILT (int)
Spectrum(w, f, df[, label]) A spectrum, made of (w, f ± df) masked-array.
Dithers() A dict of single-dither and combined spectra.
inspec.sir.warning_simple(message, category, filename, lineno, file=None, line=None)[source]
class inspec.sir.Spectrum(w, f, df, label='anonymous')[source]

A spectrum, made of (w, f ± df) masked-array.

binwidth = 13.4

Wavelength step [Å]

class inspec.sir.Dithers[source]

A dict of single-dither and combined spectra.

get_dithers(fdf=False)[source]

Return list of single-dither spectra, or stacked (f, df) arrays.

classmethod init_from_combspec(hdl, obj_idx)[source]

Initialize Dithers from object in HDU list.

combined_dither_masks()[source]

Combined mask from dithers.

combine(stat='wmean', clipping=None, **kwargs)[source]

Combine using stat and clipping, and store under key=’[clipping-]stat’.

stat: mean|wmean|median clipping: None|sigma|wsigma|robust|pull

plot(ax=None, discard=None)[source]

Plot combined and single-dither spectra.

get_pull(clip=5, side=0, maxiter=None, return_final=False)[source]

Get pull of individual dithers.

It should not be straight single-iteration pull values because of screening effects. Rather, it should iteratively correspond to pull values just before clipping.

inspec.sir.decipher_strategy(key)[source]

Decipher ‘[clipping-]stat’ strategy.

inspec.sir.plot_spectra(hdl, obj_ids, keys=['combined'], nperfig=50, **kwargs)[source]
inspec.sir.plot_pulls(hdl, obj_ids, clip=5, ax=None)[source]
inspec.sir.plot_stats(hdl, ax=None)[source]

Compute and plot stats about COMBSPEC HDU-list.