SPIRE instrument and calibration web pages
Back to main Herschel Public TWiki page
Introduction
This page provides up-to-date information about using the SPIRE instrument: from preparing observations to reducing your data. This page also provides you with the latest calibration accuracies and known SPIRE calibration issues.
Observing with SPIRE
The most up to date information on instrument calibration and performance is given in the SPIRE Observers' Manual. This is the reference document used by all the rest of the SPIRE user guides (eg data reduction guide, cookbooks etc). Sometimes it may happen that outdated values are quoted in some of the documents. In such a case use the values given in the SPIRE Observers' Manual.
- The SPIRE Observer's Manual (HTML) and PDF (15.5 MB).
- A quick guide to the SPIRE instrument: SPIRE fact sheet
- Summary papers from the A&A Special Issue:
-
Some values on the performance are now outdated. Please consult the SPIRE Observers' Manual for most up to date information.
- The Herschel-SPIRE instrument and its in-flight performance, Griffin et al., 2010, A&A, 518, L3, in-flight scientific capabilities of the SPIRE instrument.
- In-flight calibration of the Herschel-SPIRE instrument: Swinyard et al., 2010, A&A, 518, L4
- Detailed documents describing the pipeline algorithms (to be updated soon):
AOT release notes
Reducing SPIRE data
Software and documentation
- HIPE (Herschel Interactive Processing Environment): The latest User Release HCSS version that you should use for reducing SPIRE data is HIPE v10.1. It can be downloaded from: http://herschel.esac.esa.int/HIPE_download.shtml.
-
Please note that there was a bug in the destriper task included in HIPE 9.0 that may affect your final map, especially if there are bright objects in the observed field. This has been corrected since HIPE 9.1. If your observation falls in the mentioned category, you are strongly advised to update your HIPE installation.
- We also provide access to the latest stable developer build (latest stable CIB).
- Beware These developer builds do not undergo the same in-depth testing as the user releases do. The latest developer build can be found here.
- Within HIPE you can access all the SPIRE data reduction and HIPE-user documentation. The SPIRE Data Reduction Guide (SDRG) follows the user pipeline scripts and also explains the details of pipeline processing and data analysis. It is also available online here:
- SPIA: The SPIRE Photometer Interactive Analysis (SPIA) package is available as a plug-in for HIPE. SPIA provides a structured GUI based access to the more intricate parts of the scan map photometer pipeline for SPIRE without the immediate need to resort to scripts. More information can be found in the SDRG or on the SPIA web page
The SPIRE Launch Pads
- The SPIRE Launch Pads are single sheet quick entries (like a cheat sheet) into SPIRE data reduction and providing quick references to the relevant sections in the SPIRE Data Reduction Guide. There are launch pads for Data Access, SPIRE Photometer and Spectrometer data reduction.
Photometer data reduction
Overview
The best source of information for reducing SPIRE Photometer data is the
SPIRE Data Reduction Guide available through the HIPE help. This runs through the User Pipeline scripts step by step, describes several other Useful Scripts, and offers advice for specific issues that might be encountered.
New definition of Leve2 products
- For versions of the HCSS prior to HIPE 10.0, a single point source calibrated (Jy/beam) map was provided in the Level 2 product for each of the PSW, PMW, PLW bands. However, for observations processed with HIPE 10.0 or later, more than one map calibration is made available within the Level 2 product. Maps are provided for the following scenarios for post HIPE 10.0 processing:
Description |
New Name (HIPE10+) |
Original (pre-HIPE10) Name |
Point Source (standard) Maps |
psrcPSW |
PSW |
Extended Emission Maps |
extdPSW |
- |
Solar System Object Maps |
ssoPSW |
- |
Point Source Destriper Diagnostic |
psrcPSWdiag |
pddPSW |
Extended Emission Destriper Diagnostic |
extdPSWdiag |
- |
|
|
|
- psrcPxW are the previous PxW maps, calibrated for point source and in units of Jy/beam. Note that to do aperture photometry on such maps you'll first need to convert them to surface brightness (Jy/pixel, MJy/sr, etc.), although it is suggested to directly use the already extended emission calibrated extdPxW maps. Finally, bear in mind that SPIRE itself cannot measure the absolute sky flux, hence psrcPxW maps have an arbitrary offset having zero median.
- ssoPxW maps are corrected for SSO proper motion: maps are in Jy/beam and they are subject to the same photometry rules of the psrcPxW maps.
- extdPxW maps are calibrated for extended emission and provided in units of MJy/sr. These maps are provided with an estimation of the absolute offset via cross-calibration with Planck data.
- In all cases, SPIRE data is calibrated in the assumption of source having a spectral index equal to -1, i.e. νSν = const. To calibrate your data for other cases or convert to e.g. Jy/sr, please refer to section 5.7 of the SPIRE Data Reduction Guide.
- The SPIRE Photometer filter transmission curves, also known as Relative Spectral Response Functions (RSRF) are available here
. For more details, please read the .readme file in this ftp folder.
Data Processing Issues
The main issues that you might find in your data are:
undetected glitches,
thermistor or detector jumps, and
bad baseline removal.
- Stripes in PSW, PMW and/or PLW (Level 2) maps
- All SPIRE photometry pipelines now use the destriper by default, which improves the issue of stripes in Level 2 maps. There should be noticeable improvements in that respect with HIPE version 9. The destriper documentation can be found on the NHSC website
- De-glitcher masks faint sources
- For data taken in Parallel Mode in particular (sampling at 10Hz, at high speed 60"/s), the de-glitcher may flag very faint sources as glitches when it is run with standard parameters. Faint sources may have a "delta function" shape due to the low sampling rate, which looks similar to a small glitch. Try modifying the "correlation parameter" to 0.95: this will decrease the number of detected glitches - it may be better to have a limited detection rate in first level deglitching and defer to Level 2 deglitching.
- Cooler temperature variations (Cooler Burps)
- The SPIRE cooler is recycled after 48 hours. Between 6 to 7h after the cooler recycle ends, its temperature rises steeply and reaches the stable plateau. Observations taken during such times may exhibit stripes in the final maps. An option to correct for this effect is now available in the User Pipelines (See the SPIRE DRG for details).
- NaN pixels present in the PSW, PMW and/or PLW (Level 2) maps
- This effect, related to data masking or poor coverage, is more evident in single fast-scan Parallel Mode maps. To avoid NaNs, increase the pixel size (i.e., decrease the map's resolution).
- This effect can also occur with destriped maps. In this case check if increasing the sigma parameter or switching off the Level 2 deglitcher helps.
- WCS in 3-colour images
- Problems with the wrong WCS in the output RGB images in all observations reduced with HIPE 8 have been fixed in HIPE 9
- Quality flags in the quality context
- Currently, the quality flags at the quality context inside the observation context are just meant for HSC/ICC internal evaluation of the quality of the products and not for the users. In case the data had some serious quality problem, the PI of the program has been contacted about it. Otherwise, only information in the quality summary, when available, should concern the observers.
Known Issues in ODs 1304 & 1305
For (yet) unknown reasons, the three detectors PSW-B5, PSW-E9 and PSW-F8 - that use to behave well during the entire mission - were noisy during the two operational days 1304 and 1305. The result are stripes visible in the final PSW map which the current (HIPE 10) pipeline is not able to correct. The solution is to mask and exclude these detectors from the analysis. This could be done in 2 ways:
- You can use the
SpireMaskEditor
GUI as described in sec. 7.4.2 of the SPIRE Data Reduction Guide: write-click on your observation context variable and then select Level1_SpireMaskEditor
and set to Master
all samples in all scans (listed as BBID) for the detectors mentioned above.
- Alternatively, you can use the following line of codes:
##### EDIT MASK SCRIPT BEGINS #####
# List of detectors to be masked
bolos = ['PSWB5', 'PSWE9', 'PSWF8']
# Level1 of your observation, assuming the observation context varible is named 'obs'
level1 = obs.level1
# Create new level 1
new_l1 = Level1Context()
for scan in range(0, level1.getCount()):
#
# Load level 1 product, scan by scan
data = level1.refs[scan].product
#
# Change mask for selected detectors in all scans, setting it to MASTER
for bolo in bolos:
data['mask'][bolo].data[:] = 1
#
new_l1.addProduct(data)
##### EDIT MASK SCRIPT ENDS #####
- After either of those cases, you must then re-run level 1 to 2 steps on the newly modified level1 product. If your observation has been already re-reduced with HIPE 10, original and new level1s are already destriped, so you can directly run the naive map-maker on the new level1. Otherwise, you must run the destriper step: check the pipeline script for details.
Source Extraction and Photometry
- The current recommended method for photometry
sourceExtractorTimeline
task (formerly known as the Timeline Fitter) which works on the detector timelines. The Map based algorithm sourceExtractorSussex
(SUSSEXtractor) providers good results and is useful on larger maps where the sourceExtractorTimeline
will be significantly slower. sourceExtractorDaophot
(DAOphot) also provides a reasonable estimate of the source flux but may require an aperture correction.
Spectrometer data reduction
The best source of information for reducing SPIRE Spectrometer data is the
SPIRE Data Reduction Guide available through the HIPE help. This runs through the User Pipeline scripts step by step, describes several other Useful Scripts, and offers advice for specific types of sources:
- Faint (<10 Jy) and medium (<100 Jy) strength sources
- Bright sources (>500 Jy)
- Extended sources
- Observations with few repetitions
- H+L observations
For faint sources, the subtraction of instrument, telescope and background emission is particularly important. Optimum subtraction can be performed in several ways (read the
SPIRE Data Reduction Guide for details):
- Subtract the Dark Sky spectrum closest to your observation (use the "Background Subtraction" script in HIPE)
- Subtract the spectrum of surrounding detectors (use the "Background Subtraction" script in HIPE)
Dark Sky observations are observed on every SPIRE Spectrometer OD, and are all public in the Archive.
A listing of the available Dark Sky observations can be found here.
Spectral Cube Analysis
Spectral cubes are produced by the Spectrometer pipeline for mapping observations. Some tips, suggestions and examples of spectral cube analysis for SPIRE data are provided
here.
Cookbooks
Cookbooks are provided inside the SPIRE Data Reduction Guide (see above).
The standalone "Photometry Cookbook", is no longer maintained - it is being incorporated into the SPIRE DRG - please see the SDRG for photometry cookbook information, and raise a
Helpdesk ticket if you find something missing.
SPIRE calibration file versions
The available calibration trees for SPIRE are listed below (with the current operational version at the top).
SPIRE Calibration Tree |
Applicable HIPE Version |
Comment |
SPIRE_CAL_9_1 |
HIPE v9 |
Final v9 cal tree |
SPIRE_CAL_8_1 |
HIPE v8 |
Final v8 cal tree |
SPIRE_CAL_7_0 |
HIPE v7 |
Final v7 cal tree. |
SPIRE_CAL_6_1 |
HIPE v6 |
Final v6 cal tree |
(SPIRE_CAL_6_0) |
HIPE v6 |
Spec major update |
SPIRE_CAL_5_2 |
HIPE v5 |
Final v5 cal tree |
(SPIRE_CAL_5_1) |
HIPE v5 |
|
(SPIRE_CAL_5_0) |
HIPE v5 |
Phot flux conv. based on Neptune. Spec major update |
SPIRE_CAL_4_0 |
HIPE v4 |
Spec point source flux conv based on Uranus |
SPIRE_CAL_3_2 |
HIPE v3 |
|
(SPIRE_CAL_3_1) |
HIPE v3 |
|
(SPIRE_CAL_3_0) |
HIPE v3 |
|
SPIRE_CAL_2_1 |
HIPE v2 |
Spec point source flux conv based on Vesta |
(SPIRE_CAL_2_0) |
HIPE v2 |
|
SPIRE_CAL_10_1 |
HIPE v10 |
Calibration tree currently used in operations |
SPIRE_CAL_1_2 |
HIPE v1 |
Phot flux conv based on Ceres |
(SPIRE_CAL_1_1) |
HIPE v1 |
Pre-launch dummy values |
More details of the changes in each version are given
here
.
- Any of the calibration trees can be retrieved in HIPE from the HSA using (e.g.)
cal = spireCal(calTree="spire_cal_10_1")
etc. The default (applicable to the HIPE version in use) can be obtained with cal = spireCal(calTree="spire_cal")
. It can then be saved to a local pool right-clicking on the cal
variable and then selecting from the context menu Send To -> Local Pool.
- Alternatively, the latest calibration tree for SPIRE can be obtained as a jar file from Latest calibration trees
. Then, you have to possibilities to read and save:
- The jar file can be load directly into HIPE with the command:
cal = spireCal(jarFile="PATH_TO_FILE/spire_cal_10_1.jar")
. To save it to a local pool, proceed as described above, right-clicking on the cal
variable and then selecting from the context menu Send To -> Local Pool.
- The jar file can also be saved directly to a local pool without opening HIPE, running the following command in the terminal command line:
cal_import PATH_TO_FILE/spire_cal_10_1.jar
. Then, to load the calibration tree in HIPE, simply type: cal = spireCal(pool="spire_cal_10_1")
See
the SPIRE Data Reduction Guide for more details.
SPIRE calibration and performance
Photometer calibration
- SPIRE Photometer Calibration:
- Full details of the SPIRE calibration can be found in the SPIRE Observers Manual and in dedicated journal papers (Bendo et al. 2013, North et al. 2013)
- Point Source Flux Calibration ~ ± 7%
- ~ 5% absolute from Neptune model
- ~2% relative from Neptune photometry
- Extended emission Calibration ~ ± 11%
- additional ~4% due to uncertainty in the measured beam area
- Additionally repeatability of measured photometry to flux densities >100mJy ~ 2%
- SPIRE Photometer Beams:
- These are available in the SPIRE calibration context, at the standard map pixel size of (6,10,14) arcsec/pixel for (250,350,500) µm bands, and can be accessed in HIPE after a calibration context has been loaded (see above).
-
The observed beams at much finer scale of 1 arcsec/pixel, as well as the theoretical ones, are available from here
. Please read the release note for more details.
-
A new more detailed analysis of the SPIRE beam profile data was undertaken in 2012, leading to revised values for beam profile solid angles and derivation of a semi empirical wavelength dependent beam profile model. The results at a scale of 1 arcsec/pixel as well as the data needed for the model are available for download. A detailed description of the analysis is given as well.
- SPIRE Photometer filter transmission curves:
- These are also available in the SPIRE calibration context (photRsrf) and can be accessed in HIPE after a calibration context has been loaded (See above).
- Neptune and Uranus models used for the SPIRE photometer flux calibration:
- The ESA2 models currently used in the SPIRE calibration are available here
.
Spectrometer calibration
Important FTS information, including details of the calibration, point source and extended source calibration etc, is available in the SPIRE Observers' Manual, Sections 4.2 and 5.3. These two sections are a must-read for anybody processing SPIRE FTS data.
Interest groups and scripts
- The following interest groups relate to processing of observations taken with SPIRE. The links provided allow subscription to these interest groups.
- subscribe
to the SPIRE Photometer interest group
- subscribe
to the SPIRE and PACS parallel mode large map and point source extraction interest group
- subscribe
to the SPIRE Spectrometer interest group
- subscribe
to the PACS, SPIRE and HIFI spectral maps interest group
- User contributed scripts: Users are welcome to submit scripts and software that they believe could be of general interest to the community to the Herschel http://herschel.esac.esa.int/esupport/Helpdesk.
Further information
Contact the
Helpdesk