Changeset 670 for trunk/python

Timestamp:

09/02/05 14:38:51 (19 years ago)

Author:

mar637

Message:

fixed typo in import rotate_xyphase (from release-2), added auto-quotient

File:

• trunk/python/scantable.py (modified) (2 diffs)

trunk/python/scantable.py

@@ -1195 +1195 @@
         if allaxes is None: allaxes = rcParams['scantable.allaxes']
         varlist = vars()
-        from asap._asap import rotate_xyphase as _rotate_xyphase
+        from asap._asap import _rotate_xyphase
         _rotate_xyphase(self, angle, allaxes)
         self._add_history("rotate_xyphase", varlist)
@@ -1255 +1255 @@
             return
 
-
+    def auto_quotient(self, mode='suffix', preserve=True):
+        """
+        This function allows to build quotients automatically.
+        It assumes the observation to have the same numer of
+        "ons" and "offs"
+        It will support "closest off in time" in the future
+        Parameters:
+            mode:           the on/off detection mode; 'suffix' (default)
+                            'suffix' identifies 'off' scans by the
+                            trailing '_R' (Mopra/Parkes) or
+                            '_e'/'_w' (Tid)
+            preserve:       you can preserve (default) the continuum or 
+                            remove it.  The equations used are 
+                            preserve: Output = Toff * (on/off) - Toff
+                            remove:   Output = Tref * (on/off) - Ton
+        """
+        modes = ["suffix","time"]
+        print mode
+        if not mode in modes:
+            print "please provide valid mode. Valid modes are %s" % (modes)
+            return None
+        from asap._asap import quotient as _quot
+        if mode == "suffix":
+            srcs = self.get_scan("*[^_ewR]")
+            refs = self.get_scan("*[_ewR]")
+            return scantable(_quot(srcs,refs, preserve))       
+        else:
+            print "not yet implemented"
+            return None
+        
     def quotient(self, other, isreference=True, preserve=True):
         """