The SPIRE FTS part

The script prepSpireSpec.py provides the interface to combine the two SPIRE FTS spectra from the central detectors SSWD4 and SLWC3. The output of this procedure is one Spectrum1d with subsequent segment numbers for SSW (194-313 μm) and SLW (303-671 μm) arrays.

Here is the full description of the script:

    """
    NAME: 
        prepSpireSpec
    PURPOSE: 
        combine the two FTS arrays into one Spectrum1d object, 
            each FTS array will be with a different segment number
    USAGE: 
        result = prepSpireSpec(obsContext,startSegment=0,apod="unapodized",waveUnit=None,shiftS=0.0,shiftL=0.0)
    INPUTS:
        obsContext - a SPIRE FTS observation context, sparse spatial sampling only
        startSegment - the starting index of the segment, default is 0 => SSW will be segment 0, SLW will be 1.
        apod - one of "unapodized" (default) or "apodized", which variant of the spectra to use
        waveUnit - if None then use the default unit (wavenumbers in cm-1)
                 can be "micrometer" or "GHz" for the unit of the output spectrum.
        shiftS - a flux shift to be aplied to the SSW spectrum, i.e. to match with the photometry, in the same units as the flux
        shiftL - a flux shift to be applied to the SLW spectrum, i.e. to match with the photometry
    OUTPUT:
        object of class Spectrum1d with the combined SSW and SLW spectra, each with a different segment number
    HISTORY:
        Created: Apr 2011, Ivan Valtchanov, HSC, ESAC, ESA, following interactions wth Katrina Exter
    MODIFICATION HISTORY:
        v1.0 - Apr 2011

Here is the simplest call of the script

obs = getObservation(1342189124,useHsa=True) # public NGC7027 data in the HSA
allSpireSpectra = prepSpireSpec(obsContext,waveUnit="micrometer")