|
Документ взят из кэша поисковой машины. Адрес
оригинального документа
: http://xmm.vilspa.esa.es/sas/7.1.0/doc/espfilt/node15.html
Дата изменения: Mon Jul 9 04:25:10 2007 Дата индексирования: Sat Dec 22 09:51:00 2007 Кодировка: Поисковые слова: http www.badastronomy.com bad tv foxapollo.html |
subroutine espfilt
Read parameters.
Call method
if (method .eq. corner) then
call corner
subroutine corner:
task evselect (to create object image)
task evselect (to create object light curve)
task evselect (to create corners image)
task evselect (to create corners light curve)
call clean_lc
subroutine clean_lc:
extract columns TIME, RATE from object LC
extract columns TIME, RATE from corners LC
bin TIME, RATE columns from object LC
bin TIME, RATE columns from corners LC
smooth RATE column from object LC
smooth RATE column from corners LC
create histogram from object
find peak of the distibution
call gauss_fit
subroutine gauss_fit:
find best paramater fit using Marquardt routines:
MQM maximum likelyhood fitting routine
MRQ compute the matrix elements for Marquard method
GFUNCS calculate gaussian function distribution
return best fit
call write_qdp_fit
subroutine write_qdp_fit
open QDP output file
write header, title information
write out count rate histogram
overlay selection and fit limits
write out FOV lightcurve
overlay GTIs
write out Corner lightcurve
overlay GTIs
open ASCII gti.txt file
write out GTIs
task gtibuild create GTI fits file
task evselect create object filtered event list
task evselect create object filtered image
task evselect create corners filtered event list
task evselect create corners filtered image
Filtered event list, lightcurves, nice qdp plot are generated.
end espfilt