Документ взят из кэша поисковой машины. Адрес оригинального документа : http://xmm.vilspa.esa.es/sas/7.0.0/doc/dpssflag/node8.html
Дата изменения: Thu Jun 29 00:04:28 2006
Дата индексирования: Sat Dec 22 09:37:18 2007
Кодировка:

XMM-Newton Science Analysis System


dpssflag (dpssflag-7.0.2) [xmmsas_20060628_1801-7.0.0]

Meta Index / Home Page

Algorithm

Setting flag 1: MASKFRAC value <= mfraclimit

        Initialize flags to FFFFFFFFF for all sources
        Foreach source in list{
           If MASKFRAC (value <=mfraclimit) then FLAG[1:1](i)=T
        }
	# flag also summary row
	if(ID_INST(i) == 0 and FLAG[1:1](this source)=T) then FLAG[1:1](i)=T
 
Setting flags 2&4: source close to another bright point source

        Foreach point source i:
           calculate radius of 'circle_of_influence'= brcolumn(i)*brfactor
           if(distance_of_nearest_neighbour < radius(i))then
             Foreach source j/=i :
                  !set Flags 2 (and 4):
                  if(source_j = source_j-1) then
                     FLAG[2:2](j)=FLAG[2:2](j-1)
                     FLAG[4:4](j)=FLAG[4:4](j-1)
                  else 
                     calculate distance(i,j)
                     if(distance(i,j) < radius(i)) then 
                       FLAG[2:2](j)=T
                       if(EXT(j)>=extlimit .AND. SCTS(i)>brlimit) FLAG[4:4](j)='T'
                     endif	
                  endif
             }
           endif
         }

Setting flags 3&5: source in extended emission area of an extended source

         Foreach source i:
           calculate radius of 'circle_of_extension'= ext(i)*extfactor
           if(distance_of_nearest_neighbour < radius(i))then
             Foreach source j/=i :
                 !set flags 3 (and 5):
                 if(source_j = source_j-1) then
                    FLAG[3:3](j)=FLAG[3:3](j-1)
                 else 
                    calculate distance(i,j)
                    if(distance(i,j) < radius(i)) then 
                      FLAG[3:3](j)=T
                      if(EXT(j)>=extlimit .AND. brcolum(j)<brexfraction*brcolum(i) FLAG[5:5](j)='T'
                    endif
                endif
             }
           endif
         }
 

Setting flag 6 : extended source detected in one only band 

         Foreach source i:
            Calculate detmlsum(1-15 images) and detmltot(IDINST=0,IDBAND=0)
            if(EXT(i) >= extlimit)then
               if(IDBAND(i)>0 .AND. IDBAND(i)<=5 .AND. &
                 detml(i)>detmlfraction*detmlsum(ML_ID_SRC(i)))then
                 flag6(ML_ID_SRC(i))="T"
              endif
            endif
        }

Setting flag 7: likely spurious extended source
         Foreach source i:
      	    if(FLAG[4:4]="T" .or. FLAG[5:5]="T" .or. FLAG[6:6]="T") then
		FLAG[7:7]="T"
         }


Setting flag 8/9: source ON/CLOSETO  bad PN(MOS) bands 1(1,2,3) column(s)

        Foreach source i:
          calculate 'critical' distance to a bad column=oper1(badcolumn(i))oper2(badfactor)
          if ((PN and BAND=1) or (M1 and band 1,2,3)) then
             Foreach bad column:
               if(src ON bad colum):
                  FLAG[8:8](i)='T'
               elseif(distance_to_bad_column < critical_dist)then
                  FLAG[9:9](i)='T' (if PN, 'pnclosetobadon'=yes)
               endif
             }
          endif
        }
        if FLAG[8:8](src)=T then FLAG[8:8](summary_row)=''T''
        if FLAG[9:9](src)=T then FLAG[9:9](summary_row)=''T''