|
Документ взят из кэша поисковой машины. Адрес
оригинального документа
: http://xmm.vilspa.esa.es/sas/8.0.0/doc/epframes/node14.html
Дата изменения: Wed Jul 2 03:38:31 2008 Дата индексирования: Fri Sep 5 20:48:05 2008 Кодировка: Поисковые слова: п п п п п п п п п п п п п п п п п п п п п |
subroutine epframes
* Read input parameters:
*** Open the ODF in directory SAS_ODF
*** (can also be defined via --odfdir or -o; see taskmain)
call OAL_odfInfo
*** Print infos about the proposal
call OAL_proposalInfo
*** Specify which exposure and CCD one is interested in
call OAL_expandFileName(shortfileName,fileName)
OAL_selectScope(EPN, IMAGING) or TIMING or BURST
call OAL_setState(fileName)
*** Print Infos about the exposure
call OAL_exposureInfo
*** Select Data from single CCD
call OAL_selectFile(ccdNr,fileName)
* Declarations for the DAL:
*** Create empty output set and copy attributes
outfile = dataSet(stringParameter("outfile"), CREATE)
call setAttributes(attributable(outfile),attributable(infile))
*** Write standard primary header keywords
call setAttribute(outfile,"TELESCOP","XMM",...)
call setAttribute(outfile,"INSTRUME","EPN",...)
call setAttribute(outfile,"DETNAM"," ",...)
call setAttribute(outfile,"FILTER",filterId(OAL_activeFilter()),...)
if (odfInfos%isSlew) then
call setAttribute(outfile,"OBS_MODE","SLEW",...)
else
call setAttribute(outfile,"OBS_MODE","POINTING",...)
endif
call setAttribute(outfile,"OBS_ID", &
odfinfos%proposalId+odfinfos%observationId,...)
call setAttribute(outfile,"EXP_ID",expinfos%number,...)
select case (readoutMode)
case(0)
call setAttribute(outfile,"CONTENT","EPIC IMAGING MODE EVENT LIST",...)
call setAttribute(outfile,"DATAMODE","IMAGING",...)
case(1)
call setAttribute(outfile,"CONTENT","EPIC TIMING MODE EVENT LIST",...)
call setAttribute(outfile,"DATAMODE","TIMING",...)
case(8)
call setAttribute(outfile,"CONTENT","EPIC BURST MODE EVENT LIST",...)
call setAttribute(outfile,"DATAMODE","BURST",...)
end select
call setAttribute(outfile,"ORIGIN","Leicester/SSC",...)
call setAttribute(outfile,"OBJECT",proposalinfo%targetName,...)
call setAttribute(outfile,"DATE-OBS",odfinfos%observationStartTime,...)
call setAttribute(outfile,"DATE-END",odfinfos%observationEndTime,...)
call setAttribute(outfile,"OBSERVER",proposalinfo%GoName,...)
*** Write some additional keywords (to be used by attcalc)
call setAttribute(outfile,"RA_OBJ",proposalinfo%targetRa,...)
call setAttribute(outfile,"DEC_OBJ",proposalinfo%targetDec,...)
call setAttribute(outfile,"RA_NOM",proposalinfo%boresightRa,...)
call setAttribute(outfile,"DEC_NOM",proposalinfo%boresightDec,...)
*** Create events table extension (one per CCD)
outccd = addTable(outfile,"EVENTS",ntrow)
call setAttributes(attributable(outccd),attributable(inccd))
call setAttribute(outccd,"TSTART" ,...) ! start time of first frame
call setAttribute(outccd,"TSTOP" ,...) ! end time of last frame
call setAttribute(outccd,"TELAPSE",...) ! TSTOP - TSTART in [s]
call setEventAttribute(outccd,EVATT_IN_SPOILED_FRAME)
call setEventAttribute(outccd,EVATT_IN_BAD_FRAME)
call setEventAttribute(outccd,EVATT_OUT_OF_CCD_WINDOW)
*** Read table columns FRAME, RAWX, RAWY, ENERGY
* Convert frame counter to event time tags:
col = column(inccd, "FRAME", READ)
frames => int32Data(col)
call OAL_frameCounterToObt
call OAL_obtToTimeTag
* fine time correction:
if IMAGING mode then
call CAL_correctTime(timeTags)
else ! TIMING or BURST with source position
call CAL_correctTime(timeTags, rawY, srcposition)
endif
* Compute mode dependent livetime factor FC0:
* Create exposure extension
TIME = real64Data(column(inccd,"TIME",READ))
TIMEDEL = real64Attribute(inccd,"FRMTIME")/1000 * FC0
FRACEXP = 1 (real32Data, other livetime factors to be handled
by higher SAS task)
* Correct for offset shifts
call CAL_offsetCorrect
*** determine frame boundaries and prepare the loop
*** read table column RAWX, RAWY and FLAG frame by frame
do frame = firstFrame, lastFrame
ixx = RAWX(frame)
iyy = RAWY(frame)
iqual = FLAG(frame)
itime = TIME(frame)
* Search for valid events
*** Check for wrong coordinates
read window size and position
if (ixx,iyy) outside this window, than
call warning "Error in coordinates"
qual(i) = ior(qual(i), EVATT_OUT_OF_CCD_WINDOW)
! SINGLE EVENT FLAGGING:
! OUTSIDE THE READOUT WINDOW
*** Check for correct order of events
pix_code_new = iyy * 100 + ixx
if (pix_code_new <= pix_code_old) then
call warning "Events are in wrong order"
qual(i) = ior(qual(i), EVATT_IN_BAD_FRAME)
! WHOLE FRAME FLAGGING
qual(i) = ior(qual(i), EVATT_IN_SPOILED_FRAME)
! WRONG ORDER OF EVENTS
endif
pix_code_old = pix_code_new
*** Check for time order
if frame not sorted in time
call warning "Frame not sorted in time"
qual(i) = ior(qual(i), EVATT_IN_BAD_FRAME)
! WHOLE FRAME FLAGGING
qual(i) = ior(qual(i), EVATT_IN_SPOILED_FRAME)
! NOT SORTED IN TIME
enddo ! loop over frames
*** Display statistics about wrong events
wronglimit = intParameter("wrongpixlimit")
if (number_of_wrong_events > real(wronglimit)) &
call warning "More wrong events than wrongpixlimit allows"
*** Add statistics to the history and keywords
*** Convert EPIC PN ODF file into format defined in ICD
write outfile with columns
TIME, FRAME, RAWX, RAWY, PHA, CCDNR, FLAG
NOTE: RAWX/Y and attributes WINDOWX0/Y0 are increased by 1
* Write GTIs depending on start/stop time of observation
and wrong/missing frames
*** Create empty output set and copy primary standard attributes
outgti = dataSet(stringParameter("outgti"), CREATE)
call setAttributes(attributable(outgti),attributable(outfile))
*** Create GTI extension
outgticcd = addTable(outgti,"STDGTIxx",ngti) ! where xx is the overall
! CCD number (01-12)
call setAttribute(outgticcd,"HDUCLASS","OGIP",...)
call setAttribute(outgticcd,"HDUCLAS1","GTI",...)
call addAttributes(attributable(outgticcd),attributable(outccd))
for each GTI:
start(i) = time(firstEvent) - (FRMTIME/1000)/2
stop(i) = time(lastEvent) + (FRMTIME/1000)/2
end subroutine epframes