Script for Step 2

; ospex script created Sat Sep 04 05:54:19 2004
;

;Initialize ospex object and call it "obj."
obj = ospex(/no_gui) ;Type obj->gui later to get the GUI.

set_logenv,'OSPEX_NOINTERACTIVE','1' ;make everything non-interactive

;Set input spectrum data file generated in Step 1 for detector #3.
obj-> set, $
spex_specfile= 'C:\My_Documents\My_HESSI\2002\04\21\Orbit1\spectra\OSPEX\sp_20024210280_to_1360_D3.fits'

;Set srm file generated in Step 1
obj-> set, $
spex_drmfile= 'C:\MY_DOCUMENTS\MY_HESSI\2002\04\21\ORBIT1\SPECTRA\OSPEX\srm_20024210280_to_1360_D3.fits'

;Set no separate background for each energy band
obj -> set, spex_bk_sep = 0

;Set background order of polynomial, and time intervals
obj -> set, spex_bk_order=0
obj->set, spex_bk_time_interval=[ ['21-Apr-2002 00:29:00.000', '21-Apr-2002 00:31:00.000']]

;Set time intervals for spectral fits
obj-> set, spex_fit_time_interval= [['21-Apr-2002 00:32:15.999', $
'21-Apr-2002 00:33:15.999'], ['21-Apr-2002 00:33:15.999', '21-Apr-2002 00:34:15.999'], $
['21-Apr-2002 00:34:15.999', '21-Apr-2002 00:35:15.999'], ['21-Apr-2002 00:35:15.999', $
'21-Apr-2002 00:36:16.000'], ['21-Apr-2002 00:36:16.000', '21-Apr-2002 00:37:16.000'], $
['21-Apr-2002 00:37:16.000', '21-Apr-2002 00:38:15.999'], ['21-Apr-2002 00:38:15.999', $
'21-Apr-2002 00:39:15.999'], ['21-Apr-2002 00:39:15.999', '21-Apr-2002 00:40:15.999'], $
['21-Apr-2002 00:40:15.999', '21-Apr-2002 00:41:15.999'], ['21-Apr-2002 00:41:15.999', $
'21-Apr-2002 00:42:16.000'], ['21-Apr-2002 00:42:16.000', '21-Apr-2002 00:43:16.000'], $
['21-Apr-2002 00:43:16.000', '21-Apr-2002 00:44:15.999'], ['21-Apr-2002 00:44:15.999', $
'21-Apr-2002 00:45:15.999'], ['21-Apr-2002 00:45:15.999', '21-Apr-2002 00:46:15.999'], $
['21-Apr-2002 00:46:15.999', '21-Apr-2002 00:47:15.999'], ['21-Apr-2002 00:47:15.999', $
'21-Apr-2002 00:48:15.999'], ['21-Apr-2002 00:48:15.999', '21-Apr-2002 00:49:16.000'], $
['21-Apr-2002 00:49:16.000', '21-Apr-2002 00:50:15.999'], ['21-Apr-2002 00:50:15.999', $
'21-Apr-2002 00:51:15.999'], ['21-Apr-2002 00:51:15.999', '21-Apr-2002 00:52:15.999'], $
['21-Apr-2002 00:52:15.999', '21-Apr-2002 00:53:15.999'], ['21-Apr-2002 00:53:15.999', $
'21-Apr-2002 00:54:15.999'], ['21-Apr-2002 00:54:15.999', '21-Apr-2002 00:55:16.000'], $
['21-Apr-2002 00:55:16.000', '21-Apr-2002 00:56:16.000'], ['21-Apr-2002 00:56:16.000', $
'21-Apr-2002 00:57:15.999'], ['21-Apr-2002 00:57:15.999', '21-Apr-2002 00:58:15.999'], $
['21-Apr-2002 00:58:15.999', '21-Apr-2002 00:59:15.999'], ['21-Apr-2002 00:59:15.999', $
'21-Apr-2002 01:00:15.999'], ['21-Apr-2002 01:00:15.999', '21-Apr-2002 01:01:15.999'], $
['21-Apr-2002 01:01:15.999', '21-Apr-2002 01:02:16.000'], ['21-Apr-2002 01:02:16.000', $
'21-Apr-2002 01:03:15.999'], ['21-Apr-2002 01:03:15.999', '21-Apr-2002 01:04:15.999'], $
['21-Apr-2002 01:04:15.999', '21-Apr-2002 01:05:15.999'], ['21-Apr-2002 01:05:15.999', $
'21-Apr-2002 01:06:15.999'], ['21-Apr-2002 01:06:15.999', '21-Apr-2002 01:07:15.999'], $
['21-Apr-2002 01:07:15.999', '21-Apr-2002 01:08:16.000'], ['21-Apr-2002 01:08:16.000', $
'21-Apr-2002 01:09:16.000'], ['21-Apr-2002 01:09:16.000', '21-Apr-2002 01:10:15.999'], $
['21-Apr-2002 01:10:15.999', '21-Apr-2002 01:11:15.999'], ['21-Apr-2002 01:11:15.999', $
'21-Apr-2002 01:12:15.999'], ['21-Apr-2002 01:12:15.999', '21-Apr-2002 01:13:15.999'], $
['21-Apr-2002 01:13:15.999', '21-Apr-2002 01:14:15.999'], ['21-Apr-2002 01:14:15.999', $
'21-Apr-2002 01:15:16.000'], ['21-Apr-2002 01:15:16.000', '21-Apr-2002 01:16:15.999'], $
['21-Apr-2002 01:16:15.999', '21-Apr-2002 01:17:15.999'], ['21-Apr-2002 01:17:15.999', $
'21-Apr-2002 01:18:15.999'], ['21-Apr-2002 01:18:15.999', '21-Apr-2002 01:19:15.999'], $
['21-Apr-2002 01:19:15.999', '21-Apr-2002 01:20:15.999'], ['21-Apr-2002 01:20:15.999', $
'21-Apr-2002 01:21:16.000'], ['21-Apr-2002 01:21:16.000', '21-Apr-2002 01:22:16.000'], $
['21-Apr-2002 01:22:16.000', '21-Apr-2002 01:23:15.999'], ['21-Apr-2002 01:23:15.999', $
'21-Apr-2002 01:24:15.999'], ['21-Apr-2002 01:24:15.999', '21-Apr-2002 01:25:15.999'], $
['21-Apr-2002 01:25:15.999', '21-Apr-2002 01:26:15.999'], ['21-Apr-2002 01:26:15.999', $
'21-Apr-2002 01:27:16.000'], ['21-Apr-2002 01:27:16.000', '21-Apr-2002 01:28:16.000'], $
['21-Apr-2002 01:28:16.000', '21-Apr-2002 01:29:15.999'], ['21-Apr-2002 01:29:15.999', $
'21-Apr-2002 01:30:15.999'], ['21-Apr-2002 01:30:15.999', '21-Apr-2002 01:31:15.999'], $
['21-Apr-2002 01:31:15.999', '21-Apr-2002 01:32:15.999']]

