NAME

pcparpha -- generates a Plucinsky etal background PHA spectrum for the ROSAT PSPC.


USAGE

pcparpha


DESCRIPTION

IMPORTANT: If you are using US Rev0 data to create PHA files, please see KNOWN BUGS below. There is a serious error condition which you must work around.

This task generates a 256-channel PHA spectrum of the particle background in the ROSAT PSPC from the parameterization of Plucinsky etal (1993 ApJ 418, 519)

The particle background spectrum of the PSPC is composed of three components (Plucinsky etal):

The Internal Component
This is present at all times when the PSPC is turned on, and is produced by charged particle events within the detector. Its spectrum shows a steep decrease between PI channels 18-40, but then flattens above channel 40. Spatially, this component exhibits a radial gradient whereby its normalization is 19% smaller at the edge of the detector (1 degree off-axis) compared to at the center.

The Aluminium Component
This component is only present when the filter wheel is open, and most arises from fluorescent Al Kalpha within the instrument housing. Its spectrum is that of an 1.5keV Al Kalpha line, broadened by the spectral resolution of the detector. This component shows not positional dependence except for the shadows produced by the window support structure.

The External Component
This component is only present when the filter wheel is open. Its spectrum is completely flat (ie constant count rate per channel). This component shows not positional dependence beyond that expected from the shadowing produced by the window support structure, and the decrease in unobscured open-window area as one moves off-axis.

As shown by Plucinsky etal, the count rate expected from all three 3 components correlate well with the Master Veto (MV) count rate.

The basic input parameters to the task are therefore the MV dataset corresponding to an observation, and user-defined information regarding the temporal & spatial parameters for which the background is to be generated. Any temporal filters to be applied (to the MV dataset) can be specified via an Good Time Interval (GTI) dataset. The spatial information can either be specified via an off-axis histogram dataset (which also contains the spatial extent of the region for which the background spectrum is to be generated), or specified via parameters. Users are also able to include/exclude any of the above components from the computation.

The task then uses the parameterizations given in Plucinsky etal (1993) to construct the requested spectrum. The rejection efficiency of the PSPC for particle background events for PI channels 18-249 is very high (99.90%) with a typical count rate of 4x10E-6 count/s/arcmin**2 (Plucinshky etal).

Only the the OGIP recognised PHA file formats described in Arnaud etal 1992 (Legacy, 2, 65) and its appendix provided in OGIP memo OGIP/92-007a are supported for both the input and output files.


WARNINGS ON USAGE

