Sources detection method used at Saclay : MR_filter + SExtractor

Download Report

Transcript Sources detection method used at Saclay : MR_filter + SExtractor

Sources detection method used at Saclay :
MR_filter + SExtractor (current placeholder in catalog pipeline).
Already applied at DC1 (J. Ballet presentation at closeout)
http://www-glast.slac.stanford.edu/software/Workshops/Feb04DC1CloseOut/agenda.htm
Optimal filter method. developed by J. Ballet.
Already applied at DC1 (same presentation as above)
3D wavelet filtering (space + energy or time). developed by J.L. Starck.
Some ideas presented in J.L. Starck's summary talk in May 2005
http://www-glast.slac.stanford.edu/ScienceTools/workshops/May05DiffuseCat/agenda.htm
My modest contribution for DC2:
use a 2d binned Poisson likelihood
Use the Egret fortran program LIKE (1994)
J.R. Mattox, D.J. Macomb, D.L.Bertsch, J.A. Esposito
Why ?
I already used it on Egret data.
Not fancy but it works very well.
Fast.
Adjust the diffuse emission model (Gmult, Gbias).
Characterisation not as good as 3d models.
Alternative method for the catalog.
What have I done so far ?
Developped a cpp wrapper for Like:
Like *like=new Like (“counts_g12_30_g002.fits”,2,“cfgas.gal.g002b”,20);
like->Set_ROI (60,120,-90,90);
like->Set_Sources ("sources_gal_g004b_17.list");
like->Set_PSFADJ ("MA0");
like->Set_PSFADJ ("MAR");
like->do_LM (0.05,25,“y","LM_Output");
delete like;
Run Like with the 1997 (P14) and 2001 (P19-reprocessed) Egret data
Combined it with Sextractor to detect the likelihood map excess
Code to merge results from all maps (gal,eq, energies, zones, periodes) into a catalog
48
47
46
45
44
43
42
41
40
38
37
36
35
34
33
32
31
30
28
27
26
25
24
23
22
21
20
18
17
16
15
14
13
12
11
10
8
7
6
5
4
3
2
1
0
Results with Egret data and diffuse emission :
All the sources in catalog were found.
Found 12 new sources (not in catalog nor in Egret lower sign. source list)
82 (80) catalog sources are not any more significant.
61 (70) new significant sources should appear in the catalog.
To do list for DC2
Adapt Like for resolution other than 0.5 deg.
Generate:
Like-PSF file from DC2 PSF
Counts and Exposure maps (weeks, days, energy bands)
Galactic diffuse emission map
Run it (1308 cpu at In2p3 cc-Lyon )