Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.stecf.org/software/IRAFtools/stecf-iraf/specres/demo/extra/rub1.txt Дата изменения: Thu Jun 8 18:01:05 2006 Дата индексирования: Tue Oct 2 08:10:15 2012 Кодировка: Поисковые слова: comet

*** 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 ***