Please report any comments/problems or bugs to Ian M George (http://heasarc.gsfc.nasa.gov/cgi-bin/ftoolshelp).


PARAMETERS

outfil [character string]
The name of the o/p file to be written containing the derived PHA dataset. For safety, the name of the o/p file CANNOT be the same as one of the i/p files, and the task will stop if this is the case. Under unix/ultrix, an o/p filename identical to a file which already exists on disk is also considered illegal, and the task will stop. However, the existing file can be automatically removed, and the new file written if the o/p filename is preceeded by "!" at the outfil prompt.

(qintbgd = T) [boolean]
A logical flag indicating whether the INTERNAL background component of Plucinsky etal is to be included in the o/p file.

(qextbgd = T) [boolean]
A logical flag indicating whether the EXTERNAL background component of Plucinsky etal is to be included in the o/p file.

(qalbgd = T) [boolean]
A logical flag indicating whether the ALUMINUM background component of Plucinsky etal is to be included in the o/p file.

evrfil [character string]
The name of the file containing the Master Veto (MV) dataset.

(evrbin = 30) [integer]
This parameter specifies the bin size (in seconds) over which the mean MV rate is computed & used to calculate the background spectrum. This prcedure both speeds up the task, and smooths out the 'jitter' in MV measurements. The default value of evrbin=30 is suitable for most applications.

gtifil [character string]
The name of the file containing the Good Time Internal (GTI) dataset. The o/p PHA dataset will be constructed using the MV rate within the GTIs (only). Alternatively, the special value, gtifil = 'NONE', will result in the entire MV dataset being used (equivalent of a single GTI starting & ending at the times of the first & last measurement of the MV rate in the evrfil, respectively).

oahfil [character string]
This parameter is only required if qintbgd=T, and contains the name of the file containing the off-axis histogram dataset (fraction of time spent as a function of off-axis angle). The special value oahfil='NONE' will result in a single off-axis angle being prompted for via the theta parameter (below).

(theta = 0) [real]
This parameter is only required if qintbgd=T & oahfil='NONE', and contains the single off-axis angle (in arcmins) at which the user wishes the background PHA spectrum to be calculated for.

npixsou = % [pseudo-integer]
This is the number of pixels (of size specified by the pixsize parameter below) in the 'source region' for which the background PHA spectrum is to be calculated.

If an off-axis histogram dataset is specified via the
oahfil parameter, then the special value npixsou='%' can be used to indicate that the number of pixels is to be read from the NPIXSOU keyword of the off-axis histogram extension. If an integer value is entered, then this value will OVERRIDE the value read from the NPIXSOU keyword.

If an off-axis histogram dataset is NOT specified via the
oahfil parameter, an integer value must be specified.

pixsize = 0.00833 [real]
This parameter is only required if the value of the npixsou parameter is NOT '%', and specifies the length (in arcmin) of the side of a pixel.

(phaversn = '1.1.0') [character string]
The OGIP FITS format version for PHA datasets in which the o/p FITS file is to be written. The default is PHAVERS1='1.1.0' (also known as known as PHAVERSN = '1992a'), and is currently the only format allowed.

(chatter = 10) [integer]
Flag to indicate how chatty the task is at execution. A value of 10 is the default, with lower/higher values producing quieter/verbose output respectively.

(optaxisx,optaxisy)
Do not remove cover. No user serviceable parts inside.


KNOWN BUGS

IMPORTANT! US Rev0 format data (but not RDF data or German data which has been converted with GEVT2RDF) does not properly specify the size of the pixels in the Detector image. Because of this, the extractor (the program that xselect uses to make a spectrum) puts the wrong pixel size into the PHA file. The extractor will be modified to cope with this deficiency in US Rev0 data but it isn't fixed yet. So if you use US Rev0 data and you use xselect to creat a PHA (spectrum) file, you must change the DELTX and DELTY parameters to get a correct answer.

You need to change the CDELT1 and CDELT2 keywords to 2.595021e-4 times whatever the value of your WMREBIN keyword is (usually 15). You can use FKEYPRINT to find the value of WMREBIN and FPARKEY to put the new value into CDELT1 and CDELT2. If you have perl installed at your site (and you probably do since ascascreen is written in perl and if you don't, bug your sysadmin to install it) the script fixrev0pha will do these changes for you (type 'fhelp fixrev0pha' to see the help page).

We apologize for this flakiness, but it is a prime example of why we needed the RDF format in the first place.

This task should still be considered as a BETA test version. Please perform spot checks and report all potential bugs to http://heasarc.gsfc.nasa.gov/cgi-bin/ftoolshelp.


SEE ALSO

Plucinsky, P.P., Snowden, S.L., Briel, U.G., Hasinger, G. & Pferrermann, E.
1993, ApJ, 418, 519.

Snowden, S.L., Plucinsky, P.P., Briel, U.G., Hasinger, G. & Pferrermann, E.
1992, ApJ, 393, 819.


LOG OF SIGNIFICANT CHANGES

v2.1.0 (1994 August)
WMAP support

v2.0.0 (1994 April)
Public release version

v1.0.0 (1994 Feb)
Pre-release version (called PSPCBKGD)


PRIMARY AUTHOR

Ian M George

HEASARC

NASA/GFSC

http://heasarc.gsfc.nasa.gov/cgi-bin/ftoolshelp

(301) 286-6094

CATEGORY

Mar94 ftools.rosat