699

**in order**

alias basicyy="/Applications/basicyy.app/Contents/MacOS/basicyy --sam"

script_lib_dir = /Users/micah/ops/eeprom_effective/80_samlib
script_lib_file = /Users/micah/maven/mavlib_rev555.bas

class TMFile 
{
public: 
    // todo - I'm forgetting to use const in these methods 
    TMFile(std::string filename); 
    std::string filename(); 
    std::string absolute_directory(); 
    std::string directory_string();
    std::string get_mcal_filename(); 
    long long int file_size(); 
    int tid(); 
    int length();     
    bool exists(); 
    TMPacket get_pkt(int index);

    double t0_in_unix_time() const; 

private: 
    QString filename_;
    TmMeta tmMeta_;
    MomGse::MomTelemetryFilePtr momaTlmFilePtr_; 
    MomGse::MarkerCache markerCache_;
    MomGse::MomaTimestampResolver& timestampResolver_; 
};

macx {
    INCLUDEPATH += /Library/Frameworks/Python.framework/Versions/3.4/include/python3.4m 
    LIBS +=        /Library/Frameworks/Python.framework/Versions/3.4/lib/libpython3.4.dylib
    INCLUDEPATH += /Library/Frameworks/QtCore.framework/Versions/4/Headers
}
linux-g++ {
    # When Python3 is installed on wherever, it should be installed with
    # "./configure --enable-shared". This is necessary so that the library files
    # generated contain position independent (-fPIC) code. 
    INCLUDEPATH += /usr/local/include/python3.4m
    LIBS += /usr/local/lib/python3.4/config-3.4m/libpython3.4m.a
}

macx {
    QMAKE_LFLAGS_PLUGIN -= -dynamiclib
    QMAKE_LFLAGS_PLUGIN += -bundle
}
linux-g++ {
    QMAKE_LFLAGS_PLUGIN += -fPIC
}

# Bundles have no strictly enforced file extension, but swig expects the 
# libraries it uses to be of the form _<module name>{=html}.so , so rename the output
# file. (Don't forget the leading underscore. It's a common pitfall.)
macx {
    QMAKE_POST_LINK = "mv libc699util.dylib _c699util.so"
}
linux-g++ {
    QMAKE_POST_LINK = "mv libc699util.so _c699util.so"
}

[MOMA_TEST_DATA]
<other lines>
tm_database=/Users/mpallone/momagse/TMDef/MOM_TM_Database.txt

[MOMA_TEST_DATA]
<other lines>
tm_database=/Users/mpallone/momagse/TMDef/Test_MOM_TM_Database.txt

make clean
qmake -spec macx-g++ c699util.pro # For CentOS, I can just run "qmake"
make release
cat top_of_c699util_dot_py.py c699util.py > temp_c699util.py && mv temp_c699util.py c699util.py
./test_c699util.py
./copy_into_699util.py # Only do this if test_c699util.py passes all tests! 

/*****************************************************************************
 * TMPacket definition 
 *****************************************************************************/
class TMPacket 
{
public:
    TMPacket(GseLib::pbytes packetBytesPtr, int pktIndex,
             const MomGse::MarkerCache* markerCachePtr,
             const MomGse::MomaTimestampResolver& timestampResolver); 

