py/getm.py

   1 """How many times does a large mussa analysis redraw the same line.
   2 """
   3 import mussa
   4 m = mussa.Mussa()
   5 #m.load_mupa("examples/mck3test.mupa")
   6 #m.analyze()
   7 m.load("../../mussa_analyses/land/myogenin_w30_t21")
   8 path = m.paths()
   9
  10 count = 0
  11 tracks = {}
  12 for rp in path.refinedPathz:
  13   t0_id = rp.track_indexes[0]
  14   t1_id = rp.track_indexes[1]
  15   tid = (t0_id, t1_id)
  16   tracks[tid] = tracks.setdefault(tid, 0) + 1
  17   count += 1
  18
  19 print "count", count
  20 print "tracks", len(tracks.keys())