*** Start of specres.specholucy demo *** [artdata, stsdas.graphics, stsdas.ttools and stecf.specres need to be loaded] Tidying up files from any previous trials dev$pix[1:128,1:128] -> testbas.fits dev$pix[1:256,1:256] -> testbas.fits Creating three identical stellar Gaussian images (FWHM=4pix) psf1.fits -> psf2.fits psf1.fits -> psf3.fits Displaying the stellar point source z1=0.001 z2=0.02 Forming a PSF longslit spectrum from 3 stellar images Slitwidth 0.5arcsec, pixel size 0.05arcsec, centred at Y=128.0 Same pixel sampling as point source Working on PSF # 1 Working on PSF # 2 Working on PSF # 3 rub1_01.fits -> testpsf1.fits Displaying the longslit stellar point source spectrum z1=0. z2=0.1 1st PSF spectrum: Rotation 0deg; flux 5000, Centre Y=120.5 Transforming image rub1_01.fits to image rub1_02.fits xshift: 0.00 yshift: 0.00 xmag: 1.00 ymag: 1.00 xrot: 0.00 yrot: 0.00 2nd PSF spectrum: Rotation 2.5deg; flux 500, Centre Y=199.33 Transforming image rub1_01.fits to image rub1_05.fits xshift: -5.48 yshift: 5.73 xmag: 1.00 ymag: 1.00 xrot: 2.50 yrot: 2.50 3rd PSF spectrum: Rotation -2.5deg; flux 1000, Centre Y=88.67 Transforming image rub1_01.fits to image rub1_08.fits xshift: 5.73 yshift: -5.48 xmag: 1.00 ymag: 1.00 xrot: 357.50 yrot: 357.50 Displaying the 3 longslit stellar point source spectra z1=0. z2=500. Sum the 3 PSF longslit spectra with a spatially varying longslit background spectrum (sky) Transforming image rub1_14.fits to image rub1_15.fits xshift: 257.00 yshift: -0.00 xmag: 1.00 ymag: 1.00 xrot: 270.00 yrot: 270.00 Displaying the background image z1=0. z2=200. Adding Gaussian noise to the model image (RON=6 e; Gain=1) Displaying the noisy simulated longslit point source + background image z1=0. z2=200. Adding wavelength information to header for 2-D spectrum add testspec1.fits,CRPIX1 = 128. testspec1.fits updated add testspec1.fits,CRVAL1 = 8000. testspec1.fits updated add testspec1.fits,CTYPE1 = LAMBDA testspec1.fits updated add testspec1.fits,CD1_1 = 5. testspec1.fits updated Forming the statistical error image 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% - done Forming the data quality image and setting some BAD pixels test256.fits -> testspec1dq.fits Forming the noise image on the background rub1_19.fits -> modelback1.fits specholucy input table for 3 point source spectra: X (pix) Y (pix) Slope (/pix) # Table demopos.tab Tue 14:39:34 25-Jan-2005 # row xchan ychan slope # 1 128.000 120.500 0.00000 2 128.000 199.330 0.04370 3 128.000 88.890 -0.04370 PACKAGE = stecf TASK = specholucy inpima = testspec1 Name of input flux file errima = testspec1er Name of input statistical error file dqima = testspec1dq Name of input data quality file poitab = demopos.tab Name of table of X;Y pixel centres and slope of sp OUTPUT tables and images Rootab = test1t Root name for output point source extracted spectr poi_back= no Output fitted point sources+background (Y) or back bakima = testout1 Name of output fitted background file bakerr = testout1er Name of output fitted background statistical error bakdq = testout1dq Name of output fitted background data quality file POINT SPREAD FUNCTION Spectrum and Position psfima = testpsf1 Name of input Point Spread Function file psfmeth = g Method for fitting PSF peak (C=Centroid; G=Gaussia psfb = 1 Subsampling factor (Y-direction) for PSF image interpol = spline3 Interpolation method used for shifting PSFs POSITIONS of Point Source Spectra posmeth = c Method for fitting input point source peak (E=Exac icsum = 10 Number of columns of input data to sum for cross-c BACKGROUND CHANNEL parameter skernel = 5.0 Sigma of Gaussian for spatial resolution kernel ITERATIONS and Stopping Criteria niter = 50 Number of iterations for point source and 1D background restoration epspoi = 0.0005 Stopping criterion for fractional change in point epsbac1 = 0.0005 Stopping criterion for fractional change in backgr ntrial = 1 No. of Monte Carlo trials for error estimation (seed = 13) Seed for random number generator (dqlim = 10 ) Limiting data quality value above which to ignore ( verbose = yes) Report parameters of restoration during execution !--Spectrum No. 1 at Y pix = 120.49 flux 0.4974E+04 (mean per pix at X= 128) !--Spectrum No. 2 at Y pix = 199.26 flux 0.5660E+03 (mean per pix at X= 128) !--Spectrum No. 3 at Y pix = 88.96 flux 0.9599E+03 (mean per pix at X= 128) !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 1 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 17.98404 !-- Point source no. 2; % flux change 55.14075 !-- Point source no. 3; % flux change 42.55578 !-- Summed background spectrum: % flux change 98.15462 !-- 1D Background NOT converged: % flux change 17.44490 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 6110.0 52.306 6001.1 52.978 6065.8 53.614 2 1120.9 52.306 1142.7 52.978 1097.3 53.614 3 1781.0 52.306 1670.9 52.978 1654.5 53.614 !-- Completed combined iteration no. 1 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 2 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 16.42371 !-- Point source no. 2; % flux change 53.72346 !-- Point source no. 3; % flux change 40.33082 !-- All point sources: Maximum % flux change 53.72346 !-- Summed background spectrum: % flux change 11.11112 !-- 1D Background NOT converged: % flux change 8.43436 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5277.4 74.376 5174.8 74.572 5228.3 75.334 2 739.98 50.547 754.96 49.735 722.83 49.775 3 1284.1 57.045 1203.1 54.817 1191.5 54.644 !-- Completed combined iteration no. 2 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 3 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 2.87397 !-- Point source no. 2; % flux change 15.96982 !-- Point source no. 3; % flux change 10.41675 !-- All point sources: Maximum % flux change 15.96982 !-- Summed background spectrum: % flux change 2.24983 !-- 1D Background NOT converged: % flux change 1.99890 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5137.0 76.111 5035.4 76.347 5088.0 77.098 2 644.73 48.348 659.24 47.543 628.31 47.675 3 1169.1 56.435 1094.0 54.019 1084.6 53.775 !-- Completed combined iteration no. 3 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 4 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.64949 !-- Point source no. 2; % flux change 6.70538 !-- Point source no. 3; % flux change 3.57729 !-- All point sources: Maximum % flux change 6.70538 !-- Summed background spectrum: % flux change 0.68091 !-- 1D Background NOT converged: % flux change 0.99241 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5105.8 77.038 5004.4 77.243 5057.0 77.997 2 608.32 49.035 622.74 48.250 591.66 48.411 3 1132.4 57.038 1057.9 54.789 1048.8 54.638 !-- Completed combined iteration no. 4 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 5 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.22904 !-- Point source no. 2; % flux change 3.58128 !-- Point source no. 3; % flux change 1.64161 !-- All point sources: Maximum % flux change 3.58128 !-- Summed background spectrum: % flux change 0.30142 !-- 1D Background NOT converged: % flux change 0.66557 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5094.9 77.792 4993.5 77.970 5046.0 78.729 2 589.88 50.136 604.03 49.387 572.78 49.542 3 1116.4 57.950 1041.2 55.860 1032.1 55.748 !-- Completed combined iteration no. 5 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 6 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.13168 !-- Point source no. 2; % flux change 2.26050 !-- Point source no. 3; % flux change 0.98726 !-- All point sources: Maximum % flux change 2.26050 !-- Summed background spectrum: % flux change 0.18031 !-- 1D Background NOT converged: % flux change 0.50407 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5088.6 78.430 4987.2 78.592 5039.7 79.356 2 578.55 51.136 592.27 50.432 561.00 50.565 3 1106.8 58.829 1030.9 56.862 1021.9 56.743 !-- Completed combined iteration no. 6 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 7 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.09913 !-- Point source no. 2; % flux change 1.58570 !-- Point source no. 3; % flux change 0.70776 !-- All point sources: Maximum % flux change 1.58570 !-- Summed background spectrum: % flux change 0.12823 !-- 1D Background NOT converged: % flux change 0.40729 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5083.8 78.978 4982.3 79.135 5034.8 79.902 2 570.68 51.964 583.92 51.305 552.73 51.410 3 1099.9 59.601 1023.3 57.727 1014.6 57.573 !-- Completed combined iteration no. 7 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 8 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.08233 !-- Point source no. 2; % flux change 1.18565 !-- Point source no. 3; % flux change 0.55466 !-- All point sources: Maximum % flux change 1.18565 !-- Summed background spectrum: % flux change 0.09929 !-- 1D Background NOT converged: % flux change 0.34224 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5079.9 79.453 4978.3 79.611 5030.8 80.380 2 564.83 52.637 577.57 52.022 546.52 52.097 3 1094.3 60.263 1017.3 58.463 1008.8 58.262 !-- Completed combined iteration no. 8 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 9 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.07070 !-- Point source no. 2; % flux change 0.92245 !-- Point source no. 3; % flux change 0.45364 !-- All point sources: Maximum % flux change 0.92245 !-- Summed background spectrum: % flux change 0.08024 !-- 1D Background NOT converged: % flux change 0.29436 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5076.4 79.868 4974.9 80.031 5027.3 80.801 2 560.28 53.186 572.56 52.612 541.67 52.658 3 1089.7 60.830 1012.3 59.090 1004.2 58.835 !-- Completed combined iteration no. 9 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 10 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.06152 !-- Point source no. 2; % flux change 0.73700 !-- Point source no. 3; % flux change 0.37942 !-- All point sources: Maximum % flux change 0.73700 !-- Summed background spectrum: % flux change 0.06643 !-- 1D Background NOT converged: % flux change 0.25641 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5073.4 80.231 4971.8 80.403 5024.3 81.172 2 556.65 53.639 568.51 53.102 537.78 53.120 3 1085.8 61.319 1008.1 59.627 1000.4 59.317 !-- Completed combined iteration no. 10 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 11 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.05388 !-- Point source no. 2; % flux change 0.60031 !-- Point source no. 3; % flux change 0.32170 !-- All point sources: Maximum % flux change 0.60031 !-- Summed background spectrum: % flux change 0.05585 !-- 1D Background NOT converged: % flux change 0.22563 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5070.8 80.550 4969.2 80.731 5021.6 81.501 2 553.70 54.016 565.17 53.514 534.61 53.506 3 1082.5 61.741 1004.5 60.090 997.16 59.727 !-- Completed combined iteration no. 11 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 12 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.04739 !-- Point source no. 2; % flux change 0.49635 !-- Point source no. 3; % flux change 0.27540 !-- All point sources: Maximum % flux change 0.49635 !-- Summed background spectrum: % flux change 0.04749 !-- 1D Background NOT converged: % flux change 0.19914 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5068.5 80.831 4966.9 81.023 5019.3 81.793 2 551.26 54.333 562.39 53.862 531.97 53.831 3 1079.6 62.109 1001.4 60.492 994.43 60.077 !-- Completed combined iteration no. 12 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 13 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.04180 !-- Point source no. 2; % flux change 0.41544 !-- Point source no. 3; % flux change 0.23752 !-- All point sources: Maximum % flux change 0.41544 !-- Summed background spectrum: % flux change 0.04074 !-- 1D Background NOT converged: % flux change 0.17664 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5066.5 81.080 4964.8 81.282 5017.2 82.051 2 549.22 54.603 560.04 54.160 529.77 54.108 3 1077.1 62.430 998.75 60.842 992.09 60.379 !-- Completed combined iteration no. 13 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 14 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.03696 !-- Point source no. 2; % flux change 0.35130 !-- Point source no. 3; % flux change 0.20610 !-- All point sources: Maximum % flux change 0.35130 !-- Summed background spectrum: % flux change 0.03519 !-- 1D Background NOT converged: % flux change 0.15754 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5064.7 81.300 4963.0 81.512 5015.4 82.281 2 547.50 54.833 558.04 54.416 527.89 54.345 3 1074.9 62.711 996.43 61.149 990.08 60.640 !-- Completed combined iteration no. 14 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 15 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.03275 !-- Point source no. 2; % flux change 0.29970 !-- Point source no. 3; % flux change 0.17976 !-- All point sources: Maximum % flux change 0.29970 !-- Summed background spectrum: % flux change 0.03059 !-- 1D Background NOT converged: % flux change 0.14084 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5063.1 81.495 4961.4 81.717 5013.8 82.485 2 546.04 55.033 556.32 54.639 526.29 54.550 3 1073.0 62.959 994.40 61.420 988.34 60.868 !-- Completed combined iteration no. 15 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 16 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.02907 !-- Point source no. 2; % flux change 0.25765 !-- Point source no. 3; % flux change 0.15748 !-- All point sources: Maximum % flux change 0.25765 !-- Summed background spectrum: % flux change 0.02671 !-- 1D Background NOT converged: % flux change 0.12634 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5061.7 81.669 4959.9 81.899 5012.4 82.668 2 544.78 55.206 554.83 54.833 524.92 54.728 3 1071.3 63.178 992.61 61.658 986.82 61.068 !-- Completed combined iteration no. 16 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 17 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.02585 !-- Point source no. 2; % flux change 0.22301 !-- Point source no. 3; % flux change 0.13850 !-- All point sources: Maximum % flux change 0.22301 !-- Summed background spectrum: % flux change 0.02344 !-- 1D Background NOT converged: % flux change 0.11363 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5060.5 81.823 4958.7 82.062 5011.1 82.830 2 543.69 55.357 553.54 55.003 523.72 54.884 3 1069.8 63.372 991.03 61.869 985.48 61.244 !-- Completed combined iteration no. 17 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 18 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.02303 !-- Point source no. 2; % flux change 0.19421 !-- Point source no. 3; % flux change 0.12224 !-- All point sources: Maximum % flux change 0.19421 !-- Summed background spectrum: % flux change 0.02064 !-- 1D Background NOT converged: % flux change 0.10260 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5059.3 81.962 4957.5 82.208 5009.9 82.975 2 542.74 55.489 552.41 55.153 522.68 55.021 3 1068.5 63.543 989.64 62.057 984.31 61.399 !-- Completed combined iteration no. 18 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 19 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.02054 !-- Point source no. 2; % flux change 0.17006 !-- Point source no. 3; % flux change 0.10822 !-- All point sources: Maximum % flux change 0.17006 !-- Summed background spectrum: % flux change 0.01824 !-- 1D Background NOT converged: % flux change 0.09324 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5058.3 82.085 4956.5 82.338 5008.9 83.105 2 541.91 55.606 551.41 55.287 521.76 55.142 3 1067.4 63.696 988.40 62.225 983.28 61.536 !-- Completed combined iteration no. 19 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 20 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.01836 !-- Point source no. 2; % flux change 0.14966 !-- Point source no. 3; % flux change 0.09608 !-- All point sources: Maximum % flux change 0.14966 !-- Summed background spectrum: % flux change 0.01618 !-- 1D Background NOT converged: % flux change 0.08527 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5057.4 82.195 4955.6 82.455 5008.0 83.222 2 541.18 55.710 550.53 55.405 520.95 55.250 3 1066.3 63.833 987.30 62.374 982.37 61.658 !-- Completed combined iteration no. 20 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 21 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.01643 !-- Point source no. 2; % flux change 0.13232 !-- Point source no. 3; % flux change 0.08553 !-- All point sources: Maximum % flux change 0.13232 !-- Summed background spectrum: % flux change 0.01439 !-- 1D Background NOT converged: % flux change 0.07831 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5056.6 82.295 4954.8 82.560 5007.2 83.326 2 540.54 55.802 549.75 55.511 520.24 55.345 3 1065.4 63.954 986.31 62.508 981.55 61.767 !-- Completed combined iteration no. 21 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 22 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.01473 !-- Point source no. 2; % flux change 0.11748 !-- Point source no. 3; % flux change 0.07632 !-- All point sources: Maximum % flux change 0.11748 !-- Summed background spectrum: % flux change 0.01283 !-- 1D Background NOT converged: % flux change 0.07232 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5055.9 82.384 4954.1 82.654 5006.5 83.420 2 539.97 55.885 549.05 55.606 519.60 55.431 3 1064.6 64.063 985.43 62.628 980.83 61.864 !-- Completed combined iteration no. 22 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 23 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.01323 !-- Point source no. 2; % flux change 0.10473 !-- Point source no. 3; % flux change 0.06826 !-- All point sources: Maximum % flux change 0.10473 !-- Summed background spectrum: % flux change 0.01147 !-- 1D Background NOT converged: % flux change 0.06698 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5055.3 82.464 4953.4 82.739 5005.8 83.504 2 539.46 55.958 548.42 55.692 519.03 55.509 3 1063.8 64.161 984.64 62.736 980.18 61.952 !-- Completed combined iteration no. 23 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 24 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.01190 !-- Point source no. 2; % flux change 0.09371 !-- Point source no. 3; % flux change 0.06119 !-- All point sources: Maximum % flux change 0.09371 !-- Summed background spectrum: % flux change 0.01029 !-- 1D Background NOT converged: % flux change 0.06231 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5054.7 82.537 4952.9 82.815 5005.3 83.580 2 539.00 56.025 547.86 55.769 518.52 55.578 3 1063.2 64.248 983.93 62.834 979.60 62.031 !-- Completed combined iteration no. 24 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 25 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.01073 !-- Point source no. 2; % flux change 0.08414 !-- Point source no. 3; % flux change 0.05497 !-- All point sources: Maximum % flux change 0.08414 !-- Summed background spectrum: % flux change 0.00925 !-- 1D Background NOT converged: % flux change 0.05820 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5054.2 82.602 4952.3 82.884 5004.7 83.649 2 538.60 56.085 547.36 55.839 518.06 55.641 3 1062.6 64.327 983.30 62.922 979.08 62.102 !-- Completed combined iteration no. 25 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 26 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00969 !-- Point source no. 2; % flux change 0.07580 !-- Point source no. 3; % flux change 0.04949 !-- All point sources: Maximum % flux change 0.07580 !-- Summed background spectrum: % flux change 0.00833 !-- 1D Background NOT converged: % flux change 0.05449 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5053.7 82.662 4951.9 82.947 5004.3 83.712 2 538.23 56.139 546.90 55.902 517.64 55.699 3 1062.0 64.397 982.72 63.002 978.61 62.166 !-- Completed combined iteration no. 26 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 27 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00878 !-- Point source no. 2; % flux change 0.06851 !-- Point source no. 3; % flux change 0.04464 !-- All point sources: Maximum % flux change 0.06851 !-- Summed background spectrum: % flux change 0.00753 !-- 1D Background NOT converged: % flux change 0.05115 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5053.3 82.716 4951.4 83.004 5003.8 83.768 2 537.90 56.188 546.48 55.960 517.27 55.751 3 1061.6 64.461 982.20 63.075 978.18 62.224 !-- Completed combined iteration no. 27 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 28 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00797 !-- Point source no. 2; % flux change 0.06211 !-- Point source no. 3; % flux change 0.04035 !-- All point sources: Maximum % flux change 0.06211 !-- Summed background spectrum: % flux change 0.00682 !-- 1D Background converged: % flux change 0.04812 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5052.9 82.765 4951.0 83.056 5003.4 83.820 2 537.60 56.233 546.10 56.013 516.92 55.799 3 1061.1 64.518 981.72 63.140 977.79 62.278 !-- Completed combined iteration no. 28 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 29 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00725 !-- Point source no. 2; % flux change 0.05647 !-- Point source no. 3; % flux change 0.03654 !-- All point sources: Maximum % flux change 0.05647 !-- Summed background spectrum: % flux change 0.00619 !-- 1D Background converged: % flux change 0.04540 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5052.6 82.810 4950.7 83.103 5003.1 83.867 2 537.33 56.274 545.76 56.062 516.61 55.842 3 1060.7 64.570 981.29 63.201 977.44 62.326 !-- Completed combined iteration no. 29 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 30 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00661 !-- Point source no. 2; % flux change 0.05148 !-- Point source no. 3; % flux change 0.03316 !-- All point sources: Maximum % flux change 0.05148 !-- Summed background spectrum: % flux change 0.00563 !-- 1D Background converged: % flux change 0.04293 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5052.2 82.851 4950.3 83.146 5002.7 83.909 2 537.09 56.311 545.44 56.107 516.33 55.882 3 1060.4 64.617 980.90 63.255 977.12 62.371 !-- Completed combined iteration no. 30 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 31 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00605 !-- Point source no. 2; % flux change 0.04706 !-- Point source no. 3; % flux change 0.03015 !-- All point sources: Maximum % flux change 0.04706 !-- Summed background spectrum: % flux change 0.00514 !-- 1D Background converged: % flux change 0.04069 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5051.9 82.889 4950.0 83.185 5002.4 83.949 2 536.86 56.345 545.15 56.148 516.06 55.919 3 1060.0 64.659 980.54 63.306 976.83 62.412 !-- Completed combined iteration no. 31 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 32 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00555 !-- Point source no. 2; % flux change 0.04313 !-- Point source no. 3; % flux change 0.02747 !-- All point sources: Maximum % flux change 0.04313 !-- Summed background spectrum: % flux change 0.00470 !-- 1D Background converged: % flux change 0.03865 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5051.7 82.924 4949.8 83.222 5002.2 83.985 2 536.65 56.376 544.88 56.185 515.82 55.954 3 1059.8 64.697 980.22 63.352 976.56 62.449 !-- Completed combined iteration no. 32 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 33 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00510 !-- Point source no. 2; % flux change 0.03962 !-- Point source no. 3; % flux change 0.02508 !-- All point sources: Maximum % flux change 0.03962 !-- Summed background spectrum: % flux change 0.00431 !-- 1D Background converged: % flux change 0.03679 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5051.4 82.956 4949.5 83.255 5001.9 84.018 2 536.47 56.405 544.64 56.221 515.60 55.985 3 1059.5 64.732 979.92 63.394 976.31 62.484 !-- Completed combined iteration no. 33 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 34 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00471 !-- Point source no. 2; % flux change 0.03648 !-- Point source no. 3; % flux change 0.02294 !-- All point sources: Maximum % flux change 0.03648 !-- Summed background spectrum: % flux change 0.00396 !-- 1D Background converged: % flux change 0.03507 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5051.2 82.986 4949.3 83.286 5001.7 84.049 2 536.29 56.432 544.41 56.253 515.39 56.014 3 1059.3 64.764 979.64 63.433 976.08 62.516 !-- Completed combined iteration no. 34 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 35 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00435 !-- Point source no. 2; % flux change 0.03366 !-- Point source no. 3; % flux change 0.02103 !-- All point sources: Maximum % flux change 0.03366 !-- Summed background spectrum: % flux change 0.00365 !-- 1D Background converged: % flux change 0.03356 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5051.0 83.014 4949.1 83.315 5001.5 84.078 2 536.13 56.457 544.20 56.283 515.20 56.042 3 1059.0 64.792 979.39 63.469 975.87 62.546 !-- Completed combined iteration no. 35 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 36 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00404 !-- Point source no. 2; % flux change 0.03113 !-- Point source no. 3; % flux change 0.01932 !-- All point sources: Maximum % flux change 0.03113 !-- Summed background spectrum: % flux change 0.00337 !-- 1D Background converged: % flux change 0.03220 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5050.8 83.040 4948.9 83.341 5001.3 84.104 2 535.99 56.480 544.01 56.311 515.03 56.067 3 1058.8 64.819 979.16 63.502 975.68 62.573 !-- Completed combined iteration no. 36 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 37 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00376 !-- Point source no. 2; % flux change 0.02885 !-- Point source no. 3; % flux change 0.01778 !-- All point sources: Maximum % flux change 0.02885 !-- Summed background spectrum: % flux change 0.00312 !-- 1D Background converged: % flux change 0.03099 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5050.6 83.065 4948.7 83.366 5001.1 84.129 2 535.85 56.501 543.83 56.337 514.86 56.091 3 1058.7 64.843 978.94 63.533 975.50 62.599 !-- Completed combined iteration no. 37 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 38 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00351 !-- Point source no. 2; % flux change 0.02678 !-- Point source no. 3; % flux change 0.01639 !-- All point sources: Maximum % flux change 0.02678 !-- Summed background spectrum: % flux change 0.00289 !-- 1D Background converged: % flux change 0.02985 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5050.4 83.088 4948.5 83.389 5000.9 84.152 2 535.72 56.520 543.66 56.362 514.71 56.113 3 1058.5 64.864 978.74 63.561 975.33 62.623 !-- Completed combined iteration no. 38 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 39 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00328 !-- Point source no. 2; % flux change 0.02491 !-- Point source no. 3; % flux change 0.01515 !-- All point sources: Maximum % flux change 0.02491 !-- Summed background spectrum: % flux change 0.00269 !-- 1D Background converged: % flux change 0.02877 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5050.3 83.109 4948.4 83.411 5000.8 84.174 2 535.61 56.539 543.50 56.385 514.57 56.133 3 1058.3 64.884 978.56 63.588 975.17 62.646 !-- Completed combined iteration no. 39 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 40 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00308 !-- Point source no. 2; % flux change 0.02321 !-- Point source no. 3; % flux change 0.01403 !-- All point sources: Maximum % flux change 0.02321 !-- Summed background spectrum: % flux change 0.00251 !-- 1D Background converged: % flux change 0.02776 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5050.1 83.129 4948.2 83.432 5000.6 84.194 2 535.50 56.556 543.35 56.406 514.43 56.153 3 1058.2 64.903 978.39 63.613 975.03 62.667 !-- Completed combined iteration no. 40 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 41 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00290 !-- Point source no. 2; % flux change 0.02166 !-- Point source no. 3; % flux change 0.01302 !-- All point sources: Maximum % flux change 0.02166 !-- Summed background spectrum: % flux change 0.00234 !-- 1D Background converged: % flux change 0.02683 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5050.0 83.149 4948.1 83.451 5000.5 84.214 2 535.40 56.572 543.22 56.426 514.31 56.171 3 1058.1 64.919 978.23 63.636 974.89 62.687 !-- Completed combined iteration no. 41 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 42 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00274 !-- Point source no. 2; % flux change 0.02024 !-- Point source no. 3; % flux change 0.01210 !-- All point sources: Maximum % flux change 0.02024 !-- Summed background spectrum: % flux change 0.00219 !-- 1D Background converged: % flux change 0.02593 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5049.9 83.167 4947.9 83.469 5000.3 84.232 2 535.30 56.587 543.09 56.445 514.19 56.188 3 1057.9 64.935 978.08 63.658 974.76 62.705 !-- Completed combined iteration no. 42 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 43 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00259 !-- Point source no. 2; % flux change 0.01895 !-- Point source no. 3; % flux change 0.01127 !-- All point sources: Maximum % flux change 0.01895 !-- Summed background spectrum: % flux change 0.00206 !-- 1D Background converged: % flux change 0.02508 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5049.7 83.184 4947.8 83.486 5000.2 84.249 2 535.22 56.601 542.97 56.462 514.08 56.204 3 1057.8 64.949 977.94 63.678 974.64 62.723 !-- Completed combined iteration no. 43 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 44 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00246 !-- Point source no. 2; % flux change 0.01776 !-- Point source no. 3; % flux change 0.01052 !-- All point sources: Maximum % flux change 0.01776 !-- Summed background spectrum: % flux change 0.00194 !-- 1D Background converged: % flux change 0.02428 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5049.6 83.201 4947.7 83.503 5000.1 84.265 2 535.13 56.614 542.85 56.479 513.98 56.219 3 1057.7 64.962 977.81 63.697 974.52 62.740 !-- Completed combined iteration no. 44 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 45 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00234 !-- Point source no. 2; % flux change 0.01667 !-- Point source no. 3; % flux change 0.00984 !-- All point sources: Maximum % flux change 0.01667 !-- Summed background spectrum: % flux change 0.00183 !-- 1D Background converged: % flux change 0.02352 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5049.5 83.217 4947.6 83.519 5000.0 84.281 2 535.06 56.626 542.75 56.495 513.88 56.234 3 1057.6 64.974 977.68 63.716 974.42 62.756 !-- Completed combined iteration no. 45 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 46 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00223 !-- Point source no. 2; % flux change 0.01567 !-- Point source no. 3; % flux change 0.00922 !-- All point sources: Maximum % flux change 0.01567 !-- Summed background spectrum: % flux change 0.00172 !-- 1D Background converged: % flux change 0.02279 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5049.4 83.232 4947.4 83.533 4999.8 84.296 2 534.98 56.638 542.65 56.509 513.79 56.247 3 1057.6 64.985 977.57 63.733 974.31 62.771 !-- Completed combined iteration no. 46 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 47 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00214 !-- Point source no. 2; % flux change 0.01474 !-- Point source no. 3; % flux change 0.00866 !-- All point sources: Maximum % flux change 0.01474 !-- Summed background spectrum: % flux change 0.00163 !-- 1D Background converged: % flux change 0.02208 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5049.3 83.246 4947.3 83.548 4999.7 84.310 2 534.92 56.649 542.55 56.523 513.70 56.260 3 1057.5 64.996 977.46 63.749 974.22 62.785 !-- Completed combined iteration no. 47 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 48 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00205 !-- Point source no. 2; % flux change 0.01389 !-- Point source no. 3; % flux change 0.00815 !-- All point sources: Maximum % flux change 0.01389 !-- Summed background spectrum: % flux change 0.00155 !-- 1D Background converged: % flux change 0.02140 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5049.2 83.260 4947.2 83.562 4999.6 84.324 2 534.85 56.659 542.46 56.536 513.62 56.272 3 1057.4 65.005 977.36 63.764 974.13 62.799 !-- Completed combined iteration no. 48 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 49 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00197 !-- Point source no. 2; % flux change 0.01309 !-- Point source no. 3; % flux change 0.00768 !-- All point sources: Maximum % flux change 0.01309 !-- Summed background spectrum: % flux change 0.00147 !-- 1D Background converged: % flux change 0.02076 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5049.1 83.274 4947.1 83.575 4999.5 84.337 2 534.79 56.669 542.38 56.549 513.54 56.284 3 1057.3 65.014 977.26 63.778 974.04 62.812 !-- Completed combined iteration no. 49 !--------------------------------------------------------------------------------- !-- Starting combined iteration no. 50 !-- Point sources and background spectra NOT converged !-- Point source no. 1; % flux change 0.00189 !-- Point source no. 2; % flux change 0.01236 !-- Point source no. 3; % flux change 0.00725 !-- All point sources: Maximum % flux change 0.01236 !-- Summed background spectrum: % flux change 0.00140 !-- 1D Background converged: % flux change 0.02014 Star # Flux (# 64) / Back Flux (# 128) / Back Flux (# 192) / Back 1 5049.0 83.287 4947.0 83.587 4999.5 84.350 2 534.74 56.678 542.30 56.561 513.47 56.294 3 1057.3 65.023 977.17 63.792 973.95 62.824 !-- Completed combined iteration no. 50 !-- Combined restoration failed to converge after 50 iterations Point source 1 at 120.5 with 5000 counts # test1t_1.tab FLUX # nrows mean stddev median min max 256 5003.110458 95.8419 4993.81 4753.63 5678.22 Point source 2 at 199.33 with 500 counts # test1t_2.tab FLUX # nrows mean stddev median min max 256 493.6896006 38.7356 494.95 287.603 592.184 Point source 3 at 88.67 with 1000 counts # test1t_3.tab FLUX # nrows mean stddev median min max 256 992.7461538 72.3027 996.233 0. 1101.93 Displaying the fitted background image (range 0-160 counts) z1=0. z2=180. Noise on the fitted background # IMAGE NPIX MEAN STDDEV MIN MAX testout1nois.fits 65536 0.05031 4.417 -85.86 41.03 Compared to that on the input background image # IMAGE NPIX MEAN STDDEV MIN MAX testnois1.fits 65536 0.05682 9.995 -49.07 48.58 *** End of stecf.specres.specholucy demo ***