!____________________________________________________________________________________________________________________________________________________________________ ! ! GETSF • Multi-Scale Multi-Wavelength Source & Filament Extraction • Version 210920 • Alexander Men’shchikov, DAp IRFU CEA Saclay !____________________________________________________________________________________________________________________________________________________________________ ! ! Date: 20/09/2021 Time: 16:24:37 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 'l1689b.sw.det.*.an.m.fits' (pixel: 3.00000 arcsec, size: 580 x 580 pixels) ! combined of 'l1689b.163.det.*.an.m.fits' (comboexp: 0.0, combowei: 1.0) ! measured in 'l1689b.163.obs.fits' (units: H2/cm^2) (crpix1: 2.90500E+02, crpix2: 2.90500E+02, crval1: 2.4870175E+02, crval2:-2.4635108E+01) !____________________________________________________________________________________________________________________________________________________________________ ! ! CONFIGURATION PARAMETERS USED BY GETSF: ! ! l1689b 140 | prefix distance .................... ! prefix for naming files (never use ~=.?), distance (pc) to objects ! 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 ! 163 13.5 90 180 n | .................................... ! wave, beam, [-]srcmaxsize, [-]filmaxsize("FWHM), apcorr (Herschel) ! ~/Astronomy/+GETSF_VALIDATION/25_L1689B/+Images/210109/+results ../163 | ... ! paths to image and monochromatic directory ! hi.surface.density.r13p5 | ............................... ! prepared image name ! y y | sources filaments .................. ! extract sources and/or filaments {y|n} ! n 1 1 50 | separate its1 itf1 niter ........... ! separate components and backgrounds; restart iterations, niter ! n 1 1 | flatten itsd1 itfd1 ................ ! flatten background-subtracted components; restart iterations ! n | complete ........................... ! complete separation and flattening ! n | combine ............................ ! combine decomposed structures over wavelengths ! n 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 17 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 ........................................... (H2/cm^2/beam) conversion from H2/cm^2: 2.0651E+02 ! FXP_ERR ....... Peak intensity error ..................................... (H2/cm^2/beam) conversion from H2/cm^2: 2.0651E+02 ! FXT_BST ....... Total flux ............................................... (H2/cm^2*arcsec^2) ! FXT_ERR ....... Total flux error ......................................... (H2/cm^2*arcsec^2) ! FXT_ALT ....... Total flux (Gaussian) or MASS ............................ (Msun) conversion of FXT_BST to MASS (Msun): 1.0329E-26 for a distance of 1.40E+02 pc ! SCALE ......... Characteristic size (FWHM) or ALPHABE .................... ALPHABE = 2.4 * FWHMSIZE * SOUND^2 / GRAVITY / MASS ! 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 - 1.539E-05 (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 = 90.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 | ! |_______________________________________________________________________________| ! !____________________________________________________________________________________________________________________________________________________________________ ! ! 163 µ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 57 22 528.8 209.8 248.4831348 -24.7021722 0 8.900E+01 5.249E+02 0 8.900E+01 5.249E+02 1.300E+24 2.294E+22 2.332E+24 3.346E+22 1.712E-02 1.349E+01 1.574E+01 1.477E+01 2.305E+01 2.017E+01 1.714E+02 -3.528E+00 5.959E+01 5.850E+01 2 5 57 40 287.2 291.6 248.7047753 -24.6341916 0 2.986E+02 3.138E+02 0 2.986E+02 3.138E+02 3.725E+24 7.261E+23 5.002E+25 1.671E+24 5.397E-01 1.525E+00 5.898E+01 5.002E+01 4.671E+01 4.063E+01 7.968E+01 -2.311E+00 9.344E+01 9.302E+01 3 5 57 20 126.3 372.4 248.8522154 -24.5667598 0 4.233E+01 8.633E+01 0 4.233E+01 8.633E+01 5.879E+23 2.416E+22 7.051E+23 2.090E+22 7.152E-03 3.144E+01 1.759E+01 1.252E+01 1.642E+01 1.209E+01 1.584E+02 -2.311E+00 3.587E+01 3.419E+01 4 30 57 48 485.4 470.0 248.5232747 -24.4854117 0 1.876E+01 1.745E+01 0 1.876E+01 1.745E+01 2.148E+23 6.051E+22 6.341E+24 2.060E+23 7.025E-02 1.771E+01 8.705E+01 7.747E+01 6.691E+01 5.867E+01 5.453E+01 -2.311E+00 1.380E+02 1.377E+02 5 26 57 42 427.8 310.1 248.5758855 -24.6187489 0 1.327E+01 1.465E+01 0 1.327E+01 1.465E+01 2.113E+23 5.088E+22 5.855E+24 1.671E+23 6.496E-02 1.858E+01 8.330E+01 7.623E+01 6.736E+01 6.290E+01 7.855E+01 -2.989E+00 1.331E+02 1.329E+02 6 14 57 22 168.3 128.9 248.8138704 -24.7697423 0 9.895E+00 1.114E+01 0 9.895E+00 1.114E+01 1.204E+23 1.232E+22 1.901E+23 1.177E+22 1.913E-03 1.345E+02 1.795E+01 1.608E+01 1.707E+01 1.579E+01 5.004E+01 -2.311E+00 3.881E+01 3.859E+01 7 14 57 19 227.1 70.1 248.7599363 -24.8187423 0 9.161E+00 7.666E+00 0 9.161E+00 7.666E+00 1.232E+23 1.526E+22 1.513E+23 1.259E+22 1.576E-03 1.474E+02 1.663E+01 1.413E+01 1.513E+01 1.317E+01 1.237E+01 -2.311E+00 3.352E+01 3.329E+01 8 15 57 20 168.2 412.8 248.8137346 -24.5331921 0 7.998E+00 7.294E+00 0 7.998E+00 7.294E+00 1.138E+23 1.155E+22 1.838E+23 1.000E+22 1.910E-03 1.429E+02 2.188E+01 1.483E+01 1.785E+01 1.244E+01 9.814E+00 -2.311E+00 3.547E+01 3.466E+01 9 34 57 46 231.2 446.4 248.7560638 -24.5052081 0 1.066E+01 6.891E+00 0 1.066E+01 6.891E+00 8.194E+22 2.844E+22 2.053E+24 8.775E+22 2.221E-02 5.131E+01 8.477E+01 6.674E+01 6.528E+01 4.981E+01 3.445E+01 -2.311E+00 1.252E+02 1.248E+02 10 16 57 33 558.2 386.4 248.4564869 -24.5550350 0 1.005E+01 5.290E+00 0 1.005E+01 5.290E+00 2.623E+23 9.592E+22 1.790E+24 1.569E+23 1.851E-02 3.070E+01 3.863E+01 3.641E+01 3.235E+01 3.067E+01 6.626E+01 -2.311E+00 6.626E+01 6.620E+01 11 17 57 18 221.5 56.5 248.7651065 -24.8300942 0 5.652E+00 4.456E+00 0 5.652E+00 4.456E+00 8.050E+22 1.186E+22 1.115E+23 1.050E+22 1.002E-03 2.263E+02 1.553E+01 1.442E+01 1.663E+01 1.555E+01 1.483E+01 -2.604E+00 3.590E+01 3.581E+01 12 13 57 33 309.5 317.0 248.6843343 -24.6130240 0 1.852E+01 4.316E+00 11 1.852E+01 4.316E+00 5.094E+23 4.097E+23 4.321E+24 6.698E+23 4.563E-02 1.504E+01 4.992E+01 4.108E+01 3.473E+01 2.842E+01 8.882E-01 -2.311E+00 6.627E+01 6.618E+01 13 17 22 22 77.0 347.0 248.8974075 -24.5878986 0 5.573E+00 1.600E+00 0 5.573E+00 1.600E+00 1.068E+23 3.903E+22 1.599E+23 3.730E+22 1.506E-03 1.600E+02 1.737E+01 1.456E+01 1.746E+01 1.496E+01 1.421E+02 -2.311E+00 3.876E+01 3.865E+01 14 13 57 20 270.8 300.8 248.7197787 -24.6264962 0 1.170E+01 1.389E+00 30 1.170E+01 1.389E+00 4.824E+23 5.778E+23 9.797E+23 5.009E+23 9.986E-03 3.123E+01 2.136E+01 1.984E+01 1.736E+01 1.588E+01 1.103E+02 -2.311E+00 3.511E+01 3.510E+01 15 12 57 23 309.3 278.2 248.6845061 -24.6453958 1 1.194E+01 1.213E+00 31 1.194E+01 1.213E+00 3.201E+23 3.564E+23 6.698E+23 3.577E+23 9.104E-03 4.048E+01 2.749E+01 2.153E+01 2.253E+01 1.179E+01 1.659E+02 -2.311E+00 4.067E+01 4.062E+01 16 15 57 18 352.5 311.0 248.6449176 -24.6180143 0 7.138E+00 1.174E+00 10 7.138E+00 1.174E+00 2.327E+23 1.828E+23 3.758E+23 1.436E+23 3.965E-03 7.154E+01 1.973E+01 1.777E+01 1.583E+01 1.426E+01 1.604E+02 -2.311E+00 3.183E+01 3.181E+01 17 12 57 16 310.8 285.8 248.6831083 -24.6389961 1 1.186E+01 1.098E+00 30 1.186E+01 1.098E+00 2.693E+23 3.366E+23 2.966E+23 2.593E+23 3.150E-03 7.069E+01 1.494E+01 1.446E+01 1.376E+01 1.281E+01 1.528E+02 -2.311E+00 3.120E+01 3.120E+01