;Set functions to use for spectral fits
obj-> set, fit_function= 'vth_noline+line+line+line+3pow'

;Set starting, minimum, and maximum values of function parameters  for fit functions.
obj-> set, fit_comp_params= [1.02536, 2.48283, 246002., 6.98800, 0.100000, 28438.6, $
8.34892, 0.100000, 701.648, 10.9477, 0.100000, 10.6611, 1.50000, 18.5944, 4.43239, $
400.000, 2.00000]
obj-> set, fit_comp_minima= [1.00000e-020, 0.500000, 1.00000e-020, 1.00000, 0.0100000, $
1.00000e-020, 1.00000, 0.0100000, 1.00000e-020, 1.00000, 0.0100000, 1.00000e-020, $
0.500000, 1.00000, 0.500000, 10.0000, 0.500000]
obj-> set, fit_comp_maxima= [1.00000e+020, 50.0000, 1.00000e+020, 1.00000e+006, $
1.00000e+006, 1.00000e+020, 1.00000e+006, 1.00000e+006, 1.00000e+020, 1.00000e+006, $
1.00000e+006, 1.00000e+020, 12.0000, 20000.0, 12.0000, 20000.0, 12.0000]

;Set mask to all paramters to be fixed (0) or free(1).
obj-> set, fit_comp_free_mask= [1.00000, 1.00000, 1.00000, 1.00000, 0.000000, 1.00000, $
1.00000, 0.000000, 1.00000, 1.00000, 0.000000, 1.00000, 0.000000, 1.00000, 1.00000, $
0.000000, 0.000000]

;Set loop mode to forwards
obj-> set, spex_fit_reverse= 0

;Set energy ranges to include in fit.
obj-> set, spex_erange= [[6.0000000D, 20.000000D]]

;Set units of plots.
obj-> set, spex_autoplot_units= 'Flux'

;Set energy bands to use for time  history plots.
obj-> set, spex_eband= [[3.00000, 6.00000], [6.00000, 12.0000], [12.0000, 25.0000], $
[25.0000, 50.0000], [50.0000, 100.000]]

;Restore ospex results file from previous run for detector #1 to act as starting points for detector #3
obj -> restorefit, $
file='C:\My_Documents\My_HESSI\2002\04\21\Orbit1\spectra\OSPEX\ospex_results_d1.geny'


end