Changeset 1633 for branches/alma/python/asapmath.py

Timestamp:

09/16/09 17:11:16 (15 years ago)

Author:

Takeshi Nakazato

Message:

New Development: No

JIRA Issue: Yes CAS-1422

Ready to Release: Yes

Interface Changes: Yes

What Interface Changed: Defined calibrate() and almacal() in python/asapmath.py

Defined cwcal() in STMath class

Test Programs: List test programs

Put in Release Notes: Yes

Module(s): task_sdaverage

Description: Describe your changes here...

asapmath.py is changed to be able to calibrate data for APEX.
This modification is tentatively considered a calibration of ALMA
single-dish data. However, it must be updated in the future since
I don't know how raw ALMA data are provided and I have to change
code along with read ALMA data.

The calibrate() function takes calibration mode from its argument and
looks antenna name of the input scantable, and calls appropriate calibration
function depending on the calibration mode and antenna name.
If antenna name include 'APEX' or 'ALMA', newly defined calibration function
apexcal() is called. For other antenna name, one of the existing calibration
functions (calps, calnod, calfs, auto_quotient) is called.

File:

• branches/alma/python/asapmath.py (modified) (1 diff)

branches/alma/python/asapmath.py

@@ -850 +850 @@
     return s
 
-## def apexcal( scantab, calmode ):
-##     """
-##     Calibrate APEX data.
+def calibrate( scantab, scannos=[], calmode='none', verify=None ):
+    """
+    Calibrate data.
     
-##     Parameters:
-##         scantab:       scantable
-##         calmode:       calibration mode
-##     """
-##     from asap._asap import stmath
-##     stm = stmath()
-##     if ( calmode == 'ps' ):
-##         asaplog.push( 'APEX position-switch calibration' )
-##         print_log()
-##     elif ( calmode == 'fs' ):
-##         asaplog.push( 'APEX frequency-switch calibration' )
-##         print_log()
-##     elif ( calmode == 'wob' ):
-##         asaplog.push( 'APEX wobbler-switch calibration' )
-##         print_log()
-##     elif ( calmode == 'otf' ):
-##         asaplog.push( 'APEX On-The-Fly calibration' )
-##         print_log()
-##     else:
-##         asaplog.push( '%s: unknown calibration mode' % calmode )
-##         print_log('ERROR')
+    Parameters:
+        scantab:       scantable
+        scannos:       list of scan number
+        calmode:       calibration mode
+        verify:        verify calibration     
+    """
+    print 'param: %s, %s, %s' %(scannos, calmode, verify) 
+    antname = scantab.get_antennaname()
+    if ( calmode == 'nod' ):
+        asaplog.push( 'Calibrating nod data.' )
+        print_log()
+        scal = calnod( scantab, scannos=scannos, verify=verify )
+    elif ( calmode == 'quotient' ):
+        asaplog.push( 'Calibrating using quotient.' )
+        print_log()
+        scal = scantab.auto_quotient( verify=verify )
+    elif ( calmode == 'ps' ):
+        asaplog.push( 'Calibrating %s position-switched data.' % antname )
+        print_log()
+        if ( antname.find( 'APEX' ) != -1 or antname.find( 'ALMA' ) != -1 ):
+            scal = almacal( scantab, scannos, calmode, verify )
+        else:
+            scal = calps( scantab, scannos=scannos, verify=verify )
+    elif ( calmode == 'fs' or calmode == 'fsotf' ):
+        asaplog.push( 'Calibrating %s frequency-switched data.' % antname )
+        print_log()
+        if ( antname.find( 'APEX' ) != -1 or antname.find( 'ALMA' ) != -1 ):
+            scal = almacal( scantab, scannos, calmode, verify )
+        else:
+            scal = calfs( scantab, scannos=scannos, verify=verify )
+    elif ( calmode == 'otf' ):
+        asaplog.push( 'Calibrating %s On-The-Fly data.' % antname )
+        print_log()
+        scal = almacal( scantab, scannos, calmode, verify )
+    else:
+        asaplog.push( 'No calibration.' )
+        scal = scantab.copy()
+
+    return scal 
+
+def almacal( scantab, scannos=[], calmode='none', verify=False ):
+    """
+    Calibrate APEX and ALMA data
+
+    Parameters:
+        scantab:       scantable
+        scannos:       list of scan number
+        calmode:       calibration mode
+        verify:        verify calibration     
+    """
+    from asap._asap import stmath
+    stm = stmath()
+    antname = scantab.get_antennaname()
+    ssub = scantab.get_scan( scannos )
+    scal = scantable( stm.cwcal( ssub, calmode, antname ) )
+    return scal