How to apply the astrometric corrections to the event files ?

calculates equatorial sky coordinates for each photon event in a given photon event list.
Post Reply
tatsumikon
Posts: 13
Joined: Wed Feb 09, 2022 6:52 am

How to apply the astrometric corrections to the event files ?

Post by tatsumikon »

In Section 3.2 of Brunner et al. (2022), the mean linear offsets of the X-ray sources from the Gaia DR2 positions of the reference catalogue are shown in Table 1.

It said that the recalculation of the event coordinates can be done with the tasks "evatt" and "radec2xy".
But I donot understand how to apply the correction to the event files using the above tasks, because "evatt" need the input files with corrected attitude table.
TengLiu
Posts: 3
Joined: Tue Feb 22, 2022 1:05 pm

Re: How to apply the astrometric corrections to the event files ?

Post by TengLiu »

Yes, the first step is to correct the attitude table. Please see an example below:
ftcopy "300007/fm00_300007_120_EventList_001_c001.fits[2][col TIME,RA=RA-0.001389444,DEC=DEC+0.001171111,ROLL]" cor_fm00_300007_120_EventList_001_c001.fits clobber=yes
evatt EVENTFILE=cor_fm00_300007_120_EventList_001_c001.fits ATTFILE=cor_fm00_300007_120_EventList_001_c001.fits GTIFILE=cor_fm00_300007_120_EventList_001_c001.fits
radec2xy file=cor_fm00_300007_120_EventList_001_c001.fits ra0=129.55 dec0=1.5

Any task can be used to modify the ATT file. Here I am using the FTOOLS task ftcopy. Then running evatt and radec2xy can rebuild the event coordinates based on the corrected ATT table.
tatsumikon
Posts: 13
Joined: Wed Feb 09, 2022 6:52 am

Re: How to apply the astrometric corrections to the event files ?

Post by tatsumikon »

declare -a reg=("300007" "300008" "300009" "300010")
declare -a ra_corr=("-0.0013889" "-0.0013889" "-0.0014167" "-0.0013889")
declare -a de_corr=("0.0011667" "0.0008889" "0.0009722" "0.0011111")

for i in "${!reg[@]}"
do
ftcopy "fm00_${reg[$i]}_020_EventList_c001.fits[1][col *,RA=RA+${ra_corr[$i]},DEC=DEC+${de_corr[$i]}]" cor_fm00_${reg[$i]}_020_EventList_c001.fits clobber=yes
evatt EVENTFILE=cor_fm00_${reg[$i]}_020_EventList_c001.fits ATTFILE=cor_fm00_${reg[$i]}_020_EventList_c001.fits GTIFILE=cor_fm00_${reg[$i]}_020_EventList_c001.fits
done
================

The ftcopy task is OK, but the evatt task is failed:

/home/idies/sw/eSASS4EDR/bin/eva: eSASS4EDR Mar 13 05:00:20 2022
/home/idies/sw/eSASS4EDR/bin/eva/GET_INSTNUM: **ERROR3** Error: d not a valid telescope number
/home/idies/sw/eSASS4EDR/bin/eva/EVATT_IN: **ERROR3** Error reading EVENTS file
/home/idies/sw/eSASS4EDR/bin/eva/EVATT_IN: **WARNING1** General warning: With pipeline version 946 the handling of camera time jumps has
/home/idies/sw/eSASS4EDR/bin/eva/EVATT_IN: **WARNING1** In v946 time jumps are now corrected by evprep and new (v15+) INSTPAR files have
/home/idies/sw/eSASS4EDR/bin/eva/EVATT_IN: **WARNING1** According to the header the TIME values in this file are NOT corrected.
/home/idies/sw/eSASS4EDR/bin/eva/EVATT_IN: **WARNING1** The following INSTPAR file is used:
/home/idies/sw/eSASS4EDR/bin/eva/EVATT_IN: **WARNING1**
/home/idies/sw/eSASS4EDR/bin/eva: FAILED

I guess only the file containing the photons detected by a specific telescope module can be processed, how to do with it?
BinnurSokmen
Posts: 4
Joined: Sun May 07, 2023 11:02 am

Re: How to apply the astrometric corrections to the event files ?

Post by BinnurSokmen »

Hello, Did you fix the problem? I have the same issues with the evatt task. Ftcopy is OK for all telescope modules but evatt task is running only first telescope module (TM1) without any warning. I don't understand why.
Post Reply