Changeset 1145 for trunk/python


Ignore:
Timestamp:
08/16/06 11:11:45 (18 years ago)
Author:
mar637
Message:

added average_beam which is part of Ticket #45

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/python/scantable.py

    r1143 r1145  
    10211021        return s
    10221022
     1023    def average_beam(self, mask=None, weight='none'):
     1024        """
     1025        Average the Beams together.
     1026        Parameters:
     1027            mask:        An optional mask defining the region, where the
     1028                         averaging will be applied. The output will have all
     1029                         specified points masked.
     1030            weight:      Weighting scheme. 'none' (default), 'var' (1/var(spec)
     1031                         weighted), or 'tsys' (1/Tsys**2 weighted)
     1032        """
     1033        varlist = vars()
     1034        if mask is None:
     1035            mask = ()
     1036        s = scantable(self._math._averagebeams(self, mask, weight.upper()))
     1037        s._add_history("average_beam", varlist)
     1038        print_log()
     1039        return s
     1040
    10231041    def convert_pol(self, poltype=None):
    10241042        """
     
    13311349        return s
    13321350
    1333     def mx_quotient(self, mask = None, weight='median'):
     1351    def mx_quotient(self, mask = None, weight='median', preserve=True):
    13341352        """
    13351353        Form a quotient using "off" beams when observing in "MX" mode.
    13361354        Parameters:
     1355            mask:           an optional mask to be used when weight == 'stddev'
    13371356            weight:         How to average the off beams.  Default is 'median'.
     1357            preserve:       you can preserve (default) the continuum or
     1358                            remove it.  The equations used are
     1359                            preserve: Output = Toff * (on/off) - Toff
     1360                            remove:   Output = Toff * (on/off) - Ton
    13381361        """
    13391362        if mask is None: mask = ()
     
    13421365        preoff = scantable(self._math._mx_extract(self, 'off'))
    13431366        off = preoff.average_time(mask=mask, weight=weight, scanav=False)
    1344         print_log()
    13451367        from asapmath  import quotient
    1346         q = quotient(on, off)
     1368        q = quotient(on, off, preserve)
    13471369        q._add_history("mx_quotient", varlist)
     1370        print_log()
    13481371        return q
    13491372
Note: See TracChangeset for help on using the changeset viewer.