NAME

pha_summary -- Run all the Pulse Height Analysis tasks on an observation.

USAGE

pha_summary input eefile intphd difphd fitphd table text diodes

DESCRIPTION

pha_summary is a script that makes running all the individual Pulse Height Analysis (PHA) tasks a simple process. Given an observation input file name and an output file name, this script runs the following tasks in this order: phaz, pha_text, pha_center, pha_rms, pha_sig2rms, pha_diode, and pha_all. The graphics output is sent to the device (usually a printer) specified by the device parameter. The text output is placed in the file specified by the text parameter. (See the help on the individual tasks for a full description of each of their outputs).

PARAMETERS

input [file name template]

Names of files containing the input PHA observation raw data. This data must be sorted in ascending threshold order. A PHA observation will normally have its data in separate groups of a single input file.

eefile [file name template]

Names of the extracted engineering trailer files corresponding to each file or group specified by input. This file contains the data for the 12 special diodes.

intphd = "intphd" [file name]

Name of the first output image, which will contain the integral pulse height distribution. This distribution will be a two-dimensional array of size (number of thresholds observed) x 512. Each row i of the image contains the integral pulse height distribution for the corresponding diode number i. For example, row 50 will contain the integral pulse height distribution for diode number 50.

difphd = "difphd" [file name]

Name of the second output image, which will contain the differential pulse height distribution. This image takes the same form as the image created for intphd.

fitphd = "fitphd" [file name]

Name of the third output image, which will contain the fitted differential pulse height distribution. This image takes the same form as that of intphd and difphd.

table = "phaz" [file name]

Name of the output table that will contain results of the threshold analysis.

text = "phaz.txt" [file name]

The file containing the textual summary produced by the task pha_text.

diodes = "dev$null" [file name]

Text file containing the list of diodes for which detailed plots should be produced (one diode number per line). If blank, the diode numbers are taken from standard input. If this is set to "dev$null", then no plots are produced.

(t= 30.0) [real, min=1.0, max=200.0]

Initial threshold level in input (i.e., the threshold level for the first group of the image).

(delt = 1.0) [real, min=1.0, max=200.0]

Step size for incrementing the threshold levels in input.

(device = "stdgraph") [string]

Output device on which graphs will be printed.

(left = 0.) [real]

The left boundary on the graphics device, specified from 0. to 1., where the graphs will appear.

(right = 1.) [real]

The right boundary on the graphics device, specified from 0. to 1., where the graphs will appear.

(bottom = 0.) [real]

The bottom boundary on the graphics device, specified from 0. to 1., where the graphs will appear.

(top = 1.) [real]

The top boundary on the graphics device, specified from 0. to 1., where the graphs will appear.

(verbose = yes) [real]

Display number of diodes drawn as task progresses?

EXAMPLES

1. Run the pulse height analysis suite of tasks on the observation z06i0501t, using all the defaults.

cl> pha_summary z06i0501t.d0h[*] z06i0501t.x0h[*]

2. Run the summary above but use a different percentage value for the determination of the threshold. Since this is not one of the parameters of pha_summary itself, the parameter to the task phaz itself has to be changed. Change the value to 60.

cl> phaz.percent=60.
cl> pha_summary z06i0501t.d0h[*] z06i0501t.x0h[*]

BUGS

REFERENCES

Author: Jonathan Eisenhamer

HELP

For assistance using this or any other tasks, please contact help@stsci.edu or call the help desk at 410-338-1082.

SEE ALSO

phaz, pha_all, pha_center, pha_diode, pha_rms, pha_sig2rms, pha_text