Last modified: December 2013

URL: http://cxc.harvard.edu/sherpa/ahelp/load_syserror.html
Jump to: Description · Examples · See Also


AHELP for CIAO 4.11 Sherpa v1

load_syserror

Context: data

Synopsis

Load systematic errors

Syntax

load_syserror( [id=1], filename, [bkg_id], [options] )

Description

The load_syserror command loads a user-input systematic error array from a FITS table or column-based text file for 1D data, or an error map from a FITS image for 2D data. The source dataset must be loaded before the systematic errors can be associated with it. The get_syserror command is available for accessing the array(s) of systematic error values.

Systematic errors are added in quadrature with statistical errors (if specified) to produce an overall error estimate in each data bin. If the systematic errors are not loaded, it is assumed that they are zero.

Options for a FITS file: [ncols=2,colkeys=None]

Options for an ASCII file: [ncols=2, colkeys=None, sep=" ", comment='#']


Examples

Example 1

sherpa> load_syserror("src", "syserror.fits[cols channel,syserr]")

The systematic errors for the dataset with id "src" are loaded from the columns "channel" and "syserrs" of the file syserror.fits.

Example 2

sherpa> load_syserror("syserror.dat", bkg_id=1)

Systematic errors for the background with bkg_id=1 of the default dataset are read from the first two columns (the default) of the ASCII file syserror.dat.

Example 3

sherpa> load_image(3, "src_image.fits")
sherpa> load_syserror(3, "syserror_img.fits")

Source image data is loaded into the Sherpa session and assigned to data set ID 3 with the load_image command. The systematic errors map associated with the source image is loaded for data set 3 with the load_syserror command.

sherpa> load_image(3, "src_image.fits")
sherpa> load_syserror(3, "syserror_img.fits")

sherpa> show_data(3)
Data Set: 1
Filter:
name = src_image.fits
x0 = Float64[56376]
x1 = Float64[56376]
y = Float64[56376]
shape = (216, 261)
staterror = None
syserror = Float64[56376]
sky = physical
crval = [ 3798.5 4019.5]
crpix = [ 0.5 0.5]
cdelt = [ 2. 2.]
eqpos = world
crval = [ 278.386 -10.5899]
crpix = [ 4096.5 4096.5]
cdelt = [-0.0001 0.0001]
crota = 0
epoch = 2000
equinox = 2000
coord = logical

See Also

data
dataspace1d, dataspace2d, datastack, fake, load_arf, load_arrays, load_ascii, load_bkg, load_bkg_arf, load_bkg_rmf, load_data, load_grouping, load_image, load_multi_arfs, load_multi_rmfs, load_pha, load_quality, load_rmf, load_staterror, load_table, pack_image, pack_pha, pack_table, unpack_arf, unpack_arrays, unpack_ascii, unpack_bkg, unpack_data, unpack_image, unpack_pha, unpack_rmf, unpack_table
filtering
load_filter
info
get_default_id, list_bkg_ids, list_data_ids
modeling
add_model, add_user_pars, load_table_model, load_template_model, load_user_model, save_model, save_source
saving
save_arrays, save_data, save_delchi, save_error, save_filter, save_grouping, save_image, save_pha, save_quality, save_resid, save_staterror, save_syserror, save_table
statistics
load_user_stat