Документ взят из кэша поисковой машины. Адрес оригинального документа : http://xmm.vilspa.esa.es/sas/5.4.1/doc/flspec/node8.html
Дата изменения: Fri Jan 10 23:11:55 2003
Дата индексирования: Sat Dec 22 07:47:47 2007
Кодировка:

XMM-Newton Science Analysis System


flspec (flspec-2.32.8) [xmmsas_20030110_1802-5.4.1]

Meta Index / Home Page

Algorithm

flspec is a perl script (using PEDAL) which does the following:

1. Read input parameters

2. Run region on the observation source list to obtain a FITS file with the regions to be excluded

3. Run the F90 module createpseudoevt to create a pseudo event list with one event in all the X/Y and DETX/Y positions of the input image:

      module createpseudoevt
   
      read input image and source regions file
   	
      create output pseudo event file with X,Y and THETA (offaxis angle) columns
      fill in X(Y) column with X(Y) position of input image pixels
   
      calculate RA,DEC for each pixel
   
      calculate THETA for each pixel from RA,DEC values
   	
      fill in THETA column
   
      end module createpseudoevt
   

4. Filter pseudo event list with output of region and create mask with sources position.

5. Cut exposure map borders with the output mask of emask

6. Run the FTOOLS task farith to flatfield the input image with the exposure map and mask source positions.

7. Run the F90 module flmask to create masks selecting areas with intensity above and below the median:

      module flmask
   
      Read input flatfielded (source masked) image and parameters
   
      Rebin the image
   
      Compute median intensity of rebinned image 
   
      Define positive and negative rebinned pixels above and below the median
   
      Create and write out image mask for positive regions 
   
      Create and write out image mask for negative regions 
   
      end module flmask
   

8. Run the FTOOLS task fstatistic to obtain maximum value for energy channels from the input event list.

9. Filter pseudo event list with annular regions to create a mask

10. Combine positive/negative masks with the annular regions mask

11. Run evselect to filter input event list with the positive combined mask and extract the spectra of positive fluctuations.

12. Calculate solid angle of extraction area.

13. Run evselect to filter input event list with the negative combined mask and extract the spectra of negative fluctuations.

14. Calculate solid angle of extraction area.

15. Set EXPOSURE keyword (equal to 1.0) in positive and negative spectra in case evselect cannot do it properly.

16. Run the FTOOLS task mathpha to subtract both sets of spectra (one for each annular region)

17. Accumulate the spectra of the fluctuations into a single FITS spectral file, by using the task specadd once for each annular region