[Previous]
[Next]
Attempts correction for rear decimation by comparing the count rate
in the decimated times to the count rate in non-decimated times.
CALLS:
CALLED BY
hsi_obs_summ_flag__define, hsi_obs_summ_rate__define
[Previous]
[Next]
PROJECT:
HESSI
NAME:
Open Packet File
PURPOSE:
Creates and/or opens a HESSI FITS packet file. If the file
exists, it is opened (for reading) unless the keyword CREATE
is set. If the file does not exist, the a new file is created.
CATEGORY:
Utilities (hessi/util)
CALLING SEQUENCE:
HSI_PacketFileOpen, unit, filename, nPacket
CALLS: ***
DOC_MENU, FILE_EXIST [2], FXBOPEN [1], FXBOPEN [2], FXBOPEN [3]
HSI_FITSBINARYTABLECREATE, HSI_FITSPRIMARYHDUCREATE, HSI_PACKETFILEOPEN
HSI_STRUCT_DEF [1], HSI_STRUCT_DEF [2], SXPAR [1], SXPAR [2], SXPAR [3]
file_exist [1], file_exist [3]
INPUTS:
filename: The name of the FITS file containing the packets
nPacket: The number of packets that will be written by
HSI_PacketWrite.
OPTIONAL INPUTS:
OUTPUTS:
unit: the unit number of the opened file
columns: a structure containing the column numbers required by
the packet writer HSI_PacketWrite
nPacket: in case the file is opened for reading, nPackets
contains the number of packets stored in the file
OPTIONAL OUTPUTS:
none
KEYWORDS:
CREATE: if set, a new file is created, even if a file with the
same name as "filename" already exists.
ERRMSG: contains the eventual error messages in case of
abnormal termination
FILEINFO: the result of an fstat on the file unit
VERBOSE: if set, prints on the tty what it is doing
SIDE EFFECTS:
none
RESTRICTIONS:
none
SEE ALSO:
HSI_PacketFileClose, HSI_PacketRead, HSI_PacketWrite
MODIFICATION HISTORY:
Release 1, csillag@ssl.berkeley.edu, Dec 1998
Minor changes
Release 0, A.Csillaghy, csillag@ssl.berkeley.edu, Oct.1998
[Previous]
[Next]
PROJECT:
HESSI
NAME:
Open Score File
PURPOSE:
Creates and/or opens a HESSI FITS score file. If the file
exists, it is opened (for reading) unless the keyword CREATE
is set. If the file does not exist, a new file is created.
CATEGORY:
Utilities (hessi/util)
CALLING SEQUENCE:
HSI_ScoreFileOpen, unit, filename, scoreLength
CALLS: ***
DOC_MENU, FILE_EXIST [2], FXBOPEN [1], FXBOPEN [2], FXBOPEN [3]
HSI_FITSBINARYTABLECREATE, HSI_FITSPRIMARYHDUCREATE, HSI_SCOREFILEOPEN
HSI_STRUCT_DEF [1], HSI_STRUCT_DEF [2], SXPAR [1], SXPAR [2], SXPAR [3]
file_exist [1], file_exist [3]
INPUTS:
filename: The name of the FITS file containing the packets
scoreLength: The number of events that will be written by
HSI_ScoreWrite.
OPTIONAL INPUTS:
OUTPUTS:
unit: the unit number of the opened file
columns: a structure containing the column numbers required by
the score writer HSI_ScoreWrite
scoreLength: in case the file is opened for reading, scoreLength
contains the number of packets stored in the file
OPTIONAL OUTPUTS:
none
KEYWORDS:
CREATE: if set, a new file is created, even if a file with the
same name as "filename" already exists.
ERRMSG: contains the eventual error messages in case of
abnormal termination
FILEINFO: the result of an fstat on the file unit
VERBOSE: if set, prints on the tty what it is doing
SIDE EFFECTS:
none
RESTRICTIONS:
none
PROCEDURE:
none
MODIFICATION HISTORY:
Version 1, A.Csillaghy, csillag@ssl.berkeley.edu, Dec.1998
for Release 1 dev
[Previous]
[Next]
PROJECT:
SDAC
;
NAME: other_filters
PURPOSE: the transmission fraction of xrays through pmma and/or multilayer vs energy
CALLING SEQUENCE: result = other_filters( energy, /pmma ) or other_filters(energy,/multi)
INPUTS:
energy - in keV above 0.5 keV
KEYWORD INPUTS:
PANKOW - If set, then use Pankow MLI design of 22-0ct-1997 quoted
by Paul Turin: Dave Pankow tells me he thinks
the plan is more like 5 layers on the grids and 10 in the cryostat of 1/2
mil mylar, only used with the multilayer keyword.
/pmma - transmission through the pmma specified by Reid Brennen
/multi- transmission throught the multilayer specified by Reid Brennen
if both are requested, transmission is the product.
mikechoi - multilayer models specified for grids by Michael Choi, GSFC, code 724.3
-1 - 9.33 mils total Kapton, changed from 23.6/4 mils, ras, 10-aug-95
35200 Angstrom of Al
0 - 10 layers MLI fore and aft
1 - 15 layers
2 - 20 layers
OUTPUTS: the transmitted fraction, from 0 to 1.0
OPTIONAL OUTPUTS:
CALLED BY:
BATSE_SPEC_MODEL, FILTER_ATTEN, ROUTINE_NAME [1], get_slider, low_e_spec_eff [1]
low_e_spec_eff [2], out_spectra_4_designer, set_field
PROCEDURE: interpolates on a data array
CALLS: ***
FCHECK, INTERPOL, filter_cmpnd, xsec
COMMON BLOCKS:
RESTRICTIONS:
MODIFICATION HISTORY: ras, 22-feb-95
ras, 27-feb-95 added multilayer
ras, 19 apr added mikechoi models
ras, 10 aug 95, corrected multilayer inside and outside cryostat
richard.schwartz@gsfc.nasa.gov, 22-oct-1997, added pankow keyword
[Previous]
[Next]
NAME: out_spectra_4_designer
PURPOSE:
Compute the count-rate spectra for the model parameters in
hesi_windows
CALLING SEQUENCE:
out_spectra_4_designer, e_in, photon_in, e_out, count_rate
INPUTS:
through common blocks
OPTIONAL INPUTS:
OUTPUTS:
e_in - 2xn energy array for incident spectra
photon_in - 3xn photon spectra for hesi_windows models
e_out - 2xm energy array for output count rate spectra
count_rate - 3xm count rate spectra for hesi_windows models
OPTIONAL OUTPUTS:
CALLED BY:
ROUTINE_NAME [1], get_slider, set_field
PROCEDURE:
gets the parameter values from the sliders in the widget program
CALLS: ***
CW_ALOGSLIDER, CW_BGROUP, CW_FIELD, CW_FSLIDER, DATPLOT, DETECTOR_FIELD, FCHECK
FILTER_ATTEN, HESI_WINDOWS, MINMAX [1], MINMAX [2], PLOTDATE, PSPLOT [1], PSPLOT [2]
STATUS_FIELD, TOP_BASE_EVENT, XMANAGER, edge_products, get_slider, linecolors [1]
linecolors [2], mewe_spec [1], mewe_spec [2], other_filters [1]
other_filters [2], set_field, xsec
RESTRICTIONS:
only uses the valid energy range of xsec, shouldn't be used >50 keV
for mewe_spec, this must be run using the Yohkoh path!
MODIFICATION HISTORY:
ras, 23-jan-94
[Previous]
[Next]
PROJECT:
HESSI
NAME:
outlier
PURPOSE:
Given a 1-D array of values, y, OUTLIER returns the indices of values that could be considered
as outliers. If there are none, OUTLIER returns a -1.
Values of y do not need to be ordered.
Method is to identify values corresponding to the 25th and 75th percentiles, interpret the
difference as corresonding to 4/3 sigma of the central gaussian fit, and then flagging
values that are more than NSIG sigma from the median as outliers.
As an alternative, see the routine resistant_mean.
CATEGORY:
Utility
CALLING SEQUENCE:
indicies = outlier(y)
CALLS:
none
INPUTS:
y = vector of input values
OPTIONAL INPUTS:
none
OUTPUTS:
Vector of indicies, corresponding to input vector values which were more than NSIG
sigma away from their median.
Returns -1 if there were no such values
OPTIONAL OUTPUTS:
none
KEYWORDS:
NSIG = number of standard deviations outside of which value is considered an outlier.
Default = 3.
COMMON BLOCKS:
none
SIDE EFFECTS:
none
RESTRICTIONS:
MODIFICATION HISTORY:
11-Apr-02 Initial version (ghurford@ssl.berkeley.edu)
8-May-02 gh Added ssw wrapper.