5.3.2 FSP_PROC (McTiernan)


next up previous contents
Next: 5.3.3 SPFD HXS Up: 5.3 How to Fit Previous: 5.3.1 HXS_FSP (McTiernan)

5.3.2 FSP_PROC (McTiernan)

 

If you have set the keyword countfile in a call to one of the *_FSP routines, the output is a file which can be used as input to the routine FSP_PROC. You may want to do this if you want to save the count rates from a given instrument, and run multiple fits, using different kinds of spectra, and different sets of channels, for example,
IDL> hxs_fsp, index, data, countfile='test.dat'

IDL> fsp_proc, 'HXS', 'test.dat', 'test.otp', fit_pars, sc_par, ch_dta
will result in the fit of the file ``test.dat'' and put the output into the file ``test.otp''. The structures are the same, except that the interval and background times are in a structure which is not compatible with the Yohkoh time structures; you can change the times to the seven element array using the routine TFSPTOT7 (it is possible that the default time structures will be changed to the Yohkoh format in the future).
IDL> n_fits = N_ELEMENTS(fit_pars)

IDL> t_new = intarr(n_fits, 7)

IDL> FOR j = 0, n_fits - 1 DO t_new(j, *) = tfsptot7(fit_pars(j).t)
This will not be necessary in future versions.


a

ydac@mssly1.mssl.ucl.ac.uk