Source of task_appendantab.py - casa6

 
                keys = StringIO()
                continue
            continue
        
        # Create the subtable for syscal if selected in params
        if self.append_tsys: 
            self.create_syscal_subtable()
        if self.append_gc: 
            self.create_gaincurve_subtable()
​
    def get_spw_ranges(self):
        # Get ranges for each spectral window
        try:
            tb.open(f'{self.vis}/SPECTRAL_WINDOW')
            #print(self.spws)
            for spw in self.spws:
                chan_freq = tb.getcol('CHAN_FREQ')[:,spw]
                min_freq = min(chan_freq) / (1e6)
                max_freq = max(chan_freq) / (1e6)
                self.spw_ranges[spw] = (min_freq, max_freq)
            tb.close()
        except Exception as e:
            print("ERROR FAILED TO GET SPW FREQ RANGES")
​
    def check_which_spw(self, freq):
        in_spw = []
        for spw, ranges in self.spw_ranges.items():
            if ranges[0] >= freq[0] and ranges[1] <= freq[1]:
                in_spw.append(spw)
        # If there is no spw that this belongs to default to what we were setting this value to before
        if not in_spw:
            print("WARNING: NO CORRESPONDING SPW FOUND FOR FREQ {}".format(freq))
            return [-1]
        else:
            return in_spw
​
    def process_tsys(self, infp, keys):
        # Get the antenna name
        antenna_name = find_antenna(keys[0], ['SRC/SYS'])
        # Skip if no antenna name found
        if not antenna_name:
            print('ANTENNA missing from TSYS group')
            skip_values(infp)
            return
        try:
            antenna = get_antenna_index(antenna_name, self.ant_names)
        except:
            print('Antenna %s not present in the Measurement Set. %s values will be ignored.' % antenna_name)
            skip_values(infp)
            return
​
        try:
            self.get_spw_ranges()
            spw_ids = self.check_which_spw(keys['FREQ'])
        except Exception as e:
            spw_ids = [-1]
            print('ANTAB HAS NO FREQ INFORMATION DEFAULTING TO FREQ ID -1')
        
        keys = dict(keys[0])
        spws = []
        spwmap = {}
        
        if 'INDEX' in keys:
            update_map(self.pols, spws, spwmap, keys['INDEX'])
        if 'INDEX2' in keys:
            update_map(self.pols, spws, spwmap, keys['INDEX2'])
            pass
Show more
                        tsys[pol].append(value)
                        continue
                    continue
                if source != -1:
                    time_val = mjd_sec + (days * 86400)
                    self.data.time.append(time_val)
                    self.data.time_interval.append(0.0)
                    self.data.antenna_no.append(antenna)
                    self.data.tsys_1.append(tsys['R'])
                    self.data.tsys_2.append(tsys['L'])
                    self.data.freqid(spw_ids)
​
                    pass
                pass
​
            if line.strip().endswith('/'):
                break
            continue
        return
​
    def get_spw_ranges(self):
        # Get ranges for each spectral window
        try:
            tb.open(f'{self.vis}/SPECTRAL_WINDOW')
            #print(self.spws)
            for spw in self.spws:
                chan_freq = tb.getcol('CHAN_FREQ')[:,spw]
                min_freq = min(chan_freq) / (1e6)
                max_freq = max(chan_freq) / (1e6)
                self.spw_ranges[spw] = (min_freq, max_freq)
            tb.close()
        except Exception as e:
            print("ERROR FAILED TO GET SPW FREQ RANGES")
​
    def check_which_spw(self, freq):
        in_spw = []
        for spw, ranges in self.spw_ranges.items():
            if ranges[0] >= freq[0] and ranges[1] <= freq[1]:
                in_spw.append(spw)
        # If there is no spw that this belongs to default to what we were setting this value to before
        if not in_spw:
            print("WARNING: NO CORRESPONDING SPW FOUND FOR FREQ {}".format(freq))
            return [-1]
        else:
            #print('SELECTED SPWS: ', in_spw)
            return in_spw
​
    def process_gc(self, infp, keys):
        antenna_name = find_antenna(keys[0], self.gain_keys)
        if not antenna_name:
            print('Antenna missing from GAIN group')
            skip_values(infp)
            return
        try:
            antenna = get_antenna_index(antenna_name, self.ant_names)
        except:
            print('Antenna %s not present in FITS-IDI file' % antenna_name)
Show more
            # if the antenna has been visited already with the same time stamp then ignore
            if self.data.antenna_no[i] in existing_data and self.data.time[i] in existing_data[self.data.antenna_no[i]]['times'] and not self.overwrite:
                # This data is already in the subtable so ignore
                #print("antenna: {} has data present for the time {}. The antab data for this value will be ignored".format(data.antenna_no[i], data.time[i]))
                to_skip.append(True)
                continue
            else:
                to_skip.append(False)
​
            for j in self.spws:
                # Check that this antenna is in this spectral window
                if j not in self.data.freqid[i] and self.data.freqid[i] != [-1]:
                    continue
                if self.data.antenna_no[i] in existing_data and self.data.time[i] in existing_data[self.data.antenna_no[i]]['times'] and self.overwrite:
                    # Get the row and information to overwrite the row from the ms
                    row = existing_data[self.data.antenna_no[i]]['row']
                    to_skip.append(True)
​
                    self.replacements[row] = {'ANTENNA':self.data.antenna_no[i], 'FEED_ID':0, 'INTERVAL':self.data.time_interval[i],
                                        'SPW_ID':int(j), 'TIME':self.data.time[i], 'TSYS':[[],[]]}
​
                    continue
​
casa6

Source

Add shortcut
