Документ взят из кэша поисковой машины. Адрес оригинального документа : http://xmm.vilspa.esa.es/sas/5.4.1/doc/emevents/node9.html
Дата изменения: Fri Jan 10 22:52:58 2003
Дата индексирования: Sat Dec 22 09:16:15 2007
Кодировка:
Поисковые слова: asteroid

XMM-Newton Science Analysis System


emevents (emevents-6.5.4) [xmmsas_20030110_1802-5.4.1]

Meta Index / Home Page

Algorithm

MODULE em_events_module

   Read the file names of event files (input and output) and frame file

   Open files 
   Get keyword values   
 
   parameters read : choice of sub_tasks to execute
   Loop over input task parameters
      Read the task parameter / 1 to perform the procedure,
                                0 not to perform it  /
   End loop

   Get the values of files columns in memory
  
SUBROUTINE CUT_BAD 
   if cut_bad requested  then
     read bad pixels in extension
     Loop over frames
       Creation of image using inverse pattern definition 
       Connection with binodal task 
          in the case of 2 files (prime and redundant) analysis
       Set bad pixels to 0
       Loop over events
          Use of EDUSOFT code for new pattern reconstruction
          creation of new rows and destruction of old ones
       end loop
     end loop
   endif
  
SUBROUTINE DIAGO 
   if diago requested then 
     Loop over event rows for diagonal patterns
       Check if both E1 and E2 energies are different from Si K energy
       In this case :
         first correct the existing pattern    
         pattern is 0, e2 is 0, others the same   
         second create another pattern   
         pattern is 0, e1 is old e2, X & Y to be calculated 
         depending on pattern  
         e2 is 0, others the same   
         both patterns are flagged   
         ( do not correct the NVALID in the file frame )   
      end loop
      new rows have been created at the end of file
      now shift of lines in order to insert new events
      Loop over number of new lines
         Shift all rows following place of insertion
      end loop
    endif

    Check events with energy value below EMDHLOW or above EMDHHIGH

SUBROUTINE CUT_GTI 
    if cut_gti requested  then 
      Loop over frame rows
         all events whose frames were flagged are flagged for rejection 
      end loop
    endif

SUBROUTINE PUT_XY 
    if put_xy requested then
      Loop over events
         center the events, substract the underscan, add 1.
         flag the events outside the field of view
      end loop
      Modify window keywords
      Compute camera coordinates
    endif 

SUBROUTINE PUT_TI  ( Imaging mode )
    if put_ti requested then
      Loop over frame rows
         Loop over events rows
           calculate the average time for the frame 
           copy this time in the time column of events file
         end loop
       end loop
     endif
 
SUBROUTINE PUT_TT  ( Timing mode )
     if put_tt requested then
       Loop over frame rows
         Loop over events rows
           calculate the new time as reference frame time plus Y offset
           copy this time in the time column of events file
         end loop
       end loop
     endif

SUBROUTINE REJECT_ROWS
    if reject_rows requested then
       Loop over frame rows
         Loop over events rows
            Build RAWY profile
         end loop
         Loop over events rows
            Flag events with RAWY profile larger than maxeventsperrow
         end loop
       end loop
    endif

SUBROUTINE REJECT_E3
    if reject_e3 requested then
      Read e3threshold parameter for current instrument and CCD
      Loop over events rows
         Flag events with pattern <= 4 and E3 < e3threshold
      end loop
    endif

SUBROUTINE SP_GATTI 
    if sp_gatti requested   then
      Loop over events rows
        calculate the E1 spectrum, the (E1+gatti) spectrum 
        and the gatti spectrum
      end loop
      Convolution of E1 and gatti spectra
    endif

    sequence of instructions for outevent1 creation 
 
    closing files
        
end module