4137
Comment:
|
4067
|
Deletions are marked like this. | Additions are marked like this. |
Line 37: | Line 37: |
|| [[attachment:simput_test.fits|simput_test.fits]] || Source Map || |
Description of performance tests and simulations relevant to the APETOOL and SHAPELIB libraries.
Simulations are underway to test the following aspects relevant to apetool and shapelib:
- How well/accurate can the shapelet-PSF model determine the flux of point-sources: simulate fields (both survey-mode and pointed) with a grid of point sources of fixed flux and spectral. Run ermldetect to detect the sources and measure their flux and compare with the input one.
How well can APETOOL reconstruct the selection function of eROSITA observations: simulate fields (both survey-mode and pointed) with a point-sources with fluxes drawn from a logN-logS distribution and fixed spectral shape. Run ERLDETECT and ptoduce a source list, APETOOL to threshold at a given false detection probability and estimate field sensitivities. Reconstruct input logN-logS and compare.
The relevant simulations are based on modified versions of Nicolas Clerc's python scripts described here. The main modifications are: (i) functions to generate attitude files for both survey and pointed observations (ii) place sources in a grid of fixed spacing within the survey field, (iii) use of simple power-law X-ray spectra with a choice of Gammas for the point sources, (iv) functionality to produce sources of fixed flux. A tarball with the relevant python scripts can be found [[|here(to be added soon)]].
Progress so far
Simulations of sources with flux fX(0.5-2)=2*10-12cgs in grid. Pointed observations and survey fields have been simulated. The latter use the attitude file: oneyearsurvey_att.fits. The SIXTE and Simput versions 2.2.0 have been used. The EROSIM task with the XML files and the Panter PSF (256M!) was run followed by the ERO_CALEVENTS task.
The table below lists links to the input simput file, the event files for each telescope separately and the merged one (via call to evtool):
simput file |
|
Event file for CCD1 |
|
Event file for CCD2 |
|
Event file for CCD3 |
|
Event file for CCD4 |
|
Event file for CCD5 |
|
Event file for CCD6 |
|
Event file for CCD7 |
|
Merged Event file |
Issues with simulations (to be addressed in the future):
- the event files produced by ERO_CALEVENTS have colsixte
umns RA and DEC of type 'J' (integer) in the 'EVENTS' extension. Should be 'E' or 'D'
- the event files produced by EROSIM include events with negative energies.
Results:
The demoscript.py was used to analyse the data. The ERMLDET task was called with the "shapelet_flag=yes". The table below lists the Source Map in the 0.5-2keV band produced by the ERMLDET task, the corresponding SIXTE-simulated image and the source list produced by the the ERMLDET task:
Source Map |
|
Comparison of Source Map with SIXTE simulated image in the 0.5-2keV |
|
Source list produce by ERMLDET |