The problem - history is being inherited from the reference image in
wcsmosaic.
The SCUBA-2 pipeline uses wcsmosaic to mosaic incoming data from the
different arrays into images which are then sky subtracted and
extinction corrected. The SMURF extinction task has a check so that it
(sensibly) won't correct for extinction twice. The problem is that the
data seem to be inheriting history from the reference image specified
in wcsmosaic, which includes a previous call to SMURF extinction, and
so SMURF thinks that the extinction correction has already been made
and does nothing. Uh oh.
For the first set of data, the pipeline manually creates its own
reference image for wcsmosaic to use. It then goes on to mosaic all
the available data and created a single `Group' image. Clearly once a
Group image has been created in the right coordinate system it makes
sense to use that as the reference image, and the pipeline does this
for subsequent data.
I can work around it - by telling the pipeline to always use the
reference image it created as this has no history (other than hisset)
- but it clearly seems wrong that the output mosaic should inherit the
history from the reference.
So the question is why is the history in the *reference* image being
inherited by the other files in the mosaic when it is not part of the
mosaic?
Problem can be demonstrated with data here
http://www.astro.ubc.ca/~agg/wcsmoshis.tgz
and the following commands:
% wcsmosaic in=^in1.lis out=out1 ref=refimage lbnd=! ubnd=!
% cmult out1 1 out1a
% wcsmosaic in=^in2.lis out=out2 ref=out1a lbnd=! ubnd=!
Compare
% hislist out1a
with
% hislist out2
I've also appended an example history listing for an image created in the
second pass. The history is largely the same as that in the Group
image.
Am I missing something?
Andy
History listing:
History listing for NDF structure
/aspersa/scuba2/testdata/reduced/scuba2_short/20060926/s20060926_00001_0002_450
_arr.I1:
History structure created 2008 Jun 20 15:35:31.097
1: 2008 Jun 20 15:35:31.098 - HISSET (KAPPA 1.9-8)
Parameters: MODE='normal' NDF=@gs20060926_00001_450_ref
Software: /aspersa/scuba2/star/bin/kappa/ndfpack_mon
2: 2008 Jun 20 15:35:31.099 - HISSET (KAPPA 1.9-8)
Parameters: MODE='normal' NDF=@gs20060926_00001_450_ref
Software: /aspersa/scuba2/star/bin/kappa/ndfpack_mon
3: 2008 Jun 20 15:35:31.777 - WCSMOSAIC (KAPPA 1.9-8)
Parameters: ACC=0.05 FLBND=[-0.0011765670712247,-2.27587595430452D-5]
FUBND=[2.24774403547338D-5,0.002194763234023] GENVAR=FALSE ILEVEL=2
IN=@^/aspersa/scuba2/testdata/reduced/scuba2_short/20060926/oractempu
wI LBND=! LBOUND=[-33,-123] MAXPIX=1000 METHOD='BILINEAR'
[log in to unmask] REF=@gs20060926_00001_450_ref
UBND=! UBOUND=[46,24] WLIM=0.2
Software: /aspersa/scuba2/star/bin/kappa/kappa_mon
4: 2008 Jun 20 15:35:31.788 - WCSMOSAIC (KAPPA 1.9-8)
Parameters: ACC=0.05 FLBND=[-0.0011765670712247,-2.27587595430452D-5]
FUBND=[2.24774403547338D-5,0.002194763234023] GENVAR=FALSE ILEVEL=2
IN=@^/aspersa/scuba2/testdata/reduced/scuba2_short/20060926/oractempu
wI LBND=! LBOUND=[-33,-123] MAXPIX=1000 METHOD='BILINEAR'
[log in to unmask] REF=@gs20060926_00001_450_ref
UBND=! UBOUND=[46,24] WLIM=0.2
Software: /aspersa/scuba2/star/bin/kappa/kappa_mon
5: 2008 Jun 20 15:35:32.660 - CSUB (KAPPA 1.9-8)
Parameters: [log in to unmask]
[log in to unmask] SCALAR=145.08190394401 TITLE=!
Software: /aspersa/scuba2/star/bin/kappa/kappa_mon
6: 2008 Jun 20 15:35:32.674 - CSUB (KAPPA 1.9-8)
Parameters: [log in to unmask]
[log in to unmask] SCALAR=145.08190394401 TITLE=!
Software: /aspersa/scuba2/star/bin/kappa/kappa_mon
7: 2008 Jun 20 15:35:33.532 - EXTINCTION (SMURF V0.3.1)
Parameters:
IN=@^/aspersa/scuba2/testdata/reduced/scuba2_short/20060926/oractemp8
IZ
OUT=@^/aspersa/scuba2/testdata/reduced/scuba2_short/20060926/oractemp
pa6 METHOD='CSOTAU' QUICK=TRUE HASSKYREM=TRUE CSOTAU=!
MSG_FILTER='NORM'
Software: /aspersa/scuba2/star/bin/smurf/smurf_mon
8: 2008 Jun 20 15:35:33.561 - smf_correct_extinction
Extinction correction successful
Parameters:
IN=@^/aspersa/scuba2/testdata
duced/scuba2_short/20060926/o
temp8PmgIZ
OUT=@^/aspersa/scuba2/testdat
educed/scuba2_short/20060926/
ctempP8Ipa6 METHOD='CSOTAU'
QUICK=TRUE HASSKYREM=TRUE
CSOTAU=! MSG_FILTER='NORM'
Software:
/aspersa/scuba2/star/bin/smurf/s
murf_mon
9: 2008 Jun 20 15:35:33.563 - EXTINCTION
Parameters:
IN=@^/aspersa/scuba2/testdata/reduced/scuba2_short/20060926/oractemp8
IZ
OUT=@^/aspersa/scuba2/testdata/reduced/scuba2_short/20060926/oractemp
pa6 METHOD='CSOTAU' QUICK=TRUE HASSKYREM=TRUE CSOTAU=!
MSG_FILTER='NORM'
Software: /aspersa/scuba2/star/bin/smurf/smurf_mon
10: 2008 Jun 20 15:35:34.063 - MAKEMOS (CCDPACK version 4.0-22)
Parameters: GENVAR=TRUE
IN=@^/aspersa/scuba2/testdata/reduced/scuba2_short/20060926/oractempZ
6U LISTIN=TRUE LOGTO='neither' METHOD='sigma'
[log in to unmask] PRESERVE=TRUE
REF='gs20060926_00001_450_ref' SCALE=FALSE SIGMAS=4 TITLE='URANU'
USEVAR=TRUE WEIGHTS=! ZERO=FALSE
Software: /aspersa/scuba2/star/bin/ccdpack/ccdpack_reg
11: 2008 Jun 20 15:35:34.079 - MAKEMOS (CCDPACK version 4.0-22)
Parameters: GENVAR=TRUE
IN=@^/aspersa/scuba2/testdata/reduced/scuba2_short/20060926/oractempZ
6U LISTIN=TRUE LOGTO='neither' METHOD='sigma'
[log in to unmask] PRESERVE=TRUE
REF='gs20060926_00001_450_ref' SCALE=FALSE SIGMAS=4 TITLE='URANU'
USEVAR=TRUE WEIGHTS=! ZERO=FALSE
Software: /aspersa/scuba2/star/bin/ccdpack/ccdpack_reg
12: 2008 Jun 20 15:35:34.128 - CMULT (KAPPA 1.9-8)
Parameters: [log in to unmask]
OUT=@s20060926_00001_0001_450_cal SCALAR=130000 TITLE=!
Software: /aspersa/scuba2/star/bin/kappa/kappa_mon
13: 2008 Jun 20 15:35:34.139 - CMULT (KAPPA 1.9-8)
Parameters: [log in to unmask]
OUT=@s20060926_00001_0001_450_cal SCALAR=130000 TITLE=!
Software: /aspersa/scuba2/star/bin/kappa/kappa_mon
14: 2008 Jun 20 15:35:34.157 - SETUNITS (KAPPA 1.9-8)
Parameters: MODIFY=FALSE NDF=@s20060926_00001_0001_450_cal
UNITS='mJy/beam/pW'
Software: /aspersa/scuba2/star/bin/kappa/ndfpack_mon
15: 2008 Jun 20 15:35:34.160 - SETUNITS (KAPPA 1.9-8)
Parameters: MODIFY=FALSE NDF=@s20060926_00001_0001_450_cal
UNITS='mJy/beam/pW'
Software: /aspersa/scuba2/star/bin/kappa/ndfpack_mon
16: 2008 Jun 20 15:35:34.260 - NDFCOPY (KAPPA 1.9-8)
Parameters: IN=@s20060926_00001_0001_450_cal LIKE=!
OUT=@^/aspersa/scuba2/testdata/reduced/scuba2_short/20060926/adam
TITLE=! TRIM=FALSE
Software: /aspersa/scuba2/star/bin/kappa/ndfpack_mon
17: 2008 Jun 20 15:35:34.284 - NDFCOPY (KAPPA 1.9-8)
Parameters: IN=@s20060926_00001_0001_450_cal LIKE=!
OUT=@^/aspersa/scuba2/testdata/reduced/scuba2_short/20060926/adam
TITLE=! TRIM=FALSE
Software: /aspersa/scuba2/star/bin/kappa/ndfpack_mon
18: 2008 Jun 20 15:35:34.858 - WCSATTRIB (KAPPA 1.9-8)
Parameters: MODE='mset' NDF=@gs20060926_00001_450_mos REMAP=TRUE
SETTING='system=icrs,skyrefis=ignored,alignoffset=0'
Software: /aspersa/scuba2/star/bin/kappa/ndfpack_mon
19: 2008 Jun 20 15:35:34.871 - WCSATTRIB (KAPPA 1.9-8)
Parameters: MODE='mset' NDF=@gs20060926_00001_450_mos REMAP=TRUE
SETTING='system=icrs,skyrefis=ignored,alignoffset=0'
Software: /aspersa/scuba2/star/bin/kappa/ndfpack_mon
20: 2008 Jun 20 15:35:34.899 - WCSATTRIB (KAPPA 1.9-8)
Parameters: MODE='mset' NDF=@gs20060926_00001_450_mos REMAP=TRUE
SETTING='system=GAPPT,skyrefis=origin,alignoffset=1'
Software: /aspersa/scuba2/star/bin/kappa/ndfpack_mon
21: 2008 Jun 20 15:35:34.909 - WCSATTRIB (KAPPA 1.9-8)
Parameters: MODE='mset' NDF=@gs20060926_00001_450_mos REMAP=TRUE
SETTING='system=GAPPT,skyrefis=origin,alignoffset=1'
Software: /aspersa/scuba2/star/bin/kappa/ndfpack_mon
22: 2008 Jun 20 15:35:40.354 - WCSMOSAIC (KAPPA 1.9-8)
Parameters: ACC=0.05 FLBND=[-0.00066675024979101,-0.0018803140623962]
FUBND=[0.00053229506947252,0.00033720903402707] GENVAR=FALSE ILEVEL=2
IN=@^/aspersa/scuba2/testdata/reduced/scuba2_short/20060926/oractempW
kx LBND=! LBOUND=[-34,-112] MAXPIX=1000 METHOD='BILINEAR'
[log in to unmask] REF=@gs20060926_00001_450_mos
UBND=! UBOUND=[45,35] WLIM=0.2
Software: /aspersa/scuba2/star/bin/kappa/kappa_mon
23: 2008 Jun 20 15:35:40.368 - WCSMOSAIC (KAPPA 1.9-8)
Parameters: ACC=0.05 FLBND=[-0.00066675024979101,-0.0018803140623962]
FUBND=[0.00053229506947252,0.00033720903402707] GENVAR=FALSE ILEVEL=2
IN=@^/aspersa/scuba2/testdata/reduced/scuba2_short/20060926/oractempW
kx LBND=! LBOUND=[-34,-112] MAXPIX=1000 METHOD='BILINEAR'
[log in to unmask] REF=@gs20060926_00001_450_mos
UBND=! UBOUND=[45,35] WLIM=0.2
Software: /aspersa/scuba2/star/bin/kappa/kappa_mon
|