!____________________________________________________________________________________________________________________________________________________________________ ! ! GETSF • Multi-Scale Multi-Wavelength Source & Filament Extraction • Version 210915 • Alexander Men’shchikov, DAp IRFU CEA Saclay !____________________________________________________________________________________________________________________________________________________________________ ! ! Date: 18/09/2021 Time: 20:26:00 Created by FINALCAT compiled on May 17 2021 19:59 ! Alexander Men'shchikov, DAp IRFU CEA Saclay, France (alexander.menshchikov@cea.fr) ! Used CFITSIO library version 3.260 by William D Pence (william.d.pence@nasa.gov) ! and the utility XY2SKY from WCSTOOLS by Jessica Mink (jmink@cfa.harvard.edu) ! ! CATALOG OF EXTRACTED SOURCES ~> ok: resorted and renumbered ! ! detected in 'sasha.sw.det.*.an.m.fits' (pixel: 2.00000 arcsec, size: 661 x 661 pixels) ! combined of 'sasha.100.det.*.an.m.fits' (comboexp: 0.0, combowei: 1.0) ! measured in 'sasha.100.obs.fits' (units: MJy/sr) (crpix1: 3.31000E+02, crpix2: 3.31000E+02, crval1: 5.2777235E-03, crval2:-9.4443472E-03) !____________________________________________________________________________________________________________________________________________________________________ ! ! CONFIGURATION PARAMETERS USED BY GETSF: ! ! sasha 100 | prefix ............................. ! name (never use ~=.?) to be used as a prefix for naming files ! 1 | nwmax ...............{ nwmax < 100 } ! maximum number of wavelengths specified in this configuration file ! 1 1 | nwaves nw[i] ....................... ! number of waveblengths and their indices to process ! 100 6 6 15 n | .................................... ! wave, beam, [-]srcmaxsize, [-]filmaxsize("FWHM), apcorr ! /Users/amenshch/Astronomy/+GETSF_VALIDATION/34_SASHA/+Images ../100 | ... ! paths to image and monochromatic directory ! Sasha1.flip.mul1.m.m.add-55.m.bor.resamp2as | ............ ! prepared image name ! y y | sources filaments .................. ! extract sources and/or filaments {y|n} ! y 0 2 50 | separate its1 itf1 niter ........... ! separate components and backgrounds; restart iterations, niter ! y 1 1 | flatten itsd1 itfd1 ................ ! flatten background-subtracted components; restart iterations ! y | complete ........................... ! complete separation and flattening ! y | combine ............................ ! combine decomposed structures over wavelengths ! y 2 | detect skelsig ..................... ! detect sources and filaments; skeleton detection significance ! y 1 1 50 n | measure goodmin itn itx sproima .... ! measure structures; minimum goodness; iterations numbers; profile sources ! y | visualize .......................... ! visualize sources with ellipses or dots ! 5 | nsigmacombos (!expert) ............. ! number of sigmas for combining sources in single scales ! 2 | nsigmacombof (!expert) ............. ! number of sigmas for combining filaments in single scales ! 1 | nwavesdetect ....................... ! required minimum number of wavelengths for source detection ! 2 n | savespace delete ................... ! space saving level {0|1|2}: clean up more files; auto delete {?|n|y} ! 0 | verbosity .......................... ! verbosity level {0|1|2} for the screen and log file output !___________________________________________________________________________________________________________________________________ ! ! TABULATED QUANTITIES, TOTAL OF 2 SOURCES ! ! NO ............ Source running number .................................... ! S1 ............ Number of the smallest scale where source first appears .. ! S2 ............ Number of the largest scale where source still exists .... ! SX ............ Number of the footprinting scale for current source ...... ! XCO_P ......... X-coordinate relative to the lower-left pixel (1,1) ...... (pixels) ! YCO_P ......... Y-coordinate relative to the lower-left pixel (1,1) ...... (pixels) ! WCS_ACOOR ..... WCS (J2000) alpha-coordinate of current source ........... (degrees) ! WCS_DCOOR ..... WCS (J2000) delta-coordinate of current source ........... (degrees) ! FG ............ Flag describing global properties over all wavelengths ... see definitions below ! SIGN .......... Global detection significance over all wavelengths ....... SIGN = sqrt ( sum ( SIGNM^2 ) ) ! GOOD .......... Global goodness combining significance and S/N ........... GOOD = sqrt ( sum ( GOODM^2 ) ) ! FM ............ Flag describing monochromatic properties ................. see definitions below ! SIGNM ......... Detection significance from monochromatic single scales .. SIGNM = SSPEAK / SSTDEV single-scale analog of S/N ratio ! GOODM ......... Monochromatic goodness combining significance and S/N .... GOODM = (SIGNM / SIGNM_MIN) * sqrt ( SNRP / SNRP_MIN * SNRT / SNRT_MIN ) * (BFWHM / AFWHM) ! FXP_BST ....... Peak intensity ........................................... (Jy/beam) conversion from MJy/sr: 9.5877E-04 ! FXP_ERR ....... Peak intensity error ..................................... (Jy/beam) conversion from MJy/sr: 9.5877E-04 ! FXT_BST ....... Total flux ............................................... (Jy) ! FXT_ERR ....... Total flux error ......................................... (Jy) ! FXT_ALT ....... Total flux (Gaussian) or MASS ............................ (Jy) ! SCALE ......... Characteristic size (FWHM) or ALPHABE .................... (arcsec) ! AFWHM ......... Major size at half-maximum (FWHM) ........................ (arcsec) interpolated from intensities at half-maximum ! BFWHM ......... Minor size at half-maximum (FWHM) ........................ (arcsec) interpolated from intensities at half-maximum ! ASIZE ......... Major size from intensity moments ........................ (arcsec) below half-maximum for non-Gaussian sources ! BSIZE ......... Minor size from intensity moments ........................ (arcsec) below half-maximum for non-Gaussian sources ! THETA ......... Angle of the major axis (east of vertical) ............... (degrees) position angle PA = THETA - 0.000E+00 (east of north) ! FOFA .......... Footprint expansion scaling factor ....................... ! FOOA .......... Full major axis of an elliptical footprint ............... (arcsec) ! FOOB .......... Full minor axis of an elliptical footprint ............... (arcsec) ! ! DEFINITIONS AND COMMENTS: ! ! 1. Catalog entries may have a value of 0.000E+00 for sources outside the observational mask ! 2. SIGNM has a special value of 9.999E-31 if a source is not visible in clean single scales ! 3. SIGNM has a special value of 7.777E-31 if a source has deconvolved peak intensity equal zero ! 4. FX*_BST and FX*_ERR have a special value of 9.999E-31 if their measurements are < 1.0E-30 ! 5. Monochromatic S/N ratios can be computed as SNRP = FXP_BST / FXP_ERR, SNRT = FXT_BST / FXT_ERR ! 6. Meaning of the global flag FG: ! FG = 0 ..... source is isolated: not blended with any other source in any waveband ! FG = 1 ..... source is blended: footprints intersect by more than 20% in at least one waveband ! FG = 2 ..... source is sub-structured: its footprint area contains at least one other source ! FG = 3 ..... source is sub-structuring: causing a larger source to be sub-structured ! 7. Meaning of the monochromatic flag FM: <~ "negative" when deblending iterations did not converge ! FM = 0 ... source has no problems identified in waveband ! FM = 1 ... source sizes are suspicious: sqrt ( AFWHM * BFWHM ) / BEAM < 0.95 or AFWHM / ASIZE > 1.4 or BFWHM / BSIZE > 1.4 ! FM = 2 ... source is larger than the desired maximum size: sqrt ( FOOA * FOOB ) / FOFA > SRCMAXSIZE = 6.0 arcsec ! FM = 4 ... source is too much elongated: AFWHM / BFWHM > 2 or AFWHM / FOOA is too close to 1. ! FM = 10 ... source peak intensity is not measurable: FXP_BST / FXP_ERR < 2 ! FM = 20 ... source total flux is not measurable: FXT_BST / FXT_ERR < 2 ! FM = 40 ... source peak intensity FXP_BST < 0 or FXP_BST and total flux FXT_BST have opposite signs ! FM = 100 ... source is sub-structured: its footprint area contains at least one mesurable source ! _______________________________________________________________________________ ! | | ! | *!*( SOURCES FROM MULTI-WAVELENGTH CATALOGS SHOULD NEVER BE USED BLINDLY )*!* | ! | | ! | Select only sufficiently good sources in each waveband, on the basis of the | ! | quantities provided: SIGNM, GOODM, FM, FXP_BST / FXP_ERR, FXT_BST / FXT_ERR, | ! | AFWHM / BFWHM, FOOA / AFWHM, etc. | ! | | ! | *!*( RECOMMENDED OPTIMAL SELECTION CRITERIA BASED ON BENCHMARK TESTS )*!* | ! | | ! | abs ( GOODM ) > 1 & abs ( SIGNM ) > 1 & FXP_BST / FXP_ERR > 2 & | ! | FXT_BST / FXT_ERR > 2 & AFWHM / BFWHM < 2 & FOOA / AFWHM > 1.15 | ! |_______________________________________________________________________________| ! !____________________________________________________________________________________________________________________________________________________________________ ! ! 100 µm ! .____________________________________________________________________________________________________________________________________________________________________________________. ! ! 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 ! NO S1 S2 SX XCO_P YCO_P WCS_ACOOR WCS_DCOOR FG SIGN GOOD FM01 SIGNM01 GOODM01 FXP_BST01 FXP_ERR01 FXT_BST01 FXT_ERR01 FXT_ALT01 SCALE01 AFWHM01 BFWHM01 ASIZE01 BSIZE01 THETA01 FOFA01 FOOA01 FOOB01 ! 1 2 37 8 385.5 370.0 359.9749999 0.0122223 1 2.065E+01 9.957E+00 0 2.065E+01 9.957E+00 6.635E-02 1.105E-02 6.811E-02 8.554E-03 6.089E-02 6.000E+00 6.932E+00 5.193E+00 5.908E+00 4.836E+00 1.582E+02 -2.311E+00 1.400E+01 1.387E+01 2 3 37 32 389.5 370.5 359.9727777 0.0125001 1 1.518E+01 2.877E+00 21 1.518E+01 2.877E+00 7.078E-02 1.180E-02 4.981E-02 1.428E-02 3.060E-02 9.425E+00 6.522E+00 5.520E+00 2.000E+00 2.000E+00 6.569E+01 -2.311E+00 2.182E+01 2.174E+01