hxdpinxbpi @evts.list @pse.list @bkgs.list <outstem>
DESCRIPTION
hxdpinxbpi performs the following steps to produce output spectra:
- AND GTI from NXB file(s) with GTI in pin??_cl.evt file(s) (or spectrum
file) and any extra input GTI file(s) after ORing them together.
- Extract source spectrum using same GTI as step 1., if PIN cleaned event file is input.
- Extract NXB background spectrum using GTI from step 1.
- Calculate a simulated CXB spectrum, if requested
- Adjust the EXPOSURE in the NXB spectrum by 10.0 (background files are
scaled up by 10.0 to suppress poisson noise).
- Apply dead time correction to source spectrum using hxddtcorr.
- Add CXB spectrum to NXB spectrum, if CXB spectrum is created or given as input.
- Optionally group source spectrum.
OUTPUT FILES
The number of output files depends on the input parameters and types of input
files. If cleaned event files are used for the input_fname
parameter, cxb_fname != 'NONE'
and grouspec = "yes"
, then
the output files will be as follows:
<outstem>hxd_pin_sr.pi
- Source spectrum, corrected for dead time.
<outstem>hxd_pin_sr_grp.pi
- Grouped source spectrum
<outstem>hxd_pin_nxb.pi
- NXB background spectrum
<outstem>hxd_pin_cxb.pi
- CXB background spectrum
<outstem>hxd_pin_bg.pi
- Total background spectrum, NXB + CXB
The total background spectrum will only be produced if
cxb_fname='CALC'
or a CXB spectrum is supplied for the
cxb_fname
parameter.
PARAMETERS
- input_fname [string]
- Cleaned PIN event file(s) or PIN source spectrum. The GTI from this
file is ANDed with the GTI from the PIN NXB background file, and optionally
the GTI file(s) input in the gti_fname parameter (not supported for input
spectrum). The resulting GTI is used to extract a spectrum from the
cleaned PIN event file (if event files are input), the pseudo event file (see
pse_event_fname
parameter) and the PIN NXB file (see
bkg_event_fname
parameter).
A single spectrum for multiple event files can be extracted. Specify
@flist, where flist is a file containing cleaned PIN event files (one per
line) to use this feature. Extended filename syntax is supported for
individual files.
If a spectrum is entered, only a single file is allowed, and no extra GTI
filtering (see gti_fname param) is allowed. However, mutliple pseudo event
files and background files can be entered if necessary.
- pse_event_fname [string]
- Pseudo event file name or WEL unfiltered file name if pseudo event file is
unavailable. This file is used to correct the dead time of the source spectrum
extracted from the cl_event_fname event file(s).
If multiple event files are input for the cl_event_fname
parameter, then the files input for this parameter must correspond 1-1 to the
cl_event_fname
files. The same syntax is used as for the
cl_event_fname
parameter.
- bkg_event_fname [string]
- HXD PIN NXD file, either "tuned" or "quick". These can be downloaded from
either (tuned):
ftp://heasarc.gsfc.nasa.gov/suzaku/data/background/pinnxb_ver2.0_tuned/
or (quick):
ftp://heasarc.gsfc.nasa.gov/suzaku/data/background/pinnxb_ver2.0/
Alternatively, if the full ftp or http path to the file(s) is known, it can
be entered directly. There should be 1-1 correspondence between these
files and the cl_event_fname
files. The same syntax is used as for the
cl_event_fname
parameter. A warning will be issued if any of these files
are the quick background.
- outstem [string]
- Output filenames will start with this string.
- (gti_fname) [string]
- Additional GTI file(s) to AND with both the PIN NXB GTI(s) and cleaned
event GTI(s). If more than one file is specified, they will be ORed together
before the AND operation. The same syntax is used as for the
cl_event_fname
parameter. There does not have to be 1-1
correspondence between the gti_fname files and the cl_event_fname
event file(s). The files must be FITS GTI files.
- (cxb_fname) [string]
- If 'NONE', then neither a CXB spectrum nor a total background spectrum
will be calculated.
If 'CALC', then a simulated CXB spectrum will be calculated using a
standard model for the CXB with adjustable parameters. The model will be
calculated as:
CXB = cxb_norm
*(E/cxb_norm_ene
)-cxb_phot_indx
*exp(-E/cxb_cut_ene
)
where E is the spectrum energy in keV. The default values for the
cxb_norm
, cxb_norm_ene
, cxb_phot_indx
,
and cxb_cut_ene
parameters are those found by Boldt et al.
(1987). The exposure time in the simulated spectrum will be set to the value
of the cxb_exposure
parameter.
Otherwise this parameter is taken to be a CXB spectrum file, which will be
added (with addspec
) to the extracted NXB spectrum to create the
total background spectrum.
- (cxb_norm) [real]
- CXB normalization in units of photons/cm^2/s/FOV/keV. See
cxb_fname
parameter
- (cxb_norm_ene) [real]
- CXB normalization energy in keV. See
cxb_fname
parameter.
- (cxb_phot_indx) [real]
- Photon index for CXB model. See
cxb_fname
parameter.
- (cxb_cut_ene) [real]
- Cutoff energy for CXB model in keV. See
cxb_fname
parameter.
- (cxb_exposure) [real]
- Exposure (s) to use for CXB model spectrum. If set to -1, then the
exposure of the extracted NXB spectrum will be used. See
cxb_fname
parameter.
- (cxb_randomize) [bool]
- Randomize counts in simulated CXB spectrum?
- (pinflat_rsp) [string]
- PIN response matrix for flat emission. If set to 'CALDB', then the
calibration database will be queried for the correct file.
- (pinnom_rsp) [string]
- PIN nominal pointing response matrix. The response matrix file
corresponding to the nominal pointing and epoch of the input observation
should be used. The nominal pointing can be checked by examining the
'NOM_PNT' FITS keyword in the input PIN cleaned event file. See the URL:
http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/pinepochs.html
for more information on finding the correct PIN response matrix for a given
epoch.
If set to 'CALDB', then the calibration database will be queried for the
correct file.
- (nxb_scale) [real]
- The exposure time in the extracted NXB spectra will be multiplied by this
parameter. Currently, the background rates in the available NXB background
event files are scaled up by 10 to suppress poisson noise. So to get the
actual count rate one must multiply the exposure by 10.0 (which is the default
for
nxb_scale
).
NOTE: It is unlikely that you will need to change this parameter.
- (groupspec) [bool]
- Group output source spectrum?
- (groupmin) [int]
- Minimum number of counts per grouped channel in grouped output source spectrum
- (cleanup) [bool]
- Remove temporary files?
- (chatter) [0,5]
- Terminal chat level. 0 => no chatter; 5 => debugging
- (clobber) [bool]
- Overwrite existing output files?
- (history) [bool]
- Add history keywords?
NOTES
LAST MODIFIED
August 2009