Science Threads for CIAO 4.2
Introduction
Beginners should start here. The Introductory threads provide an overview of the main components (GUI applications, parameter files) and concepts (the Data Model, filtering) in the CIAO data analysis software.
Getting Started:
-
Installing CIAO 4.2 using the ciao-install script
(29 Dec 2009)
-
Installing CIAO 4.2
(30 Dec 2009)
- Introduction to CIAO
-
Starting CIAO
(16 Dec 2009)
-
Installing CIAO 4.2 using the ciao-install script
-
A Complete List of CIAO tools can be
found at the Introduction to Tools
page.
If you are looking for a general topic, try the CXC Search at the top of the page.
Data Preparation
When Chandra data goes through Standard Data Processing (SDP), the most recently available calibration is applied to it. Since this calibration is continuously being improved, one should check whether there are newer files available. Similarly, some science decisions are made during SDP; every user has the option to reprocess the data with different parameters.
-
See the current
Calibration Status Report
Analysis Guides:
General:
-
Setting the Observation-specific Bad Pixel Files
Uses: the acis_set_ardlib script - Filtering Data (S-Lang or Python)
- Filtering Lightcurves (S-Lang or Python)
Uses: the lightcurves.sl S-Lang script; the lightcurves.py script
-
Setting the Observation-specific Bad Pixel Files
Coordinates, Astrometry, & Spatial Filtering:
ACIS:
- Check the ACIS Data Caveats
- Why topic: -110 C Observations
- Why topic: Continuous Clocking Mode Data
Removing Cosmic-Ray Afterglows:
acis_process_events:
- Reprocessing Data to Create a New Level=2 Event File thread: apply the newest calibration (CTI, TGAIN); includes CC mode times of arrival calculation.
- Clean ACIS Background in VFAINT Mode
- The ACIS "Blank-Sky" Background Files (S-Lang or Python)
Uses: the acis_bkgrnd_lookup script; the lightcurves.sl S-Lang script; the lightcurves.py script - Remove the ACIS Readout Streak
HRC:
- New Observation-Specific HRC Bad Pixel File
- Reprocessing Data to Create a New Level=2 Event File thread: applies the newest calibration (degap, AMP_SF/tapringing).
- Computing Average HRC Dead Time Corrections (S-Lang or Python)
Imaging
The Imaging threads cover a wide range of topics that include source detection, creating exposure maps and normalized images, and calculating image statistics. How to create color images for publication is addressed, as well as merging data from multiple observations.
Analysis Guides:
General:
- The ACIS "Blank-Sky" Background Files (S-Lang or Python)
Uses: the acis_bkgrnd_lookup script; the lightcurves.sl S-Lang script; the lightcurves.py script -
Match the Binning of an Image
Uses: the get_sky_limits script - An Estimated Background Image
- True Color Images
- True Color Images in ds9
-
An Image of Diffuse Emission
Uses: the mkBgReg.pl script; the mkSubBgReg.pl script - Using Data Cubes
- The ACIS "Blank-Sky" Background Files (S-Lang or Python)
Reprojecting Data:
-
Merging Data from Multiple Imaging Observations
Uses: the merge_all script - Overview: Reprojecting Files
- Correcting Absolute Astrometry with reproject_aspect
- Reprojecting Images: Making an Exposure-corrected Mosaic
- Reprojecting Coordinates of a Solar System Object
-
Object-centered Aspect Solution and Exposure Map
(12 Jan 2009)
-
Merging Data from Multiple Imaging Observations
Coordinates, Astrometry, & Spatial Filtering:
Source Characteristics:
- Estimate Source Counts in an Image
- Obtain and Fit a Radial Profile (S-Lang or Python)
- Compute Net Counts, Rate, or Flux for Point Sources
- Computing the Intensity Upper Limit for an Unresolved Source
- Calculate the Flux for a Position
- Measuring Source Extent
- Search for Variability in a Source (S-Lang or Python)
PSFs:
- PSFs with ChaRT: the Chandra Ray Tracer
- Create a PSF (S-Lang or Python)
Detect:
- Overview: Detecting Sources in Imaging Observations
- An Estimated Background Image
- Running celldetect
- Running vtpdetect
- Running wavdetect
- Using the Output of Detect Tools (S-Lang or Python)
Exposure Maps:
-
Using merge_all to Compute ACIS Exposure Maps and Fluxed Images
Uses: the merge_all script
(12 Jan 2009)
- Calculating Spectral Weights for mkinstmap (S-Lang or Python)
-
Single Chip ACIS Exposure Map and Fluxed Image Step-by-Step
Uses: the get_sky_limits script
(12 Jan 2009)
-
Multiple Chip ACIS Exposure Map and Fluxed Image Step-by-Step
Uses: the get_sky_limits script
(12 Jan 2009)
-
HRC-I Exposure Map and Fluxed Image
Uses: the get_sky_limits script -
HRC-S Exposure Map and Fluxed Image
Uses: the get_sky_limits script -
Match the Binning of an Image
Uses: the get_sky_limits script
-
Using merge_all to Compute ACIS Exposure Maps and Fluxed Images
Modeling & Fitting Imaging Data with Sherpa
(from the Sherpa analysis threads):
Imaging Spectroscopy
After extracting source and background PI or PHA spectra from an imaging observation, the appropriate response files (ARF, RMF) are created so that the data may be modeled and fit. In the case of multiple or extended sources, a weighted ARF and RMF are built for the spectral analysis.
Extracting Spectra & Creating Response Files:
These threads use CIAO contributed scripts to automate the spectral extraction and the creation of ARF and RMF response files.
specextract should be used for the analysis of extended sources only. Users working with point sources should instead use the psextract script. The mkwarf tool, which is used for ARF generation by specextract, may not produce accurate results for point sources. The mkarf tool, used by psextract, is preferred for point source extraction.
In certain cases, the mkacisrmf tool will need to be run separately to create the best possible RMFs for ACIS observations; see the thread for details.
- The ACIS "Blank-Sky" Background Files (S-Lang or Python)
Uses: the acis_bkgrnd_lookup script; the lightcurves.sl S-Lang script; the lightcurves.py script -
Using specextract to Extract ACIS Spectra and Response Files
for Extended Sources
(12 Jan 2009)
-
Using psextract to Extract ACIS Spectra and Response
Files for Pointlike Sources
Uses: the psextract script; the acis_fef_lookup script
(12 Jan 2009)
-
Coadding Spectra and Weighted Responses
Uses: the acisspec script
(12 Jan 2009)
- Creating ACIS RMFs with mkacisrmf
- The ACIS "Blank-Sky" Background Files (S-Lang or Python)
Step-by-step Analyses:
These threads give step-by-step instructions for creating spectra and response files.
- The ACIS "Blank-Sky" Background Files (S-Lang or Python)
Uses: the acis_bkgrnd_lookup script; the lightcurves.sl S-Lang script; the lightcurves.py script -
Weighting ARFs and RMFs: multiple sources
(12 Jan 2009)
-
Step-by-Step Guide to Creating ACIS Spectra for Pointlike Sources
Uses: the acis_fef_lookup script
(12 Jan 2009)
-
Extract a Spectrum from the ACIS Readout Streak
Uses: the acis_fef_lookup script -
Extracting a Spectrum of a Solar System Object
Uses: the acis_fef_lookup script - A Note on HRC Spectra
- The ACIS "Blank-Sky" Background Files (S-Lang or Python)
Modeling & Fitting Spectral Data with Sherpa
(from the Sherpa analysis threads):- Introduction to Fitting PHA Spectra (S-Lang or Python)
- Changing the grouping scheme of a data set within Sherpa (S-Lang or Python)
- Introduction to Fitting ASCII Data with Errors: Single-Component Source Models (S-Lang or Python)
- Simultaneously Fitting Two Data Sets (S-Lang or Python)
- Simulating 1-D Data: the Sherpa FAKE_PHA Command (S-Lang or Python)
- Simulating Chandra ACIS-S Spectra with Sherpa (S-Lang or Python)
- Fitting PHA Data with Multi-Component Source Models (S-Lang or Python)
- Independent Background Responses (S-Lang or Python)
- Using A Pileup Model (S-Lang or Python)
Grating Spectroscopy
If new calibration has been applied to the event file, the grating spectrum should be re-extracted as well. It is then possible to build grating response files (gARF, gRMF) in order to model and fit the data.
- Chandra Grating Analysis Page
- WebGUIDE: Interactive GUIDE for ATOMDB version 1.3
-
Chandra Grating
Data Archive and Catalog (TGCat): a browsable
interface to analysis-quality spectral products.
General:
- Examining Grating Spectra and Regions: PHA2 files (S-Lang or Python)
- Create a Color Spectrum
Problems with the Zero Order:
HETG/ACIS:
- HETG/ACIS-S Grating Spectra
- Grating Spectra for Multiple Sources - ACIS
- ACIS-S Grating RMFs
-
HETG/ACIS-S Grating ARFs
Uses: the fullgarf script
(12 Jan 2009)
- Grouping a Grating Spectrum
LETG/ACIS:
- LETG/ACIS Grating Spectra
- Grating Spectra for Multiple Sources - ACIS
- ACIS-S Grating RMFs
-
LETG/ACIS-S Grating ARFs
Uses: the fullgarf script
(12 Jan 2009)
- Grouping a Grating Spectrum
LETG/HRC-S:
-
LETG/HRC-S Grating Spectra
(28 Dec 2009)
- Applying Customized Background Regions to LETG/HRC-S Observations (S-Lang or Python)
- Grating Spectra for Multiple Sources - HRC
- Higher-order Responses for HRC-S/LETG Spectra
- HRC Grating RMFs
- LETG/HRC-S Grating ARFs
- Grouping a Grating Spectrum
-
LETG/HRC-S Grating Spectra
Combining Spectra & Preparation for Fitting:
-
Extract Coadded and Grouped Nth-Order Source &
Background Spectra and ARFs
Uses: the add_grating_orders script -
Add Grating Spectra and Average ARFs
Uses: the add_grating_spectra script -
PHA Background File for XSPEC
Uses: the tg_bkg script
-
Extract Coadded and Grouped Nth-Order Source &
Background Spectra and ARFs
Modeling & Fitting Grating Data with Sherpa
(from the Sherpa analysis threads):
Timing Analysis
In order to perform absolute timing analysis on a dataset, a barycenter correction must first be applied to the data. One may then create lightcurves and phase-binned spectra to look for variability in the source. These threads also provide information on working with data taken in the ACIS continuous clocking (CC) mode.
General:
- Why topic: Continuous Clocking Mode
- Apply Barycenter Correction
- Phase-binning a Spectrum (S-Lang or Python)
- See the S-lang/ISIS Timing Analysis Routines (SITAR) package which provides a set of functions and subroutines for timing analysis within ISIS.
Data Used in ThreadsHow to Download Chandra Data from the Archive | |||
| ObsID | Object | Instrument | Threads |
|---|---|---|---|
| 3 | Trapezium Cluster | HETG/ACIS-S | Grating Spectra for Multiple Sources - ACIS |
| 29 | Alpha Cen | LETG/HRC-S | Grating Spectra for Multiple Sources - HRC |
| 133 | PSR B0540-69 | ACIS-I | Apply Barycenter Correction Phase-binning a Spectrum |
| 144 | G21.5-0.9 | HRC-I | HRC-I Exposure Map and Fluxed Image Setting the Observation-specific Bad Pixel Files |
| 198 | Cas A | ACIS-S | True Color Images in ds9 |
| 313 | M31 | ACIS-S | Compute Net Counts, Rate, or Flux for Point Sources |
| 315 | NGC 4038/NGC 4039 | ACIS-S | An Image of Diffuse Emission Calculate the Flux for a Position |
| 441 | Chandra Deep Field South | ACIS-I | Correcting Absolute Astrometry with reproject_aspect |
| 459 | 3C 273 | HETG/ACIS-S | Remove the acis_detect_afterglow Correction New ACIS Bad Pixel File: Identify ACIS Hot Pixels and Cosmic Ray Afterglows Creating ACIS Spectra, ARFs, & RMFs for Pointlike Sources Examining Grating Spectra and Regions: PHA2 files HETG/ACIS-S Grating Spectra ACIS-S Grating RMFs HETG/ACIS-S Grating ARFs Grouping a Grating Spectrum PHA Background File for XSPEC Add Grating Spectra and Average ARFs Extract Coadded and Grouped Nth-Order Spectra and ARFs |
| 460 | 3C 273 | LETG/HRC-S | Examining Grating Spectra and Regions: PHA2 files LETG/HRC-S Grating Spectra HRC Grating RMFs LETG/HRC-S Grating ARFs Grouping a Grating Spectrum |
| 461 | 3C 273 | HRC-I | Basic Lightcurves Improving the Astrometry of your Data: Correct for a Known Processing Offset |
| 578 | 3C 295 | ACIS-S | Running wavdetect Running vtpdetect Running celldetect A Note on Processing Versions |
| 581 | Chandra Deep Field South | ACIS-I | Correcting Absolute Astrometry with reproject_aspect |
| 650 | GK Persei | ACIS-S | Creating Source and Background Files |
| 884 | 0235+164 | ACIS-S | Clean ACIS Background in VFAINT Mode |
| 953 | 3C 273 | HRC-I | Basic Lightcurves |
| 990 | Vega | HRC-S | HRC-S Exposure Map and Fluxed Image Calculate the Flux for a Position |
| 1198 | 3C 273 | LETG/ACIS-S | LETG/ACIS Grating Spectra LETG/ACIS-S Grating ARFs Examining Grating Spectra and Regions: PHA2 files |
| 1451 | II Peg | HETG/ACIS-S | Create a Color Spectrum |
| 1463 | Jupiter | ACIS-S | Reprojecting Coordinates of a Solar System Object |
| 1522 | Trapezium Cluster | ACIS-I | Running celldetect |
| 1557 | G21.5-0.9 | HRC-S | Setting the Observation-specific Bad Pixel Files |
| 1703 | PKS 2149-306 | LETG/ACIS-S | Correcting a Misplaced Zero-order Source Position |
| 1712 | 3C 273 | ACIS-S | Using SAOImage ds9 The ACIS "Blank-Sky" Background Files Using CIAO Region Files Filtering Lightcurves Remove the ACIS Readout Streak |
| 1800 | PKS 2155-304 | LETG/HRC-I | Examining Grating Spectra and Regions: PHA2 files |
| 1801 | PKS 2155-304 | LETG/HRC-I | LETG/HRC-I Grating Spectra LETG/HRC-I Grating ARFs |
| 1838 | G21.5-0.9 | ACIS-S | The ACIS "Blank-Sky" Background Files Estimate Source Counts in an Image True Color Images Match the Binning of an Image Multiple Chip ACIS Exposure Map and Fluxed Image Step-by-Step Single Chip ACIS Exposure Map and Fluxed Image Step-by-Step Obtain and Fit a Radial Profile Create a PSF |
| 1842 | G21.5-0.9 | ACIS-I | Using merge_all to Compute ACIS Exposure Maps and Fluxed Images Merging Data from Multiple Imaging Observations |
| 1843 | G21.5-0.9 | ACIS-I | Introductory Threads A Note on Processing Versions Setting the Observation-specific Bad Pixel Files Using merge_all to Compute ACIS Exposure Maps and Fluxed Images Merging Data from Multiple Imaging Observations |
| 2463 | 3C 273 | HETG/ACIS-S | Add Grating Spectra and Average ARFs |
| 3207 | 3C 294 | ACIS-S | Weighting ARFs and RMFs: multiple sources |
| 3726 | Jupiter | ACIS-S | Object-centered Aspect Solution and Exposure Map |
| 5969 | Vega | HRC-S | Computing Average HRC Dead Time Corrections |
| 6476 | HZ 43 | HRC-I | New Observation-Specific HRC Bad Pixel File |
| 9100 | SNR 0104-72.3 | ACIS-S | Coadding Spectra and Weighted Responses |
| 9810 | SNR 0104-72.3 | ACIS-S | Coadding Spectra and Weighted Responses |
![[CIAO Logo]](../imgs/ciao_logo_navbar.gif)