From jhibbard@cv3.cv.nrao.edu Thu Jan  3 16:21:17 2002
Date: Thu, 03 Jan 2002 18:09:35 -0500
From: John E. Hibbard <jhibbard@cv3.cv.nrao.edu>
To: aips2-naug@cv3.cv.nrao.edu
Subject: Re: test log from last year; update to follow

ooops - wrong txt file was attached. This is the correct one.
-j
    [ Part 2: "Attached Text" ]

#JEH> Continuation of test log from aips++ data reduction
#JEH> of NGC 4038 D-array HI spectral line data.
#JEH> Log 4. 12/20/01
#
#JEH> Well last we left our fearful tester, he was frustrated
#JEH> by the strange amplitude levels of his bandpass, and
#JEH> continuum flux measurements that were a factor of ~6
#JEH> lower than he found in his previous work on this system.
#JEH> 
#JEH> In todays episode, he will map his source both with and 
#JEH> without the bandpass applied, to see if he can isolate 
#JEH> where the scale factor is coming in.
#JEH> We remind the interested reader that fluxscale gave the
#JEH> correct flux for the secondary flux calibrator, 1130-148.
#JEH> Should be 5.335+/-0.002 Jy
#
#JEH> Using aips++ weekly  (Ver 1.7, build 061)

#JEH> Tried starting from existing ms, but map looked like junk.
#JEH> Not sure if running cal.correct twice is legit. So lets start 
#JEH> from scratch. First ran script that did both G and B calibration, 
#JEH> and got same result as last time (fluxes too high by factor 
#JEH> of ~6). Now redo, but only apply gain calibration:

include 'ms.g'
fitsfilename :='/export/data_1/jhibbard/DATA/n4038/AH582/AH582.960603.LINE';
myms:=fitstoms(fitsfile=fitsfilename, msfile='n4038D.ms');
myms.done();
include 'flagger.g'
fg:=flagger(msfile="/export/data_1/jhibbard/aips++/n4038D.ms" );
fg.setchan(1:127);
fg.setflagmode('flag');
fg.filter(column='DATA', operation='range', comparison='Amplitude', range=['0.00001Jy','5Jy']);
#This doesn't work as expected:
fg.query('ANTENNA1==8 and ANTENNA2==23');
#have emails out to the NAUG to see if someone can explain to me
#how to flag a baseline.
#
fg.setantennas([4]);
fg.setids(fieldid=[1]);
fg.state();
fg.timerange('1996/06/03/00:57:30','1996/06/03/01:02:00');
#
fg.setantennas([18]);
fg.setids(fieldid=[3]);
fg.state();
fg.timerange('1996/06/03/00:59:00','1996/06/03/02:56:00');
#
fg.setantennas([10]);
fg.setids(fieldid=[1]);
fg.state();
fg.timerange('1996/06/03/00:56:00','1996/06/03/01:12:00');
#
fg.setantennas([14]);
fg.setids(fieldid=[1,2,3]);
fg.state();
#fg.timerange('1996/06/03/00:56:00','1996/06/03/02:41:30');
fg.timerange('1996/06/03/00:56:00','1996/06/03/04:23:30');
#
fg.done();
#
include 'imager.g';
im:=imager('n4038D.ms');
im.setjy(fieldid=1);
im.done();
include 'calibrater.g';
cal:=calibrater('n4038D.ms');
cal.setdata(msselect='FIELD_ID<=2');
cal.setsolve(type='G',t=300.0,preavg=300.0,refant=24,table='n4038D.gcal');
cal.solve();
cal.plotcal(tablename='n4038D.gcal');
cal.fluxscale(tablein='n4038D.gcal', tableout='n4038D.gcalout',
      reference='1331+305', transfer=['1130-148']);
cal.plotcal(tablename='n4038D.gcalout');
cal.done();

cal:=calibrater('n4038D.ms');
cal.setdata(msselect='FIELD_ID IN [1:3]');
cal.setapply(type='G',t=0.0,table='n4038D.gcalout',select='FIELD_ID==2');
cal.correct();
cal.done();

#JEH> Lets first make a map of "ch 0" (ch 11-107)
#JEH> NOTE!!! THIS IS NOW A MAP *WITHOUT* THE BANDPASS APPLIED.
#JEH> SO DIFFERENT FROM BEFORE.

#JEH> Lets first make a map of "ch 0" (ch 11-107)
include 'imager.g'
imagr:=imager('n4038D.ms');
imagr.setdata(fieldid=3, mode='channel', nchan=96, start=11, step=1,
                spwid=1);
imagr.setimage(nx=256, ny=256, cellx='10arcsec', celly='10arcsec', 
                 stokes='I', fieldid=3, mode='channel', nchan=1, 
                 start=11, step=96, spwid=1);
imagr.weight(robust=-1);
imagr.make('n4038Dch0.model');
imagr.clean(algorithm='clark', niter=3000, gain=0.1, threshold='0.0Jy',
               model='n4038Dch0.model', image='n4038Dch0.restored',
               residual='n4038Dch0.residual');
imagr.restore(model='n4038Dch0.model', 
                 image='n4038Dch0.restored2');
imagr.done();
im:=image('n4038Dch0.restored2');
im.view();   

#JEH> Looks ok, but can't get statistics to work from the viewer!
#JEH> Bugged that as defect AOCso03285
#JEH> Moving cursor around looks like peak fluxes are closer to the
#JEH> truth. Comparing with AIPS and previous results:
#JEH>           aips++, G only  aips++, G&B       AIPS
#JEH> n4038 peak: 104.1 mJy/b    18.8 mJy/b    122.6 mJy/b
#JEH> sw peak:     63.3 mJy/b    11.5 mJy/b     66.3 mJy/b
#JEH> se peak:    111.9 mJy/b    20.0 mJy/b    113.3 mJy/b
#JEH> east peak:   31.9 mJy/b     5.9 mJy/b     32.3 mJy/b
#JEH> So problem is with applying bandpass plus Gain.
#JEH> Great pate mom, gotta motor if I want to make that funeral.
#JEH> More next year.