to refresh my memory of what is
available - and there is exportfit listed. I can also type taskhelp
and I think this list is independent. For example the flagdata task
is only listing when casapy is started. If these task availability
listings do not come from the same list, they should. This may hold
for other information as well. You don't want the same information
listed in duplicate lists. This is identical to my flagger comment
above.
- 7. Image tool and imagefits:
The uvfits was in the image tool in aips++. Unfortunately, this tool
has not been pythonized yet, so there is no way of getting an image
out of casa in casapy. You have to go the aips++, I think. Perhaps,
in the short term, the uvfits reader and writer can be incorporated
into the ia (image analysis utilities tool). I successfully made the
fits file using aips++.
-
8. Catalog;filecatalog is now available, though we encourage the use of the unix shell as this is the simplest most direct route.
I like the aips++ catalog routine and some of the interaction with it.
Can this be implement at some level? Most of my deletions of incorrect
gain tables and measurement sets were done in unix which is okay.
The clean boxes work! Thanks.
Some parameter considerations: I think bmaj, bmin and bpa (clean beam
size) should be included.
It is not at all clear that if you want to clean in the same way as
IMAGR in aips, then you should choose the alg='csclean'. The other
options will only clean accurately the inner quarter area and are more
like APCLN in AIPS which hardly anyone uses. These options need more
documentation in several places.
Additional parameters needed in clean are timerange, antennas
Do you need uvfilter at all? If uvfilterbmaj = 0.0, that is the
number to decide on filtering. Also, I think users would be more
comfortable with uvtaper (the inverse relationship) rather than the
filtering.
rmode is unfamiliar to most of use. It needs a bit more information
in help clean
The need for tying calibration files to the relevant data set
is becoming more clear. For example, fluxscale has the visibility
data set as input because it needs the visibility data set to
translate fieldids into source names. As yet, plotcal does not point
to a measurment set, and this is the reason fieldid and antennaid
cannot be translated.
My general feeling is that these calibration tables should become
tables within the structure of the measurement set. These are already
many table directories under the measurement set, so why not add
calibration tables. In talking with George about various calibration
concerns, he mentioned that Robert Lucas is working on the calibration
tables for the asdm and quick-look, and the casa group is waiting to
see this implementation. I believe there is still the consideration
of how 'independent' the calibration tables should be with respect to
the measurement set. This is an important decision.
- 11. Casa parameter defaults:
I notice that the default parameters for casa parameters are listed in
the casa.py file in the appropriate directory, in this case
/usr/lib/casapy/20.0.993-003/lib/python2.5
I guess we are stuck with global defaults, rather than different defaults
for different tasks. I would also recommend changes some of the defaults
such as:
fieldid = -1 instead of fieldid = 0
datacolumn = 'data' instead of datacolum = 'corrected'
field is in the list twice
multiplot = true instead of multiplot = false
alg = 'csclean' instead of alg = 'hogbom' (ie APCLN)
opacity = false instead of opacity = true
gaincurve = false instead of gaincurve = true
plotsymbol = '.' instead of plotsymbol = ','
timebin = 0.0 instead of timebin = '-1s' (see above comments on
this)
timerange = ['01-JAN-1980/00:00:00.0','01-JAN-2010/00:00:00.0']
instead of timerange = ''
(It is difficult finding the various syntax for timerange. If you
use the
default above, it is obvious what format to use if you really want
to use
timerange.)
verbose = true instead of verbose = false
I suspect plotcal is being worked on, but here are some suggestions
which might
not yet be considered.
a. The plotted points are too large and blend into each other. Need a
parameter
to scale the plotting character size.
b. As discussed above, the input vis is needed to translate id numbers
c. A range parameter for the x- and y- axes are needed
Since fluxscale already points to the measurement set, the
fluxscale pythong script should automatically run setjy on this
source(s). At the present time, the flux that fluxscale determines is
not available, I believe, so that in running setjy on your own, you
have to write in the value(s) specifically. Perhaps a flag in
fluxscale indicating whether to apply or not apply this value to the
measurement set is useful.
I can see the text information for the inputs to any task, and
making changes is simple. But, where is the information used for the
help files? I believe that these should be an easily accessible and
changeable, so that changes to the inp and help functions for any task
can be more easily coordinated.
In the filler, please flag autocorrelations by default.
Add a parameter to the task inputs that you can fill them if you
want. (request by Mark Claussen and it is a good one but should be low priority)
---------------------------------------
setjy ('vis=src.ms',field='1331+305')
# typo created major error messages that are not helpful:
# Exception Reported: Table
#/home/sola/dss/mcmaster/spitzer.vla/vla.reduce/Aarray2.multiobs/
#vis=src.ms/FIELD
#does not exist
#Fri Jan 26 20:22:29 2007 WARN :
#No table specified, please open first
#Exception Reported: Table
#/home/sola/dss/mcmaster/spitzer.vla/vla.reduce/Aarray2.multiobs/vis=src.ms
#does not exist
# etc...
# just say that the vis file "vis=src.ms" doesn't exist and leave
it at that.
- Asthetics change request:
browsetable ' src.ms'
# why does this come up as a small box in a large white screen?
# why not bring it up so the table fills the available space?
-------------------------------
setjy(vis='src.ms',field='1331+305')
# gives an error - details not reported here because I've already
worked them with Joe:
##### Joe figured it out, selection between fieldid and field is not
handed off correctly. eventually get rid of fieldid but not
until casa has wildcard parameters in the field names. Minimum
match has too many problems getting multiple sources.
-------------------------------
plotcal gets TPPlotter error when trying to plot a gain table:
gaincal(vis='src.ms',caltable='src.gcal',mode='none',
msselect='FIELD_ID IN [0] && DATA_DESC_ID IN [0,1]',
uvrange=[0,400],refant=20)
# OK, 2 good solutions found
plotcal(tablein='src.gcal',yaxis='amp',multiplot=True)
# a bunch of code traceback, I get a blank screen.
# Fri Jan 26 22:10:23 2007 SEVERE :
# Exception Reported: TPPlotter: From Python.
- problem is repeatable.
-------------------------------
Note: plot cal is marching through the antennas in numerical
sequential order, This suggests that it is showing antenna_id
rather than antenna_name like plotxy. This should be made to be
consistent.
-
Enhancement request:Behaves as for other plots.
-------------------------------
When I stop iterating in plotcal or plotxy, the image goes away. I
would like it to stay up on the screen so I can see the last plot.
-
Problem with understanding:
-------------------------------
# why do I need a visibility file for fluxscale? Isn't it still just
operating on the cal tables???? This is disturbing.
It is used for referencing information from the original MS. This
is also planned to be used for then setting the value without requiring
an additional step
-------------------------------
fluxscale causes segmentation fault:
fluxscale(vis='src.ms',caltable=' src.gcal',fluxtable='src.fcal',
reference='1331+305', transfer=['1820-254','2007+404','1925+211','2015+371'])
# Fri Jan 26 23:15:23 2007 WARN :
# The following transfer fields have no solutions available:
# Fri Jan 26 23:15:23 2007 NORMAL : Found reference field(s): 1331+305
# Segmentation fault (core dumped)
# src.fcal not created.
# Summary shows that all field names are correct... solutions were
obtained for field-ids 0,1,3,5,7 which are the sources I list
above. So this doesn't appear to be a wrong transfer field.
# restart casa and try again:
fluxscale(vis='src.ms ',caltable='src.gcal',fluxtable='src.fcal',
reference='1331+305',transfer=['1820-254','2007+404','1925+211','2015+371'])
# Fri Jan 26 23:23:42 2007 SEVERE Calibrater::fluxscale:
# Caught Exception: Cannot find solutions for transfer field(s)
# but this time it did not cause a seg fault.
# maybe since the transfer fields were appended there is a problem?
- Review last week's work
- Dailys in CV:
- Ed:
- - clarification on antenna based flagging
- - log mouse copying
- - exportuvfits improvements
- X refant confusion (need to confirm)
- - fluxscale
- David:
- X
bandpass documentation
- X
autoflag segfault
- X
autoflag parameter type issue (fixed in tomorrow's daily)
- X
gaincal defaults Changed defaults to include calwt True for opacity,gaincurve
- - importvla bandname issue
- Debra:
- X
msplot problemsSee above
- - mplot data selection issues
- X
regression failure (used old script)
- X
plotxy problems (.matplotlib problem)
- X
setjy confused?fieldid/field selection hierarchy confusing. Consolidated to single parameter='field' which will take either an integer or list of integers or a string.
-
ipython/python versioning - use 'special' ipythonrc_casa
- Writes to local directory .ipython_casa directory in local directory. Should this be the same place as the others, i.e., $HOME/...?
-
EF: seg faults - data corruption or other? We think this was a combination of effects from possible actual data corruption between conflicting versions (running the pipeline version), rpm built without proper flags, and an actual bug in fluxscale.
Alpha Release notes
Enhancements/Defects resolved
New functionality
-- JosephMcMullin - 26 Jan 2007
* flag_casapy.html: flag_casapy.html
* flag_casapy.html: flag_casapy.html
-- JosephMcMullin - 28 Jan 2007