Oberon Version

Ian M George: Make an EVT2 from an EVT1

These pages are unofficial in all respects. They are primarily simply notes to myself. The information is not guaranteed to be correct or up-to-date.

Last Updated: 2004 Jun 20

Starting Point... OK, so you have copied the data, unpacked it, constructed a work area, have set-up the necessary software packages, started taking notes, and have located some of the necessary files etc ... Now we are ready to run some ciao tools

Set the Bad-pixel maps

This step is crucial to let the software know where the badpixel map is for a given obsid.

punlearn ardlib
foreach d ( 4 5 6 7 8 9 )
pset ardlib AXAF_ACIS${d}_BADPIX_FILE = "{bpix1}[BADPIX${d}]"
end

so again, in the example we used earlier, the actual command would be

punlearn ardlib
foreach d ( 4 5 6 7 8 9 )
pset ardlib AXAF_ACIS${d}_BADPIX_FILE = "../primary/acisf02090_001N003_bpix1.fits[BADPIX${d}]"
end

Basic Ciao Processing

So now we are ready to run a whole bunch of ciao commands.
In most cases the output from one is the input to the next... our final goal being an "evt2" file.

punlearn acis_process_events
acis_process_events infile={evt1} outfile=acis_new_evt1.fits acaofffile="@../primary/pcad_asol1.lis" eventdef=")stdlev1" apply_cti=yes clobber=yes

where {evt1} is ../secondary/acisf02090_001N003_evt1.fits in our example.

punlearn tgdetect
tgdetect infile=acis_new_evt1.fits outfile=acis_src1a.fits OBI_srclist_file=NONE clobber=yes

punlearn tg_create_mask
tg_create_mask infile=acis_new_evt1.fits outfile=acis_evt1_L1a.fits input_pos_tab=acis_src1a.fits grating_obs=header_value clobber=yes

punlearn tg_resolve_events
tg_resolve_events infile=acis_new_evt1.fits outfile=acis_evt1a.fits regionfile=acis_evt1_L1a.fits acaofffile="@../primary/pcad_asol1.lis" eventdef=")stdlev1_ACIS" clobber=yes

punlearn dmcopy
dmcopy "acis_evt1a.fits[EVENTS][grade=0,2,3,4,6,status=0]" acis_flt1_evt1a.fits opt=all clobber=yes

punlearn dmcopy
dmcopy "acis_flt1_evt1a.fits[EVENTS][@{flt1}][cols -phas]"
acis_evt2.fits opt=all clobber=yes

where {flt1} is ../secondary/acisf02090_001N003_flt1.fits in our example.

punlearn destreak
destreak infile=acis_evt2.fits outfile=cleaned.evt2 ccd_id=8 clobber=yes

et voilá we have an EVT2 file called cleaned.evt2. In many of the later stages, we'll be creating subsets of the data contained within this file.

Ian M George: Make an EVT2 from an EVT1

Set the Bad-pixel maps

Basic Ciao Processing

Onto the Next Step