    std::string description();  // this method was just added

/*****************************************************************************
 * TMPacket Implementation 
 *****************************************************************************/

TMPacket::TMPacket(GseLib::pbytes packetBytesPtr, int pktIndex,
                   const MomGse::MarkerCache* markerCachePtr,
                   const MomGse::MomaTimestampResolver& timestampResolver) 
    : packetBytesPtr_(packetBytesPtr),
      pktIndex_(pktIndex),
      markerData_(markerCachePtr->getMarker(pktIndex, &hasMarkerData_)),
      unixTimestamp_(timestampResolver.getPktUnixTime(pktIndex)),
      relativeTimestamp_(timestampResolver.getPktRelativeTime(pktIndex))
{
    // no further initialization needed 
}    

std::string description()
{
    return std::string("this is a description of a packet); 
}

/*****************************************************************************
 * TMPacket definition 
 *****************************************************************************/
class TMPacket 
{
public:
    TMPacket(GseLib::pbytes packetBytesPtr, int pktIndex,
             const MomGse::MarkerCache* markerCachePtr,
             const MomGse::MomaTimestampResolver& timestampResolver); 

    std::string description(); // this method was just added
    %pythoncode %{
        def __str__(self): 
            return self.description()
    %}

class TMFile 

    <other TMFile methods go here>{=html} 

    %pythoncode %{
        def test_name(self): 
            return self.directory_string().split('-', maxsplit=5)[5]
    %}

    %pythoncode %{
        def message_log(self): 
            from tmread import MessageLog
            return MessageLog(self)
    %}

(py699) gs699-mpallone:c699util mpallone$ otool -L _c699util.so
_c699util.so:
   /Library/Frameworks/Python.framework/Versions/3.4/Python (compatibility version 3.4.0, current version 3.4.0)
   QtGui.framework/Versions/4/QtGui (compatibility version 4.8.0, current version 4.8.4)
   QtCore.framework/Versions/4/QtCore (compatibility version 4.8.0, current version 4.8.4)
   QtNetwork.framework/Versions/4/QtNetwork (compatibility version 4.8.0, current version 4.8.4)
   /usr/lib/libstdc++.6.dylib (compatibility version 7.0.0, current version 104.1.0)
   /usr/lib/libSystem.B.dylib (compatibility version 1.0.0, current version 1226.10.1)
   /usr/lib/libgcc_s.1.dylib (compatibility version 1.0.0, current version 915.0.0)

   #!/usr/bin/env python3 

   #==============================================================================
   # Packet data extraction through 699util's tmread module: 
   #==============================================================================
   import tmread 

   tmfile = tmread.get_tmfile('/Users/mpallone/momadata/fm/2016-01-15-14.06.47-30410-MSFUNC/tm.mom.m3.s3')

   wrp_current_datapoints = []
   wrp_current_hkid = 800 
   pkt_17_count = 0 

   for pkt in tmfile(17): 
       pkt_17_count += 1 
       wrp_current_datapoints.extend(pkt[wrp_current_hkid])

   # Print out the data: 
   print("pkt_17_count:", pkt_17_count)
   print() 
   print("WRP current sample timestamp, WRP current sample science value:")
   for timestamp, value in wrp_current_datapoints: 
       print(timestamp, 't', value)

   #==============================================================================
   # Packet data extraction through c699util: 
   #==============================================================================
   import c699util 

   c_tmfile = c699util.construct_tmfile('/Users/mpallone/momadata/fm/2016-01-15-14.06.47-30410-MSFUNC/tm.mom.m3.s3')

   wrp_current_datapoints = []
   wrp_current_hkid = 800 
   pkt_17_count = 0 

   for pkt in c_tmfile(17): 
       pkt_17_count += 1 
       wrp_current_datapoints.extend(pkt[wrp_current_hkid])

   # Print out the data: 
   print("pkt_17_count:", pkt_17_count)
   print() 
   print("WRP current sample timestamp, WRP current sample science value:")

   #for timestamp, value in wrp_current_datapoints: 
   #    print(timestamp, 't', value)

   # Note the change here! 
   for datapoint in wrp_current_datapoints: 
       print(datapoint.relative_timestamp, 't', datapoint.value)

   import c699util 
   c_tmfile = c699util.get_tmfile('/Users/mpallone/momadata/fm/2016-01-15-14.06.47-30410-MSFUNC/tm.mom.m3.s3')

   import c699util 
   c_tmfile = c699util.get_tmfile() # grab the tm.mom file in the cwd 
   c_tmfile = c699util.get_tmfile(12345) # Grab TID 12345 
   c_tmfile = c699util.get_tmfile("/Users/mpallone/momadata/fm/2016-01-15-14.06.47-30410-MSFUNC") # Grab the tm.mom file in the specified directory

   import c699util 
   c_tmfile = c699util.construct_tmfile('/Users/mpallone/momadata/fm/2016-01-15-14.06.47-30410-MSFUNC/tm.mom.m3.s3')

   #!/usr/bin/env python3 

   import tmread, c699util 
   tmfile = tmread.get_tmfile(30410) # Give it a TID and let it find the file for us 
   c_tmfile = c699util.construct_tmfile(tmfile.filename)

   # (alternatively, we could just do c_tmfile = c699util.get_tmfile(30410) 

   print("Absolute filename:", c_tmfile.filename())
   print("Absolute filename:", str(c_tmfile))
   # Absolute filename: /Users/mpallone/momadata/fm/2016-01-15-14.06.47-30410-MSFUNC/tm.mom.m3.s3

   print("Absolute directory:", c_tmfile.absolute_directory())
   # Absolute directory: /Users/mpallone/momadata/fm/2016-01-15-14.06.47-30410-MSFUNC

   # Parsed from the first GSE packet in the tm.mom file: 
   print("TID Directory name:", c_tmfile.directory_string())
   # TID Directory name: 2016-01-15-14.06.47-30410-MSFUNC

   print("Mass calibration file used:", c_tmfile.get_mcal_filename())
   # Mass calibration file used: 2016-01-17-23.48.39-30410.mcal

   print("File size in bytes:", c_tmfile.file_size())
   # File size in bytes: 48384623

   print("TID:", c_tmfile.tid())
   # TID: 30410

   print("t0 in Unix time:", c_tmfile.t0_in_unix_time())
   # t0 in Unix time: 1452884878.3959727

   print("Number of packets:", c_tmfile.length())
   print("Number of packets:", len(c_tmfile))
   # Number of packets: 70232

   print("TMFile actually exists on the filesystem:", c_tmfile.exists())
   print("TMFile actually exists on the filesystem:", bool(c_tmfile))
   # TMFile actually exists on the filesystem: True

   pkt_at_index_10 = c_tmfile.get_pkt(10)
   pkt_at_index_10 = c_tmfile[10]

   # We can iterate through c_tmfile and get packet objects: 
   pkt_count = 0 
   for pkt in c_tmfile: 
       pkt_count += 1 

   # The old tm_filter routines also work for c699util: 
   pkt_17_count = 0 
   for pkt in c_tmfile(17): 
       pkt_17_count += 1 

   pkt_count = 0 
   for pkt in c_tmfile(type_filter=17,      # Only yield packets of type 17 
                       marker_filter=2500,  # Only yield packets at marker 2500 
                       start_time=1000,     # Relative time = 1000 seconds 
                       stop_time=2000):     # Relative time = 2000 seconds 
       pkt_count += 1 

   # If the directory name is
   #     '2016-01-15-14.06.47-30410-MSFUNC'
   # then the test name is 
   #     'MSFUNC'
   print("Test name:", c_tmfile.test_name())

   # tmread.MessageLog is not deprecated, and is fully compatible with c699util:
   msg_log = c_tmfile.message_log()

   meta_markers = c_tmfile.get_meta_markers()
   for meta_marker in meta_markers: 
       print(meta_marker.startRelativeTimeInSecs)
       print(meta_marker.endRelativeTimeInSecs)
       print(meta_marker.metaMarkerId)
       print(meta_marker.index)
       print(meta_marker.msg)

   print("model:", c_tmfile.model())
   # model: fm

   print("FSW timestamp at of t0 packet:", c_tmfile.t0_in_sclk())
   # FSW timestamp at of t0 packet: 19351744

   #!/usr/bin/env python3 

   import tmread, c699util 
   tmfile = tmread.get_tmfile(30410) # Give it a TID and let it find the file for us 
   c_tmfile = c699util.construct_tmfile(tmfile.filename)

   AUX_AMP_DAC_SET_HKID = 61 
   RF_AMP_MON_HKID      = 62

   fields = tmread.extract_tmfields(c_tmfile, [AUX_AMP_DAC_SET_HKID, RF_AMP_MON_HKID])

   # Note that we still get c699util.TMDatapoint objects here, which are different
   # from 699util's tmread datapoints!
   sum_ = 0 
   count = 0 
   for datapoint in fields[RF_AMP_MON_HKID]: 
       relative_timestamp = datapoint.relative_timestamp 
       unix_timestamp = datapoint.unix_timestamp
       raw_timestamp = datapoint.fsw_timestamp
       sum_ += datapoint.value 
       count += 1 

   if count != 0: 
       print("RF Amp mon avg:", sum_ / count)

## TMPacket

Just like TMFile, c699util provides TMPacket objects which are very
similar to tmread.TMPacket objects:

   #!/usr/bin/env python3 

   import tmread, c699util 
   tmfile = tmread.get_tmfile(30410) # Give it a TID and let it find the file for us 
   c_tmfile = c699util.construct_tmfile(tmfile.filename)

   gse_pkt = c_tmfile[0]       # Pkt type 8, GSE generated 
   pkt_22 = c_tmfile[65]       # Pkt type 22 
   msg_log_pkt = c_tmfile[85]  # Pkt type 8 

   print("Unmasked type of gse_pkt:", gse_pkt.full_type())
   # Unmasked type of gse_pkt: 136

   print("Unmasked type of pkt_22:", pkt_22.full_type())
   # Unmasked type of pkt_22: 22

   # This is analagous to tmread.TMPacket.type 
   print("gse_pkt with gse-generated bit masked out:", gse_pkt.get_type())
   # gse_pkt with gse-generated bit masked out: 8

   print("pkt_22 length, including header and checksum:", pkt_22.total_length())
   # pkt_22 length, including header and checksum: 224
   print("pkt_22 length, just the data:", pkt_22.length())
   # pkt_22 length, just the data: 212

   print("pkt_22 checksum:", pkt_22.checksum())
   print("pkt_22 checksum:", hash(pkt_22))
   # pkt_22 checksum: 1379474455

   print("Is pkt_22 gse created?", pkt_22.gse_created())
   # Is pkt_22 gse created? False
   print("Is gse_pkt gse created?", gse_pkt.gse_created())
   # Is gse_pkt gse created? True

   print("pkt_22 sequence number:", pkt_22.sequence())
   # pkt_22 sequence number: 146

   # Yes, I'm aware that 'sclk' is a misnomer 
   print("pkt_22 FSW timestamp:", pkt_22.sclk())
   # pkt_22 FSW timestamp: 19540225
   print("pkt_22 FSW timestamp, in seconds:", pkt_22.sclk_in_seconds())
   # pkt_22 FSW timestamp, in seconds: 1954.0225

   print("pkt_22 index:", pkt_22.index())
   # pkt_22 index: 65

   print("pkt_22 marker ID:", pkt_22.marker_id())
   # pkt_22 marker ID: 0

   # Packets before the first marker packet have a marker index of 0 
   print("pkt_22 marker index:", pkt_22.xina_marker_index())
   # pkt_22 marker index: 0

   print("pkt_22 unix timestamp, in seconds:", pkt_22.unix_timestamp())
   # pkt_22 unix timestamp, in seconds: 1452884897.2440727

   print("pkt_22 relative timestamp, in seconds:", pkt_22.relative_timestamp())
   # pkt_22 relative timestamp, in seconds: 18.848099946975708

   print("pkt_22 marker text:", pkt_22.marker_text())
   # pkt_22 marker text: <no marker text, but you get the idea>

   ###########################################################
   # 
   # (this doesn't actually work for TID 30410, but I think you get the picture)
   # 
   # HKID 100 = phase cycle readback HKID
   phase_cycle_readback_raw_data = pkt_22.all_raw_data(100)
   phase_cycle_readback_eng_data = pkt_22.all_eng_data(100)

   phase_cycle_readback_sci_data = pkt_22.all_data(100)
   phase_cycle_readback_sci_data = pkt[100]

   datapoint = phase_cycle_readback_sci_data[0]
   print("Datapoint unix timestamp, in seconds:", datapoint.unix_timestamp)
   print("Datapoint relative timestamp, in seconds:", datapoint.relative_timestamp)
   print("Datapoint value:", datapoint.value)

   first_phase_cycle_readback_science_sample = pkt_22.get_first_value(100) 
   # 
   # 
   ###########################################################    

   print("Does gse_pkt have an earlier timestamp than pkt_22?", gse_pkt < pkt_22)
   # Does gse_pkt have an earlier timestamp than pkt_22? True

   print("Raw message log:", msg_log_pkt.raw_message())
   print("Sanitized message log:", msg_log_pkt.message())
   # (output omitted for brevity)

   print("raw bytes:", pkt_22.raw_packet())
   # b'x16x92x00xd4x01*)x01xfaxf3 x02x00...(etc). 

   meta_markers = pkt_22.get_meta_markers() 
   for meta_marker in meta_markers: 
       print(meta_marker.startRelativeTimeInSecs)
       print(meta_marker.endRelativeTimeInSecs)
       print(meta_marker.metaMarkerId)
       print(meta_marker.index)
       print(meta_marker.msg)

   #!/usr/bin/env python3 

   import c699util 
   c_tmfile = c699util.get_tmfile(30410)

   scan_count = 0 
   for scan in c699util.get_scans(c_tmfile): 
       scan_count += 1 

   print("Number of scans:", scan_count)

   #!/usr/bin/env python3 

   import c699util 

   # Not needed 
   # c_tmfile = c699util.construct_tmfile(tmfile.filename)

   # Notice that this factory function takes the absolute filename, just like 
   # construct_tmfile
   # science_cache = c699util.construct_science_data_cache('/Users/mpallone/momadata/fm/2016-01-15-14.06.47-30410-MSFUNC/tm.mom.m3.s3')

   science_cache = c699util.get_science_data_cache(30410) 

   #!/usr/bin/env python3 

   import c699util 

   science_cache = c699util.get_science_data_cache(30410)

   last_scan = science_cache[-1]

   print("Last scan data:")

   print("Unix timestamp, in seconds:", last_scan.unix_timestamp())
   # Unix timestamp, in seconds: 1452890221.7956727

   print("Relative timestamp, in seconds:", last_scan.relative_timestamp())
   # Relative timestamp, in seconds: 5343.399699926376

   print("marker ID:", last_scan.marker_id())
   # marker ID: 2050

   print("marker Index:", last_scan.xina_marker_index())
   # marker Index: 136

   print("marker text:", last_scan.marker_text())
   # marker text: SCN 602 START Emission Current Source B

   list_of_ion_counts = last_scan.counts() 
   print("total ion count:", sum(list_of_ion_counts))
   print("total ion count:", last_scan.total_ion_count())
   # total ion count: 40463

   print("Highest count in the last_scan:", last_scan.highest_count())
   # Highest count in the last_scan: 708

   print("Bin number containing the highest count:", last_scan.highest_count_bin_num())
   # Bin number containing the highest count: 1710

   print("Lowest count:", last_scan.lowest_count())
   # Lowest count: 0

   # There should be one entry for every 100 bins 
   rf_status_entry_list = last_scan.rf_status_entries()
   rf_entry = rf_status_entry_list[0]

   print("First RF Entry start bin:", rf_entry.start_bin)
   # First RF Entry start bin: 0

   print("First RF Entry RF Amp V:", rf_entry.rf_amp_v)
   # First RF Entry RF Amp V: 66.67613871395588

   print("First RF Entry RF Amp DAC V:", rf_entry.rf_amp_dac_v)
   # First RF Entry RF Amp DAC V: 70.00512770935893

   print("First RF Entry Aux Amp DAC V:", rf_entry.aux_amp_dac_v)
   # First RF Entry Aux Amp DAC V: 1.0998702980577946

   print("RF Amp Start:", last_scan.rf_amp_start())
   # RF Amp Start: 66.67613871395588

   print("RF Amp End:", last_scan.rf_amp_end())
   # RF Amp End: 993.6579284607433

   print("RF Amp DAC Start:", last_scan.rf_amp_dac_start())
   # RF Amp DAC Start: 70.00512770935893

   print("RF Amp DAC End:", last_scan.rf_amp_dac_end())
   # RF Amp DAC End: 996.3826513544906

   print("Aux Amp DAC Start:", last_scan.aux_amp_dac_start())
   # Aux Amp DAC Start: 1.0998702980577946

   print("Total scan time, in seconds:", last_scan.total_scan_time())
   # Total scan time, in seconds: 0.05000000074505805

   print("Peak counts per second:", last_scan.peak_counts_per_sec())
   # Peak counts per second: 70799998.9449978

   print("Pkt source type (26 or 27):", last_scan.src_pkt_type())
   # Pkt source type (26 or 27): 27

   print("Number of scans that went into this scan (typically 10 for sum packets):")
   print(last_scan.num_scans())
   # Number of scans that went into this scan (typically 10 for sum packets):
   # 10

   print("Does this scan have a scan status packet?", last_scan.has_scan_status())
   # Does this scan have a scan status packet? True

   print("Number of scan status packets:", last_scan.num_scan_status_pkts())
   # Number of scan status packets: 10

   print("Does this scan have a preceding SEB HK packet?", last_scan.has_seb_hk_pkt())
   # Does this scan have a preceding SEB HK packet? True

   print("Is this scan a dark count scan?", last_scan.is_dark_count_scan())
   # Is this scan a dark count scan? False

   print("Sum of the ion counts for bins 500-600:")
   print(last_scan.selected_ion_count_sum_by_bin(500,600))
   # Sum of the ion counts for bins 500-600:
   # 3981

   print("Sum of the ion counts between masses 138-140:")
   print(last_scan.selected_ion_count_sum_by_mass(138,140))
   # Sum of the ion counts between masses 138-140:
   # 78

   print("Maximum ion count between bins 500-600:")
   print(last_scan.selected_ion_count_max_by_bin(500,600))
   # Maximum ion count between bins 500-600:
   # 537

   print("Maximum ion count between masses 138-140:")
   print(last_scan.selected_ion_count_max_by_mass(138,140))
   # Maximum ion count between masses 138-140:
   # 21

   print("seb_bin_time_readback:", last_scan.seb_bin_time_readback())
   # seb_bin_time_readback: 10.000000149011612

   print("seb_bin_time_readback_in_seconds:", last_scan.seb_bin_time_readback_in_seconds())
   # seb_bin_time_readback_in_seconds: 1.0000000149011612e-05

   print("seb_special_scan_id:", last_scan.seb_special_scan_id())
   # seb_special_scan_id: 602

   print("seb_raw_scan_value:", last_scan.seb_raw_scan_value())
   # seb_raw_scan_value: 127

   print("seb_scan_value:", last_scan.seb_scan_value())
   # seb_scan_value: 49.80392238497734

   print("seb_rf_freq_hz:", last_scan.seb_rf_freq_hz())
   # seb_rf_freq_hz: 988400.0

   print("seb_last_rf_amp_mon:", last_scan.seb_last_rf_amp_mon())
   # seb_last_rf_amp_mon: 993.6579284607433

   print("seb_rf_amp_dac_set:", last_scan.seb_rf_amp_dac_set())
   # seb_rf_amp_dac_set: 996.3826513544906

   print("seb_aux_amp_dac_set:", last_scan.seb_aux_amp_dac_set())
   # seb_aux_amp_dac_set: 7.297474631876097

   print("seb_counter_sel:", last_scan.seb_counter_sel())
   # seb_counter_sel: 1

   print("seb_bin_count_set:", last_scan.seb_bin_count_set())
   # seb_bin_count_set: 5000

   print("seb_bin_count_reg:", last_scan.seb_bin_count_reg())
   # seb_bin_count_reg: 20000

   print("seb_phase_cycle_rdbk:", last_scan.seb_phase_cycle_rdbk())
   # seb_phase_cycle_rdbk: 3

   print("seb_rf_step_set:", last_scan.seb_rf_step_set())
   # seb_rf_step_set: 0.1870379802200142

   print("seb_rf_step_reg:", last_scan.seb_rf_step_reg())
   # seb_rf_step_reg: -0.038366765173336245

   print("seb_aux_step_set:", last_scan.seb_aux_step_set())
   # seb_aux_step_set: 0.0012397955291065842

   print("seb_aux_step_reg:", last_scan.seb_aux_step_reg())
   # seb_aux_step_reg: 0.0006103608758678569

   print("seb_ionization_time_in_milliseconds:", last_scan.seb_ionization_time_in_milliseconds())
   # seb_ionization_time_in_milliseconds: 4.999999888241291

   print("seb_scan_em1:", last_scan.seb_scan_em1())
   # seb_scan_em1: -2069.6365661583154

   print("seb_scan_em1_raw:", last_scan.seb_scan_em1_raw())
   # seb_scan_em1_raw: 197

   print("seb_scan_em2:", last_scan.seb_scan_em2())
   # seb_scan_em2: 994.4167251074475

   print("seb_scan_em2_raw:", last_scan.seb_scan_em2_raw())
   # seb_scan_em2_raw: 0

   print("seb_src_a_foc_a_dac_set:", last_scan.seb_src_a_foc_a_dac_set())
   # seb_src_a_foc_a_dac_set: -39.6078431372549

   print("seb_src_b_foc_a_dac_set:", last_scan.seb_src_b_foc_a_dac_set())
   # seb_src_b_foc_a_dac_set: -39.6078431372549

   print("seb_src_a_foc_a_dac_reg:", last_scan.seb_src_a_foc_a_dac_reg())
   # seb_src_a_foc_a_dac_reg: -67.45098039215686

   print("seb_src_b_foc_a_dac_reg:", last_scan.seb_src_b_foc_a_dac_reg())
   # seb_src_b_foc_a_dac_reg: -67.45098039215686

   print("meta markers:", last_scan.get_meta_markers())
   # output omitted due to laziness, it's the same interface that's documented in the TMFile and TMPacket sections

   print("m/z values:", last_scan.m_over_z_values())
   # (long tuple of doubles)

   print("list of TICs:", last_scan.total_ion_counts())
   # list of TICs: (3680, 3882, 4092, 4100, 4147, 4119, 4103, 4021, 4214, 4105)

   print("cTIC:", last_scan.corrected_total_ion_count()) 
   # cTIC: 1715842

   print("cTIC Factor:", last_scan.corrected_total_ion_count_factor())
   # cTIC Factor: 42.40522135001076

   print("cTIC Factors:", last_scan.corrected_total_ion_count_factors())
   # cTIC Factors: ()

   print("Packet version:", last_scan.science_data_pkt_version())
   # Packet version: 1

   print("seb_is_fil_a_on:", last_scan.seb_is_fil_a_on())
   # seb_is_fil_a_on: False

   print("seb_is_fil_b_on:", last_scan.seb_is_fil_b_on())
   # seb_is_fil_b_on: True

   print("Scan mode:", last_scan.scan_mode_as_string())
   # Scan mode: EI 

   print("SEB sequence packet indexes:", last_scan.seb_sequence_pkt_indexes())
   # SEB sequence packet indexes: (68463, 68464, 68465)

   print("EM on times (ms):", last_scan.em_on_times_in_ms())
   # EM on times (ms): (600.1000052131712, 0.0)

   #!/usr/bin/env python3 

   import tmread, c699util 
   tmfile = tmread.get_tmfile(30410) # Give it a TID and let it find the file for us 
   c_tmfile = c699util.construct_tmfile(tmfile.filename)

   AUX_AMP_DAC_SET_HKID = 61 
   RF_AMP_MON_HKID      = 62

   #------------------------------------------------------------------------------
   # TMFile interface differences 
   #------------------------------------------------------------------------------

   # Getting a TMFile: 
   tmfile = tmread.get_tmfile(30410) # Can be passed TID or absolute filename
   c_tmfile = c699util.construct_tmfile(tmfile.filename) # requires absolute filename

   filename =   tmfile.filename 
   filename = c_tmfile.filename() 

   dir_string =   tmfile.directory_string 
   dir_string = c_tmfile.directory_string() 

   t0_in_unix_time =   tmfile.start_time().unix   # Returns TMTimestamp object
   t0_in_unix_time = c_tmfile.t0_in_unix_time() # Returns a float 

   msg_log =   tmfile.message_log 
   msg_log = c_tmfile.message_log() 

   #------------------------------------------------------------------------------
   # TMPacket interface differences 
   #------------------------------------------------------------------------------

   pkt   = tmfile[-1]
   c_pkt = c_tmfile[-1]

   # Packet type with GSE bit masked out: 
   pkt_type =   pkt.type 
   pkt_type = c_pkt.get_type() 

   # Full packet type: 
   full_pkt_type =   pkt.full_type 
   full_pkt_type = c_pkt.full_type() 

   total_length =   pkt.total_length 
   total_length = c_pkt.total_length() 

   # Doesn't included the 8 byte header and 4 byte checksum, unlike total_length 
   data_length =   pkt.length 
   data_length = c_pkt.length() 

   checksum =   pkt.checksum 
   checksum = c_pkt.checksum() 

   gse_created =   pkt.gse_created
   gse_created = c_pkt.gse_created() 

   sequence_number =   pkt.sequence 
   sequence_number = c_pkt.sequence() 

   fsw_timestamp =   pkt.raw_time 
   fsw_timestamp = c_pkt.sclk() 

   unix_timestamp =   pkt.timestamp.unix     # TMTimestamp object 
   unix_timestamp = c_pkt.unix_timestamp()   # float 

   relative_timestamp =   pkt.timestamp.relative   # TMTimestamp object 
   relative_timestamp = c_pkt.relative_timestamp() # float 

   marker_text =   pkt.marker_text
   marker_text = c_pkt.marker_text()  

   marker_id =   pkt.mkid 
   marker_id = c_pkt.marker_id() 

   msg =   pkt.message 
   msg = c_pkt.message() 

   #
   # Housekeeping data interface: 
   # 
   data_dict = tmread.extract_tmfields(tmfile, [800])
   hkid_800_datapoints = data_dict[800]
   first_datapoint = hkid_800_datapoints[0]

   timestamp = first_datapoint.ts 
   value = first_datapoint.val 

   data_dict = tmread.extract_tmfields(c_tmfile, [800])
   hkid_800_datapoints = data_dict[800]
   first_datapoint = hkid_800_datapoints[0]

   unix_timestamp = first_datapoint.unix_timestamp
   relative_timestamp = first_datapoint.relative_timestamp 
   value = first_datapoint.value 

   #------------------------------------------------------------------------------
   # Scan interface differences 
   #------------------------------------------------------------------------------

   scan =   list(tmread.get_scans(tmfile))[-1]
   c_scan = list(c699util.get_scans(c_tmfile))[-1]

   result =   scan.timestamp.unix        # Returns a TMTimestamp object 
   result = c_scan.unix_timestamp();     # Returns a float 

   result =   scan.timestamp.relative       # Returns a TMTimestamp object 
   result = c_scan.relative_timestamp();    # Returns a float 

   result =   scan.marker_text 
   result = c_scan.marker_text(); 

   result =   scan.scanpoints   # list of MomaScanPoint objects 
   result = c_scan.counts();    # list of integers 

   result =   scan.total_ion_count
   result = c_scan.total_ion_count(); 

   result =   scan.highest_count
   result = c_scan.highest_count(); 

   result =   scan.highest_count_bin_num 
   result = c_scan.highest_count_bin_num(); 

   result =   scan.rf_amp_start 
   result = c_scan.rf_amp_start(); 

   result =   scan.rf_amp_end 
   result = c_scan.rf_amp_end(); 

   result =   scan.rf_amp_dac_start 
   result = c_scan.rf_amp_dac_start();

   result =   scan.rf_amp_dac_end 
   result = c_scan.rf_amp_dac_end(); 

   result =   scan.aux_amp_dac_start 
   result = c_scan.aux_amp_dac_start(); 

   result =   scan.aux_amp_dac_end 
   result = c_scan.aux_amp_dac_end();     

   result =   scan.total_scan_time 
   result = c_scan.total_scan_time(); 

   result =   scan.peak_counts_per_sec
   result = c_scan.peak_counts_per_sec();

   result =   scan.src_pkt_type 
   result = c_scan.src_pkt_type();

   result =   scan.num_scans 
   result = c_scan.num_scans(); 

   result = bool(scan.seb_state.num_scan_status_pkt_updates)
   result = c_scan.has_scan_status(); 

   result =   scan.seb_state.num_scan_status_pkt_updates
   result = c_scan.num_scan_status_pkts();     

   result =   scan.seb_state.scan_bin_time_readback 
   result = c_scan.seb_bin_time_readback()

   result =   scan.seb_state.scan_bin_time_readback_in_seconds 
   result = c_scan.seb_bin_time_readback_in_seconds()

   result =   scan.seb_state.special_scan_id 
   result = c_scan.seb_special_scan_id()

   result =   scan.seb_state.raw_scan_value 
   result = c_scan.seb_raw_scan_value()

   result =   scan.seb_state.scan_value 
   result = c_scan.seb_scan_value()

   result =   scan.seb_state.rf_freq 
   result = c_scan.seb_rf_freq_hz()

   result =   scan.seb_state.last_rf_amp_mon
   result = c_scan.seb_last_rf_amp_mon()

   result =   scan.seb_state.rf_amp_dac_set
   result = c_scan.seb_rf_amp_dac_set()

   result =   scan.seb_state.aux_amp_dac_set
   result = c_scan.seb_aux_amp_dac_set()

   result =   scan.seb_state.counter_sel
   result = c_scan.seb_counter_sel()

   result =   scan.seb_state.bin_count_set
   result = c_scan.seb_bin_count_set()

   result =   scan.seb_state.bin_count_reg
   result = c_scan.seb_bin_count_reg()

   result =   scan.seb_state.phase_cycle_rdbk
   result = c_scan.seb_phase_cycle_rdbk()

   result =   scan.seb_state.rf_step_set
   result = c_scan.seb_rf_step_set()

   result =   scan.seb_state.rf_step_reg
   result = c_scan.seb_rf_step_reg()

   result =   scan.seb_state.rf_step_reg
   result = c_scan.seb_aux_step_set()

   result =   scan.seb_state.aux_step_reg
   result = c_scan.seb_aux_step_reg()

   result =   scan.seb_state.ionization_time
   result = c_scan.seb_ionization_time_in_milliseconds()

   result =   scan.seb_state.scan_em1
   result = c_scan.seb_scan_em1()

   result =   scan.seb_state.scan_em2
   result = c_scan.seb_scan_em2()

   result =   scan.seb_state.srcA_foc_a_dac_set
   result = c_scan.seb_src_a_foc_a_dac_set()

   result =   scan.seb_state.srcB_foc_a_dac_set
   result = c_scan.seb_src_b_foc_a_dac_set()

   result =   scan.seb_state.srcA_foc_a_dac_reg
   result = c_scan.seb_src_a_foc_a_dac_reg()

   result =   scan.seb_state.srcB_foc_a_dac_reg
   result = c_scan.seb_src_b_foc_a_dac_reg()

   result = bool(scan.seb_state.num_hk_pkt_updates)
   result = c_scan.has_seb_hk_pkt(); 

   result =   scan.seb_state.is_emon_a
   result = c_scan.seb_is_emon_a(); 

   result =   scan.seb_state.is_emon_b
   result = c_scan.seb_is_emon_b(); 

   result =   scan.count_ions(min_bin=min_bin, max_bin=max_bin)
   result = c_scan.selected_ion_count_sum_by_bin(min_bin, max_bin); 

   result =   scan.count_ions(min_mass=min_mass, max_mass=max_mass)
   result = c_scan.selected_ion_count_sum_by_mass(min_mass, max_mass); 

   result =   scan.max_ion_count(min_bin=min_bin, max_bin=max_bin)
   result = c_scan.selected_ion_count_max_by_bin(min_bin, max_bin);

   result =   scan.max_ion_count(min_mass=min_mass, max_mass=max_mass)
   result = c_scan.selected_ion_count_max_by_mass(min_mass, max_mass);

#!/usr/bin/env python3 

from tmread import get_tmfile 

tmfile = get_tmfile("/Users/mpallone/momadata/etu/tm.mom.m2")

directory = tmfile.absolute_directory() 
print("directory:", directory)
# output: directory: /Users/mpallone/momadata/etu

start_time = tmfile.start_time()
print("start_time:", start_time) # returns the timestamp of the earliest packet 
# output: start_time: 0.0

file_size = tmfile.file_size() # returns the number of bytes of the tm.mom file 
print("file_size:", file_size)
# output: file_size: 610

tid = tmfile.tid()
print("tid:", tid)
# output: tid: 7515

#!/usr/bin/env python3 

from tmread import get_tmfile 

tmfile = get_tmfile("/Users/mpallone/momadata/etu/tm.mom.m2")

num_pkts = 0 
for pkt in tmfile: 
    num_pkts += 1 

print("num_pkts:", num_pkts)
# output: num_pkts: 544

(py699)gs699-mpallone:~ mpallone$
(py699)gs699-mpallone:~ mpallone$ # Starting in the home directory: 
(py699)gs699-mpallone:~ mpallone$
(py699)gs699-mpallone:~ mpallone$ pwd
/Users/mpallone
(py699)gs699-mpallone:~ mpallone$
(py699)gs699-mpallone:~ mpallone$
(py699)gs699-mpallone:~ mpallone$ # We can use the 'tid' command to jump to a 
(py699)gs699-mpallone:~ mpallone$ # TID directory: 
(py699)gs699-mpallone:~ mpallone$
(py699)gs699-mpallone:~ mpallone$ tid 7421
(py699)gs699-mpallone:~ mpallone$ 
(py699)gs699-mpallone:~ mpallone$ 
(py699)gs699-mpallone:2014-12-12-15.33.44-07421-steves_test mpallone$ pwd
/Users/mpallone/momadata/etu/2014-12-12-15.33.44-07421-steves_test
(py699)gs699-mpallone:~ mpallone$
(py699)gs699-mpallone:~ mpallone$
(py699)gs699-mpallone:2014-12-12-15.33.44-07421-steves_test mpallone$ python
Python 3.3.2 (v3.3.2:d047928ae3f6, May 13 2013, 13:52:24) 
[GCC 4.2.1 (Apple Inc. build 5666) (dot 3)] on darwin
Type "help", "copyright", "credits" or "license" for more information.
>>>
>>> # We can use the tmread module from the Python interpreter, from any 
>>> # directory: 
>>> 
>>> from tmread import get_tmfile
>>>
>>> tmfile = get_tmfile()
>>> tmfile.tid()
7421
>>>
>>>
>>> tmfile = get_tmfile(30100)
>>> tmfile.tid()
30100
>>>
>>>
>>> tmfile = get_tmfile("30101")
>>> tmfile.tid()
30101
>>>
>>>
>>>
>>> tmfile = get_tmfile("/Users/mpallone/momadata/qsm/2014-12-12-09.48.18-20000-woa-204-wrp-load-test")
>>> tmfile.tid()
20000
>>>
>>> count = 0
>>> for pkt in tmfile:
...     count += 1
...
>>> count
702

#!/usr/bin/env python3 

from tmread import get_tmfile 

tmfile = get_tmfile("/Users/mpallone/momadata/etu/tm.mom.m2")

for pkt in tmfile: 
    print(pkt.raw_time)    # raw_time is the timestamp in the packet header 
    print(pkt.type)        # 8 for message logs, 7 for digital status packets, etc. 
    print(pkt.gse_created) # always False for genuine MOMA packets
    print(pkt.sequence)    # sequence number of the packet
    print(pkt.length)      # length of the packet (not including the 8 byte header or 4 byte checksum) 
    print(pkt.mkid)        # Marker ID of the packet 
    print(pkt.marker_text) # text of the marker packet that preceded the current packet (if such a packet exists)

#!/usr/bin/env python3 

from tmread import get_tmfile 

tmfile = get_tmfile("/Users/mpallone/momadata/fm/2015-02-10-12.26.48-30005-woa-245-wrp-test/tm.mom.m2")

for pkt in tmfile: 

    # This isn't the best way to extract data. This is just a TMPacket example.
    # See extract_tmfields() for a better way. 
    if pkt.type == 17: # MAIF packet 
        ps_temp = pkt.all_data(806) # HKID 806 is PS_TEMP 

        # ps_temp is now a list of (timestamp, value) tuples 
        first_datapoint = ps_temp[0]
        timestamp = first_datapoint.ts 
        value = first_datapoint.val

        print("timestamp:", timestamp)
        print("value:", value)
        # output: 
        # timestamp: -1.895219
        # value: 29.495738

        break

   (py699)gs699-mpallone:2015-02-25-18.49.20-30072-msfunc_rf_ramp_time mpallone$ tmmarker.py
      75.321    1000   Check RF Freq
      75.542    1010   Pressure Check
      96.347    1020   Filament on
     102.327    1040   Initialize SEB
     117.294    2500   SCN 13 START rf ramp time
     120.530    2501   SCN 13 END rf ramp time
   (py699)gs699-mpallone:2015-02-25-18.49.20-30072-msfunc_rf_ramp_time mpallone$ 
   (py699)gs699-mpallone:2015-02-25-18.49.20-30072-msfunc_rf_ramp_time mpallone$ 
   (py699)gs699-mpallone:2015-02-25-18.49.20-30072-msfunc_rf_ramp_time mpallone$ python
   Python 3.3.2 (v3.3.2:d047928ae3f6, May 13 2013, 13:52:24) 
   [GCC 4.2.1 (Apple Inc. build 5666) (dot 3)] on darwin
   Type "help", "copyright", "credits" or "license" for more information.
   >>> 
   >>> from tmread import *
   >>> tmfile = get_tmfile(30072)
   >>> markers = set()
   >>> 
   >>> 
   >>> # A string can be used as the marker filter. A dash indicates a range 
   >>> # of desired values. 
   >>> 
   >>> for pkt in tmfile(marker_filter="1010-1040"):
   ...     markers.add(pkt.mkid)
   ...
   >>>
   >>> markers
   {1040, 1010, 1020}
   >>>
   >>> 
   >>> 
   >>> # An integer can be used as a marker filter. 
   >>> 
   >>> markers = set()
   >>> for pkt in tmfile(marker_filter=1010):
   ...     markers.add(pkt.mkid)
   ...
   >>> markers
   {1010}
   >>>
   >>> 
   >>> 
   >>> # A comma separated string can be used to specify a list of marker IDs
   >>> # to filter by. Note that marker ID 1020 isn't present, unlike when 
   >>> # we filtered using "1010-1040". 
   >>> 
   >>> markers = set()
   >>> for pkt in tmfile(marker_filter="1010,1040"):
   ...     markers.add(pkt.mkid)
   ...
   >>> markers
   {1040, 1010}
   >>> 
   >>> 
   >>> 
   >>> 
   >>> # A list can be used to filter markers. 
   >>>
   >>> markers = set()
   >>> for pkt in tmfile(marker_filter=[2500, 2501]):
   ...     markers.add(pkt.mkid)
   ...
   >>> markers
   {2500, 2501}
   >>>
   >>>
   >>> 
   >>> # A generator can be used to filter markers. 
   >>>
   >>> markers = set()
   >>> for pkt in tmfile(marker_filter=range(1000,2000)):
   ...     markers.add(pkt.mkid)
   ...
   >>>
   >>> markers
   {1000, 1040, 1010, 1020}

   >>> 
   >>> 
   >>> # The type_filter can be an integer: 
   >>>
   >>> types = set()
   >>> for pkt in tmfile(type_filter=25):
   ...     types.add(pkt.type)
   ...
   >>> types
   {25}
   >>>
   >>> 
   >>> # Or the type filter can be a list (really, any iterable container): 
   >>>
   >>> types = set()
   >>> for pkt in tmfile(type_filter=[7, 8, 9]):
   ...     types.add(pkt.type)
   ...
   >>>
   >>> types
   {8, 9, 7}
   >>> 
   >>> 
   >>> 
   >>> # Or the type filter can be a generator: 
   >>>
   >>> types = set()
   >>> for pkt in tmfile(type_filter=range(7,10)):
   ...     types.add(pkt.type)
   ...

   >>>
   >>> types
   {8, 9, 7}
   >>>

   >>> # marker_filter and type_filter can be combined: 
   >>> 
   >>> for pkt in tmfile(marker_filter=2501, type_filter=10):
   ...     print(pkt.index, pkt.type, pkt.mkid)
   ...
   5646 10 2501
   5649 10 2501
   5652 10 2501
   # etc.

#!/usr/bin/env python3 

from tmread import TMArgumentParser 
import tmread 

import os 

def main():

    # python scripts are usually run inside of a TID directory, so pretend to do 
    # that here. 
    os.chdir("/Users/mpallone/momadata/fm/2015-02-25-20.51.11-30076-msfunc_gc_ldi_ec_increasing/")

    parser = TMArgumentParser(description="This title-string shows up when the -h flag is used.") 

    # This allows the user to specify --unix to see Unix time, --relative to see 
    # relative time, etc. When "tmpacketObject.timestamp" is accessed, it will 
    # default to whatever flag the user specified.
    parser.add_time_modes_group() 

    # This allows the user to specify --raw to see raw telemetry values, --eng to 
    # see telemetry values converted to engineering units, etc. Housekeeping 
    # objects will take on whatever format this flag specifies. 
    parser.add_data_modes_group() 

    # This allows the user to specify "-m 50" or "-m 50-100" to see marker 50 
    # packets or to see marker packets between marker 50 and marker 100 
    parser.add_marker_argument() 

    # After setting up the parser, let actually parse the arguments: 
    args = parser.parse_args() 

    # 
    # *This* is the preferred way to get a TMFile object: 
    # 
    tmfile = args.tmfile 
    if not tmfile: 
        return tmread.ReturnCodes.missing_tmfile 

    # Do something with the tmfile object 
    print(tmfile.absolute_directory())
    # output: /Users/mpallone/momadata/fm/2015-02-25-20.51.11-30076-msfunc_gc_ldi_ec_increasing

main()

#!/usr/bin/env python3 

from tmread import MomaPktIds 

print(MomaPktIds.fsw_time_pkt_id == 2)   # True 
print(MomaPktIds.digital_hk_pkt_id == 7) # True 
print(MomaPktIds.msg_log_pkt_id == 8)    # True 

The full definition (as of 05/26/15):

class MomaPktIds(enum.IntEnum):
    """Enumeration of packet IDs used in MOMA."""  
    memory_dump_pkt_id            = 1 
    fsw_time_pkt_id               = 2 
    digital_hk_pkt_id             = 7 
    msg_log_pkt_id                = 8 
    marker_pkt_id                 = 9 

    rsim_nominal_pkt_id           = 10 

    fill_pkt_pkt_id               = 11 
    gc_science_and_hk_pkt_id      = 14 
    laser_pulse_tlm_pkt_id        = 15 
    laser_periodic_status_pkt_id  = 16 
    maif_1_sec_pkt_id             = 17 
    maif_1_min_pkt_id             = 18 
    micro_pirani_pressure_pkt_id  = 19 
    seb_ack_pkt_id                = 21 
    seb_hk_pkt_id                 = 22 
    seb_sequence_mem_dump_pkt_id  = 23 
    seb_dds_mem_dump_pkt_id       = 24 
    seb_scan_status_pkt_id        = 25 
    seb_science_data_pkt_id       = 26 
    seb_summed_histogram_pkt_id   = 27 
    seb_combined_histogram_pkt_id = 28 

    rsim_fake_timestamp_pkt_id    = 74

#!/usr/bin/env python3 

from tmread import get_tmfile, extract_tmfields

tmfile = get_tmfile("/Users/mpallone/momadata/fm/2015-02-25-20.51.11-30076-msfunc_gc_ldi_ec_increasing/tm.mom.m2")

list_of_hkids = [802] # Arbitrarily extracting HKID 802 

dictionary_of_datapoints = extract_tmfields(tmfile, list_of_hkids) 

first_datapoint = dictionary_of_datapoints[802][0]

timestamp = first_datapoint.ts 
value = first_datapoint.val 

print(timestamp, value) # Output: 1.012687 1.498488

   >>> from tmread import moma_tmdirs
   >>> from tmread import moma_tids
   >>>
   >>> 
   >>> # moma_tids is a dictionary where the key is the TID (which can be an int
   >>> # or a string), and the value is the directory of that TID on the user's 
   >>> # computer. 
   >>> 
   >>> moma_tids[7421]
   '/Users/mpallone/momadata/etu/2014-12-12-15.33.44-07421-steves_test'
   >>>
   >>> moma_tids["7421"]
   '/Users/mpallone/momadata/etu/2014-12-12-15.33.44-07421-steves_test'
   >>>
   >>> 
   >>> 
   >>> # moma_tmdirs is simply a list of TID directories on the user's computer.
   >>>
   >>> moma_tmdirs[0]
   '/Users/mpallone/momadata/etu/2014-07-24-14.16.31-00268-'
   >>>
   >>>
   >>> moma_tmdirs[5]
   '/Users/mpallone/momadata/etu/2014-07-29-17.00.38-00273-MEB-GC-Integration-ETU'
   >>> moma_tmdirs[100]
   '/Users/mpallone/momadata/etu/2014-09-07-15.35.24-07078-firing_laser'

#!/usr/bin/env python3
"""A script display start time from a telemetry file.

Many telemetry file processing programs use time since start of file.
This program displays that time.

Author: Eric Raaen
Date:   Feb 1, 2012

"""
from time import asctime, ctime, gmtime
import tmread

def main():
    """Main function for t0.py."""
    parser = tmread.TMArgumentParser(description='Display start time of a '
                                                 'telemetry file')
    args = parser.parse_args()

    t0 = args.tmfile.start_time()
    gmt = gmtime(t0.unix)
    print("SCLK (raw):", t0.sclk)
    print("     ctime:", t0.unix)
    print("       UTC:", asctime(gmt))
    print("Local time:", ctime(t0.unix))

    return 0

if __name__ == '__main__':
    exit(main()) 

    from timeit import Timer

    t1 = Timer("for i in xrange(100): i % 2")
    t2 = Timer("for i in xrange(100): i & 1")
    # The "not" tests show what happens when interpreting
    # the result as a boolean
    t3 = Timer("for i in xrange(100): not i % 2")
    t4 = Timer("for i in xrange(100): not i & 1")

    print "Checking for odd parity with `mod`:\t%.4f" % t1.timeit()
    print "Checking for odd parity with `and`:\t%.4f" % t2.timeit()
    print "Checking for even parity with `mod`:\t%.4f" % t3.timeit()
    print "Checking for even parity with `and`:\t%.4f" % t4.timeit()

Checking for odd parity with `mod`:    6.5617
Checking for odd parity with `and`: 5.3778
Checking for even parity with `mod`:    8.4417
Checking for even parity with `and`:    7.4086

Checking for odd parity with `mod`:    0.2556
Checking for odd parity with `and`: 0.2312
Checking for even parity with `mod`:    1.7576
Checking for even parity with `and`:    0.6614

Checking for odd parity with `mod`:    3.4480
Checking for odd parity with `and`: 1.9380
Checking for even parity with `mod`:    3.6050
Checking for even parity with `and`:    2.0440

    from timeit import Timer

    index1 = Timer("x = tpl[0]", "tpl = (5,)")
    unpack1 = Timer("x, = tpl", "tpl = (5,)")
    index2 = Timer("x = tpl[1]", "tpl = (5, 6)")
    unpack2 = Timer("y, x = tpl", "tpl = (5, 6)")
    index3 = Timer("x = tpl[2]", "tpl = (5, 6, 7)")
    unpack3 = Timer("y, y, x = tpl", "tpl = (5, 6, 7)")
    index4 = Timer("x = tpl[3]", "tpl = (5, 6, 7, 8)")
    unpack4 = Timer("y, y, y, x = tpl", "tpl = (5, 6, 7, 8)")

    list_index2 = Timer("[tpl[1] for tpl in tuples]", "tuples = [(i, i * i) for i in xrange(100)]")
    list_unpack2 = Timer("[y for x, y in tuples]", "tuples = [(i, i * i) for i in xrange(100)]")
    list_map2 = Timer("map(itemgetter(1), tuples)", "tuples = [(i, i * i) for i in xrange(100)]; from operator import itemgetter")

    times = 100000000

    print "Indexing vs. unpacking a 1-tuple:\t%.4f\t%.4f" % (index1.timeit(number=times), unpack1.timeit(number=times))
    print "Indexing vs. unpacking a 2-tuple:\t%.4f\t%.4f" % (index2.timeit(number=times), unpack2.timeit(number=times))
    print "Indexing vs. unpacking a 3-tuple:\t%.4f\t%.4f" % (index3.timeit(number=times), unpack3.timeit(number=times))
    print "Indexing vs. unpacking a 4-tuple:\t%.4f\t%.4f" % (index4.timeit(number=times), unpack4.timeit(number=times))
    print "Indexing vs. unpacking a list of 2-tuples:\t%.4f\t%.4f" % (list_index2.timeit(), list_unpack2.timeit())
    print "map() and itemgetter() (just for kicks):\t%.4f" % (list_map2.timeit())

Indexing vs. unpacking a 1-tuple:  5.0712  3.3939
Indexing vs. unpacking a 2-tuple:   5.7888  6.2801
Indexing vs. unpacking a 3-tuple:   6.1820  7.5976
Indexing vs. unpacking a 4-tuple:   7.1802  7.8219
Indexing vs. unpacking a list of 2-tuples:  8.6561  8.3513
map() and itemgetter() (just for kicks):    9.1651

Indexing vs. unpacking a 1-tuple:  0.2268  0.2279
Indexing vs. unpacking a 2-tuple:   0.2301  0.2302
Indexing vs. unpacking a 3-tuple:   0.2335  0.2320
Indexing vs. unpacking a 4-tuple:   0.2332  0.2344
Indexing vs. unpacking a list of 2-tuples:  1.2610  1.2698
map() and itemgetter() (just for kicks):    5.4586

Indexing vs. unpacking a 1-tuple:  0.6510  1.1520
Indexing vs. unpacking a 2-tuple:   0.9610  0.7800
Indexing vs. unpacking a 3-tuple:   0.8930  0.8330
Indexing vs. unpacking a 4-tuple:   1.0250  0.8070
Indexing vs. unpacking a list of 2-tuples:  36.4800 40.4600
map() and itemgetter() (just for kicks):    11.0170

    from timeit import Timer

    # Use 62 to prevent slowdown from long ints
    t1 = Timer("for i in xrange(62): 1 << i")
    t2 = Timer("for i in xrange(62): 2 ** i")

    times = 1000000

    print "Bit-shifting vs. Exponentation:\t%.4f\t%.4f" % (t1.timeit(number=times), t2.timeit(number=times))

Bit-shifting vs. Exponentation:    3.8654  8.3995

Bit-shifting vs. Exponentation:    0.2184  2.0279

Bit-shifting vs. Exponentation:    2.6870  17.0960

{
 "total_pkts": 6285,
 "earliest_pkt_unix_timestamp": 1418415862.637699,
 "earliest_pkt_sclk": 74300701,
 "earliest_pkt_unix_timestamp_no_rsim_correction": 1418408433.567599,
 "earliest_pkt_sclk_no_rsim_correction": 1,
 "last_pkt_unix_timestamp": 1418420434.78211,
 "last_pkt_sclk": 4970713,
 "first_time_pkt_rover_seconds": 471680511,
 "first_time_pkt_rover_subseconds": 56741,
 "first_time_pkt_moma_ticks": 151142,
 "timestamp_override_table": []
}

:   RxCreditStat =\> swts_reg2(29 to 31), 0 - 7\
:   TxCreditStat =\> swts_reg2(26 to 28), 0 - 7\
:   char_out =\> swts_reg1(0 to 9), display as hex\
:   cleared_rx_fifo =\> swts_reg2(25), 0 - 1\
:   credit =\> swts_reg2(24), 0 - 1\
:   data_in_rdy =\> swts_reg2(23), 0 - 1\
:   data_out_rdy =\> swts_reg2(22), 0 - 1\
:   errwaitstate =\> swts_reg2(19), 0 - 1\
:   link_down =\> swts_reg2(18), 0 - 1\
:   link_version =\> swts_reg2(10 to 17), 0 - ?\
:   rx_almost_empty =\>swts_reg2(9), 0 - 1\
:   rx_fifo_ff =\> swts_reg2(8), 0 - 1\
:   status_bus =\> swts_reg5(0 to 15), display in hex \-- breakout
    below\
:   tick_out =\> swts_reg2(6), 0 - 1\
:   time_out =\> swts_reg5(16 to 23), display in hex\
:   time_sent =\> swts_reg2(5), 0 - 1\
:   tx_almost_full =\> swts_reg2(4), 0 - 1\
:   tx_full_flag =\> swts_reg2(3) 0 - 1\

:   2 - ERROR WAIT
:   3 - READY
:   4 - STARTED
:   5 - PRE CONNECTING
:   6 - CONNECTING
:   7 - RUN\

{  
  "gse":{  
     "default_mcal": "2015-01-01-00.00.00-identity.mcal",
     "mission":"moma",
     "model":"etu",
     "patch_files":[  
        "m1",
        "m2",
        "m3"
     ],
     "tid":534,
     "tid_dir":"2015-09-29-13.46.05-00534-",
     "tmfile_name":"tm.mom.m1.m2.m3",
     "version":1
  }
}

{  
  "gse":{  
     "mission":"moma",
     "model":"etu",
     "patch_files":[  
        "m1",
        "m2",
        "m3"
     ],
     "tid":534,
     "tid_dir":"2015-09-29-13.46.05-00534-",
     "tmfile_name":"tm.mom.m1.m2.m3",
     "version":1
  }
}

Script Name	Description
checksum.py	Calculate the Fletcher's checksum of input or a list of existing files.
datavol.py	Estimate the data volumen of a telemetry file.
faradaycup.py	Correlate faraday cup scans to faraday cup currents.
feiget.py	Download data from FEI.
fei2tm.py	Convert EDRs or DPs to a telemetry file.
gcms.py	Print a tab-delimited table of GCMS data organized by both mass and time.
gcsci.py	Print GC TCD data in a tab-delimited table.
gcsciplot.py	Plot GC TCD data, with or without TCD offset correction.
getcsvreports	Script for PDL to download MSL engineering reports and sort their contents by sol.
hkdump.py	Print the values in a SAM type-10 housekeeping packet.
hot.py	Create a tab-delimited CSV file showing the cumulative number of tics the heaters have been on.
htrduty.py	Plot the duty cycle and temperature of a given SAM heater.
listfei.py	List and sort the DPs or EDRs in the current working directory.
massrange.py	Plot the masses and bands being scanned at a given time.
movefei.py	Move files downloaded from FEI into the appropriate TID folder.
ontime.py	Calculate how long a given component was "on" (e.g., how long a valve open, or how long an HK value was above a certain threshold)
phd.py	Plot the pulse-height distribution at a given marker ID.
powercut.py	Cut a power profile into two independent power profiles.
ql_ega.py	Generate the plots for a quicklook.
runall.py	Run a command or python code in every telemetry directory.
sclk2scet.py	Convert an SCLK timestamp to SCET time.
t0.py	Display the start time of a telemetry file in various formats.
timestamp.py	Print the start time of the file in a human-readable format.
tm2rdr.py	Create a Reduced Data Record from a telemetry file.
tmdb.py	A method used in backtics to get the axis name for a database entry.
tmdiff.py	A diff-like program for telemetry files.
tmdump.py	Produce a hex dump of telemetry file, packet by packet.
tmexcerpt.py	Create a telemetry file containing only a certain marker or packet type.
tmfields.py	Print a tab-delimited table of housekeeping data.
tmglue.py	Concatenate two telemetry files together.
tmkeys.py	Search for and display the telemetry database.
tmmarker.py	Print a list of the markers and marker IDs in the telemetry file.
tmmsg.py	Print the contents of the message log packets.
tmplot.py	Plot housekeeping data.
tmsanitize.py	Remove malformed data from a packet.
tmscan.py	Plot a DAC scan.
tmsequence.py	Print a report detailing whether packet sequence numbers increase normally.
tmstat.py	Print the statistics of a given HK value (min, median, maximum, mean).
tmsummary.py	Print a summary of each packet in a telemetry file.
triplets.py	Print a tab-delimited table of GCMS data organized by time.
valvestat.py	Print the valve operations (like the tmvalves viewer in SAM Data View).

Mass Flow	WRP Back Pressure	Instant WRP Current
MP_PRESS < 0.5	4 Torr	0.55 A
MP_PRESS < 0.5	6 Torr	0.65 A
MP_PRESS < 0.5	8 Torr	0.75 A
MP_PRESS > 0.5	4 Torr	0.70 A
MP_PRESS > 0.5	6 Torr	0.80 A
MP_PRESS > 0.5	8 Torr	0.90 A

Mass Flow	WRP Back Pressure	1-Min Avg WRP Current
MP_PRESS < 0.5	4 Torr	0.45 A
MP_PRESS < 0.5	6 Torr	0.55 A
MP_PRESS < 0.5	8 Torr	0.65 A
MP_PRESS > 0.5	4 Torr	0.60 A
MP_PRESS > 0.5	6 Torr	0.70 A
MP_PRESS > 0.5	8 Torr	0.80 A

Field	Required	Description
SUB_SYSTEM	Yes	The sub system that the telemetry value comes from
NAME	Yes	The name of the telemetry value
DESCRIPTION	Yes	A description of the telemetry value. Useful if the name is non-descriptive.
ADC #	Yes	Deprecated but still required right now, use 0 as the value
DATA_ID	Yes	The telemetry value's unique identifier
PACKET_TYPE	Yes	The type/APID of the packet that contains the telemetry value
PACKET_SUBTYPE	Yes	Deprecated but still required right now, use 0 as the value
DECODE_TYPE	Yes	The decode type, which dictates how the value is decoded, and also how decode ID 1-4 are used
DECODE_ID1	Yes	The purpose of these depends on the decode type
RAW_FMT	Yes	How to format/display an unconverted telemetry value, i.e. the raw bits/bytes in the packet. Follows the C style printf formatting
ENG_FMT	Yes	How to format/display a telemetry value converted to its engineering representation. Follows the C style printf formatting
ENG_EQ_TYPE	Yes	The type of engineering transformation to apply to the raw value, and also dictates how the ENG_K1-ENG_K2 are used
ENG_K1	Yes	Their usage depends on the ENG_EQ_TYPE, but is usually the coefficients of a 1st order polynomial. K1 = constant, K2 = coefficient of first order
SCI_UNITS	Yes	The units of the science transformation
SCI_FMT	Yes	How to format/display a telemetry value converted to its science representation
SCI_EQ_TYPE	Yes	The type of science transformation to apply to the raw value, and also dictates how the SCI_K1-SCI_K8 are used
SCI_FROM_ENG	Yes	If Y, then engineering transformation is applied to the raw value before science transformation. Otherwise, if N, then science transformation is applied directly to the raw value.
SCI_K1	Yes	Values used in the science transformation. How they are used is determined by the SCI_EQ_TYPE.
DISCRETE_LABELS	No	Allows user to define a mapping of raw values to textual labels
MUX_CH	Yes	Deprecated but still required right now, use 0 as the value
LIMIT_TYPE_1	Yes	The type of the first limit check
YELLOW_LIMIT_1	Yes	The first limit check's yellow limit
RED_LIMIT_1	Yes	The first limit check's red limit
LIMIT_TYPE_2	Yes	The type of the second limit check
YELLOW_LIMIT_2	Yes	The second limit check's yellow limit
RED_LIMIT_2	Yes	The second limit check's red limit
EXTRA_1	Yes	An extra field that is used for various purposes.

699

699 Software

General Information

Update Mailing List

Bug Reporting

Installation

Mac OS X

System Requirements

Windows

Individual Software Pages

Miscellaneous

Mac Gnuplot Installation

Mac lxml Installation

699util

Installation Instructions

Bug Reports

List of Scripts

Version History

1.00.00

1.00.01

1.00.02

1.00.03

1.00.04

1.00.05

1.00.06

1.00.07

1.01.00

1.01.01

1.02.00

1.02.01

1.02.02

1.03.00

1.03.01

1.03.02

1.04.00

1.05.00

1.05.01

2.00.00

2.00.01

2.00.02

699util User Scripts

Alerts

MOMA Alerts

General Error Check

MAIF APV Comm

Spacewire Comm Error

MOMA Safe Event

WRP Too Fast

WRP Failed To Start

WRP Pressure Spike

Dark Counts

Laser Burst Energy Variation

Laser Pulses Mismatch

Unknown SEB sequences

Filament VMON Pegged

Expected Pulse Count Mismatch

Gap in Scans

Old Packet Format

Missing Scan Status Packets

Dropped Packets

Both EMs On

Sub-system Communication Errors

WRP Instantaneous High Current

WRP 1-Min Average High Current

Archived SAGE Versions

BASIC Script Tools

Basic Editor

Mac OS X Installation

Windows Installation

Syntax Checker (basicyy)

Mac OS X Installation

Bugzilla

Creating An Account

c699util Developer's Guide

SWIG

How c699util and the 699 C++ code work together

How SWIG and qmake work together

How to add a feature to c699util

Getting c699util Source Code

Setting up test_c699util.py

Mac `lxml` Installation